Tasas de mortalidad a lo largo del tiempo

Cynthia descarga los datos de mortalidad más recientes de Bélgica del Human Mortality Database (HMD). Este conjunto de datos está precargado como life_table.

En este ejercicio, usarás dos funciones habituales de R. with() (docs) te permite evaluar una expresión de R en un entorno local construido a partir de un data frame. Esto evita tener que escribir repetidamente life_table$ para extraer columnas. Por ejemplo, las tasas de mortalidad en logaritmos de las personas de 18 años a lo largo de los años se pueden extraer con:

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

with() es especialmente útil en combinación con subset() (docs). Por ejemplo, la tasa de mortalidad en logaritmos de las personas de 18 años de la tabla de vida del año 1999 se puede obtener con:

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

Este ejercicio forma parte del curso

Valoración de productos de seguros de vida en R

Ver curso

Instrucciones del ejercicio

Explora life_table. Imprime las primeras 6 filas con head() y calcula el range() (docs) de la variable year.
Completa el código usando subset() sobre life_table para que se representen las tasas de mortalidad de una mujer de 18 años a lo largo de los años.
Vuelve a usar subset() para seleccionar la tabla de vida de 1950 y dibuja la curva de la tasa de mortalidad en ese año.

Ejercicio interactivo práctico

Prueba este ejercicio y completa el código de muestra.

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

Editar y ejecutar código