From 9b004db6f390ad3336856b571557c1904d5b20d1 Mon Sep 17 00:00:00 2001 From: Kyle Husmann Date: Mon, 16 Mar 2020 13:38:18 -0700 Subject: [PATCH] add sex / age info to nlsy97; restrict tests to 1983 cohort like estabrook 2015 --- R/data.R | 3 +++ data-raw/nlsy97depression.R | 15 ++++++++--- data/nlsy97depression.rda | Bin 57988 -> 74098 bytes man/nlsy97depression.Rd | 3 +++ tests/testthat/test-refimpl.R | 49 +++++++++++++++++++--------------- vignettes/mmod_tutorial.Rmd | 2 +- 6 files changed, 45 insertions(+), 27 deletions(-) diff --git a/R/data.R b/R/data.R index ca12097..0913db7 100644 --- a/R/data.R +++ b/R/data.R @@ -8,6 +8,9 @@ #' @format A data frame with 26952 rows and 7 variables: #' \describe{ #' \item{pid}{Unique ID of participant} +#' \item{sex}{Sex of participant} +#' \item{birth_m}{Birth Month} +#' \item{birth_y}{Birth Year} #' \item{occasion}{Measurement occasion} #' \item{nervous}{How often participant felt 'like a nervous person'} #' \item{calm}{How often participant felt 'calm and peaceful'} diff --git a/data-raw/nlsy97depression.R b/data-raw/nlsy97depression.R index 1323137..39788fb 100644 --- a/data-raw/nlsy97depression.R +++ b/data-raw/nlsy97depression.R @@ -5,9 +5,9 @@ var_list <- list( pid = 'R0000100', - # sex = 'R0536300', - # birth_m = 'R0536401', - # birth_y = 'R0536402', + sex = 'R0536300', + birth_m = 'R0536401', + birth_y = 'R0536402', # sample_type = 'R1235800', # eth = 'R1482600', nervous0 = 'R4893600', @@ -30,12 +30,19 @@ var_list <- list( nlsy97depression <- read.csv('data-raw/nlsy97_subset.csv') nlsy97depression <- nlsy97depression[unlist(var_list)] colnames(nlsy97depression) <- names(var_list) + +# Add NAs nlsy97depression[nlsy97depression < 0] = NA + +# Wide -> Long nlsy97depression <- reshape(nlsy97depression, idvar='pid', timevar='occasion', - varying=names(var_list[-1]), + varying=5:19, direction='long', sep='') + +# Fix df attributes nlsy97depression['pid'] <- factor(nlsy97depression$pid) +nlsy97depression['sex'] <- factor(nlsy97depression$sex, labels=c('M', 'F')) attr(nlsy97depression, "reshapeLong") <- NULL rownames(nlsy97depression) <- NULL diff --git a/data/nlsy97depression.rda b/data/nlsy97depression.rda index 0d9061e88f66988ae53386ef5a0f0172e25eeb38..5258c0edd7dfd2f4f07009458799fe7bd597a025 100644 GIT binary patch literal 74098 zcmagF2|SzW{x@to+Ce+5wU)Mq2(>Q>S|oEy6KgDqgvNGCiCrrM!OWQ|wX2;38A1?C zw01$yjIBuQTPNK>{7KPw`lhCeGJ zA#zURCpp~Nuk3T}3q?c}OGO%TB}CCKvLYhtI0^O{NRsU>ty?Yx^V7eZz*GNzRsQGQZPw0C z9a>W$^67tF68Y_aJO9t=|2O-8-}>t3-y2v~bxi+J5hLZl{A(Ad%?$jX!rEHPTqUKT z=^zosbJ;@hOn+7B|5c^2$hn`szO8n>CaiK)N#xsek8e#M|Lb211cGG|0Xp`t_pcfX z76{G@UtLhU-l3uKuPK7IqajV{uH(1HzZwr1DOjITz=+WAicbG(>~~kjTBKAwdCpNm zq%`|$jc@kSt7R}>DRk{6{pGL!WD{H$FrvSR-u-%Bq*Q|rz@Q=ai-_ni){f`I#%j;p z6_F7Q`b9(JjG>5zi0D|gqXOoZG{YW`noTzr5jpv91^>GRJ+UJ|hg>nDcQx)WSQ>q; zau*{aIw$Jye?jKf97Y7-O~bm>(ReJ)|D4$9^e?5LDW3d`h{&xc3uj^e4Y+ENAhqxI z*#E;e6i?5_U>=D`8VCJq{Hv0o##bWKh9ZuJm`B!X8XBSlcO4atjHmBnva$_Dva*7I zXP>S#yEXsg_s+(U=XGui;qIS1!rhm?|Nghn9z`@Fn^qh`44@T}h;$m8NU9AWkw`>g zqb#6|#0sF1YCUK)QUEQWIG~tVOsXXjX>>XbTnZtvgpaa2!F_OnL}b(H#Y7q#+yH%w ziJ&vFn8*gh1r&z_kZM`Q9-up&7C@tcIZ5H5Z!u`nXha&Zh!y}Y)4?g|$pXhR5Tpo< z2|9thbZ|x@xe)`vH6T|}F$qipreT2rNg<(_RT#gO2jjtY`TY3_2i-OkIBcRVRF~fY zH-Vd&0y$|yMS(6tSx7)H5|Ieb!oiF{Y#>pwM-do`7C@&JyMZA+Xod4U-NhmrkQWFH z)FQcsgGs>%Zf@W%5MQVXkiXLn$l(?aG;+{^!*wSdjD#8i;laql{=gh8AUL>r9X1{Bb3<@0%P7|%p^oZkix z(9(UaK{v_P`7MI;51Z!}ymo{6MNudwyFOe}2A|zY5po>COX< znt(-g`4f)zZMto)aJY$X(<0Cqn6Wm1)(IW~cUyr-{0U4eAIv(>=Xbz>U16Xz42TPy z1dLb)O!$WZYr%s;&yk3oz!l(rQ855gG58f)naD09`SEo0Fkl#rri~xL1FAfnRsSP&a00T-l@qm4Si(4n!fHrMm ztrPY@PB1PQ2J~q&<+;N66ZR9B`3Xnx0O(>Oe5?%y(4Y(E)a9=NS%3w>#4Ds_p=X5F zB9S6)#IOLaJpfjPmIoREAOO?};ZRIk2B@c(0ngCrbGkrsLChW)8;lMQg##0YwZUOP zFIV09Hlb=fFz13je-S_nV3)Mo8v`CJ3v}SVW6H06e758gH_R&d>!l@RR%yfaNfa;2 zBk7d*r*9hn_TyFUQy1d9pHnvlyRJLW&2~5w*B($E0EZm_RILpMo(5tI7aS45;?!+lmj~l_ zPB`d*bqz2Nj6hEekWKuiHrRYC563eCek2A6S563c5ZEIG z%+m?}fEUqwwGFHm;Gx!5;Ttfh36KVuR0oU#yfuO0fdk+IKnBO>fgVDf%=0I1zj~AS z!%^?X{fv1bG(hhr;MK4WI2^_YNc;C-O<=UPu=z$4fEgjq01Cj<6PR(Ze3~YL0cZh0 z|G=zpxfZ*D1yk!Llu^iOf2>*|@e!Zja3|H)hJ*lDvVqTFel(+ngiSw*clwIVe|YBy#SycO^71^6c91O zi;4pR0Iqxv4+Wh2#gnHIPB*|BqZJ8R9q<&e62JtB*sKGj0ipn*fKP$PfCNGjcqW&@ ztT3U%wZK2Z@Wc{YwXt{)U!GDAkpbzL9P}VH~ZtMi(g3*EOz;XZ$q~0#zEAV^~jqOHcy9MaLT7@#L!VNMa zKEKIR7d(Kn6->?LR{gJzwfXE!Rgc6odZWU2r4wGTc$ za&mF`zj3c$*?)cQG+9vci#?tqRJuqVuT#lB&uDp zE_`|Snmp<3yA18oIaJ}e&XJF|UhJiPjQHy}Ps)O`|7w($Js9e_pqGvO@j~25eQ@@T^KmCX zqA=BM263EcooFd-S$7y8VXjI*VQB9@WbB0m6)Gol*}%kdGFTG={7ETQr-6QI>#QzIj~@LuNQA5l;LuEvQAPoeLd05R+>G~V6q+O!EusQd4K#H z19d&VhpLREw0Tbmwsv<$4L|1DQ4C?G(ftBlLJERL*&Qj%ot;HH?)Vz`$&RwkYLR>a zbR34IRL$^KFG!qEb}vw6{P|JW4+(2=D7f~|XQ7{(irc<<;X}5mst1v|nL}Chn9+z| zj66Ni$Au#(E8S%yGgBg8OB@efW9si)(cf@7~Yf`&Ha)R@b^0MtO0~H8WBxQm^0DS$2c- ztk(GCiBeqMWMBNO=VM0a#;|FvafvFcebp}UV;wtg^}f6~?iEb7_GhE6vG_0Tz280$ zJ`QJh8g)&o>-VQluxmOlc{nur&5s})IOelN@s!RD*%LiPe5aM$$Sb|BhG4_bvU947 zYpLm82>ZglpG}9GGd(=Qm+-4D!<{GT;dY7jlnn=pfxTkypYOw+ ztMmpJ)IVt*x0DATU;E)(^EY0Jvyh*W@wE&s-QeTJHOhU4`CQKTk)5S;Prkdk=k_Hm zd|zQ!?Y>g*-s_t$vXS>9`hnUlr6)@Yt2xH|&R0CIrkIYvjT=%^JO zKX2NB_~8<+EB8e+OP~AY$+gFffvOR~z2VLchxm_9fB#(+6V<==!-LAdddRFk*d5LA zvUBNcchz39hu@ z>3|tpCmdUOa9t2I;HH2wiGY~_qqD(S;Bg}0=3tN;1hj)76RPLv2uNBOT-d^-7>o`M zkbi@mQJCC;phY50_E+2&=AR&)0J*O)7bAsY_&gAUgn5-vsa0JN;X$$nW&(^1GzusK z8WfX49dsu^g78OA*n}g{9Ay3=ONN1f04xInT5FptU@siFE{spnyNH!0aZh zjX*gB*hCl;LEoJM#y1aludcnhcI5{_LV{Rr_yY9vB3>rqwbPfKt{-||tJHtqG35iz zOu-aDI}m(ZfndV)7Bk@>OjJ8`=Y^sR6W2i zt-}0hBAD0-q!E^1oxq_&A5Vb(Ai#mA!2|37aF-Uq29=jE`~km%RvAdoK-SXUhMS0u z0OPd5alm|Ge2_D>>FBNiIDrfSh~*06FkA=zNB(GHq6@4eY|2l{r!RzOmAsz+5J5G& zNZs5KCYtkYuCTT!cQCKeMQ*~GK=#5aDh7EVD6WJeEdy^0lO@k0E|?AuQnfbSuG1eV z{|>_>P-TD^4>CS&o*o9|+`ufrnXUYZ{nFm8E)E~0tYx4RNO#z58gPWL++b%-I64Y} z56sx11J@aUC#)^t01rY82x}+4u(FvS?}~W`&YQrP^ZaTYe;(-i#{tq%^1+ylib3Jw zc*wR~#@!U(XUeUAjd#y)zkWK&pNlyV}p7+N}T9WQRziH7L)2 zNXw~`jlTUGGBW$Y#OrIH-v7Fv{hk8*7(!X-kNapJf3olF!+ycJt!yDq3B|31?|Fcw z&vJCswzQL?{aE-p9h#hADWmUy%!WJekbBqmd)|R03I@Cu1}eoi-6pUmK$G_eE*FJt z0E&FD?0~nHLkIxc;Fk!<12~mV2MbYH8wpoICvYsN`-DYRNXUJ?#br=Cg@l8|3lI>S z)OL3=9~DPzycGkvqJs=iLnkG4wfnm)X}=YOs4$F zz}Sm`6F{EIDh?op1LgxL1pFzq6v#M9Bu7k}E1WD4X8Pb7M`Irjv<9m2|5H~BvuZ&M zr%eYeEnQH!;6Ppu%4*Ozpj-m^{@#$^Stz}+3ajkTkRy9!OHRf=$oQjn|*wtyuf%!R=^r-Pjg zAPb5nv_+m@Z@o$PwZS6D%H>AFlrN$(AJL&LHYeK<%|>&g1$kv%RXP41!9K zWJrwfSy>JF1(%e1QFT((iU-y= zN|(xpqpEwKQmGWY#Ed3)cS~w(b4>2=U{dJ}qC2R-z;y&%-usM> z4j{U%BD^?<{ELV7IAd9#!@LWm_?U=Ta|}|=8 zJve%z6L!9kv~@kWUdm;1)*Zw;$TZt}|dU!82fDk(qpwiR=?P(OTg zVQw<%(|i5WXFbbmGXBy#TuV9l*K@F0cd2|TS{p4nMmDypNY&WX@k5Q;&3^iLKW6vb z(LwL4PZ^ldi<|PbOPB!({Y9VKD8%X=eVN_Y$iv@HfBm3&)Hv_~()ilwfTS#|%x}n0 zBDz8w2egY<4gFQmqUV(+JzPUeQs%XD4}!1}Qm{b(qO@+jJkr2CZ%a02MDA;pv+cyf z#k!0NNhThV)=>L&D!fiUL=GiomhkEmOy%OC=3M_$x>;A*00}ylA}}{k7A#pbR@EMM z2mMq(pyfRe`H(Gdt~!O2#LkRnx?7R$-7c{@(u(2jKf?}JeCQ+ z+H9{?8uT!sBd~-;MQ;(-i8Zx2CdcP8N6zbDxR@bzHu|1zdZ|2vx|ckobv*e>P+Kvu zL!C0;;8-gaAXC|S$ej4;_E%R5F{^uS`HBvN+rqzKfQRJ(e4 z@1DyiJJi9cgeWa36;!e|tw(NUjR{HGQ&M5N!G48iPK+gRvjtqc$+MwB;jPBb1f_BF za$Y(|Hri@?$I@&x5I=~FUF|%g>NNBysm8~pOkMVccc4bmp4L+CSvz7)BSmO!zsWeS zHxa%TJB*ch5=zIr|ayQ|gc`9s3X z2UzGrR83wM@3rEMa{p1X8uDnXIB@tr(*YgX*k3;I~o ziug#kAfdOnzRISc%bUEz9r9C?f;B^yY}@S-X}L#S=&Sc#pw4{%-FGIRqFQfN40ewu zox>HFQzjEQlJyf~I+vm{V`v!x{bOhR)AJ5T-lPl2J=zT3uE)eT`V$$x*;HW&voen%xAXcYuhYljHQx+XiI#U`o2`Tc8 zX>3VL<@pRIBJIOtzM>=^T8Ovy9TR7VBC32>)sUfu>dNIO9Gi~pA@l1uAc6G84lS!+ zX>nAa+v<~*m?^jz6VKgaIoY_((0e5~X`)fk0cQ-Z_R;xt`=02)(9Z2P^56v5Z04&} z^IqVt z%iQU#tOE@W`rVosT)&Pw=6nNUYo1MG#$Q&>L+>Wjpz7<_V*J-K8xgvbtB3E;;gvry zFdCDlX#}5m)+Rww{c=;ssOGe$SBHVygX7L`q_x8b&?}kQNr0igJ zuQtDLzcvjmgC1u#G{JdVUHNe_!8O}&pb6b`r-syyM{1rLrOA8Bv7TT0?KMOB(5x#c^E#u@4C_Pz20eJc};zz^)0WVt&s zeU=3ly7UVc8DeecP}l^Wx#FL&R4MG8q`i1K+tsIy4cS*(nI}z|^Nv=fCB>vDH;Pz( z;H3vleQ0q+DsdOqqf>>dOqNl1`!j%CQ!4r7m#T4le2p?P)NUgOGZnYFs@ z-LvJa#>DjMb2=C_ySG@&ZT;Ow4XT?{odWH((`o;(@{B4XjiySJ^;scc6ZGa-Yo`Vf zrNZ=*vRNy}X41Ttl&|29t{{an@VVO0th%RW~dSg(;3*B)WX z|EJ*O;C-Us<@d33&*R_s_8EPnYOHN?m8@h@gf+q7GtodeG>pSf~`b2+m(S45A`0e=={b0b!Yea*Pr z4n1w9hHAJ(A}K!XSgnYz9xQVyAos&QuecLTS0U#w$c#3KyV_^qaBXU%t@bf>a<>ip z8m4P1ItC1M?3FUFX1VJ#9>5jj-6ffr!1Tqoj;f*o<`Y$zffn?5v<)?L_BGh{`uzTh zL86P;KUxVq$JD3B_aY&+3HO4QygqQEa`_O1V#g4qA1*;_@?98x8`YgRFgu%1IL~uw z{QUl}&q2t3F>{IbyIII>otTJOsI`mPuln4BX(@TLBA8s%v-*_CuPaz~!DQQAR$2_fHvPK7NZNZzo-XN{C!-_zfuM4<6G2FJDF}v>lUwdxlc;>+n;#MBwK8k3 zGddg%BM(&GhJ4fmz<4e8yVU!1x58BV%~8JHL$h8+Z~Wn#On&QZy|J9ASc4|Y5`Dql z?14%d_wC)}f?7Em_V6|G;1KUiW60rsPA-W-RQh{&Un2G{X+y+zHhM2U4^gOHi)-)P zGTzo?*$NqU)`oxs!(_==*Y91!bB8Fh$Fgy0O{bnw9>Jr{Iap1N2i;3@pv;4I_ zsb|{H?mL>VhQ+PgJEWr$yA)W>8HsP!C_0`ma*&i!*I~UbCD)dqJX~*N#i1$g#IVaZ-1uanxu*F=R(QDSiGlC&3tn)C@tosx*K3|zeST^$ zJ)2F7jZgNSj8}IuT*JRlT=5MK4`*k$zbU6mM0Vxed~r20EB;Hi=}FhM#FI{wZwH?i zvECu$PhQ^J`u^UqwvvPL_m5X1F27KXU^z>-vyNUPlAnj~OBa{N&BoS zywp6H)4l4Pzmb74U%{PZcYfLX85vhQJGV30T<2T+ish`_zCUdICFiys{^rSjrMQ!g zk+sAxBSv)-=Ccg*F6qU8!H-`iTd4)7sfhhURTP7gl0>1^CLQsEH!l77wHSR&WA_fW zv&T=8ZFt5DgF$+|>D@ z!9j1cNgtP}y!c@rBToxe-m4!mz5|a+<@ymC^HY*}S0bAW?heh)zp!+ok z(b^fbo_BC5Ep>0bi&kCE`Uz&D*IYIC@|E4&advKw&2^@6quMzE>yaGSpq26Qonkp9 zuaOwWl9ignKnwD`3{xAArA*gywXslnmW(BQR-C5nZrJZ_K$1;3|E#`8O4-J68-$*< zzKo=Tiu(MmKn_70yZ2QQJu+I_8IRyHe3ou_ZRVrW{nbsPV9uq1$`Qsryzj0jR=gtp zS^jJ^^G{G``PYjFumAGnr3B8W)5fhlDg0QY7$mfgjZUS@dpdecH1gcFO8MVm)WUBcnv2uPHqv7~{F*R>u$D$UU)(x$32%{5)|#96{-1u-aj|9&vT`0dd~V0-!Pxo0}6 z657F5*6H?#MK1Y1mcC~fEzT5LWs~#Ikv)PsW@{s;vq-V5{0zjH(!Mu2-{2!pUwI=57E5`LmS%Z@Zwn3bmTbB^>JZG0z?lTy}eQA)Hx9T^ikfzzw$D*4#_oyS}qx(@+YB^y=I%kc^_QWd3o)x{c? zzsXPN3Uf7jAGP)?Vk?vGOTb!qlkFgKv-Le$D6ZDy?jjdzk?j0SSMAykYjt?%ka-V9 zlAdSMXq%RCv#!w6$K&CK<5TA|M`=t0?LZBu)q4url=nVj4uUZe8jtBq9Zh@{Sz zmh=hw%*Hd?%u@Xf8yX}(wb?-<{MTe8pHL_}_bCIqr>d}^R9ID^1I;4X+YVLhS!~lb z)Q;eF)Rq+m<;ie1@7Tns@e(FpHrB9w@F5QORB74|w_xF9oE$=C%GK%ws#&vXs?*S0 zS(THlnM@^g>OuJOT5MT3RtzydWn0s7g^+36bDa^jq?&&ldh|F1$4M=8&*w_iy3@s_ zOHg9Rtji6ZlW3QCJ6Wcvqf@%n)>e@$%`Kws5`_Q0E}{SBhlVA8uf|Kd4Q4ubUvga% zyj2<`kmG@zk>1v}v8GzBgW?NAZS&JTd2%LQu3aH(BfV-VF9)0)Y2HnK?h$K!Y21`j zS9gu>h`eMRQGa`+Sz{*;O-UbplU5MK-oZl$3|wv1T>7+N?Rt$*%dVihBp}n3eLW2K z?8%`D223;D`TT;Ywr+i2_W5WX$gvCUu*AxZtW#UH&W7eb4(efFPuUgU6c9qJS*m#w zQ4hVAJV}AWrn8GH4DQL2+p2TXg;d;Ld0bbvc|C>I7hbozIZUvs#eTnlA9hgf?RhNE zeIE-=uat7JPT-X42!4f_&57+<61K=F+FF)XegUGd1A^~ZG6-xaGHqp+gjY{VKz? z$4SX!`R`lYr)+Yjr_+*>{_*XUcC$UEA7{3<%vRFCTK^scOCu53RM}dZ{Ai$#S$XqaP2#eTn&cD-GSr>N6#o0)gR#220b9ne*6I=c_pF?b_*~6%;tJ^a13A5gZ zYr2`JwCe~H=kwz2fkm9=toA&KqYP$)b{c*4!&Cups_Cr!?tf#rL4YWZVeLkqrIVPTtA3MpJqebD~yaWS~b$qi&82d~K3L zqXpe>ZO%X=z-MMj9no8Fa50(T<(Sa%D&~R|RZ4D3@@P?cH83pnqWG|8Lth`{NXG7K z$zdCPjI?Bn$+RywXA!2kS4LB<4hlJNkZNW)%1bm_p81?^+F2S}sNsssK=jkA-h|ks znAKAiwKJytAuW#AO<~xq@@y(!KGmXE-nX%DZOX?>66N<(r9*|z3rfXO*{^uie6v^( z!#x*83pA6$lc1N`#0?FsPqe00q4iyxRJ^qMLf$HVDfWDC`yI_^&(4Y-_HxA!XAnW-32-c6TM_Ex}fyQ2*Bx)iDcw8EZA^cQ<5hh(-52I@Z1 z;V)@_=v#&}bIPsBxX~oNPPCiO7&Ja)OmDwZhZgYV^%(~P7c{k8=*cvbx#fch?**|2gl7bnv}xCa zwobnd(|s2RX&gmpY+oZ>GWSjkO;148?$+C@RqL{76Uw|*%>+?( z@cfkRBS=0x3XRrS#^^=C@}$)U_T(t%r^k3=Nf>S1T6SGzuM8s<3hP&~2`bQ0NwXt2 zuJzTFvm}Rxu6UP!mVuaPum*$Rm8b~c%R0q10TZD3sVCUCJxw!dSMW94lciR6Xw_pT ztEh7GOeUL?gT1Vp;R_3dt8}*Ev<;LVS~@n>w^th}!h&#H>r@ZZqCxwHh2_y6+InS- zjlb7Q@Vj-xexwwo3{hEy(s9&`(zzDe`OwN+X0X)zhEr>3xN1z*ICK<|?{u6-?cROa zcN+T?zSRH3DEkdMCy*lPd=qgz<-L5t)|woCo;Q8;OTN1FIkO7>2#))fJ7pp2o(xg5 z=OK_Z`Nv%4DVK@CIy2*bUiY%GdRp!p*Y?!&kwsnjAhpZZN|DWUnv|1o)7F8juypO6 zrE5-%$IBUJ zsTWTp4_|Wbn(i2#7w~3pKy(C1%!B}`sy=|UO7Wge7OY_dZ>CayfhBRp;pM*UZd#tQ zi@nmO{SM)9_Y8H&Y^c{vH1 z$f$&-P3n6jN89^$RjfdS+UTzxsA~ zrhFy(Ab?T)lWmH|)pvR9+ZV}hFg|ntqB|>&F{ZK&6~ie~(d%iUG_@M76lKjT#nv@J zFUJZKZiw~u9i3w;5_aEIGt+4_rT3KWV=sHN2=(zO4tlg(JSrDk5yM;`f26W1{f54! z9nj#xpvq=P(o0fZA^G+th_tj8q{(j0qaUrI?;-PGrtdu^HfBroMhrK~ zqT_<(PnC21uWpSVzB-l4!q&wii7+MF*Y5r}I!@sn`=MN}&*EqWvcSb7H3d3#1(&l} z$x!remT^gUyJrdcG^KrM_YQsj9@vejT2X0- zG>2+IN?Ch!ACyL5x|C5keDq^QsW^PD(3Uv)PE}lLmD3wpqBnw1y>T0IFPZ3VwO&pe zb@|0g8_`$UkDXiA)2>dviC&LeSk~wqOs%8CHIE=j0i#E=m?i3t%}=;)6;p|pBG{#; zU|pdhlbGCgy<7E1y#XCe6U9P;oD%^{uunNF$;g8yRgL)q<3cegt~Wk=S=s`N7!jm)OR3J{P!h_tsL!jIQmQ&3?aB5_d{ ziCb8En-&3I-cW_@b1c#F58ak5LW%HMfycLb!Bvx6ySu0Gz2{y2zh##2A@y3NoYW2q zg{O|5B6EvgAWYt#otU$$T%wS@PPy`B^`tS0ut&X-w=t*Q39n9jYol{$x9&=R`tV3% zOrk+&5bp zRnp7=V!BI_EY?Fw3g#4LX2v9@jP`6ZGb1Cdv^0ZW&n~Oa=>6*r?SsRE+_Q(>J&^(8 z`<0F1BQpgp(Ek;eQk4)hU(- z1PS}@;0#>_EsIMlCGO$!E(mBo`IrO&`fpL?!O=4BOPa9y-FtHYV>!8fUxjiXe+pX=+Hzo+Yt&#l^U$vBAgF-CbDl7l_+Ue=)vz?MH zC7Rqf&v7=f=CLu`sb-QIwy3^cvTlAqd` z&3#<0bN|7-_fHpxtRio?VfyO#{Ixf&pj?aK;S185c#rFama1Fq?QRN?oqQ^hFR6T3I6e8cD$&ca7)Ji{gIqY(drOp`VNvNMEYagAW zs zB^YX|lph-&AswbWhWh8^NG9Z>;SpbXZM~W3-02Z2HEiHR#v^aB#_%40c+d8gLh3>l zT~ncvRhobctSc$I*Q$LzsGIG^Okhz-yCIOk?vTf|K5jkn1gsR?Tvgcx<-$H2hZ?(N z5E_<$r`A7@eHMjYjpFVxANMe;(CSaR68a0NNc#zIvpljD^P@}0e2EF-NJ8_2ZkE;2 z$8Vtq-^wFRzC8mIt%uay2ywqe?2dMZFhc8UH>)WEeQ#AGUR~YIp*M@qO0Acp)D7*# zp19lBe=4(_B=z{a>=s(abJv}U^yN&S9vzvjpg65oKBT~6?$Vm>_-r$uMw_rSLRot1CD`$-@uP~CoSAyr+@I4sT1 z{VByHYAm1WL|xx9?b6m{o`hDLEeoPcmu}e7CT&S7Rr`KTQvUhY-3L+PGJ^fJqxjET;|<*8 z^3Ajcwk7ShYI%lk!RiNO>XJeIwQ8NwhF8%L>&FQlOWggK29IncTBm*Z{LQ{6Sim~BW`2hIEi z%s~B$5ls$FBs^T-{VAk0wDPPnbvIwPpf?|@m0PW0hLV#|zxwIJ)28fLn2&doEZBWg zzc&3h$(*u+p2MYdpZ%xjtZ%7`7G}K`+%J5s_5Azax?Tu+6WgD$%)6d2Znxm~WiBZj z9eKok%+@m;VHlT8=wEeJO7)^Rmj~}FvqtUR=C*FGF8TDODy%MDySH7sD!bvV@Z~u7 z%l=Y}NBrD<**d%kZ<|eq8B1GSc>mSFoP>{?45j-f=J7SiKVA zPW_1;H{abX{yzB0_*XsV>5l*>TOzec?f9A;mUdME88*e~pMcq03y zQ|pBnrM=HpVdYMp(!B`eJ!c%vDE@@9Z{l{SAt)W#?;>jxmj$_0F0g?hot-RH|o2=}hG=r|kO8Ma9 z7^cbCtu6K89{=)Ub_yW`X4k=p^rrde1+{4}VC_4s$#0SlIG+B7b@plX^XH*js+a6P z%IX&AL|bVYB2Jib(T)ycH8|*p`A=-NQe=wUv^b+}V1>J4^@~^QYRqcjgQFHv!P+rz z;+(mQ3vHFNJVKU+)PPO0lD5=LyzTkPRqdla6FMwx{VXpsdc?|W_J&s`N`sEsoV334 zvSP(BC$2E2%`;|N0>S+@a`jO!QzBDbEMzBd{&?VcE~o(dp#>4ZQK?vcP%U`~wjnc% zCLp8U*Ji#U+w>2GS!GtHps|w)4f5p!s*Ve4^Hz4+Oo>*?O1nACAVjQ+z=jqPq2V$O z(9p*<#J2>gbM~r)y_^RLz0gzhErIn-qJ`MG+hutrI!}*2Z61Car)lCyc&Xr(Z^zrP z?sj-h2(AV-X8?PIy`M7|Rh6%ljcsrunuPi&N6R=zy+OkYMfdcSZIVx%D}1?@(w=8;*t(n5T_I0@?>+i9(AO!etN; zR(8kJ|OqN|bes zsXq2iOSZQ=IpEsb+}Q zwA0t`?Mk5*2?Q^DOFXp1vgDzsJ)45M-qDQHXhm^?Qc*PO(}Y(OD$II|0rg0wLA4fT zvEKH2_USlHv3mHKzg}x{wX?HlsH|ApTFOP-nR*fbXUH+?`J-Ucx@v9AvVBS%ZHK#A za^seaEpn}|3#XyK{hI3%`EfZ+EI**}`^>O{P5Ac3h!n(|s)}^ez zbi<;BkbxCKnsJ?g2*4O+x~xXGsPX2`Jb0^5fgo8D=}?y< zdDTmhC@#E5ULB4eKZ5A+=26Jm8U9cPEAq@qBocL04KY~56mMm4F!AD%Z%#}%? zp?B*&8a(35QV2d+HqNnXVft2Bl140L@-*f&`a)ACqLMjfyYD9}r@S)I*T*b$){#Zz z6A<;EPU^QxO>Er+46NKbWe!!}*P_h~zmg}1DHFB0Q(tAQHQCEWsi)p7wY$}a8HA5M zwu9zHM*BjqpD!&yzD&4|Nh%C8AR!NX2M(lV{cEDyF$-r@8j&XFWtmAuXZrF`hdHkf zW?xju8o?bZ6#BRd2{CWg1v8pGGmjXj{hZ4edpjm;!koO-QHc45`s9IB^JFut#D1sV zY>2))Rs$_vQFhU>HYEWGc2+hm@v^(h@yggO(_uxyba^W}y1>GyBZ}wjLerOXFT8!h zH~EFWQw)`R%KdnwJl#lAe*SoIcMpUNbPD*awUa4~;#7ZEs$& zcv|NugN*7V+>o_uZ{-dArq77rfXwzrDxlqcsl zV6i67cGkaYhg#@FVUarE5?<=60>x;7J))-xy23`l|V%`q{v zNcw4P33Zk_h0s=EQ9fC0Jpcato4;LMhWJV;GD07sQoXzd%c|;f+xGtDTU4^`i5A0% zP3)%AoT`pqA@5y{)D3>&Nmjo3aZ~Tj>}q`mj-j`dZ6-BIz|Rk(@Pmkdk`Hyg@O_WD zQt0isnPm@cZaVAt9YjTO*VFK$)4Kd@ulgrCG)w*x0q*ImJtHw+4{33^jWTds-b^Tw zz_X+qe_`@5qi7}U_A4xV?*MXF1erQ!3en`59m$c<8O7@OU-#$&4(tlbS7{g65 zvwmmgR06SRWajOD$B4KApJmFx9K653=S?{rQFg-0`@~N9>`YrR(3^IQ5JkC=bejmq zHI(|44sTXkm#gMl&NL9CA5Q7Le{idIgLBv|8+>(ct)Kj-tDk-cA9c;AKY_lrwcF0K z-(+15kY_~7yDd>{(~`<-xo9b=nb$sHf&&B|=6az-(-hX*xqJ2VtYV>?L(im5Oefi4 z3KAQmrEGl|T{tS6mbD{y*pr)(v34#`{Mz}9l zzR+9BU7YbAF(VbANm)~&EG|7Eyo~ccd~3{<{^{2pNh;r+5sx(-&J=uX z?KAGvcfb9{M=HO?N%WL^M58ip1jL7jXG+(IJ<%HAM2Jiu?UJ{4`(6oF^A+-tXDjW^x|**`33< zq{Wi2>-1KWHDK9NwRTE(D+Cs=`c5rMN4BFn%_^33aTeLV6U(qWhOj64y}cE(0l%ck z67yV~9Iiag8XODn7OagNzH2dK+M96({7KSwlA|h)iCLUmG`p{1zm@-1e!Um@KA3hf zUP2=MF||G~UKA^3DYtzxr$e{`o5IdySzgIUpGX2b+^PWWHPw zM9U)?pS&xJ-9kAUVMBwR8H1c=?Zo{T;G1L_@g>b#`z@v2iC;FgHf*y^_nU${7eCe- ze_SjDA5}5hA1Qosd6r@N#UXyAuuipZ>S}6>7tYtgb^A8>PU3yv{hXU68Lph?aUFw; zE0p6Q_M62uqvID!2^T!G*knXT+68Zr?Z;x?M%_eU6giFH>Z!y;caH+PzrIynEV_n@8o_jo4t`#guFEC7xRc{>EoWs!;^*(ozw9d;mw2W&bLTkv%Ix``OG=3+nURa@i-^?k z^zPGxmw1#t&*0GTrS!y-tb6m^FC!HWIo$5Vle!GgY}45l%7*iwjW;6Lb>6NG;ab#h zweW>EBhISaqD{VN~MJgOQCPs}Zlh(;XJ*R9(vbRC@TvYvMr{ChNp- zN70sx^^i+*&u7`x7wx!^q@?{>8F~r>M^k>LxOXeXT2xzNbr5Qx{+OrQ6!LIgN8=hQ zA5mHKAOoUjq-`gdh4l0=htnkH9}TpajAn{4+m{q(_G~VeRTSiGrkQnf6tAiAMo;{g z^)Sb{(nl&t%D}0&LeEaIH8v+=Yvhc4nyk^^;;5U5eGV9tIMbk=&dZl! zRB;kqTJ%??J$7TG)nyOY7}PhGGOK}HiR0(D_(LRJ2b#h>gq+Z3D2oKE^ zxw15mv+Ix|f+*E&j%*X#TlqfjbC#;?jo7Gk!NB3$|A)PIfotl@_Q$o3wmhr@4QVD( zENG|+M2!Qf<3(G6q$P=!JfvwtRZh5VUPs{0rGwPhRnU-!P#yuRMM(%INYrDKz_mUq zl7KUqx^yT8t_Iz@9 zcpi!Ksy*O3iZ<2{xhVr>1m{lq5UJ!ppT@Xb7mx~iRxL>rL^}sFY9a%^g$?#Ab@rBh zYo}FCi@dz>?3b-AgtLnL@0M`-HTl$SSrXr_{M_ACVw&{e+2(a4SL&Jwu7$&WoWkL~ z@Mx9>-Fx~ZoAXz3SHSwt)2;Am+1^i%uIqY(qZ7YZ*iwqDZ0&%DwzPs5L zQK0AJ!p5q-{IijR_Es;YV_*KIdUOf9>g>MtQx|iVPcK>2yOL76mErtcYLySDq#fIn z#TCxf11rfL`x86LNP@qv4bN9_bn60TD^r(Y`qnGE*0mg%x#;SQb!U4Bd!tW!?IVqK zj%n4?Aa5U7{>HlDo~D9J{44yCZkl_S+_A_-nYXKGpeJVT3saXAwq)gdk8o*2c!Qy) zv|u22hV($4vvD_LfRx{--6k&R!w>hAt_|m3ytujzejqa7TGZDAlOpG>>Fdmk0{a&i^({Nf-@T|$hYsY5 z2G-1tw&iAa5QgEA#LQuMJ=lS5KiRJ<3Xs+XT$G`&3>{~9xWQ$q%d3X3d?sF%G_B70 zWAu@N#$L{s*5!tQeMRs%e$UFqIf<7UeZKdlaBqw=vE=~YOKv=$EU+!=>t>v7E@~v} zQ*Bn4StmKz7|3slrYQ;s`!w%7od~Zk&#OCHlDa&`y=bHx;|L-L%cjeRmcO^H>^*^1 z>b1F+K(F1sXyfNaeTw{6w4o(@o2F=_tZoTEz?rdapKzh!Gfp?7<1bqmj`WdK;B2$MI(MedlN~}4WHxab2;5RQOcUgmNS{o&k6w@m#iDsNI6B9 z&M!Rs)%nG!p0>uece`^?zZh!n-DWMwtq8!2#N`7Wwq$umOIgf7Cu&|e(ozZ_yWQMz zAkwgCZ{|RMG{x!fI1`hOnp2n85B-(HMtJyyM!al>sl>u z&p_&z8%74et~7<_)|P#e83q2Vb%V?6%c5IH`d+p6z3TM7EGhV!adt-GNH=4|ezdR# z(0pY6*NcYxT#13)mC5V&iVGTjBM!{Gf99gVm-gY@=!JVvN&Lxl_Dy z(MU^~`0Yp-B$eE~q_8tzE9IcRfTZ1ScMkNdq|+!vgn`Spo=hkC>iL&f?kg<-H3#~) zr}nID9d`DzKP%ijV_SI(k+WvuKuq6P)ETIfe?A2u=Ixn(c)Hb)I07-ScQ;HFise1p zZKn!51$kHA6MT7W-K7f)N9O3{LtJ!dyR%Qp>2BkD88uBvSFiA_n|Ae$qM9sJ$A9DB zqK~1z_wzeG9USQ3cQi}8>6XHp1Mr=e(?{1&UFcZYmEGyHSCMK$}K1Kfi4 zm$OlGpF6QJvFZre$};+^@DBMq@CHFQeaJ-{nwcVzmu2gZ77RsO-3x2ZWanmICYzJl z=Z^?fg&i}v3j_O8_vH`tluEb-PI8%W(0w6phID%}hY(O})?I;50PLN?u7GwD^5IEv zlu>c2#~KXLEMBV39zAhj)##}^-N(I4fkq4#X3l#$>Y0m^k9M9zU- zl8ys}y~4dDj(*A3MICq-rNOT56?H^gJGLMGn$h?^=YV7@d|&`3W;sQzYod>ceIt#D zBk|~Rz@+VDom5ZL6~L_K^lj562QmY=j>t;^gztbC^(~xbseWW7*=ucyrVMpM|64xo zoz-sZ5e?+>W-c6oy=bhIf3ahF(Qx0Z$F^k#Nc^sWo~9y8{(#=>4#<_#GK1royxTrF zzv|^u@y@qro`z951yGHAQpYSSfF=a@ef~pk3-_uqcj)A)lL~$Mwz=Lze9rN z#&oiGlUr6LZ~Kef*(?5BdN6Qc#k5-U;-ap4H>K}O+TP}pMQ4>ci-rdxtB&Tp5H9hW zzoT_zj<8X+MquUiM3Xx??)QtF74S`hD}j<@CG@MiU|xGifM!Z?v`VU~S~cDfhVmS1DyUp!_(#l226an>_53r-hS9pIOwiw zA<0`vl9@Df6RFT(7<7<*l_B|#)Lgvz!>NnT`gRo#x@hRaz~!S^i#Z*(#p_yT6b{4Z z9|mkk(YXtU%}2%W7Y&!b9i1#~clHVg`<9U}y$8-Kd1mUOmcK|c2G%%xy)}Q%Z_U5b z%PB?Y)i;taFW#oN`V!rV1N+#%mV@pEkymz;yHj7awHyre>`zn=92{x+GiQyB!_I~< zwkBENFYL{1*?xo#u5KR6&AhaWZ3uLf{}@^GvBX!;UT3z>i+XyI{qR7#z`Y9qK>&EN z=<3Wxmn4$@gJ%y^4DSVQnlU`kepFgElQGzHVCTBP<<`OFG4PbTa|zm^KPq^o>hWNm@#@ZK^p zBHCMi7NxYrtVw3qB@UJkbJq>}3vxeh9d4PqrDsM#&w)BvfZlJv#2@q?pL~?_)uQ3D zooGZO#cZCJy%*o*>nt4Zc;;zL&*)pT+3)dpla~v7-z9MrQe|Nuu0>OzS_hQkeGa>@lY<@; zTz>h~Hlkp8#g(0f4jW(|x~OYsf&Oqc?2q)2I6Vi14bJu>>z{V^b&AhZTC&Rs0sTq; zzme7fp5&L+ewJL&_bOwrXt?FssYeU#ABN`7Jv_iheegi^PGlt=;L4S}ZFlO}p&3#|_+;edQ{WHnifhC=cDetX^?q3R&91fJ%XbO5xo2gPA zx{UD|&F3F(H}}N!m`8`5)%S!m#NKcl$qjcRMZLtD3;y7B^fxAp^th#bBy#UgpmNz(S z$c?*`-zL-o?AOXmIcpZ}jqGC=tXTvPYz{6LuWKoLAAGR{-u9Chv;f14=$Kx7wBS{` zR>H||IS|>ka?zEn!p8U5?-llaRXf;7$?KbecEE~R9Zey~vj}1ax)+(vE!z)weO?&A z_pUVeeF?Nos<(d4XgO7|_l=@qEt`OYtdI~K1y zyJq3Y+$Dtr9i`G`_8+Gg4DFAL9Bw}*zXD5oQxErEpj@dtSq?l8cy(X&;j`V8mNOE8 z+t;Bt1C^au5pdJjEE>!k>6G}Y&ZbMj3oca*traL<2M9S_UIp)c zymyxL_KbBmB3r%~Y`|JHeE-?pD(CLv#hkC=fJSG@2Yi0l@h3kPblJ+?JM&$PE183rBn2Ix_Il5DFUg-*H9VKrZp*E+u?xy}FT6}IlN9(W z_S($8-HS#polk_M+mI|R%XHJ$9Odk{!5brMpKiN|9#KTo=84~#wwvrcZ++F8aW*#^ z>k&g>8!F4R6b-jZZS=X;j;3{u?{E7vCwHa2>+@~ra z+93v89vRx7+5;Tfd$=a1v?THB^fg3|jWPTkP(yBhS^cGzl)j8BwT+*lw8J%@4wlsq ztXZ;l+Ny;E{ltPxD{0FkYqIiHlE||VymO}uzIv}DZgwt4bHV*eLug!Oi4q1wNAVJf=pr0Zz@9Bn~g`%wi^dLX719U^ja4VNyK5o!y{ zstk1}*DM@fwy0&l{VZ^v(?<)PReS%uu4~mI$Zk4p^mNqi>&g-bB-6PI2bJ5hiuUe5 z*>l>eg@k7@zxK*`n$Hlp?7UL9a5$omC`Yk)sV*(QQ*t1htOBoY`NSFMIS_ft+qx#Q zF;L{br~+h2@ee{~)J^YBglo~Q?p21fLDnO>6E_f=KRnHi|_yt7x-V2&7yZj)@2WYGrxS_BMpC9SW& zeB`P`5q*riczsyWNJRTFao4cp@bI!!6`Hyw^tFo)XAUmuye$&|{DX2|7-< zUNL;Vn1gGvIs5h-8GHAWX-m+d{YTM(vKj2mn$wKNUg_l{((buH3n2LQoH6vRIVSH* z=YLkv^U2B1HPK{vTMaTsMu+e$-auQHEKKa{cmuxvG~(V_IMRUSLQkaoosC@``&-?0 zsqW9$*`0&l!voPL%ibty1XP7AI2yjpB=0uM^lcSa3@wuNPcPaVO>cGA6Z=aJ-{D6|tJFFE`*{c<7%UkksE+{)=c6WMW^(J6Oc%zbv+E;-sJmdNyu zko+N+b8x6^`nLXWe|lo(p5vQtVC)+s*EruT+)z60O#P1uFS52>e`TX%+h5Cnik{h! z`d#gY`9Gu6PEG5ozj=P9`Q&$<>woCDuzQ|q;MmRUFAOaI=*1^~LVxTV?io09;>9UH zt-n6$`)9vi?!R#8#ixF%zW)2~|8#Mc{l)Mb@!!|{ZTZRL$_@UX9M}K&?KkZ&=5F}n zC)f4LA71!-+2^Z1d@eoE&Ay*{uAe#m;-a6|kDR$S^V{Q>Ho9*N-F*4Rfsx%q z$N##q>u2lFlA9er9k{;xhvQ3rKJl~m=jU(ce){5f8?XF)?#E~Uy2`lm?>}Gs`7b~3 z{W;fj{Km!~e^`4lV&=f{o(=goPhPX!u*F0?G3~p`r_`H&TXy=z?u|`1jO#!A+u85r zvn$f?`dPp0`?|Be?<|-9ZSh`v@XGQJ(`Q%yoquBR{S*7nJm0ddF8SN5?<7x49$fq%m+mS!TJiH|J1PSY{{ad7Yx}p= z6|U`nlnlFzC!IX?a4{ZrOU1^YUYzm8(3#}q=RfGED;zo#G8`d}_58*0Tkqc4y>()L-R(N>zVNvHjW6_rRxZ7hJk+u6*WG;WKYu|J zk7@44yTPHX>)3MUSL*lqJ}~GdZU5C?-6}saGkRL`+Qml-lmBs%@9xcY_^#FNyc>N# zl^E$9&NCs1ualdKchlxOXz>Dqn>^ocW9N{m2tu%%OQn%6iiUrY7b4D%=ljhm^)j8f z%4JR8O3P@iWqa`EHok#O71xBYt3!n}DqG|v%MpIHkfPd3jz@g3*qqH}^VJTTJWYTz zP$Du~IWc)oil*4g=2vwhK6?|co!J;!-4qkg?@}8m8c9rfgIUMVwD9>2)^xGoIhW92 zBx_438g^`jAw}q-m-4F`+NvCwtW-iUWABEhA$*s47S2s$r%7WJC1y$qU9Az%`dyfS zpl;<~RPYTBItB4|vVG2!+GgTuHR5>)ZIGcl1hVwHW3qUz139gi?2J5?8vENNo1UT z)z(y8ds~-_rY&jnHn4(9x7&0`tuw{vO2xTox?)=bzqX;RIz(I@!uD3DRM_`F8_uqe z7aOYEF7iT4tn(d|ctm7g5UFmY=QS{C=@Np=rWE@dFltF#W;Mp;Qs!5;b%*j}g=UPC zMOKM38&i71#3H*{#TSH#Ya3}9)#UPa8>fV(;d@=Q`Q}D~!#EdTP04_E!%A!vj{v|H1KG`0kc1aD5M4ll;$gJn3(rASQ}AC#2MMJin8 z3^#=>eKt(&!8vSNesw7FS{S=N9Xj4t>q_x=@@rjfplmKDHd5`R<-w@(eJ&dT-t?2{ z(j+ll@kO10P@U3M%lCG2aIJh{H3iI*4W5GVA{%UZj#!FDO0YBAz|tv5tqV&RR~TqY zad!yctFCT}bkWi|1hCx_tB&ojvLy&&sBCXziqLFJ;NTtf2zeHTQ$m+>=SGF_>nVJp z9a9yzc?}e*xH23$7lLFf5%>(KP06L=9kvXYwX_6$#}gI8CCf{wk?Lw2mEYA!1DoZr zFO~|Tzy*3j`4@xkzm=?mPP@tFW-#(puz89`ENXA7aFG=f0(dz10Gsc!DoU^fE+t+_ zr}8t+WVVpP=ESJP^$KBJWHpQ+oz3?-DH+aKvMOgmR5+4p0Ur`)nyG|FN``@~;?iM= z`IR6~Mkmq}E{JSzifp7QN}32@_-30*+~p+eBrythO&B7yu>H^=yQ-62ujO~OTgyAe zA`pkWAVyh2*Ye-RgLjpr63i!Yc1$M*T~I)e5HBi!XXp~|so-cBk)ujTc zUztTyY^8uo+iesftrYRv4@EjCTqM&ShpVP>gmFenMTHd~dN#bo$|<(Si@Pi!hB=;pF|Ueg zV`tG4IFW9;76JSNY*pV3i#-nrUo|Mx31a43xa4^5f=D|yMl(u-|ub}cvwJ3FI+MNuGC4thGfo+9?uf<*!j z;O&Q^!b;5Xf*20p;7;+m!OPk%%8{CI38CHgN-3WYM$(C#3lmZ_0$f!Z1`^3(m;pSI zE?0`zWiAyE8m##t+yw+fioq32nP1E1{SF|%!3w~pP2)o^y%5_!-eXf%PbthX?WmT~&EbMw1$vl9% zoVmDadVH}puTz3kSG9?PHf=6iCeZD)z4>cn&!c`}@!AO--(O^wL z>e*>*ezm#0mYp4+hE$oYvQoCl445QxU|R9HQ1OZoHZQU|rK(j3aFQZw1b{G@X_~F6 zfT^mqoCWwsin7E=(QaKps6K=@(<#Mfd2t&8Nwi#C1?b*x<_ew-&6-Pa&3zUEK9vpO z{n@Z$n}#i{71x?+Tmiv;5^qS2z`a`zPT+{08(Do4{Lhj_=KxeQ!J^2jtQ4Vxp1|i< zvFI6gn_TQyo6X7`I=8spO67aNc@ns^csAg1iCGJsz#vMNinHe%Vrh9`i9rhCccwZ4fCU z5=!Rc9h3|^rp+OP5oD|Qvw%F1OQkXK9DEfv-$hGc``S+eJXP5-s+vp1HGnw$sy0Jw zkbHqBy`tkpg-h_AS@eifNmQs<=x77w&(TOC)vhLx@?59@?<8w-$SS@Um>s(&MEYzv zBCNKSva{3q_42Lc440WJr0c*GZSh(33>OwbMe5b1SrFHmRD88LqZ%lJ2mnH!m9irQ zG1$QYF+~oAnugRm+N#We>{bpxJFm(dF~2>{NapZ+!ocPrIN%KQ`JECF-s@`1)?`r$ zcIfoU$XX!`2z;0aOp*+4;il(RQzF;^hMf?dAhJ_*Tr#(WBIlAdh{z5|ev!(L?Sdd# z-F7y7D~$tCYo;j>QPs(qbWRK#@i)=|*i>NZTspY11&Q@Uwx5Ixm7zI+ofIV-Zl>@r zYO_+K!lXEdRU^Kr;#UAnTA+1`iXD3~0bsWQ11{VNTMd7uaW ztQacdhw#D;%C{q|*vu7+3_-sJO5t9B2ls2RQ8>?pva|E5F=`g@O53}k;2ti!O&(?u}Hfj*QVV&M0=2u~*;)}0n5JNi_pGDEJvw7iCyxK)8 zw_^(aIiUX`IW!#y-)_ru(iLj4-;g45SV0G_lQ0tTFpgxdkXGK_h_8b1PFJwAHN`e= z4jIydNCS{;A*2~dY&Mlk<1QoT+1q>uAg^Y57M;Q+tF;oG%f{tbhi|3CLr5=fs|-P2 z4=*{1g9s<)!}CzJ9DE~9C$85CaJAwpM=ZVENmc;&GOGkJfjhHHfH+$-Eaa+4v|Xg8O;Qwc6JFi4%)?;@uc z8|hp`1VLE^{13<#g%?&rE(cVUOJmZD+p0si(kN_0^&x@-`1t;)5DrdlhkT(sOdRWX z(Kv|TNd}KIr^YDx^(uCiL&G-`8i26_n}LAmpyV}JRmDyCPK4ji#HeCJRZ}p@$e~0~ z*;SB`u?1mS^jAvR)#2=HRn`Jvr-I)CVgVH9xydlN;AoI~CGaoG!QVMJxFb#6rMB-U zG+=CYy(Whu=dS?6kry|D*K~&>7b%h$EinFAGQ@QQStACL$KgQ4X9H9TQ~a$E*qd;U zl&(e!wYaUTRY;+1rRmr_Lc0})FG2^D7fT0rD;Ao~Y8~6qDR?Fn$&MGtW?D+Dc|ga^ zHmVrdf|*)k<3Q#Hg4mSBSh?8@guFYPU0WT*7!CWPqJ*4|_^PZFJ`4dR&w;6cJ;p`4 zQ@~^RcA(8wuC`2bV@xpUv9mdld2FSEX9GP;09VMNP^5tBd1f;gdcV51xUJSnjkE~k z2-UH)`3+5YATL5XPz>(Y6t4vsD{_%5_+3~QExm-I6#L96HIN4bvR2tB;;O2&HsElz zRe;mkN*HCD2BHB1+NW}G1~LRGzU$DlVSrCTsWEmA6+X?M3H{s(x$_}Hl@W3=nhNX( zLa?|dT(}@Ai~~Nd5&SNk@3%8)T)x*qhltfk*T7h?GppJvY6UnK6>^1a4dQpCi0mhy z2?0V5u}$rw=ecQAq@t=3hyifda!6v~hE@pRbVvZo8*BhoN_H&>j_?R}i0t+0Y$jcW z0DjkmK^QdCBD5TW*_sELcLKbe#^*(aBGB8~T7EX0e>R*GGe;AQG1*kejNEiZ5L#+o zbjSkjSRU|fS6fwk8)W{O?X0=@b{L0xT`=#rMU*!ffl(Xjnvz2}Ggx>w*azPcOHPo+ zyrN~tc2&=f(L%EC0<<#NoAAKns1UKBr)IiB%m>N>nFZwN1_;G;7*6mas3V~m(nxYC z2bkkUB>+zjO#xO?j3uzWkS+mN!~=yo&;GG7J=sIMW#ifX0Q~rJVz4GT(tQ z&77>bNQYSkM_N4ku_ z&G4;=s2X@g1rQEFq#6?HyxJ^Ega&8?0?cSZB#;YfRCo>@vJ`bAMFkcDtP=phUrnc| zb0`YwGoi(821i_^fg%T(4J)B&5n%(CS8dgS2o2U!2zCpe30I4=%d2hk8|gYJu6ln= zdJbJKt_+dJP=M@8f}QDV2gaXLUn<49f-LHy254Htp=Tk~%Hc_?YTGDGma>H| z8W04&gBHOq!O9^{K{lU7X9q(;LK+u9ww{_Memz8jGq-isQ(+bWvZFK!%xAH)zzg2R z!L)^$l8#(V=K$VQ#G+0}O|jCn90~`RHGdYty#JXn;Ew`$&PJXFz9hg`ljmEIDvOY! z5c9wc)ZjNRTb|2EhbRz{X9q2&_zbiJ1hT(pL)pGyOn}KJh=Nt%{Uug;2|1&|s^H+N zk@bW{Rx?ExAUp^`?LQ*o_=7?|=1;Y2|G0`p9e6zVuJ0LfJ+VctnA z4D#>tT2Pdc9>JAH0s<%?Hn&{__R(Nfio2?8vQmC+JI3Wa0~16+=2M}D>ATTXJ{{7D zwt9AnEkct8W)A_}MuDgm^c$MC*qok41`F`H%qm3WNU3V1K{T+lKx395^4BZHz@db6 zIlrnEya`gr409^pkODwh1x5ui4Cc`qFwNGv&x9c>&sh5g79UyvTKD?be2ILd{gR6QH)Q?-(52r2qU$wR-l|JpfA;Ttz0rm0O)vKkjvV;(_+{j|r4hfQ|Mowg zN_+V$a|Hg}zdASmvG3Z=7t2O}U5h^U;#cqorg{#nJbLe+w*F@yq0YeSUt8y-uRePE zqhDM2-f|P20d&*7r~cbV6aAA8PM9zWdhlPYYSF;{!vp*O8{#~!H(yWm-xBA3)lLbZ z&WB4%LY=p7dU%2V*4d{&nz3oTf~n4cYW(0AcYc*%PlGZK5A6BH%I+!fq!ev=XU#|V zoV)+vv9jK;-TQC8b3C!|YvSLiC zVx4*^`za?ay)s4SV7`ea$QGw(ABo2^_9iy!Rn<{?e4>)6pzhhivYc2k#f>3J z>9TGiDcxf9kYtC%PkD5iG}2NN62Md2rAvtxj|w$Kk$9P#pU?@#iH1`=Wx9e9M$}j0 z>3XtI#dIc+GF13em?c13ZI>`BnS`d8B)+JZ9+8=_0uSx;niZiYmB>NjA>DSF#?_#= zc+#qv7CljKCvtd1YQ?#o#D8sGv4!Q$bjjS!V!bIl;`h(YSMMlJONmiz`1g*>8(tp& zIX`plJI7yp=ji!w{`BGJf4>aTsPYI+^5wZio{34+;e=JHiJHnOM$0*!cPS6A#YylL zdW#XuT(%Ukg!7twO;Z)Untckb3Lnq2HTmKrFbOL1N)#@KLWhgeM&W6DT&y{kK$Mo) z&`o5dDoPVre#|zxQZ!45N+Xa{7K@XEW9^I5BZ*3Uh1$gOJ4kqsSdU3aqJVP!B*y>I z=w#=T5TAyUr08|KO_}?cN|8sI>0$M%eRQdWsP-zRV)7(rx~RhMGl}V>%ZwX35A(H%wa=3#$Jb3sgg+WO&ThLfgvVT-8;wMYqTqwdwQ2eXj%RpP8Sxd zNRW_7T-O$nGbw^eETvL8IT#w@bRm9YvrtSkb`#x)m@*f(l!s7_bNpTnsLG7L*MDTywYaZ-j7?`~GmSY0T~=u42qd6*^=;tLm@6413c#pY!! zzfYfTTB`OND^*GW;HGPnkrPb|FENC?vXcNSAERCcD8!Vg;Bm9akDz z;UN+AIECKWMIzzV>csaFgz6P|w>{IXH}NW2nk16Q=us^dHR+U5_%0}C^!^2j4No^Z z=Qx)UK@`lcL@;Z7)MmY&XN~oHFeP!Bk)fCJ<&j%h4jj7K?L}f;&|DNr`g&NJ%$c;A z5uWMxZ`Q8%G1OiZpQIIi>@l%Csz1+H=&@#@NJG^w6*Kg&C&|*fiJoS|9FN86Qt|Y9 zr^+!DF(*+d!87Qn#sddQ z6?8=`!}O+-NiZRnIXg_OP_IIh$d}Pr%?5vf&?UkwdMwJA#3WKZM$G8NC)PF#5$7wu zc>JdlLbs+%On0`UEE7)z#{D|OU8xd>%Pd`DorEcsol?+v_}R@01|vb3mSfRuVVcBq z*jl{XmjwP#_;e=|Kga%Vg~#ctl=U*)4L*jJ;N@j9^!QkdC6U#hX_A7ucu=V<*5l9{ z>4Zr_(595fHM+YdV^I~okw!gnvj>dSyV?;qbuvlT%tP#A_Z}m|VoFnq)!sCmoGxlg z!fRJM^-{#EWHD$`#e7tQks_ol46jKFX-pR?qlo$n%&y{LsK!W}EtT1KOqGz9qLbs? zMw86b9AiurGfvpWr$l8c-Lwkh9tS=iUx}!DQ#^+dyC@)yA|Y;Xo8PDtyJMY+Vs-oLRXfR>O%d{0DjaEb7MGMqoMZHu z5Ni^$=Ma-0A>{4Bny8FeBLmZGk_dFCPL+nnB}GILlb+yZ&r$m^VkN3>S4!ncxavxg zX1Pva)icu+J7m7}Ia_#`lov5usF)q<)*CUs4!?|e0?T}^Qe-i1 z$>~%HQN7IX(C3&YPDxHx8KPVcESrp)Lmw&8U+REv069;e$b1#r+xSth<46qqp4owUc|W9W_A2$d`W$7&S$^?HjtF*BS$d8tX}l7(sx@|ssn z=6Yln6Oki3#oPgdr{8mshcStXOCyXsoD40A;4(2;98z3t?^ML-w=#6qcn$RlQ#MU2 zCF)6dJ?5Bfk`U)K(Ja1%Oc#J(f+C@|p=nN+k{BO>*LcV*wn;Hx5@As{_}y;TK_ScE zrNnP`=`mZ{C#p;aEh=KE@I52O6#4BuV<2iWW+DPiCXx6eugI4ag@qvP2$ju%_!W=P zg#<4eMfBSdh6ZHy2w4dThe`Tf5gXLxrq^(?)`8qdTVP2$j0l;wek94rULT2{NrjlF`OrJ(J3D*QjROI561 z8pgyYMKtN9%+fiXOL;1!lC+c*)G&w?2tcqryLIWmpM6jjOGW1O>(M4TN6v3T@Rn*r$q zOZge#>=CkVjk|Y}EGccLs3}3A`CwMMkh!{)mJ_RIGGz>n(c+6Iwkyd@GGFFSQXdk^ z{9aVgnjb|pq53V-Od65KG?6hq!KY=Nuqzl|Q7_|^5aBG9HEmAPxb4Df^{0F$lLa

