## Importing mortality data from the USA available on the Human Mortality Database (HMD):
data(USA)
# Example 1: -------------------------------
USA2010 = USA[USA$Year == 2010,]
ExF = USA2010$Ex.Female[1:91]
DxF = USA2010$Dx.Female[1:91]
yF = log(DxF/ExF)
fitF <- dlm(yF, M = 100)
Heatmap(fitF, x_lab = "Female expec. 2010 USA", max_age = 90)
Run the code above in your browser using DataLab