Question

我跑了一个lmer模型。我能够像这样为每个人制作正常的格子图：

xyplot(asiatot~time|ID, data=S8_earlyNA)

现在我希望看到我的模型（L4_edit）中的拟合值对于每个超出这些图的个体。注意：在我的原始模型中，我没有包含＆＃34; na.action = na.exclude＆＃34;，但是当我尝试保存拟合值时，我会收到错误＆＃34;替换有4048行，数据有4557＆＃34; （我假设因为缺少一些数据点）

L4_edit <- lmer(asiatot ~ time+ ASIMPC01_A+ SPLVL1_A+ NaStatus3+
              NaStatus3*time+ F+ (time|ID), S8, na.action=na.exclude)

S8$fittedL4edit <- fitted(L4_edit)

S8_edit <- subset(S8, S8$Sodiumblockers.early==1)

ggplot(S8_edit)+aes(x=S8_edit$time, y=S8_edit$asiatot) + 
  geom_line(aes(y=S8_edit$fittedL4edit), size=0.5) +
  geom_point() + facet_wrap("ID")

这产生了我想要的东西，除了一个问题：可以理解的是，这些人的一些拟合线被打破了......但有时候我不明白为什么这条线丢失了这么多，或者为什么呢没有创建2行。鉴于数据，似乎缺少很多行。请参阅下面的图片链接（例如ID＃378或380）1。

谢谢！

lmer trellis plot：拟合（）行的部分缺失

0 个答案: