BIBD——猫的肾功能
为确保我们真正理解 BIBD 的样子，让我们从零开始构建一个包含 BIBD 的数据集。
假设我们想测试 4 种不同湿粮对猫肾功能的影响。不过猫粮很贵，所以为了省钱，我们每个区组只测试 3 种食物。区组因素为猫的 color，因为我们对颜色本身并不感兴趣。结果变量是血肌酐水平，这是猫和人类肾功能与功能障碍的常用指标。
本练习是课程的一部分
R 中的实验设计
练习说明
- 已为您加载自定义函数
lambda()。使用t= 4、k= 3、r= 3 计算lambda，以确保 BIBD 可行。 - 运行代码以组装数据集。您可以看到每个区组中使用食物处理的顺序。
- 根据上面实验描述，使用
aov()创建cat_model，并用summary()查看结果。不同类型的湿粮是否会影响肌酐水平？
交互式实操练习
通过完成这段示例代码来试试这个练习。
# Calculate lambda
___
# Build the data.frame
creatinine <- c(1.98, 1.97, 2.35, 2.09, 1.87, 1.95, 2.08, 2.01, 1.84, 2.06, 1.97, 2.22)
food <- as.factor(c("A", "C", "D", "A", "B", "C", "B", "C", "D", "A", "B", "D"))
color <- as.factor(rep(c("Black", "White", "Orange", "Spotted"), each = 3))
cat_experiment <- as.data.frame(cbind(creatinine, food, color))
# Create cat_model and examine with summary()
___ <- ___ (___ ~ ___ + color, data = ___)
___