diff --git a/TODO/MortalityTables_USA_Annuities_1971IAM.R b/TODO/MortalityTables_USA_Annuities_1971IAM.R
deleted file mode 100644
index ba9c7fe40b555b32f04ac0a0e74c822fd2ec5674..0000000000000000000000000000000000000000
--- a/TODO/MortalityTables_USA_Annuities_1971IAM.R
+++ /dev/null
@@ -1,58 +0,0 @@
-stopifnot(require(methods), require(utils), require(MortalityTables))
-
-
-###############################################################################
-### 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.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")
-  )
-)
-
-USA1971IAM.female = mortalityTable.period(
-  name = "USA 1971 IAM, female",
-  ages = USA1971IAM.data$age,
-  deathProbs = USA1971IAM.data$qy,
-  data = list(
-      dim = list(sex = "w", collar = "Mortality", type = "Rententafel", data = "unloaded", year = "1971 IAM")
-  )
-)
-
-USA1971IAM.male.projected = mortalityTable.improvementFactors(
-  name = "USA 1971 IAM, male, projected",
-  ages = USA1971IAM.data$age,
-  baseYear = 1971,
-  deathProbs = USA1971IAM.data$qx,
-  improvement = USA1971IAM.data$B,
-  data = list(
-      dim = list(sex = "m", collar = "Mortality", type = "Rententafel", data = "unloaded", year = "1971 IAM projected")
-  )
-)
-
-USA1971IAM.female.projected = mortalityTable.improvementFactors(
-  name = "USA 1971 IAM, female, projected",
-  ages = USA1971IAM.data$age,
-  baseYear = 1971,
-  deathProbs = USA1971IAM.data$qy,
-  improvement = USA1971IAM.data$B,
-  data = list(
-      dim = list(sex = "w", collar = "Mortality", type = "Rententafel", data = "unloaded", year = "1971 IAM projected")
-  )
-)
-
-rm(USA1971IAM.data)
-
-# plot(USA1971IAM.male, USA1971IAM.male.projected, YOB = 1971)
-# plot(USA1971IAM.male, USA1971IAM.male.projected, Period = 1999)
diff --git a/TODO/USA_Annuities.xls b/TODO/USA_Annuities.xls
deleted file mode 100644
index 30ca63d2dfce66ad899391dd3152c8e14e475cae..0000000000000000000000000000000000000000
Binary files a/TODO/USA_Annuities.xls and /dev/null differ
diff --git a/data-raw/US/Annuities/USA_Annuities.xlsx b/data-raw/US/Annuities/USA_Annuities.xlsx
new file mode 100644
index 0000000000000000000000000000000000000000..9921e1512a56cf6b053b99a79918355223334bee
Binary files /dev/null and b/data-raw/US/Annuities/USA_Annuities.xlsx differ
diff --git a/data-raw/US/create_MortalityTables_US_Annuities.R b/data-raw/US/create_MortalityTables_US_Annuities.R
new file mode 100644
index 0000000000000000000000000000000000000000..f431550f7c963f7c1442ce34452e15684d60c927
--- /dev/null
+++ b/data-raw/US/create_MortalityTables_US_Annuities.R
@@ -0,0 +1,76 @@
+################################################################################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 = 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(table = "1971 IAM", sex = "m", collar = "Mortality", type = "Rententafel", data = "unloaded", year = 1971)
+  )
+)
+
+USA1971IAM.female = mortalityTable.period(
+  name = "USA 1971 IAM, female",
+  ages = USA1971IAM.data$age,
+  deathProbs = USA1971IAM.data$qy,
+  data = list(
+    dim = list(table = "1971 IAM", sex = "w", collar = "Mortality", type = "Rententafel", data = "unloaded", year = 1971)
+  )
+)
+
+USA1971IAM.male.projected = mortalityTable.improvementFactors(
+  name = "USA 1971 IAM, male, projected",
+  ages = USA1971IAM.data$age,
+  baseYear = 1971,
+  deathProbs = USA1971IAM.data$qx,
+  improvement = USA1971IAM.data$B,
+  data = list(
+    dim = list(table = "1971 IAM projected", sex = "m", collar = "Mortality", type = "Rententafel", data = "unloaded", year = 1971)
+  )
+)
+
+USA1971IAM.female.projected = mortalityTable.improvementFactors(
+  name = "USA 1971 IAM, female, projected",
+  ages = USA1971IAM.data$age,
+  baseYear = 1971,
+  deathProbs = USA1971IAM.data$qy,
+  improvement = USA1971IAM.data$B,
+  data = list(
+    dim = list(table = "1971 IAM projected", sex = "w", collar = "Mortality", type = "Rententafel", data = "unloaded", year = 1971)
+  )
+)
+
+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")
+)
+
+
+
diff --git a/data/US1971IAM.RData b/data/US1971IAM.RData
new file mode 100644
index 0000000000000000000000000000000000000000..c78d178a2a4583fdda36b280a39ceb9ea489ecea
Binary files /dev/null and b/data/US1971IAM.RData differ