Ratele de mortalitate în timp

Cynthia descarcă cele mai recente date de mortalitate pentru Belgia din Human Mortality Database (HMD). Acest set de date este preîncărcat ca life_table.

În acest exercițiu, vei folosi două funcții uzuale din R. with() (docs) îți permite să evaluezi o expresie R într-un mediu local construit dintr-un data frame. Astfel, nu mai trebuie să scrii repetat life_table$ pentru a extrage coloane. De exemplu, ratele de mortalitate logaritmate ale persoanelor de 18 ani de-a lungul anilor pot fi extrase astfel:

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

with() este deosebit de utilă în combinație cu subset() (docs). De exemplu, rata de mortalitate logaritmată a unei persoane de 18 ani din tabelul de mortalitate al anului 1999 poate fi obținută astfel:

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

Acest exercițiu face parte din cursul

Evaluarea produselor de asigurare de viață în R

Vezi cursul

Instrucțiuni pentru exercițiu

Explorează life_table. Afișează primele 6 rânduri folosind head() și calculează intervalul (range()) (docs) al variabilei year.
Completează codul folosind corect subset() pe life_table, astfel încât să fie reprezentate grafic ratele de mortalitate ale unei femei de 18 ani de-a lungul anilor.
Folosește din nou subset() pentru a selecta tabelul de mortalitate din 1950 și trasează curba ratei de mortalitate pentru acel an.

Exercițiu interactiv practic

Încearcă acest exercițiu completând acest cod de exemplu.

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

Editează și rulează codul