Tỷ lệ tử vong theo thời gian
Cynthia tải bộ dữ liệu tử vong mới nhất của Bỉ từ Human Mortality Database (HMD). Bộ dữ liệu này đã được nạp sẵn với tên life_table.
Trong bài tập này, bạn sẽ dùng hai hàm R thường gặp. with() (tài liệu) cho phép bạn đánh giá một biểu thức R trong môi trường cục bộ được tạo từ một data frame. Nhờ vậy, bạn không phải gõ lặp lại life_table$ để trích xuất cột. Ví dụ, log của tỷ lệ tử vong của người 18 tuổi qua các năm có thể được trích xuất bằng:
with(life_table, log(qx[age == 18]))
with() đặc biệt hữu ích khi kết hợp với subset() (tài liệu). Ví dụ, log của tỷ lệ tử vong của người 18 tuổi từ bảng sống của năm 1999 có thể lấy bằng:
with(subset(life_table, year == 1999), log(qx[age == 18]))
Bài tập này là một phần của khóa học
Định giá Sản phẩm Bảo hiểm Nhân thọ bằng R
Hướng dẫn bài tập
- Khám phá
life_table. In 6 dòng đầu bằnghead()và tínhrange()(tài liệu) của biếnyear. - Hoàn thiện mã với cách dùng
subset()phù hợp trênlife_tableđể vẽ tỷ lệ tử vong của nữ 18 tuổi qua các năm. - Tiếp tục dùng
subset()để chọn bảng sống của năm 1950 và vẽ đường cong tỷ lệ tử vong trong năm đó.
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.
# Explore life_table
___
___
# Plot the logarithm of the female mortality rates for (18) by year
with(___(___, ___),
plot(year, log(qx),
type = "l", main = "Log mortality rates (Belgium, females, 18-year-old)",
xlab = "Year t", ylab = expression(paste("Log mortality rate ", log(q[18])))))
# Plot the logarithm of the female mortality rates in the year 1950 by age
with(___(___, ___),
plot(age, log(qx),
type = "l", main = "Log mortality rates (Belgium, females, 1950)",
xlab = "Age x", ylab = expression(paste("Log mortality rate ", log(q[x])))))