r-使用MItools从估算数据集中提取2级随机效应

Question

我正在尝试提取我使用lmer运行的三级分层回归的三级随机效应。问题是，我的数据集是估算的(合理值(，我试图只提取L3随机效应来与MItools合并。我可以在只有一个PV 的常规lmer模型上提取它们

m0 <- lmer(PV4MATH ~ 1 + (1 | SCHOOLID) + (1 | CNT), data = pisa, REML = F)

ranef(m0($CNT

但是，当我试图从一个模型中准确提取这些时，我失败了，我在所有PV上都使用了mitools，然后试图将它们与MIcombine 组合

model0<-with(mis, lmer(MATH ~ 1 + (1 | SCHOOLID) + (1 | CNT), REML = F))

我只想提取L2随机效应(=所有国家的平均值，同时控制单个水平的变量，我知道这是一个空模型(。这就是我提取所有ranefs的方法，但我对MIcombine有问题。

betas<-MIextract(model0, , fun = ranef)

我尝试了$CNT的各种组合，但都不起作用。有人有主意吗？我收到一个错误，告诉我这是一个二进制运算符的非数字运算(不确定这是否是正确的翻译(。我对固定效果没有问题。

Answer 1

找到了一个解决方案。。。万一有人也需要这样做：

betas2 <- list()
for ( i in 1:5){ betas2[[i]] <- betas[[i]]$CNT 

}

提取列表"中的所有系数；手动"；然后你就可以使用它们了。也适用于差异。