Počet úmrtí

Cynthiin profesor ji nyní žádá, aby prozkoumala počet úmrtí v úmrtnostní tabulce. Úmrtí ve věku \(x\) jsou označena jako \(d_x\). Zkus zobrazit počet úmrtí \(d_x\) podle věku \(x\).

Cynthia se také dozvěděla, že \(d_x\) je očekávaný počet lidí, kteří zemřou ve věku \(x\), z celkové skupiny \(\ell_x\) přeživších. Vygeneruješ binomické vzorky pro počet úmrtí ve věku \(x\) a porovnáš simulovaná data se zaregistrovanými hodnotami \(d_x\).

Objekt life_table je předem načten ve tvém R pracovním prostředí spolu s extrahovanými sloupci age, qx, lx a dx.

Toto cvičení je součástí kurzu

Oceňování produktů životního pojištění v R

Zobrazit kurz

Pokyny k cvičení

Vykresli dx v závislosti na age, abys prozkoumal/a počet úmrtí podle věku. Použití type = "h" zobrazí svislé čáry místo bodů.
Simuluj počet úmrtí v každém věku \(x\) na základě počtu přeživších lx a měr úmrtnosti qx. Použij rbinom() (docs) pro vektorizované generování z binomického rozdělení.
Vykresli simulovaný počet úmrtí sims přes existující graf pomocí points() (docs). Jako symboly se použijí křížky zadáním pch = 4.

Interaktivní cvičení na vyzkoušení si v praxi

Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.

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

Upravit a spustit kód