Les greffes suivantes sont extraites de cet article . Je suis novice dans le bootstrap et j'essaie d'implémenter le bootstrap paramétrique, semi-paramétrique et non paramétrique pour le modèle mixte linéaire avec le R boot
package.
Code R
Voici mon R
code:
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
Des questions
- Comment faire du bootstrap paramétrique, semi-paramétrique et non paramétrique pour des modèles mixtes avec
boot
package? - Je suppose que je fais du bootstrap non paramétrique pour un modèle mixte dans mon code.
J'ai trouvé ces diapositives mais n'a pas pu obtenir le package R merBoot
. Toute idée où je peux obtenir ce package. Toute aide sera grandement appréciée. Merci d'avance pour votre aide et votre temps.
la source
R
.Vous voudrez peut-être jeter un œil à la
bootMer
fonction dans la version de développement de lme4,qui peut faire un bootstrap (semi-) paramétrique basé sur un modèle de modèles mixtes ...
?bootMer
la source