Wizualizacja oczekiwanego dalszego trwania życia według wieku
Cynthia postanawia zwizualizować oczekiwane dalsze trwanie życia w zależności od wieku. Czy możesz jej pomóc skonstruować taki wykres?
Aby zrobić to sprawnie w R, najpierw napisz funkcję obliczającą skrócone oczekiwane dalsze trwanie życia dla danego wieku i tablicy trwania życia. Następnie możesz zastosować tę funkcję dla wszystkich wieku w tablicy. Użycie sapply() (dokumentacja) sprowadza wynik do możliwie najprostszej struktury danych.
Wstępnie załadowany obiekt life_table zawiera tablicę trwania życia dla kobiet w Belgii z 1999 roku (dane okresowe).
To ćwiczenie jest częścią kursu
Wycena produktów ubezpieczeń na życie w R
Instrukcje do ćwiczenia
- Uzupełnij kod definiujący funkcję
curtate_future_lifetime(). - Utwórz wektor
ages, wyodrębniając kolumnęagezlife_table. - Użyj
sapply()z argumentamiages,curtate_future_lifetimeorazlife_table, aby obliczyć skrócone oczekiwane dalsze trwanie życia dla wszystkich wieku w tablicy. - Narysuj wykres
future_lifetimesw zależności odages.
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
# 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")