我如何在 R 或 glm 和 ANOVA 中做嵌套模型？

Question

我一直在努力解决这个问题，但我目前正在尝试分析我的数据，到目前为止我有这段代码：

data2 <- read.table("B.txt",header=T)

attach(data2)

Sex <- as.factor(Sex)

Class <- as.factor(Class)

Order <- as.factor(Order)

Envi <- as.factor(Environment)

model1 <- glm(LS~Sex*Class*Envi)

model2 <- glm(LS~Sex*Class)

model3 <- glm(LS~Sex)

model4 <- glm(LS~Class)

model5 <- glm(LS~Envi)

但是有人告诉我，我必须将 Order 嵌套在类中，因为对于我的数据集，我的自由度将是我的类，并且交互测试可能不是最好的方法。 该类是一个固定效果，顺序是随机的。 我将如何在课堂上嵌套顺序，以便我可以针对课堂上的性测试寿命？ 有没有更好的方法来测试这个？

Answer 1

我从问题的描述中了解到order是嵌套在Class中的随机效应。 混合模型可以从几个包中实现，例如：

library(lme4)
glmer(LS~Sex*Class +(1|Order:Class), data=data2)

我们在同一模型中估计Sex:Class交互和顺序随机效应