Numărul de decese

Profesorul Cynthiei îi cere acum să exploreze numărul de decese din tabela de mortalitate. Decesele la vârsta \(x\) sunt notate cu \(d_x\). Poți reprezenta grafic numărul de decese \(d_x\) în funcție de vârsta \(x\)?

În plus, Cynthia a aflat că \(d_x\) reprezintă numărul așteptat de persoane care decedează la vârsta \(x\) dintr-un grup de \(\ell_x\) supraviețuitori. Vei genera eșantioane din distribuția binomială pentru numărul de decese la vârsta \(x\) și vei compara datele simulate cu valorile înregistrate ale \(d_x\).

Obiectul life_table este preîncărcat în spațiul de lucru R, împreună cu coloanele extrase age, qx, lx și dx.

Acest exercițiu face parte din cursul

Evaluarea produselor de asigurare de viață în R

Vezi cursul

Instrucțiuni pentru exercițiu

Reprezintă grafic dx față de age pentru a examina numărul de decese pe vârste. Folosind type = "h" obții linii verticale în loc de puncte.
Simulează numărul de decese la fiecare vârstă \(x\) pe baza numărului de supraviețuitori lx și a ratelor de mortalitate qx. Folosește rbinom() (docs) pentru a genera valori vectorizat dintr-o distribuție binomială.
Reprezintă grafic numărul simulat de decese sims peste graficul existent folosind points() (docs). Specifică pch = 4 pentru a folosi cruci ca simboluri.

Exercițiu interactiv practic

Încearcă acest exercițiu completând acest cod de exemplu.

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

Editează și rulează codul