R中的N - 方差分析

Question

我正在从拉丁方块实验中进行ANOVA测试。出于某种原因，aov（）函数仅识别所使用的3个变量中的2个。

batch = c(1,2,3,4,5,1,2,3,4,5,1,2,3,4,5,1,2,3,4,5,1,2,3,4,5)
day = c(1,1,1,1,1,2,2,2,2,2,3,3,3,3,3,4,4,4,4,4,5,5,5,5,5)
ingredient = c('A','B','C','D','E','A','B','C','D','E','A','B','C','D','E','A','B','C','D','E','A','B','C','D','E')
rxn.time = c(8,11,4,6,4,7,2,9,8,2,1,7,10,6,3,7,3,1,6,8,3,8,5,10,8)
rxn.exp = data.frame(batch, day, ingredient, rxn.time)
test.10 = aov(rxn.exp$rxn.time~as.factor(rxn.exp$batch)+as.factor(rxn.exp$day)+as.factor(ingredient))
summary(test.10)

摘要仅为前两个因子生成输出。我尝试使用factor（）函数而不是as.factor以及其他几种尝试让R根据需要运行3个不同的变异源加上任何噪声。有人可以解释为什么R不合作以及如何解决？

Answer 1

这只是因为“批次”和“成分”的组合是相同的（如果批次= 1，成分= A等）。如果你使用随机成分，它的工作原理。 E.g：

UIView

另请注意，我引用了lm函数中的数据帧。