提取迴歸係數
我們常常想知道模型估計出的迴歸係數為何。雖然 summary() 函式會提供模型輸出,但有時你也會想直接存取這些輸出。
固定效果的估計值可以用 fixef() 函式直接取得。隨機效果的估計值可以用 ranef() 函式直接取得。我們也可以用 confint() 函式為固定效果提取信賴區間。
broom.mixed 套件也提供用於從 lmer() 模型提取結果的 tidy 方法,也就是 tidy() 函式。不過,由於混合效果模型本身的複雜度,這些結果比起一般的 tidy 輸出更複雜、也沒那麼「整齊」。
已經幫你配適好一個 lmer 模型,並儲存為 out。
本練習屬於課程
R 的階層式與混合效果模型
練習說明
- 使用已儲存的模型
out,用fixef()提取固定效果係數。 - 使用已儲存的模型
out,用ranef()提取隨機效果係數。 - 使用已儲存的模型
out,用confint()估計 95% 信賴區間。 - 使用
tidy()配合out與conf.int = TRUE,以一個「tidy」指令重複前面三個呼叫。
動手互動練習
試著完成這個範例程式碼,體驗一下這個練習。
# Extract the fixed-effect coefficients
___
# Extract the random-effect coefficients
___
# Estimate the confidence intervals
___
# Use the tidy function
___