我如何在 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交互和順序隨機效應