在R中复制权重、lmer和PISA数据

Question

我试图以最“准确”的方式理解和使用PISA数据。

最后，我需要的是国家方法(但我想控制个人水平的变量，所以我有国家净效应，因为我想分析教育系统的特征)。我目前正在使用ranef从lmer多级模型中获取它们(使用(1|CNT/SCHOOLID)和(1 |CNT)嵌套)。

我首先想到使用调查包，所以我可以包括复制权重(svyrepdesign工作)，但然后我无法控制单个级别的变量(没有多级可能).I尝试svynlm -没有工作，它确实花费了太长的时间(48小时，然后不收敛)。我现在明白了，我应该使用参议院权重，而不是学生的最终权重(到目前为止，这不是问题，我可以将其包括在lmer中)。

如果我使用mitools和lmer，我有一个多层模型(具有参议院权重)和合理的值，但不是repweight。由于我没有在我的进一步分析中包括se，我可以只使用参议院权重，但我需要一个好的解释，而不是“技术考虑”。

有没有人知道走出困境的方法？

Answer 1

看看BIFIEsurvey吧。对于两级回归，BIFIE.twolevelreg允许使用PV和复制权重。

还要考虑如何衡量你的体重:请看这里: Rabe-Hesketh，S.，& Skrondal，A。(2006)。复杂测量数据的多级建模。皇家统计学会杂志。A系列(一般)，169,805-827。