Antalet dödsfall

Cynthias professor ber henne nu att undersöka antalet dödsfall i livstabellen. Dödsfallen vid ålder \(x\) betecknas \(d_x\). Kan du visualisera antalet dödsfall \(d_x\) per ålder \(x\)?

Cynthia har dessutom lärt sig att \(d_x\) är det förväntade antalet personer som dör vid ålder \(x\) ur en grupp av \(\ell_x\) överlevande. Du kommer att generera binomialstickprov för antalet dödsfall vid ålder \(x\) och jämföra simulerade data med de registrerade \(d_x\).

Objektet life_table är förladdat i din R-miljö, liksom de extraherade kolumnerna age, qx, lx och dx.

Den här övningen är en del av kursen

Värdering av livförsäkringsprodukter i R

Visa kurs

Övningsinstruktioner

Plotta dx mot age för att undersöka antalet dödsfall per ålder. Med type = "h" får du vertikala linjer i stället för punkter.
Simulera antalet dödsfall vid varje ålder \(x\) baserat på antalet överlevande lx och dödlighetstalen qx. Använd rbinom() (docs) för att generera värden på ett vektoriserat sätt från en binomialfördelning.
Plotta det simulerade antalet dödsfall sims ovanpå det befintliga diagrammet med points() (docs). Kryss används som symboler genom att ange pch = 4.

Interaktiv övning med praktiskt arbete

Testa den här övningen genom att slutföra den här exempelkoden.

# 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")

Redigera och kör kod