ANOVA
到目前為止，你已經建立了輸出物件，藉此儲存並擷取元素，以便客製化並更詳細地呈現結果。不過，有些輸出物件（例如用 aov() 執行變異數分析（ANOVA）所得到的物件）也能作為其他函式的輸入，例如 summary() 與 plot()，以產生客製化的結果報告與視覺化。此外，aov() 的輸出也能作為 TukeyHSD() 的輸入，用來進行事後的成對比較檢定。
已為你載入 abaloneKeep 資料集與 dplyr 套件。
本練習屬於課程
給 SAS 使用者的 R
練習說明
- 計算依性別分組的
age的n、平均數、標準差與變異數。 - 以
sex分組對age進行 ANOVA，並將輸出存為abaov。 - 執行
abaov的summary()。 - 對
abaov進行TukeyHSD()事後成對比較檢定。
動手互動練習
試著完成這個範例程式碼，體驗一下這個練習。
# Compute n(), mean(), sd() and var() of age by sex
abaloneKeep %>% group_by(___) %>% select(___, ___) %>%
summarise(across(everything(), list(mean = ___,
___ = ~sd(.x),
var = ___)),
N = ___)
# Run aov() of age by sex, save as abaov
abaov <- aov(___ ~ ___, data = ___)
# Run summary() of abaov
___
# Perform TukeyHSD posthoc pairwise tests on abaov
TukeyHSD(___)