Dự đoán từ mô hình thất nghiệp
Trong bài tập này, bạn sẽ dùng mô hình thất nghiệp unemployment_model để dự đoán từ dữ liệu unemployment, và so sánh tỷ lệ thất nghiệp nữ được dự đoán với tỷ lệ thực tế quan sát được trên dữ liệu huấn luyện unemployment. Bạn cũng sẽ dùng mô hình để dự đoán trên dữ liệu mới trong newrates, dữ liệu này chỉ gồm một quan sát, với tỷ lệ thất nghiệp nam là 5%.
Giao diện predict() (docs) cho các mô hình lm có dạng
predict(model, newdata)
Bạn sẽ dùng gói ggplot2 để vẽ biểu đồ, nên bạn sẽ thêm cột dự đoán vào data frame unemployment. Bạn sẽ vẽ kết quả thực tế so với dự đoán và so sánh chúng với đường biểu diễn các dự đoán hoàn hảo (tức là khi giá trị thực tế bằng đúng giá trị dự đoán).
Lệnh ggplot2 để vẽ biểu đồ scatter của dframe$outcome so với dframe$pred (pred trên trục x, outcome trên trục y), kèm một đường màu xanh nơi outcome == pred như sau:
ggplot(dframe, aes(x = pred, y = outcome)) +
geom_point() +
geom_abline(color = "blue")
unemployment, unemployment_model, và newrates đã được nạp sẵn cho bạ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
- Dùng
predict()để dự đoán tỷ lệ thất nghiệp nữ từ dữ liệuunemployment. Gán vào một cột mới:prediction. - Dùng lệnh
library()để nạp góiggplot2. - Dùng
ggplot()để so sánh dự đoán với tỷ lệ thất nghiệp thực tế. Đặt dự đoán trên trục x. Kết quả gần với đường dự đoán hoàn hảo đến mức nào? - Dùng data frame
newratesđể dự đoán tỷ lệ thất nghiệp nữ kỳ vọng khi tỷ lệ thất nghiệp nam là 5%. Gán câu trả lời vào biếnpredvà in ra.
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.
# unemployment is available
summary(unemployment)
# newrates is available
newrates
# Predict female unemployment in the unemployment dataset
unemployment$prediction <- ___
# Load the ggplot2 package
___
# Make a plot to compare predictions to actual (prediction on x axis).
ggplot(___, aes(x = ___, y = ___)) +
___ +
geom_abline(color = "blue")
# Predict female unemployment rate when male unemployment is 5%
pred <- ___
pred