以下移植物取自这篇文章。我是引导程序的新手，并试图为带有R boot包的线性混合模型实现参数、半参数和非参数引导程序。
代码
这是我的R代码：
library(SASmixed)
library(lme4)
library(boot)
fm1Cult <- lmer(drywt ~ Inoc + Cult + (1|Block) + (1|Cult), data=Cultivation)
fixef(fm1Cult)
boot.fn <- function(data, indices){
data <- data[indices, ]
mod <- lmer(drywt ~ Inoc + Cult + (1|Block) + (1|Cult), data=data)
fixef(mod)
}
set.seed(12345)
Out <- boot(data=Cultivation, statistic=boot.fn, R=99)
Out
问题
- 如何对带有
boot包的混合模型进行参数、半参数和非参数引导？ - 我想我正在我的代码中为混合模型进行非参数引导。
我找到了这些幻灯片，但找不到 R 包merBoot。任何想法我可以在哪里得到这个包裹。任何帮助将不胜感激。提前感谢您的帮助和时间。