有谁知道如何引导由 lme 函数(nlme 包)构建的线性混合效应模型?
半年前,我可以通过“lmeresampler”包中的“bootstrap”函数来做到这一点,但它目前不能使用相同的数据和代码运行。这些包裹怎么了?
mod<-lme(PC1~SA+Y+CMIave+MATave,
random=~1|PLOT_ID,data)
mixed.bootmod<-lmeresampler::bootstrap(mod,function(.)fixef(.),type="parametric",B=1000)
#This doesn't run, showing an error that
> "some bootstrap runs failed", producing NAs
#However, lmer works
mod1<-lmer(PC1~SA+Y+CMIave+MATave+(1|PLOT_ID),data)
mixed.bootmod<-bootMer(mod,function(.)fixef(.),nsim=1000)
mixed.bootmod<-lmeresampler::bootstrap(mod,function(.)fixef(.),type="parametric",B=1000)