stopifnot(require(methods), require(utils), require(MortalityTables))
###############################################################################
### 1983 Table "a" (individual) and GAM (group annuities), period tables
###############################################################################
USA1983a.data = utils::read.csv(
system.file("extdata",
"USA_Annuities_1983a_GAM.csv",
package = "MortalityTables"),
col.names = c("age","qx", "qy", "qxG", "qyG"),
skip = 3)
USA1983a.male = mortalityTable.period(
name = "USA 1983 Table a, male",
ages = USA1983a.data$age,
deathProbs = USA1983a.data$qx,
data = list(
dim = list(sex = "m", collar = "Mortality", type = "Rententafel", data = "loaded", year = "1983a")
)
)
USA1983a.female = mortalityTable.period(
name = "USA 1983 Table a, female",
ages = USA1983a.data$age,
deathProbs = USA1983a.data$qy,
data = list(
dim = list(sex = "w", collar = "Mortality", type = "Rententafel", data = "loaded", year = "1983a")
)
)
USA1983GAM.male = mortalityTable.period(
name = "USA 1983 GAM, male",
ages = USA1983a.data$age,
deathProbs = USA1983a.data$qxG,
data = list(
dim = list(sex = "m", collar = "group Mortality", type = "Rententafel", data = "loaded", year = "1983 GAM")
)
)
USA1983GAM.female = mortalityTable.period(
name = "USA 1983 GAM, female",
ages = USA1983a.data$age,
deathProbs = USA1983a.data$qyG,
data = list(
dim = list(sex = "w", collar = "group Mortality", type = "Rententafel", data = "loaded", year = "1983 GAM")
)
)
rm(USA1983a.data)
# plot(USA1971IAM.male, USA1971IAM.male.projected, YOB = 1971)
# plot(USA1971IAM.male, USA1971IAM.male.projected, Period = 1999)
# plot(USA1971IAM.male, USA1983a.male, USA1983GAM.male, YOB = 1971)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.