顯示披衣菌結果
在前一個練習中，你已對伊利諾州的披衣菌（chlamydia）資料擬合了 GLMER。這個練習要帶你走過一些用來呈現結果的方法。你可以用這些方法為客戶或你正在撰寫的文件產生模型摘要來說明你的結果。不過，我也鼓勵你學會自行操作並探索模型輸出，發展你自己的結果呈現方式。發展獨特的方法能讓你在資料科學領域更顯眼！
你將會做以下事情：
- 檢視模型估計值。
- 繪製資料，並對每個年齡階層擬合一個
glm。雖然與glmer()的輸出不完全相同，但這種近似方法有助於用更直覺、好理解的視覺方式來呈現結果。
本練習屬於課程
R 的階層式與混合效果模型
練習說明
- 使用
fixef()從model_out萃取固定效果估計值。 - 使用
ranef()從model_out萃取隨機效果估計值。 - 執行程式碼，使用
ggplot2方法繪製資料。
動手互動練習
試著完成這個範例程式碼，體驗一下這個練習。
# Extract out fixed effects
___
# Extract out random effects
___
# Run code to see one method for plotting the data
ggplot(data = il_data_2,
aes(x = year, y = count, group = county)) +
geom_line() +
facet_grid(age ~ . ) +
stat_smooth(method = "glm",
method.args = list(family = "poisson"),
se = FALSE,
alpha = 0.5) +
theme_minimal()