Số ca tử vong

Giảng viên của Cynthia đề nghị cô ấy khám phá số ca tử vong trong bảng tử vong. Số ca tử vong ở tuổi \(x\) được ký hiệu là \(d_x\). Bạn có thể hình dung số ca tử vong \(d_x\) theo tuổi \(x\) không?

Ngoài ra, Cynthia biết rằng \(d_x\) là số người kỳ vọng sẽ qua đời ở tuổi \(x\) trong một nhóm có \(\ell_x\) người còn sống. Bạn sẽ tạo mẫu nhị thức cho số ca tử vong ở tuổi \(x\) và so sánh dữ liệu mô phỏng với \(d_x\) đã ghi nhận.

Đối tượng life_table đã được nạp sẵn trong không gian làm việc R của bạn, cùng với các cột đã trích xuất age, qx, lx và dx.

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

Xem khóa học

Hướng dẫn bài tập

Vẽ dx theo age để xem số ca tử vong theo độ tuổi. Dùng type = "h" để hiện các đường thẳng đứng thay vì các điểm.
Mô phỏng số ca tử vong tại mỗi tuổi \(x\) dựa trên số người còn sống lx và xác suất tử vong qx. Dùng rbinom() (docs) để sinh theo cách vector hóa từ phân phối nhị thức.
Vẽ chồng số ca tử vong mô phỏng sims lên đồ thị hiện có bằng points() (docs). Dùng dấu gạch chéo làm ký hiệu bằng cách đặt pch = 4.

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.

# Plot the number of deaths dx by age
plot(___, ___, 
    type = "h", 
    pch = 20, 
    xlab = "Age x", 
    ylab = expression("d"[x]),
    main = "Number of deaths (Belgium, females, 1999)")

# Simulate the number of deaths using a binomial distribution
sims <- ___(n = ___, size = ___, prob = ___)
  
# Plot the simulated number of deaths on top of the previous graph
points(___, ___, 
    pch = 4, 
    col = "red")

Chỉnh sửa và Chạy Mã