Đánh giá mô hình cuối cùng

R dành cho người dùng SAS

Melinda Higgins, PhD

Research Professor/Senior Biostatistician Emory University

Bài tập cuối

Tổng kết khóa học

Chạy mô hình hồi quy
- với bộ dự báo khác nhau
- cho các nhóm khác nhau
Chọn mô hình tốt nhất
- lưu kết quả mô hình
- trích xuất và hiển thị thống kê độ khớp

Thể hiện kỹ năng của bạn

Đánh giá và so sánh mô hình
- bằng trực quan hóa đồ thị
Báo cáo liên hệ tốt nhất
- giữa các biến
- tổng thể và theo nhóm

So sánh mô hình

# Chạy lm() cho diffht theo bmi, lưu mô hình
lmdiffhtbmi <- lm(diffht ~ bmi,
                  data = daviskeep)

# Chạy lm() cho diffht theo weight, lưu mô hình
lmdiffhtwt <- lm(diffht ~ weight,
                 data = daviskeep)

# Chạy summary() cho từng mô hình, lưu kết quả
smrylmdiffhtbmi <- summary(lmdiffhtbmi)
smrylmdiffhtwt <- summary(lmdiffhtwt)

So sánh mô hình

# Hiển thị r.squared cho mô hình weight
smrylmdiffhtwt$r.squared

# Hiển thị r.squared cho mô hình bmi
smrylmdiffhtbmi$r.squared

# So sánh AIC của cả hai mô hình
AIC(lmdiffhtbmi, lmdiffhtwt)


[1] 0.003281645


[1] 0.00121824

            df      AIC
lmdiffhtbmi  3 788.0816
lmdiffhtwt   3 787.7052

Mô hình theo nhóm - nam vs nữ

# Vẽ diffht theo weight phân theo giới
ggplot(daviskeep,
       aes(diffht, weight)) +
  geom_point() +
  geom_smooth(method = "lm") +
  facet_wrap(vars(sex)) +
  ggtitle("Chênh lệch chiều cao
          được dự đoán bởi cân nặng,
          mô hình khớp theo giới")

biểu đồ diffht theo weight theo giới

Hồi quy trên tập con

tùy chọn where của sas proc reg giống tùy chọn subset cho hàm R lm

Khớp mô hình cho các tập con

# lm() diffht theo weight cho nữ
lmdiffhtwtF <- lm(diffht ~ weight,
                  subset = (sex == "F"),
                  data = daviskeep)

# lm() diffht theo weight cho nam
lmdiffhtwtM <- lm(diffht ~ weight,
                  subset = (sex == "M"),
                  data = daviskeep)

# Chạy summary() cho từng mô hình, lưu kết quả
smrylmdiffhtwtF <- summary(lmdiffhtwtF)
smrylmdiffhtwtM <- summary(lmdiffhtwtM)

Khớp mô hình cho các tập con

# r.squared cho mô hình chỉ nữ
smrylmdiffhtwtF$r.squared

# r.squared cho mô hình chỉ nam
smrylmdiffhtwtM$r.squared


[1] 4.00807e-05


[1] 0.00804139

Hãy tổng kết bằng cách xây vài mô hình dự đoán tuổi bào ngư!

R dành cho người dùng SAS