Jumlah kematian

Profesor Cynthia kini memintanya menelaah jumlah kematian dalam tabel mortalitas. Kematian pada usia \(x\) diberi label \(d_x\). Dapatkah Anda memvisualisasikan jumlah kematian \(d_x\) menurut usia \(x\)?

Selain itu, Cynthia mempelajari bahwa \(d_x\) adalah nilai harapan jumlah orang yang meninggal pada usia \(x\) dari kelompok beranggotakan \(\ell_x\) penyintas. Anda akan menghasilkan sampel binomial untuk jumlah kematian pada usia \(x\) dan membandingkan data simulasi dengan \(d_x\) yang tercatat.

Objek life_table sudah dimuat di ruang kerja R Anda, begitu juga kolom hasil ekstraksi age, qx, lx, dan dx.

Latihan ini merupakan bagian dari kursus

Penilaian Produk Asuransi Jiwa di R

Lihat Kursus

Instruksi latihan

Plot dx terhadap age untuk menelaah jumlah kematian menurut usia. Dengan type = "h" Anda akan memperoleh garis vertikal alih-alih titik.
Simulasikan jumlah kematian pada setiap usia \(x\) berdasarkan jumlah penyintas lx dan laju mortalitas qx. Gunakan rbinom() (docs) untuk menghasilkan secara vektorisasi dari distribusi binomial.
Plot jumlah kematian hasil simulasi sims di atas grafik yang ada menggunakan points() (docs). Gunakan tanda silang sebagai simbol dengan menentukan pch = 4.

Latihan interaktif langsung praktik

Cobalah latihan ini dengan melengkapi kode contoh ini.

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

Edit dan Jalankan Kode