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)

Answer 1

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

首先，从构建无条件模型开始 - 没有任何预测变量或自变量的模型。在您的情况下，它看起来像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级组。

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