Vizualizace střední délky života podle věku
Cynthia se rozhodla zobrazit střední délku života jako funkci věku. Pomůžeš jí takový graf sestavit?
Nejefektivnější přístup v R spočívá v tom, že nejprve napíšeš funkci, která pro daný věk a úmrtnostní tabulku vypočítá zkrácenou střední délku života. Tuto funkci pak aplikuješ na všechny věky v tabulce. Pomocí sapply() (docs) se výstup zjednodušší na co nejzákladnější datovou strukturu.
Předpřipravený objekt life_table obsahuje periodickou úmrtnostní tabulku pro ženy v Belgii za rok 1999.
Toto cvičení je součástí kurzu
Oceňování produktů životního pojištění v R
Pokyny k cvičení
- Dokonči kód definující funkci
curtate_future_lifetime(). - Vytvoř vektor
agesextrahováním sloupceagezlife_table. - Pomocí
sapply()s argumentyages,curtate_future_lifetimealife_tablevypočítej zkrácenou střední délku života pro všechny věky v úmrtnostní tabulce. - Vykresli
future_lifetimesv závislosti naages.
Interaktivní cvičení na vyzkoušení si v praxi
Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.
# Function to compute the curtate expected future lifetime for a given age and life table
curtate_future_lifetime <- function(age, life_table) {
px <- ___
kpx <- ___(px[(___):length(px)])
___(___)
}
# Vector of ages
ages <- ___
# Curtate future lifetimes for all ages
future_lifetimes <- ___(___, ___, ___)
# Future lifetime by age
plot(___, ___, type = 'l', lwd = 2, col = "green", xlab = "Age x", ylab = "Future lifetime", main = "Future lifetime by age")