Tingkat mortalitas dari waktu ke waktu
Cynthia mengunduh data mortalitas terbaru untuk Belgia dari Human Mortality Database (HMD). Himpunan data ini sudah dimuat sebagai life_table.
Dalam latihan ini, Anda akan menggunakan dua fungsi R yang umum. with() (docs) memungkinkan Anda mengevaluasi sebuah ekspresi R dalam lingkungan lokal yang dibentuk dari sebuah data frame. Ini menghindarkan Anda dari keharusan berulang kali mengetik life_table$ untuk mengambil kolom. Misalnya, log tingkat mortalitas usia 18 tahun sepanjang tahun dapat diambil dengan:
with(life_table, log(qx[age == 18]))
with() sangat berguna bila dikombinasikan dengan subset() (docs). Misalnya, log tingkat mortalitas usia 18 tahun dari life table tahun 1999 dapat diperoleh dengan:
with(subset(life_table, year == 1999), log(qx[age == 18]))
Latihan ini merupakan bagian dari kursus
Penilaian Produk Asuransi Jiwa di R
Instruksi latihan
- Eksplorasi
life_table. Cetak 6 baris pertama menggunakanhead()dan hitungrange()(docs) dari variabelyear. - Lengkapi kode dengan penggunaan
subset()yang tepat padalife_tablesehingga tingkat mortalitas perempuan usia 18 tahun sepanjang tahun dapat dipetakan. - Gunakan kembali
subset()untuk memilih life table tahun 1950 dan plot kurva tingkat mortalitas pada tahun tersebut.
Latihan interaktif langsung praktik
Cobalah latihan ini dengan melengkapi kode contoh ini.
# 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])))))