aov() 的基本用法

方差分析(aov)用於確定兩個或更多個組的平均值是否彼此顯著不同。假設響應彼此獨立,正態分佈(在每個組內),並且假設組內差異相等。

為了完成分析,資料必須是長格式(請參閱重塑資料主題)。aov()lm() 函式的包裝器,使用 Wilkinson-Rogers 公式表示 y~f,其中 y 是響應(獨立)變數,f 是表示組成員資格的因子(分類)變數。如果 f 是數字而不是因子變數,則 aov() 將以 ANOVA 格式報告線性迴歸的結果,這可能會讓沒有經驗的使用者感到驚訝。

aov() 函式使用型別 I(順序)平方和。這種型別的平方和按順序測試所有(主要和互動)效果。結果是,測試的第一個效果也被賦予它與模型中的其他效果之間的共享方差。為了使這種模型的結果可靠,應平衡資料(所有組的大小相同)。

當 I 型平方和的假設不成立時,可以適用 II 型或 III 型平方和。型別 II 平方和在每個其他主效應之後測試每個主效應,從而控制任何重疊方差。然而,型別 II 平方和假設主效應之間沒有相互作用。

最後,Type III of Squares 測試每個主要效果每次互動之後的每個主要效果。當存在互動時,這使得 III 型平方和成為必需。

II 型和 III 型平方和在 Anova() 函式中實現。

mtcars 資料集為例。

mtCarsAnovaModel <- aov(wt ~ factor(cyl), data=mtcars)

要檢視 ANOVA 模型的摘要:

summary(mtCarsAnovaModel)

還可以提取基礎 lm() 模型的係數:

coefficients(mtCarsAnovaModel)