R單向方差分析（僅獲得1個DF並期望2個DF）

Question

我正在研究UCLA網站http://www.ats.ucla.edu/stat/r/faq/posthoc.htm上的單向方差分析。 當我運行命令a1 <-aov(write ~ ses) ，我的輸出與示例輸出不同。 當我運行命令summary(a1) ，我特別煩惱，我在ses上的DF是1，並且有三個ses類別（1,2,3），所以我期望有2個DF，這就是示例在網站上顯示。 我已經檢查了“寫”列和“ ses”列的數據，其計數和平均值似乎與示例匹配，但是aov(write ~ ses)的結果aov(write ~ ses)匹配。 有什么改變嗎？ 為什么我只得到1 DF。

hsb2 <- read.table("http://www.ats.ucla.edu/stat/data/hsb2.csv", sep=",", header=TRUE)
a1 <- aov(write ~ ses, data = hsb2)
summary(a1)
#              Df Sum Sq Mean Sq F value Pr(>F)   
# ses           1    770   769.8   8.908 0.0032 **
# Residuals   198  17109    86.4         

Answer 1

您正在學習的頁面有錯誤，因為它沒有告訴您如何正確輸入數據。 ses變量應該是一個因素，正如我們從它們提供給我們的數據中看到的那樣，它被讀為數字形式：

str(hsb2$ses)

如果將其轉換為因子，則得到與示例相同的答案：

hsb2$ses <- as.factor(hsb2$ses)
a1 <- aov(write ~ ses, data=hsb2)
summary(a1)

             Df Sum Sq Mean Sq F value  Pr(>F)   
ses           2    859   429.4    4.97 0.00784 **
Residuals   197  17020    86.4                   
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

此外，大多數R用戶強烈建議不要使用attach。