Provide USA 1971 IAM individual annuity mortality table as data set

TODO/MortalityTables_USA_Annuities_1971IAM.R → data-raw/US/create_MortalityTables_US_Annuities.R

-stopifnot(require(methods), require(utils), require(MortalityTables))
+################################################################################m#
+#  Skript to generate US Basic (select and ultimate) valuation table objects  ####
+################################################################################m#
+library(MortalityTables)
+library(here)
+library(readxl)
+library(dplyr)
+# library(purrr)
+
+US.annuities.file = here::here("data-raw", "US", "Annuities", "USA_Annuities.xlsx")
+
 
 
 ###############################################################################
 ### 1971 IAM individual annuity table; with optional projection scale B_x
 ###############################################################################
 
-USA1971IAM.data = utils::read.csv(
-    system.file("extdata",
-                "USA_Annuities_1971IAM.csv",
-                package = "MortalityTables"),
-    col.names = c("age","qx", "qy", "B"),
-    skip = 3)
+USA1971IAM.data = read_excel(
+  US.annuities.file, sheet = "1971 IAM",
+  col_names = c("age","qx", "qy", "B"),
+  skip = 4)
 
 USA1971IAM.male = mortalityTable.period(
   name = "USA 1971 IAM, male",
   ages = USA1971IAM.data$age,
   deathProbs = USA1971IAM.data$qx,
   data = list(
-      dim = list(sex = "m", collar = "Mortality", type = "Rententafel", data = "unloaded", year = "1971 IAM")
+    dim = list(table = "1971 IAM", sex = "m", collar = "Mortality", type = "Rententafel", data = "unloaded", year = 1971)
   )
 )
 
@@ -26,7 +34,7 @@ USA1971IAM.female = mortalityTable.period(
   ages = USA1971IAM.data$age,
   deathProbs = USA1971IAM.data$qy,
   data = list(
-      dim = list(sex = "w", collar = "Mortality", type = "Rententafel", data = "unloaded", year = "1971 IAM")
+    dim = list(table = "1971 IAM", sex = "w", collar = "Mortality", type = "Rententafel", data = "unloaded", year = 1971)
   )
 )
 
@@ -37,7 +45,7 @@ USA1971IAM.male.projected = mortalityTable.improvementFactors(
   deathProbs = USA1971IAM.data$qx,
   improvement = USA1971IAM.data$B,
   data = list(
-      dim = list(sex = "m", collar = "Mortality", type = "Rententafel", data = "unloaded", year = "1971 IAM projected")
+    dim = list(table = "1971 IAM projected", sex = "m", collar = "Mortality", type = "Rententafel", data = "unloaded", year = 1971)
   )
 )
 
@@ -48,11 +56,21 @@ USA1971IAM.female.projected = mortalityTable.improvementFactors(
   deathProbs = USA1971IAM.data$qy,
   improvement = USA1971IAM.data$B,
   data = list(
-      dim = list(sex = "w", collar = "Mortality", type = "Rententafel", data = "unloaded", year = "1971 IAM projected")
+    dim = list(table = "1971 IAM projected", sex = "w", collar = "Mortality", type = "Rententafel", data = "unloaded", year = 1971)
   )
 )
 
-rm(USA1971IAM.data)
+USA1971IAM = array(
+  data = c(
+    USA1971IAM.male, USA1971IAM.female, USA1971IAM.male.projected, USA1971IAM.female.projected
+  ),
+  dim = c(2, 2),
+  dimnames = list(Sex = c("m", "f"), type = c("Period", "Generational"))
+)
+
+save(USA1971IAM, USA1971IAM.male, USA1971IAM.female, USA1971IAM.male.projected, USA1971IAM.female.projected,
+     file = here::here("data", "US1971IAM.RData")
+)
+
+
 
-# plot(USA1971IAM.male, USA1971IAM.male.projected, YOB = 1971)
-# plot(USA1971IAM.male, USA1971IAM.male.projected, Period = 1999)