Míry úmrtnosti v čase

Cynthia si stáhne nejaktuálnější data o úmrtnosti pro Belgii z Human Mortality Database (HMD). Tato datová sada je předem načtena jako life_table.

V tomto cvičení využiješ dvě běžné funkce R. with() (dokumentace) umožňuje vyhodnotit výraz v R v lokálním prostředí vytvořeném z datového rámce. Díky tomu nemusíš opakovaně psát life_table$ při přístupu ke sloupcům. Například logaritmické míry úmrtnosti 18letých v průběhu let lze získat takto:

with(life_table, log(qx[age == 18]))

with() se hodí zejména v kombinaci s subset() (dokumentace). Například logaritmickou míru úmrtnosti 18letých z tabulky života pro rok 1999 získáš takto:

with(subset(life_table, year == 1999), log(qx[age == 18]))

Toto cvičení je součástí kurzu

Oceňování produktů životního pojištění v R

Zobrazit kurz

Pokyny k cvičení

Prozkoumej life_table. Vypiš prvních 6 řádků pomocí head() a zjisti rozsah proměnné year pomocí range() (dokumentace).
Doplň kód se správným použitím subset() na life_table tak, aby se zobrazil graf měr úmrtnosti 18leté ženy v průběhu let.
Opět použij subset() k výběru tabulky života pro rok 1950 a vykresli křivku měr úmrtnosti v daném roce.

Interaktivní cvičení na vyzkoušení si v praxi

Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.

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

Upravit a spustit kód