클라미디아 결과 표시하기
이전 연습 문제에서 일리노이주의 클라미디아 데이터를 대상으로 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()