Os seguintes enxertos são retirados deste artigo . Eu sou novato no bootstrap e estou tentando implementar o bootstrap paramétrico, semiparamétrico e não paramétrico para o modelo misto linear com o R boot
pacote.
Código R
Aqui está o meu 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
Questões
- Como executar bootstrap paramétrico, semiparamétrico e não paramétrico para modelos mistos com
boot
pacote? - Eu acho que estou fazendo bootstrapping não paramétrico para modelo misto no meu código.
Eu encontrei esses slides , mas não conseguiu obter o pacote R merBoot
. Alguma idéia de onde posso obter esse pacote. Qualquer ajuda será muito apreciada. Agradecemos antecipadamente a sua ajuda e tempo.
fonte
R
.Você pode dar uma olhada na
bootMer
função na versão de desenvolvimento do lme4,que pode executar bootstrap (semi-) paramétrico baseado em modelo de modelos mistos ... Basta verificar
?bootMer
fonte