如何从R中的merMod对象获取潜在特质分数?

时间:2015-02-06 14:23:02

标签: r lme4

这可能是一个"愚蠢的"问题,但我需要从 merMod 对象( lme4 包)中获取潜在特质分数。另外,我似乎没有找到merMod对象中值的任何解释。对我来说,知道mu,wtres,eta,u,LUtx,Utx,Utr,V和Xwts应该是什么(通常)是有帮助的。

我猜(ed)潜在特质分数不在对象中,但需要计算。 Irtoys提供了一个功能(例如,dpv),但要求我有一个响应矩阵,只能是0或1而没有NA。我的数据包括NA设计。此外,我将使用值不是0和1的模型。有关解决方法的任何想法吗? Irtoys包装似乎提供了我所需要的一切,但如果我因为这些限制而无法使用它,那将是一个遗憾。

提前谢谢你, KH

1 个答案:

答案 0 :(得分:2)

我将继续推广我的评论以回答。

  • 可能存在关于“潜在特质分数”的确切含义的术语问题,但我99%肯定您需要ranef():请参阅?ranef.merMod。 (在线性混合模型世界中,这些被称为BLUP; Doug Bates更喜欢将它们称为条件模式,以便术语扩展到GLMM(它们不再是最佳,线性或无偏的)。)

  • 我很遗憾地给你一个链接而不是一个明确的答案,但是找到merMod对象内容的最佳来源可能是this ArXiv preprint,正在按< em>统计软件期刊(希望,现在任何一天)。 ?getME也可能有用 - 它的优势在于,您在那里找到的任何内容都可以安全使用,而不必担心merMod对象的内容在将来的版本中会发生变化。 (使用@ - 访问者后,所有投注均已关闭。)