Dự đoán khả năng sống sót của chim sẻ
Trong bài này, bạn sẽ dự đoán xác suất sống sót bằng mô hình sống sót của chim sẻ từ bài trước.
Hãy nhớ rằng khi gọi predict() (docs) để lấy xác suất dự đoán từ mô hình glm(), bạn phải chỉ định muốn lấy phản hồi:
predict(model, type = "response")
Nếu không, predict() trên mô hình logistic regression sẽ trả về log-odds dự đoán của biến cố, chứ không phải xác suất.
Bạn cũng sẽ dùng hàm GainCurvePlot() (docs) để vẽ đường cong lợi ích (gain curve) từ các dự đoán của mô hình. Nếu đường cong lợi ích của mô hình gần với đường cong lý tưởng ("wizard"), thì mô hình đã sắp xếp các cá thể chim sẻ tốt: tức là mô hình dự đoán những con thực sự sống sót có xác suất sống sót cao hơn. Các tham số đưa vào hàm GainCurvePlot() gồm:
frame: data frame có cột dự đoán và cột nhãn thực tếxvar: tên cột dự đoán (dưới dạng chuỗi)truthVar: tên cột nhãn kết cục thực tế (dưới dạng chuỗi)title: tiêu đề của biểu đồ (dưới dạng chuỗi)
GainCurvePlot(frame, xvar, truthVar, title)
Data frame sparrow và mô hình sparrow_model đã được nạp sẵn.
Bài tập này là một phần của khóa học
Học có giám sát với R: Hồi quy
Hướng dẫn bài tập
- Tạo một cột mới trong
sparrowtên làpredchứa các dự đoán trên dữ liệu huấn luyện. - Gọi
GainCurvePlot()để tạo đường cong lợi ích của các dự đoán. Mô hình có sắp xếp chim sẻ theo việc chúng thực sự sống sót hay không một cách tốt không?
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.
# sparrow is available
summary(sparrow)
# sparrow_model is available
summary(sparrow_model)
# Make predictions
sparrow$pred <- ___
# Look at gain curve
___(___, ___, ___, "sparrow survival model")