From 8649822f0e81cccb46377b8ec371610c9c717f57 Mon Sep 17 00:00:00 2001
From: clblatte <clemens.c.blattert@jyu.fi>
Date: Wed, 11 Dec 2019 17:08:14 +0200
Subject: [PATCH] add xlsx output overview table and modification s on Main
 skript

---
 Main.R                                        |  40 +++++++++++++++---
 loadDB.R                                      |  24 +----------
 params/Overview_outcomes_SIMO_2019-10-28.xlsx | Bin 0 -> 24757 bytes
 quick_summary.Rmd                             |   2 +
 quick_summary.html                            |   2 +-
 5 files changed, 39 insertions(+), 29 deletions(-)
 create mode 100644 params/Overview_outcomes_SIMO_2019-10-28.xlsx

diff --git a/Main.R b/Main.R
index dc1831a..8d50497 100644
--- a/Main.R
+++ b/Main.R
@@ -31,10 +31,40 @@ db_names <- c("MV_Hartola",
 
 ### Restructure the SQL database: querry that creates the final UNIT table, which contains indicators over time under regimes
 #
-# !!! Only needed if the DB is read/loaded for the first time !!!
-#
-# source(paste0(path, "structure_SIMO_rslDB_FBE.R"))
+# !!! Only needed if the DB is read/loaded for the first time.  !!!
+
+# source(paste0(path, "structure_SIMO_rslDB_FBE.R")) Afterwards, you can exclude this line with: #
+
 
 ### Import the restructured SIMO data (from UNIT table) in the R-environment
-# It gives a single dataframe for each database named by "db_names" AND an overall dataframe calles rslt
-source(paste0(path, "loadDB.R"))
+# It gives a single dataframe for each database named by "rsl_db_names.csv" AND an overall dataframe called "rslt_all.csv"
+# Select columns that should be importat from the SQL table UNIT (created before by script structure_SIMO_rslDB)
+# An overview on all available SIMO outcomes can be found under: params/Overview_outcomes_SIMO.xlsx
+
+columns <-  paste0("id,
+                   year,
+                   branch,
+                   branch_desc,
+                   branching_group,
+                   Age,
+                   area,
+                   cash_flow,
+                   V_total_deadwood,
+                   BA,
+                   V,
+                   N,
+                   H_dom,
+                   D_gm,
+                   Harvested_V,
+                   Biomass,
+                   income_biomass,
+                   CARBON_STORAGE")
+
+### !!! Needed if the DB is read/loaded for the first time, since dataframes are also stored as csv-file under "output".  !!!
+
+# source(paste0(path, "loadDB.R")) # Afterwards, you can also exclude this line with: #
+
+
+### The single csv can later also be loaded into the R-environment with:
+
+rslt <- read.csv(paste0(path, "output/rslt_all.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
diff --git a/loadDB.R b/loadDB.R
index f4196b4..e899767 100644
--- a/loadDB.R
+++ b/loadDB.R
@@ -15,28 +15,6 @@ library(dplyr)
 
 
 
-### Select the columns that should be importat from the table UNIT (which was created before via a querry - script structure_SIMO_rslDB)
-# An overview on all available outcomes can be found under: params/Overview_outcomes_SIMO
-columns <-  paste0("id, 
-                   year, 
-                   branch, 
-                   branch_desc, 
-                   branching_group,
-                   Age,
-                   area,
-                   cash_flow, 
-                   V_total_deadwood, 
-                   BA, 
-                   V, 
-                   N, 
-                   H_dom, 
-                   D_gm, 
-                   Harvested_V,
-                   Biomass,
-                   income_biomass,
-                   CARBON_STORAGE")
-
-
 # Variant 1: create a single dataframe for each SQL Database (selected columns of table UNIT)
 # Dataframes are given the names of the DB and column "test" is added with the name of the DB
 
@@ -52,7 +30,7 @@ for (name in db_names){
   
   write.table(rsl, paste0(path, "output/rsl_",name,".csv" ), sep = ";", row.names = F, col.names = TRUE)
   
-  # assign( paste("rsl",name, sep="_"), rsl)
+  assign( paste("rsl",name, sep="_"), rsl)
   rm(db, rsl)
   
 }
diff --git a/params/Overview_outcomes_SIMO_2019-10-28.xlsx b/params/Overview_outcomes_SIMO_2019-10-28.xlsx
new file mode 100644
index 0000000000000000000000000000000000000000..cf37eb1226edda814089f6b4a4399148ab987e57
GIT binary patch
literal 24757
zcmWIWW@Zs#U}NB5U|>*WxXsaIIh&DzftQ_uL70JqDcU(duOv0EBtE3FAhkGFuOc@`
zVG#?H#sUT=hEXu0Ltsm&xA$QOfqmTn+Bbw}Zqw@U2<Erz(^~PM^sL|PrN5KXSH0a|
zzG<b>rS`2~CrvnB)BiJV(!r-^Za(4_n)N`tFy@;p<37GyKMGan&CjcTzQ`dcPFvb0
zB#BQ$uf3hUd3hy^&W>Bjs#1$86jiVBTC}EiJotH1LuGk~{^ba_Q2#3uPdxY6tG#}}
zY4df(QGvq;O?$3f2`V&xe0%Ps@W-x&SNC2I3UxY~;k&k@^_2Pg`UzhSWu#<XSUAtB
zP0vs*KuVgSXj{)ogXAg}C9^gsy{A42h1Qo^S5MO3Y>?w~M0wl3riHbeZe3mUZr)=L
zH|LKFgzjn{zxVXasqYii)+_p)`l@hyo7(mjC8vJ)-#z_f=X!okmnScTa@IfC?4Mhz
zZzP~%^x^8ttvz#opPRN&R)@DL`dw9*$enXvWY6fI6c2nJ6Y%Wv$9cy2?wZW&yXTY&
ztT!$>epfYR18<8*G53!5oaYzJo*_Ne!Ct6d<5;!$2ifYmeR+OybEdW(c-ne<@B8=#
zy7!*#jb8LH#h?4Yv43kGez7+_@yA%Yj)8&U|9@r%)I2R5w=b26k%8ep69WS`0|!%l
zQEE=Hz8;7i&C4S?3&TCXW|FTrv!OuS`=<YjXZBj()>68Zqwo5{H{~y*;mp4Y<})WO
z+4c9kn`6nsZE{s})|bz#ocT3v)8dN70@Y<g(_A^5=A>QT9d^k3?~Wbpf|HBfdf#{)
zYE;Tu^?LPJ<6f<&3vYb2>ss?@QH!4H{|1Rt_pqZ)wmV}oEh45CNko0|;$Nj2Gxe!#
z+kL5kvnEH+d2m=4GcU6HvFGO=GjrednWwArPcBu~p0Ld2+=R#<6K*=a`fR!=Y~fOs
zgB$8jz3e#kTPc}YPCe@H{1@-DF9rp<skY8yp8M?Q8w=JgdzQEs3)kP7=~naHV2=_1
z@`ybBlZ(C^ZnRB47k=Jc=-j<czp!SR!rq<TADI{!4zV#X@G~$lROINF=NDxs<>zOE
zs=zfN(b=~RMC#<<*h}o)Jx%6*(yS>zjxj}R#$LDMyBa1Wc-v@G(c&-r-+S!bx~1om
zugav-zArDWub*C>8yuajR;_nSP)S@Oa*eL*(c3>c{_fI#!x6c6QLbzD#Z7{@Hp<2J
z?yuW<hht}rNKz)F!|nbHPAm1w=Io2v(DqGyQL=7Env+?`!LW04_f75z``cuA^-5oN
z;RBN?+cWwFPWLBt-0YdIBm8t#;@Vry9$~?(Uw*T`|D_xiT*tC*VbRW4Z-Pu@9c(_$
zT0NK1H7xkp<w|Xp+BFB?=zd(gwk4jSzA3S&U%&Cy*<fD%r7w>9#4>Eztp4{{U6ZU+
z+|GpG*VyHqKgpcD7IXZ=5ihp+DWY>WG#N`A{NLlIu+?YXkq3MH7yVSv)02unS8?{;
z_m$Zf*X>@`+ZWQiX4W^Y^*fTTa<py`l`a!&2ras*svXySt<})s+>a^e%r7xT9DUOJ
zBs=Qe0<|qw+K(5v=TF)4pn6tqdwAU?(>*VRQbOD3%wf_}+<#c%laqNtb5g_);Yq(6
z?RhRGO%lm7@GP)6c~?o7X@%a2l{-XAmWr`Qo4n(w{AayNzNGq!fXw=e(q|t`J}fGk
z<kvAbEiSJ*aK;VKfR=63f4fdmuj^7jxOtmv_*>ib+5d93-ix};elIM^*hyf;-R}3j
zg7zQ0bPsG%=xO*8_pSHiPo^Dh&gU=Z#QJ8o%=?gdV?$cip)$t9KJIfDKHBpjNs76^
zX!X@u31XMCGwgRIOKf*~lJvJ>PMDh3xsQIHuZ5JBix(GoY~a@UYS(&e-Aqr`HEYwO
zzwQlB+Zr~BIr`B<&yMzay`LNT`=oc3J>iY9{%U0Y{jd%5bKC#tv)W6P>S}Tp%=+Hs
zRjQu3vt?WL=`_x_imSh=G2Ne@B5d#VfkC5+RWqlbXX9Lz#`NEjh98_KEfsic9Qn54
zrh2w+(V1Mgt6ol**)OHZzusx#;eY66Tmi@J)yqWhr!8yxw*7zEh8@@0=a%0u?2-SH
zv*@|ZUh{v<*eaC!OVTx@7#SE+SQr?j7&sU!a`Yi}3UVHS)h~=13m6#T1<eR$?3v=5
zeb|7f?S1oK_Q$^$sodRV^y-Td+ul6MH(4@In~#-kxSIR(yRoEqu*OZ*FG&+x?tkCd
zvgdB?>uI%%KKGh7sRy2sluGSB7Sgo;xqi9yIUCR56W63f_PC!fzxnp;H_>R%hm+RN
zKDT*6(w+yVf!|pkg!V@LI#HJ4v~-J~H`l_|3nsr(&~7p5yZ&IAl#TV}l;2JT+x5fd
zJgpU<^Z)v`)SDVL+g`1>c!8Ds!<Fq18Fi+wa}xSv-p8cee)?O#_pQ~7)mPg$-S!Z0
zHAq{icxAQatP)QtF6TqSCj|}_s1@^^oWv7<O}3`jL%q-NY2u1Aht{mOomaBXl>hp&
znQzbjT{Jt`yZxy2&JTt@;u*94W%JJY%X<@4s$w=(3d;3Gvqc#g9#k?gNFd7H;*8YP
zl45-@We6&O=fp<$M_cyQ#V7oCKKOp_roQ)%TFvW<eqOr4$Yxz#S*GO|^kP-f+S8uj
z6E^;@)s&U*eDYm$+wzN^=N>N9=o0db)%bNd{$EYqu3sDP?zepxmsj;Lc7OTS+R7WX
zbu~Y$F8@CNKjM$k+t|IbX8QBZ{dZ4)>-%2*_V#=6cK=@QulbR@{r}&$x67(`2yyA>
zoA1qUKm7Xl>(|@c^KU;tdwciY*WX<C)!ckL-<-YjZvEc6>bv>o^Me;YU;cX8effC%
ze7X24-PZ5Fi!<KrtEsHHaju;I{`~Fw?bA-w{oZ29@jw5*T>O2jtn>WOCsf~iQ~&?p
zuhWO${S5V}sj8{6`J*iOHrCdDPkept&9~<B=PkKg-}dm?!&m3`{<(i*zGJv=le6#n
ze*WKgZ@)YJaBjf9yvk$G&aX?qKD*tA^LD-1OoKi5^6PDm@#!p?Ro}PN;MA5`PB|a^
zg{0Gy0%LxEzb|ip|I5Ga^?%-r|1>|c|JT#??{5G9YFGE;kA|&Gbm`Xnbvs}FmAdy|
zY}w>{Uyt)=vi!aN_*=jE*P7BBHNSp;%lut%!%KHzdFhS4|89Rf{V?zM7nP;Qjrtx>
z_`BuyU%BpetSj}m?pdm|aoX9#4?iB?=K4QZ?QH$;`5}ch)4bLful=?2al*v5pZC9o
z&)*(zc7HzK|9!ucx7Pmo^;TgG&wXC^iH8oq^Sto%{p#n(S7aL<F{@uz|G(l-&5yFn
z*=g^WExuS_Gq<NIod5oPzTMCN=k0Gv{<&?*hiP72A^ob$pNq=Kai7|MylHcJ;thH8
zE&E<idwAvD&u1Hc{z%;a`>~hX%7TC2c6D0G-{W5SBma25b~{T)Zhlk855;LeO~tZ)
zKV9}+?Sj1gZ~1lg9Xj$)`UTEiSsL0^6>j6yI<F)<+~&!nl}*t{GIWl~Ch_PTlQC?Y
zbMnR;rKUd5#<Xv4eOvySJ#$$q&Joy{wn+GT+TYw)T37C;?Kxz)HeiA4iKUaHG-vm(
z<9_xmdE&m8_c+%EFIeh#U3>e+6TbexmMxf)=b3Z2R3q^0qk`GeO-uE8v@dQ=ROi}q
zCAex;?!KQds*>~UroA&LJ(n?;Nw;*`xvbY$<fX667k}RPrG2GwN3c3??}M1^$5WnO
ztNYrOv_EW$N`&)7T}P91bH&X4{PU6*e!6<HK5zdtz721Ck4ambxqR^T*G_so<=M1{
z7IN>5{N!>3gDTc(L_cWLZPX2CUSGnpoN-;$ZBw6lN-G_7rtXxwa*S<zilcWXL-OGh
z0>|1kqxapidvr`MW3_2PR9*J|D(>w4VL3Su%=DYO*X)_OHGc2Ht?}7=&%cekP?p8<
zH*a56cJ}G8Gr1ew-mPGG`(%;D>ov0kZ}S~nxcyk?BIykmL|gpMdK7d^Y)O=En%L1g
zN6Xsoy%Ouei_BRZ&yHLPGF0z1ov-dPW5ZmIy;Bt}nDw39mK$wLnf6pYWl_^sJHH##
z4o;M>;5p83-xe&{(rbERxk^V_)f4u_MHNW`LURIB1<dZ4sTasN&0V^m^~BsZ3GV9!
zr)M4T`%_n8ER+A<*zcQ@V9>JpTMwRjn-{}q(JFOtg{sAY(DvXRm$)oe^(N#^oz1ar
zlhok=*Exa(S09S3+p<?aye2t0_uI<E+vgn5+T6E3=W@~X#BTq#VCik*lHb{SKS{XB
z9b{zV^gDXn?Ang`lN6%9ZY)3A`ACyf$ZG1#9}kpn&)X3Dv~JtEf`plp*;`B0wqE_<
zx~=xYx=mjt<{Yq+EIf8VVPWvbx6R!-)rZP-zi*rSK%x1CvYM!>WXA1z?rvwT)pCw6
z2p4L4?0as`#g!T-9jv=9^K7%$Th_gIl1kUP<3CpNADnq*ip%ew%RGx_&uEx5r%d%T
z<0GGJuSIu{XR=*P_Hp#|tDgF^;lfJI&TlHqC0+cDSXCWmWSU$2d_A|QTb@;@YS7f+
ze5BCa;@9i=)-$a#ys%+e-K+089DHX4jQr*)8l3Vw|8C(GCI#<?H4K}s2^hD_WPJ~b
zntLE@|MI$RZUuWp_nx-cd9`O#@V#n{tqdoG6to-Gm>$^{_1orZR*9Ffy|j#w%<|ji
zlFPTuno<<J9YxWtUn}=4zw`F>l%0RiNH;IM_mVL_bKkiv>08%0pWfO!e`cku{H&Af
z=DW;Tf3<3D{>iVi?q;vQSTrwaN#k47)a<ug^^~u#zRqg=KD<5cDg#L9Yn1<l*H^7=
zjW4AB6}_E(>e{QVY<AaIe_xY+&zJve_?qie)?PjR_FF*yo))_m8~^o)zFxnz(e%+W
z<EE=G+rIyPW*m~UFHHN{wd^I|w>sZn&vQ0?O@#Qmur=wAj^^|gtdjY4uSr1VgP8lA
zP0~jF0;*?D2}l$ko_63`XY0Q01}vOEvvij{J3UAAYxeP6@w2<#-|mTi%gXqKPocQ6
zB<@+9t^SqILAL9!eBM=SW|vfQ(Z?of^}EKyhv!FLwAP&Scd_!5Gne-|#jbe1v?^}p
z^U~ir(q$KU>Xzu{zNp?TYQOtwDaXGDXD)G^6S#Bw`MRInU-yKZ*Z#UEZ+`A<y;)!S
z3pU9yMR0C9{^YKJYKUd4WUR?s{~LbE&Gw5GJu=P}+>>{DI-%@{tOrwd+#j1pg@s>V
zgbAc>c|1#r#s2I0C8sq$woGf-%yc5k`OFK+=>MfjuXPfu3x8y?I5w?Y*pl^-`P${V
zNrLJ(56Y+C?>m2Ej*`W~l{pXl*81BVO?s_0wb%dYSIz^n4$qiQn#Cb`benAb#k#8!
z)7cr%ec}`<vr!b$N=vSp_}u44UXEEElhcF=u6KT|n$>F<AFi_D`O(=1O)P23k5e`u
z^k}}>ki^}f!?+3LrZ+Eko|xFzxxA%=r>Q0Oa7^CWoR}x#$B*d$de9!-)S_eLWV*R<
zVtmiW)|X2x<%OBRVJ@;caGp5VK3FInkG8KobxBbE`VY=W5)%Vf^#r8KR~Ak5_P!Ap
z`ZuIsfw8-10@p^vpgA9oMl)H*g}C@iFL#}AASY!`j+sX0{+}#I#TX))QdlQ!ZM#tx
z|3qCq#nbYG*j$Cg>jq9!Q=dFh+CO7MZU4`!s*X-#o$)duVN)h||4BMzc9E6o<*}o4
ztLCkKUj21Xz<KYldxFlt{&utMo9pTm5<PFigPHTSIE1_=2?Q}#Ogq|r=SEA#{Kolz
zg;Y`|3a*l|^8Mi#$Mn0AJ?-@ut9N?BTFaU$4CH5*^CLND-Odx{ekX$XAKl>)D&tjL
z!o;H6E&lTmXJMlAenXE56I`8Eq-iY^|LJt_Tyb~EB^^%&qwpW!UsW!gH~m%R(s|e4
zZl3r&r243c?;N&@w>-8Btwkk%t>tNAVOrXvl|5<Yw}kW$9yaGM=(7m<-WGT=QE(OO
zoNRG5o3h;>f*&Q<ta+sC*!0r%*h#bIN#ABm2IMeKP<04nh|KPgb+~0X@1kYOLh~=*
zYPvjbSoQ>@MjBm>`^Y2LY2Ivc%H{JFhNKp)*OOM3b*F!bFf+^#4$L--N?-o1l2u*b
z^U7x@Tm7q_y=>RtPQU&7mh-%Wx`p}52GZ$o8cyjN&e*uZAno+xZ)xctJoX5%EcX-S
zbXWV_ZMu2k-J6m&Lgq`D;`_xT&qdt!oteDm$kcON?1Pu-9-ON$@w3vxk?HcZ!cCdU
zOWN%s7`t>CHnE)GQpj%jQ`Vv~)zCTX*`Cz&r^0u-59LP+tBAO6_l$CJ{1)c^&qI0g
zx#(+<RAA~*XYJIXF(>%kq8D@0=3QOFJw4E(l96SCrNcYMP0*5%@lqVa-j@^ALTi5X
z{P-|2UM*1N`Hgpd9;fu&JuFYGW(_<)XY$!YYJLKTRs<*7TQnNivsv{Yy^_1L>-H?m
z>$5DC*X6&_t<$QRw0YX3X)`VVgjl7NUbVcUoW9-oz4UvNBa<$~<;{3rV*IG(SY6Y}
z2{m{A#=2>}TBmH+#Pw)fPv-Sy63X+^dD1$5T1xQA1||Bmx3Gxr$(6kP&&Sa1&e<F1
zf2Ez-$(a04BTue-Wnt9rC0Ezqo$`ih+vdHKW{7ToS)J~%?BK%nu}+zV+~vv9u1(w5
z%}Q+ae=%VSn@T|D7A}_ZGdDTj)rhb04>J?qqt+GCVp+ezY&wU;UFFHAR{1G~YL$m~
z?hF-tF!zZ-w*mV{<&|NoJ5TNLITWgOM)J1t$5U6QeZKpyrj=*;Ijgj<UJ6IUQ|H@$
z*&2HJaBoh&V@_H0FOzav-J^3kPp@FjZnra^p?l{lcjHZw?50hx*Q|eZPU39l+kGkX
zqt8xwFwc7bbKi{|<p~zS2QO{Dd`fP|kzFBfGKSp(lkQIYmiuA)Qm$!}<{Hkwd`ilp
zDaWEb_3WZqFLbM;f2owu410Jk#*JS~=+>>BuS!ddD(C%@y*<He{xiO1&sHX6ZtIFD
zyA!CbbN%7omwSvA+>gwgI7RyVk*ITj#L^Y6pJC>^r#!99<k_;y8COHJ<Fw9y|K#0}
z<WZm=)ls-CqQ!c~qGe61n6p-iuJ6}=*H_-Y+I*?-u5_t)n#MZkx2TvJ{7ZUf^2x65
ze74!&O;_}_|8BC|x8q^sa;>1^y`5RVHBRV$-(zxF|IprBrkC`&|0=B!GfCXKhhz1Q
zQ<EB2pRPJAQO<G4@5vR355mWO#NHJ4K9jGRk&>CLshO*zf9cpKrCcBVb=Ge@y)K5$
zsqx;VGyRvU%=?Ga7b><b)Lr2&99OCS{rl{F_q82fUAoj)rLLK|p51e`t;el(4t`AE
z7A*5%xEsl=cv<oOhj8Z2x)n}LJ-dpZE-I1LTv{ckBJ@FFVgq|8i^vJD11jPtYC=ON
zrLboFsteUjEZtysa=|8^Lkk^OTF)-)zaq{NrIgXSU@N2L!5EPk4aNp5)^4t!cI#OB
zZgOqemU6-0HR<KA$&(}UCFNhavIOk&xO}?X#A@%3XK^A8ii&j`j^6!}=vI2fsQ2A6
zVf)-IakEYB_C>ilbXH8+r>}qM#V4g+J^f2B%(`ct`y%Rdi`7)o;#tUrYVThD2h&1q
zbAp(s*(Qtr_B+|k+R}9C%DQ7ojEpL0zeQZCU@iI^H&b|`hyZU(GkcThpSOx9!>{Br
zXbC=;m#}z|?c_@y_P?8ynD{5zeOC}{Wbfn<InjB*Z_g&-6PB8%mT@;MDcs@3Het(=
zWKVA!kElnEmxTByP1Bf=eKPg4_N04@8Y_M|Y>1G5Tx275rg)y)(d*@VSFf`;5_M(4
zwP4e)uDAS_WvThLT}!&dn9KC0bwY{wn#vvX`dLFiXJucS@^z6n+p?XPqo?{h-1byU
zJm|}s`em;1V<nH?uA1`X6-KH`o%&LuU+?X-wEr$3c7CRi;F4qQ6T8J8PT_D^<Z7YH
zG@)n7t0^a*JPELMGs-p*aF>%1`*vIK_4CgzO;=2|K27%8>v&Ak^W?>ZOS_e)eY+z0
z#n|6>f|RIiZu>OPs|#7Z)r19jWIUq!+Bm+=`FV8WlxcOr-7GwM_3mZ-Du>-JAJ}$Q
z@HEG^+wcFbOb}eTl9iWz`oYB;Jof0c@NQrF){jXzXZyyF!RrK91`0OiS%m!2^W1aM
zp(H;cBE@3MzE3Ul6oZ1f-*PkWHVjDd-P7y8s_gbOvp@rmS?U3Wx?2VI8LU`yU-(X`
zhRd04OJ%j|-s-H8Fkjjf6&oM6hwWt2MV<i93)wQ<uVVLcXaxm2&Y3oK%7X?8r%AIG
zHX59E*ifT?PjrdVLfdJ^#Z02bFCw;dueE7wJ0g(u&$-Dp;)@{PXMxk2Kd;Jo`dmLy
z<|%xqSIeizV*^KVU8~ZWmjcahxx5x9L#JsQUT^k4Ve#c6e`A8+vx2QjMxG{GLKh5|
z^UkvR*7Z`U!7_lU$#urzt8Jcc7gq>49FKYGqxWpehciDf%)RZffunc_n}MXZl&J6X
z8LMSZ@^%Kc7HkIzTQv7S?Q*(OWZX8RwQ1>R--6`=0#VndZHs8Jw3%SMQ(tPT!3j@4
z={VbGpW^0aDQdAK2|lw@NqyU+878Pa&2locv*GjyaYYLRSh#1N_Gq=7thzXGqFdii
zFCUjPEe}PVQZsDyN;t$4`TQUI%<FYym9~8LBJ<Ar#wRB;YFrwZXu8=u&2Z*V{FG^;
z6?MWeyt7l`q><E#4-z0Pwk|s1$2N6I6rIWBUwB;Rqril_hBC$(9C_=y4BF+CrUjq5
z%B1Y`kLRq+M}g!Jjw4EE7PcF_G@cc;?D&LBVH37IRhjF`94UMI((*SkYD>3VeIvB;
zhurnGqWq4Vd`+}wEWaG~McJk3g`(2eOKZ3C^dDTZ@xooNZv|=vZr(;=0?#fuA1UhE
zd1~%~1fhtGXC3Vy+z%&to%<*-AuMU$Da~hDT#1`}G@f}#b+9S7a)^29DNi*NcWm96
zwIIOji$j=@3ZqMNqZ3Ov%&AS1F=?H9+AN;68+_x~WX{&)I>VV?@U-T$tn7{{9~wRR
z`hBE$gq2<@aTLE$OtL&0<Et3+twYK3kz&Y-ix;O%a1m^={35V%PSUzn0#khjC*56m
zP(cxDi0HFbhEAHo%&)^XwH^2%H|cnzL3?&bgr8w}Xve8zJSWdyNGRby2nwH$ry|9n
ztt;O29}!_+W$}Y;5>w+GzZ9L+M>-E@>Ft=p5h<<|^O2RCecELP{(nYc4^IoOTiJNJ
zxj<>QOUH|@*cp}^RoLg<-MQ50gP{KRkJXmh;riZxQ<uhHJ1%ZL-E$^qhm>o@%(7oI
zX52n-Nlswl;;UsX*-}r8WTqTk@b(o$(RbEquNywJvRR9T@UGMJ-*xBwlcj5~fy5_?
zFqB#;uQbeh8CH6H-Zrtq%xcpYi{CBSs*<x%wuh5t%9d5Df;TpQblR)cXH|4gMshNf
zN6WHhs}HH%S+G7Q?D+2Qt-ABs{n>2KE$rM<Q1<3U*zua92FHJY>(b@eqW1kGXY`KZ
zM=!3m{Ykz4+1VxhPLs6dy$k<uZYq1%=^(OV-wFANrz`jFJGFlQntuH<l^M7HR8|K5
z_f*Ke&2G2&YOcccFSY(2w{--|Z`H}~7r(qNZrZQYUi&?6dq_+@$yDHW+pA@g$&x*H
zUWz}AORcyJt{z`KQ2(m5X3~S@C6W6gt_JH(i+Ljd<Fi>N^UMdk8&%Fas`#Gjop!3v
zerKJ1LrGm@!bAP(sY@O)om;54_?Z9arOq}UyN))e$W%@!naE_x88VIKNzs1>Wg`x;
zi41z28m05awNiAJ+RmHiaiH0SNlry5Kta8+Yq{++5rZ_|dXEVSf}9KXFTJAZqo1tF
zSrMhg$S=fMF~P&3aSqFx;Qe0qZ)Wki-@f0PTQC2$^mFO&>u3MHeti1&)yaRpt@(50
zjkV^E{CW;E+tceWsjWEi>+(B={U2`at(`pI?cOQ<@Nf4!%j#~v<-6OxMf=j<-`_s`
z{hGDCe#(4Xzx`!3m#go-<^O&5v13++e0@~;{(pC79r<(njQoFjY(xL0b~7_4voSFA
zNP*UefG5UEGE#F>^}!@~Z1zoLboXsHkvjb^_8fby_jzvGdqg_trr)hMYKpIyJSoiS
zUwg98#Lf6#FK6Q<o3Bd)JiOz%YF|5R@|s+@F#X?*))uqBg0{w!*3SHDG5hqgEvnfO
zOD>k5t0_Kr(Of)sPl^7rbAo=BiJrID#;w@F=XWf}X;S3vy{0l#X01ES>ZR7PiDmm;
z)o1UP-tD*F)lufslKiHj_E34i?cSCd@28xu;W^t8pHMvKz-CSFWR(YIbEm3nTCwdn
z)V;1Xb%BiXq#b@uNd?ulR;tNwc$RvXmp|R(|Ikb@#;hl^N+CA+uv4OXduZnaZ@EiX
zES<Kzb+Q*OvpMp;<#oUjR-267qF<j+E7zZP(J#;Cs95ybyaoS?cHc=TiU_V<w)x_&
z;PapNF4+8SJwr&&qhGys7x=IJ6Jh8*Zg8om``zXD+jiO7AJ1Cw=F+3x*I!&+%@+4L
z{{DkIAD^B-pdIk+{_O+&Y>F-y4`>Ipeb+V2yy65BZsRxm(CfAQxhc;cCf6w{?b0v0
z<gXmQnsI%C=Hl=*AEv$MmD+lzvau;EB7FLrq*Fp0*Ir91nZA^n!9&5dNu@91+((h#
zg<Z}rCpA>h$<1^Ps3@Ldz;||PRBFJ=!}UKuUzys)aqFPR=4pExIPYgns*$ntI~lX(
zuK&T+s~zXrSC@VI@bcv0cz-!ttI8k6_qG?M-iUQP>9f81=i%`B+V49*Ki9v%zvfTr
zwXc~r*FWml{d#mXfB*hdwmmyIj(*Ocx2NLQlb_q;|9yO%-2L6J%gw6BBk9cC>>oFF
zo=lQh&}uKYG(XL-N5D`gfMuWVy5CAh4cDi2+$xC8Y5P#drIpMa7k9^Whn&kxCjsyC
z;frrDyfb#&t+SQw!Kydj-5&}jo_RBOIk#-N@IvM2`-)yYEZ@b8Cxj;zOS;`yP}#i5
zWr3vMhJq)?WgQz+OUgM6AMNr#rujHWYmW7imr6;zCnxmuIBE5)n;`LN(-y&&5b2FK
z_rEQX+}gXVz28S=k61_434TM<Llu&`EM2*J?HbFf4t$@`k)y&t?PbLAuR9xj&Ksoq
z)i+sgzZv77xF_n;`ldrSF7NZ7zOmDX)pE<vJI*~XPsmyGD2tq5JI7P#<qiKIic>U-
z@(dE2-y47Co%nqM`;;8t@3~FACG&MzwmZa>dbu{9Q18An>GbItGv7a+<R=)W#e6H@
zqsnKJ@z=joHvRtJv+1_!FBy}&?|#Hy*)v(_wRP#!DQ)?O|FTcr-7a>Lmq-8fvol{E
zci3}m|Gi)QUGDo&KTnq~yr6LJXjJUlJCE0Lf7e|2#X(2&-s7yd=Mv*)@hk|;)w*)7
zXz#???gcaRg?aAixSl%ZIrB$P#|s6g;)n;~SItU}6<;@gcSf?kbbnTubW+W!?YC@o
z9@Uzi{I>e)synmO9vz&haZCH(F4=EVO3v)PN&MT&)|YKH-Wnb!IAL49lhf-nDuLT4
zWUP@=<X*3|(5*~QCrV7r_hsGGRiR?n7cykMUUT5k$;PK)x4-UvBg3#Qq~lTNncy;&
z4X4)moZFdJ6>YFGj&Va*M5x2Wm(!l?Dy=9giCnqYZ5gjn)XteITa0=xuspxIQhv^y
z^mWm$Uc6OCBCEGATyQaCUw@3Zpyg6tlZ#WE)Xo*&=6z*1r?>awp0|F9_xY!Pl*wv8
z&8ol1*uqA-+oby8eYePGn&P%c7VA9lDlszc=sd8}+4D}>(+}6LItVXS?26vX*7RY$
z>UE~bpW7!kFPZIm?6RgDx9!`g?w*n|Zm#*6tzTCf)mLxWlEUrr-O<Fd%&lbNJ(qx_
zy2V<|7K$=`^{XOH)LG6MUv~<-a(`u~=JeRV4uMw#uJygJar)V_sJoIaTxr*`#qX{e
zUU>J((AmM<*<#j7&Tn(m?`~e6C28qv(rtM;K-4n(*3mD2eLqR9P7*SWNWHXn8}FVA
zyRLb?UJ$T$^AX9-@3MV3q_!?HZl0C0NAqyd=Fef1Z&~VbZ1G@;eza|#@|uYiw`a<~
z@e5x&M<MLBSLlydhui<<Hub)1s`&Ciq{5Lkd}WJ1zh;wbd|BU|lWcpGzPs#mYtc@h
zIB8q;{zv6|Hau%z>UH`3^gX-6HBJBe?Y;8xZ`IcK;o5ui-ZJ{=D0Qu_`*tyvfjzKf
zW_?xQTj}}nkFG}EZZ+&a?Nb}4sr~5HXYZ}s{(kOub|`trzSVK{y;p2Y!}c0$RPN2&
ze1)qyZ>4xk@4W@$zs_EJ(D+(LbHnwO?Qs|Kr!{3sOxv?S>7{w}(zfu`rF-0#-~4XH
z|0$KH$7XWBjhH~((d_T}4A|<hoY*KKX=Vn7AT9<5UdSY4aY<!PYB6YLaZgxm{%ix0
z+WA}lC;XiI>EeM&95bdZ^O(A8RgiV&dN!Y|v+`$VbmiQtzi&FBW!>FdTklpd{7#x*
zz3%<Tt&hLYeeh!1t2WJ<3I9(V?l|Q%)hAEr&r7S~XA-7uCngoBW;Si#a`64-?fH+6
z<yh|BweX5U!vx=^{GOH5*WDH|51H`H<7&vE`CgfhJSAU`<!1SM-n`-2w7c=S>pBIW
zM~sP@UCSFd7VY1%HCQ81_;As&nK!>hNOUGuZky;ec@Fc{bMq%MmDVj`zVdSJ<a<@;
z`x9EjdRLthe6VZPstCvTmt%KuPZwcVaj27<B)==JK>eGay>RBUtg6#FzY3(ny=JkA
zu9~WL@M4C?<7d*lRAnx`PB!o5K04KnZ4*oPm1q^g5{W`5i_@l?4{Iq&+~<3CK<d|#
zzQ-)hsoUJH?~;+Rf8(D2i6#9*#fKY`>!0dupYpFf;o|?81NGa3^xfZmzUS~r@_F@9
z*090^kK`M(+OlRkviH7yoy%J_Ilz_e?w7}Ml6x<vyk>ux(sbL}=W)zSp;g!XUfWg3
z`DC8SOm&c+m7DkG`U<IgJB`@N%A}_E<V;bylT~h>oa11@`t;AKux;nB%D!3Wv`KWP
zL3(P_Lbs@^?`H{|xFPeQt?%gkMZ7nJZ&(KWtB$x5TK(t#Sv$5m{>it_yDeGw)2vzJ
zn9v97H4m0X)P0$9^;hc+SM8!zoBzG2XnF8PwE1hi;OVX>8&zs@r!Jg6D{M`*%BkfG
zzH0tGaMO)z*ZJ)AjsbUFZk4y}ovyLTRPmsA*;z}O!o_C-?iuCHsByT`5Lwoqdr7=J
zx&K}r^V_`)UnhkvdHSsCzwyELZ|+f>!WzyU3!MHYY^I#j!S-)E@19+tdLwL$@pIle
zmLFW#rg8qO+f}taVt!hq@%OOe^fDjkbisAk{4A^^*IvJQh<ndu`N`Y|xcmP9+s7n*
zeZH<V-*gsRHtR|It!@a}yyW@v+lp<=@;l;=Pii@8&!6nGy+6}6<=vUZk1w2bvX5}z
zd^K}MtMI#r#_Kufz5ac;_3$MA>1~{+X5M+4HGSbq(Mc?K?O7haGkPee*Zk5SYcA3@
z-sgT+h=E~4AOnL4G#6zg7Nw>HmlS2@rGqomo0z*rw=JgayT7WEeah;p+cT%-ob<c&
zbMx|)NY4=6=seHD{PVorixM~mJvaPRcCX#qe@T4h-g~Pn8=m-n^L%E_@YQ|7HlgoT
zX6xp@D$EhEv9S1WVZ*-u-?xL`Z}_-u3^{nwuXg>m!!8e2P2@S|ZvCEZ{p#~&5pzBy
z9I+McSUh)UNBO_9Khx{0MUQhH-tAbk#n^Mwhok@M-uShge{3t-#o)g#<jH52dbh`a
zEqYkwj|wXa-M?+W;_AG2zhl^}&xB@9ka=rU_phXVY2zPVQD5Vtih>=*f3&6@Uq9)w
z-iO(qg=epbPCn?d%O)^z&-zP8Licq@&M$oVaiRK$DPM1=$8I~lUPy4`zw@$e6E6O_
zxccvtA6M-fFDr0<7B4&z`-Wr1(FXmpKIy}WDru9{7qFh*TlhzB!pX3JhOilGXP(ZR
zUDYWOIl<5MjhoSu&`;eb0_&FeZ%Kc9`A@~->se2ZtL^!b%l9byeZOm6hF*n;_?cwo
zLRIFRmg71zIgTA~vM>zV$5S!E!A5D`l6_mtjO{l?ueUp0{AkBXVZ$)N8C;f6CaTyS
z(-8_+>QZRY3FO(Yc4_j-IXmV0lWo}QH|*&6;J0Q8yI?`A=)?<~=3g=qlS_BmX7@cK
z>Ok4+^YiWMkA%47w>7ED$wgP!zE1G3T(Y+D^3i?jPJQa0s{V5E9~x$+yDU=?a<G``
zQO?!rqgGL2bxyzNRnOYenCvGl$8Q)EPL6%O^`LNkkBY;ic726Ul?T#$zwNxdU*o}%
z&+8wuXs*2{*zNq#>k&&C$IJ_Ak5pH0=9XMy*!}#TkAq%HOd`8PaifuZQeCU;CgF)o
z5(L&9;rJhR<MC4utCp+|hRyRfaEmUJb7-vFH}BgNQ=2AT=IQJeZ;!uuzj3;nQGj9F
zJVVwLrHH8xP8kf(E6t3Uwk-}>Dt6vqNB0x|H*OKlubHb4i?WxAt`TBe^itW#I*L<A
z*E-YU;)zJTzWYBsxLxKcX0JJFG?776D4%az)6vt{dvzk$r)Af$K5+`ET5G;mY35C*
zK+bzXMhn`MFWvL#VO^P}6gFwnede7!%TM^cy|(JBzW5c+_RvWY0ZtB#QLlTJ7Oy=W
zc0998>Tg|R)s^h{#q!FA!t=Ol3(nur5|Uagt<R(6oiw3aQA9h3uiZRtqkwY%o9Xjn
z6KhN!H*#}5lb`W&gSumYiS2&|J&RtWziToQX8xNLfAE0I0r!J?ndefIrB`gop6L@I
z)iy))@cXWvGv9pN*`ZZ%>xL1Bjatad+S83Dv)9(`uf6nc3;X`4mUW^0``XSKd{^pw
zdZD7#SGt5Z*gCQE%hu^i%hLp!-T7_=Y~XC)-2ZS!{2})Q<7{EaZ+AF^<zocZa;9XZ
zn><k#)BMxWD(U&@!3IOEhGlW57aaH#w?(S?&%Ew08JF+t%3R?NOLmv;TgSv*y`?0f
z^2C0Vy3$?i=Wg%4m%2KCZS}poY?J4_-utXr=wVNqe#3J9-OqJpq<42sna}C2=A@!1
zTBfReq@6o)&6>Y%af&H<vkqkRYhMsby70<BH`H`R&e=IV>)aHVx(S6lKbibMaWiA~
z?)O&l2YLE#Jod`kA$@AM_r<w$4z5XU%0Kz>z$R^F!(|x?{Wly`pImQNIBMVdVPV$H
zYlmKEeDzl22vR>Ky-|MN2LJzeYaD|P%3aSAefF7u&A%U>y8I!+VkKw$ug}<H6mr~I
zY(>jQgD(%#J+@twPd!{+$sj)c=fp6nl}z(11y=^XYI0+0oI5?M^h}Q8_oS9mj^nTH
zhw|@l{@?P>L0z(vv)Z<O(bhW!g)H$~RQCQkG@&qA)~uOxUH$P@Z;ro@i*o#Rbm?Yi
zj=m?NN}QXXM;_Un{4ji*@7{v~AJttfdV<9N^QAuV`^0chk7rVH;5B=le%EPbJ(XE&
z4#%?iJlq<jEK_g8a!k!>^O2%lOZNO{8(+@$NEiM!-BHyhkJWqf4e{jHOAdOt@U_N9
zEB|}RU!l2YY4g^%na&(_CyP$=)*h@1h)YkKS$NRrjfvxrCI+uQg}P5!ez!ayu66j*
zFQ2%iAmr2I*4Bo})AwhqMzcMs&PjOo_(@yC?L8Z2+4o&Bt=YZa?NGDovrHdBRvvE4
zy%ER1eT`cCOy0;taps?ypMv7wY>buoCc*9?$bMhlU^@q&0252G{IUIqZf*D{f5PH%
zOoQ{BE*-&=N}IxUHl?|dPF7wLj91isI}H7vC^dXs@OjSaM=IB{j`mg@Hn==L#cE6H
zNdx90i`JtYHEtyj4yV+1ebcboJmFv^QzVP`tg!qGoy^8J^bD%E9M>xdGfkQHPax>)
z^&XvD>a!GOf<mVRh|FNmb`X6Pd$Q@+QK!J<(?aP7xURl+%lP6ZX|m?Y4vVF?em3rZ
zWW#GTL5KI@vc$6w96EfjTU2d+f7biHzLCk4>4|AQE?kLXg4wQ;9m+2xvIF;LrmLPv
zsh0O@oih1ppzk%8&(+P<Z8L6quZc=M%n|slrJ%Lq=!v8j!zJ>ZyyA7X=l4uLbhOL%
z;;A~lFym9skAB(z{~4e2eaT(h!xs-{mRZjYE57h<p0)Vicdpyr_lHH)e|Y!i(913H
z{I{OxUz=yM$tJXI>f{?QRd_^y@b|DidopkDF0Rh+-xl17?~a@oGtE4Ene&U)3!6=E
z8rofwYtPjznLXWa2H)v<2VUQrU=g!uY7w)<<jf7vU)MAzf8tH`pMI0Qx~XWBjf>;0
z#JM7~J-eGuPVY0aWZ5aZ<WkCU)fSe8%lw=tyxw|dJ)g^6Q?r*Y@{`1lu57;?_V~)r
zrMGQM{ccZR8@K;Q3%h)h<lW$=m3#kld{o?Yo|BK4>052zWJ$kcQ$Cay>-1=d*UUWG
zd2V%#bsBqx)Bo59#|z{&@7>pY;9|CWU-SM&*SRjvmdhye=Tr7v9DZlI#kQ9ZnWZk=
z%$8j#_f}y?m>}=F>=`EuWj^Qx91=QhVx$<d=I``dfv*nl@B3r4^l$v}V&%e=i_iR*
zdw+b!8{a!i+0h|W;YUy4tKSy6t91<hUM|p?z57V<)cNrne@eaFoSp9a-<9kCf4R5o
zS?33x^4oEuCuP;vV;}M-{`hC^T;tWOoLA#<WPYvjk#@%Vh@$iRB`-Z@7C+~r^8FOQ
zsL$K&*<b#*zOxNDIjcLe{9`s_`$w1A!PDiZg!7%gaVC1X<lB9vvnC#Slqq>o<M5ZY
zF3v5l{~vdXm*LJ7yRR#>%l(l35(kgqe|m`*YQ6qAb2XVfd@sawSN);9AH$N9b<;mQ
zw>$aoL(}}<-=CVd+hzAHw+V^dJaN0vF5MTAllMtG&kJ?r4ZU8&9JwHr^SrI<8o2<a
zjH7RMyEuG5*Vrz?x;CO<dccc&ksl^JS}4l>PDtdgHMg8!%9f-9v4<8{ryn@nJ@KFG
zx&O}`G~X{d#k_Ej?%lL!xrM=TW~US8eLmvII={Gc{zhYyjVH9)Ta#`dexzS|b9?rx
zhGkr8f4AP561ASiDtqw>u6yOacXz9=3rJWt+wg*O!zVK>tHrxkDXcDVUbuO~wL3>X
z-}$v?{eHQ11xE3GftrhcTktT=|99ejo^ZzfRRO0Mud>gKcF3{{df(xD`l|Q=UCnPV
z`>iMDw08OGz0Q7eY}==jD925iTbfcPA8#x6J?hR^B`x%P!?TJXKab2xv77R0X3s8#
z8~dwIo?d)dJ?hAr0~VDh`PTkA=AAD4N0{mVjd>wF{|a^a@44m3{t0VnTxR>;X8%LZ
z_18Wp*IudFq;QqhC?w5RFs~`8!n?sUyL;(c7N**UxDUHxc(yLP7x?Gg^c^=JZD%c-
z$tM0yd5y#EPdivMbXTZ0e617~>p4HSH~*rv+Q-h^c|5BfVt8YiZ~U3Dm#JXiitrtU
zGe7srO_cs%8Lzk3XwUXF`&-54ZC)R6;B&6Xg|<Dd?{XdVF0ACu;l8@==8L%t7Ob*=
z@l7f}s?v2<ANQ(7vuvb-{Mn*pE}SXIHuNgmv$O4T?fHwDX)C|&HojJSw#@qR%;qBR
zD+{fRzZ_QCb^L;(zckl;S&_6?PMJFYnP2Y^eARZLv%TMn>;6)%=P{b+6=#-}Bra~Y
zSDm0YZ}s0Vzj9*tPFECOv7qYm?euWLUqx*i8(nS+lrl;!nI>BDCZxV{OSOfW#Dh88
ztIs(~=H8Ro6{Gi(NsA|=IPJB@g>xQvt+E$SW~t$wVjMB`q_tL2PqGvL#TT_w>d^}}
zn;O~Q>F~Pg#U1kbbV`s~P3Je!lKvS@E?GYocu8vQf9q}d`_IH#VysH){jyhH^?#2z
zankV3BC)gHrL!Mb&1Sf&mv&L$iK5^~<!{g4O0zABJ}td%jj+(s)QOW{cWdac+bR8z
zZSJK-65*8{Jc~bm$arr0ul>^FTk}H>gvPAcVr#kGNoMH=yI-CEw=CMsE0<R{#Vn^l
zs;~O3<rT*4PUB+%dcJ8c(~}pLUb|s_*}|)l(cN8TSNh&5Qp;D?FwE6+RXu$(%5`<h
zrG{T6S`U;LbSJ%Gn&NlP=WC)|=;gmgPMjyXwm*_P$>Wgu^U+H-p4?|MHWpu)&oTA!
z>9Uw&B}enLq|)RyE2|gYzp=c{P((B8IeXQ_RJClCfS$&8c6S+OE(}Qh8<xFr2FGcq
zDIL$lH~1uW{yj0x(R8k%>@MEbho8DguHMMLInX*~_q?FswzLyIQ#OU_F235xePaEW
zTR*1EE?;upRNeO4ubHy5Y^B0?99r_vf?wv{g-H<s{(q*;diMKQqhU|9!=-m;o7T%}
zpGxYwv0$Rv^6<NTdk_78x$TGa`tr|WeB19eWpB;Z=i02e!FTZy5$+=kJ}TZ(JUP)<
z(&od(we5Vz3y)kBbK*O^eU6??mE7B{M^7AHeDI>q*EfCMbM?-%u(rQl6F)a((~Bjw
z%Qt_>df0U3;H^WhmCZJVVH?(zoO~&@LYM!i<%^?vFLU?5{bJ>}KX2-<{jU7BGB!~u
zAN;TVS*3r(!ScUeT=kM*Z9{j5pNa1lUX$ZB2=x!ztUZ~rXk(|tVP>s-&7>Vpw+-*k
zPq=j?Zd2nNKh2ibd#^Su{`S~=T8L#oljYw(h6n#^zhTJ|UMSDm?v>SewN|Li|8ElW
zv=6FRIgh8xI-XJgQ~T;P*YEzF>K8b6z4<)7gK>@CFV;PWZHhVKzRX*Ev(4Yqs_=e=
z*Q$wqkCZLL3N5pyxq7X7fAOe=Z=o>jI|~8F*1xhdsyfszDo<LJxJYK&RBMk_!S8&R
zyKFlz@4P7b@chNPa}tcro%VaZ6q|o5dXAQ*_wo#F_8l7k{s{h0Sl#T=Dyi(YgNL>C
z?YB~~nKFMPSd6)zhEy8I%qlSSe_;0ExarQ1K^o_|7F-CGZanVzN78EJ<0qewf8TQZ
zna!MnJ3k&AoUz>SSBm%c-*?zvee9g|DPZCTSFUi+9g`g{_@&mp3-VvH#XsoirGP0r
z#oOlp4@}G#jHt?Z?8Yy%%h7w|%IKMr43m;psV&$u?f)_d`R6P~5B(NZ_AE2l@@6`>
z+g=Ao8`J2?-Rsy5TVoiMc0F|YDaP$`-~VE?c)r2JbvHcnLKJU{s$ZSp|2V|^pI1V4
zk&e4$THca^$LDx&>jX3e_!h|dJxjD`_h_hz;*Ocm-O9IGVvj>{Op9>9&I4DuOB9vz
zShA~F6?UxVnBHCCocd+UoZYTf2ULBId-;k#xG<Y7B>m|It}To{J||k77AsA4Fk5t`
z`A6^r+vj0M`I>uJwHR{MQqSic*IZ|OV5)HfL${!O<%LT|z7`r6SXveuyl!2_(SPB0
z+2l*|A{i|Ax-OY)IH_iH=hdS}ns3zWbX_+&XeVVIoRl3AccOJ~&rZKJcM@juaZO?V
z9=c-Xg!h-ZtdCx7(VO(~z=fhrYmE;aJNjPiZTitvF>hnGbK)WG@Frf*pJLo5OyUs_
zx2ZI6J$>GP?C_sM$7Rg97c@_=(_9+fD0k-B{QrVp7dU1uk@aamEW24-?(w1~>5#7r
z9N0O;4(Q07&TI)`F5<}G?z<q)({HsYU*|xsoch|h_u<x4zFc_PnS3LLi{G5<&ZUd3
zDpB93J3084{!ML=n%h@=?Zw|KQd|7yZZ7}3eQ)0PGrJldhC6x7*IkT!v$#`cBcE*j
zUz2$XZ@xZ^(4F!y`zMQIHy`&kK2x#i1I$a_F8p*YGv#B0;w`T4JcV&}zmmOosDFrM
zwK>~1^MJ0`1V!JEZyU{ar>nm$zQEO{(Vrju)4BHH^(7YfMC7?nd=YBnbLSHe-*_-R
z>_YG_>vx<BH|O{4f5F5Z=qth$`0AJQna<lAw0Od6A9lN$@`(y9-07|*5iabo*m;Xf
z70>m(3LFPy^}_GP+^_i@t+m)Ag8e$*y4jgu3tonnJ8W`z+bO>GnZK^KWSEd|H@{Vv
zp6l5~VUHJG%{IF=*YEZzllyy)SRMQnP?+qu`g-0R>-E<AgV+oHKQDh9xBlwJH7=6f
zX95$vXYt<q7NwhhVK4JuF3X(6+dr(|#Z+B>z3+9s;@WFRPhMKD$6><2-tcch?%j_?
z-1E(!o_W;#;e#n(BbQx_vLdhjn;ZAvzp*UabbDTM<(z%C>1VgTeSGZb<9a4+gUHso
zM&1V*7#N%w85ooy!(X7?Pv9*o#rlx#P+$&tP87Tus-^#6*C7KD*YC{#oGYp)Uu)%Y
z4mMZX5&E1t>DMi`+tVj}7x`PgStj+s)29XJpPJ47ZL?B3jc=mPm0e7&Pm;wmHhieQ
zrMus8qs6N`BAmK56A~}$?VH8y=jYXVd_sgM_r}&hmA?{)>yF&q*vGkYg59&5o}HBj
zCx3QJUs$s$EB8}8`xbq_t?zZ7)P^@(*#<@ad(&4dZ<Dup|2dCO=UY52CnVS{@~}MN
z`u0;plH=?0Z}YNi`B;|ao;m&b8z;69&_5NBGmnvhL7ka_L6CuiF(p4apeVnfSU)+x
zC>1o6uCRz@bn_Pj14B*FLBGca0(;kQ`7iKu&&s)OsohuEI0e)vge_X_dS}nKfITYP
zQvEr9{w{qpX&Xo3=QGux|NOC8-MxLCZrnlk+_SSnjNP~lE@|a%%l=ed9IWRl7$Bh<
zm7U<r#_?p)&wsh9veT~}*3Xx%xFXQCB86#sp6@(C+2tlD(qFjDaBb^+F!}4;0};P;
zOP%(0=w7f)PT#^T;<aGG;)IuKuA8QZ6iq6zzOTlps=4Jt<3Wejk*@qPEcaXe*@cRZ
zyfNX^+n1UDJtF!fkKdc^3)GYDf8V^%&ur2EF1tH-?=4?G?N(2ur@~n;&pxT>jIb{S
zJ<B)z{>`uQ>|2P7LzO^AE!#4sV}F>O^Y?NITCRN5)1<JsOu}CO{9F^I-J*=|ALZL+
z226J@zRA_2`Jn%iR;&1=C+wGgY?X`qx<lyaIp?J*T!yn*LytsM@O_m185|#*pXXtm
zesxhwx@vUJD%JIMpDe3o7F*1ePd?R~m;WjI=;uvIPiuB2KHn4fVV7M{jd;eTG_IDl
zU!Aw<FW1stdvJ&8ca4JW{UNvfqmHk>`-dG{qL!Z>+x?z_fuW0$fk6N{Q70A@jHYSC
z?#($9eEq>28a@8BYkW7A^bnqP$>*_mjG9Hy-NjelDn(g)9lG+bnn#$m@JZV5t@8V~
z-RyDJ_nT0pu;%HeOF;}(SxjYHYwp_X<t;gJH1Tnj(~74ok$G>PwpMzIC{AUFG+-*$
z?O%FiYtT*k#Dbkkl98F3yz2^PS-Iv)+uWI4FVE1u>r&Oy1f2$6lWq&<r!#lhP5tjA
z{b~NcPl@)fW;tigj@ea&|M1$)Z5Oy~TSxmvXZ?8_`8nMFzB9CLU)P+{J8QAE)S|<k
z-o;;z{5Wzsf4d$tcYTqX+`08zPlqmDzxCg(PqD_Yo_Bu3y8rcZ5%+aIMh1pSjG#p@
z$R%HCaY=seXgWboC#Un84Fp`jPyNL|?fa2tVIc<xmKh#*ODC4~A5ii5p?G^?=c<3z
zy}wS_?UuYdbFuN{x9gLfjLpkJjE+C#GI+Ep_^*b)Wu3KLlZ5i?sGHj-PqH~~p7(KI
zMct{m-b;6OOpJaac<I@XH*@?rH{PAueQjNhPOXuKa@6zaVx=sJV|~U~V_By61vKcs
z)VroB_@Z9s-@~Pwt3__E&fB>8{qfn~cDa7z-oMas$Hwy?)w{EfDRnU2QYjU1ieNoy
zb(2+q`9jyV_9Z#{3@)8EIO?(O(EK><mwUc3)xNYXEql4devL6}%t2#@S+dMx$2?uG
zZ_zTJo%Q<lu7#owA!?>6PxeioXun!X`Fpz9eFki0jm`d<A3Yfu8167IFsMVC5Cuh<
zc_pbu!Ko!BnR)5O`e=MZy`;>%q=W@L2Ad2v-({3KxaRx`hNks9K3N7enAkhB2n$c>
zOGsjyWXH3TN7Z4*sQ{OSwF+;00*oC@gs1$M=eT3r>^q^C9jj~VyPXdoU|?WS0w-z>
zhGcMBj>yf?&n!vJ1-D`p7O603EKoqiE)G`DdHd|c4m|A-*C*wjKbg5e_gGomI!=*Y
zC8FtbYbQ+Zn;i7~@piqx=4F>Z+1X!_&;NO8PvBJ*n;m=eW1VGEylnX2WfgvZ(U%;*
zD`;uJR+n2fB5PtL?<9p;G4OJyhRhFMu}J7F`$+}EWex?`w|$>*vcSptN9Y~#;>%x{
z*S>yry?MG-{!O{K^UJ6ES5J@NwY*Rsdg%W(VRy}p)6Z-V7>6aV@;?;zbWz?Lb<XQP
zIq4~LGIuAudH8<&V{5Ddvslq{-xdZ2hLa4SJ;~?+18(w=6goV;r@0P+vP9E=fwcI^
z$5ySCaOG{PX5FqhH*4X{iE3R6|8_^-6~D9jX6+$`dltv$opi7aR+wP3-S#eTjQ;XP
zY=V;)U*q-3Xpkto{W9h9vC`L*x;x!%)OnP|Yp-oK{a|}#dFoQ_(1=={?p5o;T<^Yk
z^Q<af@5-K>HGdmVWSla;ctCmS^QcoV<RxO}d=%E~m+e2eROa~ZC+l9@Mtxh~Fl9+t
zmgZmW$A@)!6t>r64UxSqNy-x#7#MCdfXYAg5HW(3f5at~oIY#5Lk0q@>zgi3)_rl|
zkx3(0;8Uh!7F|o%T>qpee_xK>v2jY(yfnEol|7|0SKVd3KFvKSoGxXv^yaGRXS<(i
zw5J96Rx=b%R@c^LogP_n<(grmkJ;}didr^hAGT(0Js1&tN#ad&NW7u#ZtJU8EMqpG
zP4!g#qac^>cwJq-N&J7zac+kn)_M1^dUV<HwvHVPc%z;q-}3Y#^#Fd)$U5$$?Y6NO
z9QU<Xoj%Vq31mA`Z@&8UUiNo4-~2nQ5od9<@#KM|xCVz$cX&&T3oPYRS(Vh>=St7K
zz+m`h+vfSc;@9VT>#5EAFvXz5hR4$X+JWlS(n#}M&DYnmey;w`9Vy?}oAWmQ>XcoY
ztD65coX7%2zwOfFQ7d1(mx!73nOXDjyyXvq)?M~3yw5ME`#asiD=ShY>Sy}3_G}5q
z-2GTXCFCHpkOVse!$oCm(Qiyj^q-0H&cAITuuuF&y}*{|5!2*5bFXMlvbf88Vavf|
z5pUKm@HlnNO(=Ad*PQzO&aZT%w=GE#&*ke)Jy&q=^I47gcJG<)S1gOE(21G#Eg?h7
z<Va7&!?es}muzRxoUd=6xbuL{6q6~huXorz(&e4N$X~DHdzmBPU;WDq<{}rC?o6`a
zYw($`Ra?D1YS!Vq7ux1Oe0wx&$ENhpFEea*?VIJZ?)<~=vu|JgT41qv*8K*T0_V<e
z>wRv0{$;MMW$;~TW$Dg`BA!_a%1iXwJfBKj@su;{cJTXaVcI;6xA(*!nRtt0r#F^I
zzg?N>m%gr9{mM-L6IV<Q<SmRy*lDso$c_87&@HVbfr=T%2V*P`=xmtrx@tz$;)fsZ
zaIJ7*@l2UI<)j4b3jW0%&0k`wLNcC(2Nj*~R26;p^f>?gvL+p&l{}q4+Pt0$3o7z(
zm5cA%;Vs>5%(ly0F8J@ouPgjpjAYBwlGBs=KeDP#mYh4GILZHHhSgl_6qA{j%LKau
zg$z4+9nZwfn?8N|!P6B7ez?s%v+du$|Fc}8)0WT8s5h%xJ*%@u;oYS<3qQ&%=CpOb
zvd?06&&fYI&u;xVbJW>pwzY9PPi9_;euu<#w}=_LMR!iM?EMv@+MBk)#JlZ0pF@hP
z-QATbXJmcrk4@P5^SM^P&cFA89cHsC_1-YXo%sHGQNb%c-={mOy==aC{4{R)S}S!y
zfHCfe?Bf~D<!#UA$+36tS5UlvB!TtiUN?V}H>H0=b#6-Ky68W;pXp&4RGRjdRf}=g
z>iXY2a{rh+-|{N6TSjQr2K&FbpWS!qe2~6&{>$@=L|ET;mOQlzHH>QAFp=v^>&FDw
zSdrfWDs{(Bd|%)&L8<fq;vRL!KsmqbkN+;&ppqie;xg%nU|zVxo(rjhlcGa6>8)Wu
zqUv&-_lx(rKZ}Iki??1g-4S@f=JGE?d%X&0%Pb9sjEAi)?-m3(HGZ#=x3t?bbC11%
zzl&X!Z-Aq6_gZ%jR+&kYD;caRI{pZkTzGZz`-OU6&CdaBn}2MHQMDAkptSK+%wb0+
zm6i#SHOZ6ILbW{?aYfpm4~Yuv&W@ZL6?%HD(fhQiwq~+Ao85oEj#wVEsrd4Zh7Y@3
ze5QJ`9}eNXpE1?+?#c45nL*lGJxTM;yl(9i@eQ6967k<k>iyUAPVWP^Cs-U<>Fgr2
z;9qG^h&xAX@3Ff%Th6a&?)=vI-=ACc>kd)7m8HQFJx&>q^%w+ITIXpbe5z#nvfI3;
zuqsS0&(A?(3U8!nk8r!{l*d1K?4N#kZq87sx#Y%S;lE4Ym|9mZZ(nY#FCuhMVd6rC
zpalywp6%PFaxg4QHSs3%ZLM=oidu;qj2i!`*`H*%JjWqCGIa9{=f0k>mG}EkOj9_x
zo+tary!&2B$G<0BKe+Om=Fx63iLeNJM`wrP<ZO$LZOYvb7`OaQ+K|~2-7Zqs|9b1k
zRYoC(9YK#3@3YNlW7_2XIqmP_Z?1f=ye_*~rN5IuR}ga9#^VFiYF+NfTLn4H%y#F<
zZIOQF*t9i>ligH0q-{l+jF{X3*Oh7AeL^cP%7lp?^?Uu|>q^tFcB~JQcU;u&7kRJt
zdnMPt&>Qh*-*Z)(W$m9Z@5}4S@u9g}e|Y83e4(wqeU0hghl$cQ%1O`VyRJT4Ycprt
z7Hv(l%)1gQZMBQm>-S9lRS<hwuk+%hQ0dZ-fjc@EmVBScx@lgXME%<r58Ga<-xb`)
zZ|pbUc5B7k6T1GV4N`)xZ+YWsrjf?q+gAU2t$I)TljBz!eteqJ*`m1p3diFE3-2VX
zQ@mholYHkyiwJWvuX=#U&n^ewfCDUSt2$Z_HD{i0h|4}Am}8;n+z?&ITBLZHx9fvh
zy4jB(4<}@t6k5E_^^-{G7F*k?KeP6;bmwfX*ipCn`&TY)?MTN{+rll>wmuY;5&wRw
zKdgG;lIJRkAD78~_$YdCWvM^+nlJNKPb<m4oaH{}*0Ni!@(->uElK`$_g-2_VBGq*
zF8kOn929(bbX7#r6{p8RM+?gKm6f{Hwi<RjSBB@e%4U=vUMcp+HcQnb*OwzdT<iCf
zw8+Vy;w~Dme#E^svhiAl&z+V%OEU5wxm<qt<ofE_+KLk27Cf04eSq)N*6q6+!%r>T
zdHmQyW<S-2d_OCuKa-wwi+J7G?(Cr|{zN}7Nbl!PlPP@{x2jdvh)!ZF*ZdP6zD@0{
z50~%4PQRm9P8w)`@laa2u9ok}s)N^Z`mAR2NUWC8Xn(EA$DVUJc-pR&QJYp?TdDNv
z*R)kfHsxzaaQ5pY{9B^;^c8PN|E1|gN*lhcT5#2Dss7^Dqk;KR5i7nd3q7%?XT$Cj
zJntv=U*6}u(eA;^cd45%Ufmn18UJ0KE&aPXTmE--w({@lZ1v^pajU0#y;obh;{5~G
ztAcB8Ep*?y>e1DxsoWQ~hiJ0xT3vG9xc0lJ%~w9puoE}V@8|k5cgdBaXS)MlUyy$q
z^N$-_N5MLjC4xXRgrr&*-VAB_FL2LJ&F^y5l9ea6uzojwV76p^RK~r@YNrl7*gyN~
z(|gR{O_uN9&>_CAFm2O@`<VhB^N!EV6RWFdi)`*VdGUpq+7<_%<z>4wk5$OT78!o*
zn%?60B=KD3OPxRA_n*bQSTv<JGUU>fQr+{;_fI;<%eCcQ*lPHUS8288eHZTUQ{>#2
zF4~`5P<Whw36FYQv-i$*=0`7^RKK}l&m0sK9~FMeezRT7VWw4ou{J|W>@BMfF)%P_
zFfuU6F>o*<IxLXUS9rHYuP8OASYdQ965LyBIe)N=$&iQb!R`mvQ|`<?S0OO*Wksae
zp3WT)ww4$fAI%J&rSa>Vn7_8p8$-^Y|CZdfKeXoEjGMtG7eY)mditv5-I?v5zrQLb
z(%rjvuI8dfb)z@CpDwa}n6p&p;}hZ4i(JeMI+x77_C}3onzon3>68UBYo-|WTSb;{
zRen_(n_u!T*Wr8OD(8F!`Cn{Te$Mpz^xurPx%`LgJ28V(@`AVSMrgzwsyQ=d#|g83
z4YoZA729X6y2HII`|I)gyIAu^%rteS!vs<RxQ#xV3c&FXAzIEK%sOnq!}g$h;roP=
z_&MJUygDy?OFl7k<K8Rf?&*_!I>IaSe%-R=Vr<r4PWSaq+pA|><u$JEVlxZn<yyS>
z-tmhEet!A7k@;v+a(;SO#)C$l_OiWMbp_jkrd53E3=LG0p5bySF1t*a?X;+BPpVM>
zUzpcS_ql6!-xB<~ZJm8-o#}#F<IwBUCt~+AmR2p9^k={C;RmvECv_PY{}Ehx_mI;=
zp=a-XI=gc|cO)I)ZcEROo?7j=y|?c7`oFB$ipc8EOp6I618^BRnhYQ*08Ap1K_WOA
zlz@`KX^kglZtT0I+zZbnosM`F0Z9h7pk(m1{r-+s@;Wjak}qb7iKw{;d)zzl{N;~5
z%tw<hr-XZDJZR!+f0%i>{z2AKy^lght5>O*ZRmKi_u3n=o@vsjX1q#SCbag-ri!_1
ze{Ylfx-B=pbf;;;4&$p9k{XY%G+M<L`PSuJwqbA2e{xQ6@*m-acW1Uf3Xi$Zow{z5
zuG~=vfu8Hvyo2Af+J5=<ZvR7mSTbN_(q+JPfHl~cNCcMStOLAJ4M06D7fCO}6=nto
zNYw&PMCcmP&kRLqWoHL(lmyKW!cGuH*Nl3kB0?(zgDMwTGsI$KBd|K60PO%hgh};0
zU~@nwfafKVkJCeU4f-xmgeeM2U{k=mJ(2c&q8os|#Rg%(6AiEd#qbR`=vvYD6(F=~
z=|Hri?=nC)0exx~VS=VT*aW07fCMWzK){1|SOZlbWu_Bh90P;40|NubyeGPQQ0J78
zwKF&)X@|}*q3cGSY(Unn>WZWr()t07y`UR_KFo|TBF7yZ9B6~gSO=Za4MiO-Mi|M!
z@XMV6+t@L>0gF*aRgn!4_Cz-TF~*8+3hGECvMH)w=%#>M0q6#xj_@HHz~_x_0N5E=
zO+g(ILpDXl2i+8KErjk2)D}Ln0aJs~4FEd>t0|~0Vq{ZRhaj5*N%W}gV{{jxwyTiM
s+8Bo2EO5wR%|ED>GqPDbBCwkU_8!zMRyL3hF$OV)2doSX`=dZS0AUI4>i_@%

literal 0
HcmV?d00001

diff --git a/quick_summary.Rmd b/quick_summary.Rmd
index 0576616..6a0e525 100644
--- a/quick_summary.Rmd
+++ b/quick_summary.Rmd
@@ -16,6 +16,8 @@ path <- paste0(getwd(),"/")
 
 source(paste0(path,"main.R"))
 
+library(tidyr)
+library(dplyr)
 library(ggplot2)
 library(knitr)
 library(kableExtra)
diff --git a/quick_summary.html b/quick_summary.html
index 44d6e6b..831f982 100644
--- a/quick_summary.html
+++ b/quick_summary.html
@@ -385,7 +385,7 @@ $(document).ready(function () {
 
 <h1 class="title toc-ignore">Overview simulated data</h1>
 <h4 class="author">CB</h4>
-<h4 class="date">09 12, 2019</h4>
+<h4 class="date">11 12, 2019</h4>
 
 </div>
 
-- 
GitLab