From a5b4bc3301101505b00e7c97a9f405601574b273 Mon Sep 17 00:00:00 2001 From: Stiofain <17852477+sdfordham@users.noreply.github.com> Date: Sat, 11 May 2024 21:00:24 +0100 Subject: [PATCH] remove R code --- R/basque.R | 54 ------------------------------------- R/germany.R | 72 -------------------------------------------------- R/germany.rda | Bin 18990 -> 0 bytes 3 files changed, 126 deletions(-) delete mode 100644 R/basque.R delete mode 100644 R/germany.R delete mode 100644 R/germany.rda diff --git a/R/basque.R b/R/basque.R deleted file mode 100644 index e64e178..0000000 --- a/R/basque.R +++ /dev/null @@ -1,54 +0,0 @@ -require(Synth) - -data(basque) - -treated.regions <- 'Basque Country (Pais Vasco)' -control.regions <- c('Spain (Espana)', 'Andalucia', 'Aragon', - 'Principado De Asturias', 'Baleares (Islas)', 'Canarias', - 'Cantabria', 'Castilla Y Leon', 'Castilla-La Mancha', - 'Cataluna', 'Comunidad Valenciana', 'Extremadura', - 'Galicia', 'Madrid (Comunidad De)', 'Murcia (Region de)', - 'Navarra (Comunidad Foral De)', 'Rioja (La)') - -dataprep.out <- dataprep( - foo = basque, - predictors = c('school.illit', 'school.prim', 'school.med', - 'school.high', 'school.post.high', 'invest'), - predictors.op = 'mean', - time.predictors.prior = 1964:1969, - special.predictors = list( - list('gdpcap', 1960:1969 ,'mean'), - list('sec.agriculture', seq(1961, 1969, 2), 'mean'), - list('sec.energy', seq(1961, 1969, 2), 'mean'), - list('sec.industry', seq(1961, 1969, 2), 'mean'), - list('sec.construction', seq(1961, 1969, 2), 'mean'), - list('sec.services.venta', seq(1961, 1969, 2), 'mean'), - list('sec.services.nonventa', seq(1961, 1969, 2), 'mean'), - list('popdens', 1969, 'mean') - ), - dependent = 'gdpcap', - unit.variable = 'regionno', - unit.names.variable = 'regionname', - time.variable = 'year', - treatment.identifier = 17, - controls.identifier = control.regions, - time.optimize.ssr = 1960:1969, - time.plot = 1955:1997 -) - -synth.out <- synth( - data.prep.obj = dataprep.out, - method = 'BFGS' -) - -path.plot( - dataprep.res = dataprep.out, - synth.res = synth.out, - tr.intake = 1998 -) - -state.weights <- setNames( - control.regions, - round(synth.out$solution.w, digits = 2) -) -print(state.weights) \ No newline at end of file diff --git a/R/germany.R b/R/germany.R deleted file mode 100644 index e135dbe..0000000 --- a/R/germany.R +++ /dev/null @@ -1,72 +0,0 @@ -require(Synth) - -load('germany.rda') - -treated.countries <- 'West Germany' -control.countries <- c('USA', 'UK', 'Austria', 'Belgium', 'Denmark', 'France', - 'Italy', 'Netherlands', 'Norway', 'Switzerland', 'Japan', - 'Greece', 'Portugal', 'Spain', 'Australia', 'New Zealand' - ) - -training.dataprep.out <- dataprep( - foo = germany, - predictors = c('gdp', 'trade', 'infrate'), - special.predictors = list( - list('industry', 1971:1980, 'mean'), - list('schooling', c(1970,1975), 'mean'), - list('invest70', 1980, 'mean') - ), - dependent = 'gdp', - unit.variable = 'code', - unit.names.variable = 'country', - time.variable = 'year', - treatment.identifier = 7, - controls.identifier = control.countries, - time.predictors.prior = 1971:1980, - time.optimize.ssr = 1981:1990, - time.plot = 1960:2003 -) - -training.synth.out <- synth( - data.prep.obj = training.dataprep.out, - Margin.ipop = 0.005, - Sigf.ipop = 7, - Bound.ipop = 6 -) - -dataprep.out <- dataprep( - foo = germany, - predictors = c('gdp', 'trade', 'infrate'), - predictors.op = 'mean', - special.predictors = list( - list('industry', 1981:1990, 'mean'), - list('schooling', c(1980, 1985), 'mean'), - list('invest80', 1980, 'mean') - ), - dependent = 'gdp', - unit.variable = 'code', - unit.names.variable = 'country', - time.variable = 'year', - treatment.identifier = 7, - controls.identifier = control.countries, - time.predictors.prior = 1981:1990, - time.optimize.ssr = 1960:1989, - time.plot = 1960:2003 -) - -synth.out <- synth( - dataprep.out, - custom.v = as.numeric(training.synth.out$solution.v) -) - -path.plot( - dataprep.res = dataprep.out, - synth.res = synth.out, - tr.intake = 1990 -) - -state.weights <- setNames( - control.countries, - round(synth.out$solution.w, digits = 2) -) -print(state.weights) \ No newline at end of file diff --git a/R/germany.rda b/R/germany.rda deleted file mode 100644 index 059630880c3c01ac098c9b8eda08c942dbf33861..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 18990 zcmXVWWmp_d({&OEgy0ff0t5)I!6iU&0>RyayF2U#hu{t&i@W>c?(T~_yTIbGxO{n^ z`PfqL0*9ST%YB{B~-d!#m78}h`*Z&B$oaB z6V>CjhbMD%IP349)5f$k)UrX|$4c9s&a9`GzFFrFYrZF!K#Gnkq=Bi)bla-->t)x; zFED+0-v!}q#>oAs3qLt$OjGe^8O1>V;rKHO=;B5~Z%TFGJDI<6Jo z=rof|U);BxzGZP&32eLw-nzb1LIpAv2;HhNZggkK?_oWwm%5)7)%yOs>n?b!Rs!kH z8p))1EIFF}Hvs$_@SJQa?`8kP{qZ;}a&h{H8yCFg9)2DMSJlm$dvO|3ISh_&JQNgH z*$$SygjD4>xmt{4Wf8AZdxXE3ely%ps*Z3i^a!Weej!%D5wy2)y4}U9k$lXekrxa_ zRVE`_rgoPV(>UZ2c3ZHcV7S*Qjj8=V$VOs z@U6>9&ZV1lD9TF~cS=m|KcO-Cz3=GykvDVpIcyOrDeQBsN}+0Ij<9mNc)Fy{SMTWB zxZ1SZsM>_u(At2D*L(OQu%s%+4K+Za9evc!kKbzYY96^sMe2p>c}Geml1quT*ydvL zrX$@+DM^h<*#tiX)vyYlqKb~D*iFI!XhK&C+KNj#Mm<k^|6KeZ;Ga40;-%7PBG&yn2si)H7=Me}x2G2x zEK_V4OYuoq?jg(OyO$0A*1>47bs5@Z8=WFMP21N8*AhMxi;pH)?Me=;CeM)K)3%xD zAuUJT^)%@b#~Jwg#V3}Nhjx?IocqtX!_VFm{d8BN={ZJ0&7?Q_u82gcS2VK;MqfU$H6| zUVELfZRL;!%j~-fN@q=?*)Wvs4r)3rCv!x{#Ma@XO**5oM|~X2_DEZcxowNv zo*lQp1vAdROaUCxhnGh(J(oX&iEeN$W?8WXo!klJ>=Cqk1=i}oOGb8_6#T)kR_<_pT&274>I-3kmv!czL$!RscVLgM zkg8$Px5LvSpWo>c#Nbd40IG=uk_k=2Wxi&QnN1PnoD%k|M+_CU;yj)32jw)61P;*h zxwlGBq}Ds`pJ;UlPf#W5r?H})T^GFgs#1JIHM%`WUW8v>(iv1u6)>7H1jfcC>>Th+ zDaLEgj~*&OjWCoS9riZ=5WIIKuzvhKW}|}ErQfhe{t!}HELM8Wa6C``@fd=FOOKW7 zlkv^#*Rmy7G%DkqBUSF6jKsUx886qLZ{}&kxDwBo>8psKQ3y3Y=v>5n1nm&q3IqmC?EDfSe zj>l2<`jfEi@yQRtt!aZ`hN0M28h5U(64zjBy?2wExS4sATe)vcQ05hQqkWsv(JtF% zf$R_Xh7@1D>(_#Y3MO_RXO{BY##mrN zqvpfsE`O-jSg8?A#r|5ApAe(~R67@ZUpuuRr%Y^>#gVbI%jg&Yl6w_5&I{c86iEmN-x~q4vLux4jNk zA`ES~(%P~J*Rqr|Udah|N{N#jEJQCISutfs$@*Nlex@wRR~@;DRx)vI&bTD4Bl=^_ z`@{=hlwb|YB&qy08^II6?M{2)M_zAXvX`t+`fDPadf6XEy7y7k@X zNy3Lpif4n-Q1&0EhOz{X+}H#$S?>}CMs_vYy*<_4whK3x$4XiOL-8kbDU$e8l^4HL z+W0c1xE&2_S)udMV#*RK+3X!rXOwx`anHn(4Sw0dX=+_S@I3^7?xkqvanp^%WEtB( zhE7=daye2um}IW9&bnNzC+;2J4f-rJoXJiJy?f`q54|`&{rt%uhS$?DP0e?Chh!eB9X;w!hko zOAIuK=;CvY_ zg$p!=*$B$fYn;2vcO^#R$9I(TdJR|thRG5mwrg_UruR2URx{Zyw8^LgKp3tdW$g5W zh*FhnU`(0%pD<9a^O?HHJpMe4L{Mo~2tN}j=M`}%5Hzw>57IZ4C(XLB`pXcTJimcg z+DX%sVermz`timCSjw(Y%1#PVu+nGQ!)Z>M;2SLloX*@y_LI?QT@0O%MCZ5FM%8;d za+$L-c6x@)RZ)J^o(&_f%5kpd&6q2LOUFQF@1&6sunce<+^*p6*u$nss04JJa~aY; zH#t!LF{;R$O-@0lQyk_TX8kBAkMK#z1lovfS^KP$v{c3`T^?V$0lXq6Tx(@;p%b5; z#*E~5N0+~K%%l!!S-(`x-fdv|fO-J=hSPjBL~8!Y$O#YZYjWDsgBm0x#h)!@+gtH= zy`-pA4aEl~rM?IlmMD0!N6ocD$ODUhsvxnTct7+pqMLGot{5Vu0O3RAZoIp z?P3eQc>7JFU!R;nJ`+|D-tpIzfbA`ZfH0w3=2Bqk`O>#vyyg(_rq+^!lUG0{JOJCf5_h0M!zaO>N$0H>2NrB!-F$0<#RsA zPCw&b6U6f(%UJkW;~);R9|YqR&cOVa1J5uT!rW|5VK7|7v+Vqxt(mVq7^FsZ2bV(y z-ZB7;cy=LN!wi4`Ly1Fte@o)mTISP(9nn7!73GN$+V`JgvP{g zu!rHfb*ve9(T_V0kPLawze9<#aK_8|ccbTWz+vnFyly@iw(r_p{4dd=*y;68mMy z#;5>rMv8z#_;Bx}fixT8G#B3mf$aI`NAwMOJ$buc(-J*D|HKl2@z&ENuQYPM?c$N; zf@ki6%^iaGGJ`$c<=WqKJ@dPPL@wMP6iU89qYw0!bh1J@97~HXA}i>*P79e>E)0_= z+JF_?Q3Fm;xb=AlZ4vU$zXK(j8{#+#LRp-o7uYN|WPg!+1qMskqOs6&N73rSTgn zws*$k*T*MpSK_Mky_@vph=98lQls&akcoPGXPMCCp0m9%_cA7*uHDPWAPqjQDIv#z z8GLj0eJo-XFhl>l-)f!*if0rk1wT_icP0q@*l92;nm-L(-Cdmqv;pEa@To2G+^OWz z-j24TPRL;IRIjQtf7x%3tc@LF3_Z+$TLIbOnQo@fw$L`>gGRc6A6k=RB#*Am&kJS+ z;>cXqyC>a+`#jobGpB**W-Y%@SA@#myFBBk`rc((it=fo8YSGKl*n1rqJqsiSbZ!k z4gWq5`F$C>l1tP^X2iB{9g&oKnW@@n`4y!@lwEyje>7S`AIIO~-I*Z&s-=yU%9_<4 zvjBYv-FN9t{O;~vk=jDGt%-@p59?2#P<;l-Dt-iC1@sC$v%UGu(yV{!B7pka6GKWG zlhGUj!oTG4!l|FfigXk5BfNLv2< zj<(d|fUpmQ4_oVq=GGLYARGF~-)(K<_sHs|Sx;B|;d*%=eRBs2{T!8x$(g1zicV3| zKj1dT(5R(3LXCFkSZC8?cRF7?%mH&&}sG$|{S*HP;qwMm_?EC`Sd)$`=SO zmeG$PQZ3naL*2JxbEuMsX*oeFJ4HBq7+9n#mr>RcJdwT|T*-IF_fN6j&q`Ipov?R> zH5D?Czc>kT+LKv1ucfqP`HyFtJd|UtpNWn9R|UiGezckUTZ87ap10=zrY2TT1|+i8y2bl&;<>WqYp`KsEp(&HzFi98VpCV z98Zx&*efTB$uZrsn#+A_!6XEvrwmuZyQm?dTPjul3|~z*41FaT$*kC^hBN#;2;3~` zpePkvUxoV3<%X}_ieoEc)EkVCy#Pt?s$@Q1KL7l9(>KUI2D|T&acnahZd(I+N_sz0 z%WC<}+4wi3xHxxD7jRX?36%D8^Hv8hMl}m$rrUmGtVmW{)5G|QmrD3ubUyWt$7t&b zRfp3FN~N=R=k@whLzDBdrHt;^xSO-14rM zLEzkup+0oMM0!i*4zp;=eifBFgif6Ak?;W9IXX|vALO+;Kh*h<=3F5Zzd#hVlEk|+ zj;rej=+!;Wuh(Rhk)VG+J6Gc$Rwg3EZe5i*kQiE3WZYA8i*B!Non-?%%rOtu@unD- zADmy1{jb&!4)uM0zVX4@C&1_SLgvA5p=%t!M>>>rfa`8pK1tO_tyR`APu;n^_*|+q z>@V{d$ivO4KK}__7SbaVm#bpi>3K&bs|~e1_*Cn0V^Neg{3N+apcO0HP9Nya&Sr4u zXNfg?onLbBI6u9qzoF)PI7{qAIRukohzS`=sX`z*z&Th^O*>eG6VPk8UXHjA0CDgA zTdPAk13u5kE`FOw^$}I>YOaREYZ(-t=c>|DznNCTC_jTS<`C`qZwrJUC`D+*Bnuz& zEliIMWdF8yoW^E4Qg`fCUdtn)%2ooG3EDg`3C06!w%S2w@||F6&wDCMVXp?Xoa1LS6CvbW^c8-XPi@k=TxN%~#^qS=KbDAbPahX_^0yR7-rTs*wapJ75+|fNv9&}o z%LNO%cyL;cJBA%EPM@tYFTV`oqyLKbUS}_V;gt0 zppRet5wm2#I$cuF>|;w#y3TShl1q>oqG+-yNC`J;;LCK)Sn8BPUL;FYIz%Sss99M0 zl~NxXQiF)X$5F73I`s?Ii0v$1iqwrX_(xrEUw>;3NfVFR#aDkT!PSUm7irAWaV~e% zWFK*p{^@%nUYfWEKd{QFi(|MDaV#j06U9K-{c;BbFZ=O?S8Fy5+r2ZCWgORMztLNF z0n%R>S3(TvTnVzqY44detV{Z?yY^P0LRiaeK;4iy{Y*5R?#759CA_*Kdx=`Tkcj3 zKMRYhfjx)0hdi~Uvrxx>Jll}r;O@p&@+lnO-`5ORb=FzMKZ+LG2OfM-+Ta_5(AAB4 zwt+C5!LP9m*LEakbaH+5ZyCZtF9{#;i_s%C&!^dU6`!LVH)hl@^wkrn1@nh0%c^+H z>X?2^*^{b2EJ&{k55_z^8*x?j+5bB1S{MtxwHp%g@^m&dCKU`TFbf_nkYx66SYYe7 zaLcS=xl!?;uYZmke_Bopr3@7q{ZMoq5VPw#Q6E_@R&BMzBXcq3Ix-73Y!=Ec`chGO z*o;0$#kOr>*%8nNv7LyJ3ZpivA!ty<=P?x3t8LJEOIo1@7myJX{~|nDX%5$vJ<*38 zT^Wm+MMq>xIUKPX8{oT|bb{u-Ksc9Tv3%Yh35x`JukxVkvg$*lVTa(bd-Blo%hz(} zhW2qPallQrnGal-9Z?5+Td&SdO{TRBw7$A*mqJ?}j`FFEvTiEC*!(^|t2!nhaJv6+ zs-7mLH%96$kRW@6$d*c$1#Sz6lNI$hq9_LE0$|19+*s<=ni9yDJ2R(DL-Y&1`)#9V zPvBbZZ#Ql?wZ|!*_PFTCBJcPQnU>iunT~vP`V{ZEJpW9N&xe&;749hLZnDYvn z;0jIoCIr`x7p+f`P5*7x^<}zJpG$Ee-_F@-j3)VXi%hpvS(9r?nFEO?O)s7mCKe$5 z?xRnM5>y_(J1H3}S!Cw1cKqD`1Hz=MZ)fp3xfjgEcfucavNV-Cq}vZA8AJ+O%opFG z$Dh)Hp#M%uy)^Xl3wbOrZ_!gRj6+|Xqrd4Q%_U)c%^2B7y{ybL1-^ZI$!QgzzV)8B zirJ7^F-qS@3KN;4PZyYX=woN|K_k~LLb+*`J)QAzv)WlE(ZaTv~0Bf#x|DC?GLuamEq(DoUUBE zn611vj#=-P0CE~Ryl=(K<8leN# z<)b7=l9)0ctA>{O&RtUIOp(+Nv5POKPF+@?n+apQ6S9KW*ixt7M$zHr(>KEQr^ z$o=>%*eCqZC?tl6uaej1HJqq4h>|}>1}4w7agy^vuy7JrG$}!Do$h>!;re-`S+R`w z{3Z7Z<-q<4#VJY)rD3f2=Lu1#{Ng7-#PC*>kYCx?e+PcYr_u8i!EeV{EO$SLpc&{& z*FJhW80nTv*qROb+#7f!_O_?I9H$Uv1B5<8;L@tfJKVm-#UsYt;dC9EsaP4SF$6nl z7z5Xpq>8$cj0s!qS*(=kTBI`Y*&5XOgPxqOAuW5kVRclK5A~k}$d}mv#+rfb46+PZ z_*0*{k9)JOSMGVKANQF{jz*|Gr7+cgQtOC(F$G4lUDB8s8#h{g%-?Wjw~y95DV6^bQ7X!(tAYJq&;}Dth|1I$ zJc_iK&RI^?%4d0x(tD;FIvwNt%S6AoMXM+_chc#cSCo&?e{ZEXgI;m(H&GS89{T9L zXCyC|P4!=V+=YTHJQT4V0YCl@vY_KbjK-LqpFeI#tCp`BOb&U4&|GXys5v`seuhc# zVmWa9Sndf>fz3X(-bVC>kHI+jqlArS--0e=_}F^86nc+i|FD1S@aD4Ns>DE)QL+c2V_IBJ%m0K@p5`e+;KDO8q7l7C z0kTEf3%GGPDVCL#0tgqC-TM^_RFxWe9OBpKJ{%a{Eu;VTX{qRiU7hj_BR&sdt_612&XSy3Af*}kC5zRK;A#GG(L6*EUNqFhC%>S_C>&04SA$DU z`D)w``Kw4~dtM}6XH0{0#7+Qnb9Z6W5a{BVQjupj2lIDEg;)e!!%lDwPUn_!9)$KMt0-WKz^vmo&|pdZs`1f^o1yY{dr4A;`Q+JskVYY}a& z&mNID{w=39V|R39TspvCrUqbs)%`h9OjB~!$w841o7HNLjLsTB@^*fV*%b35+I&5v zVq4)Tdw;*JZ`hKP06yzJ0l5y^J^-tv5?)SOuY(qxsaDe(zieEx7xqt{iI?hfd?Rwb zX}swD$QRa^qeKL@T~{dyy86>*Jlbd>i_RON2akDxJW5n>wOrzZj)dr#bCrFGI5CZA zNtMGsOI$q{o{)_OR2F^DzFNgqm71WQ4I-tkQ|V+hzQvT)Mo`pnTOj32Yl;$oniAKd z7u|w7h|V?H)?3_qeRWr&hi#4$P5Cx+dY!NM6B`Rmh`*SALb)pUdg^ib5GAv#DZ0uKCSV93Vg=FUoim?~B=E?$@jOlV*=% zW6qY7r$0vv${3`lRhJ`;!hxCbY7>9>2D2~c8%@-)1>AsJWKVAMdR8mzw)#6=P);Vk zAg*sj5>VZuf>LjxG%d29rc|RLoUFO=Xy@ShG|j3Rj_B2#AFKS#GYkXd+p%`3A>+&Sncghp!a6)LcLbtrCJC1;pM{q{?2EuY5NHBsL5N!Pd^%c9#rJ7%8+ zWMH>yO|n)>eJ)UI@>k_>hKZl)K_*`T`xNk3Li{N+#Z`whR9~!0gyRnmsQTIAouYbO z6$ac56G=(H&AJ&tG9WEL3jIEr{iUDgvyk!Yx!-1b3laYPzu-&L)My7Be>>!9i^s5$ zh{lOr%h7)MsIwD)uMb@q%>IU3B>E$s#rO>vh0|zK3}3C4Bh6c6{&Zjmf^o_D@~hU= z0Ke}@P+?*+{wvH9>4vDGUfk`kep$3O{^d0^%5g&kr8x8! z6}8gS=KZH{O5d(d!Gw5;N}srBAzb4I36rEm$uHWLmd87CNBweq*`8TrD~Qs_pxC|3 z?3+-0o=#iLD@zIteXktSdso@#joIBA7o5&s-hTbu4Te(J8f-z*I zRCbr?u%#B7L2A9i9y_&6arH5rw84v;OtyUt1(~P*FMJj~E|xd?6M5>|$C;HKYwd4S zN+9kO_z5AIR6>14OQ6+TlE62WYt4LhA353?2_Nn;k%RWWfmuPQuO0o|L#Dr<t z;){|fTV6Z1zY6oxo$AqpaS2F&9m7&gh6*ko9agBib@Hme^BsTs8~b5y531SK-=I?k zi;Q)ew>BT!7B${q$*2JjN_QW7Hx6Y7v~yl9RrWKe{@q$`;e$fY%q8Fu9mCFekau^) zDFx&sfKo~56eeH%JxDE5zU}lcofGeDMKo3X8yL;?fe-6ooT!Qjbf~qv`RDg$o7$dqJvgWdQjb# zVN<`TlrT~zTUD-KCn%(7WjCB_y!nirxbR77fx=L2!ra#JeW<5dsE6Jv>8)%rs{gSj zW|bvaXKudF4H;7{a-nVpTXkddD^F*;w z%qU9kCp9ELkFH}yKG*l{b-gL-&X(G_$_OR@DZa6(9jMnKT zFBSGJh@-G98WF7afJ4khx4xi3C@u}Z-e4P*p2}@p6S2gWJd3Z~Nfgnd-90iG)u7c6 zS*Wt~f()(PWo-exouI3!bz2TRQ;jd@mri_=tg}sIC!Xm2glX?K0Zg$*YDvgot%F?V ztQFJ!$m%wYtXR{1`0cL|c4?{)@-4Fvinpc%$ z$>?e>*6A1DCKIyEEsFSTe}iMuw%HP+FSZFr}TlP`i=LL zZ||3~u-n5d13qrKCvgfVGN4ZSE>jr#82GTHyHjL7M6~7% zTa0Z-AD-zuzpenEbd?doSqDWQ2d{0P?0&zD;{PjnpU$c|bP%&p$tk2u1$WDv2{b~> zLp4@=zMe!VSC7zp*tBO$oaM}6FAQA6910DB5WJa5tDzQ)7WX#an`$W(!R&aN8(6=F znctU%3%|zTBIgFqeRe{)1`V1~h4y}I&J^Tk`^9(jnEMAZu-VyA3%so*3}oDZfVw7i z^N#axUj;&#YrEj@W!m5VY*|y=noA!)Nv_>YGfk_4XradFNrM-T*+SLiQm#orubq~0 z7mKHSQIH20b=Hs--Ziq0b;N1BM+L}At|p>6l^0`Hw|ZoXu=_h_9l*DmdU$n2?rqUC*=a22+_DL7eoy7UL_BoWj7yimF8@^2cQt2LJw=g&4;=g2Oj zPKT>rZV_tY+;Ll;seo){rFwgOvX{efqpiWQAp#v2yg*}s$E_kdKQ<}ffd95Fzs@O~ zv9OJFf=rGQc66lrqQj1)WE-9OP-(iOs`)~>)@kr-ROR3VNzur(cq9AvywR28>a({} zWx7^__V+GUw|u4J5AIDc)``pkmS1p-k{UF<{4FjSS~lAg_og-m>Xp_T0|0hB7WK7aVu40z@KqC7o#;mB z__+t)A(Bvung9)`u-C!rH~^7zUu~P_)8!qGlqSDE?>IQ|O3iGRjsAaq#f#m?sp;kIn`PcK~{?3fqjO%pl;{-%h9w zH95G-6wmWh-cKZ|2`1GFy4=>KCecjyZ-+11od8(Xex@xyO6SBm zbuBDu)Aw(mz($pLg_3V;FQ2pI)Ht|^^Jry;Q_U170?p1=;zSjDv&vTQPfCIJnpz zlaC}_qr6g7-LvmrNhe(x>XWLVzA(^X-t>bw9-Ep^AW#JHdSUUCRSef7T)9t_czi+T zp-1wgwC-{GyoVy33J&qTh#APg+1X-!qWp!DzQ)acg89Ku%aFz+~LpuF(itlr?g60_4HfZ~je>=}N&@D4a z%dgl;6uMRtqs+V550x=%yY~|!VKy2gH*Z9bEy2W;5R zcZ9As_b&M`jvcu(omH-xu(LB_p4x*}#BM*<0}Qr%T`4UWd*%`S+F{fOPUee6VxHvs zOv64Y3QN2XD$J6bEv8miDSesT=08;Z5m68-L};dWimR*cYAu?)--sjeLh zE5l(w#wIe0xy})kRr!@5OM=rC4ue~HUDhSU&7iht>=RCMn=vEWusJ7KLA-0?&y%upq~_61-Tg8WyA zJWNXS5mrVzVrgs(!dZ?@@u1j~Q)f4b9)?mYM=g75nD{){rYHht?8Gqr7wO#8_!4Sf zCUZ%tZU0xlg3}$ZT!gX`j}tp|Zy^2q0|S|f({9EH<;9K(qpRfMPW^nfdkblOS}S%V zDki01Nokiy(zM~yg0Xg4vOq9n?eTf(w9($1@`Z|r@}F^`;yQi63_H?MjVG30!{0BS z8xH1Yf39WKKFKRV6fuP(#8cI8zA>q)jz23 z0z=?e(+m$X-7_5#r(ifP%jCHB8Al}*3s#gIxsN#s3lBYeY7-v{;blu!Y`QvjWb-T* zrt9y({{C?Ry1OzH5sqwxUxDbyGi+7Tjko%F*VZfi|GlhmJDHgzJ zXS3{Kyi3_V@<>%Tg4!ZAW#i=LEn6Kg&ziRyMk>dJKrOYF^jYIXeCTDbV({9@#K00o zDkcAtNw@{yu>cgrx!WGPaOW+vrK(3!y`Y0u$!&TCPLb*ERnV#8hz6y4W_54bF7$47 zwEG+FNIB^`)9pmWV*VwS)pyBYFep2AeMhO#M%&gp9#k~pY*1Iab4Kp;1=dakXJ{Ry zH|B?9&n5<~iDm_ANvdY{95S2{Q`KnKa*`TSN=gN`iEN437}_bOwx3N1jnx20?jth9 zR~rZ=w;HNt{tn%7cwtd{{|%e3NAj_~>2ei5=Q9+!bcOhmc8~hwT-I(Fdad!HDHR(& z#UgCS+8{Ro!gq=JU&;_a!Ch{h@se4h9vmXnrK1_o#1^=q!Oo7gcb9c_n!`sGA^lIV z7y~AAlaw{lanRdQ$;xp@Mm$>O(sR9-epLbjlA~v4BkK4Z_oF*G;^CqLkhvIUMCfb5EW(H7Mumv@(zZ zz!HyGB)}JMr8fm>DW>?^?O>~2EDCmX)a2kz)T{wXr|rL1KRPRVVAwH=_Ev#N`cN?? zRWCeb*g0Z(#{XjLqRmtBpkEUM&EFX4_zjPI(A}mfTN`Q6F+^~Wx@L2O;nHu=gep~r z%_ex(E_JWycU5;YyH2Xxm||2W*pch+1g;p&57b6=-J(bRAlVrlq|2Ebx9w46tRi5M-|YfIhJ? zSrk%k`LO@+unj!T;_qodrACalgZXj_R`NR`sKKplFCZt#p5v|UlVB`PrI8Zw>F=|K z^3b!;`=!j*xSh~P#hwi#&*C^A*AExo@m0f*VpZs5i1wq6j!y7srvuCg5lmimEfFJT zQRVa$6nBiKH)*)|u=m&i012XbiVTu(Rl?p;!_fpzD1!YAu_O#3Ir5>W9NiF-)~lZP9~akt_S%t6#}s7-UoFPJ$WxQTv%`c~zb z{dt|OJ-DqBR&$!A_qgIic{>X4Q5@noWjpc)U(TN5T5UhGA_Q^OC?b(9xj=AX>`_SwzW zXN(rFpO=+_7c#bQPd7aOtwE_oV-c;HZ>Yp!@3$cz&42uZ({BE4;?u@1wq$S=_k_d* z``?f#8Hq_T*WfZ3%aQ*(J>?;ACilJ__F;PltfVO0&u${};e&)yd$gd|}fvO_|co*#8u0 zjGPrHt6%Jy+!Qs0R?P+R_-wC+33L|&XmNQuSnR7atv=5`-5LxhohsFrSykTtPf@Q2 z`gy@Npc@pAlS4_XiZzj%toQ%H9EauHRO=uAQzq9N+irJ9jMic|k`+AF|Ai*3FJ5of;JIQPXcn>>2b6 z>%9y##_tXLJnZ5hb)fvjON#sT$`g2#2s&%q1{H{r}*iDwOU=U1HB7q5qyN z%zv^)>Bqz}K9B)|TT#t2iW8+!)B%YQVOy&Ccb5RPK0 zE8Qv4a7scQ6?D-z+uey;w0$Z^Y1E$ESRB%p zOa3#I(ze_ zlx}R7R;{1^muT7*icM$c{Ind zdFo9lG6JGADtiyNdur+m=X#l19LGpAK52W_xtH?@zOzl8B(%gs?LG8?ux}}0Rli5` zDV~RaF7H#g&`NImwHLAzIo{9#LOhbDQ9>fe;*U)zm1$aB+3yeHuSuY_2fr$e`grBdHh~+ zOIiU6j?ESFj?F`wQ9lrawp(dMr^cHgEB3RFwiBLSx4`=m4q#W zEBT50ukBYDx*U{j=e_zrRIkBT@w5af8qrt4u=CLIH?>)1?xV{gCd1Oz>lJ=m&|qq@ z@1229SqL@}vG?K~YQ24%e*9LxHaEKGud!M(aXE!fzJ64`sSdW_UUWM>bWKGxnnx|hgSWo|2)ru)e*ZBH1$P_v*ilW#7`T2bCn!~ob(ZvFtBju

xgVEz%QU3}!iREIya~|mK5?mg71HK^e zAzdB5R{i78dd4dr`!yG<>WiXT4MuaO!bb}PVfLFCK?aim?1AItKnn1&%3bJq^>;zX zhwE#$#Wt#PK-Mf~Ba`w{3|kQ94M2WAfN(K3x?%+- z%T#uEPg9vBKeV+jBu!TPfP;lr07T_`yphbPSwUY0f4IesA#3Z2%1!&!Q}H4_{cEmi z=~-S_`aNX8GCb#uiNQIRcnL^M6O|2hgIJ{geGhn2odsC`=EQ7 zFTK7@wFyjhw>(da9oQUi8p~?)7W2~Gl7mCmP;YdiiWg~7glj@u623^4Zyb>zc9%Jw zg!)%nMug3G4DC;L3>1Eik6y_%*X{U|FHt2myZdJCc0W`D{khH&L+RC%e;Vp@?~#sJ z6$a+yuR;v|pxIgtSL-c7Of*L*7ZEkf8`z%4+_nYR-t5=j7-tjpU$lNj(=PE1`gD#; zjS8$7v$meEtnE@%`>fbESf0Gqma;0yFSHMu&kYEy@Q3s;*|fYQ9JHG3t$*$qKjLN~ zTI9ieXBBMv?P2nI`03_ABk@T2aO_}o@7Wno30=`OvT&$N?(6lbYdMYc%GRePs<9j^ z1Bxar4z=nVc{SM##_j_|5cdvU;x+=?^-d&;<61Sn9&mPieBJ+c6rV-1*w~E?r!x$Srg&QJIJ|P}cy!xCkEz5F&(L|;clxUzk#^qt`lVNGg8?dJGO z0e812ant_`ObN60*QIn;o>~cR|bl`sXzAga1N{s*NNZJQHM1flsPF^U{%`m=;=zq zPq$y@&D}+b7Da$oME>tL}Rd3%Y65QjxYci|%fLGJcHuiZA___x(r>*-8JZWoBsKKo! z`zbms_{t`J+aQ?!Fz3QvCGg%mMA_cqftRPOz3;=t|9-pGQ_~jv#E#3>VHvN}&w*RD zwS8kU)_)d>^*_bxb-Qf-p`cpUkB(v|1nawz+|4=|>W z*4gHm4fuPp59A3z_O}2-e)k;AGZL2|;uPbyoXgG^?tYoOKn%>ByojA9zF?yJO0Hei z0wy)-@f|Y-FvGvU?*#iBWdAw47FLqAgWZmWjZP<&VL#zR-saJ494sTU?d{QO18%&6 ztuZ;@W}0&Qfj?=TNA`2I8S(x1SYdtjL@#dT6IcuPd%yp&o*j2qwQp<GQ_z>}N6pV5`<(Ni-r*kHz6AWjZt;uEjKE(?o)>?991>SYg1^aq zxkiZ+_$w2&lSY_>efSF#i^@rENDmK6l01?ph1`k7(aG753vLCTQPs{ptRJ z*Hz2Dy!S)!`cBHc|1A?d!TiSU{in0@Cu`lYEq8%6WFK3sOc`B#X*y6#=JA-etU*2X zD7c;9ww!Od56n+#-^dmi(l6`4vs!ju(4!zXl%2nNEMMfJZw>poruBv~A+VqQRYKY3 z6ByNA`E`$%!S18j(j7S}unT;*Xd|Nr+qK47onQ2Y?Vy0gr|(I^dNtYC&+4$_n6-L7 zuv$v?d$tVd(sI3^8Wt+w_vM6z3*_ayb1z_~Pe}ESDljihtlsyQ$R?o}NOl)B{Zr;Aho~K}& z+AqQR;)6a`g6%xy0scG{5*{skJ2A89dGBiN1YX!ObN8iR;h0t4|f| zcp>u%#;9x7g$LVVzu|SIpqTZ;-{M+={_uePoXMew4@tm&MqA3W$?su5XoS=~MPJr$ zlXW(GlU?itvI2H5uB}oT$c~qjLjw*T5U}%915I}1{5=1VykqR`G!yW~FSiX6jbrnI z>va-hbHJNP)*pBq$^NZ8U$QSSx9w@`qi;jOE!~_nciR9q{^#U2Zyp5R@=5hGpSFSf z(KB=P*(q$EW*i;a@)+ESsY`8-E(7=AD4xRycAl}8&QD zeWP{QeQ4!}1=Pa!Q-d*oHDf~c`eQTM`Lt5| zMz}xSez{`#!AkG|_K|&YVMX#FZk~?V(N(O!;Gc`x64wLVQTj`}tQi3;93<%RbSVIm zS7KP7g!&6sfa;ZRCA02=D{?O~B%RfxQ=e>`r`v!uvTrokRlfZMAkjPM!vifgpJ1oM z@7KWI`eICvAlC2h3eL7TUj{ql2j>^nyoSwlHZ6{5gpFT$*8TEU*l4>OltuM`jl#KX z?M=b3DY&HZF2N1fE$!Rq>cTpcJQwSw87;1__rl8YmQ|HXFszn)AL2)?ftAbA1v_6V z{lE0U;fRZ}FGM)v92u4};nU|0m=-)WbQUXgX8(ysvH!KTwJ85vMTsAkU-0A5g3#b_`dTb(+>9T!CQch49z6TUNfNWega2$nX4cpV)2C0F zI4z78NKc$LH+0r)bL0OiTl|}B`G3gc=FDpE!+-6&_^cVBW7%#fwfBX+MWyezKxEFp z-{;sV!LxqcrzkCyA3J@@l+dxWr-uiQn?5IG%AY;{vz!%L|7(j9ZU5GY(!UxRJ0Unc zcCrgLfNtWf4}YT=lmG;BK`hN z4`Tn)!;~K_NR0~}A3SHu>>qbHE_BMossGlA9}P;g=LHhC^Fzn~IyNJGI{P@LwtE>B z+L_jWYe$^cm1}P!aCXR)z;WYgk>r1o#8`>GlpPm5J9q&55T=IyX!ZXCTe;K5Z5#mr D9)$^t