From e2da404867b592d7c93d6b2dde694ec101fc048e Mon Sep 17 00:00:00 2001
From: Reinhold Kainhofer <reinhold@kainhofer.com>
Date: Sun, 21 Feb 2021 13:27:06 +0100
Subject: [PATCH] Provide USA 1994 GAR/GAM annuity tables as data sets

---
 TODO/MortalityTables_USA_Annuities_1994GAR.R  |  56 --------
 TODO/USA_Annuities_1994GAR.csv                | 124 ------------------
 .../US/create_MortalityTables_US_Annuities.R  |  76 +++++++++++
 data/USA1994GAM.RData                         | Bin 0 -> 2733 bytes
 data/USA1994GAR.RData                         | Bin 0 -> 3059 bytes
 5 files changed, 76 insertions(+), 180 deletions(-)
 delete mode 100644 TODO/MortalityTables_USA_Annuities_1994GAR.R
 delete mode 100644 TODO/USA_Annuities_1994GAR.csv
 create mode 100644 data/USA1994GAM.RData
 create mode 100644 data/USA1994GAR.RData

diff --git a/TODO/MortalityTables_USA_Annuities_1994GAR.R b/TODO/MortalityTables_USA_Annuities_1994GAR.R
deleted file mode 100644
index 5793dd3..0000000
--- a/TODO/MortalityTables_USA_Annuities_1994GAR.R
+++ /dev/null
@@ -1,56 +0,0 @@
-stopifnot(require(methods), require(utils), require(MortalityTables))
-
-###############################################################################
-### 1994 GAR/GAM group annuity tables, with improvement factors AA_x
-###############################################################################
-
-USA1994GAM.data = utils::read.csv(
-    system.file("extdata",
-                "USA_Annuities_1994GAR.csv",
-                package = "MortalityTables"),
-    col.names = c("age","qx", "AAx", "qy", "AAy", "qxBasic", "qyBasic"),
-    skip = 3)
-
-USA1994GAM.male.basic = mortalityTable.period(
-    name = "USA 1994 GAM basic (unloaded), male",
-    ages = USA1994GAM.data$age,
-    deathProbs = USA1994GAM.data$qxBasic,
-    data = list(
-        dim = list(sex = "m", collar = "group Mortality", type = "Rententafel", data = "unloaded", year = "1994 GAM")
-    )
-)
-
-USA1994GAM.female.basic = mortalityTable.period(
-    name = "USA 1994 GAM basic (unloaded), female",
-    ages = USA1994GAM.data$age,
-    deathProbs = USA1994GAM.data$qyBasic,
-    data = list(
-        dim = list(sex = "m", collar = "group Mortality", type = "Rententafel", data = "unloaded", year = "1994 GAM")
-    )
-)
-
-
-USA1994GAR.male = mortalityTable.improvementFactors(
-    name = "USA 1994 GAR, male",
-    ages = USA1994GAM.data$age,
-    deathProbs = USA1994GAM.data$qx,
-    improvement = USA1994GAM.data$AAx,
-    data = list(
-        dim = list(sex = "m", collar = "group Reserving", type = "Rententafel", data = "loaded", year = "1994 GAR")
-    )
-)
-
-USA1994GAR.female = mortalityTable.improvementFactors(
-    name = "USA 1994 GAR, female",
-    ages = USA1994GAM.data$age,
-    deathProbs = USA1994GAM.data$qy,
-    improvement = USA1994GAM.data$AAy,
-    data = list(
-        dim = list(sex = "m", collar = "group Reserving", type = "Rententafel", data = "loaded", year = "1994 GAR")
-    )
-)
-
-rm(USA1994GAM.data)
-
-# plot(USA1994GAM.male.basic, USA1994GAR.male, Period = 1994)
-
diff --git a/TODO/USA_Annuities_1994GAR.csv b/TODO/USA_Annuities_1994GAR.csv
deleted file mode 100644
index 0f6eb80..0000000
--- a/TODO/USA_Annuities_1994GAR.csv
+++ /dev/null
@@ -1,124 +0,0 @@
-USA 1994 GAR,,,,,,
-,,,,,,
-,1994 GAR Male,,1994 GAR Female,,1994 GAM Basic (unloaded),
-Age,qx1994,AAx,qy1994,AAy,Male,Female
-1,0.000592,0.02,0.000531,0.02,0.000637,0.000571
-2,0.0004,0.02,0.000346,0.02,0.00043,0.000372
-3,0.000332,0.02,0.000258,0.02,0.000357,0.000278
-4,0.000259,0.02,0.000194,0.02,0.000278,0.000208
-5,0.000237,0.02,0.000175,0.02,0.000255,0.000188
-6,0.000227,0.02,0.000163,0.02,0.000244,0.000176
-7,0.000217,0.02,0.000153,0.02,0.000234,0.000165
-8,0.000201,0.02,0.000137,0.02,0.000216,0.000147
-9,0.000194,0.02,0.00013,0.02,0.000209,0.00014
-10,0.000197,0.02,0.000131,0.02,0.000212,0.000141
-11,0.000208,0.02,0.000138,0.02,0.000223,0.000148
-12,0.000226,0.02,0.000148,0.02,0.000243,0.000159
-13,0.000255,0.02,0.000164,0.02,0.000275,0.000177
-14,0.000297,0.019,0.000189,0.018,0.00032,0.000203
-15,0.000345,0.019,0.000216,0.016,0.000371,0.000233
-16,0.000391,0.019,0.000242,0.015,0.000421,0.000261
-17,0.00043,0.019,0.000262,0.014,0.000463,0.000281
-18,0.00046,0.019,0.000273,0.014,0.000495,0.000293
-19,0.000484,0.019,0.00028,0.015,0.000521,0.000301
-20,0.000507,0.019,0.000284,0.016,0.000545,0.000305
-21,0.00053,0.018,0.000286,0.017,0.00057,0.000308
-22,0.000556,0.017,0.000289,0.017,0.000598,0.000311
-23,0.000589,0.015,0.000292,0.016,0.000633,0.000313
-24,0.000624,0.013,0.000291,0.015,0.000671,0.000313
-25,0.000661,0.01,0.000291,0.014,0.000711,0.000313
-26,0.000696,0.006,0.000294,0.012,0.000749,0.000316
-27,0.000727,0.005,0.000302,0.012,0.000782,0.000324
-28,0.000754,0.005,0.000314,0.012,0.000811,0.000338
-29,0.000779,0.005,0.000331,0.012,0.000838,0.000356
-30,0.000801,0.005,0.000351,0.01,0.000862,0.000377
-31,0.000821,0.005,0.000373,0.008,0.000883,0.000401
-32,0.000839,0.005,0.000397,0.008,0.000902,0.000427
-33,0.000848,0.005,0.000422,0.009,0.000912,0.000454
-34,0.000849,0.005,0.000449,0.01,0.000913,0.000482
-35,0.000851,0.005,0.000478,0.011,0.000915,0.000514
-36,0.000862,0.005,0.000512,0.012,0.000927,0.00055
-37,0.000891,0.005,0.000551,0.013,0.000958,0.000593
-38,0.000939,0.006,0.000598,0.014,0.00101,0.000643
-39,0.000999,0.007,0.000652,0.015,0.001075,0.000701
-40,0.001072,0.008,0.000709,0.015,0.001153,0.000763
-41,0.001156,0.009,0.000768,0.015,0.001243,0.000826
-42,0.001252,0.01,0.000825,0.015,0.001346,0.000888
-43,0.001352,0.011,0.000877,0.015,0.001454,0.000943
-44,0.001458,0.012,0.000923,0.015,0.001568,0.000992
-45,0.001578,0.013,0.000973,0.016,0.001697,0.001046
-46,0.001722,0.014,0.001033,0.017,0.001852,0.001111
-47,0.001899,0.015,0.001112,0.018,0.002042,0.001196
-48,0.002102,0.016,0.001206,0.018,0.00226,0.001297
-49,0.002326,0.017,0.00131,0.018,0.002501,0.001408
-50,0.002579,0.018,0.001428,0.017,0.002773,0.001536
-51,0.002872,0.019,0.001568,0.016,0.003088,0.001686
-52,0.003213,0.02,0.001734,0.014,0.003455,0.001864
-53,0.003584,0.02,0.001907,0.012,0.003854,0.002051
-54,0.003979,0.02,0.002084,0.01,0.004278,0.002241
-55,0.004425,0.019,0.002294,0.008,0.004758,0.002466
-56,0.004949,0.018,0.002563,0.006,0.005322,0.002755
-57,0.005581,0.017,0.002919,0.005,0.006001,0.003139
-58,0.0063,0.016,0.003359,0.005,0.006774,0.003612
-59,0.00709,0.016,0.003863,0.005,0.007623,0.004154
-60,0.007976,0.016,0.004439,0.005,0.008576,0.004773
-61,0.008986,0.015,0.005093,0.005,0.009663,0.005476
-62,0.010147,0.015,0.005832,0.005,0.010911,0.006271
-63,0.011471,0.014,0.006677,0.005,0.012335,0.007179
-64,0.01294,0.014,0.007621,0.005,0.013914,0.008194
-65,0.014535,0.014,0.008636,0.005,0.015629,0.009286
-66,0.016239,0.013,0.009694,0.005,0.017462,0.010423
-67,0.018034,0.013,0.010764,0.005,0.019391,0.011574
-68,0.019859,0.014,0.011763,0.005,0.021354,0.012648
-69,0.021729,0.014,0.012709,0.005,0.023364,0.013665
-70,0.02373,0.015,0.01373,0.005,0.025516,0.014763
-71,0.025951,0.015,0.014953,0.006,0.027905,0.016079
-72,0.028481,0.015,0.016506,0.006,0.030625,0.017748
-73,0.031201,0.015,0.018344,0.007,0.033549,0.019724
-74,0.034051,0.015,0.020381,0.007,0.036614,0.021915
-75,0.037211,0.014,0.022686,0.008,0.040012,0.024393
-76,0.040858,0.014,0.025325,0.008,0.043933,0.027231
-77,0.045171,0.013,0.028366,0.007,0.04857,0.030501
-78,0.050211,0.012,0.031727,0.007,0.053991,0.034115
-79,0.055861,0.011,0.035362,0.007,0.060066,0.038024
-80,0.062027,0.01,0.039396,0.007,0.066696,0.042361
-81,0.068615,0.009,0.043952,0.007,0.07378,0.04726
-82,0.075532,0.008,0.049153,0.007,0.081217,0.052853
-83,0.08251,0.008,0.054857,0.007,0.088721,0.058986
-84,0.089613,0.007,0.060979,0.007,0.096358,0.065569
-85,0.09724,0.007,0.067738,0.006,0.104559,0.072836
-86,0.105792,0.007,0.075347,0.005,0.113755,0.081018
-87,0.115671,0.006,0.084023,0.004,0.124377,0.090348
-88,0.12698,0.005,0.09382,0.004,0.136537,0.100882
-89,0.139452,0.005,0.104594,0.003,0.149949,0.112467
-90,0.152931,0.004,0.116265,0.003,0.164442,0.125016
-91,0.16726,0.004,0.128751,0.003,0.179849,0.138442
-92,0.182281,0.003,0.141973,0.003,0.196001,0.15266
-93,0.198392,0.003,0.155931,0.002,0.213325,0.167668
-94,0.2157,0.003,0.170677,0.002,0.231936,0.183524
-95,0.233606,0.002,0.186213,0.002,0.251189,0.200229
-96,0.25151,0.002,0.202538,0.002,0.270441,0.217783
-97,0.268815,0.002,0.219655,0.001,0.289048,0.236188
-98,0.285277,0.001,0.237713,0.001,0.30675,0.255605
-99,0.301298,0.001,0.256712,0.001,0.323976,0.276035
-100,0.317238,0.001,0.276427,0.001,0.341116,0.297233
-101,0.333461,0,0.296629,0,0.35856,0.318956
-102,0.35033,0,0.317093,0,0.376699,0.34096
-103,0.368542,0,0.338505,0,0.396884,0.364586
-104,0.387885,0,0.361016,0,0.418855,0.389996
-105,0.407224,0,0.383597,0,0.440585,0.41518
-106,0.425599,0,0.405217,0,0.460043,0.438126
-107,0.441935,0,0.424846,0,0.4752,0.456824
-108,0.457553,0,0.444368,0,0.48567,0.471493
-109,0.47315,0,0.464469,0,0.492807,0.483473
-110,0.486745,0,0.482325,0,0.497189,0.492436
-111,0.496356,0,0.49511,0,0.499394,0.498054
-112,0.5,0,0.5,0,0.5,0.5
-113,0.5,0,0.5,0,0.5,0.5
-114,0.5,0,0.5,0,0.5,0.5
-115,0.5,0,0.5,0,0.5,0.5
-116,0.5,0,0.5,0,0.5,0.5
-117,0.5,0,0.5,0,0.5,0.5
-118,0.5,0,0.5,0,0.5,0.5
-119,0.5,0,0.5,0,0.5,0.5
-120,1,0,1,0,1,1
diff --git a/data-raw/US/create_MortalityTables_US_Annuities.R b/data-raw/US/create_MortalityTables_US_Annuities.R
index a4a91db..c63eff6 100644
--- a/data-raw/US/create_MortalityTables_US_Annuities.R
+++ b/data-raw/US/create_MortalityTables_US_Annuities.R
@@ -150,3 +150,79 @@ save(USA1983GAM, USA1983GAM.male, USA1983GAM.female,
 
 
 
+###############################################################################
+### 1994 GAR/GAM group annuity tables, with improvement factors AA_x
+###############################################################################
+
+USA1994GAM.data = read_excel(
+  US.annuities.file, sheet = "1994 GAR",
+  col_names = c("age","qx", "AAx", "qy", "AAy", "qxBasic", "qyBasic"),
+  skip = 4)
+
+USA1994GAM.male.basic = mortalityTable.period(
+  name = "USA 1994 GAM basic (unloaded), male",
+  ages = USA1994GAM.data$age,
+  deathProbs = USA1994GAM.data$qxBasic,
+  data = list(
+    dim = list(sex = "m", collar = "group Mortality", type = "Rententafel", data = "unloaded", year = "1994 GAM")
+  )
+)
+
+USA1994GAM.female.basic = mortalityTable.period(
+  name = "USA 1994 GAM basic (unloaded), female",
+  ages = USA1994GAM.data$age,
+  deathProbs = USA1994GAM.data$qyBasic,
+  data = list(
+    dim = list(sex = "m", collar = "group Mortality", type = "Rententafel", data = "unloaded", year = "1994 GAM")
+  )
+)
+
+
+USA1994GAR.male = mortalityTable.improvementFactors(
+  name = "USA 1994 GAR, male",
+  ages = USA1994GAM.data$age,
+  deathProbs = USA1994GAM.data$qx,
+  improvement = USA1994GAM.data$AAx,
+  data = list(
+    dim = list(sex = "m", collar = "group Reserving", type = "Rententafel", data = "loaded", year = "1994 GAR")
+  )
+)
+
+USA1994GAR.female = mortalityTable.improvementFactors(
+  name = "USA 1994 GAR, female",
+  ages = USA1994GAM.data$age,
+  deathProbs = USA1994GAM.data$qy,
+  improvement = USA1994GAM.data$AAy,
+  data = list(
+    dim = list(sex = "m", collar = "group Reserving", type = "Rententafel", data = "loaded", year = "1994 GAR")
+  )
+)
+
+
+USA1994GAM = array(
+  data = c(
+    USA1994GAM.male.basic, USA1994GAM.female.basic
+  ),
+  dim = c(2),
+  dimnames = list(Sex = c("m", "f"))
+)
+
+USA1994GAR = array(
+  data = c(
+    USA1994GAR.male, USA1994GAR.female
+  ),
+  dim = c(2),
+  dimnames = list(Sex = c("m", "f"))
+)
+
+save(USA1994GAM, USA1994GAM.male.basic, USA1994GAM.female.basic,
+     file = here::here("data", "USA1994GAM.RData")
+)
+save(USA1994GAR, USA1994GAR.male, USA1994GAR.female,
+     file = here::here("data", "USA1994GAR.RData")
+)
+
+
+
+
+
diff --git a/data/USA1994GAM.RData b/data/USA1994GAM.RData
new file mode 100644
index 0000000000000000000000000000000000000000..29dbbb73edca9b5676203845428ef4fcc1f902bc
GIT binary patch
literal 2733
zcmb2|=3oE==C?C=Ew0>@JD%R7ImKg3^fq6umX5$HNAmZmyw^`QxOZy4<SZ$-v=$F9
zuC0@z_P!H)b?R2`OGkaJl<k+Md2zZfPP$#1?fpG-y~jHqd+p8Z&)0slJAbeC+4&#;
z_CGIvHplXR@xS`}mOkuezO8GWlqPDZxSr&1)v{C#tc~pZyw*d-m9<mRyKAcIsz<Mr
zv$s`jJX#c8F65cmwZx^-Wtq?9l}kNbidu!{s#K~T;J$R>pu$3}g{u~4E!x-ri*Xg-
zD!o@SA+{mRA<iMmA;MQ!SEQ|Yw!$QY`%3GI?f~7uwE^Ob^&FQj)^pwJw$+g}<XFct
zzSX)_(O=45yuEOH@pkw34fh!T8htr^!Fo}-%YK=^MPK}0%)MZ}xZLf1(>>O`+Ew9S
z@?PA%uzPX2`}@Xw%zM?V{J+G#*ekI@`KjsQSG{+ncbXUFKk<P=yA!n&ek*);Jm&;~
z@4D}Z@6<1fe^Phi?}Xoq-<{v}-;v*GUsV4j2MQTaFi&7t;CGPguo0+K_~h`!@kH|k
zc13<CxgHydN`p@cPZUo$PiR-*cih+1a{PumL%VlJSmjZ!LvHU2*gux7jrisgk>yl<
z@3SGR&d!f^#;0Tt-tqdo`^((_ZLCwdySD1|H)QUg?Db%)nt`Tec$w~b!5-V}qb8S_
z4IY(js942+aD|<B`Mz?UGx0$usz0wcGsv}FZL7P5|IEtgCt5EqOSl+4@nc0Oe?;ZR
z+w(;%8EW4<?OL5_c;m<Rd-`Eo3~iwWOMYxrypgr$b@;VNrfuKXKYOdi#kgDFtoGvE
zoqTV;w?<D@?>_MGPs6>Y<!Kw<_`co7|8sld&e^knd_H7QaLaCHVFDlP9R6~>uI&B1
zdb^)Tu9l8vyS_uo{OK;Xg3T|hC%wMJoOk~D3eoGcnySBrwVhiQQc!<n_NxVE?00sr
zf10&!9^bpayN$I=rW{bda?9SZ#=fFxPVH&s$p^C8Zv-Ss9rzvKey5*#*MT=7Nz=B^
zVGZ;66>=>$_e1bmkCoTfGX6iN+wNy1t*~Irk(5{Oq&RlPU0pc!6HnL!3+MA~t~>=>
z6?JrvuY39^@PW$rm!cB7Z5_cAPt~qj-}u~(FSTkW)BU|OQZCzav+VpcMc6SsA-C%A
z;jdhO748^(j`j~zJ`(GFe@DX(CeCUz$<j1dhoG&mo`yd9yrFPo_N0Ty^)_f8FJ0NY
zLgdoAS+<#XzYFH1=AZihbq~YLZ|`gWUH$BHW8x*@c$ZU6GM`)y{^Hb5yngu2zf0@#
z6TjD1+5Z+3w#&=6nP}V0T|GZ6_m7K#B;SiI@w1LSI<)Yv)!OQ&9|aaS9$Y;-_uvHk
z``bTw+^W-9G4D^*`%m@{Gd6u_U$=AShyE!Wwq+a+vG{jjgJ1bPz9}+mKb*d2!|8YB
z%d(62a+-~L89#d$?|%3sdG^7~I1`QiQ7gU&XftSOUl;gzSN_1u^Sgd)t+6clw)#T)
z%AcMEb<bS7Z#|yxT)O$9x%Tci3(Nkti0W>??lAYJnbKk(;eFHil3lf9mhFnMxgpH=
zHMl%2Kf{{se#M3L`Znj<UN4=s?8s(?9hRz@e|{)P?|Az0$?C&re_q(J@!S(92a!AH
z_D4-#rr2Qpe1-3^;G-8xmLGNc$*j8L_oo|yahqAUWZa)Cd8L#;O1bQX?@7gjd)rE)
zZdM3xnVH{z>(`gsjQG=L)ro64`77dYT)k?~zH!>qhmWRcDxA9%y;5ui-=k<-^HtxT
zH2!esZR>c<E_?X<^ZK>=bszb{pQ#<;Xer$2{Qr<U*YSzAHFxivFbcie8e?Z3b41w1
zTIivY<C-V3R}RLndn|KS_#?-@T$S)&0+NsH-$w<zw0zuIc89}J=X&~*m+=|FhO;^*
zZCm`$?`Dpt$;F9j!m^$g7de$%_7uoA@5tEK@Swo*b5H3a4!iJWlB%A%c?{D#j;c$}
z{BqfF*}aNP<vlM>++S&C=yl+s>xV;~5B=AAOq}xd<P}AAkAGjh<MLJYb#r_|Z%x|5
z_)_)2w53lE=c}yb<Ga3)=c>2rNx5XJ1#O=>7Y01LwZeN@h|zqflh1{$+>3W`9+ppd
zWBOD;iOE@3a;ny4K1;7<4!TBR0^*Aw>vf0<*38ogWb;obTI@P0Ua9HybVo1cVC85#
zqt;E2C;9z6v13JqhoDzm6Sqmsf#=?7LhQ<3JrDaNCFZH~oj2o)s|{FI&ebk%+NJKf
z;n3~~?W~jqIxn3}J0BkV#3FO^s*rQbdPSww>&s{Tp7EWlkav&l+4Ee@DwSz&zOrGJ
zX>MyJL8R-y>F3T&x%m6$x<%RdZl2TFpp<5o+A`zsl%C!<Gh_ZTX-UW|o$kM(g6Y@f
z?;e|F&e%@ca-pDG^h}ut@9Ccp)n2gG+^Si*U$A6ykG{6+yoQ^CllMq7dnfp9{`c<>
zqwa;}FF(E9DShC>gTBkM(fSRZ(@q|J+m?MG=hxoSwVspyg%#(!OuNd^TbA~fd-ZV#
z>&uJ(zq0mUz`LT)d&R=0R+Z~4b~y~(i$3+4rs_v5vAnXBBdLC!=x=-OyIWbrF6Fzr
z8FVE~ymxXscZU1HNug3x{!HTy6XjA>z3aOnrSjwA;`Pg!_r|YTqqd~E^lWWytCM=e
z9{=gjQmq&k|2$;3|B%Xt!zuald%smVx-#m_)~P$}mBV!W-qlkxqr?x)+%nV3BH4Jy
zD*h!p-doweXGXnVZFJ(n$J{y3j%hwV&>i@9$CWg(iu2p*B5STQ%rM#bP(AVFgJ#eB
z(;a8N7F^5yuOgP2r#1Vkb)-N^#eL1UYm_%Itq*=tYW+H9v$ud)$*N6W%Ip?f?y7#v
zIw!DkwZ@OMsSFiU*RR``Q^y-vc|Uw_MtA$xDcb!_<|`Wtj_-Z3!;$Hl>*P)RQFmL9
z=-clt&}eSpotU1kzCE_}igwqlWqbV<eGNAoR&G>Xvbme1)usMM*gD<4MiYJ>S+nu+
zMx~tk!>KXeFWTNLeRVJTeYLZnO>4c!Q{@w;=YACIbW?w0c;9XP@)GB`GmE1Ri{;+9
zzTIxFj14!RUDWId-FqDEFI@}&s>x1}pD$N^!GV|SvHZNcWp}xcd;N`-G&Hfu+viuJ
zU86W<#oGrwQ{)<iP43Nn)_67`%Ko(Be+fxHIpsUmr~fSc^Jm|ag1_<)uB_p_f4h?9
zz|7kdFK6FXdLX0rTIy2R>Cad67k(?L>D(zS{AyFv?o-d#w5-2=_G)3Jt+ZbLiZ_eC
zL_5{~`?9yOcJBw@_UTL0^I1$Tgy-(xo-TNYC+kwhu4im_e(CJCxns$h7ks<?SzV>T
zyAOJ<k(WF#TzPTlW#tj+S+&mFWiIbx)LnmSU4fBAf&TU2@@6NoEh}$`30w`ku<wh5
z=+1jAauWM1?nb=j=C4|Fu5*5Pb47PsaB;g_hjp6A{>guLaz+*wm~{13Kd;&I#cq1%
z7O(dOsW+z>elK5?#ZqpvW>t@kOG?o5mUzwi78Yiee67d-J3g4BV!h)2i5XeV^*dvJ
z27mgk?RN6yrrf_-WiqpEHhN#ZcvI$gse6;5YyX*J0s35vH#evY?l)7LKQZ%f<o&3{
z)`!(DUClPyvgYZ-X)65NKNVQ7slQ&JegEIo{Qsx6*Z)7iz5f4}g5Sr#CqBIxCYSp^
zV@vJ(xfPawzIw<RJY^K#Ym>nD`o3wxsdE3PyO)2u8~pR#%9`&Xa)$q}ba+<%mVdsb
Ky6=w`0|NkZ;AbWP

literal 0
HcmV?d00001

diff --git a/data/USA1994GAR.RData b/data/USA1994GAR.RData
new file mode 100644
index 0000000000000000000000000000000000000000..d90ee6dfd05945bbe0eae72ece932b3e5b5685e2
GIT binary patch
literal 3059
zcmb2|=3oE==C?DsBYb{JA73sVyMm+JRha9g-L9bbA-nv<eZNofF^x4bojFq}a;d1c
zr}Uw`6!Gp*>98rc-Bzi(?sN&7cI~f_cF68qGoRV6+WhPG-Sofm^V0c#&inrT@7?<6
zDR!TY(tlkP`75F+(#7K%_QT!Hufluc=Wm;?n=TGWS=!n-WkZDHimrR}uC+aTc3Ec4
zOp&fFvo^TQQc^YU3h}mUIhy`^bEM6~y_;6dnJ4Zp7+?EJ)Ai_)k2^TuNR+K%SgqGm
z+Fe)pBQfGq_4Tbyi$71>694=7EZN<OTdO~p{4kodeR`FaxKX#Q+>&1Hwl8}gbK1)r
zw5Q(qf9M^XtipjG`@<IfKl<}iThw<UMb#r)-@SZyUs3eA`<djQT`Mm&&*3{P9x;DL
zI@>noXY*U%c9sghR=M_fU9HM{?u+t~ogX<=pE-sl&3p25((}K1Z`a2s*?0U^{MLQL
zKFOX__1nxff9Dp<Y*RK7dOfjJ?^nv)uVqVXK7`K?kC`TV@L6z2`R>q|=04MxpH{3|
zaX7$zQJ>?|MSbpbJ65$?Nrp;?vWN0tWw<JEm1Am~74KDvt1MSVuJT-!xyn>(^@=CN
zJ48FgVX5#erLgRtwT9ac-e60P?#VTr*ZI+>=dHoE!#CuTt9xV($`0HROm^>?YfyIR
zhGcSjkEDdLLD|6@tjXa$v4&-bZ}29U_uMs{*ZNV%N-Ho<D~$h(?S=0S_nB+ee))Wv
z`NH_(cbEIZwf<jhFMfBr&t0qkW$%UWZug~Y^S|u9_}%e7d#(JJ+6&)Z?~B*Qf2qCr
z-T6L$t^Jq37rwjSm#?k=^7kUYgFVAvhA+$)=J)=~`I5faUiz=Y7v_unF7^U{6}~iI
z;CHg;_-iS$=U1HsNAhNtiJL3#OqgeW$1LXNV!Z_Zh#7Zk4333~_xNl1^yGyyH<wKp
zl@6_Oh&Q#de4S<cVp`m6mAz*LOP=k%&S9at;H~&|fe;P;;G}Ci*4wV)2!3_u_^+~B
z-r(GsZ#?gxYRsFr`X%RO2?L9a4b?kdhAmKyDSEUzk7@E9&h~4s{TM#)5u5zuwRl3C
z+}o3vzZxYh-T9{S=v%giq`51<=Du>cZ*!nx_cg%{d{qx+PwZhfOIuZb`Pv?q=;_hE
zI?KIgocz7!Rn<wO4Y#ju*}gbZ{!O-6UgU#s0y3|rJ$;vcts!~hmwC3|^&2mLnq+iy
zbN@lfH;3i(W-^>BXIpvEYf<AymyMeeV=N8)miJ%(`haoqw~Np0KP_%3{Cwb-uA<$8
zq^pJMw|{ExkUE$=w|9Tj?UmNErde-Iyf(-1u&9_mYroDh9_@p13ET3{btJxQxc=Pu
z^t6J?2aC7<pUIeO{o`KJn$J`0x&HC)oUO*5;=JHw+Y4cTD}@T5^v>W8;|*tg$|R@O
za4u}V-5>fSm0|WQ%b2VwGA?T-Z01~+b--rXw5OLUpDU<cxNu+Vhrp4Y7GgGbM!XfT
z&PN$tu2k62Eg@B6$1Wu0ZvAFoC67=x8(Z<~6(43MxpaKG8|Cp)?_SEfiF_*szfAdN
zdOIQVX}(`to!W+^lgnprwO4hl>I$w1`XCWen^oDEch%s;tzUcgEijYVB(AeGzJ2zF
zbHD#3ukx{ba$%+SWbNmN+@>?#Uwo4L;p#WHrP_ASJCK(9)u;Os*NNY`m%dLt&wt1~
zrdK?DmHVM>2j$kg{XX0NmcPd8`j34rbFa)*@Q`QR-M!}ed|&UjKfcdbhAsTQ;DgH2
zEuniQzNJnTI}vl2*Z$WEv5?od94G#|<8L-+YeVx|^D@Jzl!ZH%TCCi7m7$O4`uv5_
z0X)Zx%-MTlp71RF694SE>iRa}@XgEqJn(DYv^DJ5Q}%BTl_{>@x_1XJxU_on!h6Z0
zbNF^wuUMMN7V=a}KL4Iy`>WemgRdN$+wfVrztuk~hkxtk9DTR!sZMiC-J6!LWze(G
z>yQyF<<r}hcUD(AlyQ2ErRhPY$Q?4VlihBc3%n}k`5HU#eS@~R-CwIcbJ=*ca=%P<
z^eQ+Tq?TS<BWGgu&h8icW6hZIGc%&U=GdE9_GQR!*u^qytKPB|$5%Uj);#uk)kNus
zZ9o40xqj={HS?Ipy#{@!`=6Ko7xxq4n^5$*tZ<@IXz0_;mTv+aMa-|w6V-a~sYREq
zZ13K_;G8K6T}L<eC3M$ymi(@lt(wswdMNSshQm8;wy9tI8vScZZ|c#GUx9B!%|bjh
zLj;!vw9Z-N+?!Z<V@VjZWlUH^NDj-tdy^LJ+4(C(TI4WqH&;JbYA5^jse;o=C%p2|
z75R5r^G~JjAIZMCykcKfUO8+LdJ`F<m*)}qPU!l?m+PiGDaO~bP2hgI;YZCSq0m-!
zWgD@51#vaYZy0=Y`pEJt&sW0xfRoQ;r>V+cl|`Mk0y-ZQK7PEy$87#tZ-Ms#`O%%4
zDW(tO!z^YyZ({5CHbvG=YpVLW#ywgWcM7habl6F9UxDI{XpQNI1%*`fTqj5zf45m@
zkH1K;nqIAkcxK7W4J&5IMZ0Erf8S=Bw8Mt6TUgcOi;!(@)!vFPXM^v2-`Sbm>b6<w
z1pl0_?D=&+ewgfemzexg?^JYTM4QI{|Bd>LSt$yy`8S+M@L;nks%BdAE0}4^iH5KF
z2fx10s&DWP`}FI?Yo+u5Yn%V=-p8O&FX-^k=|L>R@A9d8&h6*k@t)~!Im7$otot_E
zzAs+AV@lsTmNLU>2lTXl3EgmY{omXayK~<6n^t%Bhb>>g>GSEmxYT_qCN`ypFa0&m
zYoZSEKRep4?P#oVBR#!N%SWc>;nx*E?mIlLoV}=IqiKf5D%IIjZ#NyglAk0lx$5u_
z!If`~)%zMbJI>7Vmb-aZT9Re*Y5)A<n^A05@6=~yPEoty-eP*BtWCGd*?Gy{11IXe
zYhBZ2C7#OJy(mdkm5?pJvyto20S>`QX8$Db+ivB}a9*zVH$eCPV&h<^O>*<uZ_P~l
zI!AY!LGIj3f~98@4DO{b4o+Rjx@^iTlm9!E3})*mNtbSzaCXYp3mnZat?SLJKD=hz
z@MO=dt2Q<44e76Bw${#cSYNSo>x+w78v?`r_tZxT=a{$6l8w+|u)f3I{xY+CsvPG5
z&)N3zTaR)}e2Vuxo_(wN^u^-I+Mf>{JiE%^(wC=&1^X?{p5L3e+~D;?#p}<W>K<@$
z`EK)a>)J!t<n{=M%@Yg~`>{#9qR%#g|6ld@7rhbDT#Xa?XWnj~;<(}Jso0mMvkpkj
zy0Fq=Yr_QN>$7vU^_?TcT|X~p%I(;pp0w7iBkq}OnZXm&g7O^|xBbKz{PSPSUtV6r
z-Q-rbbWNeQj>8qdPv3p7pIg<I61@0MrgzxkONWA{hb~|J@H5-JbIE<78)7<Vhnspx
z?fcPUHqY2r`V#xAzPf}+Zpp_&%IEHgJEX-qCa%5pd_w8(n|FhbGo1LhbIzlE`6_F&
ztqLD~+G%nl;o{coZ~ZKDs$#BM?(txlJ=IF&&8^vqt3#s>-HeYvw0&`WT)>_at#5<*
z-CoE~YkMD*w}JP`v_pAsOAm=ob^6C~`AGHi_s$D%Zd-IT@7IR&*BvYO7n}Tfapl&f
z1qEVVLHj32Pd?tBbaPI}qUH7nA9mIB-i$xnmYlje<@=f5wo9(xzxibD<WYC+zrMpy
zBm49Odz0tZ8vBdpO#gg{SxY}QBW7!j{G7UovTTup2XAi7t*-xLB%Rg5e|^tn<p(uq
zR<C?D(KBH8_QjF!Wkg?D#a|4HmMEC$*l&DGgL(SC@V$n0H&j;0*KR&Idm>}FkyO&-
z-4O*}J7+zv%TkIl*4U6A9mwBXcv?I*)3E8aYVWTv34A7|uA(=l%x!#jC|+vPlr#UM
zx5!?&dU(R8z3)U_-dlFMidO&hO4%fFZMxy%o3>LdUz|;^<y!jvM)U4pN-5^7``C*@
z6wiOHKd>n4Iy+lOzQ8-@^g|1u)ctNfu#cnjF#kKt>-%;FFYL|mndALk<C5+3M)QY%
z*Pfqj@nC<?+Dnu7U42~_XYXP+RV7WeGtx8Zn{AAc-P<s$g_=oIKWJvHj8LuV>iT#5
z(DF%FR;zr~^4weZWaah$J^sJXHhR4Xaa#WO^2-TvEgknWYdU%~P9}O4&6p(fta^t^
z&yBSuPQgdC5?!tYZS;~56#d)ZDyVAd<H;91>F|^%2C9}moT~qfy_W3xC(Oh6s8QXD
GfdK&at^OJS

literal 0
HcmV?d00001

-- 
GitLab