datacamp-logo
Session Ready
Exercise

The number of deaths

Cynthia's professor now asks her to explore the number of deaths in the life table. The deaths at age \(x\) are labeled \(d_x\). Can you picture the number of deaths \(d_x\) by age \(x\)?

Moreover, Cynthia learned that \(d_x\) is the expected number of people who die at age \(x\) out of a group of \(\ell_x\) survivors. You will generate binomial samples for the number of deaths at age \(x\) and compare simulated data with the registered \(d_x\).

The object life_table is preloaded in your R workspace as well as the extracted columns age, qx, lx and dx.

Instructions
100 XP
  • Plot dx versus age to examine the number of deaths by age. Using type = "h" you obtain vertical lines instead of points.
  • Simulate the number of deaths at each age \(x\) based on the number of survivors lx and the mortality rates qx. Use rbinom() to generate in a vectorized way from a binomial distribution.
  • Plot the simulated number of deaths sims on top of the existing graph using points(). Crosses are used as symbols by specifying pch = 4.