#con los datos de mortality.org

#Spain, Death rates (period 1x1) Last modified: 22-Feb-2008, MPv5 (May07)
spain<-read.demogdata("Mx_1x1SPAIN.txt", "Exposures_1x1SPAIN.txt",type="mortality", label="Spain", skip = 2,
popskip = 2)


#Seleccionando tasas de 1975 a  2006
spain.desde1975 <- extract.years(spain,1975:2006)
 #Seleccionando poblaciòn con edad mayor o igual a 1965 y tasas de 1975 a 2006
spain.mas65 <- extract.ages(spain,65:110,FALSE)
spain.mas65.desde1975 <- extract.years(spain.mas65,1975:2006)


#library(demography)


#Modelo Lee Carter  hombres mayores a 65 años, años 1975 - 2006
spain.mas65.desde1975.lcaM <- lca(spain.mas65.desde1975, series="male",adjust="dt")
summary(spain.mas65.desde1975.lcaM)

#Proyecciones a 50 años
spain.mas65.desde1975.fcastM <- forecast(spain.mas65.desde1975.lcaM, 50)

#Grafico de las proyecciones
x11()
plot(spain.mas65.desde1975.fcastM)


#Grafico de   alpha, kappa
x11()
plot(spain.mas65.desde1975.fcastM,'c')


#method based on life expectancy
spain.LChom <- lca(spain.desde1975,series="male",adjust="e0",years=1975:2006)
summary(spain.LChom)

#Predicción LC de la esperanza de vida
spain.LChom.fcast <- forecast(spain.LChom,jumpchoice="actual")
spain.LChom.e0.f <- life.expectancy(spain.LChom.fcast)
par(mfrow=c(1,1))
plot(spain.LChom.e0.f)


#Predicicón e intervalo LC de la esperanza de vida
spain.LChom.fdm <- fdm(spain.desde1975,series="male")
spain.hom.fcast <- forecast(spain.LChom.fdm, 50)
e0.fcast <- e0(spain.hom.fcast,PI=TRUE,nsim=100)
plot(e0.fcast)


#predicción LC de la esperanza de vida restante
#modelo lC hombres        0-110 años, 1908-2006
spainH.LC <- lca(spain, series="male")

#proyecciones a 50 años
spainH.LC.fcast <- forecast(spainH.LC)
spainH.LC.fcast.f <- lifetable(spainH.LC.fcast)

par(mfrow=c(1,1))
plot(spainH.LC.fcast.f )