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

How do I obtain latent trait scores from a merMod object in R?

这可能是一个 "stupid" 问题,但我需要从 merMod 对象( lme4 包)。另外,我似乎没有在 merMod 对象中找到任何值的解释。了解 mu、wtres、eta、u、LUtx、Utx、Utr、V 和 Xwts 应该是什么(通常)对我很有帮助。

我猜(编辑)潜在特征分数不在对象中,但需要计算。 Irtoys 提供了执行此操作的功能(例如 dpv),但要求我有一个只能为 0 或 1 且没有 NA 的响应矩阵。不过,我的数据在设计上包括了 NA。此外,我将使用具有非 0 和 1 值的模型。关于变通方法有什么想法吗? Irtoys 包似乎提供了我所需要的一切,但如果因为这些限制而无法使用它,那将是一个遗憾。

提前谢谢你, KH

我将继续将我的评论提升为答案。

  • 关于 "latent trait scores" 的确切含义可能存在术语问题,但我 99% 确定您想要 ranef():请参阅 ?ranef.merMod . (在线性混合模型世界中,这些被称为 BLUP;Doug Bates 更喜欢称它们为 条件模式,以便该术语扩展到 GLMM(它们不再一定是最佳的、线性的或无偏的).)

  • 很抱歉给你一个 link 而不是一个明确的答案,但是了解 merMod 对象内脏的最佳来源可能是 this ArXiv preprint,在 Journal of Statistical Software 上出版(希望现在任何一天都能出版)。 ?getME 也可能有用——它的优点是您可以安全地使用在那里找到的任何东西,而不必担心 merMod 对象的内容会在未来的版本中发生变化。 (一旦你使用 @-accessor,所有的赌注都关闭了。)