基本上,我正在运行nlme::lme
一个缺失值。该模型非常适合na.action=na.omit
,但是为什么拟合/残差/系数的名称似乎都移动了一行?
## Generate data ---------------------
X1=gl(2,4)
X2=gl(2,2,8)
Y=rnorm(8)
dat=data.frame(Y=Y,X1=X1,X2=X2)
dat
## missing value -------------
mis.dat=dat
mis.dat[3,"Y"]=NA
mis.dat
> mis.dat
Y X1 X2
1 -0.06845332 1 1
2 0.89169085 1 1
3 NA 1 2
4 1.88997449 1 2
5 0.95912879 2 1
6 -0.64049400 2 1
7 -0.23354948 2 2
8 -0.66869350 2 2
## Fit model -----------------------
model=nlme::lme(Y~1,random=~1|X1/X2,data=mis.dat,na.action=na.omit)
summary(model)
## Notie the names -------------------
fitted(model)
> fitted(model)
1/1 1/1 2/1 2/1 2/2 2/2 <NA>
0.67179438 0.67179438 0.67179439 0.02855517 0.02855517 0.02855517 0.02855517
attr(,"label")
[1] "Fitted values"
#model$coef$random
#resid(model)
注意拟合值的名称?不应该在第三个位置有一个 1/2,之后的名称向右移动一个位置,从而消除 NA 吗?