如何在 r 中使用 lme4 获得组内关联?

How do I get the within-group association using lme4 in r?

设置: 我正在测试成对个体之间的特征 (BMI) 关联是否随时间变化。我有重复测量,一对中的每个人在 7 个时间点提供 BMI 数据。下面是一个长格式的简化数据框,其中包含 Pair ID(赋予每对个体的标识符)、两个个体在每个时间点(BMI_1 和 BMI_2)的 BMI 测量值,以及一个时间具有七个间隔的变量,编码为连续。

Pair_ID BMI_1 BMI_2 Time
1 25 22 1
1 23 24 2
1 22 31 3
1 20 27 4
1 30 26 5
1 31 21 6
1 19 18 7
2 21 17 1
2 22 27 2
2 24 22 3
2 25 20 4

首先,我主要感兴趣的是测试within-pair association(下面BMI_2的回归系数)以及它是否随时间变化(BMI_2和Time之间的交互) .我想排除任何对间效应,因此我只测试成对内随时间推移的相关性。

我正计划拟合以下形式的线性混合模型:

    lmer(BMI_1 ~ BMI_2 * Time + (BMI_2 | Pair_ID), Data)

我了解模型的参数(例如,随机 slopes/intercepts),并且 BMI_2 * 时间交互测试是否 BMI_1 和 BMI_2 之间的关系由时间主持。

但是,我不确定如何识别(均值)对内回归系数,以及我的方法是否适用于此。

其次,我有兴趣了解 BMI_2 * 时间交互中的对之间是否存在变化(即对之间斜率的差异) - 例如,[=37 之间的关联是否存在=] 和 BMI_2 在某些对中会随着时间的推移而增加,而在其他对中则不会?

为此,我正在考虑拟合这样的模型:

    lmer(BMI_1 ~ BMI_2 * Time + (BMI_2 : Time | Pair_ID), Data)

然后查看 BMI_2 中的方差:时间随机效应。据我了解,大方差意味着这种相互作用效应在两对之间变化很大。

如能就这些问题(尤其是第一个问题)提供任何帮助,我们将不胜感激。

P.s.,如果问题格式不正确,请见谅。这是我的第一次尝试。

回答的完整性。 @benimwolfspelz 的评论是正确的。这在应用工作的某些领域被称为“背景效应”。这个想法是通过对每个组进行均值居中并拟合以均值为中心的变量(这将估计组内组件)组均值(这将估计组件之间)。