From 3f2326bb815537e287e73b2f2fdc2d53267799aa Mon Sep 17 00:00:00 2001
From: Reinhold Kainhofer <reinhold@kainhofer.com>
Date: Sun, 21 Feb 2021 11:28:23 +0100
Subject: [PATCH] Provide USA 1983a and 1983GAM annuity tables as data sets

---
 TODO/MortalityTables_USA_Annuities_1983a.R    |  57 ---------
 TODO/USA_Annuities_1971IAM.csv                | 115 ------------------
 TODO/USA_Annuities_1983a_GAM.csv              | 115 ------------------
 .../US/create_MortalityTables_US_Annuities.R  |  76 ++++++++++++
 data/USA1983GAM.RData                         | Bin 0 -> 2582 bytes
 data/USA1983a.RData                           | Bin 0 -> 2680 bytes
 6 files changed, 76 insertions(+), 287 deletions(-)
 delete mode 100644 TODO/MortalityTables_USA_Annuities_1983a.R
 delete mode 100644 TODO/USA_Annuities_1971IAM.csv
 delete mode 100644 TODO/USA_Annuities_1983a_GAM.csv
 create mode 100644 data/USA1983GAM.RData
 create mode 100644 data/USA1983a.RData

diff --git a/TODO/MortalityTables_USA_Annuities_1983a.R b/TODO/MortalityTables_USA_Annuities_1983a.R
deleted file mode 100644
index 0fa20c3..0000000
--- a/TODO/MortalityTables_USA_Annuities_1983a.R
+++ /dev/null
@@ -1,57 +0,0 @@
-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)
-
diff --git a/TODO/USA_Annuities_1971IAM.csv b/TODO/USA_Annuities_1971IAM.csv
deleted file mode 100644
index 66506c1..0000000
--- a/TODO/USA_Annuities_1971IAM.csv
+++ /dev/null
@@ -1,115 +0,0 @@
-USA 1971 IAM,,,
-,,,
-,1971 IAM,,Projection Scale B
-Age,Males,Females,
-5,0.000456,0.000234,0
-6,0.000424,0.000193,0
-7,0.000403,0.000162,0
-8,0.000392,0.000143,0
-9,0.000389,0.000134,0
-10,0.00039,0.000132,0.0125
-11,0.000397,0.000143,0.0125
-12,0.000405,0.000155,0.0125
-13,0.000413,0.000167,0.0125
-14,0.000422,0.00018,0.0125
-15,0.000433,0.000193,0.0125
-16,0.000444,0.000205,0.0125
-17,0.000457,0.000218,0.0125
-18,0.000471,0.000231,0.0125
-19,0.000486,0.000245,0.0125
-20,0.000503,0.00026,0.0125
-21,0.000522,0.000275,0.0125
-22,0.000544,0.000292,0.0125
-23,0.000566,0.000309,0.0125
-24,0.000591,0.000327,0.0125
-25,0.000619,0.000347,0.0125
-26,0.00065,0.000368,0.0125
-27,0.000684,0.00039,0.0125
-28,0.000722,0.000414,0.0125
-29,0.000763,0.00044,0.0125
-30,0.000809,0.000469,0.0125
-31,0.00086,0.000499,0.0125
-32,0.000916,0.000533,0.0125
-33,0.000978,0.000569,0.0125
-34,0.001046,0.000608,0.0125
-35,0.001122,0.000651,0.0125
-36,0.001204,0.000698,0.0125
-37,0.001295,0.00075,0.0125
-38,0.001397,0.000807,0.0125
-39,0.001509,0.000869,0.0125
-40,0.001633,0.000938,0.0125
-41,0.001789,0.001013,0.0125
-42,0.002,0.001094,0.0125
-43,0.00226,0.001186,0.0125
-44,0.002569,0.001286,0.0125
-45,0.002922,0.001397,0.0125
-46,0.003318,0.001519,0.0125
-47,0.003754,0.001654,0.0125
-48,0.004228,0.001802,0.0125
-49,0.00474,0.001967,0.0125
-50,0.005285,0.002151,0.0125
-51,0.00586,0.002371,0.0125
-52,0.006461,0.002641,0.0125
-53,0.007088,0.002966,0.0125
-54,0.00774,0.003351,0.0125
-55,0.008417,0.003791,0.0125
-56,0.009119,0.004284,0.01245
-57,0.00985,0.004826,0.0124
-58,0.010613,0.005409,0.01235
-59,0.011411,0.006017,0.0123
-60,0.012249,0.006628,0.01225
-61,0.013133,0.007219,0.0122
-62,0.014073,0.007773,0.01215
-63,0.015083,0.008285,0.0121
-64,0.016185,0.008775,0.01205
-65,0.017405,0.00929,0.012
-66,0.018767,0.009888,0.0118
-67,0.02029,0.010622,0.0116
-68,0.021992,0.011536,0.0114
-69,0.02389,0.012664,0.0112
-70,0.026,0.014029,0.011
-71,0.028341,0.015651,0.0107
-72,0.030933,0.017548,0.0104
-73,0.033801,0.019742,0.0101
-74,0.036976,0.022256,0.0098
-75,0.040494,0.02512,0.0095
-76,0.044393,0.028369,0.0091
-77,0.048715,0.03205,0.0087
-78,0.0535,0.036225,0.0083
-79,0.058787,0.040975,0.0079
-80,0.064599,0.046386,0.0075
-81,0.070902,0.052513,0.007
-82,0.077668,0.059409,0.0065
-83,0.084941,0.06716,0.006
-84,0.092874,0.075899,0.0055
-85,0.101689,0.08577,0.005
-86,0.111652,0.096898,0.0045
-87,0.123048,0.109338,0.004
-88,0.136123,0.122978,0.0035
-89,0.15107,0.137508,0.003
-90,0.16804,0.152472,0.0025
-91,0.187147,0.16737,0.002
-92,0.208457,0.181776,0.0015
-93,0.231885,0.195386,0.001
-94,0.257146,0.208071,0.0005
-95,0.283841,0.219896,0
-96,0.311565,0.231097,0
-97,0.340214,0.242211,0
-98,0.369769,0.253823,0
-99,0.400194,0.266452,0
-100,0.431413,0.280535,0
-101,0.463312,0.296449,0
-102,0.495756,0.314535,0
-103,0.528599,0.335121,0
-104,0.561692,0.358537,0
-105,0.594884,0.385122,0
-106,0.628022,0.415238,0
-107,0.660949,0.449274,0
-108,0.693503,0.487649,0
-109,0.725521,0.530787,0
-110,0.756852,0.57904,0
-111,0.78739,0.632529,0
-112,0.817125,0.690903,0
-113,0.846198,0.753081,0
-114,0.874915,0.817218,0
-115,1,1,0
diff --git a/TODO/USA_Annuities_1983a_GAM.csv b/TODO/USA_Annuities_1983a_GAM.csv
deleted file mode 100644
index eb7c1dc..0000000
--- a/TODO/USA_Annuities_1983a_GAM.csv
+++ /dev/null
@@ -1,115 +0,0 @@
-USA 1983 Table a / GAM,,,,
-,,,,
-,1983 Table a (indiv.),,1983 GAM (Group Annuity Mortality),
-Age,Males,Females,Males,Females
-5,0.000377,0.000194,0.000342,0.000171
-6,0.00035,0.00016,0.000318,0.00014
-7,0.000333,0.000134,0.000302,0.000118
-8,0.000352,0.000134,0.000294,0.000104
-9,0.000368,0.000136,0.000292,0.000097
-10,0.000382,0.000141,0.000293,0.000096
-11,0.000394,0.000147,0.000298,0.000104
-12,0.000405,0.000155,0.000304,0.000113
-13,0.000415,0.000165,0.00031,0.000122
-14,0.000425,0.000175,0.000317,0.000131
-15,0.000435,0.000188,0.000325,0.00014
-16,0.000446,0.000201,0.000333,0.000149
-17,0.000458,0.000214,0.000343,0.000159
-18,0.000472,0.000229,0.000353,0.000168
-19,0.000488,0.000244,0.000365,0.000179
-20,0.000505,0.00026,0.000377,0.000189
-21,0.000525,0.000276,0.000392,0.000201
-22,0.000546,0.000293,0.000408,0.000212
-23,0.00057,0.000311,0.000424,0.000225
-24,0.000596,0.00033,0.000444,0.000239
-25,0.000622,0.000349,0.000464,0.000253
-26,0.00065,0.000368,0.000488,0.000268
-27,0.000677,0.000387,0.000513,0.000284
-28,0.000704,0.000405,0.000542,0.000302
-29,0.000731,0.000423,0.000572,0.00032
-30,0.000759,0.000441,0.000607,0.000342
-31,0.000786,0.00046,0.000645,0.000364
-32,0.000814,0.000479,0.000687,0.000388
-33,0.000843,0.000499,0.000734,0.000414
-34,0.000876,0.000521,0.000785,0.000443
-35,0.000917,0.000545,0.00086,0.000476
-36,0.000968,0.000574,0.000907,0.000502
-37,0.001032,0.000607,0.000966,0.000536
-38,0.001114,0.000646,0.001039,0.000573
-39,0.001206,0.000691,0.001128,0.000617
-40,0.001341,0.000742,0.001238,0.000665
-41,0.001492,0.000801,0.00137,0.000716
-42,0.001673,0.000867,0.001527,0.000775
-43,0.001886,0.000942,0.001715,0.000842
-44,0.002129,0.001026,0.001932,0.000919
-45,0.002399,0.001122,0.002183,0.00101
-46,0.002693,0.001231,0.002471,0.001117
-47,0.003009,0.001356,0.00279,0.001237
-48,0.003343,0.001499,0.003138,0.001366
-49,0.003694,0.001657,0.003513,0.001505
-50,0.004057,0.00183,0.003909,0.001647
-51,0.004431,0.002016,0.004324,0.001793
-52,0.004812,0.002215,0.004755,0.001949
-53,0.005198,0.002426,0.0052,0.00212
-54,0.005591,0.00265,0.00566,0.002315
-55,0.005994,0.002891,0.006131,0.002541
-56,0.006409,0.003151,0.006618,0.002803
-57,0.006839,0.003432,0.007139,0.003103
-58,0.00729,0.003739,0.007719,0.003443
-59,0.007782,0.004081,0.008384,0.003821
-60,0.008338,0.004467,0.009158,0.004241
-61,0.008983,0.004908,0.010064,0.004703
-62,0.00974,0.005413,0.011133,0.00521
-63,0.01063,0.00599,0.012391,0.005769
-64,0.011664,0.006633,0.013868,0.006386
-65,0.012851,0.007336,0.015592,0.007064
-66,0.014199,0.00809,0.017579,0.007817
-67,0.015717,0.008888,0.019804,0.008681
-68,0.017414,0.009731,0.022229,0.009702
-69,0.019296,0.010653,0.024817,0.010922
-70,0.021371,0.011697,0.02753,0.012385
-71,0.023647,0.012905,0.030354,0.014128
-72,0.026131,0.014319,0.03337,0.01616
-73,0.028835,0.01598,0.03668,0.018481
-74,0.031794,0.017909,0.040388,0.021092
-75,0.035046,0.020127,0.044597,0.023992
-76,0.038631,0.022654,0.049388,0.027185
-77,0.042587,0.025509,0.054758,0.030672
-78,0.046951,0.028717,0.060678,0.034459
-79,0.051755,0.032328,0.067125,0.038549
-80,0.057026,0.036395,0.07407,0.042945
-81,0.062791,0.040975,0.081484,0.047655
-82,0.069081,0.046121,0.08932,0.052691
-83,0.075908,0.051889,0.097525,0.058071
-84,0.08323,0.058336,0.106047,0.063807
-85,0.090987,0.065518,0.114836,0.069918
-86,0.099122,0.073493,0.12417,0.07657
-87,0.107577,0.082318,0.13387,0.08387
-88,0.116316,0.092017,0.144073,0.091935
-89,0.125394,0.102491,0.154859,0.101354
-90,0.134887,0.113605,0.166307,0.11175
-91,0.144873,0.125227,0.178214,0.123076
-92,0.155429,0.137222,0.19046,0.13563
-93,0.166629,0.146462,0.203007,0.149577
-94,0.178537,0.161834,0.217904,0.165103
-95,0.191214,0.174228,0.234086,0.182419
-96,0.204721,0.186535,0.248436,0.201757
-97,0.21912,0.198646,0.263954,0.222044
-98,0.234735,0.211102,0.280803,0.243899
-99,0.251889,0.224445,0.299154,0.268185
-100,0.270906,0.239215,0.319185,0.295187
-101,0.292111,0.255953,0.341086,0.325225
-102,0.315826,0.275201,0.365052,0.358897
-103,0.342377,0.2975,0.393102,0.395843
-104,0.372086,0.32339,0.427255,0.43836
-105,0.405278,0.353414,0.469531,0.487816
-106,0.442277,0.388111,0.521945,0.545886
-107,0.483406,0.428023,0.586518,0.614309
-108,0.528989,0.473692,0.665268,0.694855
-109,0.579351,0.525658,0.760215,0.789474
-110,0.634814,0.584462,1,1
-111,0.695704,0.650646,,
-112,0.762343,0.72475,,
-113,0.835056,0.807316,,
-114,0.914167,0.898885,,
-115,1,1,,
diff --git a/data-raw/US/create_MortalityTables_US_Annuities.R b/data-raw/US/create_MortalityTables_US_Annuities.R
index f431550..a4a91db 100644
--- a/data-raw/US/create_MortalityTables_US_Annuities.R
+++ b/data-raw/US/create_MortalityTables_US_Annuities.R
@@ -74,3 +74,79 @@ save(USA1971IAM, USA1971IAM.male, USA1971IAM.female, USA1971IAM.male.projected,
 
 
 
+###############################################################################
+### 1983 Table "a" (individual) and GAM (group annuities), period tables
+###############################################################################
+
+USA1983a.data = read_excel(
+  US.annuities.file, sheet = "1983a - GAM",
+  col_names = c("age","qx", "qy", "qxG", "qyG"),
+  skip = 4)
+
+USA1983a.male = mortalityTable.period(
+  name = "USA 1983 Table a, male",
+  ages = USA1983a.data$age,
+  deathProbs = USA1983a.data$qx,
+  baseYear = 1983,
+  data = list(
+    dim = list(table = "1983 table a", sex = "m", collar = "Mortality", type = "Rententafel", data = "loaded", year = 1983)
+  )
+)
+
+USA1983a.female = mortalityTable.period(
+  name = "USA 1983 Table a, female",
+  ages = USA1983a.data$age,
+  deathProbs = USA1983a.data$qy,
+  baseYear = 1983,
+  data = list(
+    dim = list(table = "1983 table a", sex = "w", collar = "Mortality", type = "Rententafel", data = "loaded", year = 1983)
+  )
+)
+
+
+USA1983GAM.male = mortalityTable.period(
+  name = "USA 1983 GAM, male",
+  ages = USA1983a.data$age,
+  deathProbs = USA1983a.data$qxG,
+  baseYear = 1983,
+  data = list(
+    dim = list(table = "1983 GAM", sex = "m", collar = "group Mortality", type = "Rententafel", data = "loaded", year = 1983)
+  )
+)
+
+USA1983GAM.female = mortalityTable.period(
+  name = "USA 1983 GAM, female",
+  ages = USA1983a.data$age,
+  deathProbs = USA1983a.data$qyG,
+  baseYear = 1983,
+  data = list(
+    dim = list(table = "1983 GAM", sex = "w", collar = "group Mortality", type = "Rententafel", data = "loaded", year = 1983)
+  )
+)
+
+USA1983a = array(
+  data = c(
+    USA1983a.male, USA1983a.female
+  ),
+  dim = c(2),
+  dimnames = list(Sex = c("m", "f"))
+)
+
+USA1983GAM = array(
+  data = c(
+    USA1983GAM.male, USA1983GAM.female
+  ),
+  dim = c(2),
+  dimnames = list(Sex = c("m", "f"))
+)
+
+save(USA1983a, USA1983a.male, USA1983a.female,
+     file = here::here("data", "USA1983a.RData")
+)
+save(USA1983GAM, USA1983GAM.male, USA1983GAM.female,
+     file = here::here("data", "USA1983GAM.RData")
+)
+
+
+
+
diff --git a/data/USA1983GAM.RData b/data/USA1983GAM.RData
new file mode 100644
index 0000000000000000000000000000000000000000..423cc0c08b31fb8d111d85424cec3256dd2573fa
GIT binary patch
literal 2582
zcmb2|=3oE==C?CALxPXX9dl;V(CA>jBkmFq!lJ67(P)1{N-WaJKw#U1I}Tn>6I@)^
z^2hY86`Q(ifySvt5w&>=dJ%<=tS&JVj;_73YR7BmDeCVoasOeS^y;<cyz6%U@p>)#
z_y0UQx7WB>y=EqZIA84S#l8yzGDJ#0wJ&qrw=`n|^W9^&f;3)*Tv{>lXw;Id&5?eO
z)o)(9JGVwdWR+c{)6x|ee`#M0dhRN06jXY0`jv|fm&z=@Im@<WwPf{VHCZuT<ta5;
zWx0y=mBp2hE5cTkgk)dgTp=8o8;~2QzR=HYt8?mtTZ^q0Fs)kbH2L;b&b`9FG+!)V
zQ114==N{``?yB%F>t5WwuzO*-`}@9o%zN3Z)T-2~{J+e5vG+pt!tZYHd+)LD<^Q$v
z#qvepUEg=#<KH`9dDWsn6}x{ae|LY^e@DKO;S<9X#uLmFHbWu162FUFhmC-xGz3-(
zd{TJgc%pfNzq1|hPm?E?Clo^=yE4C<T#t=JCC4X$CyFPWClpVFz;-2mSGg{mEpxTk
z$1on8eEf&1$6~86H6M|=>vdinWPN>F`SRxbl1zVn9$whu`^}-;v(GDc1<S6F#xvhH
zJ$G0x|7+L9ZE`Q}r&+w5uE+K(+5E2dt-l9uz4wf5J72azulGz~v3WD+z1K=6v%?R!
ztUBdryM5KcBbM2Io4afe9GjsXA3vSTtpD2!-skUGZ~wg_Yb|%}(2_l}Vzpti?I+#Y
za-_fM$i6u>U8?`9Kf~{TZN+@`%h~u0Bd?xWypkc`W{dN=YfX<?-z_*Q*L~on%e8{V
zTU=%oT~YpCEznT>cWvn1Z#>N5kJtL1*Y$6>-4w2;*34j2BB9&FeY&C2qx|>Xw&n-9
zj`g2PczJ8M(%)|Q#KCZ1o-<oog8!eCdH3II^V%AAK1;p3xZHVxyI0E5&3gr=Bs|?4
z@Phx~D^~UQmHV16{41NB{OP(g%e>>eQ$42JcriCN&r93><h#P3$md;8pB4zbU0wZU
z-weh_M(bNO_p&?Le4Tc+%O`wBox^X(o^@&`()PzyU(T{WbKv8w)QdZ2y9OkFThN_)
z{jb=eJB^#yXW7<reSYIFeoSg#(~{ovwsP?+n|~&6di<<Mo{PIn)pvcx$qgFLRv$f!
z1a$1GpSLdlTN)vxx!W$~O|6Gu*y(S3E;fkrR%;i>{MazbXs1Eedz*;2^S)g9FUNh!
zFUs^A@0sS3z-#N4?roi5xGu_A-CoX!xBbkF70#_WhvPZRI$57g5nXoenOu|Z?yeTQ
zGT#kb&8A%{IxD!%c=3#zUxYUVpSJ7We?scpmu+jz`nwuxPBYGQPoH-9UYh&AV+H;<
z!qfL0-m*Ym=0|G!;ofCj8a0l4qO0DjTjW%^*SvI5mQl=ecdNeBmi+X!{0tl8=A=3M
zg+J!Y<tEKKd~($dmiF213tDC-UB6t)IVH|QqO?`9@$$s?cQ3O~bNCsZJh|Sqo>g@F
zUCFIRjB~R8nbimF)Q_>KzLPe!kvrtxxi8|Xk_V4Y-|e&b`clWOnHP@s_P=-NP5JnE
zU;j6UTwf2Dn|EsqEPIw6*Yjy#otAAgdD7R1>3y2hkF9NeZL8ha&MI=aw<+k<isc&%
zUWXV?w0h6J<M1BCAJbOw-${Sl$7wCW-TU^n?n>X?PybBRTI-{$WEiJ4d*X#}`oV!K
zR+>$H`*!WJkWhcO4$gfsD!s*zKR&(0@`<N+TL1OLW1UZw=H1Yol@a?V_s@5EF1Dh$
zDKlm`aV+Iz@xJ)Q{Z>f-W6!AD7k@8!qiw_`u{DiR*@$O`v}Z+<)5|sQ0&OG(eLe_%
zV0J%Sc=>jJtx2h{cO_F`%$7w`OCP@oTD#muOMuh0ChoWIs>>%fUH2-zB6QL||HSp`
zPpKWfFSl&Fa?2sRf8*0@pS}reTIs$m&B*;aN%NHZbnc7~M{Q>MWNLglA*r(5-|gP?
z9!Yftq1*WjUJGzl6v?HjnTGf-E!(Bt7Q8aXcTva%Z7WAB<$EXGXDK}Ibp3SGYll07
zpt^c}rIXW|7L(}XvjqMe?8(@2Zu_F?v$_s#YiH&4b^9hX`-%6_q@qb1=YP8S(Pg^b
z<KH*m@-y{y$8}%oy~v{EseVn>@T+r<M&K9I4PW(i9(A%z3|a6;;ppcd8~0>9Ro^<@
ze|mB1kBU9-eb!wO>(W}gXi@mm(2ysqwsdv<@%LQ4rLUAj*<iL5d;YGNU(2}KOkzC0
zPTkYmc#iSN2^roO`jHQJukw)o)j!Sm`fk}9fp!*`k1xC7#^h{h%y;<Al8khVv<Trd
zdz4LlR`RAkH7Huv{QTnjmkD7(`)*plR_S?mB%b^BW!Ym1YO>bvXSfF)_1ci4%W?gD
zJm2H6=Q69-X*7OZ8OOt3DHoCYr=%fr?-FI_=>2#8$JU;?+}D#^{I+gC>u#?Ep?ll!
zZU{Dh=O@Ga?e20tnSiDI2LJ9ChOT3MpkZvW+2|3YMdI#YuX4V}>z+4t`*vpVEM|#0
z-?3RD?$(38@26Fch0K5QvibTg&a${^w-aAINbEd&|9W#gm&Jug=ban+uc^gGs<7;8
zd2(g?4F>^zvkzhO7%Xn9T%Nt?hrp+>za@cgye+Pue)mPVwR4>5*|Q<7EaQM}N7g?3
zu9@Bq3pTyW^-PrInEls%ecr+bjk2_wv~}hkYyVFydZ)#_B=_*SwwHoef~Oz4>+nn`
z`&F>LY?_a;MTp?xmdKkt!TxOw&t^o1dz$ZW&Yw5q%dNN^&XXbUUa$IMYOufOa=%Pf
zL`5z4MC;44V(q&#Pfg!_c+!TI+MlZCw2G{m-}N`@XT^8te~qde)b1WxbIgkS+TR0z
zQ`2~u&0@M{m1=F_FiL**Dt%2I<IQNJGwiSSHlKC<nSS}s8HMk;dl#O|W51!6a^gNq
z+ClsH-)p@-^EAe`KR@=Yf02{i@8ep_pGck9pM6R^?p$lr*IL!<=gbfMJMxp6eZ5Xg
zvv-Npj?BdmE*O_h{@*60ace^K`Km(3zO(!1f4;GXq3fih^lV!hzoctAhVPDRa_O~f
zdhoOO`ja(nMMj$K%k3LBZ|wc`<hu5bglvlgmKR0zCNXX<&lC;0`|QFwNuGvTr}D+#
zuGen$)<4eEJ2RUh{O0Gpm`Sf3=H}i1rymz~VEwj<FQv6D3W~4IZMVMSyRa@gXM^J$
z^E;^)deL!mT>Ntm$A5O3Rv%{;n*WlgahCj^_t$pWO+A$Ioqxx)PxYnE>l6b$N{+eY
z9*qj;Qna-?r77#4U3q@mx|!>BxBcY|^bg<3|6lR$*Q}p+T+F*aO;_2~H{<Kk_8s?5
z<<=T)Ej?RP;W_tQ=eu9BEk`8MM1`-GiDsp5Ic9#y|7+*c86WS*?hb#Rx5D*y^IEPc
zs~lExO?mrkDXZ7jwB=cwgSX!NX@Bd@pY^P#{NGvqt3Q5ne@63%_wwr1^#U(H-0sr4
zveId;)|IWVr;EJ`+1;`#*6FI=!mQ~V7XNuK`Z{EH%Y6CQ8t?0wS2s*jyP?6r000t)
B4M6|^

literal 0
HcmV?d00001

diff --git a/data/USA1983a.RData b/data/USA1983a.RData
new file mode 100644
index 0000000000000000000000000000000000000000..c5dd8c607871d1f834ca9709b00db11628b1a8fd
GIT binary patch
literal 2680
zcmb2|=3oE==C?E0Lqfvkj<stkw5YXc2<l8=;qjWd=#cy5d+iH1TOC;Z;($Vro<>td
zXP~1BXJ<wTOOA%y23KJLR?{wz19B!=9Rem=VvF3$*#33&s8yePS7oN2&erz6_IvsL
z+U<Le*Az3)aXvHEX-b67G#>A;AIi$=HYq0ynPTm>q-{)D+Un*Nu|^|=mGAb!E5~ll
zIi2L0&9&FW+v}O^MypMCEYw0*<a6Gy`JK4t-t@hEEA>{|z0A0@aY@^f_mib2_9>ex
z&UMK0%5ux{%W|}ucxBR+t}C;yw3P_GOu3}FWZROmB??>LUuUg-<#4ZYk8-8^C%-3i
zPgqYbSF@k`C*+gflh_lxCzY$dce~fS$GWooliU;A6V;QxtGsu;*S*KQGXIm^liCx%
zCw*6a?|$$6ye<3p?DGHH`CItA`a6G{h6;v{438L(+=N2rPIe)F6*-SN&JbwRQ6cbA
z;gRE!=1z8Del<BCo0bZWj{=Vrk2rU>3-PPUdD-+-NPIMSlz8NDXS*=Jx}2X)Qw7UM
zjz@yEF>hvubS3U+s@N}5ccoih_*&?t+uO^x9Qd^E&#AxcuMIb@3=H2H!WH(4Rnniw
z{9x6RGs!zuKPN2y5d3ZCmXaGLH}^>2*(^4Bn*3_VWM#G4wR3iC31*+H|8xIkgA}P{
zI#<?RmWpM&?Rk{*{H}Ee{7Qb^k-Kqs=Y;38fACbiS|(DqiuIdzi^fhyKCw-*a@t!S
z6kK*u_OcS_dwTP!to+%=%bS$g&+b}U@a$dgwy8gPg7<PQKWVSkaJu5H%C%Kn3zq)d
zX?Nkrj@Pr-7G&2hsI4fMyptWo^KIw*mg9V7uHUA}e)zMg;C{zb=?9^f#mAGMr~J@p
zPc4p;J5>JP;_+wxuMf(zJ>NaPP{DG8JIO+}{D7TMM9lLiOn&9>IK5evQ^Ib)<Tv>$
zFmJ)ahxV3C|1)OW>Uams-FZ@c?E9Bbg0aiqotNDxxp9@W^7k7lO8ZO<_8;HOV!ZfO
zI=5WE<D|#^ZNJ#eJmS^|8ePh<2;#}zoNkiio^^G*aR{GT@;3KtOueh08!mr)+1Imi
z^Mq}eJo2wDkvg-L-$HzMsiIMe&ej>qzBB4mQ>yNq<V%U2eE0i?BQiH<b*JeV{gdII
z`KjC8$o}A>-<7wuYyWI~GQ;zfYA<I*+2PF%YYR9||201$oBCx#QN$CWw+D((Tz(lF
zuqN-{+GW;34W=gd{MYS%#xn0-ztidX8ZOh~4@U*pGOI>+ZOF@fetZ9$pu=q8?-$K1
zO1n{!e{_1a>E6H?`_m`Wf>gcy?6rK1=Qiqg^&jJTcrr2i+A-erk|hCo+dMK~bR1rB
z;fle$?VnA`-hbK=BYwHZ)^qCHo7aj1pH7Z8{Ic8Z&F2>}l8Q5WnJ<Ux{t@!}@U>|1
zp38k@7BS_u+}df<tES#A4S$=k@Yu`Fi{IZaE_jl>xvgKqF0201%U3o}yOJB<yU6gr
zbGoa#K02!^&0^lZ<A-<dV^Y`Mv*wl@|8dRKy^mv)FI`mFGb8W3eum7=rB9N#ilsMC
z>pOFPyYuCDzvt-vS2BrCv-**;ZL#Z(-Cq{Jvbg`Epey%DUw65w`+fhfhehR|bMC(U
zZNb-d^FEy0*^s^cIosA*``cfiZtt`B@<1%_m*L9A%jP=GzU<QKp7(Xf(LLE4)&@27
z{rxH^*S+sh?(fs@J;irF)O%g(zf5WJv@`y5=U)0f@n*b+{?o%BKTQ%^y=vO24>K|x
zne0!tO6xvY&T)Ri>&jZ$?@L%5nWjG%INl-mFYe!~rR|(M0&R{>?0mZ1;LX2f@t2q7
z2$=_7>AU5fdh3wF0(~#0UyY00zU+MbMR=!;n~k1JVerg=w%JbqOP0L3ZeM*@h3DbN
zoravAg@3O`=*Km$(p3MHER-a?XlJeL)Tx5gYz0><PgVQ#Ipjy3UQKpjveu$$F1#IX
zzPjC~e*2}Y{PN}M#FuNgH>m!TDq-B2A!lwHnqzgPWwl1kp^(J=L3bx7&DgS!k7wR7
ziFFI?%{?cc>fXBej`h-}eVU@lL7hpgUHAToFixMY;33ZJ%sfM-@xyfK(;tpZ*psnr
z=8XEyaz%?o-^(wqsZY$`e&WOuuHzoQZl5b{9tvI7V10S=%mTs7g;n=+9X2nzrhLtC
z=8TpjjDMyHc<=AoD)jsKt?K<VRg5(`|44kitua&DcCxg3-TGBk{8r{u^ZK4M-wfT=
zu6lCyk|(N_Q>W}uSwBgrUb!;#+0o(-exc&p>oP%yZoFGHNBZM;1FiGLU2_*CIqdJ3
zj(aTHeev_LtJ+829+7y?8dzB_Wq9rVb#2#s(TN)p@3|~Kdd*--$?t#9PEK{!t!oSP
z__XvyESJRcm{Y5MJgA)!;3ez-@=A!O@#$YPgwJdd|M$!$O#In1_Ov@&_c3yNEf!uG
zdwTVc5chAF9K6I*1uO41mh9PoE#v5|qcbjr?R}QC#`F!3nEi_E-7H6T6l*<xcF=B%
za>bPOGBPI5Y(lThXK6pJZXX*f#(CU$)$1PlUx!~F`)cZ3)_W)<S@qw>r|b_lRnD}W
zY%G&kr)d8t<{1zFjw5C3<eAt1IM;CY)z%N0q55}r#*5aaA1a-`OV`a|iplTI3xqjL
z<{8*6oz$za`pdUt-)dT(SlPwDnBQ5@?EZ}R;?j<qUC$P4)w24;m;CJ5d9Bh#a>lCr
zYri^ux!NVR_f3#N=t|k@M`lWS&-r;Jdxa-nSdlMYyIOQldF9=@@U0WKT(pZw%(0vF
z{pPDDR^JqkelDNfRH4byrXDf3_FO`rhW6)86Q*50s~bAAtMZuG9Zpy737(lX2Vea(
zp1C%9!-YSqF78^vqqD0vas4X}XC2FI?F`Kw{i*z1+1h<=NtKV*e?8jzE8=9C?#dT|
z?r}@9Jm0Okc!J+5BuR857pv=Iw-tZqb)G-drDGIs@_^fxRdVJloi*LRWglKPOFZm#
z{?kuM`J~IsPE38-yY$ANl;cOE%lb5g&u)`vza2KajE{F#@29ZU6F$gqU-0NgXVQA!
ztv^;x{GWX+!a5{#&Yy>?E*HLup1LA*-Lc~bHO0N%+~W@<>aNz`o|d(!%45^h&R1(!
zoSeP?mcj~7_kE?`0&gyBT3lUU9Txak?$w@l<4=dT1lH&ND!zWcv{0$g=zT#V+rG)&
z`_@XS9$4A@`pdHYu^nw+Tf_CD=ej7Hr++s%>T3~QvU$xX_4u0gDYIj0;yBO7pM5iL
z^FFcX*NRuo550H1wK#p+iaXqgpU!AzjlW;wx;lRQ-g6$2G1;>U9-OMEj=6t1?pAIE
zyZP*U=X(FA+0>bt=J;pV^?jQ*$0PW&PRX@A@Bh2@q!n0AeJk0~_(@-uU(j}bQen_T
z)pZ*GR<8S;<<pbO*0d<-|7Sfhm$rQx=^NwTO)onC;L3IBJNF!JRD0LR)K7XJQvGkX
z&v*Z8(&tL_o36#%9b4JG{cHQq+=5q&y;t$*e17m=!O5p^Tm8(c^-8t^t6wU7>pynn
zUHtp-b!#u}m(RLX9qi4tR3*TRDRut8t4vE%b_Z{|t9kR?AN$<$f9IK=&VOg^_2`N8
zeWgFYk9%GfpV+r=rH9j8l`Ec3vGf03b()(pJz$f!#!b~JZ-h3*y}N#2|C(nR=l{tO
R$E*J7oBmq!;+6&j0{~#CS?2%%

literal 0
HcmV?d00001

-- 
GitLab