Следующие прививки взяты из этой статьи . Я новичок в начальной загрузке и пытаюсь реализовать параметрическую, полупараметрическую и непараметрическую загрузку начальной загрузки для линейной смешанной модели с R boot
пакетом.
Код R
Вот мой 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
. Любая идея, где я могу получить этот пакет. Любая помощь будет высоко оценена. Заранее спасибо за вашу помощь и время.
источник
R
.Возможно, вы захотите взглянуть на
bootMer
функцию в версии для разработчиков lme4,который может делать основанную на модели (полу) параметрическую загрузку смешанных моделей ... Просто проверьте
?bootMer
источник