NHANES RCBD
回顾一下,我们的分区实验包含一种处理:医生要求受试者在饮食中减少脂肪或热量,我们要检测这对体重(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 = ___(___))