在 R 中定义嵌套设计答案

【问题标题】：define nested design in R在 R 中定义嵌套设计
【发布时间】：2014-12-14 22:38:50
【问题描述】：

我有三个因素 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)

【问题讨论】：

【解决方案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 级组。

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

【讨论】：