El número de fallecimientos

La profesora de Cynthia ahora le pide que explore el número de fallecimientos en la tabla de vida. Los fallecimientos a la edad \(x\) se denotan como \(d_x\). ¿Puedes representar el número de fallecimientos \(d_x\) por edad \(x\)?

Además, Cynthia aprendió que \(d_x\) es el número esperado de personas que fallecen a la edad \(x\) de un grupo de \(\ell_x\) supervivientes. Vas a generar muestras binomiales para el número de fallecimientos a la edad \(x\) y comparar los datos simulados con el \(d_x\) registrado.

El objeto life_table está precargado en tu espacio de trabajo de R, así como las columnas extraídas age, qx, lx y dx.

Este ejercicio forma parte del curso

Valoración de productos de seguros de vida en R

Ver curso

Instrucciones del ejercicio

Representa dx frente a age para examinar el número de fallecimientos por edad. Usando type = "h" obtienes líneas verticales en lugar de puntos.
Simula el número de fallecimientos en cada edad \(x\) en función del número de supervivientes lx y de las tasas de mortalidad qx. Usa rbinom() (docs) para generar de forma vectorizada a partir de una distribución binomial.
Dibuja el número simulado de fallecimientos sims sobre el gráfico existente usando points() (docs). Se utilizan cruces como símbolos especificando pch = 4.

Ejercicio interactivo práctico

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

# Plot the number of deaths dx by age
plot(___, ___, 
    type = "h", 
    pch = 20, 
    xlab = "Age x", 
    ylab = expression("d"[x]),
    main = "Number of deaths (Belgium, females, 1999)")

# Simulate the number of deaths using a binomial distribution
sims <- ___(n = ___, size = ___, prob = ___)
  
# Plot the simulated number of deaths on top of the previous graph
points(___, ___, 
    pch = 4, 
    col = "red")

Editar y ejecutar código