Sterftecijfers door de tijd

Cynthia downloadt de meest recente sterftedata voor België van de Human Mortality Database (HMD). Deze gegevensset is vooraf geladen als life_table.

In deze oefening gebruik je twee veelvoorkomende R-functies. Met with() (docs) kun je een R-expressie evalueren in een lokale omgeving die is opgebouwd uit een data frame. Zo hoef je niet steeds life_table$ te typen om kolommen op te vragen. Zo kun je bijvoorbeeld de log-sterftekansen van 18-jarigen door de jaren heen extraheren met:

with(life_table, log(qx[age == 18]))

with() is vooral handig in combinatie met subset() (docs). Zo kun je bijvoorbeeld de log-sterftekans van 18-jarigen uit de levens- of sterftetafel van het jaar 1999 ophalen met:

with(subset(life_table, year == 1999), log(qx[age == 18]))

Deze oefening maakt deel uit van de cursus

Waardering van levensverzekeringsproducten in R

Bekijk cursus

Oefeninstructies

Verken life_table. Print de eerste 6 rijen met head() en bereken het range() (docs) van de variabele year.
Maak de code af met correct gebruik van subset() op life_table zodat de sterftekansen van een 18-jarige vrouw door de jaren heen worden geplot.
Gebruik opnieuw subset() om de levens- of sterftetafel van 1950 te selecteren en plot in dat jaar de curve van de sterftekans.

Interactieve oefening met praktijkervaring

Probeer deze oefening door deze voorbeeldcode aan te vullen.

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

Code bewerken en uitvoeren