Tassi di mortalità nel tempo
Cynthia scarica i dati di mortalità più aggiornati per il Belgio dall'Human Mortality Database (HMD). Questo insieme di dati è già caricato come life_table.
In questo esercizio userai due funzioni R molto comuni. with() (docs) ti permette di valutare un'espressione R in un ambiente locale costruito a partire da un data frame. In questo modo eviti di digitare ripetutamente life_table$ per estrarre le colonne. Per esempio, i log dei tassi di mortalità dei 18enni nel corso degli anni si possono estrarre con:
with(life_table, log(qx[age == 18]))
with() è particolarmente comoda in combinazione con subset() (docs). Per esempio, il log del tasso di mortalità dei 18enni dalla tavola di mortalità dell'anno 1999 può essere ottenuto con:
with(subset(life_table, year == 1999), log(qx[age == 18]))
Questo esercizio fa parte del corso
Valutazione dei prodotti assicurativi sulla vita in R
Istruzioni dell'esercizio
- Esplora
life_table. Stampa le prime 6 righe conhead()e calcola l'intervallo (range()docs) della variabileyear. - Completa il codice usando correttamente
subset()sulife_tablein modo da tracciare i tassi di mortalità di una donna di 18 anni nel corso degli anni. - Usa di nuovo
subset()per selezionare la tavola di mortalità del 1950 e traccia la curva dei tassi di mortalità in quell'anno.
esercizio interattivo pratico
Prova questo esercizio completando questo codice di esempio.
# 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])))))