r - 在 R 中定义嵌套设计

Question

我有三个因素A，B，C。B 和 C 是随机的。B与A嵌套。C和B，C和A之间存在交互。我尝试使用lme函数，但我不知道如何定义嵌套以及C和B的交互。

lme(dades~A*C,random=list(~1|B, ~1|C),data=mydata)

score 1 · Accepted Answer

我将尽力在这里提供帮助，但如果其他人有更好的答案，当然请听从他们的建议。

首先，从构建无条件模型开始——一个没有任何预测变量或自变量的模型。在您的情况下，看起来 A 代表 2 级组。这将是使用您似乎正在使用的 nlme 包的无条件模型：

model1 <- lme(fixed = dades~1, random =~1|A)

然后，您可以对 B 和 C 的影响进行建模（允许截距而不是 B 和 C 的斜率在级别 2 是随机的）：

model2 <- lme(fixed = dades~B+C, random =~1|A)

您可以允许斜率在 A 上是随机的：

model3 <- lme(fixed = dades~B+C, random = ~B+C|A)

要对交互进行建模，您可以将预测变量相乘（就像您所做的那样）：

model4 <- lme(fixed = dades~B+C+B*C, random = ~1|A)

我不确定如何在 C 和 A 之间添加交互，因为 A 不是预测变量，而是 2 级组。

同样，如果另一个答案更有帮助，请使用它。

1 回答 1