NHANES 的 RCBD
回想一下,我們的分區(blocked)實驗包含一個處理:醫師請病患在飲食中降低脂肪或熱量,而我們要檢驗這對體重(bmxwt)的影響。我們計畫以性別分區,NHANES 資料集中性別儲存在 riagendr。請記得,分區的目的在於建立彼此盡可能相似的實驗組。以性別進行分區表示:若我們觀察到處理對 bmxwt 有影響,則更可能是處理本身造成的效果,而不是個體的性別所致。
在 R 程式碼中,你可以在適當的建模函式裡,使用類似 outcome ~ treatment + blocking_factor 的公式來表示分區實驗。
nhanes_final 已可使用。
本練習屬於課程
R 的實驗設計
練習說明
- 使用
aov()建立nhanes_rcbd。回想處理變數儲存在mcq365d,你要檢驗的結果是bmxwt,分區因子是riagendr。 - 用
summary()檢視nhanes_rcbd的結果。 - 使用
dplyr的函式,依mcq365d與riagendr檢視平均體重。
動手互動練習
試著完成這個範例程式碼,體驗一下這個練習。
# Use aov() to create nhanes_rcbd
___
# Check results of nhanes_rcbd with summary()
___
# Print mean weights by mcq365d and riagendr
___ %>%
group_by(___, ___) %>%
summarize(mean_wt = ___(___))