사망자 수

이제 교수님은 Cynthia에게 생명표에서 사망자 수를 살펴보라고 합니다. 나이 $x$에서의 사망자 수는 $d_x$로 표기합니다. 나이 $x$별 사망자 수 $d_x$가 어떻게 보이는지 그려 볼 수 있을까요?

또한 Cynthia는 $d_x$가, 생존자 수가 $\ell_x$명인 집단에서 나이 $x$에 사망하는 사람의 기대값이라는 것을 배웠습니다. 나이 $x$에서의 사망자 수에 대해 이항분포 표본을 생성하고, 시뮬레이션한 데이터와 기록된 $d_x$를 비교해 보세요.

life_table 객체는 R 작업 공간에 미리 로드되어 있으며, 추출된 열 age, qx, lx, dx도 함께 제공됩니다.

이 연습은 강의의 일부입니다

R로 배우는 생명보험 상품 평가

강의 보기

연습 안내

나이별 사망자 수를 확인하려면 age에 대한 dx를 그래프로 그리세요. 점 대신 수직선을 사용하려면 type = "h"를 지정합니다.
각 나이 $x$에서의 사망자 수를, 생존자 수 lx와 사망률 qx를 바탕으로 시뮬레이션하세요. 이항분포에서 벡터화된 방식으로 표본을 생성하려면 rbinom()(docs)을 사용합니다.
points()(docs)를 사용해 기존 그래프 위에 시뮬레이션된 사망자 수 sims를 표시하세요. 기호는 pch = 4로 지정해 X 표시(십자)로 사용합니다.

실습형 인터랙티브 연습

이 예제를 이 샘플 코드를 완성하여 풀어보세요.

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

코드 편집 및 실행