如何绘制lme模型的各个轨迹

Question

我有示例数据和模型

x<-rep(seq(0,100,by=1),10)
y<-15+2*rnorm(1010,10,4)*x+rnorm(1010,20,100)
id<-NULL
for (i in 1:10){
  id<-c(id, rep(i,101))}
dtfr<-data.frame(x=x,y=y, id=id)
library(nlme)
with (dtfr, summary(lme((y)~x,random=~1+x|id, na.action=na.omit  )))
model.mx<-with (dtfr, (lme((y)~x,random=~1+x|id, na.action=na.omit  )))
pd<-predict(model.mx, newdata=data.frame(x=0:100),level=0)
with (dtfr, plot(x, y))
lines(0:100,predict(model.mx, newdata=data.frame(x=0:100),level=0), col="darkred", lwd=7)

如何提取每个ID的建模截距和斜率并绘制每个ID的各个轨迹？

Answer 1

不确定你想要做什么，因为你的系数几乎相同：

> coef(model.mx)
   (Intercept)        x
1     54.88302 19.18001
2     54.88298 19.18000
3     54.88299 19.18000
4     54.88299 19.18000
5     54.88302 19.18001
6     54.88300 19.18000
7     54.88301 19.18000
8     54.88300 19.18000
9     54.88299 19.18000
10    54.88300 19.18000

也许您的真实数据会给您带来更多不同的结果。如果是这种情况，我会在abline电话中使用mapply：

with (dtfr, plot(x, y))
mapply(abline,a=coef(model.mx)[,1],b=coef(model.mx)[,2], col=1:10)

这是结果。由于所有系数几乎相同，因此将线条绘制在彼此之上。你只能看到最后一个。