毒理學資料
一位毒理學家想知道,隨著劑量增加,接觸某種化學物是否會提高死亡率。她讓測試生物接觸不同劑量,並記錄每個水槽中死亡的數量。每個劑量在 3 個不同的日子各重複一次。
首先,建立一個 GLM。先從簡單模型開始有助於除錯。之後你也會把結果與 GLMER 比較。這個模型將 dose 和 replicate 視為固定效果(fixed effects)。
接著,建立一個 GLMER。這次把 dose 視為固定效果,而把 replicate 視為隨機截距(random-effects intercept)。
最後,使用各模型的 coef() 檢視係數估計。注意兩個模型的截距估計顯示方式不同。
本練習屬於課程
R 的階層式與混合效果模型
動手互動練習
試著完成這個範例程式碼,體驗一下這個練習。
# Fit glm_out and look at its coefficient estimates
glm_out <- glm(___ ~ ___ + ___,
family = "___", data = df)
coef(___)