Кто-нибудь знает, как bootstrap линейная модель со смешанным эффектом, построенная с помощью функции lme (пакет nlme)?
Полугодовой go Я смог сделать это с помощью "bootstrap" функция из пакета "lmeresampler", но в настоящее время она не работает с теми же данными и кодом. Что случилось с этими пакетами?
mod<-lme(PC1~SA+Y+CMIave+MATave,
random=~1|PLOT_ID,data)
mixed.bootmod<-lmeresampler::bootstrap(mod,function(.)fixef(.),type="parametric",B=1000)
Это не запускается, показывая ошибку, что
"некоторые bootstrap запускаются неудачно", приводя к NA
Впрочем, lmer работает
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)