Zaman içinde ölüm oranları
Cynthia, Belçika için en güncel ölüm oranı verilerini Human Mortality Database (HMD) sitesinden indiriyor. Bu veri kümesi life_table olarak önceden yüklenmiş durumda.
Bu egzersizde iki yaygın R fonksiyonunu kullanacaksın. with() (docs) bir veri çerçevesinden oluşturulan yerel bir ortamda bir R ifadesini değerlendirmeni sağlar. Bu sayede sütunları çağırmak için sürekli life_table$ yazmana gerek kalmaz. Örneğin, yıllar boyunca 18 yaş grubunun log ölüm oranlarını şöyle çıkarabilirsin:
with(life_table, log(qx[age == 18]))
with() özellikle subset() (docs) ile birlikte çok kullanışlıdır. Örneğin, 1999 yılı yaşam tablosundan 18 yaş grubunun log ölüm oranını şu şekilde döndürebilirsin:
with(subset(life_table, year == 1999), log(qx[age == 18]))
Bu egzersiz, kursun bir parçasıdır
R ile Hayat Sigortası Ürünlerinin Değerlemesi
Egzersiz talimatları
life_table'ı incele. İlk 6 satırıhead()ile yazdır veyeardeğişkenininrange()(docs) değerini hesapla.- Kodda
life_tableüzerinde uygun birsubset()kullanarak, yıllar boyunca 18 yaşındaki bir kadının ölüm oranlarının grafiğini çiz. - Yine
subset()kullanarak 1950 yılına ait yaşam tablosunu seç ve o yılın ölüm oranı eğrisini çiz.
Uygulamalı etkileşimli egzersiz
Bu egzersizi bu örnek kodu tamamlayarak deneyin.
# 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])))))