Aan de slagBegin gratis

Eén jaar aan sterftedata visualiseren

Nu je wat gevoel hebt gekregen bij de life table van 1999, kun je de evolutie in vrouwelijke sterftekansen \(q_x\) uitzetten als functie van de leeftijd \(x\). Cynthia's docent zet in zijn slides meestal de log van de sterftekansen uit. Kun jij haar overtuigen dat dit inderdaad een goede aanpak is?

Het object life_table en de vectoren age en qx uit de vorige oefening zijn alvast ingeladen.

Deze oefening maakt deel uit van de cursus

Waardering van levensverzekeringsproducten in R

Bekijk cursus

Oefeninstructies

  • Gebruik plot() (docs) om de sterftekansen voor vrouwen in 1999 te plotten. Het eerste argument is age en het tweede zijn de sterftekansen qx.
  • Gebruik een tweede aanroep van plot() om de sterftekansen voor vrouwen in 1999 op een logaritmische schaal weer te geven. In vergelijking met het vorige commando moet je het tweede argument transformeren door eerst de log() (docs) te nemen.

Interactieve oefening met praktijkervaring

Probeer deze oefening door deze voorbeeldcode aan te vullen.

# Plot the female mortality rates in the year 1999
plot(___, ___, 
    main = "Mortality rates (Belgium, females, 1999)", 
    xlab = "Age x", 
    ylab = expression(paste("Mortality rate ", q[x])), 
    type = "l")

# Plot the logarithm of the female mortality rates in the year 1999
plot(___, ___, 
    main = "Log mortality rates (Belgium, females, 1999)", 
    xlab = "Age x", 
    ylab = expression(paste("Log mortality rate ", log(q[x]))), 
    type = "l")
Code bewerken en uitvoeren