我使用 nlme 包創建了一個多級回歸模型，現在我想說明為某些患者獲得的回歸線（不幸的是，我不能將 geom_smooth 與 nlme 一起使用）。
因此，使用該模型，我在不同時間（date_day）和這里為兩個患者（ID1 和 ID2）獲得了以下預測值（predicted_value）。
df <- data.frame (ID = c (rep (1, 10), rep(2, 10)),
date_day = c (7:16, 7:16),
predicted_value = c (33, 33, 33, 33, 33, NA, 34, NA, NA, NA,
55, NA, NA, 53.3, NA, NA, 51.6, NA, 50.5, NA))
ID date_day predicted_value
1 1 7 33.0
2 1 8 33.0
3 1 9 33.0
4 1 10 33.0
5 1 11 33.0
6 1 12 NA
7 1 13 34.0
8 1 14 NA
9 1 15 NA
10 1 16 NA
11 2 7 55.0
12 2 8 NA
13 2 9 NA
14 2 10 53.3
15 2 11 NA
16 2 12 NA
17 2 13 51.6
18 2 14 NA
19 2 15 50.5
20 2 16 NA
現在我想為這些患者中的每一個畫出回歸線。所以我嘗試了以下
ggplot(df%>% filter(ID %in% c("1", "2")))
aes(x = date_day, y = predicted_value)
geom_point(shape = "circle", size = 1.5, colour = "#112446", na.rm = T)
geom_line(aes(y = predicted_value), na.rm = T, size = 1)
theme_minimal()
facet_wrap(vars(ID))
scale_x_continuous(name="days", limits=c(7, 16))
scale_y_continuous(name="predicted values", limits=c(0, 60))
但我以以下情節結束： 患者 1 ：線路中斷，患者 2 根本沒有線路。我該如何解決？
非常感謝
uj5u.com熱心網友回復：
謝謝@BenBolker，確實改變了第一行
ggplot(df%>% filter(ID %in% c("1", "2")))
到
ggplot(na.omit(df)%>% filter(ID %in% c("1", "2")))
允許解決作業
轉載請註明出處，本文鏈接：https://www.uj5u.com/caozuo/451747.html