Wskaźniki śmiertelności w czasie
Cynthia pobiera najbardziej aktualne dane dotyczące śmiertelności dla Belgii z Human Mortality Database (HMD). Ten zbiór danych jest wstępnie załadowany jako life_table.
W tym ćwiczeniu skorzystasz z dwóch popularnych funkcji języka R. Funkcja with() (dokumentacja) pozwala obliczyć wyrażenie R w lokalnym środowisku zbudowanym na podstawie ramki danych. Dzięki temu nie trzeba wielokrotnie wpisywać life_table$, aby odwołać się do kolumn. Na przykład logarytmiczne wskaźniki śmiertelności 18-latków na przestrzeni lat można uzyskać w następujący sposób:
with(life_table, log(qx[age == 18]))
Funkcja with() jest szczególnie przydatna w połączeniu z subset() (dokumentacja). Na przykład logarytmiczny wskaźnik śmiertelności 18-latka z tablicy trwania życia dla roku 1999 można uzyskać tak:
with(subset(life_table, year == 1999), log(qx[age == 18]))
To ćwiczenie jest częścią kursu
Wycena produktów ubezpieczeń na życie w R
Instrukcje do ćwiczenia
- Zbadaj
life_table. Wyświetl pierwsze 6 wierszy za pomocąhead()i obliczrange()(dokumentacja) zmiennejyear. - Uzupełnij kod, używając odpowiednio
subset()nalife_table, tak aby wykreślić wskaźniki śmiertelności 18-letniej kobiety na przestrzeni lat. - Ponownie użyj
subset(), aby wybrać tablicę trwania życia z roku 1950, i wykreśl krzywą wskaźnika śmiertelności dla tego roku.
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
# Explore life_table
___
___
# Plot the logarithm of the female mortality rates for (18) by year
with(___(___, ___),
plot(year, log(qx),
type = "l", main = "Log mortality rates (Belgium, females, 18-year-old)",
xlab = "Year t", ylab = expression(paste("Log mortality rate ", log(q[18])))))
# Plot the logarithm of the female mortality rates in the year 1950 by age
with(___(___, ___),
plot(age, log(qx),
type = "l", main = "Log mortality rates (Belgium, females, 1950)",
xlab = "Age x", ylab = expression(paste("Log mortality rate ", log(q[x])))))