Taux de mortalité dans le temps

Cynthia télécharge les données de mortalité les plus récentes pour la Belgique depuis la Human Mortality Database (HMD). Cet ensemble de données est préchargé sous le nom life_table.

Dans cet exercice, vous utiliserez deux fonctions R courantes. with() (docs) vous permet d'évaluer une expression R dans un environnement local construit à partir d'une trame de données. Cela évite de devoir taper life_table$ à répétition pour extraire des colonnes. Par exemple, les taux de mortalité en logarithme des personnes de 18 ans au fil des années peuvent être extraits avec :

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

with() est particulièrement pratique en combinaison avec subset() (docs). Par exemple, le taux de mortalité en logarithme des personnes de 18 ans tiré de la table de mortalité de l'année 1999 peut être obtenu avec :

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

Cette activité fait partie du cours

Évaluation des produits d'assurance vie en R

Voir le cours

Instructions de l’exercice

Explorez life_table. Affichez les 6 premières lignes avec head() et calculez l'étendue (range() docs) de la variable year.
Complétez le code en utilisant correctement subset() sur life_table afin de tracer les taux de mortalité d'une femme de 18 ans au fil des années.
Utilisez de nouveau subset() pour sélectionner la table de mortalité de 1950 et tracer la courbe des taux de mortalité pour cette année.

Exercice interactif pratique

Essayez cet exercice en complétant ce code d’exemple.

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

Modifier et exécuter le code