?k5fbS^+GPHjZuNlT;1%qV<A{H3;bktf*jM27l_cf|NhTEUcHwv{qDT;65XNTzr!t=@4I&{bbw=qJyZ$F7ngCIy}g^gr%mXP6|WApkv}hl*mYs zIlH2CnR}>Vcp8FXR{HOGm{Z?0i)Y`-^%`Qs(n;Ag#O=XIacsYhM9^EB4OkO>6GkFp z&RB*E{BNmVOAJSu>8MGOgqegY&yOuMV#@3QEl`HC9Gz|-)-6cVZ+2sUIsw9 z$Nn1#rM-GRS;zA7kj*Al6u6E?pPut6)1ifqHc^wN&bPnW5Oa8H)v0cf_M#7`A^9tHsh z@iZ7~RlFv9ErTv)ZWYrAcGPG*p|5a6#p!{9S=4^Lu`NR-(r-?Jl!HfPY!R|%yY-lK zDc-Nen`SF0)g~ENJk^d8_2i{^_f$<{Cbb*`D%vjEvQ#CNDv1_-+HY68GBF@BjO)Lvw8;++*kcc*!n7?P85Ysi43dUwDgOQpE z&a%{=(~aW`RZ0ZvVmAL)=C?S(Y;AbGOegv@0?YA8Mc!(bWsfDm@_P-36lzZtPQY** z<3V&;lUDVv^Ta;qli12c`JS|g!~HCr(1eMpR$V%2~|0Eh~8j;E=*);suY>TsjD5c zcPz-PPe4vM;H*#q`344PrFY=OZZA;eIZihwXSh!Y5)l_q%XB8OyogE20|fopm$;MA z6#vqD`iQh@ZInh&Vix-?^D_`kL4R(kDTKso&fJ3T2`EXdL=}RuW`S!ZL39#p5h`Y+ z5lu3S?>Q6CL3IB<0gYoPC`mkv-$R_rnv=w{6CetZlq+(#{_Wcf_-kwKWmTuwy>N2$ zr|ss5!m;8&4Kgc(C#L22;nDU02|tzl5~jsPH6j+j+r%=u&o%lRrA(a2L7*;8Vj5i< zKr5y;gcWPl(6pp9^fwEs_D5!d(ix8?)r5{27oHgvgj9&enM9g)kttn@rxA^=xKmU7c4<1Xa)>p$+4#^9D2Lh=)sI` z23=+XDVe~Wm@L;J7(M@U%|J#G2d`<)9>5vUg^j3Xd6`Tmi8Rq}h@B`6fJQ}=u-^MFHGU~4vnqrb`QRG}8>S>?qk3(`>Tb-~n}pjbMfEV7g) zz+w{ld$~fQ<|U_|s0C64>Q=i=WF9XglOE@B=u~*Ul1NP=^X%_hJm9N)fPF&NiIIN8 zS5PaoahMUPTg5qsNf}}1D;SA6Nm!*cs5}JRi7Lfx4^Q@{UdnSWSdk9Uj7_Q7 zxshjrdi;E@v6+w`Qk7s za!AEnp|N0=Bs+ko20V|Ih-ysC^k(rgeQAPJB}HfRcZ%#JCXQ#9X`>>fQB0w%0dsD~ zR!n9!zpErPXEHo0gWD6XUZyrulMu|tj3UlQuszMZN*THZdQ+hwIY zPqCIk&X3pZaq$$i-wP2Jan^G?r9AbiEi8Y7FP@CaoD6rdeFZ@SkDRN-G*EqSg;AqZ zco6C)YEqO%b8af$!?QHFG$%4&F-a%uMHVb7E)iu)nGqxz!+kD+FI8B2r&gR|h0=(L zD>QCzrk*KPB&B7}LgVzTIlVD+=6JNE2xf$_!r@R}{4fZO=(iW9-bq78_s(p+mb>+u zQXEhms9Ay6EVEd;h+CLsd^`!SQCs#n7s&j*F#Z5J4VAOsZJMgDXo!atg~`qcQ&b3z6>DN?S(;G&QkYPPGvO%-wRYL3=oSUT zWqfl|npu+s6BP?&?#=o=kg4(%e<*sDobY}WP1zXZ_dgx5_+Wzf-yz-xSq#fV;m*BO!hwAV2pWx4?zafS_~ko05@-gR7%Mkn zY~d`p!q#-(w(pi$qxFu8Hi3jX?<+KVc$dI;?s>4_=sMnC$7&k=J+S}5LQS??{RrGv z@K{$IcRDb>mwDIP``9`wEVUN&j9(A#j4rSZN^Twr}^!^VQ*8gx}-B{y{fc#DWJ)*od;SV2;7k>Ir!LsYG-WKw% zO5xPPiKoX)PI#yUlgr0C`4|aqb$_&G5}E)y>^hCMNXL3NT77>|gdZ+EsIuTaA8eZb z#kdep{>!M!Mjoc$3FKW9`l00b*i9cTd;Ri#<8Z5sYd4U-J-KG@OeO_85$8g1jyxJAD;VeuWCdvd}b9yKJxV_F&W zou?=K+oLRPVAHtP7EV0%5EsgJ2IOzOu)64r@p5<9Hv4XWTR8C-Yl8;9cyj-ERTI8= z@=-dBZXGW*X;d2zx?n`V=SRukzv*Iwd{(a&%zc7%&vePGAJ~)X3 z=25zOV$Lmpx@8!@PW2;oGv`msnlNGNDDcl8ZU6RYzV-1Y^3K)w0OGC*x2$Na@h@+I z&ON3y)&9Xw^uxPta$*a4C)$E^-VH#v(JI*Zt!?5*Pmclgy&ahla=(S&*V9o@4}#99 zY`5d!?Or@!z`ss0S8u#<@yLPdN5j!=q#D_9{``g1?_vAl+uJ?)rT=d~Y;(KZLow+7 zLbnUwXLzVHuzJFz3tx<3?ANzTO?a|(+S)ha8@BRM1#Cf`26WRM^?tDFS9nw8=?RZ; z@_StL*C}))Yv$}xqqyBK(jQUh-)`f}2VYJg?^N9t9`CJm%Hv@au^*_hf0>D{T-*)q$jKTgq(C%$nmn@(3YK)b`?;CdncB-kKUXB7PV>o{B28mKnahxq}wBUFC;vanmmq7bWXSf z9gl@S2&bduc$8^9w&eKO-08rmBJM|nDKIw}7zH={s{pBBN2_lstIIaov5?@ zE<-{Szz73l76dXKgj;SAvaU~@0py;bJHc?<*6u1vlD`>H-HP8~&yD9oV}yIh$ZLt>(U{GSkM141dgs)CaW5$Q{Ra8YjTxg!$wTbu9t}V?^T|P=3g`dR z87LZW_gAq=#%BTa-0eto*8qgvXHHXF7=gF{i@bT1ooVkV|MoBiC7cMNuDO#!^@0&( z-ywZ>Y-+*;u$VuN<{uBVj(;f#NK0Ox@PEcL#NfC9gn*5pQ@nN zaQG9q*x4=eHJX2nI`Uf^R+aqlr_Xn1e?RRJ(*xR#2fqAh-u9 zeyvH3YJc)WX!Jd`LAo=l?|Te>0-U$M$IyP$pi@hZxozW!34&5#2iPlyGi62 zG>?+*<8imC^`mD4hJrL{bbAa}x=l3in8PCo&@Cgn_1|dmyZ3M19|Wg6RgLZ~zWMZM z+JEe$@eA*4dgm`ehWhrc(!nNCzzuXvUlS&b;@i|#IJ;^5rs*%eu;I>8(I}Eqn{J(1 zxM|_vX;6Ip4v-CP!QU;-lwmk^Vb6Gd(RV7FF#VVs za%aZm$vso=sB`+7r#UA&uPto-!@MU+Tl6Q~U-aeWw=Vqbp&9&mR{y|S!ToCo|KRbp z-G|hh;SWRq^jz@IL*H1X@!fhr_1f`gFWtL>^AEOg-)+0=?ydFq0@98h!|*rD?%9^O z?z?x#?%L6%nkD_Wmuc>;>z=#+IR{JrviN=>{z=D2hWC?c~K@*ctlj%{Slj|FWL%EAc|z+L4WK-K+3h_WCPK==J`nZ0F-M_Iq7nOcy`>%Px>TxNWE%U(&n&;>RQZ=a;u1jy+U&I|d)tQ2!sk-@X3& z4<{cc+k;{cf4_0vwO^_HnM)oj#;)=1+}=5};h_b8g{pVI^YAYHid+BdMRu2582eXt zbLlni%_!%El9saR|ChZx4N7b47Y3ggV;mDvoD;7iqT)Qq>D-8bf&b=+lW&BREfS#6BkC5FHV-y*shf2oBi9b5QK#DQ(-$&dYlMB*mhPd6wH3{dMS2U{1F4{sV&L` z#sc~h*C&~zxF%auX)dn8{3e{Yi=arVv}|6*>lp+)OvEFoX_+neLt9*OCP_`LBFqRviRl%zYOo}B+}0JuO>c8uG=0{9!ox1*|}vuDz-Nl z9MpFH<}Go%eZW9MI)C}|-Nm0Io6+yQLdK_9agHNp{}R#u!t_@|F@@w&=>YhU|2pU& zvysIa-!}tl9-=ZvzGh8tW@#at%blNQzov|yEi9)>uP_T++q2_eP29Mx}NhqXZLhJeSm*(c0U>QH~yZR0{w00SQAUdV}Gu_ys1%rkAe*CoiMy(dnGG7Ga?gh*+dVFhvY?j9FFh@^G9CI znx1a2)+H_5)V6_pJd>kum6x9K@IO5T;daXl^#E$S`JLcSA;P$ngk^Oh=t9-WJ!f~mO7nhYH2$!uQC|1^{0>)H5h>Q=Qgc<^`yG;;71W)td0%O@>n zGD->-lFOd051;qw4=j|tPPh3lyA zR_45ozyYPm$+-QbT&UVs_1ElwYti1@0VAk`n{-X8`_T(~GuE9HdQMpTa`FA`yIs&j zuW$Xuh4oUMe{6KvB!LH9e5q-Uyvnyl{iS_>V}3+F4@w=?CpHZvk;BVSn!dO|DG3{f1o>~wxw07jPBm^S}D-sS7u zjXX&X`G;$gZF$huv6z&9$@twTu>a;$^)C}pi=rQ=nmrSzd=0mXLZKkl?aE}wVPv<({ zHL_-%!qBwYg{CB3(?_;L2?&@FX-@y)@iUw!K$+K_e0g$fS`ejoHnRP#q-YIc{IHj} zzvdj<@0h8^<#c8gY_uRpXO~)VA6&ZsL_vCNR2FK}O}o0ECT#zUi_O>cC6{m9ae$lq z^sCUXIY%qio3GnvnxhL-{fJh>@smuw@0kl;39%H9AK!ECYE8=1(if_hIu$%! zFw*-)xU#l(`ZZFIv^&mZAQH0H1jR1B>RJFDgJvIb5*Bgr}x(f z-*&YwUyn8!`KmN{`e|Z7KksW%!FP?hr0o&GukeyG-QrDCgJ09P^L|ab%)IIPH3_=< z)AQLpR_O-`bs4*02bvV~EasYB-9)-=S@(O2=MySHpZ8u>@->6H_SLw)th0fcqq(O+ z8PCs}ZL-0ClEhz&zVAEy@cF3$`v=**mr-r=N^y(ipkTew*8h3@ua;#$CrfX44<_Es zN_guP(VAnX^dDEVD;7UKRLY3<5cjB>L=T*(muhc`0S}c1y$6ao( z4Lx)0J?u|ys2ig z*>b6-9pn&`Vbk<}WUqANtgAZY+fmoo_`S2IF2yExUS9_DEy8B+=`KN5MB9Exn>`&} zAyrI+0XJK1>W_T)oM-sVzv=v&pP(-#HKTM2Q~vJNBgj40FO91yfbU=L8F{mDK2Fm2V=!vfIhRW1 zm(IBl%C5!(3?jC(|3n8^?_cLXye}<3DefK`1?=KevR3j>)BuVTiOT3&wp95JZzv30M zfPMe?^gnu~|8qH`NKyIw>HGip9maqD$Nz5;Xqx-K?c9Gp@V{aNzL}DC{&MI3pRe2e z&pP+N#xFDb3c>CF(_5YE`ma^v$G;~1_-~kBG2dhV%bI`x2UT6{(kO=eJe<7p^Pm4) z;2A?PsG+K)7#R8Ur+-&0RL?1jD=M{eKzP-;OGw&9Ovm}&i%Jj9-N`^l_g@GNCo{Mg zF-R(gQxiRaRabhTyk_$}s>^+#gj3pDrBaC;+|cWMV0})hY~Wmu%O*)}J zL_w*uyCcqp+rU7a^H3!-bzrOBjdQ^B2qo${r3Waij&&Q13q&DEKXeoyLGe*AHx$LW zbLxK;)_)b}VGsUVoNN7$@;{3ISa~k|zi;$c3H zJJh?>f!8jmMOcRe^|1EkjEGWj4cSdA$~vQ)6+RbH)Zh|&e)p-(XB8GS3R86#qDFpm z-e$K0OVg|jf#{w)Z`F;C(o&(hX+9X6RK`BzxXnAexLg2~syBAGn~k1(bI#&lO>(mC z{9Hfv+yCA)G%WS~eWB`|gmVU#+BwDT9>tAV#9xWG(;FPSV(uH@SM60w?KHRYd+Y}1 zwudhSxz$Q^``h3GpiW96y2A+{MsQK|-JB`TadMwL3cDz;s+=zU8Yk3|&mKMbL>vd4 zKlomn@0uR|5F5Sx=$Y4>*MzRw;3xYjOFXf@@`Hml1CHV*xfOr)6RR+IyKA}uFOlAr ziVTZB=!)hk@}rAt#vR?8QD`v~RY!}kJJIxQbE;~3@?Q5@*uma{LH*O*)%fdgXJ+HC z`_DPB0==v)v6-d`x66zf7>5{6t>%y20dc7SfcSXCMX{njKh(=8ddx=PK65g|W+Yti z^r}Fi4rwUjB7}s+0E@f5QII-wU_&u?YKVRb=+K^S^^(daR}dDoTO|@lsLR-rfpxD^ zi7jgbA6u)Ec~tGE!B(%;V>0Io=EcyN={buj$#?v7Y8KDcoe&Fl$#X%?pMbx6-aAPa z%z4J+r6!QVFtmAc|C4BH;QImI2%u_C7DsbSwE+GJ^3)Yb%<=97TcsbsKX}k_~!PSg9<=H zO^@#(lw~N)bY)6D8+|;gl!O_0LQm>?K~%^hx8fZ-RCL9b9?^SoY2Ddk5F7Fbd_=R$ z)7sHecht8#Mo;H+WKhrKH}RwAJ=GL=bf@U*P{MT!Ej^ccZfEZ3$w`#+?GJRG8-4LR zHqqpzOjF>WkZf#fz~=b}mYPFg7qv;IbK2&OMW!8g;#y-BLUR?DmzE<6fPoZu3@#X0 zua1aI#}e}Q1Y|DMcj*799wnvLkvXZd{k5`Pggb1j3hqOkjnq>X>Mae~4&|br#o;@y zrMz>zL$sX*psc&7ZPgiwgz3^bu&kgBpB|Uw+T0+v98Ir76%AU7p!L>))*}cp#4hvW z&NZ*Tvs*zDYn)3;e7vOL+G42qWq$RBg${ROe!lj zPWVkm@qiPTCMpDwOki{o;t+Spn$F_kug701z{z(W6${QY_JW-opb0vOI$fGS@{NUq zRs$dcubLJ2&0<|GoqEJbK|}#r+}J}))7k^wn)HDjerst|cq(8oapz-2NiqHK-JOig z&Uop%E#GHg?JTak{(_2QI85NJjm`WNove|BPs@L z6j%ck*!+CBQPjADU?5UEy;R6hRro9VJIxk7X*fy(l=i1a$-AJO22t$WF@bVJIAaHo zy7twGBC;JqBR2BNvim8M73<#cs@b>~`QFKCTI`*9sV=dd8%6@hdXbXB4EXIMuMlNg z90i?3JXWK)jm|`ZD(rgPu^{kfS^TW_Xc5zsLH}quTitmEU@X1YUj$H(s z&Wj-F<$IYEoi(=7GBeGQlznQ3L8ne50?xP*!?HKpBd`kkhBd}woOPZ@*{HJF!~Jo~ zrQZW#TER=I+!3hIhc{UgD)7Eufj<(wx~T$?i?4lXeo*rG?zzi|wp=z4-hebssEVB> zR+C%&C)A%(2$sYgS=f{Q_HItdftew z$bdY1TyyBn#&sI0-k(exJs_&K34CHk-geDi@H1ooKD?^E+`6iR!a}gsdRx;T8WaY_+i!NvH*eLISq&ac>Ue8@J|LL}(tUFBHWT=gy68SB-Tj*oyArB}A;+orlyoO|^ zLulv?JZ>nHp_M5`aB@+ltO;BT6}zlLL6{YkIfF2Hl3O6(apnsOBGhy7x==sxqp=N- zSb3~h8x%lH^9hT8Vg>s-fBTE&m=Tya(+1JAqzSLmBq*-y*h%;IK*J6(r~-)gvz!U{ zwTV&MTfo*wbg3wL#uL24K6}gG3V2oZRsKPnP;O`N03W7ZNw-|iDo`(zT}*k_kQ+`) zV3Mwd#DuiMJlRfj|=XJb)SDD3W*1^JR24gh!X5#6X>R{CI*Q+a71)T1N!2_ z1Zu>6v`xcfuFY|UKZ+`Na+yNRyn*L^A!Iq>lV<}UWS2!e`vs^&Kb3>+*O(RT_HeJg zwC)Annh45yb64}KrgE~qvOS1t8K^hO{x(^nh8=UVkhJ)Ldtk0hYn$kHm)nPkSONZ~ z+0m21LD#lI z0riFFy(DJyGN}Vbx-I8y=P@voRIDg2nFuTWp{V1~9*JWHJI(0@Vv=ga$RA8NIK8)T zP2We6xAwi=LAKQ$4eIP#!3a|%r_fzRx zLyb2^WpA$vENsPc(gJZ{`NmRQ<-la=oWGMjHtw12D0sfAM}oiHW=mPPFxjoQEGbR5gQ;&Ky^8>*krkagB|BsuK36Up zy0$o{+`jBVh{wy5tJ|XBNQi5n^CIx&iSzT%Mw>z=aKWo~-wM_jBKA6NZQ;}P%pI64 zBnb#hX6Wb~+L#*)K4;Zt!l|SpIF`wZrl&_ST5FX#LZB79wyRn4kd8VNO`9~EtUv9*MPcOZxYKAQV_d?P7ZBGYRDoiVQNOxz-Xh4fER3_ z803@=EE+<=Fm@bZsNDefW+T_-9KtMbexB^DWsqv5N4if(spN)-=ck+abPY{4)e`)2 zr9woI5eVEv*Gq$3tyw=VIoxlc53$hY9Y!BLEX?*(Gvw6R{nl=3R8(NOQ7URF@slgV z!+r2E>sp*5q3k|+izwfkNSVzV+14|1ydGGnLdYky>08w&DW4I~eV`#4H`;PYhye>7>=|hWuDnXVaAVIzo%uh}B$d!r-J?SW(J! zh@ZKiOK`qNniKh8<+ZP#ORm@0J518xq1;%+R~uuDlSl%TTQ{qDQ^) z7_D@u#y|G`({h z1Dvp3BkRMCy72p`hian0w8|&ho#ze-ZTF}?qM=nRCWOmYZQSU&6M}9OhQx0@Zu~A! zzCb({@x*TSS9xRY#_AWn?jfmVI7CDqh2V>d_~vK?w0i{eTiNX6(1Mlve_Yq{@|pyw zN~pu1?#AXbgn=xKX;bYbc8K=@ghjM4CKcx0E0!Ld5=-$g%0=kmeQ$u*IhQZ8FO%c) zulcin4nbu60d+@DyGoREy@)3Z9gGjgkn6g3i%_W~4EQiAIg1_{?(MIq);&dZNNreI zH_m_UzUqn4GSKXWEoN>Uxu{wML&+7(LFUbS)$#8Q3sIX*6bGI;6#&PEhib4UYf}-< z+E{wW(NN_+NcEkyk7IGXt5hrycx?#@j>(BbHU8Kl*fI)lLEVWM;xHj?`z<{NX2T_? zZ5)Y#gIy(L-5g}ZaI2Gl+!Dm?)x4F+g)Q?BsV;_$5Ymin#)Lz_=7LvkfDxtQkR|2g2 zrf|bYvlf1>He3plcXmEuRt9>i5QCr8PyB)qT)1`m#&{m6Yieo9Ww;uMKFQDX z)`-Xc2=2E-jPiqSqWl=|9E0ef;Z>bwE6^42tDAVS$C?b&&gMv}N|-EznrV0S2j?GE zHusdb#xbhfDQ_2pJt{nj)#bq+^Sx$uz1Sm{gOs`=?V6Zsb8}rnP{nKq6I>q-a`T-f z5ju81|1K&PT8l3@;Nt|MmUud)4()J!w3wFrLwY1ARX?{Dw8T;5Nxto^y`eddhpKXg z9y0DmlopV2`sRgt33tyH7g(>Em4!6*bb9&LeeL?&xl2)j*MAeQU!VK|i*mO&`r|K0#KRRfp3+R#LXyu5~Rr*l}k);Fy zHk_7wiQK17pRP+;I$NuumK5`1!9&(jH7&y&jzH-7?;?6Dm?ja;M$QjpT>4dY*LI|4 zA4;pKyl*&(yYb1`{C>!|+Tds})Cx)qqTjTM`%^GHhi)8&H<~a=mVvb&t^LtuP1aTYsifpDt;0`knaWslLWB-lfFc*HyCk_(EhCKA zSciLDI0~N5{rBR*hiMd6J?nc$i78ax1%$_qhIq4>5X7O%5ImLOa|di|4kGVT)Bwd1 zFEqd*wo5u3zpY8BqmL98*eu4~c%|WLVCUPjsQpJ&M93r1x?GNdPbMIJ*omoyQ@ZB$Tbif5S5g=xO3D<|CVVBM;1Nqt4-WWeFU3oVxImv`C)0Ez= zeIgg1AItt^W$!=P^TzJeOHfwF4tmI(x_7u5S0Nc;m|JuMO^ah}wno-j)3$WL)zSt8 z;9_%Qz7VmD3JmTL?04dWs3%;Lkm*2A{SjU&L?X=)jo(!ueplIs=mQdd@1Ftxz1=;3QI90 zq5siP^YFFRK8(40KlJ59bP&cqg*3kG5(%)s=@_)o{H;`|hAH%F-gS{JI^;}hZ+-ch zNE&h2PrK$eRi7+oirrU@%?(P?AK{9OshaqcPIN|kI3klp$MA}ik{Pt1f%yta3Fw=p zzbe+zRq3g3pZ+5+RH{lD7W&zP$Y2fI(L#A(5`=rxuPB1H+ETi>>Pk zPfE_KCN$X>%vr-~g19?|L!rrBI7YlaN}j$+{85TZwZwQFb`>qQlu5oVMH$<7CG5ky zr{*xeCGC@=bfqZU&AbbvoQQ4>35A(z>ahlI>iptv=^4z&JCw=J2V?}IjB*6#tyMa~ z*OYvQM|R_(4(@$%J?O>m)u2BbNCEKS{x^ebZWD?7rG&*R`H_pD=jDq>lxT#_p^AZb z{Hj4fDP?9_zDG$nH)X_{M;U>mg$)U%-~=2nR*(1^^uh)??}i^8@j!*bIw@M(S-E5* zkh<@6nOM_r>MFHn)p&>+c7GBT{sQWiy+bbTGS!0x_xDxRt5RpWc3G^zT9^1NzdMDu zp>$|aThxA4ApOgz*8xvlW_RxP7HEQA)cY@XZrzkbwfXj zaUZyfelD?Ldx#v0nA&Bj#hwu8#Zb*eJI7$2`Be*JZ|4WQQ#3}F%Q-`D&^1n)SqnSG z4K=4)8s7cLvaWLY@x?o>PU!aPN0|KBFVnb;Hq0&@(=`|9aSxVr`Ld>-alox1+7rFF z$>Aq=(W=W8Q|)c?{JS(5Rq2Gv{x*ACujopF$El^8$;gqoj~st*6Xuv`&NPH}HK&?z(@S@FjH12P+Ye2n&ZiV9l0I`; zZ8(%MtNK}eCyDpSi7mT7z>)wiMs%iJsw_e|Li4mn`)uT3p>3;i zNA%i+_yk;^zd)(U-Oq(n?H+JTz9LB zMc)u6kcd)cTMwrx!#QsY^Y`6aY0|H-E54s^l%g$-Jl2b2RHds8haIJY;!RP1ZjD(8 z)zR@zihdEPU;?bKhZ2|Aa@QE7hUgL+YxELX?|C1*}*gKn#;416u}JdsF_#IspBw~_p;h!RNEz# zkgHb`|1*mr3dRJJOoV%Y2mkYh{ z79h32#y#_di1@*@^29D>?1JjR(&&#wJ@#i~ui`os+*+u2kic!g2KoTjtE}5uBL#vt z6~_D5*1jF57Im~IU64?@-IoPO*+Bw{80CQSkEP;B!Y{bI)%w!Zgm}mN6=2{iYbZ?z zqkFEkXG>tanp7s8ev%Z2R(5;WVj+r)Eap{lDUI}Z8RGylp&bWG@<+SO`ECpvcNLahGc8yk z#gK?uO$~P+Iox4atyhci=hUnla*y*K7KZLRcviPE^z`1@pNFFQlM@oux-yE?$^sC6 z=*h)U{jz)L7wL2v#=nVktbaI!4*FnJ?~AUTw-xtt!UH2dW6(Aq=F~&1z3X$vcEhVF zYF1^3L(?%~*^=^PpC-+xG($GpQOWX?HdQsT_UPIgQ>h4)b6ef z^;;lPLESg8+7&C~t0J+Qt(4fx2_yi>RQ*9Me1p;Ix$GxS#2;p0$AX=r&&RXFL!9SJ z0ZLE7RonpU0Z@2_GlkX$A#s@zYMh?FL`{Bgr*3VyzaBuxFpr#f!^Ia8!0Pw-%`JJR zLy{>q2VwyXH^~TL2^cmArgMR?&dpRAFsNh~9v`xJSMufB!t!|G_tZ< z57zM~;E}1x6C+z3KlxXO(>X80M>fu=LVRFZtW|fkfBHs%;+T{8NPTc6vdyNxTzRLO zTWZ8*`TnK+wWLc0BDJk2TP-VMH(bh|m&xpY+(Dc2f|ZYC@8n2IaOGu?!ZBQgp3$B~ zdXLify)x$;Hw#XeA9;#I+#N08UC_I8G$;5-4D7^KLaev(w6ZTmT1qG2^52BHkjL@c z@f0JFIlX23``UO_kMMA+YM~mK{z%1Za{kF}>_x&Dv^%{KRCc&4yZR*!cYwis`4+;{ zynEi_L04_$gNX~3T}$`~+$qOqW>${&DtODYjA@=I2@i|yQGC?6#msX_8S9l?&#<$byrufwe{D9i=`7D zRw8mbXHjLsDdJ}fI-H?XD(4jvaSTq|laFyv-Tmxz{QTS<+7jUJ-in0bVBy!FbKf-) zZP1VWj7v$TB$;V3+`8TF9tr}I3n8IX4bH9YPfmlW^FDmMD;s# zm&y*8Pxf3|R<%={4-1kVidU8(r~Jz?;00#f>l?fIKD-b@e}xi(uNS4W5dfOTsKL3G zI&{;#Vu5lSu)Z+T#T^BGUZD=1N2~50Ouulmuq1+fU45%Te{`AV;(|uD_O&#@!^c3> zc@>oJQB!*hxRAJ5?Jg}47-(6EwhuX}g+vQcwH0W`a8kZ|G`uw&6u+u#LntuzC?q!2 zg~jd3AmoWFg>h#*8z*LGy@N3><6*JDh5|SmJ~75j)w#RSQw+8a4dGF zwBf8dZ^RfdU2#SD`Ie*_4B9^GnVEDh-u}IkaU=hItS0ZbOpHPNwH|?<)m)m(fAeYJ zRD9x?ARFlP-x}`25eqc+6aY$-v8YB>ld5!3FT=kWL%dF%17x~L+p=+ZcbNQ#Sk7|WJmFZ?^<3}H~YlLQ&CFxTE0;1^h(1%%RexDZ`x*{@gg_^%Owko(9sI3=~#QPav9tGqb;+M z-xf+x1@mw8_e6m3R5f*@*(hBXceRO&kAm096A#~huqHpaOntboKJCgyz>%b`vbA0{ z4BeD@3ptuOnafU#^>QWqwkR$nW(j|G0lSjy5Zk9l9mR zRn7xZEIo1+9HnAL+Dy3~+r%@+>NzpDFc!?KZ=Gx2YMo;Pq`&C}@iKUNv5H_r3= zoXPG??#zrO7bVvsler|lwYRaQvt5s^RNq_pxHx9Dr}=J|ovUh~$#<*@NATLyZ?1h! zc$1JX!;{H1L5ksCYgfg^XIPJ~c@)~CPrU<8J*l}Qq8qY=hkYs{VT4m){#T!u@dckl z%suH{G4-yak&`rR3*Aq1SH-BiZbLsfyOhjYc7i|6sD+aHbZzfC#z+Uc(E zNC6g8&4?Jq5cWMb`9(;3h@GE-X3xr`M+nXkX3#f0g8c*MX-#HhjJ!khRF{ulZr*xt za-9fz>a64a*+q_wE;?*Et7wwYqF%**bymsI3@IePKWt>uFdv zN4S()fAcdgm#sSsw)YbVIyOb5Fwz~@0^vnn^&4nw zK&Q^MfjUj4>(T~F^1e;k;z33g?V2yh*7=H?J&bzdnr=Ju&R2M)FtG1P)bPv)zuU0y z>NeG7K-^)L$P#xG<`+Xouw`-BxIXnyul;dxp8@L0yn;7K_+xL~8PZ{b{H#F9O)g}l z53jPSM53U))oYs~9)+z9pnEMM&tgcUOyHf_)Qrpc#hVV?YDG$MaO#kU_(4}>mIKT` z3ZSPGqB(cB;MEfGovkF;T$RcmK;0YWwT0fBr*t7^*hS=ois~W4FpaC+_$D* zTyqxl`^tzAPW~f2V<32v7IDm+tU0L6r`1K}UDoaDtH>*rZoXnp2n!Rw%YIYIJ?LpZ zgm}%)UXNEZh+LYn*J+CiRYSy+zhy+2`@D_?GkNlpouF*8 zq*#bQYEhj|tiO(SXrdR_#FjR}Et<;qAJgH;qKKrbC^(mW2qG1k0_(GdbuivxL;z$j zWK4{(mGZ#llFUf3E@c8tW8#W=^ow5F3pNty?e$r_t(&L8)x$>0z8Ph|crJm8gQ13H z+BOvFf^u* zOvp_S-3M8SN=e0=FXda~QhdNUWKi|&_L1N%EC5e%Nb1_n#9;IY{T6d>tG=1<+FKYY zHE+&h^tI1zCSrGk<7>OZuq1$gyy$}sWIxtAR0F&d?<0U@coA9*^6q&bj~NSdZ&ZQ( zT^&Au_5EY?5q!Bvl4B(EDtMBoUg2tJr8o^uXTY20e5c-;SNkrX`BqF;R>Nihvj^7< zU9w%xCyzvo>iHL{y*eO<(OT|jfi~hJwk&fpuCaE-BWN{k;xIR{Fiy9grJ6kF!|*Sz z216Y!4U!Eja<lM@mA>`x^!00KpN*8xD_VK+68|Ya zu=J`(lZs@l<}bieLC6$HGQdE@2S>>W@zN62@?z<4-uMV()ve!7q!6!$ zjHzc~v<>axlNtsNDt0)Bumhv^q{4VUr4Sz>z#AD01Spw!U|kl)FpYv_|&=YtW2vJy^-RQ?pX+6<1{x# z4Byjv!HZ%(A3|A~(|rdYB@YkG|6X$Py>+X^f@}ZMc6P`TbJ1uB25Zw}!II$sFgT$h zvoIIJh})Ad5RCS-q(r)^{2)?F6c^o#*azKoh?k#CfEW5zIB4FubrgQCS%zf#fj9!v-m~QEEwtU7NF>0x2!rvq*k?oZK`I$gtaW~Lw-E|E+Jq*Eh zq| zCoan$iz-iAhU8!BeOyr#M5b(LxqDXt_z-iWwTn1l>2eytFQf(YeL3rRlr`u4qe zatD>58x|^eSZ9IFGa90MJ|gXU;h`){rXt@nWS8l*VRUO^)e?-l(I3NNtvHrli0Bkv zVmYIgjVuD!R+(P}k|5xkm!@a8h8Y1U`2nSAwk~KP*;RPhJ|dORcM6u0H(%wskk;k- zVh`)RRhP5);+`sI{1)>eR3+R@e?LXeONcL%-YEk!T+Mk|Oio7$XCVW$wL`W6naX{? z#1!hpOWUQ(lNDwNTH{yZS)lNRU&8~D$nZopV65M{AV@o8>wtbgV#zJxh;2VR1sUEd zGsH+2bVwC!l>S-qb_~aG%#0Ox1$rj*8k?@n&OLgmXzg70A@eMiF=)z1{~m5{qW+~O zj@Y?%-9x^Y5!_zZC0a#AMguBu{Mv8E$p_tV#K3Lhd>eoW#iXIM|T ziiq?E!Pk7vHO;e;B35qpEhK44zae>CPev##$OzK)m3n-89GBtf{cRU{F@CEituQ1x z#CpcCLQgvW1TzCMa>yqX_>XKpB?xpB9$?e#W8sDqMhQeN@uSAgft2h3svPn&;!#EZ zDulV$fh7Cu0bO;~$*!xSF5M)^zN?K@9{yn?SW`B6Bh#AodG=6F%l$Mkqb4bI(@az= zo9d9xgAB=S&N}I~^$~APTEOiJWMFpEAs2H4zhgK*Gu80avzS+3c&Ax{Fct}6&fM2p zlc1`{1u_#8(X)a*H-}?MsGWt*UEn8ym+&le?>Z42Jk^*gj$3t$(IEI01{ZBUJ!SC%~adfne+!|gk+3+aABz*pOLBt6P?=wEfm%cM&jH_u3Ngd)Hh^LdMgc@Wj zUS!DLl1)(vZF~V@UVaE-$>GT$lu^6SYd1Ov^zMC_Y^p~^t1!NDrucPWvwHCW5pA!x zTaYuGY8>pfshrLkw>14)%=EaB5a&+XNbg#}Sl2qLcLymywX_?aQ^VOf$ZICEGl3Pg z^HqYLS(2KzPPU{WU0gGU5E2D%fk8a~B6N!J=aMi;SDoM2Wb!r<1nN{mI?)p?Fqm{3 zIWHVFSeT2`b_SpTmTE8J3)Xmd>InQNxwHBXNtP(`xXPYvs@`2(7ZIX-`SUK8!9kd& zM}3T9YsXZ98$Y%~#L^ABh9_V0k?Yr1#>6GmVZkE;K2-LmWs5AE67hW==eOVJ2B{9% zhTE`Iu`#r$xXr5v@8C`5$x8%AXzTU3zB@BpyH$y>q{I;^$RUm@7nKok3*4%$4H>(n z0)QypxS+z=^K@#If5}d=IX7kHfS6fyK@M-V@=L}mM7^jU>03cTgJwN}2Wz6^WlJqn zCDzY`hsmsCp0Zr(Gen`Hp2*oW(GXb>H$hhiuKIl~zQ#ccW)QR48DAjR zLruF?{uH4OVhF0wAG)iUIuIbmGd%0hEyhbx_`{^DuXRdCnIeIL%HqU7d*20UicB7UfPf*>Z;aQ3Qu#eA>12bo=5D4 z5w~R3r0t_u<}jq?;(1{D6Ms3@4Q3w6$j?{deCw&j%a8GH+5Qg}B0WN8F8K~u8^6;? zkhe&|rCo)&$<=x@qpo%ipZ6Sco6g_eS~X|^HtJ(g+Sy8Hn0t%gPI+`T{5}>#3?d7x zo@9x#h6>j1>tPKdBJc+hDe}W9rlx#lO7zt8^4ZS;mTFdB0Aq$O#E^hzSog@bwLonA zaQB!vzPJE?YOe$4-MfV&l*^{6I7(qWJ2gvGL*9x(j6sm_cTJ@W{-U9-Gg7IPf zHSDSN6lAXUz-&q`uW4^D^^2ql~r#NU7 zpJo*Xk*yY*EtF?Z*4Jxig1K~bS+wwe6>A48dlDzepCt#iBQ2UAA>I?y7vuJvEu|l> zxNjd%0Za)kC=jjs2G$XO^84rI-;*tk9n#g{R?p~OfM`l~xxLZj#xI9LsZV{1^ht@f zL`fvnOlA?YZ#wL2ig^2R-_J*kQt>(_actic*9=opqsFYRRIF2B(~gns73ORC2-+cn zu-{r|Ebhq+@~2ufq|;S+&j|ler7KYElOfl7x`kc*vOLx7A2IX){;&V|w;!81tw)_a z(1TWX==qDvF`+*Y82IJni+TVhFyI^pehA(kS*S^3nmbHP|10UaEc79@m)dT7o<@Z! zXuAQV5>4T^;R4b(oTC+J69%OKuQJh22uxBE7Q=NnM}$z9B`G@4mv!Oc;Wo8QtYsfF zv?uazKV(6%Ar$1CzMyb!os8Y#NWPxFEy@IwzJRPiv>QH;QqDPFnWVCykh2wpHFAMU zQ!ZCvbPxCdV`~La7J<0Se~xNTGOoN6vNjIP;VW#xid!2C{G{aP801){`oxEjoVyB! z@9u}c#A*b;Kl(3E(2zE0z2R9x+iM|5Eq9rOU9V)qBUqLW(m(vyMV0JBnU4h$x+ST!s31l(Bn-KGhC*7wWEj1j$dg zugnZ>WwY5YZ)=z!*;>k(!3sbeQUMsGRo>!STlZlVv=dw7mf{&E+!ccP0-uQiVzQBI z&{g1DuRN8Eazbt8SDb!|D*4hNzs%u<9gmFf|S(RZ$lb{?ok zh2M!(_;P=7$Bz)n?g|5sO=Z&%5bL#ZoClk$P{^?9(UR#-Xi1qmuI$DpJO`MKgjN7? zR9d&kIE`wr0HT+h6nrL)ZQ;+Yu(frgsz(CdUBi31D+ojWbA>53vF>i{ZlW+<6+6_f z77oNM*V4M&^#Frz$W3Dpd*GdvY`*NP& zS}CDXJs9qH2N2wecK9u>m_}2?USXG4nD^?@Oi)N}Hxz8WUJH*hTw{BnL9{Lp1JE4~Asp-1)E^MI{=bPG})%3KT56u_CiF2~~l?m}5Xm`pemw3KUk1tGuPq zm)f)II179Q{6tIA#&F1*2n;`+%MH-R!pn>)T5NY??x4rG@mMCUYR&l;hozEk;W8I=qA%t&`*5uQ+M5*VgQj3rQyF zz4}`13c7N`EsD)`PVbLU#FSy5q@!R2U)I?IO@b7$Gyx>{;5?R*%jLFu?Wx9N3Q3I? ze#eln;0&{)^l19_CP|p|&7@50<)rkri9sjin)P6HQglF$iOxlR5IeHJZdYv0e!w#*kEkNvb{S63vW+hb0go1nYqb z>X1i$v|0hXmgS*gk`YLJg4!@aI$%%+L)ub2w^TAEh?Zc)7`nK&)>^F|Yki%bv(8#~ z_a8_yd-neA?{80%&-UXKnJyv|H9Z_onb+hfO9r!zoG$}}EjJR%YsDi>c*FdY-7ZW{ zuB(T9 zwIgiNf)bu1&%maHrN9M`0g^YeTq4tC73G4Z4XKp0*AURe2E46DedrJ-7`M{EYO`y* zX6BWmp*jFh!m!Vbdi3#x$S#CA4lNlu5549Ur^a`lgN_cL^ro3-Y5ll+3cNPJkp;95 z@Oh8}vO_>B3&R8lN=0&PdN~Fbs6=WhjBy8ndTFD2t&pqDY%cwv2(M%H!S3eUExvO` zHCmT5HLnOPQzqWQHuHgWOGjeg3Xc^@4Y`maoua~~K#b3F zC|8{UrZo~iE@mj7Nh~0fE%_F$Ps;|EqB>eGYEi)WLpGPP`h*a`XQZP1| z`g)1jcZ?VEalJZ|!lKdk5qd94 zl-dB$MM$nODY^2}yaMggo@3kqS`TxH44fd98HF)yD~q5#TCxx3O3@%phI+W;c~}6} z0GvKWAK!(^V0GvMU@c~N3kgpK2BtOGIOu*zVTKa4WrUrLMU#z~z`xX`1NFRadm$jU zLj=HQsc2>Z2)uZ*nHTbV^wD@3PiYIMkOb+D$~Z6po`JSZ0Y%bAC;3$}D3wctt+GzU zxVt4Xf&+UDDuP|shb>o@k#bTc!zNXUm2&EQMHrC-oEggH{`L8=K8;$Z#li^bD^bn2 ztFNW>bK|M-F+JV?I+i5k5NJpZ(J*48?L2q5jYWy4IM597P0muT3?J8m>nmhi`D!#o z6TOha>I8i6p`5&`*^}oOCQr1>ED-ij)pcUu|Db z_iZ>eQ|n)@j8&L2BSa?x;By&((sN_F;>oDum{XXp<@s9BFr-Q@^ID7=AK%C&EJ?5^ z9Tc`SuEN@`OsBgu8Bc6|P~WdMh-I(?m6jGr0CnpChP90eENbDzWcW-1NoGGgk3|l( zi^r|Jc))emnX;i;uw*JG*H7B|?HceP3_ZnFLdj)O#tzp~kp()QPcTnE2EgLjP+utr zZzjk<{YXB*<{fQXz2D}STJ0|BhqV-M2>jPTu%(U)s`F^1ovDO z$kX}-6WS7$$_48sS)ra;7BkDMza&YnR1p+T;*mxrD|Vx$C{3h_%4Zu=5!wT2T!QEU zS~fSpRmosmWO5EVgmkIW4GJR?mFpAN>qAC3UzCh7{#=%x9IBgVnmrbEXkci_pGeq^ zDPX6=dd)%;h^$CN=HTF}aRU>85kNT8fRc1r^pQxnq^?|*tc&%*#5^rsOLl;ZQS#ZG z6py+HB#lxaXR`IGF_pLvXzfrU8PKRrHZUut7PApK$jDN-0fZsQ_#ALZfz*bakthzD z0w(}>D?AlGi`L54>dK~Kl_^-WM^D*#ld}m7XX-5Ius&JS7aNL!Y6-u9$=2x|S1UrA zZix$AY1ptn(t9!)B7^-Y4oiuEig-k+oLa9%&qSu+E`i4i$as%d+iMU?RG|!N3g$&7 z#UP#0K7EfJ1Eyu2#UwKtVJxBICZ7*i&WR=109PmH7-9Y72yX!|3Q`##1yU>7Eyf*I zPDM0TL;%FS2L~Ys88;?zm3A~j>qb4WO>xp=0YD~78{}%t=0z14O-^yaTIAAzekE$% z3=`^>X<4NZ8I_vSuCaO>fX5CKqKMao%SlN#)*jDz$^t$!S6W19ZBZtS=6WiXMzu*b z1oQJ^BTSk}f!&0WtoQ)g5%3voAXCiOd-L=@y{{;T8zA~Nu=Sv~QF$!2Y%`RNgbYd# z7MWuk4k1!RRvxz`FwB=S0>R;3EswT&}DH-*X67Hv@?{JQBc6@CW3 zF~PfDrwfWB0I9bxP#9GK@RNz!fGrtlYy#+7+J*c3!Xlp>msfWnF?ZN9=3;BykXW3`NoZBj8| zuqXSFNDm8uscl>)61gosNlLZ72+T6 z8wo0(wIBNr2E&le0Vm zD?;!>z2#yWX9$rZ&zRAN*;zg_Pwdl!NPM1MDVZNODB>KVUhHOU02%7=1CACE2O)7p z371Rc`IYgr{L95M+q~Oi#$k!|S@94iOl21l*c9RnYoTV645N0w5xCt<5OLR%0LA7& z^|gUmJVe7-6d$5(Q89HJls;)$EwJt90Z$f4%t*5rm0OT3Y1w3sRzH-4FgXU~4Ypi`fB&HH8ycW1D66xfHTYMe)@F_e* z%+SJ6xen~{A#3c^eEtS2kCG7xSffVWD8tl7fZQrksVIj@70NRs0l?PcwH<^(megyt z1N{~a>+%Y57uIWq9s0-ut&ziWf!NqC2U(ur*rJL|LvE75YKJz$T$&s3SaB+Bz!PyM zt-bqGU8=Y!59DaUFo%<8Szl3SP!g!jV*&YVg54^`z%~2lz;+`{qh!cywcS>5nW)ql z>(NsE;H_1ECAFI(u-!C+;r2Uo*{GAS^>8w1(0g=&h2%lqp02RbwkCR>RGwi~S)m~^TF zZpOvuusQQ|+I~A5AqlURib}+ja$E(vJQvFcNlgkzH@Ndmz-^060neV1fxHDq8@2P~ zYzyo49stvOk+4PMHFH>80R_e~tQxK71X(ez< zN!r1&5J;Ic;soY#x426IZiH1N_JUHVA~Lihajm0BGGa06yCo_QW}8$ZG%5yB&!b;1 zs4MeevRXpn(bme1n8#>dj<;u&bRBCzd)qml3{n`rs-*$_t;|zfD#h*)BKE7=RZ?~- zTg?@W8cX9W$7*FP!e>(CL1P4*!8U^Zz!(DoMwt}Q6yrGN2AN2unk`+g_vC}fWnFK` z0XwKWLJ~1lgbbk`t`cROax5w;gaUkRzcWM%q*p6!)D+J)=8ihBT<|gxDiYQMg&P(7 z@(hYo3Ij9A5D>{QQ|Xtat+o)WyUt=MbhIFo1j{7J^@095KutlUB%lG5H7|q&nkSWf zWSj#^XGu0&nQDp;$Tv^`%mz!whb^L}5FhZZy(YP6+@i+=uu+XL5hV?dDfFPguyk0U zD3t@@8bwjAk(;VEfkU}ogh`}f?nJ54AT#D+ekDM>kq}7e?p~A{qIfjNEH%~=-gc3NSqMCIcWltju$ zRZ5Kur3ot9XPQm2!{!t(Hl^(s5q#*1!h-aIXJWQVF^w^lQ{|LuL2}9jk=`H!!33_b zqd|3BUYgYbx1|W9DAjmlZG>H2A^~%mV6)X) zzsi%30HR(2&589v4m9Fzar+Z-Wt~A4wzP1lb~ru|ADGdv#H4y3hRbX;NHPwS(qSpf zVuWO)!lcw#jd}Ad$3ZossHhayYN@Hb<_yZ8tXGl|JZObfcHT|4Ny!FwyS6U_>y~1E z#?C=~ntp(}`&uBJu;?&w3wmw`m zn-QqoQ$h>48>49kNusC1FQz7n94!?L$Q6PFG^F?N*)3WEVU(b8bVkKD)!e!~oy80Z zqY6X1X?yQciy>sD4oNE*x^0MUz_0Xs@J6Z&j^MyC#9{yNwL2LvhLC_?Fp#Zy_V z{G#lbIyG+rthNH4T{Q;=lzfJ-MaX6FpiXykDTvlB4Ms?g_VZa@3#(hMtBjY^B??1r zuu=hf%c1U^O7Wy41scYf0yqhi8f*%Kst%MDT>~zjVFRtuMjq9@G$OTvbHj+Y)nSq$ z<3hw8fT4~o1)@@zR2{umP7(@2^a9wYLjfBvAh@8-D5Ja;-~mIr&L9qJgvdmZGY6T& zhanY#IGoxHT#gT+Q9+hCtnYzn!dAzELM19$nnak*2P-hqS;yxF^vvu4smje875MZM zhTMd1Njf$IjHg5*Nq!w1$uf8+29_Q;Milhw&^lfnS7Bt4&HgMV(S384_Gm{m&2p$0 zVnWnpT9F+|bjB~22xmv`|SHF62-hOJq#`@!CYsoI=1l0K{I}?j zEo<&3y`%oZ^R8Z~t**+s{&TS<8^1sIy#z}kqxsoie5~ zhxxOA9+_)6ToZj-GJHz@@J46QvwwQ<+xFz(^|^w3H+HXQl|SfgvR-aDjjk~aMdnd=sGrtwBePSG|`eb@y z_g($60D5uFf6gXcy~sVaAtPIVdt!UBsB_!O=2!j|mhu?TMyK>X{zd!kPZJl#{5R(O z__yi356?cGy!4r|2+G;F?Anj7F5hV?e!~>MKKf`iw{Y{}(9FhM{ zj+G}1JH~r1kNj!ex~r)4`cr!5gFBBZ@8$opZEwyFzV4fVBJJOo$!V#v`s&@tmA1IV zhg&46CvH_f%#Bg3xi{^C_AOc1_CfvDD@!j=?})1@-_!W(;T4yR)5+D(ekn?6_`C3$ z-rTW}T6yMi{Er9ju5qr&uGyEf+;Do`m&ca>+Tfag8M-^&I5cnb!TXg*!p7#R`^LpN z&h8!6KTSvHpFeS2AW55>$iBU4^Q@7(v2Sk1$2|1mJ%Uz;;p>EnoHs2=+gpuq%oCSv zN_MSY_{QS5`o2_LdHY0x@B{E2U29I}^K#RHKNn7Sr{wZA`9H37^?Y;z>D-XHqd}BX zJ>PkK+3NRdb|0x3%{;c}`Oytz^5_ju{6#b+Nxc=$xfuWUVp@s5JLibK;K?TYpt>*J z{ANkCV(d;wJ9XgE&55HY$mz%r!H=4C@Z%?0{~okwnYO2N-9q*x<9a-uXVBs|7quW{JJmU4J^ z^6rlc*7;j(d5udN?wxG^q})2OsJXk^S-E{+QS}1D2kN~O9p8Uk_0}63R?q+2SB?kQ z;p}1s`MX!U;f3Ol0hjpl%8e_7bF;Tds9Vx$@S`_W-BI;+=Yws(wsy8Ub>DA1erW1Y z+8>K58mpGp9?5GRFw$Aymv8CvZkL_ju=x#3@?i7X=YL+gfIME2IWt&ibgVn|^vR;% zx4hVWc<@H`mP@-ogI3v2J((9*vtPP6dHi7Ov%AOvy2Ye~P6plk7ZlFj8whD>LlsxbVW=o&<7hQvEOW%{~XODl{Q8R0Ly=skHir3%ttsDHR zt^LIz*5{k36$bfT)yv)IM(7K+?ATedjdPD4uA%5_dp|#C)j6Fckut|+y--zI^n7R5 zlC4LN8kWDhSkrMTo968?PS30vxzXCWDM=?A3YuaEU_d?0?Gb^prQ0>=E2NZsB`YtJ(=Q#tc+{zQgcr{~+}F>Y!l#%_Ye7 zWzMAq>P6>DstaM(vXO1$pBILH>j~~ncvgKovvAc|_b*2qgO4Zvfxfy}jSN)jQOWa- z$MFZQshx4BYd+9TNDJoN?|3#`eRp)Bop-BFb!>F!<{=6DKf~0(*J+31w|-DkEr3}? z&o#Ct^E)dGN|NL?T*;Cbt6%l}yu_V-ZqH8Z{fXPj!K5FHzB=H#KHFL+(!BgK+@Cu7 zt$0wFBu=dfrPqthV{P8NKr-k2ouyR4E$ZFL2zFjc3Y@gs@R_8aZ$Wtj#BDF9n}$U6U;O@PgJE zabT{h%MW`JeU(UkizdN2*d!QeETp?t11n{w#-}yYYn%2L7GG_+;cFV4uKBL+t)O9I z&=A6ZgS@(Y|2yX^a@X)X-Z?4JdvABOq&K;ioxC>wTWp8ru4mbb z*S5a+n%N}yPptXIw*2!=@!E0iTlFK~R!vE=xe_2Xmd^{0A0 zH#Ro2e^TL{J-Dm-qjk_PQxea`>OWq!I6L*^)>DO3tA2fN?djyj;d7dYKS&hy?_PRs z&CH?0*RKs`ep6RHeedJJhc3EjQQN{bhRvbFfz`}n{x^l^^oy^wo)#2orWUU9UwIg= zK6+*NSf6WQ!IFs!n;Kuny*2aKx9?pmhPQZgp5=%)KCTF!a5x2BSqX8{6`98a*KaMk z@!-^ZNP64u87Fft)TICMWXp#a-aomE!oO$c>nP~|Ng%CU)PXN z2>br!S86oZXn~druKw!9qpFYTH#)8!4P}j<+56jiUHtE>F7FGz8(#R?qAy*;|G0ha zmRh;z?B)AE>~|CvoVnNZzH0&N$S2#b>Z)8Dch5fcwN*W_GiLlP&tc^j+vaTrE6Nww z`6^1F806=(4+ok`_s$WviDIh`o#S(yzrMU(@a?vlYghNqJU^iN^8TUT)UzF3m+$YG gs3Dd010#R`;@R*26Z^mI9}5c_;$FXVJ9r>i&x{i{1MqJWwR zii!hIaYBsbM57TfRp%I!i}RQm2hztJ5l3#K8UFX`b>H^f`@jEnb$3Q1xC-V^x4B#OshdxXELo}JkAxoriz!|&|pdD;*mThDL2A{#h4P@YUXG{oIK<>kr z2|@VC&%`kU@&glt?q5#r_5TM_fBk=zI{R;bvWW(zAO6I42B)U`6W}s~KPLhJKKtO$ zqCu^}PqhXHwx<=Rg;U<2>_5pr;UBpDe@Om=ajN9Zsm1@S^7>Qc6l_mt=;oilQ@Z{1 ze+p+G0RD+@8vQA6^N)hVg$D)(4yWYLr}z|3tNsE1X!s}9{*U6Hef}pv{~iBR*Z*wr z?>#?N{7-=gi6&43;Mto7CI)9v#{V4mzc>AVqMzIUFFfm0`pm>X82{sC1_r0pPyh0b zN|vSnx9E?up5>*AM`_JjfB*YC>c6OE|K%0s?<=z2!I=L=mS*--=*@-2o@Sqg#Ke^H z1zTG$gJp-W6;3z3($wQ#X$#9?7G@@8WCX-2hWc7Z3-*|^Fk@fu_ zG6Mtq_W$%>Wt9Ku(N8!3N&5T;$FNr*q;&o}`#-qAZXiyZoUNjDYt-(~6TX_tvVh1HynI`rnL z=)#}!e?EV@d1n9IgEQ@Ei|s$}Ka8PLdc?I&CUpzViRnExKbdq!6Mv07N2sw)OfOIG zAv|n9^YD^Mok=GFMhI(8?`a}L+9svf0Xx?_ILl6vw*T@j(xchK?Bf4U4f-j2D7)o< ztEqTBXjbw1zr6kJ>$9@|Z2X^1v=_@`eNx-rGkXJ^sj~{p6db8k*i#Qp=mKE|3)0FH zcd2s4Zrt6@Yz~9I~+np zzo_x}Mj}r}t<5~~bvHu`ei3HBz%h6E;e6V10MaYMIoN3`2$1GP5~dMJ*23-h9b@wS zthJhKj~>e=xnL+YiCB*w(T?aGfjDWHH;TE#lP)ppeowS(&yNJ+< z7yP2`Jd9~bCWHqEARuN|?fm3P?U=^$hQ@Bd>jYW z4Fs?e$;m0*jfzs>`Pro0Zoni7Dv$TYb{!9`ywX&)%tt(=(J_;40hX_^t-TOH(zMdL zFAVMcx{IboDz}LRCk12GUIPU(gs+Vvq>$M<^Mpdvo$X(9{<%w?y~rO{-C*aiMApMd z%bwdMDJ8bb?b;@4M@fRB`61cWwmd^%hR4?;WbMjbX=V2%9PJC>_3#@{fDU<-AE-db zm>6l0Tj3t=w)I3Y9~d>62$j>*)spX6=B1>U@jdrRn8~wOyU)MwYpF}y3p@+!T7g{i z`igd)abwPuN$}5o-b|ZTC{tcK0X0uAF#H^@G$&AA0hK1e6%1n0~7Cvp`ULr8T@n~c_B_DL>i8@kzEY*OYVkUd75m2b#gKG*!cFSTB-S9M>ns==_(_cxaV?`o_g?Ej$() z^b$0I$50F5>SM1ZZtD#R{zdlcd7jsV^G1<{^Bm!E%|Nh6AI8ax$FzChdr1&w+e955 zN+B1vG25l>wT^OIulx@}|!Z{V;Aabc!XItj>Z$c1Xh%;T%#m>8ecM^O4wrE8L zmjwkaE_};&_N3_iEzu;RcpbyQfb$72z~jNt9Oa&OAjz{&yWo2#75-vy472d6o)w=G zo|W30mf$$Du1{Q%Bg30YWDHD(l^ZeBrskfl=Uw%aX#BHK58<2x;bYOq+>bdYPl{=O zn5@RsR!4%il?^t6jJKwT!NiHQnfIS8E+RTzodFRbpKN4Gg);#&ADDTeMiLn~GHpmQ zAE$U*1dL9^OmBpRfQQoeV(#5(cry%j_piW-sLeKv60K?f2rizzSd)( zqx$A{ltCE#;i0UA*zT%OXZ&2>;wajds;YiCcdV`97J^uSABx0Zjw*jWUd zK7Nzsv|JoZb?ne!Q;saW%OeC(VJ4rR3QHi24KqKNy)#<|rqFHntOfziMFRj;>lam~ zaZhL3i{a&}$l@fj4eBYW+#deWyIkGP{3>}YXP{SPuY!XM)a{+GpYorPJWG{hk#o6M zr&V3ta(6rBOsNXOI{cnnKajh49e#wwkjcq$dzA|3D)NPtdmfS@WBHS4&-;h^nd%<} z4}Q7?o`jaedc~NZ8|!DOjg0L(wc_UO@=g_NKlon@+*Mou9eD9q3@Cj^d5QXOR>hF9 zK+sHC)$C>B%V0SDrc#e<+ zrb*kQ@7NUZy-nwo;6wD)<2p_4oTbMUPZu!~A=qxstV78tmuGhXlC67_%L^;R-ZM${6Pw+6L9$o$-0JCj*@g6r?uZ1Dwkr*F`RUqytE-+I zDmo_UYtoeKhLi!~rq}1pN3k0$wE+o-%c+g3UO~OUZ&Qe9Ca29!n~x!Qub^LU*;i93 z=ocB_7CKXAQsM6tettQrhXQOWOzo}rA-lplHj)aLyfDv&t03yqn>VPcmgna4SEUm# zHjkwg$}>Rw4mEa=_RepFlLMoG?VHejN!}j{!a7BMP+v&Cypj#fM?{Ct2kco&oRtim z0vzAGuPBlJp1ZBxV83a$EBZZOQkoc(t@3cXj-WI<`p+9gW%7RHe#q7zbL-`yYJ+ybdI^Kgs7CueVl9Nd+L(q8XtE0aa>9E>-XV8vY9BNTKNNXQEX=#l# zewKJ`L?maYaoOO7rL+cA23V}*;3L}mT|%cMd)=hnZX zBpvp5^fvM>-~EPhfm}u4ByKlO6-325GWx5;N(&g6wvomA(wW_*YYaD^k5CUrJ56E0 z;VGqFl6_2oh5xqT&2#(`x^8W-XtWhvo!@msHCWy4MZJhvxLS=4;u$*%BC|BPZ|02` zzTW~|eMav|8X`R`xw6B|stNjmkjeT)0rH}a?)&Qu`uXKG7hY?_%7<3}KbS%5J0CEv zkn(cn?Q_~3S@I-#vu+}P$*1>Ph*vao-UrL6Ua?}PW+^^(gaDnpLM$VTW?e;>J^3m2 z>xWS1&Hj#x^td!}r4wi}ut6Iq>i{AO^VOZ>94S|JH?;}$`c1!3BW9V6PS+5Sr^GKq=UP-S4dnPSWtN^Ek{NvCe(W8?zNzyXtmgSS z&K|}$0Ghx}Mq(^&oajz>%%Awg*acs$E#L;+*6&?zbF_R-xBF`X4-}%;`0lkm_97{e zDW2JY&g*&p-kZ%PKmp!aYiTLCPt9yV%+|?s69Sdwts3U>mQLIZe4x_o1l@?;1H2w{b;C+I21+;{#-N}I}dk1km@4|weG-1PDOrc(-=`aIN z+HEkM=u7#kCZF`Cw{wMa>DGndKIkz~13#HMKfQb*pHUoe{Tz(GG1x>?&^BcbeDltw z%mmD81Of!vwVh+mPL8N$#>>vL4jE^*Nh{0?ALbEtYkOY#4=xs$4lmlRp3_%dQeqy@ z93S&igM}-xD7e#z>F2T3-tt9we^KL4a>2+b##(%Hyx`b-*{sD2vT7#Yj4E|P_dQ2~ zsno2XXD&6Y{^^!-nEg?~NsrIBcWRypM+}Fnr7Nzpm+4c9!mu*?nX08OAwYH`$0+`> zOSAZJ9vNPP$qRy5TpI1~_+;~LuIv_0N-5>zL5jGFZXW$7RN1m&>q(_VNwOrDVbJ}R z*KJH6R??YO(-*NTWpA$5mJ0yAkTDA4xWm6SWrWrdGnF#q3guBG95xx-8}uD#N!q_1 zd)@T9vsKpFuu8m zLB30AD#>XZVQXt;NzTa`E{#H#8;1gvp2LO*r@O9e(A5G8ScAAozV1TlXsNU)RVH7& zj(v6^Ml~N<<6i(!4uVb`FZ!xiPzYXfx2}Gf?wHCANn4L?Rew(o^-IRHZ$D7RZvhtr zcG%7xEc+MC@8at3ir|f6<=%9cmGL95HjSt1<1^s1yv3ix^FSCn`RmiZsb>&*H-a#J zd|Wf>Zd|RN$IjE!^SpW6ETgkZZ7F(x8x6ESowRiq4{o37w+3hXZf&`Emw^#!r3?nb z?s+D!@}`F54u|+a>Dm7F)Np!D$YKtEpUK3XmFEGy(6y&Z;K?oa09V>eP*Un!FH~?qzr6s&X4m-_$D^MIs}e zcW&Xh8B$N%@z#u*M&@mH?THd8awwrr4e0lchJW7o_BMvRe=}8NksevVQDxzRAc(q! zPIaBfta-1qdUWZ46BUQNm27EoJ)q~c?S}`%r&|<)MqnLgH;;4lE3=toN5$a) zpMJPnf$ng$V*Yj|7r5`h&|a295s3KVrDM)_Y?3%881e&q*$M^391>w1lAb}^Ja$+S zn%vwZf!Solw`f!Nxh&sef}>Rcgg^12P5!Qmb!uU7O)OUWt)D2cvQ!l=Msvf5gI zN$=UQ?)!u>&#*ny$lZp9?O<7n2dz!aDgn8~YiY%4ni$dcSVxP}6PlXH1!rLlM!(8i zNGI8aNI!k)+w85xlZpHn{f45d-52cOeWLRh>#+AWPC&g|fC+e@jF?%;c<{tR4t zZvXiDwBSl(QqdJpK-JW~j{kJaUPlie7DM|RyGR&@!Sn1`cm(* ziS#1hC0C~CYOSK{HF1)RKQ71J@GTK~8HE=@cW0f`#7@5+MK-cDkLiIjVH zZ(F*sd4A5X0q(>ujwq2f^|(EyXhc#X_m~MTUw1Bs&wQDlcOJH~{lH-GO}D{YcAh4c3-%j66)sRt+)za56ZAP>gRThvppB>?-NY}BZ5(n%KOxMfvlbJkp7 z!MC~hi#(@3OZBw4Xh-`SUxSiLyt3VzXg{(i{3}o-f_f0AmZBKT!U@4%i20}r zPgyz1a+?t@TV#IU5rp@kaN7dQO+DjiqP}#r`6^n2G%})wI1hc$M&bh+d5$)R=~lKG zLMgca7Ul(dQRg(RQl^=es$xnN)hU*KL*sCse~)g^^EtxSMeHMsoKMQ_2&#UV;uF=! zHe|s7j`c~78C4##b)OH4m8yyHP~`_p{ZE=z_n`8VA`%LFtc-BxD3*?YD9L{*Tq|=J z^Icr#5psSge_Q~E$h-4<6cd4R`d?etL_L^2>plpM_?0obrZWC4^i85@&>-^&|4t&iie^+yp|Gb`t?K-dbV}!{lpvF8mjk# z+lw5M@_=9GQ0$W57(l~0NGvgouvIo-U0i+vF=bq%`tV_u#zk6 zsmF~yxVVB#`9bcf>DmI+ag;>?FFPTi!dgty9zCszHvT-|w0e~frEt+JW3Mc(?7jhST4$Ee>=dW&18Zmz28 z6{9#<%WCeL)0JHbVxVO>d)2jX(&ldez{pRKgKQ2j5h_7K*A)Msz9xl3t zGXU*K*9jTE85aV_k3DS#@ z;WFcw)z3oRW*7anb^b5Y*1&1>4sYf{oHGuIVZ@??juY*>-ZMGUt5ZlU)ht%?obsHnu{Y&1HQa3xg*cxASv0o(qf);u`6^^|TN4 zuJCdZ&?0+-4J57{yJl-tP>K^B!g=a=(H*Zk7vpu_e@?Fd1CVv)<)yITO*X5z)et9Q zt?B;=F8v4(TAIwTMDJHGrFt)m*}lf5z0_G4QkOHf0M>_P7f59R3j3L$qhI<4WQ(2g z1C@8u-wuXq$^{rpbci4YwtjSa425-3P$fqmz87ycrB}8Jv9e(2i`Qs1p;oAE24hA3 zw9y?&Xiizjm_ZRyk<-(`@UrdX+Srg_7jLt9s+swKdS|*ATa%jj$#g*7qco4Ltdylt z7Vbq~-PF#~H8ZJnWk{)AK5tLND- zI@0EIbmc6$R(VIRnHJ5pd*oiB=`C%`@>|RiDO;5tVFRgYbNR#BT~0A%TKsJiq7+sja4_AMW-Cy6{C)xD=e{%!E8z) zIO2IxB)tN2?J_=^xnreD#5RlLT)LL2ZG-PUM6tJ-ptDb|!3DEU!y93STgJMZkH zsk#XVNmX;n$aKuSD_L6oO#*MyBosZr9`sp2>vvYK*lm_SI*Tsbv)Ah`5ITwb*rd$* z@xgiQ?0gWs*4Df5t5N-p1uTwRC=~lrPLa^|^$wY@DnD_l%0EGrrlFkSv_Kk>^|%1b2u=7;eFLUo7_P1lWnx+RKx;e)h7ASvWX z8AHqrox4k1ksV_>v>7k@tW9x@=2A-jL)4G_xyU)kmp7i2jPPuCB`G7-U6Wy~R5fYxPXZGv*`JbaAjwsaM7qOSfc@;x|8S;Cqi zf3wnw?|f`&2j8>R7HZ!%XLgv!Hl|Hwi1zn%LJZaBws^cKuO%xx;N1W^uEQ~8JBBEE z)dH0`^vFxYX|iWIw~`qW%hH(2o*3gS8salAz?s3@vc2P*aQKt_LmV-jySJ^*zHqJL zV2}oX7TVgkR!+S@x4ehfaVVJmK?-l6KW($QA;EW!mQ&s+b7fuPiYo-kuC{S}TLiyb zU_py%3?Zz1Y3MDI-;{)~w3Z=fGC<9k*RWH=`E{87-@2%J?k*A!lP%MN#wq9~uz3Z9!GCj885yztUrrIMv_lgRa zE2?Kru6Pb^zjP+2{+^)%bB&`MVhw{o!NL8$ImvK;Ms%m4}X`soD4S^j|qVS3?-~z3`R&uVq z7DN(w;$ZInW4)Dv9C^J3`ijhdue6NyWB7W)lL^K z$8Ub_6W0{!>a8P98*ddqQavSY>+|hAjabT)G98u#CiHx|7^GS)N{0wb9o@ln4qqJe z=|7q318o+3e|ejt9ZK@Md?6Nj-Tmy=T<+d@?oM$LzV_dIFO*?m_D|XMr{B%qCgt~2 z?mykg2+QKV5rQ^~o!ZN19LBg;Elf%i8 znpt)vU7!_xN^%>yyhu6gQ!m8}qHXqO&!Qmei5XrHq3jFc{CTgSo+2a|h}MZ&TOAd< z82^lPJXZ`|DJrqtx`n!;LnefP`+Q>7>cE%r0AvKYT@r)wI@TBE_oNo)*|<*nY_iS2 zzI#!@A|zV`R;4tG#YL7zWgU$b9vQy_L@$oa#~{R0C+?5H^X!{AYYsL$k}w(gRj zGbICNXt&tb!|9UZY|7$AKM&sAA21v{hJFwAuNCeJVn`$72(DzWwe}d^MfRSCJg)MK zfF8hXv(n-~NFZr0+jwlWHF_CX6K;hrE=E_3zA;W*C*J9^Wha`B_PuM$bt__#$TKYl&Ip| zBm!#&0?d0+UUxAY!a*waEc;^CD~&1*OBkB&VlW9O&u<*`YgXzf-t1^iGOkrJYtvdG z<-|sbEut&Bjd!ZlAl8@M<#=wm>J!C1$p#3Y>D}ul|xq2=C8ry|A6GQ5sq@l zB#A@(6_U%u@`Vn=w73#;F;QYI6j5ljZGA;xT$MyKJ$flLAL$?K)4=0 z$hVP`q5~^ZccN=K>b?-rvAwcQybQB8TQrH8sGpf`<2b&sO1CVpLd$~tyXf4=w(^$P zB&yIFxn2UVWsoSb6_Hg4!1t;~_8-v8Szl0`wL!@>;0StBblolZXv6*oqvH5?i`aJD zYg55YOoLLNGKU_$IrwpEY-WBoGq;02KtqilD`!! zo}c$JbBbhgsupn$uTsJhgcH;f58-`yTlph*+^$N!d2nNiW*@^$p z7G0~vuWFRT4A0?*xHhRVYW%z3^QAB@P+!~Q^9xVL^!ICR3~(%b|2_&frtO9xTi zFIp5p8)SpRT+Pj(v_%HmZX4N@Z{yk9v92-j3OZBxYvT3RM=p)SYKPf%Qy5T7y_HD6 z6d5`xIfN~CR49P@-AJg3<_c^ED{x_Mk6U5@c4QR*gq#vuz!V?BL{ zw-4=CqzIF#OeJBwJvn07MCrg!vQ}19boR2n+o*u(!jCxxe_VR_@%z~EHb=2uC#l|! ztxw7Ko@gC*!*r%K-4xE{FG+ar5J}lLwOMLnR7p~L37*XmBq2FK;a)%nV?N1Sa)sZ| z-75%bO1LoP!SqBWzvq+ivvd?4%eEkPbaxTi*0fg}r(bMU!9lCO(;LY#@BSM2&?W+arLTtZ zb&Xa_(g8HI(~B!!2ZAYNMEXHj4Jxg*p=*WTVgKv3eEE18Z#02r&gcEyTdMRfrj2Rz zJ7loNp?WiTSfRZ&+8FeUM9PsiZv;!{7~Tj(w2j%CW4N^ZJ+lllIlfH##XXv^Z+PdS2%wl>6;2eOm-++~E?32fXL0 z)y!^dq}>2*Sfpu@I%mSJRWpyrx^%A4H@55H<>7r8A|}6X$n?Rtn<}fuN>n2O{wUCS z`2*+cN7ctJbCCa%iJ0zb+S?fPV`#MWrKg}Rzoq&?NA`zU#mq5+@U3ZTe**Ptja{@+ zpPk|9_pcgV?~C#4EyV$qQ70WHd-c?gyAffr!kupO#XyxjO_|lJdeIn&_acl1VV}hK zn2l)7YAU+DTY19EaWmghG(MDwKOPJE`g(eKvl_0B_p7H!W&jSy(rX}pThOJLPP$xy zo-QtWr7_6xNt5T43c7Xd!iAZ~TuFC zF{Urzg0eo_(V3r^?$-<<6GMU+;ER>&*CB3b+0i5BQq5&t6^8EZ9^9XOxmWnEGp9Z; z3;&eD5eZ|$6dWA{mUhuP+HDa1^j$cpN)t*cn8!P;JVOLY+%NjQ>5986rL0J*Cf5}b zL)}oqEoFvw&zgHQ;#<=r4j+7VA%Qi_H-y4M9Wxud;@NKLhYn*Lo^YW%68p;64jqAQ zcJg>M?=}97qJBNxXWdvHhAFp}mZSS(ijz*d>ibO2S{rWyX0An&oEJnwdseCc)WOIm zf`{N%4xJ_r?N>H9#8FmM_ae&=%(244b0Z;FCj0JFE$ZmAMJUtE;pc1CE{F@cJqH5h z5M^wCFS%ql9$7v!?2fQKY~jDhO|Q+JyH;kO``EK*Z?XmtKaBH|O+L*Wao$Fs3k}li z#Rv0=E~=iqjT$enA@#3S%Y>IOrn~Jg(>ZmyC(~r)2i#%GIhcK2TAwb4>gqA`&FAwI zimxWvOXd@3v$B8aST`ATAlR>N)&JgUeI0Vr^n+Qmrjz&u&uxWn!Xu>Z^;$5BV=aL- z?+Le~^3%KnJ0d#J$yLo<kDTc9} z`|K7a*Icf!5YFl-ZpkZfwAR@IZKWCUtXtVG=lzD!A2!N?^1I_lM1lV$@%T`59o#@U-!B5&@vH zwj9yH{zR_N?uovmlWbqHniuRGrSxP<>N^ar zv9F5Y$Tk5pP~-PuZ%Lz@40RI3jCJR+4lC*XXlB(*)RWF`B|EzjxT%%{{TSEL8|V|I&3emHC&57BCO83ikP0 z1-%w|Zb_~zt>+f}&S3t?)#ffj$I7WI(9Po(pe$Al)Xrpt^SAn_+S@e>zwubwv2)bH zqqaFhirVLN=(O_M?<4y;mT|qZb<-CIElOyG!B4B084xXfQ8%eI{1%83vBHN_SeEFR z;LMRXOy%BzG^W0sgds0VD;tM@!BaY$2+)0d>it_qf83xP%YqLD%T0RhrtA}#fPH9 z5q=56ZWg@wl78J2FQKa2`uvINkOq*fg-ed9oh2U!?9K!juK-+#cDe+IYMXP# zIUh4MA7S9Cbm(n5e&g0;%g9>K^AbJGUawX6n69=>S070B{7Jeh<|xWszuDIWaYFAN z-E*2pHj7np9+O2KQqrVT{<#nH-btevOLCXx*sNI+k z`aP{Pib-tMXbLr~F;A)O_QB>M;k>&#LgN+2`s2*b5z^o%-};K1N5+MAl_LT06M=XBGSZc`qI>UMu=X~AcO+q2Z$)jr>n#fNPi_ui{&*WrW z!LJb_YN^Lj4)}^_Lhv-Iy8b#!*Ihy;LmZ?{p}~~{t-sVk82P5J=IEOChkAL8STr+E zOGPY&<-fiLMlm5yJ-!Ly^wt`bsG(Qd&^t0Xc`rR~EB9-#?Ot7m7H2&qq%IQN*1R{` zI*yP7@)u6(GPv|n(*zE_MQq(+-7LP!`NFBhOzk2q8N<~**^ab)~E z3zJsS>;HCZMH~w$9c8N&T{e8b524WJgP8vOrI?Iy`2ky_78&{#i;s=RCN}HiO%69Z zoz;qS3XF?e_19-`iSx^ZpPZD;0Fhtbo=e5xyZ~{@h5#TK}%(6!%my~K6wPv z^|oT>6FZfYW*Ouasc-A@FR0??w|=AtnwTYN_{1a@x-QU%?HL;rFuksoD#Z01G1{{5 zET@?bZdRIm!%_t1ICe>T4jHtpEXCGtiyU2C5?q8OQOLth1T1sZ##OUKEbX5jAdLenCtdyMG({n871I8 zt@c{=a8RfE8M50jiskpx_dMPk*gCcyml~W?sbI0%*J?ei{5bngH-#uT4Zl35Plo-)jsJYM z;OW;;sN2~gr_I$cfzV>81Kljp;UHtmBz!@E%yiG#l zPirJjaZV~v{3M{sW0jvif3i9Cr&P#qVHtVF(eQ#vN8^RXI|bk6AP4c>eHzd3S;2A8 zOZc?vW#NokTh-VipiQr_Lz={j*KTi{ZQ_iRSt*1=e1D{iPjt|;{7|>#)eI`DFCt6U zH;)9dSdphOXUBreLz&%MCCs8QHk|8~aS-B&^)8@oI>8=-P6U~n+|b#jSotS{$inX3 z1S;yFplh^P#SGR!Y=7c!QLJxlq`EC*UqC)nx-s)u+Bm*E81S{uX*`ufT@9vvMR&)0 z3u4n3I>R=KqtjeXW{eSP4ZxCMoKZL-{rnsxgv_|r_IJ@D0~*DIFAxe~=YDv_JwDXh z;-yc)@q9E~hCr+Y7X8-Wm=Qu>sCuD@-&+j$eflY8iXbnFaKc}1yp3POkds^5F0nKL zoTyoTbF5OkqiPi)4~CiBW%}X}@SPr=@m?bmqM%PdhbunpCGI4pNbycLN@MYkf(2I- z?aIC4!sW{1%ob({RzjGPLY}zJw%Qzij+j?=BdxFA5#6X$C7y5TIT2@TR-x?#fPNLKL$wa0!~88|j80T?b^EzY(3Tj`ld!POshh57+wQuwAoHp3J69CUeF zdz((3T03ih#+gtzsj8ix0}90Su1W9GQFnLuD1C)qqGemO84FTW%I%zre7eZ=M{Mwl z{0)1>Bq25U1*bmaO~SgX!ZJMSE*(94wSS>&gm7F#gEtgy&n^L8rs?WQn$=eKP?*mD z(Rt;bb>Zamz|Qh;r~Jlq7q2s7y&B*@ioEAGS2ODEKnI<*k}bp7`pH(#JOw{CP5Wj@ z0O383)esWE-O{u=bQv%v)qcFMlH2}d|6{EEM^3hR%Kq=+Y}Yb=A_F7FpgNM8LU7Kp zViFybhd6#5FMsG1x94cSGU>JmZ?QS*wbh8`7`)^WqTub+)rW;KBO1gp1i=Q6^BjYN z5Cd_rj@e&EkMgYv5(w}3@c__Y&zMhHrM2B2-~3*ZnY`P?+5_KKN*R-N1Y^b8hy%v3 zB-LYoE_Ns_e3$~qt9w0RQ|Cqtw6pq?wR0z<;V{+}|GnCjfDrr1!Czg2Tfe9BXSsEO zzkD9+>Qo-QyeyY{Gqslz;{+koUEC(F7yeLZGoxgy`~AJkT~k!uKr~CE0@bNFUiVZS zv#I01(a`hJ?nzm;_Vk?_nVHuv;%kbMSmEVp%uf!8mN2Ci7v~~!5@+~co>hyFrgZ4Z zXT|CE^hlOt^7(tWxi zZ57he)~e5Z*vZas8HB}LRgLNkm#4YfHK+;g6ZBYR>w*2Kmjpp%(ME{-E38Zby}O#Qqb=h zo7mogH*^1b{+CG~sEtcwRBZfZ&0^kW({_Y*yYe-QhbG5^xr1x*c=~r|RX+>(iRhO6XRGrzYjSdhO zlb>sB7NPmIwjKM#uE^p-L3+?`wL<@+4S!Uylf1>M>q4`kA2a=3ZY@d|R8ST6_D6C>1?< z@5=3%fjrSCSpzk`atQwA7H;2IJhuoe`Yk|AwH4s4v{Wb`2>jHy zpF_c%=5v1UBnEJy_6%T5!+iw_Qj4AuM78(VgUOj~bpRbn_3EJGOlhBKH>@b5Vd=B< zQG`R6yW9HS&fH1_ZOme`qWb%xKKi;#;cI5w!kxy;a2-HVM~{>+Hi|~Bdu7EiXO}{q zKwpED@10}FHD-jtXBLIib(S9*(zn}%uo%df$3>s=yP>7e+!Ng?(rs0+dptBvJWmb2 zi0g|4od&96PWO)?72tTE@dkfBj*l-Sb8e;pWkpOM)m33atRchszOT(L0)}Ru4&&y1UwrQdf@9LWoC_*?X;;dyDvho}KN?x^)`N3_Yq{!b zD#ViDSASjKQ%j^chy!mv+PA2EMYgMei^aq487ZYu%$(LEE}wp(r4Lu@sEM>(Q}ZB_ zW`*1lJY%zmw@Uy5MoF=?^n#!+3@7oCa&IYIxkK|CdFzyk9OT}~3(o7$LZy|&Z`8Cj zvUl(u7H%ejSpVF^$vQhtxni02v4pj4DZk1Yf8u+tBOKjcHN_*DL=-I4Moh0z5H+Za zL90OWSl25H`uXedzKC%|-}Qvf%i>jNZR<2JFg&;~AkDF*3xbR*T%MGd*WP}UVhQ4v zcsOfDd+NE!oEGk5{7TS@&R^!eBFwflGonO0f^IJ+`qQT15FW40iOAz*CoJ->2%HI- z8L8YRlD@5A7ooM)Ifs&>-4O_P%Hxo}U1>LuD??j?dmN!92(EZF>95glbl>xk ze6N!wtW*9nYTB!QXw558lDuUPv_z*&VzaWMW`jmEwkP~H2z4!x0T-;PCYx;!Uhy`T zSgyjW0AV?Lwq|tVNx+XQji;_@&4gx^z=dM|Q%9#4FQbdY^W>^NB+3iKdk4uL#J+JhgJj6_Q29x5fE9> zaoXbdym% zoV{{e$El-$qZxSa4p+nb@p25++9v`sR^YV__hGXWRqHtL;?i-u{4w8fBV9U|*>$!| z|9A-wu~xz!vsfYa8JJJ0GUpmDa6p;A<~25e8~JNeZ5S)Rdps0iEYF-T9Sw6S_=`KS zYR)=9m&JKf+_C*r$6}x?0Egj%vmig{udaJ~dt&T91Tc58n|r>SIl8dooJFtnz~Lp) z+*Doj?@`kEg z(qD>cTCXpfIXofVrMc}*nT276d(qz9h{V5vcblXW3QkWg-igxcHQ+;|w0aBC(EFm( z>mM^|<(OD@csY*oA+02XG`en09CisBLl_T1y0IMG$o=uocL>OsAR19VeJh!6wCtQ` zCss^e19;{>qZ`8cAo%iC(I<*@uh*Z2+XYj=EL&$SDXhD#oq8=<9*Imd5#LgV*PwWeLcP4(%*1mnf9zF*c8NI4w{9q0zh3O1`uL?`R=KrZ(y=u% z??L;V_jEXlAh()(qtaL@)&&s`>!>-mh zX}8-Fm(eXH8$$)a!iYe zP(NkrSHERP;Qs8yUsfrd9h2wO1;S|R^t+`#=lP&X;~<%nI#}H~<(3bj`d*IcsZU2y zpX%6l1)Nv0ut??EORCP|)}%zY<&RryyGe~~aAsBJMQ7~N_HUfLqT$g(4ZHksZU&e(~V^sDSB=QM2f zuzYc9$652%(D}Ubkf(DW(iZg<#lJHOp1BUPzO#nkkPywOQo*Ct# zAxYT~Uhm$Y>qYw-7Q|IDHPz_Zl}`c|#txen z-^)Cgtect28b!*OkzrKjn+GAgyCN497X)U$>-}0ArwuIM^zz)@!9%W1@OIu*E!k=D zor9uxe0A*Vji_G-%`5_T4J?<~);3&n2B>oahz%HD>O{!hp~UD)^2H|Xq5B9nN%S}= zI@WUK*%cDc^?27^-A7}0*6n?9{24z3qZj!vjJSOL<(K;jBx3WWMg18ijqa6ihK~>I z0dR5UnKhxYmWGB6O+9VaIi{r($zJ1JAMHkgz*}CHm7CG1M@QNAV5&0v`IoJa75K`T zxt{!jfE7o2t&|M>wFqZ%_@u%yQ`k`SEWSXgh&xVUxh6u=nGw#C9>L{^x66y1nA*jw zuYhntZLOg1I9!}}f5oZ|se^;wIY&DCz5BA$Zp9Cxw7$w&IMO6IrHlFS@!aC+o?9&h zU>A(JqqrJ);JC+`Q`oGCS9X80<&(9iKe~}M?TT4qxtIx%FgJQ*n+=uePK$-SK*H)h z=&pVTz6y;Ra^1f135~RJL{h#g)I*ZnGaiqQD!+8Hlfhpw;M~R=y3~6rkbKCwtfl(H zwK~}COk3Sv%I>-qeC1)(oSZo)^dcq|gkrXSP90d!KhmJh_HoU-6z~a5{Nu6troZ-ezKc{wz-1B_YD`IhT886Qc zamxee){H#}x1840uU&H^F8K7l%}(}_D;jJJH{)eRLrZOQ^TWB~U4CSb=ZA(|UzWlL zor6~I>zS2}bEJ&auA}kXdNEo^SMh-_iZ_>k-R)U5N7`Gl_}SUs?6cTJ((ZcTF?b?n zeSKv)t>hpj)AjsrSBmpyWEJI}CsLOm97%mv?&0meVk>gwa^i-Me{$tL*ROw8y!wL! z%i*s{F~5ATx%w5UW+kmMUBXzsF7}>hW@A9skuZme()1W*-A?JFg zJ)71+J*(mc%KpT67z+v5>>8z`0HnNyK^Bjm$uCI ziWW$VKZ{!f*MgNipCJqeJWk3xi|Lq?y%1Gdf>gCE{Qbs9n3=^%b4__2@B8{EM|npJ zUY;)GsaifDn-8oDY8+6$xnH3FaQo(651vqNXbqcc&-qsF=*lWtif#jgE?zvT?u3K~ zYoS$XJ6FKqM-5vBMM>q1&U`}3qumXykgvQeU-l&Uz}3dOxp>zfwJTO8Z>`R>m?Yp7^1={I#yLj;x-YQNb; z8;QPfGnX{~<*AH&F;|~e;^%6U6>;&wb!&XrpR`ue?p7Ads|yP!cHOte=;uOu&!iVk zEV^w_HQ#5#S7#>DNUMt1S6ATIZj3&5D*IhgN%3Y&72W4pzN+WQ#D=EXz>L>l(Z%it|^R~xpLaM?BI#;K(?Jyu76wkI=c7t^sB0u%jp%nS2C6c#AsO8 z=nd=<&9+$WmWr9iKq|HLOv<^#y_iTW@Cd&CO^gPg;(iX6)x1%Bh5MJ`^2b=C)Ynmq1GNKZQ(bvwGQJ0fvjduaWI;$p8W{(+%+`y;j; zX#hOBH)r;JyPL6X9j%f~*F}6owWil?3ay3nW4G!D?w%APl!g&JB|HmZ?`64fJy72_49EFwV&sb%dxv_eOlWer(KHv zwjn=p?lW~$8N2YA^hixEo75erT)v~XIp7{`9tqWzzs4IN&a{(xvp8*F_?qZyZu|ZAaWfVH^6~ZC%afjJ5_si&f%#?k z>e9vBSK{TnWA9ul+yviE&~F(IJVL;0J7!y3doW2KeUyx3vb#3twGPJh?~d4XEZq>$ z;&NtX^lxh5-@Qk^f4c^veI63ckCh=Ep5#*q?i_p}e)vr+fzz{w=(ZtTmlCEZd+!zq?UR+Yj+K>IDf9T#_F2W)f6Po?K{|WIXuiy!|zYnp2&(IQ@-nKs0 zH1B_@4sRFJvilo0HTIlqQtQ6xQ+5;t2IDeI%^9h6&C%s=+6~_vn>NWUfsUKo(h}bG zJ-Y-ipwPI9><3fo*MFf*1;3B&k>dyof_ME*v z|29r>@?FUXY!6eK{(N({C*1#j=r_EfeTVv~1xFi-A^Ln;UpytsG{FrOJ>R`Pf_i?j z6xm*VWZUq>wnaA>=20v;pF|tEd#8icd3em8C>DngLxEn=F!ZqCcVY4 zI%+0=&#`bvXwqnb!SxX`cDn;U2rUTr-4@-NYk&R4L~+Uz7@rZ~eX6nYL}m5!Awc3$ z@x7X&lZyL`W?K)abr3uk#`>hV6-wOKAn9vjHIbB2=vJ>Aa$*~>swa>?*qBzXq+NY; zhm1U+A77YuQJ;0DbV=yjpy)H5D*z2}NJ;|>y7glNao;T880S%NImCX#g1jrRG^W$L z2gEBWAu^DWgXe6T?xHEW0>{U1w3}O*jh;J-Vmo_=1d03V23?y6AEp;nd;4tPc!#xA zofbiTJ)#Uzyxt+;cHGQ7Q8j#lR9RJD*@!NAyi<)koAe^J;Vf0WIl@md7V&bkaOzU9 zW`Yx^P+mElT$Z}5>s}hUD&^4lk(9`^zEqrylI}frw{z|^(V6C9i*rq%iBBMV-$~SL z!FKVa`S1wBcD+j+y z*V=IgI^2s#!uufRjJ=?gp0L08;Fz6~7^l#le7&*E7L`(+Z$MW{0Cw-K^Avj2Hz||; zzV0D3T5(`>Pyf|tflaYP1BKi>9kprtP1m=OL^12{e6~>;Fk3LLL+7H_Y)ID^;1()< zw;uzN=IU^~#8%p+`17qF?W)@w6d;Y$Jo7#qot+Z6qKp!Df_$SRqC;yGC{Glw#BR!( z(%>G=omu7rYEOZ&YJ+lKmRMO|p?>Ad^qI@K+}!!@u*bf_evi-hordo;1@fLYf>$1I zZ&-e4JWw>G>M0rZn|HlD`SN7mQC`Gzy@k7NItt$U33*%0o7>kSJw0)K-xbbl1BT== zBO_@8Un^;V=OHo>vhPQp4UE_ce(Xd84@CVq120r#V3lFxg<#VP@d?BW$-%I244HBgbAik?47+-v5WG_Dh4&;8%j5+KmTjZ zI!vdfF0qpxwyt1vp)$C$nhPXE(lTi$s-yFgR#mLhNfNVZsuOjq>hs!#4GjsOlLNGe z(Wh2W_UztpZ58&NWiV@kx{l^|zB45;dNg2ZTHE~mCo37zgV?S#Wd6H@yWRTP#!HLi z7br^KqAQqj1zWK9ysz{EYuBP?@Znb;ORj|Ig1igipnH^tkV|1hI-#TPE-6A7vgXh| z(pBrYxOMj~bz0bl@Cwq6f;88&msg@xmabZAPRMKPbd+%?6w&((cyz$2$KXzZBJxmQ zah|KQZAad+cb|*4QFJTr9Xf<7+;+o5Us32+xm)@HH`oKih3_5yD!8q<$Yb`&Cf8-M z{&*b4Enu0r^S~U@l5Vk_uGad77o~kw^or~ppLr1L4R!4+F3AVrlHskFoHI{9iXjPl zV|cB~XC(m5d|TC?`NDFR*V>;?T26hKPOo4)2ezV%de!HSG3E~MZm3*Wy!lGi(wECB zW_g~>$x1(b*>))~Hve7GMAU<4n1LlzB6TOOO+xBT)>gN3+tyYP94c}4 z_5@OU+TE12tf%+3(=_E+c*`b{b9ymmC#`)8Y4qxw&h8}%r*rHYJ&>BUeqA`W72S()RB|Zm3VMPXPr7_;^OpNXOKNkQHecECdg}D5(~mY^+A0d} z?OjiKK<3}vu!VN%_Jy>XzMU~0Z}iZ?f+#rN@VfBS_VVKySRm&4+^#z($4=UqLvss` z^nRGW?77qb_J#XVOPq~2mrMJ)`WI6=I+m}rC&Fo!2`?*O6AKkro0iVt_=Svd2QF&< zndS51{n!&52>zi3y8OrQuFW2+UK*-)y@;84z#BUGDlf8najW`x!YN*LdjhaMd};Y~ zR{g$>Lx(6&dRG;%NRLbNS-FZ9e>-4(4auh_)b-*`QNg{iIAg5k+#I>ZW8LA`v+;|@ zF0Fn{t=myDGXt4lMfAYb4-ehE$3#x1vt7lMp}B*`wM9>Du4&4Z$Cwwj(TB4MK~!E- zwn5}i=AN5Pw-mIvT&~9_Tg1hdl4iG{MX}B^81HSH7t_AMdX!R0F$RSyIVv;_BB!?) zCimtYI`I5>0-W~Np7XQ>8nQYIUA?=S@9hoUdwk{Ho&&bQnFlLS=59j^CTOoteSF4k zcci36!?0XQxAreyaEmcJagqAE=*BLRvL?Ud4EJkQ0H)+Y#00rzu5wGs)qaY_vTCDv zT`|}s+T4Bp@NR7B)7M|<)}Ppr_Oz{ucQ5=CYP@sOayhc3XHnC(En|+|#bKfEshkfM z!K(WQtj>P_MIN7Bd$?x$`JF`foz<`n=e` z$|Pa}cV>`9nz;=PvCivVJ2#+VX58n8^26`2{o%W{cS6$gBT8S)EQ)(|FuzlksGq;) z%qZIB-W;Jyv(9l`a~q>c>k6)*>v_1A%luIF`n370*|APbahh%Wp{b#DAv`d6sMQ5- zjmk4HHkhCF zqx-`WA~Q*$_8i(SDy)QylemTd(9f>Q8i;E)|h3ovckUDPPPxqs~SQNRt^_9`n9*hjE3#RzZ zZCurXeQICsNjupIJ#G^!%QJwP;LgJkh?t*oj79QObztB?2wHNkZ z5k6Ch!X_Ww+a6)KL89H{TJ)DGS*trwY=}OQemn78)==K-{PYPH;7}usIQ5m z-&hw{y^8tyet zcvI}&?EU=y>Q_BIH!p|h;aHuSkgx$5_ujDh4$AZwq!|P5=hM2n#;ux`o(Ha;Qi{jD z^j^#5?HlMf(r(m+GY)mug%2zUuUK1LP@O&J9c~S=G+R1d=$67`r%I#!G_2bSUL|RJ zNN&TCPD^KIqJ#(b{_C!O}(i^TxO>TUk@UTBtR0q(;@-K4uhQf4T$wBskG4D= z!tw$*&IaH{L}|92Jjkniwxa}Bow;{ycks};tvVf8U|iMvs2X=Zwr`$CYWr>#dFagQ z`<|QAEujnBldu`NlCd$`C7Q9Y1r5!V-(ez}Tp_!{-Xj*Ad_Dbqh!&6~=_%=~&sOtRm$donTV{sS5(s#&JzY}ZO$m!%@$jYY`ZZ@k@#$r*lpa39P4PT;c<_s> zy^F&__F_jl>LvZ5hYGI&qv>OO#pk4Jnj>5L<`u~2>G7j=>k23y^INLQq5X%P(|tPb zEsMUi%zu&sepc~i(&Ng&#snav$TrEUxO21{8+Rw7_Vn=aIKN{~M^Lutu_-FJac0rX zVIHZLi}jcjhmrOl8kz6h^zLv##gTUBhL+teZDk8q6%}u-rD{UtH*4n#^)ca@fvf7) zo@q^a+UR@8IkCy7n{IPhl3yZ2~5U5xO_WrP1m8_R zK{mPvf}IHkMNvAaHaj|yuy1ol-qut1S(#U#?>uzSZ*Lhd|J~KExf?dZ>vk-CckRf& z>YKitu$KGI`;@a(j!#;z7$`||PPK2o%q{nd2sv_c1bd&N`^K}aQnUUBRX@zxI-thbViU z=`m56<+(#W!?m+7(WwCquRW;Xd<&)WuC8$?B;1KDW0h-s^Pl)1om^2FalX!Z-$`&@ zDT0^IY_1?KO^ozt+*{C=?(}`X0+T|)0G&07S6iAIDS)MDa0RzOSM58L_p;vQyu9gX z+7)Q)#+T*dZ5u*FJ2xD;%)4t$n=MFJ2W&qyq5$K0C^*MNZ?ojX8Wce~EvMlnmI6}F z#ihKXyH3s)9CD2f9C%h5FgrW{>H&MU2e~4@bJ$hS@^ej9Xra` zmnU5>&?LHVztX$Lx%qPB<}1_{fkB%)gND2ovPP~W(>4pW_}~ZOB$dTf&uR3k3XEDY zb7OPwx6YaclpsZ(>x&Q5H*8GUZt3c6otf3)?Iv$@IZ`*>v*7eZVl?+mLAuv_;H_fP z!+ejG%rnQ))orTjyoTV)p3(XFh?@Kv@2t@D>D$z^n=&;}AkLLhrS!-^=ZRN3hnXIt z7_ceKQV>aqqDK_E(t0nRVxW^1r>;$DhupS|-6cV$-PI2+ho8!!&&j+cIz;nH@wI{G z1>b1+Pc|&~>?i;d^#QS4(>k8hhqn7V`@E5$wVQS^8HD4Mv>AW@)8pt8fJSdW%7g0c zfTKP4gQ^sl3Da-0AYeN}Y|GG8)k|NL?g(tJP2U8~B<;SG$pz0X0M}RVB|~}>_S)E; z#hjtn1KSp2u+@~>z-y=lJA+PZ-=3ep#~tGrR#^jdeNbdiPrl35G7&((AM*N(W6m_C zrZ}H_x7X`HaXK|i)v)=BzQ3!uNEJx7tbxqii+Tg&urB24we@*9m{*$fgV9KTK$2ass<_NSi79*^crV@$nW&jXT=z%?v$lh$#Fd@Yg9V9^H7B+>5pqGD>lxXXS2uXbIovy!H&DrITxnUou_2b)@64Z`QCSX_7Uy0r znpP7n#%kxmPd3Ix*W5khX6dx_YGy|IPo~yTxqavIT#=FsSaC^Dj-p9Z24<&Ec(XF{ zBz4Z)A<*5Hhj}^Gs=>nj#n(nQ>~QA?D(uDS%FNYIa~o>W2UG14o_jIVP!sA~dF1jD zU z{HiO--mhNF_Er)qJ}P;-ZA;L3bl0#nY~POb{3lL_K9t{le>d}Th^t3o-&1;Y|DnRR zJGii*gbb9kjZ3XNunj_)^;VaLTeq%gHm)`ez7a+bc`=Kw_jrHL^u(&u*e1M0$8^6A zTdtfbCSE>t=VU2ZIoFC?Rcn`&gS?FygiQs(HdQ%z>EfODVtZfNp@6l`?XDNGr>$xI z5lWJ&u(U4ZaPq`q-}Qkz@+mcmZhiwzgw@k$T|FzL*&lDXnQXpbe)-7wb>ZrnyPYjJ z>F00Lv>$t>Ron}Cz*?qC?BN#f4sW?taOvKb0}I0zZJ_G4iFMvlBCirilhgh{B#FDDLOkx5++R#sai=` zCU?DkSDbU^fM|67>a(Gx-y(0%&V0~~K7^VVVlF>)FQ9Tldg=b;;{lT5YMsY`<`UQQ{D`=L z_1(*odjHVdBd%WHqp;zLG{t~kI31yQLs2DmP@n2Eo?!?vtu>$JXEw$SH|)sXT{N@N zRZ_A@o_dCxzJj`(RwG3g@s10^8o~=gUObodWbmG~wgPV4j(z8#k|b?7d&#t)ec|nw zX93wcf7~3)HPcv5n0@u|>?qwt%~hET;-{kQyk-f2qen@f9ufgurdJ6huJcDkx0XW|>Z=O3>tPO~5Y4JUm5IBG|0z>bWRio2bKQ2>Ft`Y8su z{YpiMZPBhscdj^JqebHj#XPV3d^brEr1Xr=(p;UhAiBHLFGk-{&qYCq`_Wo};{OP8 zag&4r7lLR6^0MoppLlwL#(=tKZh6w;@cU}=n$ebo-P_Znh}31L*K?F2R`%QtI{RhH z;@P&{71EC~xG(%Vy_|OU9O4Pmvff8M*VPXT&+vR(cXZB76c<_#=3k|Ge#c8KJl$A= z!8PmM6}>U7eDq1v%TE-4G5&mW?d6@%uB|Ib{+H8tOq+gs;Xd))S9dPQEyRl;G_%`Lna`Foy>uT{R@l6Q*nJk2A*K-ZEd^HP!m?-3M2$ z^}F^L9y#x9b-g=r_1g2V_J4fj<^O8^w_kst{CLTipSu5g`;C9{{rW4#$FYxF8IK3n z|3rJ?%gbdCSALZK=<8GiR`z<(?6Ygx6|^ACtv;i zYR#|bRv-Al`>&+_wVHjZcyeO(!CN1!KJe*#2d`i6`(Ezj z=af(9evnVi^MhYZOKeZT4Fu*DJ8VV9K14p?UL8`&HAIy?+CvL`u(L} zY&-e;B7Q9JBl#bn*8J}WUl#4j9%3!r7uGm1TM+I#(mego-S2)L-1H_bZ+Grvm#qRyZKelAFup$eB+~A%e?=yI@kQJ@IQ+Cqxc^Wezx+)e|FE0Z`k?#>G$ni z^IZ4il)vY09_6?HdDZ?Rv!DJ8Me_^RUtIoUb^qqUuTI?hOZEI~dHn`0>F=xipNswU zX8E6A_#4^$I_PgC{*9!6TebLar2a8!?L*!7RoVIJAFuo)2fzB&-_`THckTa3u4sLJ z@cGwI7XR7#57YkN9H;#^=fBf=@5KM#e*^Ep|C_e{LC5}ihCeI%birRu{80z3o#_A3 zmA~`!?drc%@c&xuzq^?N{=NbCe_ek{!~baNf1!r|(*^&fcKy#N_^+DizuV{hx7+Lg z83q4U7yp0P1$i5Y@ZcMD_OnmdPaAxXuVFp^SA$r7f;P2KHaL{Jlqf#Aw{`T3P z-@gU!zwL`*o|63a%-^5SC}^1rA3E?uk!43W{_yl4I9?(*PyPeqf0yX8_r9-9s@(jC z*X8|P_k1tpj;>XB_IFGFrQqLaxYh6dM|=HPXu*`+eDp8#La+RHG|it1`ry6Qe>%ZW z$xVMc`>&3E{1g6Dd7{N)Y5Px8{_4eLOXOG0U;owQ-zS)-AEiv9Pj|8af$^E&<`uPhh_e; zkUtza(Z4^o$zRBQpW=l|*uVJTkN5qD(8BNj_saRtE8p(B^zNT8{5t%boY3sdfyPf6 zr^I`6a@^gYH_F@%xw+C|xrmWD47l5Qu)t~(+9NSUiQcFp`rMXdg+UA)PAtbOVKL9B z(VN+O`iX8B*d$9Ou;j3jz}6-Rw?&%zMo*cSKwSd%MJZdG91L|6*eU|TCI$d)*Np*s zt6AAsr=&7XBI5MUE(bum9B!ivk!*D?0GOJ*#CYmBb5ajujl`y3=_g zqt5_Y*hfJ026z1!#%#S<0fTK;qi$R)8U}gY@C7$;GBym*64do} zxqq7>ku4`nvzP*;(kq1ko06Ac0l&rAwg$uNqh!Nrjj@5|tqmrd^ zf({l;vfW^8fC1`8cLEVXZ4Oz|0UHexvS3Wl%!FY+w9$Jw9JVsnO{bsqk|3%~QzoYA%=%P}5u~xC9Ghw{i_hl+VSryQ zq!Kv1V1RdqQ{$s|7rj+V zGJQH?aVo4(4P6LO7$gq5#0wd12JO*kUvOAS0ZN9C2}cCY3s?o=ju~lG5R#L5~I>_N}7;hx# zO|pLhEz)Ovp>?aoDCy@W+a~EiB0o8pK!d@+WQHoo&knWeRiIsO*QoGtrru#l5VG-X zHvy~`*iaqoYl@Fd&6bNHGf9O5RU(#%rs1#wzD&r(+qF4#H=qWBA&fSeK`4_Gm|33n zH#CkjDZcrX#~pStcmY=?nxvW)@i5x^kZ3fpA02N93x^m@c$fhDz@jODC!GusLtzU*(S?~_P6%A0y@-T2b9!H2J^148!-hM%;!nbK0HoH)L z3>;?%*7#o(86*ZNDzZKahu20W*MM3nn}9GS185RDV_ODi1y& z=%8`faa|B3=1aZXjCwMG46@boE@?0hR+iHtcec-nNQLk7C4)wcd4am#FxeED5=S-I z+7P=>Hx7WiD1yF9Q7jHY6djy?(90yUbEU9%lNmZsIF*n*V$H@gsu);w;A;7vnOzmx z!CAgsV|=O@df&|$u*WJET})=fL~bsq>BI24q%a_HI8cvLhG+Yj`%<}*aWbeP)Bq~L zkgOtNBxa_WW5Q{L3Q;mFF_XbAAxCMm@>0P`n@!nowi#0cVG;FW02Z6flGdbfV51vo zIHe@9MMw;%Mo!3zjEdiuJaD1O2vAdmbdXcZiB_udYK(*~=!ONhNm>|F2ANP(3`1fm z$;eL>n75%G0qPD^ac0jI@VAu~rAuL8ArxGc#X zEtR2_V?Y&AV$Ymp;cb%~j@^jo8X zG&}wxfRvh=2?$)nmo4ZrPbvr7E2N@SY&EJ8`Y}IugNViyn8Z?s3LF(nZ0(`Bkvx-} zmI{NO zqgK#n>l*~=bp$w58zp6-tA7ODVUkJZgNm+H(BT;6VQM6Htr%xSYc5g73Z}dT)&PEE zic&7vl^oD5N_DsMbJ$9jP%W?nB6KV>#8|mnjswIpLAQ^Vh-tH#X)x?YB#OZi$f!Q9 z&`&KAkzl1=H_q`g3OdRS0==<=DMKw{w_*w_>c0@^aUJzfMv30OXA}#ez5{3EM?6w|h`CwpJ!5 zLVBydpFlzoyVO%*7nVt>j`1?kj;%+bkP&eh`9jiOP|IurysWYQNib|oCXEPCy@dhPhl41ly55b^5;zLmXnnj`$T7&lFJy9~SSCk| z+61CE&gA}p6ztQlCrjG;P5k55DhG!auaIGS3Iv>|i&Mz$zCRJAw&qW~TRI1Mi$Ph~@!{ z`s!IZ*evW|$~LGrn%7iXr#_RKuuJM_(?Xf?J{?|xrX-vS;!8#eY{YIv6zT*iEnZ}8 zXMk*#SaRKH5+VvMpf_=aoB_Q#h@`N(+eXo*OGSqXAR4FhV4RXVsdVT$!^x~i$OATQ_h3h#)KLRf^9i3#z7=yA_@z22r^K~mOnLel`hk?9>BCeEC$%s^orrZ^JRK>Dmj1ss@1|qjB zR1j8<*X9yvv|u;Q5>W>Wf*P@MKd_ytA*awIA(6M;@IGWTw&}G1Ai6G%p0uiE5}U$` zXCU8upD?KxR#as)E5_NC^y6Bvox|xDj4S&k_DG)OK8=iz#qfn>2%N%u!91CdezI(j zC#>`j2I4gI$Q(8S_D9wWs%-uGNzww+__?;Wb8YA|gz7Y-GDxY^jB^-RAi&tj>xlJH z%NKAlcqP?tR*pgvBLpR+;yvT!(Meb-lzBRolDqp4{QsB5)Uk3KW zkv``EA7iF@(kRD@IHS3d(j`(dKT$}c#lQh#F`-{d)ptS04iCpUH?bTm808B+y^N@N zhRiBRXH=?Wpu$`x%{H8a;wWu(yFBTTz|@=w0BuHzRm0SuC&0Q8^M?L&lDH>VY7`z(37{Chj@jj#SHctP!Tz9RgOK<_hPbI+WTu~QBVbjBazvN<4 z;3yvH62M#`t<)qL7P3|ROuYfGFjoc#qTZd2#R@sW;G_zVMnG_=5me`9>M!zB-2gOj z2DL}OJR(tBUFOMA&{9aq@p@RxJjPADUdL>fE1!{N~03N~= zL+cFyte<3nOd{y1!Vwk>*O!UI@HCEExYua55+R>a)UM+>$gq^gUTmmG0|JBeI?&b! zNXD>|QHe>&3Pe$Y#EQE3(Zzx>-+D&}wHvQgjPfy3RxUWIyn~{8YzH}K&=gKhq+m2RN_f;U zDu+Wy!Haw;o0|j}MRt*wXq=`3M0Ub9T-%shUM7x|ILsu(VMoSrSTu;%8vz_@WBXBn zaFM{|r;u?lF;L$qWGnj7*u$5K3=2#G2El$o1*-YUUId#_-;FbxsW$8P=Srs83zs0( za>Wo5DwzykvC`Wly&;V<$%Nke4kltP%Vd<2v7&K?#%_?wl{cj9+77B27SToIM8TAs z!N*u?G}|0wgi;2O!eO-OOaNZvu!8mrWe{Yx8xzF|XA=@?f6xr09VQsSZk6w>mfta0;X7w3mnA?QR+k`vDp_v7Xcy|OazT)Kb9Dy z36+c~#)C21B%IG=aDt?L37#wnmF1$LoKY+o6$5QGqBfEzwPU@60~;{z6nnS2UJGFO zB4ahD4C3OceFDVM?XJZsVo>rx8e=ndyQ6U-iYkeCMOQNFw5*!Z+<-&@86_3uVxUWo zH=@oLZyOc)CVR#kBvy<`$e~qDVf0hrXfDBusO5N@)rc}3=3oNHW<%pK)Q7VJ`%ZOZ zWMYV)00_EOpbA2|l1cGA!(Bp~#>QbAMMfNfG%6_rU=d7cYa2rWm{f@<9GB=^Jo?sO zi4tA)l}Sd$*t$ftnb-slsZWMtD>zh-9L+*Q zAEQjlMEW5D6T<1KV+4g3WHpVTmO#mFJdalqU}z5=jjBg;SU?omhl4e+mUkWplEDst z4h!YGRHy^y#1iPIRD^Q`G>+NR7m3i~0j2HxV6TjTSi2;MElv({f>DFrp-zqyCnU)A zV_F|$_f}a7NpBRgBn~To7Zxk^9Hk7VO2tqLGg)ItTc(U))A!YtqENQmU0o$HNjU6i z4{4udvL7nd+jWyv#QwciHKDFXEV7%LMMj`3g@ZP*T}oml4_8#un<0a6p)@F%pitOl z%`!oIs0?K(@ODBYpr2$TMo?$tC>ezE3rfUj5TIUCX%h6+wVem;@<0_f$$b$x1kO>)Npu-#ss)Z^*oXL{!VGuzjNomDszqj*b$N;&Q1t48`Tbqt77+2Kbt;B$% z##H=`2ZH2e$60x;>|8=!7ZIYMFvpId0Nbh};*FqPF|J_A37-0MkeDYEj)*u5L5>Qw zlMa&{z`;_E3VaZZZzGL!V7(nRL0q9ahlwJ^esrUQesXD@ilC&i&*P1$@@+}$>%Fm# z3qH(TY$^^wx8(0Z^!kn{WSB2vUI0+jgs#OG3s3{0WEwNu#uTPHY?Ut?Bc%1iN<}RbzM)MgiPP|%9TDts(oqEwh+aVj2G z*u^-|D8~`X2(SdD=GbM>;&W2A`+YKG4i!=zO23X^9|tUAaRuYsq^cGC0JzH5rshkj zMuSY!HcsP=NUD@1hbnLk4{K2jgZkCXk%T6unO>VK7$w*qP_>Ymp^W8oYx#156*TH6 zsi+Yt+cQZ=-~8a{_O0i90i~{x<3|7jHB&I=rVqvn>0&8Ms&YsC5C$N|7=;{9gyiT4 z^mYgmA@&Ox4NIyRUkY6ykZc?dRX-k20ZSEjz5YB7-b-a`7~n-Bit&_Gt(Yfp7~dC% zDz5{2eHUtTj3}e>{Q%bjC@~tZm7sJh>JZFIyJoW1aRE=$f+(%6N)XwoMp-jqp_Ii> z1(ap-vOQLEbmUC$-r_ez-%tRmQS zhH7yVwq7ML;_cr%sbB)i1}K<{vG_fxFxa6gBjM+6C@+~jQB3f2jDBx5RRKnUAs0{uh*4s|kc-pVJZ&h& zbpeJHW{yEJZnI4?;6y?@8XswW?0|mU+X^BVQv(eUhA(y}15y?WVIyVI{#>8-Mn)-B zM;&t$EtR5X6kBap5h2kOH+ck*K>$#iDGeryA%A2z)lDX3$q_Rd%da0>f$|_InwG{x za5n*`49=|B3t-v7)Jc`#ftSc(l*Yq|9))`xbUZ7qD0Rw40zsPT|SWaJV*6V(*gI$na>BRH-1X4CUyQoFzow1*JDj$iTTX+B!IFLOaK< zv-PLqYJ@86Adzq}NhK4KGKDBzf>VJGg_(&(=?V$bN5r7>P$nGyVyy@rwe1`)@OrJG z&1B$pUF_gh%W8>u#11ubg^J1RqBzKiMp3!=evM)h#>la59*H7YEoA_aCdoLR4;njY zD}BB23dCkN5hYW274Ez|Qjjc)i%tqP(!TfZXv3Tn&sQ8Rg3|&u|z@C0Q_<3P5tj6qkt8O#v8zNm332m=0epR&d@27E3xf_BMmSM9@OI z(Oi_(Wct{#cQeO8N0mv)MAQI5wI~?|*Vi`nqx6!|%wed@QN_i|dEy*8IthoB*N8?4 zc3qd;-Qnnnv1<9EQmX}XV< zE`RmzhpY06E^YoNb#_=-*v}WMcrao^8`KR(i4R^fiZR$kf|VPo1tfZOBE$q$lSTSQ z%~F9jLBPt*PN7cLU;qH#V`>gG>eY2VwyjZU*fx+%htBbjp*q}ly*uc*kY*Lx(H{Lj z?R^PYQ(3YwF3@hMEQyeGqq46mki-S>iIQMggs>S_HDL)Mfp!N1Gd7}PLR3}>2x5W~ zFf4%(S~0#>NfbdKh!9AVZrcSFkX8g*{b$}g0YnycoO$!!_ulu<`SRVhZryY0)Tydd zRp*9m&6Qbm+ru%4Gr%8O1ze^_y3^?>_8CK6yub&~FXDFGgfUNMh7&HqVO3C6`WBQm z#+r|VNzo~ll)2P7$wfj-Jn$D3mZ`!aI_VO(Tmg|tgqoYXVma`+AhNUlcmvt`jFm+&0@JL*i___gp=P-#RI3+YE3^iWs>bFJ1bEbqR3AqQCBV$me;tmQ5vXGgCZyQX zMXPuu7gxv+=1kZPVKSc&o>@h~*WUvUMr~gH?zcIuBCA z@vboBLiBElxYCo5>Wy+QqHn`4bfhzk)8}rNl_Iy0+Ov#hE-V}eM4F)b;v^Y@h%rJ! z(5g6XAmw!~!_nxX3@_jzM{44&=*JK_0>eC9_!1ejiij_iAp#c?Q2LTeU{xG7!eaQ^ zkRqYF(_M@0lt2x^%|dn+#m6TQZl*z;OR&<$GU;?lJ%){N0lFa1UmM5@;7_jf>rMk+ z1``d)ftg3db4y*A0$)yODwYyO zU`CoFHGm2LHolB24b+iBQsEe8Awtapi%TU4fbB|UctMTaN|Cx-la5i!I7E^lPdhS~ z$uCQ#0GDxFlnxtA)?jQY{U{Eub`!&gVD6*^5^fNnVPrIb{%DmIkg2 zu2ZVd4_u8mFRG{8_}6Hs1wxZAYf>3#HL(Us3MP@Lbwpf15__JS0ToX4!2n0OD%!~h z1S}uSNs?F-3RAP-8ep2I?^pk!W&qO7%c70*bSm2SiiY9iYy|L+L=HbG7Kzn)cJU!xtIs!1EJ88Bv;+d z0RkWJ`fdZ@P@&(@_O~PhhaU#Gt7HNfp{J=|cTQp9#f2dI!c@YXPi`oblDJ^5&XpNp zxIndI_*&! zLQAwUQx#TEriJ>F>#L#tAr-gh$ZX4A|P0YaZ0=3L zfw~Av+I(LowAz9Q=A0l^B!vl;2A~ia@U0S=5cLxfjK-WukYnM<#M*Fhz@^OO;~!{T z!eD5)nUNqBe; z1-3v6WHg?!o8xXOAXqc$VkpDU3hu~#0Gv$iz$J(#0Uea#mBqmUc7v1kt-T^mix4;v zbu5wztrqxDw2fu56ic+bY1xHE%&JH$xCWmBLQfdES|#v&wPEQXM5Du|8%m^jA>tVD z9V7#rhet*ME>0JK&H;7}0eIq4fNxhz8(**CL#hm1?L&42&NZ$}JNVdwIGnl=cB3pw zI!B1oU)7OPgC;Nt;tYRNuB_X}7b=EHjGc&Rn3#yGnct{W145lbRFSMkmCw=8>B!S&K{(aYR|D8MjL@UuU`8hhgo6dSJYNXE z6yZY!e_jjRYET4@<81}85qm{pft-QBI2yj1WXmWBy**SSX4a%gpo&$+8Db zzD{>aw-~Uk6y)q=lB(6g#2c7+>cFe34krkcFT(`Dk6gBp%|@Bi>1X5w8H$NhCA7hD zne8}@tfWc~i2Vc3oR3b;E$raH31XN5n&aeWZcX$9Q4UHjo877=hHn9`EgdkQVq22B zTb&4*AU_F+q)0@dDlp1RQ6RQQ02xh|*@kpfJrPYL;-XYxQXvE42-n7RqfF`SqHcBI z)5i&L0f1*vbas1=P3|#ORRYV4eD@ALiOB;KUziog2fUIb0Qn;VDsY86lFcn4+HPW< zpYRgE`3M~`gsCG3ZZI)GL-fJHBoYCxFex3#*+}57IO+?N04W7@G?CPeLJ)vs{|twc zXWSl$t#1cnZnX*BT5w87XM) zFuew(*LmSDYk*KUaCEAf%6E%Fr-NiQWKmfvS^q}a98!0BI+%c5o(~CdcRB`j(@ZP1 z!j>_a5CMfs1R)~XQeaI|-`!GtC$H4mCWHXvfGQUo~gIbM2FI}GR~x|t&7N}%NV}U zBdI>wT)KqL2To7@_G|`Dh=jJ@Nx4D6m6}R;5PpCG-hyB)z)QG5itxci2f_)xK`I5d zIEZKRshWNgaZ2}gCWaHGgF+CU3NJzg$9W~zm`jFtjp&5T1~Meml{!}wkCGsbpq73& zlBCR3Kg&pz6_f{Tcs4_a$-@bJe4;?~mx*cfg*)YN2+5TJ9G{bAM%>bRs+u^7&McI2 ztx3XlQBbiH@L_s^93E^5>IN2Md~rZ#7D!x!4h-Yd?!K|PzOgx=m zVz#(QicyDJbHUWF10@FY08QWsL$T7qOcF1LldO!KNWk4mO0UQ9OHAd!H-NF$psICQ zdx0RF3vM@`=;?=Y=cCaXeu%VC;7zbDGt5P}N!e96jwgn$&JV|dfEp!23qzP?mKhvNbi(r? z%Y7(1tvIi6JfDL*O0qI$K#ywQt1^tV#BwMc{w`Wby@2>;hJEe26 z#t|_bl?FSrsn=4Qny_{@Qzp#$muzjVr$QL=Q^WU6OJp!SVk>%nYQv*b8=f?tT@!R= z-;CYuZ})~1ZEbHJXqq<$1psINz>DjrE$eIw4<}5iM>@8<+M1P%>*TcBhlhtxJI_#Q{>2U3HT6xHnMid17Lqj9jC#1u0~x|tyOJrdI=puVP=;Gm zlcc6BDxhgv*$c>Z3-~9d7w}1dZ@-@gr#u0mx4)eZ2+pjPmzT@}{vSrX_r`5NxBt`S ze93=AeTtz!a=HJk_6_zveS1A7Z)_ug68`H4GuNC&R%?2urfmX_%`JV z4sZ`F{rS>w%f|=)ZeS~+9DFrXS>yvF?-Q^85)BqN%`Q_`NGVr8eN6!?q^2wqQ7#1e zn)0zJD(lCp^|v=uny*|;@-<~jbLRc+7mPt105t=|npp7`{WTN-Fh@SNN4b>fZPcJc=!m{5~ z{Pf4SU~%96Ba`U$t|`w8s#ndiQPG;c33;z_7%=gN)nfqH8zp?Tk&0X9C!R3bVmIDL z$`PQV@=aVbSWIa5Nz5u*lfcT#EFDs=Iy@ArV!UN;C{#{qpiW2sdEl`prfzl|Ib1l| zzVDRPJ=cFn#qG%a$&K3fiPL^h^uQdu}{GlyDj1c<=0qh9+C43CNO zk2Y> z%ioWTjC`o`>^*2#;q0rQz49gSDU7C~9L)ep!tIuyUOl>fTF*ZAE{Rn5SfC(%evilDu2y8GfYVvu{ovnQYv%Z^Fkq zVD@GbYSKbiowHmJugRTpEwx$rRC3?R~8rFdZR4EuP2-S zrGqgY`cIK@;{xuRis~^2Ix_iEhb<0&vnQ3n$_`&_1F&SQnv*`GJHMtU0awxD6P7D` zG8izp=U4?k^isosid%N;bQwAW>6O_67SsaTl_8YJPGyfIMN_ssW>QK`kzDf|T95OPi+N{5Ya!f<68+*oWTviV^XFj$=b>*5>uw!Gd>(Ynx|j!P;rk>in|%_wPT^ zD*)yCmdWP)jf2TCJ`Oik+AG=@B~=^}$N1>nB+#4}^EyI-0Qv1RYG_7x zDy7T}E*d-i=Yr9CRxYqtnlT@i*g{9-)JK!=rv0=qUA|?C%GxyNQQ>5ZHmUYEP=$&y zuGF1TU6U;U?KmZr@be14n}TJk3d|_f+cd=}{A#ZIjs^HCb{5hlK)Fsrqi3D*9MRtAG1uH_J|`pWOsO5=UPh*BOtnPI@!J^*~x zSrH;iT>2N$kx54Q&^tyLM0gD#th6Z+OZv@(CjR zbB%+~Vgqju{iszw3#PPL0MCTG0w%KwPzua3x$~jD-+Z3$l;UZuj1hH~cv>8dr3#bB zXn+MuaVp6s)Y!Py4p09QDNbouuU_G9uDkbpGh>b%abv8m2a zlh-~ipa4ljg-2MMRA0%dwUfYK(_;W}AdOmk;%lt3#W7Uj?|-+U3svfAn~#?J6mh1a zvKH`K^Q2`ow<}#vWsj?L6<~0OPFnXWO7#~5`~(MQ&-xn1G~niMwsgjSAQP}5q(D9a zE)h;*Q%){Cj;OG#gxxbxkDHsbETow=Q1;D)i29(Q$5?8dul6_7s(a3Bx5U&=%=c!q zlTrcx@u%x!<=-v8(k7p1(!>B5KO5(g@k`K;bN&yl|MDL6#)ZEAvB~zYA0X$-W6ElO z=mSH3hUGDK`YADB1sJY-`s91Micqlv>g3kfE7K{X^J6V$=-f|P<;t#@0XE~M{(z|j zBr<4+J0Wsn1nmoSp4^>%c_s5C2&zZp?~OjO)RMzc@^kgd<9HpDbf?W!jx-y!7$T zL<=T%W7nr1PXLbd-p8lF(I>3)Ipe%^1RwxoaQd$noR8CN!ht#GS7nX9@cRSY>^(C1 zPCccc5ImrC8-DV+gg~lr=%-2eWh&4wAO{jt3}lP_B)H1&N|Jry_czfieV0lEN^hhR z&nvKL`H0_M`AMuVyxwcuJLT6e7r*+ep5Drh_ltL&4^Ol! z-jTma`S`Tm+huK|w>$oL^z7B}zpoI#O4#+uO<-9^-@&Ur(X|r??;Eg&1)d-K)@>H$ z-h_}M}5~?jF(+IUVpCq z#rc%O-rm!jCm_GD3;y%sX`}`z1j^X(j9Z8**1yd>ucb&Xz;KJRkM|lw!YnAE z5A;bv=54Irt@5BD;Sb+!MZ|ZNY`=Br&2>?25${U!@~i0%Cm-B>H>fvsQMaHoL*aaE z`-%4F0fqL8oDsLd*qc{+g(&BGw^yhGdF$@j)t-jO)in7#^qDZVee8Eyz6jzV7@Zbn zk~ZVkJNMS6H;ZBw34Pld9%oww?Xh29do6Sb-SB)pEfy)dU|Af0RR`C{)Dn$c8 zz{Mvz_nJjPT#8}@MSV7UD)sK#D7hVCo%fD%$W8mEO;zNkEPd^Pf(GvXJRevMYu(E( z*yy!5*?G%c*~oXgYl9*RI__;WUwBD&?d=09%+DvE_Z-zKiolgfUNJf6d(`YqvYZquh^W4A zuAaBCOV0r)Y~pP$UbnRrqpm%XT;_dZCoe+1nJcy zIqSvXnwIGJ6%Sf!MK98xTh(x^cihuS5B)urXuhCVk(lZl@4di__husgC{`GuX)SzfcP&apJmS5~ zBmSaO{1Eegr-vx^6}4+rn6=h>1KV3a_ovK9VTqd2!%XAMb9b4M+jlt(qOcC@Z4L6% z$sOhH!b7)3r!Bp;k!6E13*r(tHx@mgU+xw^nrUcakry1L7@&({LRvidl@$33OYo1H z)SAnO`WkwDHY!diMjW1uik=U=I2x%z3d=b`xs;HTD2TiEURuCu4~_H-Py8*ROx)i` zStDR|NlxC%yk}f|?5F zypEuoZzXvYMq@EF*vy-%Mrov`WbCgJ+4f09E7sdbp_YWzQIne$KZ-V64LpijR#q%i zoHq_}c5Se);N|YJYdP)ob|}N7rYnR%zAqZ|apOLh3(3QSdv#I*-h{ol|J1{TySJuK zaeO4S9*wKo8jR?qQ6KI1KTAo|tLzkxD86-D!8n$>Lifx8cGG6P+O;&Iv!%E%cXV^N z#5rXT#6dQ4bTQ> zHb>t425FOO+G6?p&zXN*y*gZpu)i(lX36l zNVS)uPLexfX2}fA_c|zvF9;6EzFkW@8{ZbxP~zjP|89-v&A{E~i7|E^E8`(Oni$zZ zJKdp&9TD4#_3GnKQ(y~+rTNUXH97mjD8U78S$*ZSpYC2Su#kpdpoNSg2IF4y1aDH_ zmxT5H^uF`j?NT}H^?+3x*+h-Gzh7^k*1+PlvhpYG=EIN&_snlrxws1@hgQHTluDOI zzaJY?Tzp-d#i)xTBct(ESJuV2dQx0MkW2Pl6%9Q4`65d4bK0^D!k+Mm;*z3a-Nfg; zx6e2qEc1T$rsXthR3z$j zhv#PQj|d^`6C|g-J|^Jni9B;tSy|;kZ)f$74yg~1%9~%gNK;Kkb{QKN*XHAUyf}kY z6D;dIOyaH!iQA>FaIrXTZ5b0Ea=?Qx2=&sg4Pw!pUbFCt%RBw|KqCFq=vskGdYto( zk4FWn)-COlZ)jr-;yJR+tCtRb|Ki0{(k|NGQl0j2!&j+SVuRj_`hB~aT`)tBQN!hl zQ@?!1h;06y9c{42R}SCH9S^=-HF?>-Dn+vS-&;ue3_@iC9Oq+WadV(IpjH@JgZ z=Ztqob}lR#@RvwWw(uS&^0Hc=hNSE4!LzobHJLnygu%dv?mx)1zBYWTfZ&vHGuw)} zKUX`@>3L*&V!4M=ZSBPr*M<4iNV7C-Ds#bzCY}+vJp21+JN0gy_qN$LG6jla9?{%>i0 z(3F10qz}I!j8TOmj zSNfM<+&Q_(D72={2i?9i=oWfs#GbAC#Q$u!diBV+fMiIqhrJHTF_`shUBq&#D)x+EdFs!Z^~p`I z@O6t?_5`(aga^yI7@K9F*LCc znd9~d5z+lRzJj&}S~w&cq}`2-cR8Nn5)%~gx)i?7^jUTvsr`ain#)Q1p3om?MrXsZ z5Ie&xlEGQi2QD{tC1*|3QmFQ)Fw`9Fk4g0I=)HJVz zKIq*YNk#|C!%iElI3HV{(E856=^)x^mFOJD3wek|bb1jfCq!nUMk@{r=ISL@Y(Uo= zr_{*w3j7n@FZRDZ5n)Nq5sV6PD=jq9eVg5`hByu`z3idS-lTI_}9*wlMuHRzEl^ zuR3+oqNU~;=BkIyh?(KGje)^lrw-lfm1GtSt#K7$63-Hd&EGXG)*se(X|q}0S{U^tqwsZ7>)&F&uDFokI=o@ zKF_0Ea?w#gF|4zo`{{4K0$m%AbE9OhXIW1O^sBMH0Yw#*k+o4x*v)6DHCH5U8-Iv1 zJsxVfTJ1?@0Rlsi$7S4WcXDuZd1JfNgVzgbY~0g>OOCH~Yj7Kj$FkGrRHF)w6WkDdO0w07yXZA;xkbF$(C zol^jPK*JJ3|c7AnY20RpNBK1pAwrUY3t;PJ{>=niF~a_Cg2R-*^f1t@Pg*9rTvzGMe9cwala~5&weu*!uwpw0C*?QC#tS>AzrQ0#U`BfG`MM*YK!?U)i&oZ*YO_7g|JK7V_r8b4|0=I2hfi;M~nOZ%5OvtR**a zTrE@ef24(2JsH-z*HPK{s(L!x55 zc6cwK5kx0F_Sr;4hW+L?5Px?tI+fInH60Sg$padw)mP8jgkN|T7r)J65u?G0{bPBi za~g$2p|*(^6xj?^mXzE|`++7Nj2*Fzan9(fqDCfOyL$P)S2eOv5@i%65%rxM)VpfE z@0hbee}sGAGGmV=<-OUDrFQJ$hlB1dB~jN3?!OR;UY^(EneaH7ExI+0W!K&}cHKf9 zOfUVBG%_T2jZYNCCO*PE3sa0$Ft@*c)76NG>Y5V)LEs~ESMnVMB3u2P zYn?<-k-Uwg!q>r`-j<7)A-8K3LR?}LTj9aRbF538wPnu(Qefu;{nc|Q-Cd>iy~i8V(_G~ZG{x)Q zGU|qs)wsGPzd?LgE)?)%oSAF+MovDPYK1GT%ncnDY3(cTPW`n29ns=(E$E3+EIl(h zQyP7yL}+x+IC}WWlb?BuwxNCwNa)vXxEh|JaELt`Z~s=+Erb@h{F25tl34D=p*^78 zHQc%^N)Z|VhRXALYkmB&Z(~f9d$-Hy?>jBojXG>0rF~ z_hyz^QsMX4u}xC1VVkWdYW*s6%giqnq(NAhSP$WhT&!D0hb&uGBwXWg^3AQyS2Nap z+gPP~vs|;j^5OBYmi*uR@!@*doyPY}9Xua+H}*yQnjS-y`K{NwH9zh=ZOwLaaYO#V zm4lkPceY3Hx7I8)Yz&^ax+Xo`39WtX>4k7ft>@JN5#>Em3vx)Db7JAK6^e!58YVAE z)jVy`qECMvRI^@GyO!FHF(a15p(qq#HdxdW8e9$ubiLE0O_ z2|hoE{s+Mb=ulExF%Dbojw3>){%qHA_1#bzS80AAb~Gu+B5#AM{#u zA7>L)ZFkOMD{;}>$4l7spHVm49Z$!rA4sd^pVtqur*$-0b%unzj)`4{5s7{e|6~1) z>(D;}AI2O5JGYsh(0`r(k5gbul&R+7<4MmtIG7rMs`+~N^1mSe|D^y34o{DDOvD=g z_Xc49RkQyt3hb3kk7xet?Ef4E#xp)DR*$87fGFqW@AdVacmKyV)qm}UKYzdBQR+W{ k;@<2%0C#{Bzn4ZNudrbFUxyPuebjvH$=8 diff --git a/man/nlsy97depression.Rd b/man/nlsy97depression.Rd index 87c9cc9..1687d07 100644 --- a/man/nlsy97depression.Rd +++ b/man/nlsy97depression.Rd @@ -7,6 +7,9 @@ \format{A data frame with 26952 rows and 7 variables: \describe{ \item{pid}{Unique ID of participant} + \item{sex}{Sex of participant} + \item{birth_m}{Birth Month} + \item{birth_y}{Birth Year} \item{occasion}{Measurement occasion} \item{nervous}{How often participant felt 'like a nervous person'} \item{calm}{How often participant felt 'calm and peaceful'} diff --git a/tests/testthat/test-refimpl.R b/tests/testthat/test-refimpl.R index 1d5eb8d..344a0a7 100644 --- a/tests/testthat/test-refimpl.R +++ b/tests/testthat/test-refimpl.R @@ -1,5 +1,10 @@ context('Reference implementation comparison') +nlsy97_1983_cohort <- function() { + data(nlsy97depression) + nlsy97depression[nlsy97depression$birth_y==1983,c('pid', 'occasion', 'nervous', 'calm', 'down', 'happy', 'depressed')] +} + mxmmod_ref <- function(df, do_fiml=F) { require(OpenMx) require(mxmmod) @@ -47,7 +52,7 @@ estabrook_ref <- function(df) { longWeight <- as.vector(t(solve(weight))) mmod2 <- mxModel("One Factor MMOD 2", - mxData(covData, type="cov", numObs=nrow(na.omit(rawData))), + mxData(covData, type="cov", numObs=1397), type="RAM", manifestVars=colnames(covData), latentVars=c("F0", "F1", "F2", latDeriv), @@ -95,17 +100,18 @@ estabrook_ref <- function(df) { } test_that('Same results as Estabrook 2015', { - data(nlsy97depression) - a <- mxmmod_ref(nlsy97depression) - b <- estabrook_ref(nlsy97depression) + df <- nlsy97_1983_cohort() + a <- mxmmod_ref(df) + b <- estabrook_ref(df) expect_equal(a$parameters$Estimate, b$parameters$Estimate, tolerance = .001) expect_equal(a$parameters$Std.Error, b$parameters$Std.Error, tolerance = .001) }) test_that('Floating point occasions', { - data(nlsy97depression) + df <- nlsy97_1983_cohort() + scale_val <- 2 - scaled_nlsy <- nlsy97depression + scaled_nlsy <- df scaled_nlsy$occasion <- scaled_nlsy$occasion / scale_val trans <- rep(1, 33) @@ -115,7 +121,7 @@ test_that('Floating point occasions', { trans[29:33] <- scale_val^4 a <- mxmmod_ref(scaled_nlsy) - b <- estabrook_ref(nlsy97depression) + b <- estabrook_ref(df) a$parameters$Estimate <- a$parameters$Estimate / trans a$parameters$Std.Error <- a$parameters$Std.Error / trans @@ -125,28 +131,27 @@ test_that('Floating point occasions', { }) test_that('FIML', { - data(nlsy97depression) - nlsy97depression <- na.omit(nlsy97depression) - a <- mxmmod_ref(nlsy97depression, do_fiml=T) - b <- estabrook_ref(nlsy97depression) - expect_equal(a$parameters$Estimate[1:33], b$parameters$Estimate, tolerance = .02) - expect_equal(a$parameters$Std.Error[1:33], b$parameters$Std.Error, tolerance = .001) + df <- na.omit(nlsy97_1983_cohort()) + a <- mxmmod_ref(df, do_fiml=T) + b <- estabrook_ref(df) + expect_equal(a$parameters$Estimate[1:33], b$parameters$Estimate, tolerance = .05) + expect_equal(a$parameters$Std.Error[1:33], b$parameters$Std.Error, tolerance = .01) }) test_that('Argument check timevar numeric', { - data(nlsy97depression) - nlsy97depression$occasion <- as.factor(nlsy97depression$occasion) - expect_error(mxmmod_ref(nlsy97depression), 'timevar must be a numeric type') + df <- nlsy97_1983_cohort() + df$occasion <- as.factor(df$occasion) + expect_error(mxmmod_ref(df), 'timevar must be a numeric type') }) test_that('Argument check timevar', { - data(nlsy97depression) - nlsy97depression$occasion <- NULL - expect_error(mxmmod_ref(nlsy97depression), "Column 'occasion' not found") + df <- nlsy97_1983_cohort() + df$occasion <- NULL + expect_error(mxmmod_ref(df), "Column 'occasion' not found") }) test_that('Argument check idvar', { - data(nlsy97depression) - nlsy97depression$pid <- NULL - expect_error(mxmmod_ref(nlsy97depression), "Column 'pid' not found") + df <- nlsy97_1983_cohort() + df$pid <- NULL + expect_error(mxmmod_ref(df), "Column 'pid' not found") }) diff --git a/vignettes/mmod_tutorial.Rmd b/vignettes/mmod_tutorial.Rmd index a1f0f4d..fc7dbd3 100644 --- a/vignettes/mmod_tutorial.Rmd +++ b/vignettes/mmod_tutorial.Rmd @@ -60,7 +60,7 @@ The five items are all on a 4-point Likert scale. Participants were asked how often they felt "like a nervous person", "calm and peaceful", "down or blue", "like a happy person", and "depressed" in the last month. These example data are included in the `mxmmod` package: ```{r} -data("nlsy97depression") +data(nlsy97depression) summary(nlsy97depression) ```