Cangkok berikut diambil dari artikel ini . Saya pemula untuk bootstrap dan mencoba mengimplementasikan bootstrap parametrik, semiparametrik, dan nonparametrik untuk model campuran linier dengan R boot
paket.
Kode R
Ini R
kode saya :
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
Pertanyaan
- Bagaimana melakukan bootstrap parametrik, semiparametrik, dan nonparametrik untuk model campuran dengan
boot
paket? - Saya kira saya sedang melakukan bootstrap nonparametrik untuk model campuran dalam kode saya.
Saya menemukan slide ini tapi tidak bisa mendapatkan paket R merBoot
. Tahu di mana saya bisa mendapatkan paket ini. Bantuan apa pun akan sangat dihargai. Terima kasih sebelumnya atas bantuan dan waktu Anda.
sumber
R
.Anda mungkin ingin melihat
bootMer
fungsi dalam versi pengembangan lme4,yang dapat melakukan bootstrap parametrik berbasis model (semi-) dari model campuran ... Cukup periksa
?bootMer
sumber