阻隔(Blocking)
雖然事實並非如此,但先假設補充品類型其實是我們想要透過阻隔來控制的干擾因子,而我們真正感興趣的是維生素 C 劑量對天竺鼠牙齒生長的影響。
如果以補充品類型做為阻隔,我們會建立更相似的群組,也就是群組內有相同的補充品類型,讓我們能只檢視劑量對牙齒長度的影響。
我們將使用 aov() 函式來檢視這件事。aov() 透過呼叫 lm() 並以 anova() 檢視結果,在一次函式呼叫中建立線性迴歸模型。要使用 aov(),我們一樣需要公式表示法,就像隨機化練習那樣,但這次的公式應為 len ~ dose + supp,表示我們已依補充品類型做了阻隔。(我們會在下一章更詳細介紹 aov() 與 anova()。)
ggplot2 已為你載入。
本練習屬於課程
R 的實驗設計
練習說明
- 繪製盒狀圖,視覺化檢視牙齒長度是否因
dose不同而有差異。dose已替你轉為類別變數。 - 使用
aov()檢測dose與supp對len的影響。將結果儲存為名為ToothGrowth_aov的模型物件。 - 使用
summary()檢視ToothGrowth_aov,判斷劑量是否對牙齒長度有顯著影響。
動手互動練習
試著完成這個範例程式碼,體驗一下這個練習。
# Create a boxplot with geom_boxplot()
ggplot(___, aes(x = ___, y = ___)) +
___()
# Create ToothGrowth_aov
___ <- aov(___, data = ___)
# Examine ToothGrowth_aov with summary()
___