Skip to content
Snippets Groups Projects
Commit 89cab4d4 authored by Reinhold Kainhofer's avatar Reinhold Kainhofer
Browse files

Austria Population MCMC: Extrapolate to age 120 (both base table and trend)

parent 363c34b8
No related branches found
No related tags found
No related merge requests found
stopifnot(require(methods), require(utils), require(MortalityTables), require(tidyverse), require(reshape2), require(pracma)) # MortalityTable classes; new; Excel reader stopifnot(require(methods), require(utils), require(MortalityTables), require(tidyverse), require(reshape2), require(pracma), require(MortalityLaws)) # MortalityTable classes; new; Excel reader
############################################################################### ###############################################################################
...@@ -38,7 +38,10 @@ mort.AT.MCMC.load = function() { ...@@ -38,7 +38,10 @@ mort.AT.MCMC.load = function() {
data = list( data = list(
dim = list(sex = "m", collar = "Gesamtbevölkerung", type = "MCMC-Fit 1980-2017", data = "MCMC", year = "1980-2017", Tafel = "MCMC-Zerlegung Bevölkerungssterblichkeit") dim = list(sex = "m", collar = "Gesamtbevölkerung", type = "MCMC-Fit 1980-2017", data = "MCMC", year = "1980-2017", Tafel = "MCMC-Zerlegung Bevölkerungssterblichkeit")
) )
) ) %>%
mT.fitExtrapolationLaw(law = "HP2", method = "LF2", fit = 80:98, extrapolate = 90:120, fadeIn = 90:99) %>%
mT.extrapolateTrendExp(idx = 92, up = TRUE)
mort.AT.MCMC[["w"]] = mortalityTable.trendProjection( mort.AT.MCMC[["w"]] = mortalityTable.trendProjection(
name = "Österreich MCMC Frauen", name = "Österreich MCMC Frauen",
ages = as.integer(dimnames(data.array)[[1]]), ages = as.integer(dimnames(data.array)[[1]]),
...@@ -49,7 +52,10 @@ mort.AT.MCMC.load = function() { ...@@ -49,7 +52,10 @@ mort.AT.MCMC.load = function() {
data = list( data = list(
dim = list(sex = "w", collar = "Gesamtbevölkerung", type = "MCMC-Fit 1980-2017", data = "MCMC", year = "1980-2017", Tafel = "MCMC-Zerlegung Bevölkerungssterblichkeit") dim = list(sex = "w", collar = "Gesamtbevölkerung", type = "MCMC-Fit 1980-2017", data = "MCMC", year = "1980-2017", Tafel = "MCMC-Zerlegung Bevölkerungssterblichkeit")
) )
) ) %>%
mT.fitExtrapolationLaw(law = "HP2", method = "LF2", fit = 80:98, extrapolate = 90:120, fadeIn = 90:99) %>%
mT.extrapolateTrendExp(idx = 94, up = TRUE)
mort.AT.MCMC[["u"]] = mortalityTable.trendProjection( mort.AT.MCMC[["u"]] = mortalityTable.trendProjection(
name = "Österreich MCMC Unisex", name = "Österreich MCMC Unisex",
ages = as.integer(dimnames(data.array)[[1]]), ages = as.integer(dimnames(data.array)[[1]]),
...@@ -60,11 +66,15 @@ mort.AT.MCMC.load = function() { ...@@ -60,11 +66,15 @@ mort.AT.MCMC.load = function() {
data = list( data = list(
dim = list(sex = "u", collar = "Gesamtbevölkerung", type = "MCMC-Fit 1980-2017", data = "MCMC", year = "1980-2017", Tafel = "MCMC-Zerlegung Bevölkerungssterblichkeit") dim = list(sex = "u", collar = "Gesamtbevölkerung", type = "MCMC-Fit 1980-2017", data = "MCMC", year = "1980-2017", Tafel = "MCMC-Zerlegung Bevölkerungssterblichkeit")
) )
) ) %>%
mT.fitExtrapolationLaw(law = "HP2", method = "LF2",fit = 80:98, extrapolate = 90:120, fadeIn = 90:99) %>%
mT.extrapolateTrendExp(idx = 98, up = TRUE)
mort.AT.MCMC mort.AT.MCMC
} }
mort.AT.MCMC = mort.AT.MCMC.load() mort.AT.MCMC = mort.AT.MCMC.load()
...@@ -73,5 +83,5 @@ rm(mort.AT.MCMC.load) ...@@ -73,5 +83,5 @@ rm(mort.AT.MCMC.load)
############################################################################### ###############################################################################
# mortalityTables.load("Austria*") # mortalityTables.load("Austria*")
# plotMortalityTables(mort.AT.MCMC, mort.AT.forecast, mort.AT.census[,"2011"], Period = 2011) + facet_grid(.~sex) # plotMortalityTables(mort.AT.MCMC, Period = 2008)
# plotMortalityTrend(mort.AT.MCMC, mort.AT.forecast, mort.AT.census[,"2011"], Period = 2011) + facet_grid(.~sex) # plotMortalityTrend(mort.AT.MCMC, Period = 2008)
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment