1 Way Anova Test in R不正確df

Question

我想在R中運行anova測試以確定每個品牌的週數是否有差異。我的代碼如下：

weeks = c(100,96,92,96,92,76,80,75,84,82,108,100,96,98,100) 
brand = c('a','a','a','a','a','b','b','b','b','b','c','c','c','c','c') 
battery.exp = data.frame(brand, weeks) 
test.1 = aov(battery.exp$weeks~battery.exp$brand) 
summary(test.1) 

我在摘要中得到不正確的結果。品牌的df應該是2，但是我的輸出是3.有人知道爲什麼我會在我的anova測試中得到不正確的結果嗎？我也嘗試用下列替換倒數第二行：

test.1 = aov(battery.exp$weeks~factor(battery.exp$brand)) 

但這也沒有解決問題。有關如何解決的任何建議？

Answer 1

默認情況下，使用帶截距的模型，這意味着您需要爲分類變量使用參考等級，因此您只有兩個自由參數。如果你不希望使用攔截，你可以明確地請求與

summary(aov(weeks~brand-1, battery.exp)) 

（注意公式中的-1）