绘制分段增长曲线

如何解决绘制分段增长曲线

我正在尝试绘制类似于 this first plot 的分段增长曲线。我使用了单独的斜率编码方案并在时间 2 处放置了一个断点

|时间 | 0 | 1 | 2 | 5 | 10 | 15 | 20|

|时间1 | 0 | 1 | 2 | 2 | 2 | 2 | 2 |

|时间2 | 0 | 0 | 0 | 1 | 2 | 3 | 4 |

我使用以下代码来创建我的增长模型

m1 <- lmer(sdmtwr ~ time1 + time2 + (time1 | id) + (0 + time2 | id),data = SDMT,REML = FALSE)

我也在探索与使用以下代码的 2 级分类预测器的交互

m2 <- lmer(sdmtwr ~ (time1 + time2)*edu + (time1 | id) + (0 + time2 | id),REML = FALSE)

我尝试使用 ggplot2、sjPlot 和效果包创建绘图均无济于事，由于编程经验有限，我不知所措。我只能分别为基线模型和交互模型绘制段。

如果有人可以提供有关适当代码的帮助，我将不胜感激！

编辑：这是 dput 摘要（编辑长度以显示 edu、time1 和 time2）

> dput(sdmt)
structure(list(id = c(3L,3L,4L,5L,6L,28L,62L,108L,119L,120L,132L,148L,175L,178L,201L,253L,327L,336L,343L,360L,366L,366L),time = c(0L,2L,10L,15L,20L,0L,1L,15L
),sdmtwr = c(20L,24L,18L,19L,9L,17L,41L,33L,27L,29L,31L,26L,32L,40L,42L,46L,38L,14L,25L,45L,34L,30L,71L,52L,51L,50L,55L,39L,35L,37L,36L,44L,43L,49L,48L,NA,47L),time_seg1 = c(0,2,1,2),time_seg2 = c(0,3,4,3),ed_dich = structure(c(2L,1L),.Label = c("< HS",">= HS"),class = "factor")),row.names = c(NA,-80L),class = "data.frame")

解决方法

我认为您想要的是分段线性样条。您可以使用截断的幂基函数来完成此操作。在您的模型中，您将包括时间和一个时间为 2 的函数，如果时间大于 2，否则为 0。这使得分段线性函数在时间 = 2 处相遇。您可以按如下方式在模型中执行此操作：

library(lme4)
mod <- lmer(sdmtwr ~ time + I(ifelse(time > 2,time-2,0)) + 
              (1 |id),data=tmp,REML=TRUE)

然后，您可以使用 ggpredict() 包中的 ggeffects 函数来生成绘图：

library(ggeffects)
g <- ggpredict(mod,"time")
plot(g)

注意：我无法在对时间变量产生随机影响的情况下运行它，但如果有更多数据，也许你就能让它运行起来。