使用 agricolae 绘制 BIBD
我们也可以使用 agricolae 来绘制 BIBD。design.bib() 至少需要输入处理(treats)、一个表示区组水平数量的整数 k,以及一个随机 seed。
BIBD 的一个主要特点是:在输出的每个区组(列)中,不一定会包含所有处理。
不过,正如您在视频和上一个练习中所见,BIBD 有时并不有效,您需要做一些计算来确认设计是否可行。若设计不合法,design.bib() 会返回错误消息提醒您。
让我们用 agricolae 画几个 BIBD,以查看该包给出的不同警告与错误信息。
本练习是课程的一部分
R 中的实验设计
交互式实操练习
通过完成这段示例代码来试试这个练习。
# Create my_design_bibd_1
___ <- design.bib(LETTERS[1:3], k = ___, seed = ___)