From 60500b5428a379bd8a750e344b62676bccc6d504 Mon Sep 17 00:00:00 2001 From: noah Date: Thu, 29 Sep 2022 14:43:04 -0500 Subject: [PATCH] Part A --- .~lock.Schrick-Noah_CS-6643_Lab-4.docx# | 1 + Schrick-Noah_CS-6643_Lab-4.R | 32 ++++++++++++++++++++++++ Schrick-Noah_CS-6643_Lab-4.docx | Bin 0 -> 15320 bytes 3 files changed, 33 insertions(+) create mode 100644 .~lock.Schrick-Noah_CS-6643_Lab-4.docx# create mode 100644 Schrick-Noah_CS-6643_Lab-4.R create mode 100644 Schrick-Noah_CS-6643_Lab-4.docx diff --git a/.~lock.Schrick-Noah_CS-6643_Lab-4.docx# b/.~lock.Schrick-Noah_CS-6643_Lab-4.docx# new file mode 100644 index 0000000..111e6d0 --- /dev/null +++ b/.~lock.Schrick-Noah_CS-6643_Lab-4.docx# @@ -0,0 +1 @@ +,noah,NovaArchSys,29.09.2022 14:42,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 new file mode 100644 index 0000000..bda527f --- /dev/null +++ b/Schrick-Noah_CS-6643_Lab-4.R @@ -0,0 +1,32 @@ +# Lab 4 for the University of Tulsa's CS-6643 Bioinformatics Course +# Differential Expression +# Professor: Dr. McKinney, Fall 2022 +# Noah L. Schrick - 1492657 + +## Set Working Directory to file directory - RStudio approach +setwd(dirname(rstudioapi::getActiveDocumentContext()$path)) + +#### Part A: Preparing Data +load("sense.filtered.cpm.Rdata") +# load phenotype (mdd/hc) data +subject.attrs <- read.csv("Demographic_symptom.csv", + stringsAsFactors = FALSE) + +if (!require("dplyr")) install.packages("dplyr") +library(dplyr) +# grab intersecting X (subject ids) and Diag (Diagnosis) from columns +phenos.df <- subject.attrs %>% + filter(X %in% colnames(sense.filtered.cpm)) %>% + dplyr::select(X, Diag) +mddPheno <- as.factor(phenos.df$Diag) + +# Normalized and transform +library(preprocessCore) +mddExprData_quantile <- normalize.quantiles(sense.filtered.cpm) +mddExprData_quantileLog2 <- log2(mddExprData_quantile) +# attach phenotype names and gene names to data +colnames(mddExprData_quantileLog2) <- mddPheno +rownames(mddExprData_quantileLog2) <- rownames(sense.filtered.cpm) + +length(rownames(sense.filtered.cpm)) +#### Part B: Filter noise genes \ No newline at end of file diff --git a/Schrick-Noah_CS-6643_Lab-4.docx b/Schrick-Noah_CS-6643_Lab-4.docx new file mode 100644 index 0000000000000000000000000000000000000000..3d88b9c5cb9037d32ed3ab56aceb33f300cc9c60 GIT binary patch literal 15320 zcmaib19+v$wsy>pI!?#7la7;)ZQJVX*fu-1ZQHgxR>!ujKRx%HJ9FmVJOA2Geb26{ z{jS2QTJ_bgdSxX*!BBvpprC-%)47#^{vxm+-*s$_EbQoLe?QBoyTw5n5Fs|)GWC<) zoEnq_4Qtzjf2N9hct?$$L8MBG7O1bU$-=bADtF^DJ?`%l<*j1BcDNSC;Xs0yw)lkl z>K*rFZ8+T5wFOBvn(CNV<_lgO=vRzL55WjWKX@mNi)!zu5Cz!I*j=fQFu-BdyUc)Y zwwRsKzn*woPOAaRee3F}iX1be-xk6%1T|=|;Kv2wKtk=>-U%LjH2)_E5RmNu*+EDjzc92ikhQh4wxct!vNfW023UB= zOIh{Me`-HhQEX1E^T|f~7QNVKB#;7TI%BD3j<1m_9P;>dF9CD980c90s-5wQjW8Hg z6YQv)UUCzhm+mi^WfV`l>M`roDz=67EEMdX=C3Db*6F&tSo6dEoyVaDVeVN1K6Ekg% z@Ob#r7`~cK7erbbD1WXikeAttZ`v*Q2w@b%sy#7mn0;ne(Rb_yTw-g1Hw(}B&Y@JP zQq=guOpir%fd(55vRS(uc}ftkk;W_!xf|E~2eR4HtJW*RuE0N@TduG}L;IbY3^@(2fj zr;SKG>~q!9{kGYn^(}@2Fc8quKiUlOuQt=Ow*GCh^q5i09(qKf``1u%eTg(MC!r&f(6Xr(!K@@zoB&}e43FKZI``?G)~W-yZOa_OsZt{BfZ3-y+Cn4^A!iV zexcCys$sBADh+=x$4Ahx-i+b?=~N1*O4ki8Docw%bw2KkNJeSe6eg=&gk1+{8vKbH z!h@aXZ9*drM=mqXHUXu{RrGlzyBV}MqKWPKonL~c;<_0yg|)+=&SzkL=F_gtd0aHx zk-d|VA8^-3YY0l6E-Bj-bTgHt4s&+tx~T4{=G^)iZN?lCLyl9H&!KgfS!XaavgRxM)K_{~Qh=kp1)&jR`rooYrm&;DFno&V% z*?h6sU7!{U(JmcjvVBf%Rul z7Ta6|(V0qL!ZT3~-fycMq}HaZepp5QAFYD?SF0E}*x6eF{xHl$HB1jZO2Fo7A=J`} zB4qS*T^Jju;(gRw1F)UdkW`q58L=Ya2^bGf_KpEl^l6f;76$KCt96T9~)oVVR zSoiRTF0OypE(CR}i`Gcz}L5-K83)PG=^BU=n^(6 z2+CwDjR!c+xctdYf}K%_D7rfOz7;IdU{;O|#Zv$cwh?V3p5s!-XZbR2+6Wk)D&k+a zuH_7qZ38V+2UMp$BW!&1naz{gh?tRI0)y`~+85rl?(x5db&n_`x@j|>L`8I^1h!Ez z8X$pcK)1;~H)2}0q-B2^vP+$_PTL&4of%wrY;o$MzPD>}ft-~2EUYE2Tg>g^m)`6I z50B+)VCOg+>-4+V)lqd4sXs;&EZ9Hw`ftDdN3Z|2zfZMv>{i%fygFrkY-1>jl?yD$ zlG88CVy+f#NR>hcQB7#=SibK1ii!&If5rxiD_cz|N7dgm4KCcFp7 z?1B}C8pR?mzBuGgt(r2w-MtV-40((73r#OA-#RuxNw_0jUG!EjY22^t-b%vVU4xm${>MIpTfFsl&C^vuYXYR zlt(S1#B9uRK}X2bF>|fVuf4`^aouZ9#y@WrO8B!Z|HUoTP1Efs2t^ezq8G{WNCh*V zi5qf#!zSJ=Ig?|s?68Jq=Y2mN&09TPdb3LL=lZDoUZ}U@94QdLg$7wYAapk_-_pSg zD@bAYDY?F)!@_i@S=z&>_=(o9XC_URndwkiF$=it*8;<@k%p-!Jdd8|D}`SD%F7M) zGHpB%;L$3G8YJ5+r!i#6k!1~)4L3ujf)YrEJV#hL;D@#TNRH9sanwyHvQono@hK;# zu|#Xf4A5m8%!G6!!niI7d7$FN3O^2_qLVxkW5C85Fz7#9gM;)Co3*iXn8!4dL9Upb zAgpkbYv%}#i6&B><&w^8nw?SQ1wlSEtqmbUsPLCS$1iRs#hEI~j)6=YxelrFPe+mP zN6HI7u^FS9X@(A2XUY%^wRJe*0NJ#il*T>bABzITf;QZiDSZa694FWHt;Ar&Il6SHUr56!V9h{?zmP$KP+6-A<*TV)+PEl0g=XUD*iE{FcTHaHm zoD9^XxDoh*`P}!-&Ws-S@f#3{y27pAhWy{YRSPm8Sb$WCsBS!(owN%TovxA+AOssS zgtcC9fzR-5I9G8YU_0L)1Xs(=2K!~pb()(^AMYfdfqojDOB%sURblB;fCeBqptBmd)3elT3uCDR`(4khM#4)tp`dz zirV&DFriD4mNrj0by1OLHwQCgj)9+BF2I+y(Ym1PvX;+6oRWBoXSThCis&%MiaZFO z;W68o6eNd758vMnBM(E~O{(TsGA#q$1#fpxNJ(qP^_W@U<+6>PS|-g6PjfHZo_UDp z+8QRpj)x9UHXm~p8aJk8kk;)m_DPm$NLIVQM-AiaWcCDL=C2b=H)l^Y zr>r#jG$dSnkH^8<6xC1aZF1Tl%cVX6Vt|n#XG@lG`^3H-A0|MMZ0Bt5i+)2>5xuoM z>(UHj=wMY1-`BP^pou(UqwkRB91d%h@TCc2A!%WNK{x(x2oxlO*?h_7v1_=r`h~+g$QB7jl-fUv98XdW)sx%D9^FqOY8YPu9ujz=hF7sQe>j=d_|f{ z;;o!6fctz6qPxZhECxRbFv0F;R4<&|f;dWG~paI01;%L>tV}hps$$>8XxzlR& zAkX3cx@9gqY_ZlV0~~I zENS~I3TO<(a|0_KBoYt-Dpp*I+`pf}eQE@ZZgT}0;Ou-b_G zS})hEsef7EN|(>aSVdgBAeV)9ta&{%mt#oCDv%A33^$M9cG#g9k3N$cXYjvSnwtpQc28shPY(IE;0iV~iivN|M3baWG3I-hJ0 zM&s9;w>KM?l0CR-Cv8g}zF-ydM%j`ZR9sC6yPB7;rxz{w8Fp;ZJ z3R(woSaRN?NQ8PfX(nnIyK!;BG`D7ll_{!&Jn81e&^LV*5s-x7Br>!hL;M`zTGKp! z!LpVvBi0bqB63i!f`sGs~YTE$*C zdUXsrzbWDzovzof+|%s@-iF|v<}`FxXMCZ4_Xyqpj!%_LI?C0WK03xh{xFKF(5FIH zcCF-xJYQ~d38F{C(lQmtbZ;wBYVHe&({ekVpu0u)G}J9+9DAr~_AlF-ZMd-QvEWFd z!)$4YTg0iDFFX+*d(hq($wJ^zAo%<%B8U4yGmd<$4J7CWkbw(LbI(9E`s*h9y(6lC zM|Y{!h6w~d1AB%P4grb??)Z-7)KvmsQe30Qx)sNr4r{55#T*;mI?+u{P_lR|0QcSA zaPu-BP`Oz_HzqWyFbK0hr&vQBGdfaNvr7f8jg@1|Zlq(72X8kL103thC9Z{7hieL# z^*H~YEDyBHS~W~xtYM=G`~W|*e3Z>n(!`-xUdrR0e6`?zy1+6?X9N-Hbmi?c8rYDB zbwQ@XNRi-Jfac-MzqSPdauzh82tg!0iiG6Y8l+FZxwleLn(B*nDH^J0QmJa+6*lf+ za>i~~DJ^B-PyJQ+uFBswQUY-a1wn4;-c=7Y1MH+Hufvy*lST*qF1WJ0L((Iq?#+|w zL`@5=2W#L3C)+iXUgz-q#PlC;NNxgT03NE4G} zH=5~3_y!>Qka>nw(-$~fk^-4E`D`r1b9!9m{*<%>TrV?z=z0yi<5X(&+U|6meDw># z%FU1*y#LZK`mR_l_b<_30wq_p5ZI$I3|de86L%-2Y*N;J8cPh7j`Y%7y!Tswtw)60 zku}^ZU5`lhLo@ZOQ1^-DXVE?u?OsLSgjJ-K{hci#266KFred zB^VP`(M~;dpfG%mj+O;?>9sIwm?2R^-LF3?&_n#_0&?%@vGv4Pj!Ht=4;?qO)-j|@ z{kn$@+sTA(o}tb&YK@L1?CO$79+9?-D#ElX;kP+d^tqY7CWx#gcYa!fYg$<`+QRQcT(9}8lZ=%(c6M-#?YUL9hu!q;c%J<}@a9b1HDM`L(7l0HqB=jhBmPIL- zcS@+)8&v3MkflzPO05Eucp3i!q9}m=aZ!>&?gkR_qPqAPZMzKzg^nKo)B?-T)xK)o z6JGqnkfc2Yp*qyGGTOF;4m5ENV@X!<=S=YlV?`&PkfzW#=Pz$GzBsRkWYIe86272&xVh zzzSJitP@l2<3k4>V$_=?Ol`Je!e9p=B-Cug;H4=g>poS;zy^$xKcxt5G4iT1rVr=I z6QW0Ubyva5JUK$P(yWG(j^32AA3nC^4^E=9tM`rZ)A*Vf9vT|PC%qqTw2#hK2+`S* zoUbL@o@3C{foR(uc^DdbJeY$wv9XFgujRq4HjWxHOssH>(@qCp_2IgO-b377U=NR> zBe&bzpK%;inC>2*2?&O_lT6^OEj}Hx{76aG5561Qx|EMo>*%RX{h`<+)r?bxE2CO< zcJ&NSm#Ax4{N+TSphLYH$ON8x_w;A*-lo^QnY3PJ6(@Rv>b`-w&?WSk<)2I(P#MF;>h~ zYacO~2{`W6nk}k`ek=s7IV;&W04=Fu8rpH*18>3EDW%WdJ4HyDBAu`lX6S6Vk~!6T@w!S;vmCn@WW)6!n6UYHl|d7Et-$)ZqM*F0`Xz5llvV8U z+$X#-o}rv8Za2>rU`(_%=?h(2%_t%Ht`TO_mR7{wgIs~Pibv&{J9aI#z|TSoI`-De z%ZO~^yf@E27mgX06;Sx=fgpZ%(7Kx&r&@ZSGrBW*DAlj8dA+SMHV)Z^>?+oZFF#K# zljTRzS6GcD=y6Iy!ud~c04Fb}!dp&09O%nvUB(DCAE%`C;zU6G$yQ>DnO$p;)@9W3 zr^=n+6?JHbe^$)|^+FtB| zKOXk4pdYi=RtTLKLyn)PCh2SyY2ByQb2PM-O))m@^%evTp6*}tL+0*pKfu^+p(eEW z&;q`tm2`BrwN2iH(c)jQs1K`gm|O=E&DHxbhU|G zT->#x2_*Ur9H-iykWK8=P`Dv|k!_R5c2#bJ#o%Lj1sVs54XJ$a-efn~+jzpF1h(>j z2+WNe4kp73`L+MeM4|)cuE98NOz>Q|>~n!mZ)3x<9s4!)EQiwb>76y(&F9Y1cH+=> z{~&FhPs>o2Zdjc$Q7av~Zt7VTxVg8Jz_l<9xSi#prxq?nJ(@BkGWG<$D{h=FSI&KN zS8k&7ZzH9i8^Ysn_an()fl;C}o|Dnq36sI{00K0nZXxC6pG>mLD-M}yLo~x5r^nVEf!V%w zKeJxQj%P3`?O@15c%lg%de=nN>(a8O`_-WJ}Tw zlU#YK9q%8E)x>^bu} z`ljUFvq^2xi6CFjB zu|v{>+4cGkA23*e+SJG3FSSid?SUIAt}M}?=pwBSIwhZJnpdXdeIc5tqQ6Ix$P9V5 zR$>4K=|%l4Q%OU*gv`$UQi`+N(&%dSh0ADk4kn1y8w*dSFJJ9aGF{q^;|q(voR^OB z6}!>DTy7A=JMhx$#3tPfz%|Oghf0d@tzKVwUu6h*Id58w-KsdrY-fOh$)2 zB@UBBoW%TbTa`WJqU0sIy0|*gymfNfvaxGJWNkX{Z`>TSgaF z4=L)ahzY6U)o7}pg1#b4ECwesH(|>|`&LzAW^7wa0eg*5AE{&@Spj(rx*`@5O<1p~ zeD8em%V$u>ZbXk(R_{5z%?J@dpipN~r|{zWtN$yd%iANfpS3!1#O zxF+GBgm|21^oz?Syg$2Uoe-y$zj@l(&%Pc0%+}DT1_o>Js~5N~mtPHj!n;NlQ;<72 z*&4F4s+GecmNQl@@PK{kwo`%Y%qmp!>Fg6LK{WI$qMU#J#p zz52&@KwS2He|JZ9>sdyB9GOcGM2ffqyijznPvNw(Iyl!Fjb*s8-c~ z^tQFcL9HgKK-&1bJeA`_rh?x0mbhHdwgy*Qg0nbM@&Q$iios{%d?=2M6Z4@tqdjBs zu*K2KQ;fQH@7tnMr#8YUJnn8d7*89cHHkhirGDTNXpbCDLQybQ4s zMzI*5kUE(nX0>8$^EaR$Ie$#Y974DLEUcS6NiBS35;Ry`T3%got1;m7U(fWsH!L*0 z2;4fW4Zdam<*dJw6wbWO`Y`~zYakG?uG_0}Kb>oE-6gn%cL*Q^*?XKmh{ z{sE6J4;?s))=gK)l$%|8S~< z@bB{lJ9`%kBfH-s!zO zUD{WP$J`$K`7LT)diSpnu^`l*CosaO#(4ZVL%=oKDJMF-VwMj^);qilXSg_ku7R8-?~}YXfHSWjp{;zZ^N-zv@0jb{*`jOT#XQ|2YtAyWsziC&vlk6Kz0_zDz~ zcOiKyPa-Zkk$I^#G$SPv{ByoJ?l+p}Zj-%C_<@<4osa2wt>&zgaCtE|_~k=_Pt)g&!QN$r-uR1$!ZK+* zgEPR0SI6vx5=s*zTx}p}6-Kqm5FXLVY`Hz}x+2R!Z0Q zVbA_#7&()V!_03)1L46<^Di&>oB+IR1ZJ0FP$``wz92e>tR)GYgM)l{`xkjgj>l=p z&)I{SS9~XYtAm)xHl%`S6iE9~!AwM3*hn@JA#7Uh(_>d63Y5oPQ*wA=`8xNcDmHu>gP7F=s!op@)rKtILlQBS$_XektY+55NpNB>bF_ za^oyGEeayZQ-g|ZW0f*`W3oY@Ws*>F?UEzUD?*!|-y~Bjndj8+@{s(jhIT|hsgMkZ zL%$`46=2VUb4e0ebCapS!HP4BPTV|P%A*<9G%KZNo5rl=;O}@_m-i;K>96YlfGh$V_-n2t2%~ z5C%&w2$eaj*@Z;Qup2_);?Ug0E} zqmf$w9!aATyFL*vW$~HKD)0W}m)w{cYKXa3rPdDD_j;(AH=GusT$qYb5fp9INb)X9 z?RP}COq$p^kke4o%;dtbJOv@9|=_OzW9y5#VkwFZF0W-i~TVZ()7XIol=+o`Ok9vinBExPAA*)X~{VT z=B}@dVupMJ)Q%wUPa$+~!V2BB)X?ymA?j|zb=-M+`@Xs{!(!7!ghr#1NH8(%ir1by zld!{tk%5k1E?=;s{gILOsa2m9tK-ANB9|Tns<5WpbX19FI}nv*SEIYkR|0KX9^zYH zgIc~EMxA_mjU#n^m!ordF67r}dK>kp(R?TIRc1=XjeGyh>e~QSW?s%>xUG}IcZg)c3 z-YX(>q~s^KTzUT9BEz5blyJ;H6usY64rX;<>1K+%pL9xNdxH~zSTN!hqU+?&VHy6qJt;ly_BScq2zkx98!M}vhX1&11WI!(a7kc z6o0~v`~)x`Pu%~$2h-%RT;^5i@$8qvsIXi*2mUOjS({SDO(4!qCvh|wluUFbirzLU9AI&pE z$@aXl^(-Sv&(+t~2Ez8ZPGi9npi!7pX2gssxq-w?`PYkvAIW)RP&cOFCSboz!B1^+ zQR^Bp%kle->yD~n=!rP!AjhEWv-rU{WZlY8#-J3lK8-JN(%QM4x}x z83Z*+1DrDN-*pGPsioC2YdP`uoV#MtZHoMQh*dcfcwIXy=bk+)s!=vpzE4G76MfuZ zUJGGMwhsyppvzm!@r1zel_2!ws9e7?xYdlTkfN}R_(>+$jkZ@;O0uYs+|ca4`kB~e zhE=9g!&hm)4scVh;cI5S%zALtL|rLATtKAzi*;=FhX}-V!&Hkgd}|ZnM+IgNuH~2H z?}1f*D+Ajpo7uS#`&!4n5wFapv~1x1tWIS}5o7+WjsEJG_jjCO#-#R2ladFNxyjzs zGr`@k--WYLS#wYF0`&spmw9Bz!OCW%bV*pc7A4^1=H&(A(r2x8`^4d(Rf|s{O_;jH z&AqD=v2x1t1IamyHx&cpu|8MH3xFseWgyoe8PKr^s(s)b<$L>al+8uE{=gCZC+J(k zkA#1L%tfnWuyjv8bbT)-D3LO657W=2nH#eSt4blTV7HUweZ5T^12OG73 zHaAM7nG>oVO*lGQGLmq>Jgp(gD%5J2B`^DU6jf<87$#sP(y)jD@QNIlEQZNh zXf@2^KT_+(!$%Xw=WE_?nZmb~?ay5tUAxy*~juvIOt-C_$M=hg7IE)b9na?RZRSLcJj>KV&Qn&l%o@9f${QXn@M? zytOMEHcgmZlY96|L}BqlGpdo$i)rBs#XMLt@95D}tq!i<>@f7XU=;6ShjsUh-5Tbj z*W@;i78z5~ckOwFj1akukZ`WpO7)D(#26FzSwu&X)g7B!Tm#F;wl?{UEwP2U^yFQ; z3XMRm-rj0BeU$XNWSXmWz&NUG#rQGOr4%KucGlVTf_4tXIJ2~2R$*2F+8mzKh6UI6 zL6R@6Fl?r$2g2`MW+%wCovyjP0?H4i-N_XSD$Zk%1p0b!_9H#Te9b9(4X{q`7WSV! zsOf4^^X_iIS7};T5i%yck3SZjZjIWlOWD_zLWG8b_*Y!DaV?M4>J%J`fssms471Y> z=o&vVK^(5f5soC!0$CZcxyTN1Ej4bGd)WTdj$x3M^nJ!~!L=4!n+}SlS~jZ(qwN$o zl?XjQ>~@H_Bk?sY_7b|M61jQIxdnbN#AQ4{Y5^^51b@bXB(k$c{Di!VhMybY>-8e% zm9W9}3H{hxbDN^n&*Eh1HyOy%)6jz5%DGZ$_cf8HNWDITL$5vUW`(L1u~Z$lU6~@+ z(Gxfe4x9-tyRYMcK~@NpVQi`tz4;gMy%DsubhHq=lbV_Ppf!N|7ef7L>&m&XjhR~= z%ai3Z7AKCBaOah2gQ5gn&WU=9iq81h`aM8Fh@I6mT1maVTm_oLQ9Il~c?#sj!dDi* ztNuNyFU;F!9%06c;W?+n^YyA1XMU+)!Pi=x zO#pNjw{LZb9cVdNlsap-DcdzB7(I4KTUnp8X;sdpeyYh1-nUu}g`42U%1!&?YUuHH z+Ms8$B6D<0C6kqSC11nS!`$9+{;)DL`SALh8xP)~E@QhZ7dfR}@&K8_4}YUJg9NF8&cP60h$T zkGQ4bw7UWz8Zw(~sZ*K4Aa>ly4vv-|?@3G*PaPM~esqVKHj5@EGZ>m7lnfpD^a(ILo>$IX|ty)OFFWRu0eW9eck8v_DlzF4=G8 zEzM_-cOw<=zj{7K)m3V9!eDaX^ULGbmmTef28iz1Mhbh5ws0P3nX50b@-~HjRpjIM z1;QDI`&9#yzXi9FKr3nBk~~Q4>qJD$U#3g$zJmOKSm!@hp+x!JG1$y78!9NzXcaTO zNa82HoEbnm-U-uPP@EK*7TL|xTrN?xdd-vNk&LXS7mTB;%*Y2P5)XbTx(2aVzB?>I zN$j24?Ycw@ilJUmIrcc;t67tF<8(oueG>t;+q|;c)2Q{Q&gc!@R?Dhop$5%OKF~=E zqL~Bm$^lF1X|Xh195?OJnPQrF?vI6=1@w#KQ^yrGy3J!IcLU({vNDzia}sxhKqs_U zD=P<%?<8$Zy_GZ@ZYuK&SVrzJFc{%?=s(4{!$we|kW9Y*xFGFHoqKWa2QYPsCvo-3!LLlpc#QCbc$9*n9SKFfY_i2G$JA48JFb9;`=l0EvTIzvDEWj$M>;(rrr1@j?|mHm8xVf*BhMNq_u`p`enmUOw(N&cz>o6aIV_8@8FvQjLfWwA%o~k0<#8 z(%5shMl&AV`AI!3%z5>f(p+dQ^+rsHs%U+HAk+JIX`Q7e)Ks~z=MX7ECl&BR57y6 zC8SngU=^t##6F2~Hp2GIOAP4`7c#4o*VPl$YVStHg zmTh!z;ist13e5caSb!3%?zrYSw#~R&%k`L5J;N^?*(=~asen=A!5QFlfm^xiNg~re ziE;Nw)L4Ow=3?-~qNhgbZ{sAl9YeLb2`5i6hM5*;FT;+XnjNL3>)V?I0jIbz5on&w zTfLCFZePyW4HCwOFoc(SD?~wYZLMN>8I$Gs6M`x7Yl~!rgIz^tY9@@B-<1j$$DWgY zi0}BTJX>@lG694E)~EKH{Q9m-RQ@v>X=1>sY7-Y~@vx!PjQu%Z#gQo#5X2xOh^xO+~#Xiy+D62D`fbTLnshHKn^(nw5)vo_p)Ma zWofUdr*H8G&1GES56(o}Jr(tO7-&#QC{$67IbE?f zSrNk)Xz~FArHDt0P%)@jS<$6a!e#bre8H8TRO4-$r& zYLp?9(yR~wwovC;M17DFMJe}rnZe~*-w*+-xNrJY7E7s2x}3gcUC6j>oF8|fKGU|i z>Bg;PN_&W$d#i)1rB66*pzdHHM6soKFln7jo!w+D>?~``L^E2jV1%g)-3XNCk4!lk z)v~xMIaXx#Qs$7QS9*kH4-YKR>H(8RmUUxQSk4Q*@Gd2WD}q<4;C)54n(69c(ay#@ z$_j#~`9=A@qyz1E(TnZ0b7kXMOS}*MBK7v>W1RWh8aH>Sq|zUE+k3$O$r}HCx82Uj z-rmg8lL~!SAKJjs6MSCbE;LteM`^b1_M-Ef~M)44CH(%k0B4XXSzQgWP?Mq z;G+|=`uvEChewgF^8(-D+%p!T4=|n?LBY_dq2-y~^K4IXX8?y4gv}tFH|AqSA4S-2 z1Cu~Hk31^n7QgD?Dw6C!GTdpTuf_U9lJaw z_u%A2@k)na`q9KhF5TI}@S^GKQjd-g8RwDq?ly(z22H)#m{+LaxN8@KcANGp_{U`hVX1Ztye9WW| z8i8rbc5SeBjNB_0$4ea$$L(+HA*-!&v3YnnsKo^ z4dY=&zVzR}536<2hs!OO+Q~k1_tc}Umq@gZFe{_o zNTdLnXITG|h7sKx3lJ%w_vz@9Hg(+dd=l-sXl9zwAchU^SNpDPw3vi{biCs@aW zYdRut%NH5XF9u=r6Qr-&qT{N&DE6{^+~)#eTd;8{9k}3qFyYIM%7K`c$XAPGzk6po zUC7p8eR+iW&E5tGuQy8d^>6N7-(i;bcXkM|Qok97c!Rr@dtLg#v9#CX;&S?Xr*vYt z@BfOYA2GZ3O7aBxJ=snzRbeeYj3xcy{$wS9K~RAHc8RtBCY)${Z~TL??LvbeE%TJ{5Sefa@F4z z{FPhto453*2z+>?e_}EHUDIDr(f+d^nm@Yxe=u|Y4*% zbxr%dA^KCQX#NHM`}XMXD*oz-e=mN23MBKtRQ$(M_;>tY9l?J#Kc3~^_5X|W{yY4y hL-Riksr2=K^&PSj;2+)kL5lkLbbX8z5cc1({|C3huf+fW literal 0 HcmV?d00001