stopifnot(require(methods), require(utils), require(MortalityTables))
###############################################################################
# DAV 1994R exact (Male, Female), 1st-order only
###############################################################################
DAV1994R.exakt.data = utils::read.csv(
system.file("extdata",
"Germany_Annuities_DAV1994R.csv",
package = "MortalityTables"),
col.names = c("age",
"qx2000", "qy2000",
"trendM", "trendF",
"qxAVbase1955", "qyAVbase1955",
"", "", "", "", "", "", ""),
row.names = NULL,
skip = 2)
DAV1994R.av.data = utils::read.csv(
system.file("extdata",
"Germany_Annuities_DAV1994R_AV.csv",
package = "MortalityTables"),
row.names = 1,
skip = 4
)
DAV1994R.male = mortalityTable.trendProjection(
name = "DAV 1994R male",
ages = DAV1994R.exakt.data$age,
baseYear = 2000,
deathProbs = DAV1994R.exakt.data$qx2000,
trend = DAV1994R.exakt.data$trendM,
data = list(
dim = list(sex = "m", collar = "Rententafel", type = "Rententafel Deutschland", data = "loaded", year = "DAV 1994R")
)
);
DAV1994R.female = mortalityTable.trendProjection(
name = "DAV 1994R female",
ages = DAV1994R.exakt.data$age,
baseYear = 2000,
deathProbs = DAV1994R.exakt.data$qxAVbase1955,
trend = DAV1994R.exakt.data$trendF,
data = list(
dim = list(sex = "w", collar = "Rententafel", type = "Rententafel Deutschland", data = "loaded", year = "DAV 1994R")
)
);
DAV1994R.male.av = mortalityTable.ageShift(
name = "DAV 1994R male, age-shifted",
ages = DAV1994R.exakt.data$age,
deathProbs = DAV1994R.exakt.data$qxAVbase1955,
ageShifts = DAV1994R.av.data["AS.Male"],
data = list(
dim = list(sex = "m", collar = "Rententafel", type = "Rententafel Deutschland", data = "age-shifted", year = "DAV 1994R")
)
);
DAV1994R.female.av = mortalityTable.ageShift(
name = "DAV 1994R female, age-shifted",
ages = DAV1994R.exakt.data$age,
deathProbs = DAV1994R.exakt.data$qyAVbase1955,
ageShifts = DAV1994R.av.data["AS.Female"],
data = list(
dim = list(sex = "w", collar = "Rententafel", type = "Rententafel Deutschland", data = "age-shifted", year = "DAV 1994R")
)
);
rm(DAV1994R.exakt.data, DAV1994R.av.data)
# plot(DAV1994R.male, DAV1994R.male.av, DAV1994T.male, DAV1994T.male.2Ord, DAV2008T.male)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.