2

サンプルデータとモデルがあります

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 の個々の軌跡をプロットするにはどうすればよいですか?

4

1 に答える 1

3

すべての係数がほぼ同じであるため、何をしたいのかわからない:

> 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)

これが結果です。すべての係数がほぼ同じであるため、線は互いに重ねてプロットされます。最後のものしか見えません。 ここに画像の説明を入力

于 2013-05-26T20:02:19.330 に答える