r - 如何从 R lme 对象中提取特定组的随机效应？

Question

我已经成功地拟合了一个线性混合效应模型，我正在寻找为各个组提取随机效应分量。我知道可以使用提取随机效应的完整列表

random.effects(model)

然后 print(random.effects(model)) 给出了一个两列的组名和随机效应列表，即使数据本身似乎只有一列。我的问题是是否可以通过组名“查找”与特定组关联的随机效应，或者，如果没有，我如何找到与数据中随机效应相同顺序的组名列表由 random.effect() 输出的帧。

谢谢Mark Ch。

Answer 1

事实证明，问题在于我对组进行索引的特殊方式。 ranef(lme) 返回一个数据框，其中行名称是组名称。在我的数据中，我使用了一个很长的数字来区分组，R 将其四舍五入到少数小数位。这意味着不可能按名称准确引用单个组。

我通过将每个索引转换为 base-62 数字来解决这个问题。我使用数字和大小写字母作为数字中的字符集。 (也就是说，数字匹配 [a-zA-Z0-9]*)这大大减少了组名的长度，并使 R 无法对组名进行四舍五入——使用的字符越多，它变得越短.

现在，如果我这样做:

M3.ranef <- ranef(M3)
x <- M3.ranef[group_ID,1]

x 是名为 group_ID 的组的随机效应，这是数据帧应该如何工作。