这可能是一个"愚蠢的"问题,但我需要从 merMod 对象( lme4 包)中获取潜在特质分数。另外,我似乎没有找到merMod对象中值的任何解释。对我来说,知道mu,wtres,eta,u,LUtx,Utx,Utr,V和Xwts应该是什么(通常)是有帮助的。
我猜(ed)潜在特质分数不在对象中,但需要计算。 Irtoys提供了一个功能(例如,dpv),但要求我有一个响应矩阵,只能是0或1而没有NA。我的数据包括NA设计。此外,我将使用值不是0和1的模型。有关解决方法的任何想法吗? Irtoys包装似乎提供了我所需要的一切,但如果我因为这些限制而无法使用它,那将是一个遗憾。
提前谢谢你, KH
答案 0 :(得分:2)
我将继续推广我的评论以回答。
可能存在关于“潜在特质分数”的确切含义的术语问题,但我99%肯定您需要ranef()
:请参阅?ranef.merMod
。 (在线性混合模型世界中,这些被称为BLUP; Doug Bates更喜欢将它们称为条件模式,以便术语扩展到GLMM(它们不再是最佳,线性或无偏的)。)
我很遗憾地给你一个链接而不是一个明确的答案,但是找到merMod
对象内容的最佳来源可能是this ArXiv preprint,正在按< em>统计软件期刊(希望,现在任何一天)。 ?getME
也可能有用 - 它的优势在于,您在那里找到的任何内容都可以安全使用,而不必担心merMod
对象的内容在将来的版本中会发生变化。 (使用@
- 访问者后,所有投注均已关闭。)