1 Простой тест Anova in R неправильный df

Question

Я хочу запустить тест anova в R, чтобы узнать, есть ли разница в количестве недель каждой марки. Мой код выглядит следующим образом:

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 в формуле)