模型後驗證圖與變異數
在上一個練習中,我們發現 int_rate 會隨 grade 不同而改變。現在要驗證這個模型。對於線性迴歸,這表示要檢視 Residuals vs. Fitted 與 Normal Q-Q 圖。
如果你在 R 中對模型物件呼叫 plot(),它會自動繪出上述兩個圖,再加上另外兩個。你會解讀這些圖來評估模型配適度。我們已在影片中說明如何進行。
ANOVA 與線性建模的另一個假設是變異數同質性。所謂同質性就是「相同」,在此表示各個 grade 水準下的 int_rate 變異數應該相同。我們可以使用 bartlett.test() 來檢定變異數同質性,此函式以公式與資料集作為輸入。
本練習屬於課程
R 的實驗設計
練習說明
- 先執行含有
par()的第一行程式碼,讓圖以 2x2 方格輸出。 - 對
grade_aov(已為你建立)呼叫plot(),產生模型診斷圖。 - 使用
bartlett.test()檢定變異數同質性。
動手互動練習
試著完成這個範例程式碼,體驗一下這個練習。
# For a 2x2 grid of plots:
par(mfrow=c(___, ___))
# Plot grade_aov
___
# Bartlett's test for homogeneity of variance
___