在混合模型中嵌套术语

机器算法验证 混合模式 术语 嵌套数据 lme4-nlme
2022-04-05 09:42:39

假设我有以下嵌套lmer结构:

lmer(Y ~ X1 + X2 + X1:X2 + (1 | A) + (1 | A:B), data=d)

这与以下内容相同:

lmer(Y ~ X1 * X2 + (1 | A/B), data=d)

现在,如果我将其写在报告中,我会说:

Y由和的交互建模X1并包括(作为固定效应)并嵌套在其中(作为随机效应)。X2X1X2BA

+ (1 | A/B)+ (1 | A) + (1 | A:B)象征嵌套 和 是等价的。但是当我不想A估计的效果时如何:

lmer(Y ~ X1 + X2 + X1:X2 + (1 | A:B), data=d)

我还能简单地说:

1. B嵌套在A(即使A现在省略了随机效果?)

或者应该是:

2. B嵌套在A没有效果A

或者由于:固定效果部分中的 表示交互,我可以说:

3.模型的随机效应部分作为 和 的交互作用A给出B

1个回答

我会将最终模型描述为简单地具有

每个观察到的 A 和 B 组合的随机截距。

虽然它看起来像一个交互,因为你只是使用 R 的交互语法,但你实际上只是以一种比以前的规范稍微不那么结构化的方式重新定义了一个分组变量,确保随机效应现在可以比在之前的安排。

出于说明性目的,您可能希望为“每个观察到的 A 和 B 组合”起一个对您正在分析的研究而言对读者有意义的名称。