From d9e072f1f2a5848c201bbbe0c840c8e5bb42fef3 Mon Sep 17 00:00:00 2001 From: noah Date: Fri, 30 Sep 2022 12:12:24 -0500 Subject: [PATCH] Part D --- .~lock.Schrick-Noah_CS-6643_Lab-4.docx# | 2 +- Schrick-Noah_CS-6643_Lab-4.R | 49 ++++++++++++++++++++++++ Schrick-Noah_CS-6643_Lab-4.docx | Bin 27599 -> 40441 bytes 3 files changed, 50 insertions(+), 1 deletion(-) diff --git a/.~lock.Schrick-Noah_CS-6643_Lab-4.docx# b/.~lock.Schrick-Noah_CS-6643_Lab-4.docx# index dc2814b..20caa52 100644 --- a/.~lock.Schrick-Noah_CS-6643_Lab-4.docx# +++ b/.~lock.Schrick-Noah_CS-6643_Lab-4.docx# @@ -1 +1 @@ -,noah,NovaArchSys,30.09.2022 11:57,file:///home/noah/.config/libreoffice/4; \ No newline at end of file +,noah,NovaArchSys,30.09.2022 12:12,file:///home/noah/.config/libreoffice/4; \ No newline at end of file diff --git a/Schrick-Noah_CS-6643_Lab-4.R b/Schrick-Noah_CS-6643_Lab-4.R index 5038cf7..51cab6c 100644 --- a/Schrick-Noah_CS-6643_Lab-4.R +++ b/Schrick-Noah_CS-6643_Lab-4.R @@ -81,3 +81,52 @@ p <- ggplot(mygene.data.df, aes(x=phenotype, y=gene, fill=phenotype)) + stat_boxplot(geom ='errorbar') + geom_boxplot() p <- p + xlab("MDD versus HC") + ylab(mygene) p + + +#### Part D: t-test for all filtered genes + +# Put it all together into a function to run in loop. +# First write a function that computes t-test for one gene. +# i is the data row for the gene +ttest_fn <- function(i){ + mygene <- rownames(GxS.covfilter)[i] + t.result <- t.test(GxS.covfilter[i,] ~ pheno.factor) + tstat <- t.result$statistic + pval <- t.result$p.value + # return vector of three things for each gene + c(mygene, tstat, pval) +} + +# Testing on the second gene +ttest_fn(2) + +## Testing on the rest: +# initialize an empty matrix to store the results +ttest_allgene.mat <- matrix(0,nrow=nrow(GxS.covfilter), ncol=3) +# run analysis on all gene rows +for (i in 1:nrow(GxS.covfilter)){ + ttest_allgene.mat[i,] <- ttest_fn(i) +} +# convert matrix to data frame and colnames +ttest_allgene.df <- data.frame(ttest_allgene.mat) +colnames(ttest_allgene.df) <- c("gene ", "t.stat", "p.val") + +# sort based on p-value +ttest_allgene.sorted <- ttest_allgene.df %>% + mutate_at("p.val", as.character) %>% + mutate_at("p.val", as.numeric) %>% + arrange(p.val) # sort +ttest_allgene.sorted[1:10,] # look at top 10 + +## Plot the result of the top gene +# create data frame for gene +myrow <- which(ttest_allgene.df$gene==ttest_allgene.sorted[1,1]) +mygene<-rownames(GxS.covfilter)[myrow] +mytopgene.data.df <- data.frame(mygene=GxS.covfilter[myrow,], + phenotype=pheno.factor) + +# boxplot +p <- ggplot(mytopgene.data.df, aes(x=phenotype, y=mygene, fill=phenotype)) + + stat_boxplot(geom ='errorbar') + geom_boxplot() +p <- p + xlab("MDD versus HC") + ylab(mygene) +p diff --git a/Schrick-Noah_CS-6643_Lab-4.docx b/Schrick-Noah_CS-6643_Lab-4.docx index 0c5c630d6e7ea73943c7cfb93413497e6d84d3f4..78c7611ea2b7f37041b9421c53929508d8b54f87 100644 GIT binary patch delta 23565 zcmZ6yQ*>ZW8*Urh?ATVv>DabybZm97V%ttSM#r|(v2ELSc7J>P=i=;jQ}^?!^^U4% z%$iyPeQ5_nP?Q0OzyN`PfdL5)uSFmN1q%*OWB_FX)(QzS^0av@^YjeSW+MdGg->n$ zIik_7&bAU8tZ}}R92j;#-P%y)$RfL9p7&_3ltvApl7UjsBq>+2-s+C`%+vB1qTxvnq37a+BeSW-D6_LRs07l_XituqIB%0Ac$#k--Su4W9Vb14G&vHeaG-LJ7Kkf#n{7_j$i^XM> z_Ht4r`hy?z_QAgzt6g>`P+T4P*Nn%CIvJn^CC|;ml?n!pQG}W#z8HNNA3MJ6N=xK5 zf`wnP{bqoD!zGuHRND;u;kt$A=4f)84mw~6kPsFsJKV4O-`1i1`L0~>Hq$`?fG}3K zXxnex9p&{MnbF}8f7N0(w>Z$%Bq&bCOzwd#K569Ed_l{;`;n;+NS!Kpa~RrBt8>&O z`l8>NuB#I`Z8dFl0e(Gp#Uzf*(tq*-->#>A9U3UgK*PLtG_1LRf`FWYfq?wqM4_Z=bR>J3wEFuB<59|R%1Ir*O=+Ws4gDYzMBa{1)M&ri-L>JJ@R=TD79 zN0E%xu9&(#++_=-viV#Hg0Om$uG&Fs@2+XV{&%{V9kXy)JBYe1rpPWPji6@~Eg-O? z41Kd+tW*QxXjYe`4Vc^#w$_)-ja2!u5;WUlK%_I`XrdB+2}Bj4D(K*}OHU2OB-3a2 zv~Z+3V3)duaQkiIHwWfZ&t}FXEp0^1FxY^rBD~{ChYgk34&noT$CtJ1zA=2T>`a#W?qi!M%CvAy=kn zVD301x2G0jfg|6`79L|3R@zopgQ8j2!^Wuo+)-Rg;W*vvWDL1qtY4;%?PBSH7vJ(_ z#7V_ikHZ$q{fFMks)MVE4#m~i;1-w-~S| z|3{u*4dj1cVi3gmf4_w?Br_ln#*E^3#XDH+u7iacZnTq6iQnFJV6q5gEB#w+{E`?t(1ABQ?=ahk+xDEG7^j2OxMstvM zd)xr)Hx_mDRtXYerNDT99q4Zwv}e(U-+obK(PJt5f4hGZ*;zy2co_s_&CQ11|Hx{w zY;ixsuc;?9ZjWHefJ;(IhY7Ah36o1#`>`~|dhUOBY{!uys6haqX*+5wS<4 z=ZN5Yhrx}YePCrQHp z9)OCs6Y}@-BNjtZjNeScf)rL0{)>W0(J&20Vm7(w?Y#8PYJFJM%DUCby40x4xPQIR zdvyLm=VhPIl##o)#n6QpKS)AF0E9t&>P$8|uP3Boc9xMM9V-%NL`sKk1i7J$PGvwtP6D@q~dkguq zmE=j)1v3SyX8Smq{E*GNfTsP(la_qI%> zM&DRp!~>p!ws4Ih_0x+i$A{Baif}?YXsDlc27%cjb8}(c_ty5Q@0ufs*-k5OOu6PT z0Kvf|_{PE6K_GGM*G|)CB^K>+Da9x;!fv5UVE)bh&GcZ7Mu_GQ%kkyaz#dI7$~h0#B|>)sN|fH0eDKnledjThtuZ@H|*- z%HB`29LMN~J^aE4T&pr?fAua|#l z<%(__b(q=tcyms2;gn!@W_2!@5%zDs)QO~lr*w&3og-a3y-fu7;0zXhq6oRWGEdg+ zTiMM3zUIJGXB(q3rX{u)Q7#kL2mg7oKs-1P+3WqP)}yIhJXqu4Ep1a=zinZvlO0px zqs_O8&cDbuECL&Qc8t-Qbidp`?f{=gZPez%O7?QKxlP+?8v9lcxYHPwS3Fq0`Pr_^ z)NhZ)q2?XjGx`s1m^2ad0~OxC|EMV1pCn;9m#wJ2^=sC z|7!xlxx9RE#aSovt%9}trmg1s=ZB0fpyBRNm+RazNfe&OC84vwi*R4)N0CzD1vfb; z9thn|^!0QX^VNC{SdvtY7LT(#72)|Gk55-;`z9yBdw99>DGxYO>C-NMI!~fTf55s0vorLO!5P)1#8k2b;Ha7~B zToclLdc4560niZ1S6dQ-zP3KGsZ}2wt`n2(yu@x!27R@&ZOB#N3;mkuFiN1EMG7qp zxMzw1>KxHWeQ_EBMUK0rzQ^1|!-5n^+RCf4Kv&o&5!|D!uol_j@Im#D%cxm(%1KkL)(Kg zkTh4W0D`LPChOCK0R9I)jsB6$k`e^b+Rt1 zO2=1E24Wvh6hH5ZE`liGMJZ*<_J)$#EhPtiB?_cQr5oISNbt|CgY4Ek6|bqa!{Y*t zuc~n$R>I%&u1WF~Dg5IrP#b229JXIinaW>DKyVONUiXeycsvJTpQXoIzAGr%qaRmy zZ9j_wmEkXB1@yQ8GdI`r((*INJSXz)LA2LnKYbC%YVkuRv;H9%^(Z!Emos~I zx!=iUcD=Pb5qi-DoSF5tYx~C+j#zRFrV{#bpns>~6ldc%XIkjS%Z9GFh2Vb9M7rk? z0NH4ZU*J{>UrSwzQ^9^PuLBi$g~c}XsaXqefBGy}+An!8rK4&qS8eF$Px;kZJD(Uk zw9)hl)R)!{&*!p10|Td+JRchu_)6wm7DRI_N6%gL+-bAj_#a#uf>^+PF*bX$fr5fq zbSzSk_~Nil9D{F?1mDs2R~crPZU5%60~fx{8Z$fN$1^JN!aWpIOs;$-)xnjt<%6WAhjP)gR7z;L|%& zsN4M0{My$-IfV&SsKV>vRPAUa>^hC56Z;`g1>#pGLkgnlM%L z1*>w>+N`p0P9u-08ce$}fgs2xfN5L37m;{ia9lg9jb{p_NpkQ~2}%PNNx}xi-YQ)x zQ;o*P)k<;|;wH0#ffy*v@`Lf3q)5IAAAaORhDu8DRGxM2=DtOkbrj@16DW~CB$CqJ zMd)Jzc{w^X9}HV(Dn@qeZ1lwb_Yz6^*Py&h;pE6gVN9GBPO^JpKA$X3V1(VBzBpgaa;c|sQFz^rmTF_vOT-y^1B zp&3L2z`wc0pNx81?=@2(z=Vds{rkQL6fbry9EJn#VQ(|0STH|O5~4uB*|6f3>Rq^e zx@=yefh12-ove-ZKGWisl^_^~Wv3$57>Z2OW@GN2abs{vhqlG8N7?GrNarqjpr#pW zCo0r$8X1ivDqe^iKpod2-DK)LYmLO!XZxMn%HGR7jhqveGBq{}?b`|4^g#ty?DflU z(N|*nnJNOyn#=E-?rh;y2vsbxYvDHBetwanofKdJf)g z;0lFi81kXHIDOZ3c>Y96r?b@g`e)%2=c+|tL6a|ev3d;+(99(|C(hR*7EILJb6RH+ zMneBoe%!A)4nW7};M1}YOcGq1&Vo_#dFzr}7xAc{QeA1cD8IV#q7aJ1!1lBXz$Pd| zbxgdxkdiR(n2U_$xwJpJ=$DJYf}BLaSnE9J@G3+97X6bqm0FIJiWJ|%a-G*3xEcgW zzi+D>PgYC=+M0QE!wsU^5U4gQ&l=YChNqX^v4t_!w-Ad^HyHeOE zm)ehw+g^}wBt&XA|FjTTXh2ZvL8=ht0Y!Mem-QKTqv5E$N*=$ugDelG;)VAv2L0kc zw%eSi`*ZJU_3foXz{3|a@qW)g@8wOEBp=^=b1Zxb2=B~Vw|_&Y$PChUo3HnzqRZ=E z-CMFIM}tL^e4=Oo&HCt?9kAb=kk$dCFtMGIx$D_CVK^3wpjYs zC<}g!S_^XaH4W>0iEO=?wBEher=HOWdhd_8N1F##BNO~MT};l=-nnG=0UO|`75IN!7H|FO@bI9CBv|T?0L0TU92>KcDB8> z=Y$QORk~G#jLk<|9CAxUK*|IKIuj))&n?CQk)~wu;x=U}-vnm)onkvwZfB0#qSs2VWmR&U3B*b|9rkYwCzQGm zH}IeL4wC!eS-N<1r8L5-*{m9I>}~lliC#+%r`T26;b7qr-;Atce%fca2P68cv2(5f z(ME3?A1zaF#zSA$-iBdqU@=i%+Jp530RtDO@(tJ;F?&h*{P-9z)^{!35rpoQ%a3pB zn!iMvF?vqZ7ft34{Ai1_7yd{nyf*NzTL03e=MFIfCOQk+PHVDEDs;s=OY0;SdQmYT z9)7|u3Ao;G72u!7W!c7g)AA)c+yyh=VKKl`ti&Lg(=+=UeP@#L~1!F)gOO|AZn&J9JYAMjd1 z@LzwQ{qfT!bP(3AwB)KC!rx7wSqa#{g*xd1cl-GoAVPu$vt*I#y?`~=hG^~r3>gko zITkrz#Ce~K@RY*{S?x4@Q}t%_XnlTovgC4!aa)3Y4=W1(ynA@m<_KKrW$EE2pzht}d3tW<@K8}*07Z~&I6wmVsL&$+a4uMyh2uH2>QNs+AW^jvBp=2Q8 zD0B(Th>s~z)>B*PDYWwrmW>zy?ITVH`!6f`%AhTIjb6Bzvwt+1m+!Isa?3-U=lapV zaj@d8y`_?NF4)|*J`^T!$@N=rUkp}gEZLxWkicNO>mXL=N_Y8C8J^t5yyJv;L(zKq zgi3rFg)RsPGLHMLWAuQi+JeO~R~&J_T_DOnp57%}-H{~+G)5x+He>$_oY>W&LHtpT zb6t?t*Ut5OrYM%}mrk|oO3Faat3vOEFOccUu>3CXD71*0Hi1;q#Xi9k%D_TS07~;4 z-F96cs_E@#N`W&+ZdM}4%Qytoal5VMh6Q1Gs4&(q+JqhuZgANHdsA;FTQFH3-VZ|_ z7oU2usC-a-n=9|42Vv4^;MOf{-wgGt^WF%f=ID1%?zR?(Y6iNq@_^?+dqS?Z?ull< z7SYWF)3)D+Zj6ri68-e=t;NTN&u!$2kb@&`EE#<=%+3B+#>P9$4)j04+n*Kdf+hn| zFgu}#I?f!SozfG3a;Wa5DYrslAEVulUv*D+bToa#u`%2DkA$2GfUI1a(0GrrB2f66 z4E~M1N9~>(4Hih&{uDXIyMuw^cHQ+|J1LZVB^gTNeyzyP{~~-e49)Z0mHA?{J<9nr zM~HF3XfgqAJ@i==^C0&-@7Vy2CUcsoiyn zlJ!N!Ab6CNRi%CixyN-Sl}G;ebAa@ZfiWI`mvkNb9$rbYl_udRqQ%-o(w&2+{DYys zn;fFu^t`s-q?}#ou3sRm;*P9$1_?ZItQHRE)6%d}BDSZ3SeE%H(C$cI-O^ zqbpaf!l^AVk?a}88v?6`Gk_o7Z?zvL3M;uXMBFUICud|WCjb892n{j`FWo)l0(mG* z)4nR%xKO&B(&RF?zRoKziHyRDbs`G-Rp}AcW7&UAKB4!r_7)!nJGvt@YR|3DEPV4Y ze|@eqq0*lU3$Mvrkq$zxU)1uWEAfW?5R4euuJZxp17a-La`1PNHT-7q^Q`H^0MX1% z_cK|~I>NLX_w%(0{1!rqKCN+s*9 zN_s$iSDPEDDgAS0k>m#w@emndPi(r`t8P~Ff@mrX0=GlwebGW-&`vC&FmB)yQA7Bq zhuee(;#U;w3Y&(RW!B}fOVOp%U*U`VhE}W-5(GJ!^GdOal`SpACAkJa;lfiQH(50| z1PTj26s}brKN$$A4>kToCJfQ9t|?c@NFzX{jgWQdah@AweFxms?${39PbZF1DA2IE zo*qWamiJT;io6av3^Cy0;5g|&Uj-Ee=P_zkrDE_3e^8|F9SYor!&{Mj`6e{R|mFaC$Il%}V|_ zeLX$G=I?jyv%f>h{g%zSjKd(krz8o}tF#Sx!~vQZxyy5qTM;B02?!S8-I?{ULRY6C zTs@~->9gshQ(TbuYM?$U4W54}r{@7$-Zl}}gb>xDrY58hD|>t7Wp|UL7RmlGu8nV= zP)}={?c31e*?!vRf8V?@t6u{DK7!}H^M)68*-(A`8-57%UUxwfJ&YzrQv&t(YV}UP z?ka3z*Y0HR&+)oJ@tSoX+>)702ZKYNAZtWJ9}>Wu7`?4OlEsSsLvI_DxDtVo&@l^& zg|rimWooB@%YgEr)C1p`xw6ycG8Kkm`zd@&#HY!!S$amX@6kyR{s4IgIw{Aqx^*;a%z)2&XjNEq1|Uc?YG3=1{j_*IJ*y0JPwWt=E|c-md)cX zIX& zre%Z6s;;KAVE8L7Kgj^HkiAx%Qfs&6X=%=t<^m*&9OvxL5?(OPJzhJAeiKUHunJLr zlNCnmGyL#74Mp9pn}w{H)BaF>LXcKVa-;Jn?oQM~tn&e)XY4gY#4;e}z{njwqui&i zD@Eq+B(xtQZ}USQT_&4DGcs-5_Bh;uTR@P3swMI8m9LUT>VkuJ{fCvX)QL{7yg1* z`~LC7L5kZ?YZ+^?0S;i0E+{eRr6he?CbCWc+dE@^dG-55@wrhB1CP4m11lh`^7kb_ zJfXcj3?PVD;J1H)N+S*-CUmZavH{p9G)q*ErXy`p~2C6KZoaKEVhGr595u#fvj}dtAUW| z{!FeAFUBpQOW`uVnwVj1Ro>3jy%|NxTteY|(3;Wt!_W@s98O}dcpx?%#d-D7Nwukk z_f*qb+U{K62WEBchkS7Nz{m?cgSTz23O52R+09ZCi3eKI-;14wnjovE5N{;XK!3Eo zmHgzX6xXE)#NQ3p?TtZdm{I=TV+ftXl_bh`4$Wm$uKKaWD-Ai^;HN_)M!_y|#Q&H{ z=d&O1i(3cSEsCkrpw*O;;MSyVf|xY69*FFzTrBspxL|!eR4p)29oxL_kEO?Ru<=9D zgt(>8+Z`&6Y>GYL-|G3fu!((H11DIlCWf^pKD*z)&+Mxl+`BOcK1fD&mW(ff4RU@N z!f;!IgHv2U+^!8tlvtB9$AfV~m+m0m;cMRI70L#jzJdovx*o3jWsjf#716R5wYl~P zd*@}~!=k%xD=pI4Y9jd;&Re#>3etR^@qOE8WBm#cfvT zPW`^HrQKl9YiykZL7m;rpJcabEl`{ARy9`zFKO40Rcj!P#9n=p+@B-n+fVoszYeAi z#K-~9>~~M8xuaFIELAc%;JME(LW;U_Lb1%*tb4X;zdJsfA1+vLjxc?xdHj29NLMeh zchS19gjUsl8{3ZY<e?LbW)4^=3coyC^VH@|puA0`ACxoObjy)OTi0V!;}3|C;j7 z zwd}YF3Vihc)i(+MxwLW1#kkIwD(Vx>9Q;cOjK%&SJs9+MJ6DSK3pD2I>kSA#DuV4X zx#_PLUEDm_n%YJ$l-2UEJDwXcA+&z$@i!@QbyvKv-99|o#BX6{tN80n{Illn zNp@!Z%;KnVQbbYcW~@39kL5fQdUM?Q&Zjr zA==GxZtW*2m6)LLUL~B{Zgq;PxEdg1N@R0xOJwMm z)F*x$0TcKW*$;ltJi~-Y=vM<`s#o~3`YYF;oE9&E6E37)?F$2AMVZ?91nr`QRuVN$ z<~YMQg~9NDr`x_{9W-Jy6;EFZ*3N$EJ4o>BthqSrjMsDV)sJ@smz@S1<(rU=L~B0_ zAa=}l?Q4}@j!fsz#w%d$#=cIQe(C8-PGNCHvKLC$6}>aFKk^gdnT=?z6rIqg`pG?h zZP%0N`WoBL4POn{Kl9h1Xl8^8V?AFrLdlzxoC;fqV_6jldkqHH1n7^0Ovx4^Xs1V* zUq{VZnVqh-yZXo?*yVRLsNRYEChQTCUdw4IVUFr3Sr+GS%mjdTqK;&(b5C4m`C8yv z*L^RaP-C*WS!jbLc+#A^H$|%vJ2oez_eCGiwnXjDnaE4hQ!+Yg+oPWY@1c<74gib7H-C_c~RomHCR_yEF>ct+~Lz2Q5WH6&M-91@WXAT?Y%`^HH_|a z*UQ*9OUL8gbeq~lhqbGXFa=w3#W&G;K& z>>-Vd2yu#FHKIQ7bplgS&5nuNYY@m*>QT>GyPey)U zt9=7JHh^-DxZ`9jQRIWZg>K>Z`Zp}TV}Nf4B`2Uq6h*W2;LPZ5<$7EihnUf6rXm&fuZ3$N$6P8v7)(q;gk)RYMorO9Cfxwi(UNk_bji>)=9dWJ0PY( z&#Oi#>nR#kG;0p5PR_2(ytdAulZx#xODLAX=D+X_E@d!nth*Z$YL zcst_Qounuv%zd=Wsums!mBt%qr*XQeC_Y+KE@CHcFA5|I*S7NXpc`^L`^ zbZvj>SWzEIwiYFJ#@9wpJc^P-=-nsD1eN^onJh2QuDdGDZr$t|XX-qE@=B#sOR(}`(m-wjm(!T7OEd1 ziShgQ|1xj4e>l@QV0ZgtUA?u1S>XMCQOw`u-+mV8NgkGDjsCrafr~$kWY2~WM2FaS zG#2iI-XUqsJ^G{MtcY>375G7m_C7q~FsV9(q~ovO6+<*tiSJIr(XGvPk)W_;(z9|t z1V@@ibx1W$wwWAX^bKo4QgJ#1TB`%0BVvMCm$D@&r>{b+>3KKS=SYfu9&Ju6GL`$Z z+ZxRQiZqQ}C|s>GCoxvH($mm1NRMZd-28SaReukHz`C;iRwZwzMqQttNM_%Jny^ve;n3F0PeU9Y=}hs< zwbVYi?f#^I)1!%~dj6EicEDlV(cHt-5a3+*DE-0x&yjO?jQF+3{V4$^9-?qBy7-(F ze=0Fz!@QFW%`*N66d0O9?8V9BGe}v-AS0}}CXT&c9->x9Ae=8~gdyl{J>GGD1wqk5 z#;-tZ`mW|FggmlZm`oa@;wRXc8Qle+R0gsRlzVo@$t4U7*!L!4MkE(dfgeI`2_zg` zULSML?JwcO*UV=y44;@N`ce8YB<^3Wrr(n1@0sXva+08S04Y?G#@4AG5B&IDEygWz z!hxAK<5KZM+1p6?mj@2rTz|8=2W;ZRMT2^C8lEt+dp4E|8ik_oHm%A z5)O)1zL4{>u>1jrB`!e@53tq&u#|n%W%E##K_)478PkSZ{c`inZUt?%@^T&h(j9vc$Ivdsfwh^{+giOH1Ih^i9 z0HBm_A6aY}Qy3^ZJNdOfQTg1%H#*E8wP;mNrJIhUauVL*=JL)s2tQCQW(L~Bc{GWE z`Cxao`VxV&vn$R~z}r9ROG19O^pzz1Ydw~bZ`HI4;OWlo^p~KX5kHB>UsV1uZn@=U z!Ov2&n|WI_UbO$KbotUzs7U5vbg%a|lOH-iD<0U2Odx$)gv6wKZrp@rm3#zpz{T|S z&VA=XlEOSsa3dT}R{8#zK#HgIaaU$Jcc@w%W$ivBL-EQbs3VY9e&?YPTM(O=_hlWmu0frzYaP_>J8XB9L?66LC?ZH{E**_^W?ID3Skj->YHzJ6}1|jwtT#{0W>~*l_4PJE%qy;)V_7mnZUBxpk?x;51F?MUO$NO)8 z*U(>4d+q3RSZ1~PcMI|fWh&%#W2kk4Gw6?}y1^;8JIvSiUFv17d+ICtfSC#C>$@9C zp+KdG8p{^GJ@md7od@@g4vV8eO9k4<~sE^@W(1%o%S>GF5=Gw)k!H%LmH#ou4xjk^tE_d zqICyrvFU$}x33DHrG-KLd^!Ga;h*+ie-vU=F1#4eoEAytL}r9=J6fwfL%*%x>7NeO zs9}B(nQZDR#byNti26$pRu z^yyeckb08-s!sj6{_<<5?K7kw)KHw}C-OnlwVe1|Z68 z3>#sZrsBg0@Ukp}jsBGYQZ3HZhEET^nhOGD;2h5sAGxuJZ{HG?kbUMXk-!I< zeyf4Zir(X2KNb0hF>3J2t6&z<{+gVc{Q9H{V0|-X3Ts*YNkOe%~Q)2^GkMIw(?UKHa%sSPkk`jA< zYeUaGgz)}T#!-X^;bk$ZKZMJ^!mhDv`cp*Dvg3WBNg-y-3LEdd6!nB&B5#&_-hb}ux zYCtJ9Zx4~Bd8r9y+WW+fGd?NN_*6aMBt<6Os7sD+GNwH_5(GVR% zSPbHQ-<1OQe%1+afWGuk3>sAF86rsSZa66+HP+iMi|8=uD?&t20M?uzUq2H}FCdv| zny2&kRTvKRB(1WKQ!H_GK|%J*87PcV*>I=r%}1bTEe8Ubnyd>hxbaNWP13s^JuGWe zR|xu3NNCqzf<#$T{eJkOgk(?4jINlL9^)O;skS0YedG<= zN>C;4*+w`$90qNIvu~xkBE6muJCp~kiqHh4$!m*Sdqs1rkNc$(8?}F>InK;=QC%yD z1?2D_&`?Rrf$esYpHq2RNz%+d^@kpR-+HY4CYiobca>!%VE4{d*GRLS%}|Q{y;i2cef*7j znEwEuKRd&yvzx*3Pj9GeTrd}@P`_4*n(p3IivO{m_1&Gv{>C{vNN77YJpN@RfpzKX z3;REL`d~ipzL6dTgox+=k*DM2dyxLg)Ah=FYEb{JWhUAG_esYF6n3J48Lr^UG@?*; zJ`aRWtsQt)zQs^JUSk4vm<#LE*3t!B(X7oB^O?dUvdy>cMak(#SIl*0^t<373)wOeFR>mYUGj|$n=Sd-;E zTPQ39x`;lLX=CWlRCXKO--Qhl*Nrc{z@Z@1y{)896!GarTI%{9=(MsVF;mImwuC=w z05_lL#*imfKf4xN#%cZAUke9Y%3}FqX_~S#-N*tJXN+vXhdbGuL6rORoW2}AHX4EO zi4yNADD|I(N-s|rL%BKK`-MA=T6W8%xAaT0ML3TCEMrSzd{!YL+5KeoK=#F}q3!eu z?F`2rlegNOXbI<^j6*KGv#-^m)XeK6o9Gf<$EQ9 zXQB4eI(b9P!|(<6Ukcc;Io|&N6!t{$|9=XJFc#E+eE)Y8v8RuSVV$Iw3ow$T26#f5 zHt;M&VWDO;Nj44WKZT=ZO!r#fSM~lXszN zXYaZH;yQr(Ao!)eG=r#fxXsrEzue1_pElp=N{x~*bBGpImjmw_Hu7wWPQ=0fJT_JE za6us21C73Vb$m{Wo2hRa?)Ss1LZRieyG6F3UgQ6c5>nxz}K!h|NmYn75@pkBG z*iULuY?^U?bOFrBpG=qVel<9kU+HrKR(eV?$&mt>^|d1T;Xco z{MFhSPl`0R?ro90_^Rp~E+R@FM0nv4D?zq#?O0?r@nFRsxzc}|nU|}5=8H35kKlc@ zUxy7N#X}fHUJBcLh&o(*4F;tu(L(l)_>fCi22>Tq{{@DxO$9sMcEE*?FlnC;-Kl&| zdwL9Kd6N(8`L7=Mo{5Sq!M=4+2A{TH*NpLb`~v+d|GqmH#VhX0pG;mLt&aR7-t%bl z1wkPnn+W0}(cR4H;-132E6zvRns5nOJ2Z`!#ZoT3!X!@pHfB^_eUIQG?oc&8@q|hx zWyB0cY6loXy2V98inqV^YyEm7#wU(Dm&5Vd^K;+4JJz5^vib2V zaqx@B!slK|GVhkp{n``gMXk7?xwT86t{`05cZ!3;2@XAYw~-Xvif5}Z_-_nN*+=)u z-~pWJ-x8SLiH5q~AC3STu$Sl`EJ{Ai2!2Ss%)BJP00sy~=76b`*swp0^)Fi`zn-e#-KIDNMFqR=xSrZR zEJ!#|$a)5^zQV-r1C9NS&EjAM361PlU_Wjpf>t{=$oBC_)2H5%l1dD59c+pIY(Zmn zczX74L*rO54v-o8tz4`y$rqbFNpNQrel8L=UZq$&Toe5BFH7vPF3|3sn7&fY@HlmFJB*WiW0NT zy*A@jVrwXA*}z5=tDZ~ z=3>O-BU(f)#2q{O_`B6v`oq2`GxD>{%@xJO{Zp73Vog0q=< z30~{;YFq_l=4I^6cZ0;aDQmw?sOzYYy9Sz>B+*5Zr>R^bP)l2x3_EfDV3*Pw9kk-# zm|fdDIn357W@5&KL!XXe`9?-{Nj~Vneq0~4mPKpV%;rA8sl4r+H!{gqG%auT&Z98K zd39yos7{9-3)xv}pet@T1=uj^LS5u3jg6Q;Q!bjYle%uGMQg5C&Cc^`C4dBlF z1)#cDd>s}trlnr$t(JEaT2LowsfcL%Agtc3CS*_^xo8JuVSe<}gJFg~F3OUFqS@9(38R|3DR2(nK1~m0H!RP65 zW~~E12!5D6Oz$%|NvH__CR)NX7~Wxe3K>yY76=x!UIXA08*R_es91!ALeWBE0QJP7 z9FYa_L4KJ@A#GG{dJ8Q-{Q&O4m7YmK&%|LBlTa`KQ`B!1*vG~XYCFIO3=+=besn~- zm2SecULR^Fo~RhyOH$E zi$4*ZPa1z2p8G5KhG;O$FpRN*YnDo-SW!g2FGi2JDWjb9VMLo?h?BjZhDt+*Uq?wKC5FmAj2K=5F)OWjhU^9?u9bB?+KR3A^VmQXm_#r$5aU{g{|=T|o_m zc(pllzw~&1cpUE{NVecpf9HGXCz9kyno)0UOU%$d;gmUjAwF^H<`P?(vZ93|fVTS3Vgjvb2EYzu$wsa*9mtE@5O$(377wYzeiMQEQ_fi=x`jFsn|C}L?DV`*hV zfdk6Kekp{ZaE26ydHK|nEY1mDtT-`K5|Mxw+V(4!8w~;ihh-78*nN_JoD>_jP%jKCoFsm}4I9R-j1`wj zNhU>F#-O2Cg4!5K?5oPzpxAb7El16aQ;C=pNkR{GVyo9?L$;KVx$>^Rm=kk?PXi#o z&qqi;`SKV=4m;RmiUER|5jZ$He34f~t1N_FUvvEePzWv0-3Wm{C)V;3LRdxgzwJ|L4Q6y|! zzEqqF7pm6!iQO)*CLp)cR3~7}Iss5%jm{{Z)18Lus}oHW?t!5H67PV)H&o$bIGs3U z_%p7s>w0(Hlv%9EVtA|zx-b?VOdgAvmdZBPLQ!>6R#8PjLL+heUD*q&i}*3UFij#r zdsS{SX_h;kkJs1gXHiu`$qsemoGC%Sng6v>@e&sV_zoi4120COi3%ykYd9dSq=f5% zQR*^8F-P~arz%W~{KlfP)0yz??0X7D%#;SxtN`kj<$QO3LNU@A5pFby{v*7~5B(Fc zyqr%PKDuV*6P=DUnNXS;rEc(o9K%(q@=R}DA-7NA9TYm06RTm{xUBwO}k1 z^ZCF`$k(hpyAucBNMo$mkUh&q;M;L&-poaR5IoTPMlaxlh;%r(8}uPxp!*s6dB*?g z^T#6$O6^mfaVK6YmodicZ4LYBV&BgL;4WO?QO&k%Y6s6}KPplbYCZzFcy#S2HJ{Fe z3nZS9w!DV*1nKB;W~X@ga9YiKzT0l*s$gk&<|Z(Bj8OSj>u}MwJakYt5GrBX>zhMl zPY1gdXLk7Ar0%*1$0w!Pr1;4ZH5bjU#*{rgtQ?;3GTMbUy`ZCaFApZIaSGk@p3G;! za)K_+-7mX=RWA~bz3}l)FqK|rY@?+G$IJ7BH`-@Eja)s=EuuTjVih~9ZBG!LG}r!& zLSO2M1X*m1m0^EIhsKDAe%Njp!j7&q8^#G1c~$~hal{|#2;Q#UA~clsp{O(Tx`p~kkdXCyHUwosHZ#S32COuBo1>{ z_E?(^@gVI;e^qF@CU{?fe_+AEk<-8*P8xne&?tVof*`oKf354AlxZmWOIQ95Fc>;b z9QK9|y3G#^<(sm514VMHZKiijBeHFCn(MR;fj?|W-w0@YH;GUeDXOfO*`)T&gEzsn zw1>|}g^63#heKWP?d_t(vORVcE{>+Fam>m79K7#-rvCZu*sVyVAV#7;q|e~nbzo}d zw5VdqnLyc1SBB{&;bl^uuITPGV2xuItw>m4C474jRCkdju3R@-C@#TR_pL>ylHU>Q z2^GA#bUJrS>h{4CDo|_omtlWwLP)v|G8=*l~#ojEmZ}U z9GU?0r<~!hFyac7H5XXcWVZ>f3V!>7-vL{^wm{~14L=UfzAye7#dPBQ3JX#jE?PPd zr5BOX-{(fg2fB^JoAlePIY$R+u}&ehqFzc&i#j*d$X7b7s;lG(@^#gJ=y3>=lFP8M zhP_4-G0?F>c~`}CIhj+aL`xfB7LNRO@WA8}gFIgGHQ!UR4;47>@B4CF@mhY*y8!hX z6MVRohgZ;MV1xMss0$R(zu+*q`z+72vh(F-&ngo1CopVUbv06nGjS^BE{_u*PLp-e+9U z2wou-NG^8r32^=-EP`flJgM2yDK`Lo=$|NdUFwvdU+IfpG`$|-DA=rb>$#FoseN)E ze?SsTGLfzn>c67SQPLZoj0;!G%NPQAkuqKU$i%B0-EF!dZPFWyG94N{DzM%ih9)mX9}$^FeNRH27=@h1b0!>gG|D4_+dP_5MhZ(1 z1C>Bj8;%!unLWGzMOm1w2D76@OdpQ}D=Ay9rUi+YU4LO#^Kvt-z*B@o&vQ0VDPfu_1r#M_~`0 z^+Ko{$rn&~la-!f;b{R@XVPbuVHT=eJ$5WMWw*LCbcXw4;Zl^})Or^{bH`I#-t3;( z!7dx!co9#f2R_S=9#uMt2}Jc#^K1SffBcAY&7zn%gtsa1-X`biZ4^zOAdV;PkH-%T zCBZk{)kMh!=7Nbp%UWM=c*iUT{)PSL%CZC%Jhc3b+R=on296F&Qs=<*9+P_`l-{5a z)KMnISn80J8ZkQ;w^9ABpw9;2B&%8)eLcyI7;2$Pt= za}A{Cj>@gdJkk+wHG8m`s@!Vs393i&5%*`(v(cAY9wTM{e9t0wTAO+l2& z8nbKe;@b*;h8dvX&k7+ofRNbyhkD>dxa|q{p%FkPOB5{3r3rtZ{$8f33`&w~npy9$ zZg6wLe8is-xv-Eg=TF}Qt5;yf?gX4{$o;p>>$jP0fr{!fL9Q*bL03UdJ}jsX6EMY* zJdoM!K}EJb({BSr9eGdgav4KMB=g9zs${Y1LmBi*T}FT2Xwn&LFg5A$HySzhyMNk5 z6{!I}o|(TJBu{0eUl_Vxu&*185g8)9C3II`Mx1R{b5Jkaxj2XZiIn8GD4xmM8Mvco z(BfVm?;U-nPNThheSyV0Ek)m+i*3fjFsr?t0qr2k*L&AvjtJ>Fb8MLzG@L$}sD_RZ zvpsAojTA;Nh*-~#4aHY7FOHkYL|ZS(cd!QH^IaImocvz+`e88@)`^k0P}FJ{BddgN z-fa(;1Iuq5`@asJgzU~sv=Zc*&1X>P#xf1znIZpxTMDofZOFFA{nqh2;&{pMiC?uq zJuo5zQ8=GHDm9`Ajh1Z342q1OBuA7hEhTTwg9t5oO(_fh5dURi>T_&tvA9=|c`62A zscaF!}eHCF!w1_Okw<2z4N_Ne$sXVuY3cgJrsWTa@m&W^$beGzH0w? z-a_!#N;xs&LF;E63P{Zov!GnS%8w-qpQ)_Aa$64U0gVdPZbJe@g0yL5D(xXG-N)4J zNb6_ux>YP?Xj2GUjqR4|)I~5m=ZOkX(Aj?~u7f&e9H30380*pJHy#4di0vX%eB3WX z))?U9ge^Z=4wM$lGsAq;$M&@=rLyQu$J&WEh;ty};fr%XKD>Y;^Bxe*wMG@{$A%+_ z!sw7>`uT}+@m4JY!V{)Sp_GwKCyp@x=bgO=o+R6Q%=l;#e(I4pfw_9uoUnI*or1!3 zbExxA&5rIa!0>x%3;(Nu57Qat2JdYrrKWQb2}w9VmLPU2RtO1}wPfEAzKOzU$R@XB zs-b2w=n?-Bc2UV?ylTLYTGz!WI@(`Bo*fxOT>e;n1eD{1y=~|_5ij}3>4?kqaUaYr#$GmpS>^7N>TR%0W7vPTpvv4fN4Y)Y2 z;znE_q|0Z<3om}OQ2S)Pad?<;#YfUx;@^sCGbPycLhP@Ii`XwY}!4VjC`F1pG z?~hySh+hcAiPzGXtrkyuu#}zmnhu`%?22--_W5m!*2YM8x0v7i-f`LwKoX?C(t{za zdx1R>E-RP7D8bz`q~;OmZnF;+q*k6#&?mjjUwwc-?m${zanPJ33y@RnKPeJI{XP!9 zsb1ho0o$f;E=CiDM&<1}S7dhhT11*`O;3#17`REM>U`N$uI@AIPv_wi6;4i_Bsn?5 zA`>6*GzSw~G&XQAoJ5Uw{mkc8%VAMBwb`+{;xPS?;Gs^liH8b=;pha37Cq}mGWUsc zSlP5zS2)cQE;Q>G*KGiYG3}{}Dm5B0er8#f4$h4_Jb94jV%S{%#jyq170cip?ir^p zx#qJUmA722!SlJ^BO^pgNNXm>v_8P%or&7KOnZT9N9Z6biP zH8_vFHOL7TnU8=oBzmav%b6ol>OcB@Q&9(wDE5$sRWVuY5u))u{|Ey2`ooQ&`Ye49 zaGyCHawg@eqaPn{n)nrv9min&g&+!Hh&z6Rf=)RuCDVXuey>>?E?rS~>tQLW*1toG z3GToO^JPry$2aLsgsQdbH%QH>P<1k-#8&W2*Wm-=Nt=Lc03lrmm6Ugz2YYCF#lw1t z=f_40&R=Say5-8>Eolev4`t~aS8~0h#sD3Pt27yrD>7E4x&i39pjrjrFd469f%%Ue zh-@~QFXOtaeMpMD##u8i&JEbN@di#27EH@2mXJ}EOX9@oF7wCKh9ZvSSB~)2Wi|Eb zyF=13MmB(6wdB+zSWgV5ZsuO+ejF%oG`H0mR!J&2stAt?6E+fl$Ei4mb@x+#9arTy z8WU8>JxP8uZcxY5HEtW z)x*~P0mYQFZwsw3QK@v7r`7)VTaw%F&fF5{XbstM7#n zHuwO_m)}n(E45(UiwaJ3G@YnCi#&}AvU!=ZRlZ^kVktkz(e2+9KE!D-K20r66#2?R zm6|cpMDX*=qIn*AboDNI2Eq>WP#cB)ccS-Ph`i5U7kk%8_hNZIuy#4iHe_>tVKLF5 zjb4i~^(t*2@<_+8v0rEJ>IbIZ&73KpWP<@cd(yhno5yab&CfnwKkdkLUqNS3l*!~Q zM^%;U){hljF=?Uo%g(ikIPyooY7WEx%S}~-!!?emU4pD1K(~U zuQo1MU;lRUi9zsozpeaJsj+wu?a5;z>tADc}kZa>`|X{OCNfd$~@0K$uS)o(n2f1&*#M<_W>ZEbZp>f zMq=)uj*VU@Prc0-_{6b7ZJ^)O7GmG<&llt+!vUJNu?}8(r>j$)B0{3-BJDRC%p-gY zS4&d%qVuK#jvXLNK1=OPb}OM)gV*_xd;tB*RwaM z7L&3?NXn@~I3B~-QhzKbjv0N`c?yKx0QVeH=o=w<%U}9KfVb`$)5NYkAk_Et+L6@Y z9uHVK;Wa+so)1whg@G7Kd9j8hyvd%?tq>|NsS@@9KLc?-j~5Yb-p^Z{4^EGNIK6*A z?m#N)4m%RpPE5%^MUCO75ZP?}NV~Z(AP^S2(#npMVDGE&sxt?g1s_-W(xhu37b?%yVg8 z%GbPSRv$ntI1||NF|KURr^b{@?k4Od>V}EP-HL&ED@fFUTn`;LmWN4R-H3u|gP#Vc zpTQ{QlGjPRwGr_J7pb6kqy3dHf_QO==-YmxBUz`Bx9a!j;{P58mxWuv{$G5NxV~4T zzJ+lo>jwY$Z4W^l8c?^Btt0F(rwg|&;;K_&y1FP&)Ze1NCtjtX1EZe*BKzFP;qV88 z!Sz`x@=wkACHC$HvzOYfUvUCWrBfx%))d@L0N(I&uw2$H!D_QSZ>)LuoEqI+77zIw z$Lhq@=B~Wr|KW^7UvAm?R^}q#ggTmx;35rzDl!tMf6}VvJs;ev;nqv&l6g- zE8u;|RV};6+o2g}7B@R6w<|%Oa~+`=+g>bh8SaVHM0|*{5jgRa*+^V=&6!*OaHs0P z!%_>XKd0VJPkkUOG&{S)IJ1|l0-S2?H1%PW#1YX zh0eam6QCfUFwnFWq)9Y=e5EgY7|Bio!aci2R7AL zSs9oNcREqd;^kUL{6PPb!!jF#r1?C&<0g~+S#)bC{o zkq=yDsGg4fKr;%yb#rBtM^DqHrhqL~1hExej4VEWIxFFwl#p!838J`?lo7L3-A;|k z!h70Wt(0>5m#ncvRD=sg?-ZdP@`VoiElCk=mZ>e&gU!RJlOiU*;m@G9_Y*@`UpsRd z*s(b$I^TLdUe`VUNxPdcjLT&a*y|%p6u&SKQ9e`$!y}ST9H;W?)r3`xxBx|GJC{;> z8$a&%9&ZlG#_mMamuoIDzDZ!Q?h*NmWH;IBb*r3W9<|U9A5EQv>9h$$EtZ$V*2n|bdHgR^cqk=%NdJnVlPXPlyu zl}J{}Mi^>G7MeRLALbU%m+g1=V8)#F5b zyba9@`Y>*3d0(kNq~_P7zkH(dT;Trc*<8a>W-cijXHLR$LjYQ51{qKj?6>L<6+QA` z_<#WW2&R9j?D26=MfL1kSUnebF|hPNv$G%G4+b1p*y+X zsP+T#G)~}2CE&0|O#vdPaZzVHFxBx9>;A559Fzrk$E1{&9(L37Bo2g!zlz7!Nb)&| znW?lIW9+$89HjM!ON7nT%1`LpO@3XpZT|W)U2nA)WMQ#xOgSEh3iR--+s}BmX-PN-3W30pZ+-8G zXau@p|Aku~bE2~CuJAT~;W zo!X0L+|sWWbvz|Gg}6-m{uk^A^5wdRKD(H)!}8poGa`hPO?{>-z_uR%KbmIB~5uQT+a z8f5z%Ms_??^f<>3MWwPt!gr;0l?h$n+q0T0mf!*9uS0J&K=0faQxlA3tpn&LEP;GU*}m5Qn5# zGyV(tUoy23xk2k}8W1S59O~8a777{%g!A9_!hKfUh{54PuVU_31L6c&YcBG?cK`Im z#rQ89+`q4L3EI}=f4Q@Qtx2G46Oyf|{to!3nJdY^0>l!2TED0GD`^r5*Bfs-5a^r} z1j79nX6{Xbi48UI7wDhD#(#rinL!{=M;9{|dvh~uV-{Ha5$tDCo-xhu1mz1=^8DPJA1h1BUN7>00O6=UwC|GP_Bn}lr} z>c2-(c{k}~_$n<^e+6Rwt3-}d0S=Dxq(2s zVvY`O<_>O#YTi!fuKNG0;lKO+=gs|hl@9m+Ol#(7>Tds<<^PAnz4iX@_w#<|t?xN?*7&h% zu9_Hi-*Z&Wx$1gqgP2?ZM^cf8gu()YgM$NO4y!{V1&3e`OJoNB1ZXQ-TO?cNx60Eu z!obKBT@^dF4d9N!xH#EJbg{u_bp*Z{|5?^8X&@uOG9F%qj$Aqm;!5gu6%_Y&&|Xr@ zQ$C9%AQu@xupIfC3rO*&z@@W9i&m_ZQA%j}PCY(|^GlS49!n)qN;*+Vtg8jC(FK=u zb{HS1?!oF1zwFQKeW5PDdghj`-&FGag!#N$7OSW{$373ikB}U>2?)QCx z?9*D?%Ttmmy}b6A-A)mY^0S~KpNYNr^NBB-;9&AN4GvmiI>1Lor;pr0m+rTd!!$$C zyEKsCRrS)7{MLGBAvfR$?j%<+b4V<6a9_GeKyjdt(E1KlPTO^&ae3l9vDt=)~1>-CQm#F)&E6g@S9rMcP z-3x&@?vYk)&;512h1lP7zZilLztaBF36;eP1EsEV*iy4%~^ABf|>n$ITBfP-YQQ%DS7oKR2CxV5EP?y;%$Q_Yue*b}jTnqK1- zPS4$f($o;nulyFJ{Ior~U=yr3ekU`FqEO$0Dqi7^1l$2On$)_;pSGvVp*X= zX;U}_pPpCbAfS~S+gSMo$sGw5E<5+bEnzYjv#AiY%n8Eq`HjPoaxE+mhNx?;j{lO@ zu`XB-YUc3@lTus7us##Rm^2#xs<@jI3@3>!1!G5Sq*q?WlVE(9>`e%|4H+sF? zKaVaZ6U?=q%og?+`7 zGr*VBQ2WKBQ%Kd<(3P34Nj)`>L$!eIrFr7-efsQc0=H7zhePb|5Pt=qlQc>C{73Lu zsDC#5e+5q*ais$28#=G>CIGsW1DwoDigdISumuu5G?hG0Un}L=J9&XXUFUKv2{8;B zsvdArP)2n}#*lYsQ}yzj(~WEmAsM|~is@2;>5^21u`CrFx1S|vYk<^4ZmwZ4I6iB6 zbNI4=_uW3qw|H6~lDlOA5kH@n7nd+PMaodd?ibCvq{Cx<<>0V_l|2%95w~etQ}5DRSD!`8BGfs_F^c8|Z1y4o zzY&8zVW*7F$|-B^Tuv&)G8iw7+#!Za+njeL@g4j=ekNKE1&#Ewjg`KzgmKw~Ctt*R5^x^C5BP zk)+tP5*Q^s9r>al1lKuHkN+_oAUL<~$(d$+sD%Uo9gTzv%V~MhE`1o74YU3HP-wep00X6||6|;!Dcn<_V@{*xN-D{a97tulhxqvbrpOh8!{3y+(iX zQiZeww?TP*9@Mgsy80sK23mgy&L;-h*CjQYoe`!Ekx^bF7)P157ZcZuT#Xn%8H4ye zF|7dful?fM;tga%IBdNv#t~W2sXi6-e1&R>^W=PcYb=Ssv9n-yB@N1y{2->msT1g= zTyo>%-YKZA(wW(e`VSEBf@Re_-Mw_EN!Lj)a)5@mXJj4`803`170qD@U!99vO5;e* zcTQR<#qsoxz9^RUBMcQm6-FEHC!)-V4^jXy9<-@$>CB75Of=(1^dv&Vd=NfsUd*L} zqYaJbN0!63|A1oUGZK*gl?7+XUFZ+bT{mf&_d>Z1A1qil58`hraInW@a`$m=C~-*h zQBIf&GiEWGfjn2J$7|WbdbGqwSVuS}9302P6c|JF=oi{_=_}T!7sKct_R^yz{W=El zM}@2qePq4gi{1DrhnwX|8J4&kV1*Ed>%WY9%}_B?2NPhYyV)3f}wXODSZM~Id)h$4@kqa)OZqmD&FTiicb8PC7G!Wc{lgi!a)YI0YSpx}c`fA}6A z*a{#jAJr1#@&w5}Z`3a|WGp3Pa(M!5ufm6|(rZ5Wcl9%IPIIKYCAvF-R1Xzz2V{+C zl--oVCdrqfo7c$9|24G5vr`@#BY3{C~tDg!lT#LaSUPtVu;SIN06Sy`aag$#o< zwg`&UbiU3q78FGdt$I{iZC;}(doJ_RC%v<*M`J5rjWgW{7dqvrs}T@ zr)ZP3urXF3h@0bNE8iG=ICaF(M(R*20+ZzgikPXBBk8i-rL??GI!?XU~ev?GcOcQ?r{A5;|1f7+seBo1-9HB`soHsDj3}F4u-ALb}s*b}+ z-fsPCc}WoU0);*xi7LAYj|JA(z#;L=z+Otqk4GH?4^oeM4iJ6A3?5*j{`&0g#H;#j zPPV%s=UfYiZBZx_5f|`T-~uQ(oGWZ2)T!*Te&OhT`?i|=i}7GiqkY97;|g{VLTk*+ zvRZ6tNI2STXIta)k-@_|#+UI7F@#D_4(BKMm_%7{dSJPD;dF!u#9PemXYI*?3%knOP%2aU@@`gl_S^NeGFU#b9Kft##M+nkz?ZuPqWgl~#phg(w{2kYX@86sHNqrukdKy7ebbZ{>zpb?2G&wk4(g0&d_fBt+iB+2`55O?XY$(1@| z&p_2Xjq2tBL~`_eu`L(ol!pWkiB+crnd#9cbA@65(mJG#Ei^yeWCodW!87iZBUn^; zIe!>cN+^^*K1kRS@A8z6T{mfb^thm@TEh}}cPV3A;czzUyC|U!_}hNjSk3iQdFPtr za`Mb^`*NdGx?kC3;y2-4w`S+(MEK2SJ_$iy!aSQDP^5E>81~Gs3UGBac<3h!5#D&p z<|Faz6a?B2?rJMy0mS1gRwN9F^k6PGrW>g(S+0@ggIE|L4k~X|UG4K4+77>;E7e{lh3PgFN`OHVCpqVHMjAG9{XpE{ z+P(pQ_~Ou}$Dqhhc|(~;a2t{#j+q%NDCN5WI90XzEE-R4LTKcq^=fWC%ruhveO&P_ zhWa+Otx#hPg4V9bLke1kJ=!P&1)VK9FLg+vpbd^e>Z__m0-kGOew-ou5I)ZuQ}-+j z3}Zr#4ob_pauu$Ys4&f7&){J;3aCXUFi+v&Tzkf?26{<0LK+41By&F{5dYwj1anXc z&{4!NfJ{_{(z8UxVXSN-%yG5Uv5)^w25YpG$A%eC!?HLhLpxPOi~mc$dcQSucuk_y zE``XOTH06!1trHX^En`;+sD(s+I8;KS~wdo&>4Cv+1MUi>zdS=^lJy=V4>=Y5P=8$ zBA4@EY-f@0an6cIrfEC5bD>_bzf(h^I9mTb%!d8PNcO}nHcEkRY z?S0!O#eCCl)@F!|)R2%Of8?)m+@qn?>l|OSMxn&ia(3%9H4qJxf)&)NC?%OsXzY6} znKe;UF9>wB%XVdXk9T?Xhc30L#L|wCQKm$pf@TF{L*N#J5kedb*wvH3Nw68d0365f z5%_;H8sp0LTG?MD481c1ntk59*)zl^iF8+gm@F&O?tacRr}WQ!%Z?zJ97FVQC& zrA$HOyKkKezQ0aCz1Nd|BmMIy8o;4_8V`?}H6cCm2)7C5?x|v;$@&_eNeev1`{Tw( z=CQ^?;sQg9=WmsF+8108+M%Zgkv3O;ltb=bv)Qy)QgvD!KQ{s$ zgGjHK-8{0~A3dLaD<*kFH8~YiJl!ex(iK%MOBqYTVTY=Kl0hc(*@X_^K)^o6n#_p& zgPLMxS)gsQL7U5-^I6CR{$T|Ys;_xexk+JTXey<-a*>xgGvX;L;qJt8Qno+-3Zi>e zH%vzo70lR=D#1m)BsxJKc$lNG+hEQl!TvW0=b4~qZ?Oz4a?i8Du}i+MPGc;$+nTlm z8NM1PLYs)r6fu{}Esv_xi$I(JYK`n`wU67|?!$dKyHOkw(! zb;6SAf&5KkaV0HcI$ObN-E`PVRoIYRA$X#fpxnaXH&i4YNK1mZ8#T_~f*7+!rMn$Y z$*VY*#3jYW5@lG9u46iy;VK+^8;0K(>gk(bC-yYFoo;SB6V}p5Hvm5foxk1GCNrdl z?cPPKjRF1~C3{KoX6UrIng@NMO0Wo*Ux~@!r^j2qV>()U$4l=*K5IGF+DEgq;x-VW ze=|X(hpEsf*a6}NRABrl!fVuBx=*<}<8=uQ;b0ocnn6k|wzfySxzm>3z~i_t2^jbZ z7VG2oZJP;9e4~`_!VOR-?Sp*(#J3_VuwyYJbRuxw7{CUQFIzVoZj@0lsP4ed-!P4#557>%##lqq3Lze9Dvt)XkiztP?~VWp5PUy;a6%<&2%_oQypRYt*B-G z#nxR-1cQ*xE?~V0$Iov3uxXxCL}cMkUhF|pNK=;Zg!NBBPBJySYX0`Wur(8Dkbnbyv4UV_TJgHEiX3B{x(ePSFv+ zQ89{#oBxzF8!JkV>1hdzrTGjJiLlFIP4KxLaDCN=oAjcolNzrNa0%K83bM~_AgW-@ z&ww#Ak-I^Qtr1|+ho9sFCtMDExsLSd6CvB+?qDRrslL9{&g%!$qLgHqqKjOb@V4w4r#U_w+T! z!JE!Wccgh8PVsgJcpXP`SyZhNyBsJix0wLP(;UkA4t~1Ptln?7bO-N1_Xv4{Q#S;4 z=)Ju&5P!@g3D`{ReK>uwCl)YhnqquN3eu{FS9Ug=04-o2Fe+$A(ziTgO~&?%ITni8 zWR?#?#e-k+OYN8m73wX3fI+n&5Fut5X5zaGB+1a&8Fk(&Z;#5u19a@NdLfo`)gi8)xYEHPEgT*NvAS)9Eek6#$k z1D1oIH8Edig;)BL*PMZ{tIS~?0Jl3zC5G8Dnr06HwHOXA8MR~{HGC>CHI*do4=o&B zV!*5npFR&@NLQjt$+N1b_wUn*9)~3Cn%#}fvKpNQ5vs&}t2i+Nn;%BRAnmI*PEzq1 zF0x|_pA=HLDB2DsR-I=en@y7xT)^R*@$5^sHBXe#?BKxL?q{{8+ zXyzv7N?<&iQehMvTDIV+U#b9R@$RLYw?BMb-#q0@$1^WMwn`c9j(WDA47i&m^oWcr zou7jw**A%0_sXVE$P215dOr((r6V;TMNf(e9^eIX=3ydzFj8%9aK$Hf_F^PD;(ff#^Zu?ZA7gLOt! z@`t5n(yK*Zotx=0G~;51oaPaFBTmO%Xwv|QaSm2XIpD+SPfO^-oaR0VqZ965KR${) zoHls}{rp5f*cV!RMIO^GOs7XYe-IB~Xun4x+i@-y+281#?CKcr@@3sff)DhiNqWMD zyFk1-(?XDH_QT9(ErMlHiSBvp^9T74rdg-R;G~_^E+Ac5UI?zU?e20VAl|IPeQ*QD znA<$*3l^S$pNiq~1OD6+Eom8V#G5M9guUuN9*;O~R~QY({Fc3Ftlvh}tJ#tpXsuRr zYdmW2PK^NH1^{fzKM`$T*Yb#8sx^VXV5_h4_+BP6lu=RYrRwZbJ)c9EZStr0YdHY| zJNu{sX*Y(LtM7far1|gnN}^`5^w|KMs}Oe@_y^=}I+|qxd$DDc%TbGwAO?4);*vbi zes2+vkm}ua`B|hZ2BA~zn)VPS4XaS+R)x@;q?TSXdAkWpGO4Uvt67twv>R26me5*{ zj|zKcb=Sy!=tFv*GKT#eB`UZyk-nxrURuP?YCT;}-klQKOIXo5iHkkzAymK;1D%?w z_X+#*j(pxO(o6I+!$~kmhah$M7)sC&-WYnre z^%Sg!cOTZUJts|Nzac^q{t?m!d~R_8Vs6*r2%3mwlK?dW9ZGe`P6eBvV2WIq7&U4a z7)Ug}JO^dHBB`-TMlLJTnFz4iMXzd9ux&ruqv=m@l`%TkCM>y{ey6;!86$Dkyc$>E zYDwsv#MFklYw+)AD91H>Ac+%snc7IJV9r!Ol?Y}_)}Bxobfvq(qm8;;+=3&npBf-~ zYXrJ9p9Y2sj^sSn4w<83Ak4unJhEhND;(4;P6Jd#Bepxj-CLv!XY)HfP`I2^c_y~R z<2@sBZRi@wZNPC=%k5SlLeF?^3K;QxBAJzGy{vJ_ofj%`%|n#fpRriYx=W<{WHWHn z6c0WH=Zl`EfC2IS=NO>g`(^qsp#sSg2AG6GS?~D08S>*So28Pv8~VkGurJFNaK*Ic z+5qb57LHV#tg)wj(W)sC@oxS8iWZtxnAf2yOY+nFEPK*h+^G5b{F<#w>k`4w+>R?y z7=+`7@pc%VC2>}WM_1rEIMPqXq4a{0VQj4Nl!7scaNl5B93{W37R6DesJK9dGOPYMJm*h@X#o_BZ0Omfr3Uexp9q|3ggF@j-kNrygzJ7? zPNyOi77>FW)O6jb9K~aJag=mU3*#ijh%K|xOPigsusLe~cX{VWTe_BbEt76lCeOTe zDl(eLUI-W7%*J`w&)WLkK+l2Zl6O+RNh&VnDwQ z)}>a?W7nC_3BCB(jr$hM%ORv$tl9OwV=8d6PJV>t@FM+Y#CKw6zuDNcCVcCe^GkD4 z1Zq0<{!$Nj4lF)vmgFgavT1WwUttbA4S3;_cK=Zc-?``F>P|5%wsrrQ7+0Wi>9L1t zcG^c@k9TJIYU)yE`ieWLBi|@f8o&b4px2I(Za-q!&N5EI?8Rlr%vCv2{eaHVg90};JD3z>h|J<+V50qPzOk`bJ$aGAp}s z;<}g}@_ww!!G+~RM2Sbov=Pa^qUA^wjilA9*m44_i7;sl0J&%y*8Gld0et_N8M4^$ z9de@InnTx5PDTI9u#O0(~gXh~)rru4BIYnev?m z<4`KvG(b0&@SNNX&n(3JM#f6Y7s#Qrh~x^qW}X;95gc8O^%vmdU%A|HK+XNhI!}qA z#{2{=hXKt^!Z?CF^vi5w3qW=>Mq`#%JFI=kIIa^+?sGXgD+Y*>$bnmqU3<2*5(jHK z9eQU}753dWh%$LhFdB@g?VaK>zHoriZCe}X`h>NXEmKm4M=62$yQPH?x|d4 zC7<^m@To4&ZT)fbk+?ACR`_=OeAtOMl)b&dgg1=??&UlZ8G!o6iw8R9r5~kiTTiLA zDwfy&QOQ%nP4VvUt7^kEhB5QOT|iVjXl-cHP+5EN1E)*l_^DAoA0r!_<)ckN5R)bk z>*u7p7K;^Bnr`X&9iY63gavnX3jM^i)8y%-&gai<(y7+Vf|h>h`o~`BG>&{kzyd?) zTxELAX_`4PBZVgNyXS>Nt@aC)3apq9`C+t)+a~dTleO8WWJSByv2g>p0r?dCMrAF& zTtGAk-XipkMbr7$<){)7#fK7K7rsbt4Hh{cUs>saWTO1RI6#qpb@4~~P+K8rcP&snuDfv48Aaq4QD;Pnu#nymtui`?WhRw`YS`otPvgQ8rp zi{gMOQES8POR{Xy?;LazFAjuCS_r&8D@fLGC7>s;q|@W$Y5?~7Dvaa$c;E$h_wA|; zclRy|^%HlsIe=O6)a_Ho1{Zg(Nm+UjV(7yGaSvDD&KG=X1MV_>rltI12p+s+n2H8)$H;3xKabp_P!_UFW^0hFU^&rcc)Q5cJ2LpwwHzFowWD|monIrJ z{%S(-5|%@|yd~#|F`EmoOgqT~t6oY3Xs=FBHi5|q5)i31Z7JrdD+-Qq^hssRg1dhn z4AM2AZj(b>Y^15$XycvEoPBc&mjbl~dBRrgSQ#NbVi{RiQlTrP)Qy0{dmj!mtFlkZ zgqLWZvbOI;6>Xk!z4oC$FJP0q-&Z!!?*|S0# zn@hi>1D&|eW0%)UZRZ>>^pw+NbX*k$z$kd%CGJ>THH~}nwUD$^m|d~{j|k1 z+w5Gvy@C&F-S^$iz;t*nH1Kjg+*#l?Oe+lX0zQ_>JQ8Q_+nm=0-g4B%&{?BL=@pH> zMn)JUX;*!oOsmN5!=4G2Oz`%NcHndCuWjcCIO>Hwku;BFfe~Jl4foj(s^zuXOjfCo z82QN0n%Y`Uko+Dn&Gs+(fx9jNiN@R>*27@9B`}!|V;-w0pfD>tgZ4Y&hMJPgZax!k z5&-o3#NyxgGWO#1>rIWwk`ym=1zy`WL_?X7IH|8v`Foi_0?iM#h5FC;8I@VeZFP+jLbf?yB>l))6rkzc2n)y@6+es31oP_`y+k>3m&-ptFHf>EZlc1VocZ3a zUb^6|$LErf6S!*om8$r*B4A>@>$N>U9-e=jXP>U@it@@T-r-n0r>5aN-z+y`1 zyDN2yg>HEo`Y+0C%W&V~p0iMHr)kL^;O<*tWcVKuQvxIlAoI@o`1VHRMWIk%xcM2P5dfQ1 z&d<&EN*8iC8JHt@zXYp6St^YbWDp?Rih(IAXa!16HF2yRT^6vc5@T*ypAQKcD=cs2 z;qn31+lT?q2y;G$ho}%IZctvh=9$R?aw8#ouGU@UY;wn1gI1QX0ZzMM)g#Jy!tXJ6 zHKHQw$gNns*lNXRgE%_<{E8oexBxOYyws>GEpX^PI^9wB9u4gCfRDK?Jv9}%#O3wC zv$O?b_%^eh<9$T87RR}f1U9zo^V?OXPcav3LujycHEIDh6=wp`yF*GQtLeWfJc1(5!YadVA^P99&D6~HMR)(LTE z-qU-)?b5|Zj-2GhX+%hV*V?b~O7Jvo*W&rB>JW97=ZZL3UAwV;@B}l=&Z-9nZIOEr zMKM;+72?=e)}2aGBL5~z3(vRP(WN}?qlzu6;Ck~45ieXMJ{;4?1HLx;{^0 zUg=zIrgq`waoCBdw{<#uckK3%%2xA)Hq6mhnIlaCzoU8u<~i(T%--&DGrRl;l)Sbv zF?%5pt8<&d81%P55A7b~C zP!q&yboklqP_09Y7eH@}v~+LSLsi?Iafj7Rxm=@+g5pO={Se~J41Glg4BpGp9Oz|d zRlm|MyKyX_|7#@585HGjEk8wGno$DwXXsq~+Szf(7PZUJ2INBCi`#JhLyl#W#y!Ui zqA9#a_j{e{w-RvlPEx@&dvoY@RY6~zpEaR`FOAWBG92$q@&L}?Q|on;9LEPDM#iML zCkSjT^XLT8VYlYb6%VqHQwht+`Chr6f)*+&A$=SX2fp70S~*p$C^O}!8M1n^_JBndnN z{$c^yMDmo)n52KW008)32?Udy?(QECo<4=+lA5T{+TbfwncsqCn*o2i`eK zwIN&EY{VC$55|_`r_oncERrA`fy0ifZ!>VzdqM%V;|EWOuwMBmqWi2KKPO-3} z^5g-pLR_%o2+?S>$K=P{69V|>!v$07xoZHaAR7)gS9T3KS^WbYqUwyoL%`{&`n4hX z^pk~!YPQFB)5~U_rDluv==(+sM;%Abt56gJ8&T<^kquYF6Y?#D&(Brfexz*j8#0>f z#R91udzqH53m4S$9PjFv?HRs$7M7&7GQqtdphuyd=X10~`Oa}??_w`pA8(8OJrz`n zyqijgl=-&itxKz~=oFGXbI_ zOn~C#5j$q5t-UF|O;Ml?y7DpUpRl*@(6p#}#^fwmfz~R(H*ug@4JeIS@|1hQg zeX?~K>Qo4M7#Dg7_)N0{OSf2$0NNoC;n7P<-0eJyD$8QEVFYP+v&U#GmpiLGLG~u& z2+S0pnc_2?kaWe(8=ldA)jtO-iUe&x*bRh**vaaPjw7-(B|;?NJ8z}af6udLgCA$f zDvhE-%H=Z=vWp5jz3e63lytS8g42t0AVRNJhRrC5Bl>LuC|61NtY-rc>zue(9IossD#Si}1s%pj=fjj{ zHcDnZtxN{^7|IPcpbk|F-@&;e_@-P}6^2 z^NFIFq<>55TP6{>U7}wm#Xkf7H=#rHuK`ktU73tze<%H?T>qz9(f`$F3MTP)Ha60~ z5ft