O número de mortes

A professora da Cynthia agora pede que ela explore o número de mortes na tábua de vida. As mortes na idade \(x\) são indicadas por \(d_x\). Você consegue visualizar o número de mortes \(d_x\) por idade \(x\)?

Além disso, Cynthia aprendeu que \(d_x\) é o número esperado de pessoas que morrem na idade \(x\) dentre um grupo de \(\ell_x\) sobreviventes. Você vai gerar amostras binomiais para o número de mortes na idade \(x\) e comparar os dados simulados com o \(d_x\) registrado.

O objeto life_table já está carregado no seu espaço de trabalho do R, assim como as colunas extraídas age, qx, lx e dx.

Este exercicio faz parte do curso

Avaliação de Produtos de Seguro de Vida em R

Ver curso

Instruções do exercicio

Faça o gráfico de dx versus age para examinar o número de mortes por idade. Usando type = "h", você obtém linhas verticais em vez de pontos.
Simule o número de mortes em cada idade \(x\) com base no número de sobreviventes lx e nas taxas de mortalidade qx. Use rbinom() (docs) para gerar, de forma vetorizada, a partir de uma distribuição binomial.
Plote o número simulado de mortes sims em cima do gráfico existente usando points() (docs). Cruzes são usadas como símbolos especificando pch = 4.

exercicio interativo prático

Tente este exercicio completando este código de exemplo.

# 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 e Executar Código