1-Wege-Anova-Test in R inkorrekt df

Question

Ich möchte einen Anova-Test in R durchführen, um herauszufinden, ob es einen Unterschied in der Anzahl der Wochen jeder Marke gibt. Mein Code ist wie folgt:

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) 

Ich bekomme falsche Ergebnisse in der Zusammenfassung. Das df der Marke sollte 2 sein, aber meine Ausgabe sagt 3. Weiß jemand, warum ich falsche Ergebnisse in meinem Anova-Test bekommen würde? Ich habe auch versucht, die vorletzte Zeile durch folgende zu ersetzen:

aber das hat auch nicht das Problem gelöst. Irgendwelche Vorschläge zur Behebung?

Answer 1

Standardmäßig wird ein Modell mit einem Achsenabschnitt verwendet. Dies bedeutet, dass Sie einen Referenzpegel für Ihre kategoriale Variable verwenden müssen, damit Sie nur zwei freie Parameter haben. Wenn Sie möchten, einen Schnittpunkt nicht verwenden, können Sie verlangen ausdrücklich einen Schnittpunkt frei Modell mit

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

(beachten Sie die -1 in der Formel)