展示沙眼衣原体结果
在上一个练习中，您为伊利诺伊州的沙眼衣原体数据拟合了一个 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()