Hiển thị kết quả bệnh chlamydia
Trong bài trước, bạn đã khớp một GLMER cho dữ liệu chlamydia ở Illinois. Ở bài này, chúng ta sẽ điểm qua một số cách hiển thị kết quả. Bạn có thể dùng các cách này để tạo phần tóm tắt mô hình cho khách hàng hoặc tài liệu bạn đang viết mô tả kết quả. Tuy vậy, mình khuyến khích bạn học cách tự thao tác và khám phá đầu ra của mô hình để tự xây dựng phương pháp trình bày kết quả của riêng bạn. Phát triển cách làm độc đáo của riêng mình sẽ giúp bạn nổi bật hơn với vai trò data scientist!
Bạn sẽ làm những việc sau:
- Xem xét các ước lượng của mô hình.
- Vẽ dữ liệu và khớp một
glmcho từng nhóm tuổi. Dù không hoàn toàn giống với đầu ra củaglmer(), phép xấp xỉ này giúp trình bày kết quả theo cách trực quan, dễ hiểu.
Bài tập này là một phần của khóa học
Mô hình phân cấp và hiệu ứng hỗn hợp trong R
Hướng dẫn bài tập
- Trích xuất các ước lượng hiệu ứng cố định từ
model_outbằngfixef(). - Trích xuất các ước lượng hiệu ứng ngẫu nhiên từ
model_outbằngranef(). - Chạy mã để vẽ dữ liệu bằng các phương thức của
ggplot2.
Bài tập tương tác thực hành trực tiếp
Hãy thử làm bài tập này bằng cách hoàn thành đoạn mã mẫu này.
# 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()