使用 agricolae 繪製 BIBD
我們也可以用 agricolae 來繪製 BIBD。design.bib() 至少需要提供 treatments(treats)、代表區組層級數量的整數 k,以及 seed 作為輸入。
BIBD 的重點是:在輸出結果的每個區組(欄)中,並不會使用到所有的處理(treatment)。
不過,從影片與上一題你已經知道,有時候 BIBD 其實無法成立,你需要做一點數學驗證來確認設計是否可行。若設計不合法,design.bib() 會回傳錯誤訊息提醒你。
我們用 agricolae 畫幾個 BIBD,來看看這個套件可能給出的不同警告與錯誤訊息。
本練習屬於課程
R 的實驗設計
動手互動練習
試著完成這個範例程式碼,體驗一下這個練習。
# Create my_design_bibd_1
___ <- design.bib(LETTERS[1:3], k = ___, seed = ___)