Los siguientes injertos se toman de este artículo . Soy novato en bootstrap e intento implementar el bootstrapping paramétrico, semiparamétrico y no paramétrico para el modelo mixto lineal con R boot
paquete.
Código R
Aquí está mi R
código:
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
Preguntas
- ¿Cómo hacer bootstrapping paramétrico, semiparamétrico y no paramétrico para modelos mixtos con
boot
paquete? - Supongo que estoy haciendo bootstrapping no paramétrico para el modelo mixto en mi código.
He encontrado estas diapositivas , pero no pude conseguir el paquete R merBoot
. Alguna idea de dónde puedo obtener este paquete. Cualquier ayuda será muy apreciada. Gracias de antemano por su ayuda y tiempo.
R
.