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 = ___)
___