From ef24a718c86bd9c4ad5b9b15e2f64703bc3ffc51 Mon Sep 17 00:00:00 2001 From: Princy Parsana Date: Mon, 22 Apr 2019 12:58:08 -0400 Subject: [PATCH] combine fig2 and fig3 --- publication_figures/fig2.pdf | Bin 13763 -> 23247 bytes publication_figures/src/makefig2.R | 104 +++++++++++++++++++++++++++-- 2 files changed, 100 insertions(+), 4 deletions(-) mode change 100755 => 100644 publication_figures/fig2.pdf mode change 100755 => 100644 publication_figures/src/makefig2.R diff --git a/publication_figures/fig2.pdf b/publication_figures/fig2.pdf old mode 100755 new mode 100644 index 6fb80474f0f025b63f2be6513283116213138043..caa87444b8036d32696162755db9f5c421da7987 GIT binary patch delta 19432 zcmZs?cUV(R^zSPoAkupgA_x+SgpQPe0@9lx(tGbkY6J{|bdVBy7m!{ALJuHKLI(vz zszB%^5=u}M?)LqibMASbd;iGJo>{Zk%$`iXYiG?T{ce-Z1j&*Q#6_e<#3Up|#Q5$A zX?xi}wefS{;C(71Dl8%`{C{jIW=ljsS>h4m8q9!rMnDC(NMavqF2J4%NG9LC^JdF0 zGyq?lpZ6)sPKuGh>UMXsx@oae9HqAORF7_LQq>X5A&t}HU@Ga7G3n%r^P|D5x>to& zXLWa?SMT+$-HX0H8N6P&lGZZ%YJbgd(?0U`{g!geN{`H+&7b|Go!{$Wm7V#Di9-P$ zGpif2o1#{m@SiF95d0&#!&Br}t`)wtpme9-k3NLFcz&6+RqY`6Hn~pU1CH!^xw`-O z>TLV^IJZ9@9x7r{>+1@4tNV?Wx_5O}eEny!KIzY&Gq3CUt8db5S6HUiKWD=1mal0j zd%UjyY{VCBpFi%o7(WrZ+IYO^!d`Iwmzru481SrUD_jnInS>mH&A7DZeCt1Hjdt?t zhC;7ZFVB8m3SC{EBq?uuGibzHZYF&#`T7F6gAChwnf~zAFh_@6&8vz-$a^uX`Z(q7 zpjUpRMzh75H>P)83nV(i#e|fvYNvm|T8%=QcNmsV7A)0cDUImbj?OJ22a9{<>VNS= zUTygjO~YSU>UFvO(5|h$YjNDLeV$|Nw$`kv}h&1Fud;k}1e(z!PWnl;D<)jePGPsgK# zE*9X3;xJ2DnO`qD>_TwTBbUW%Lo%&pE{kAJ4KV+z)(FSrYbzqTBVwlT`%iz|HJ@^h zAfvgsCCyRaW=~Dib9YToO^M$qh zMV~ABmD@CaXdim_7R#C==jhD?qvO3L5YEx?_^_1>3W(ZJ58$Iiw(SN!BqBb$D7Ihu z$&+tu_hU@u2_Zb*laJbBZd3($D@@scQYO%Kk02ApG)9mqA_yV0tQXAQ1cNHk@!a*$ z==k#qt~ZC?1P+!!d6hgmo>AwmVKmt*93|7w_Q19n=L5%2M9R4o5On_Pk8q35emD{I z7bO!-b1EUE82Xr!Ne@FEO~!iUP0;_SE$iJfOeuQT)iXyqO35fXKKFV*f33KNWmq{P z{>$&k*xMnIWU-c?=p@^cox}eUGtp{yh(^RW8cJT(d+V9PrBa*$`^1k93dqQUsj}#H zvA5pbiw>^aK*OW<-pD*@nU+3yB-`EhvAxO&j#3i8Hs@(ffvmaK8h^Os%$lvO{_L8q zhWk%=DqgsZ=@LbSk8-SGiKFoYg!{x=4#OB3`sLQMtP)lBZyeenug{S}Oqw($gM<_ll591=kK>B_IajDO_+HBhcpoL;ACJ2uWa4jTtYMIm%eyyIPB;u*znQJ#j<2z&o^#$(T%{xrN0pNCFEbv+(oQR_+u|@ zzQ@t)B$Fkj#o)@^M=PTwSL%_(V_%NGpO1U(p`zPHie8c`BhiZu(v&j;ZO(aD!K-XN zkCq?C^<1w!8p6K+TVCX~YUAz|j5IXCVCPG`t`BS{-J!)6i{f~9IH^tNZe8JX7o!5N zr$e%j^uUmD)@msTNkKkS_UBSP8)gL`S^6N`jZn`&aVE8jD*X> z%g0Q|Bhe0qcY?DVe|_BZ)}kPyD2$pjQBPGQd%TDG{mmqbS|X2Ju%`Nf zf|Ix;b-Al>YAG^{_iEn>0gN4|Nq}ZKy4|6%?>gn9IfzqG(!kF>3%( zq@JUsrNC;dJ@}wrTSJcYmCD)&umT;|E3LAufrpAK9Z^+efjon3uV}?zZTn{YF+Lj( z%Z|>fo_d?(2S=M_GWDv5oQk?6^7$YTB;Pu&*Hyu-b@tf@z4eJ&#b#2j{G-JKTtJr| zRa?V*)5F6?{<0OiF(KF&4n&?V`wzzxS|oQc{v&T>B-hzW6H7hQrNEV!Pb@TiAfk2( z%ltkC4zEuQzxJ#4H@~In_Sa1AqiBn;Y(&y6sHeEWy~5?&q®8`W}}K>14@{UQIuNQENJSkGQP+tq&KRw^ILy^V(wah&nw$@9 zG>7*~V-)nAh+*8422;ZkQB*h%6V(VPEzXV{r*?qY#rSA?WOk;iDw4jAfO0?rabe$9 zh6(ppbwH%y@d(w_>?b;)SP8>?D964RMnNO&5KUN#)d8iPeawg67u~RogZ`*pOd;vx z0aLxg0ri#fI-tXy+w6JJ`}_`~=m>ZL3PHlQx`$Ed&xbmeF zX%|;g438i&tVrJ_AIkrd28$*P&Ht~Zs`7~(Fg0Jp?!t6YQ{MXEq_JD>V5*?QeCW$bINBbQGA%s(4mt&H zN+FT?l?XMy;GnwYkq`YLzV;4^_FdhA#;7`}f=Z*85hNYw3ZR{dKpYzVQ5a2FM&yoB zm>Kk(h@cNx9ETPCTGj#mtXAW`HBmxXgk||YAG0cScunT6s@EKvgk}9p$Nn| ze|z_e8%)LKi4ke|5|Ixrdn-S$0}@{vE(w+1hL6G2n#IKTtpTIy6<~~2$W4hxo$xN9 zmG*8<@LFgLEA{EwR;5H2Ku6B>aZ|XEJSYYetG76;^zI#PMgnvItI3}!4MW!|acQZ5 zz9^Jz^`T3|3GdLjd;-vQr@C~Y~E9b_e5|`KW z-N?Yl4+Gvzp@cRlwBPFCyNhxSKyT0g8XGb^pX$E%p=#k2esORqIy3e3hte_jTU9d{ z?m&E%Rs3rrI^vq@@zu}e?VdVShh?f4ZktDAmvGrLO8I}kR0x;(v|Y8@#uV``42r=; zTYqP`$2y#))z=Xk?R&y@WMIT=J|USgUM?Z)PvCTQN0ZTHz zYJJnbnIQzEPJI39TPdYi`w5LgC*piskcaymx(frYQYKki17VSft^1k}WUA$Jhrif2 z#w8dO=1=8Ma?UH75Y&Ayx^l@*QE=-}Z#k+WytKc1^k|%W^T6nQV!MbgOeR)+rGweM z`C*$i++QoS{h?@WvR(QylS$K^E`Of}q9(=`!@wN(!5S8p^G14Fik6$SuzmdudS6oT z@>88zLrngZy`fHgwiu5%NzMuP#vN82dE@HMN;2M*7KoFW)2FTblc7GpODGq~a_}Uu zTV+pVVIK!~NMKLX=)}=n#b*(L7=@f!cG%MaIG|x_)lhWAY&0{bmyv}XBT{mZzz*YJ zpBau2tdm6BDen#fGPjpP3^PE=jBzBfd{+y=yutW_lSpIR6hJjw-`_#idgIt(1~2y| zaU5k63Lq93C(y*W2cfFWzEA)~%l4#@n0^W3qSC+blRPwHU<{^e>GK3rsXtSIs|Mnx z1a&|(66t{Io8AjbR#hNP06MJV?*LOxnM=WTs_v>?%EE-fpE;5eTtmrm6>6XZG+`I= zE-YCuLk?DixQ9^~*#N5K+7h8Lu@1T*O#Ui5!jDcCCskD;iTW&xh{y-Ks%Ezvc24 z!6mafiLcsWrmu{rN+W?VVYO`mgezItQv+-qZ3kmRx{w8RG{n{-qS^8*yue9} z@45rQR3`a!m0ny_+57ju9brW3E7WYjNi45LbwH`CNo8=YJm^g}d?4<&)~5g*dI#|x z2qJhu^U==~z*N#Q`u`$bGY~mBkX^YaPjIl@6YCiXbcA!R$Tor1f)40mxs58w`R6Mj z04>f8+bIYhzJpSqNW*SPLR+|~zI-`CM?56!6|`V;@fJ`&1DfT*j0q?XYC)D77a=>q5ar`!MR-fb!W?Ejtu{+IoS_kYHJffF#Q zs=Md)|EG#m-v5dH|8-H=+J8-k#;PEG|Megy@P1~v1QH`93x~}l01G9W+pdIr5etzt-dq^0I$9w;e3682v>P5j%#6~H8ek-Y-D!hJbCM%1r`DGW(#y|d> z)Fb%x&?AYb3Qu(7^Z7iNYOSpJIrAcX5d`ZtykYHbwSX#2wcrrPjs_`NeuP;~)$*yi zu1s=Xq5M8Ak6)VX)r}`jqJ)$vv~S5%pS*T`2`#nChgn5EwXX@x!uP<=|6B@DY5PCE z-0k}N>*#D$==xRA)#YUH`=pWaLi^zaWZ+3p+u;j7n-rU+q}8mU?u`aGXkQOr5PtieSU1QznZjUr6FKpflW=#D)6=g3@RqX7j_hMV8 zW00FlabNTNOXc!Zd9`Nk;>mYqtm14#zU)$@-9*(g3tN-11qL!>##1Boipf50TPJYM zIAeQq;>W9wZRI0H%XdQ0L#8{Y;8HdnUoE60H&C(>cF!l7+R2VHN)N?Yi%V#&K4&3# zF2(}d?!kkOkN#d73CW{g!D{3>e){ZxAtSr~6$t-xYiM+KzxegZ!C_|f)YH8rFL~6} z-e6~o9_8fK@>{pdjYS^op47u{tFMEaGe6?}O?@8rT$ZoCCa(vnUAK4AufjiUZMHzp z^{@e*q9IV#poK%E5k9pOEwvB3_k|RmggqU-9emug>r!5x4?5QAF*iE0gFiQYFwL1> z?66s;oGJm*GHI_F_1a9GQmix4T-OlS2>m9_kjRsLqk4uR**Ue{GQ5pZf|Xl@Ui&)P zxp6>OPc~O2;(!V^eo&Vp0}pdxa7$E)pbxLA9#9j{?ypx5(kLFbtEFBn4`ZT9X6X!&AX`_7W7Tf-;1&>d6CcjWGOE<5nf@_0sB0IP z7-FZH%2w`PIaaQ^zwSb-J$KK{FOn*jXApa<`LHfUJXZn@b%NBuJqx8=coKPZwep1Z zO4J)vO?)73;y#OY)n}@7dqd$hXDzN$(r3ah#Eug+Q})SUY}ET>8FwgIMEfp7s#U_v*O6zR8d>M!srWXlCZdMznYe61h&!jw7Z!C-Rz5}Q z;@2C}B-~+&U-Ly1;rZa1^2mw@_Z?^pv#^RYNNnC6{%h(l+*Hu=FvH|Kd~r%DU!>KE zJ=9DVwwC0BZG9=#ze0lA>g>4N@2Mh^gp4~q5O`jt5OlEpi@P|`gWr@e-O?c*tUz*m zE($#RFc4Re_lpNhYU+*yg?yET=@tt;!1l^vaRsEO+ADN$rW2hNAUBA>pP1kZYKB^d zgBzKbfcBp~aRsPKz(th}0*|rS$H|}4vvFey`D+@7gZ0bY5%^nfk}%yHsL@MOn*$zf zwjv1k`o6vs0&kCYKwm~y4rs5a{_H-)gmSJ~tTgGTDXlPlcCc7+A~pny$5r_#uzyPw z3hs0k*@@Dv$@!IxCo)NoLKQxf7=mCAS=&&P3F_J_;7UHVZv15io=G5iZd2}VAuy0e zS8@oJOCVwmCtuvLVX;DDaPeU@*zjXHa0vP^dKqR~3c(WI9Snd?enxko3aRP+VUORl zu0ya6Arwo#nRuAP{qLAA@(CZBQ8IcNKo%3Ti<9Rz*zXz*uDp(lTGd{$dvjQVcUPy*#OsGd zrcwTCL-iOlh-TtzWuI(dd9V(y zF$8O4w2tdmM!@mF-P=kSTdc(YjTpUT)b^g(>>7KTYL|=3+~X)el@a>W z5bU7h1j2VLOlPwzEsaz38zxj^hs+=|> zn%f2@0?0~jcF}9f+MT|(jksa4qI*_@$veN$b4(cy){9kJ&p~DW4z~}Q*$g!`*^3+v zX1M{5c%-zl+SvBPVpFcPPyUsaiSHWjN;33y8LG!^;kw>)NX(7;P@CE-^!Y(zrqC%v z8eFcF`e((7KRz4^LXjK2F#yX68f!f7P%^pmqL4Z7-b^du7;@-)bhr^%oX`Z_nyVABmB9eDFFojX1(^T6?D*f zuW?QjcNm|sni^>7Y7MS_HT!7j4aR5Y5RUJbew|4Phb&Us1uvFxh6cB{%Kbi{glvt~ z8*ho5fjRnrc8R~}-v9PdjiPS9ArLYLj;r>!{1%Qi_Yy+loBd&(Ke!&8K(3<-ysML5 z+|6d{O+_jdd#zp_Znk!wbYHl*ldQhza}~-Rjk4l>d$xPo-Dy52MBi7%6TU4%4?oYu zu6DuvFUE5Nv6io|sxQ-I^S!x}jCxG%C6=MOa>qxQOFZ)HR~;%pFrL|Jz=v&)!)k=E znJv26_>Q{@Xy)TY24U`WmRoW!Z}msz)ZfMAD#~OegrYtg8N|8M)&_Zg{!(YYVIqlA ziez&?{$fywjL&NJ0H?}E5azI^oUCy>ZWL=4Sn9|rDafG zFw@#|Gsq=`Z6mfk^$L~7OuxKyH2GZ9{CjryiKoRJqiK`)^so_&|6BbH)^dk3787Otkj;{-2HtD^3EK+sG^ z{>t|}_|;{Fs*8=8)EyCa77sNgn=)D&v}27~2c>f;ASD`fJ`-CJqt1G#QbuBZ4Xk&4 zEXPQP8ZXvG&&*9umolB3;JZZQ-v%DU^S_r?_L+9*_}Vx3A|W2~DiBCLQ4Go+zuRgK zSZ8;+)ym^nn2_bO0@U-T?!8`vjW6RYnltV!sLAvM)xq)E&T5Trmvrw;f_R4r|6xX`naRE4n zl7Gb|XOS8UC$sBLio?f$j|^OAcbM^554}GeV}bfQ?ZkwwG1%2X8F3dz2dP?Z^Q`DC zevAzihX99}3c%T$P$q1OM70HQW0lP7?C&D|k3kPocG-?~QfUBRiHe@a-qZSd{`a`=Amtj$8(JTHQ~al$x|Kc zm`^F8z%xS>$CTxyok7t}qZTMJk_&1qOFJ|wNT@=;j~W}j|;L*y5sTC z-jdv)$r_v`6IPH-+rBd>)4LR|jVe#yeBBo0wKkBsG#LCm)tP%@u=BQF8q9WGZoK9( zF!he-y&XYJVhkD!BF8IzCh)n$*gz4>q5%dvtcwGkBxGe1Gt9KT-o?%?7x;f%>k@Bv zWdS{K7iyzg&7v1%n6NI6d`gg&#ZlENL%^$X!w6*M`IC%E$jV^OyLX=Ja(%5p&-0_x z)4&y5pioi3AguN7XkE@gB@ePv&fp5C!_rmXflm!y4!B0O2kC0BymxzaTrfvE+q&-} zcIc(xm7j!Cs>TitTylwQ$%}Vx_~UyQyF+yyS511uJ3{4H7eWPdegp5UujPJ&*n4bD zK8W1QSiLfTH0iZklKc7+l)L{8&)l+l1#$0Q?N)tsbuK5l4VcW}Lzkz_sGdiUqr~Cb zkD#qgV%ML~@A*58voCHuD~?|Tw;fKmw#qH2vTk@Lb$1;$Fb*dypYJSYO5Pwt>XP~0 z7%Jr=8oSBnHf2d_5@r`{mU|Im&-i2Bz$xM#$qxcgMg6Fg*A6^CvJ46i%|F2}j-#>i z*Zy++mtGZ!)MLjRul@JfLl%H%YLX84^PEQd^v;VJL9 z?X12knrKZQ z4a({j0s4lviwSJp+O%v$NNx7l&h-HrJ{V2M>gRMcYuz1_g;ph{Xj>97e zuMQi%O_^=NC0V-vajJ$ClZYd+VW& z)r|hf6O7M+C@@RHvMwh~)Io>kT^a}uzA6Zz(Im&A62I?*hY&nu2q^fTV1PBr!3cxw$bWn)WADk;(x&Axra>Wkz0qgk*u|QP# z(qn-lsrU;VX5LdAPb~daLg^+zJAw+Cu^mB~n-Mgyte{D^*3 zJ%~yyc;6h=DytSn932$BusV1tUNW*|fojEO_@yYJ#3D-K@yLm1mz_bA?!`wX>wj&D ze?V5oEu#iI^MqO39v#D1WW;Cgb^?mXt-spr6kFo+gbCXiLoeMCA82Ai`JfJvJqgF@ zHuzDuXxc~2x?DpJtpy4!wY>(MHriXDMsJ=tC7@c*DR!*@xy(EL40P_jUc0w0*Xsi; zu3EPzpZ3<-_lkgfmf{PaAC#pCJch8`c8Uo!QKn10gy#l%fpx3={f{AFT}xBIMF{6c z4`gK{Au0yYX^I#Ou%c*ucmXxHzk`eqY?&IQ#=0js#o#|3g)m|18h->FgA9>dZcNxi zsxQg`$16HFMGE8b)YDIwsj&u^lcD)`|2~ISqe7%V@Buo9(lcSd>=bSF4RlJjZO6a^ z4Rx*xzB+fvVmI*H z9cHJ1wpL&L26Z~XJ(4_@t>8=nWY@nzcAHkKcCDxNSdqD_L+n>O;Tx-em{CtoLYc#! zYcDVcx*g4fqBW|^IR52>8lvH+*ZTVJu_7(4a`vJ&-vQkf^0=XMn4NZ$;*$$s2S`Hn zMtquUUm;7FLrB^2IhWWCO|to-jgz}!&!+|uR=EVmnt5sY1Q)qZlReis>hk`k!j;8B ztYyz1Tz&L*&ozqacLiIoi@xsE1kE4(qGdWQQs(QB*O=Yv@%-bYGhVo^>`KQ5cx4xr zxiL1OK1Mu<^A0uG>DRpu@h{bt(|){+6vX!;_)xXE)hT4`#b?}+H`BY#oxca< zb>MxAb;xa3c0qmEG`#bfU&sMv;M%Zf$sa{~5Yx zi3+~HnvJKbw*X40!*XV6TF%hY->nnnQoLzAp+*uUvoidaPndReHT-5eLWWpS(072w9)XFblidWOVu-WV%`^TRz|e3$JZGuFvi_pf0Q7YSNn*Xd-c`qAVkfj z()xA(JUXGKlfD?@d)H18_WpgozomaIdE@!B+Z{4db?F3-Kb59t;(4DCr`}TxUe>={ z;rg?(*sAYdOK$5-cBo{PTATi+E;Kbk$3%c|`Ik#d%}kafJT04rHGXc+xLMLoLG_#a z0~6N>c8ko$gxW41$#`KR7NM9%W?HFV{^>eATP)J?hWSE?`V-_Kb=9V_wxqm0p5v)g z{QLIF^i40w@XEX26=a*M>ZLOn9KwWhv{>8}ui}*Nr7uY|rcM?UvDsRrNZU;h+qZ9X zKX!}Nmr^PjhBru18O{xH@C`I&kR*wj`MVX=+qxAvZV|iMe9tK7*aSODNDW1%ljgi= z+S{YW%p>AHo%f}3f2oq*Xb5$Es`A1g>!uyxz4bz%ZfNca4a-+~ZJMxJAMRnCzOS6K zxJ!22xuU*#A$PsTk#WP0UEpE(Y?F-j+qf|eLSe@lRyg@5;}nr5LY2H{1#Tqd)44ks z3Sn*uw&@fC6?;eyC38wG%?tt8gqpYCrx-0UzfZ}j!aF;}%zW?bs|i&`+}AWn5^g>PgzVD=u5ye zJ4u4d9ed_J0;QKQ1-i|!=>bRssf=|E9nML{`XSj(U^I-0KRy%KDK{su31iYX!lsW{ zD+4ye8X6^Tun8~tBk)TEb~d72grDDjsv;y7fb-20*ePzlp%A_`8gspLL#b!V4*To) z;QfZWQeb}X6e$%sbGV{2W$%rWNUrvcFo3HGu?-UuA+aY(h!s>^AW-_#w6j1VJn%U> zYK=lT=jaLlcLJcP%*#*CYo%)Lz<6}{JN-X>`sNu1Y^|6}TWi%gzXGnU75GA6SNt4K zP0YN$paS@u9uW2E2>g^aU*t)WP!mQMINEyZy^jDTt75a=;~nZS;JrxutC{u~eg zPg#!QIY(h4g3f8v1pRrx>4kYDN!G$>c#y?0geunBHH0cLYW%X+$ohQbBqeigCFiCd zyRf2iGnif2R3Q5od6f9{)3UqP@Bz|Td6EQX`%r^B%%TERD-t9LS<@#0pPch$`0%@7 zGy;fSVrKovNn;Ak)<|bjvHh^=V)g2sFd8QDvE48LrV{py@3v+TB}w=QslbF25>r0< z^keM~Gt`_QB3i)ur#BHfs4nbHKZS5m+f8CBqyO9Z!bBW8 zqn;uN=n2UadG1>QH9m1^DjW2Cu_{5dop-qHz{Je7+8I1>(N6$8=RnMC&3(JwKFImSwnjM{wz`6@YS}A zgSM;?;+tzY^)bThT&EVf|S!yV8qd9pJzdz@LuQVCz@7z+UhNXKwe-ugi=S(# z%`Lp#?0HKEns2!ONLxKBV8Q>8Xx{SA#acDis_W|LlInUP|Igs0gKUHagCt*A>t0YX znEP*Le`n7<%dS61v%A|k*wLxbjIEo+FfA2d!mRlJkcFQef3$#Nwg``=E(j?OoaqwVY!hdL;E_Rmk|SX-tMvWW-DPoHL2GNYAn^@3GL~FuQKqTdAHsax~r3N!&MpwsvRCApIWKwSqTafbou)x zOXP7kg=Ogv#F~;157SOI?D1=u*4>EQe-@qxC7ax8cZwGrm}_Cwtme>^jxDzxc2(r+ zrWEnJMF?xui8L<%TAa;L(at;gnZnWFp#n;{+JUzyp!+PWSD!wHKvzZVr5+QU#e_)H zT+P?uE_268`-rk;`%S-X-<1Bir1(-s&HS6grlTwl&y@MXdiDf4=1fKEwP!8malJei zP;E>3)1^haYJ?<=+gIB?BDF&jY!7?!wK!LBmj<=nr@WS3;DqXA(oF3U4X*EUp7|kWZ`E{5}!P^Ii){+O1oQj8Nj->FhiwNkQ zPJbN`#X@pO0M1!eXlF+S^xS4<6Ql5EcDTesR%FT_mt`He2W&M@*TQ{q$$Q^!n*SgCaOWG+8-f-L9nPxngJELWWYZ36DGaP6jj^VY8hB* zHmNZ6Q&cS+p{nwNG!XNqdjdej9xMacm=<9G&YFRR1ZL1qtON4? z`~koSOnM@qjun3sp|=*$wZ6D4g)6;8sivFn|6jAG%~=Z&~0klTuWTpLG@JfNZ;z z)|ZF5K5q35b1BcMOT$WIn}I3hD{(!~htgKB%?+1y$+<7BqDw-927CdKKtiS>-Vaw% zqv8Ti>LWMLhmM5D6$89eK82zRm}=KOIp`5Ex{7x1|Glfe(ejgw71*J!s084OXn}tl z({QT2T0h(m0HuPO19jy#tJ~iVXO?y`3hWs=|5LYZv#Q!|=`)YlYfB}Oh>`hs85QX0 zvmk3hbm+avA91;l_7ZvO6H}t25Zj)Z_baZqAeDE2O(Wl= zM87SWdB#Q6Y;pg%4(JNg%L0p$U>lcYfvHKcrbz&cZ^oP%IOqE{JFLWCU#w2XMRlxS z%#!(qq*M=Jr-uAi@<=EPe}|n;|5RO%fIHW%??>gfhF-#7qt0xOZ|XKr4r$e2?UPzA zIy`*8tuE0Jx%wTu(64#Gld9F1Sy_&s3%yXz`zB*%XX_}q`1Pj;jI}$utYyf9YIOyx zN59@8=N{anx`B+6Sp>EG{p9iUk(KeV`S_V-$j{W&>dmpGaQT*@_x96}_MQI#sbqhN z0Lk-)F8E;z8~bliXz(J}&P(@2Q0V0s%OI9V$#pER)*|0O9aj`Wp{HTJjzPQMH@g)c zp#rKlkH4=u{50Fye{->RH&T;j!K$3o(p~UyV_ACr`sW4}*sT5facf^RYiadD_FB*5 z{}-7u=^IlX`&n6I==zg0XU%sma^3GXobaRbgdKfYspDpz$IBRs`}QvA_(($936or5 z3NX0PkZ++@iP?G{Qf@-Lv~yccXOJg)p1(`lA(ze8&2FpA6Oj+i<`^yxNTpa!Zk_Op zeT$IBpRXBL&NEQ3ZW+73j1jz(ULq6zAP=+XD4k+1YQVb~3kVPeY1@{@-}-i0qGgVC7Z(Nw?C zl>-s(5}Vr1kARjotH+n&cKRymPxh~jpb#W2f<<7jDU`~=k|7QNrex-*3bIcxmxFN}%L76-20%1vP5_zXVj6%^*sBJZS9_Vq%s7r0 z((&GB_zO`S4}3YB&%iDuC1zPEG9vCGIva!Z%sf5$%lOTqr~$ z0isGHiVf2n@+AU_JKkCkcf9!T?%-%biX|`1n^UX`$ax9CeZ;Q-YO{d{xW!n!0*G7T zH9!~r3FM-}`2doBF@BDD^T`3gN@bh@tQ54Q0&*XDI;_+l4Fy8NCP>^x1RY^;u#!R& zex?KR>-#($p`VEuddt5^j?0>Ip5aJ}J@y1AeJUG~f!zw&mWADdpAY)uZo5-AAyl98 z?L|P1YtafIqQp819IyLnJD~hxH2@UKl3gqd(-NZx21zGyOBl^n)X=2^dNUipih z)%Z`n0aWQ9Z@fkY2UCG8I~eYL7-xn}{igXJ$!84n#uX890A!f?t-j%i3~QiBhWw|| zH!~~ioy~(<(la=Jd^gN_n^Z=jks0Y+Hja+iYLtRK6L<}#Vp+W(xH4=v;v9OZ0}}0t zZvzJcJ?nN!9st^FYXNs&`Cr=mdh9139M({_VD$+x_7eo*t1Gak*70e=d-u{=9p%5x z{4(T3N9uK$^e7l?avXdiWB7uPsdXJZY|ZW)>eU*k)VgW$Q_g|MzgOx%kUh)Q8~kkMvm8(EI$y>dkM%5C0Tv zW4aY&;rtyzze<`;;&T{YhcE|R(a&)F5@M4D|m8jar2ZeZHE@ zJN)yv-YfWp#?;H$dp#map+k1p(76)mVbK1&!l3dK%Xeb3Rpx%}&+}h3*Mn#$^Ut5}EasCpy3iIet;Rngub6laQCJ{n zrmcvb_%ZK(3P4kB4lf%2p;1}-PTC86Tql_SZ(3A>;o4ZRGct2<2c$_@Lg(P9OMU#`ts_N{z#iH&Jl{0X2@361{?jr#v1 zqq+k)3#VP+qeA`rMAsSR76k#Y@E;-SH4Fwovh|+HM>_IanPRmL5|6Jxg`>iYM2#I!3uEcM>vAFY63*D7- z>hNHtQCNBRu2V~&N3!Wlgg7HZOzR z5p}9PJyXlG(ByQJZ%U4mS}sYA44=nU|Ar(-*xv-tG8lklz7^$1StH{?+q>2t(N9OTU81l>oy)B}Ho znZtuDC8Id7K~6Cv^p#?_s^eQlZneAIg)33Fj(w&c@cP1Ht1`$g8TBM^oLMm=SlP*w zTtrM;nmOD6GtV48Q@_L<9#RO5MLcJQ#ny&inE7VPH|yJe%MqM0eC~F)^=n**s4;Ta zXQBgxdu`H6S^`x@+vFnS6)h;Sr~3pkyW=djwrA{#oTb#vEMrD*lWE^@8Ax7dP*C9;eui4_7cH5a+)<$C+wmp`4diM#>ZGwd8VXHJ**M1QFiVg9BKOWR z`7W>#b&e$EVU?UQ)0e_NrI2_(`;mf`U*OVZ?AbiOK@@pp=%fnygObdFnD|{ z_I{vIa(wLojw&M zNS>E*(;<&#AojL(c1@Uh6JVRXprm#pAcEeJ)j!lA*Nzr36|uE$4wDrzkZjr>c*e3w zY0EL<{%AkU*@2p{{V0)nDq)C&SlsD2EOyQI>{|k}=v_`RYcxZ%;%BnjH87lTK8j5G z{p@4{BZd)S{efs#Ol>aZiz21-$3bU>&$4Lq$S;ayJcaM46)!<*@m|QOaP@vCB*$O2 z=>rY2+FUMnd(jj}+E+c3l{uMQSlSTHsGePB+0-3Xp z)6~toAMPQ;=!46N#4|o_@^W@vrdbXv-%19f5?OB*RueM)Bx-losNhTT62FD25Cs!MXo)kYAs1j`WXu1YjT%i0a zjzj%*HCwEI$<%+|^B{l4Rm-JF|#~zq!Bn3->%;FwsbrQ}0S% z5Yl&Rq7ljF%urNl98t@g{+MX+WI|b$7JVatyvfizk!D;HtT^E|tAg>MsyO4@hshjw&^zI}*JeQb^ z_E*Rz5_zC^zj22Uwy_wP17=AEyXfnwtd(eSCo`Penr_B1wa641behJE{T{O})f?%T zC@nA{4oo}9{KZl2@qUNDBs6_}ZZ}f_4l7{<1@JB;TjFvi0AM-j|^N0ey;0^ORzu z{&4YnYg*lL!P_t8^)+442a*>CX)=#1?!Q+yA$g&7Kp^|#HX}!wPZe4V;#YOs`KtSQ zWI`Sn_sc3?XIPN}Ly9_ad7O*GlPv?0hTUNW2AU5j=@;;tv4*W@4Wg@lPH|zl!#uu> zY@ZsD$eO@66Xw~rj@1`3&&P;NIAoqbH82~`E^hq5FZ�iDDE{enP426j%0g(k3n( zNA=H6S1vaAwh@uXCwej*SeBI%T_UFZrlVVo2>1iNmV)Sandi64-jO`l+7=4& zdq7!Mfq@J7lo=-*`P$t{O>D8dLob&V{?toLp)OdMXjp7HD;)DePGqo4cyJK)HG+=^ zQ*o%p>F=RNcAz!Q&BJQCt3~tdtb2mq#_|4W_I>*~uYfXZ$LR@rQB&Q%d_MbOy_bji zd^bme{K|}-By*h98XRwPtI$6?WY$9(Qk*TS&^N+)+4R^{=(T4HMKoTg5$zalsfd0f?P*WIp2xUdpu20;9D)n~=cf1O|UtKC-$b>K*mSbq2N zt-*{+6X6IA>ATdss@xUH30(G$pk0^En1K;4xU{O;1D~fiM`&zMRLD!%IOTW>v%?r} z=X241px(S)S@dp8;{Pe*+T)pC|M+n`MmlICeulXXk;{CIkuqmB{S+qZM=~OpXv+1d zSZ+16xu!L@&dL!Lig7F7Ddq>n}tnmXY0;>Hge=w_VYOZHBNK$sjx4R&?HR8%rR6cy^;PB`v4u`$X@S_d$cOBM^ zkB}<6I#pt{Ip^kR@Oj~51>>=Z*Jj5x8!f4*+h)x_`r@;y*mx<__QiR@Ow|1m(vd}> z)^&8Lx}wtjI|x=N*|lvO6V@pjX*G6Up&}(7RF?JWb8zszhDodLXC0HK*RYRGcJ7UG zUDhhL|6bs>j$Pmi8g{2dW3g`yz5rnz!7}iZ?FLmr`GOWq;Dx+%P+bS_4=1GTgHwTkNK~m8|f<{_$^3`*pnJ8oHM04>jU_dB-IJ z{NF2u36R;eJ@&1hgjoDLR~Q0BEIM0_2s$A5k{83w*;ffoj525c;ZW#&E4TP7rg+X0 zHLH0i<)4MHfl`4YxiH$9n-u|Fj^*7SDotUjgy;J)MIEcAm7Yha&bNp6tz2p?ejjGY zq@#oFvb~vU&HI{A7h@wWFJAoszfiNHluQUYZDBPoH}h@`_({-c7sf7nk;I%2^QlEko}+g=QilTbDt zA&CwO%d597>j=@HpRNJ7x0NYI z(F0?R^~ZRa%8I~{_INv&MLCL|Q+deTq0W?K1@dKDU{Ahb)+{i%eX=je5nM%5E0NN7 zNCjRpF$cRo&Hy;E1BYgXlQwAPZ6-#)JjH~Tf${=+o&uRHmo#BLO?N`1IK`7C^Js`D zZZN_LO1L973&;;n1(&4U=pz)ssh@j->>xH1uwto7QF-zVQ29)rhX(ucUps=$5Xm6` z+GC0gAK-M|PmwyRa?-T}od_Q|$DYx3UZWPlmi{>gei*3)#QWzUX|}12L*mL~tiXgf zWEmRFX(0u^$QtJZ8Ubl*b*&g8<5WnTtvh|4`miZ7M=lR2$`fu*Q3=Q(9!URLi!R*n zn0n@oFpNj0@*v8tM}iLMPlM*>ak)8RNlFESeK&Rkt~MhqX-aeXPzgrL9;K-->rAn~ z4DYfORdkP_*q7($hbP^R(L9K+hwdq@;6p_m-dbZxjdse`*+jP9@q-pT&|Gz9Za?7x z#@(z#-j$8QHfNsaflm9Eg=3ByG(}UqbuL?BZQcHc;#leq2kfO~>enzf(^*)l_uUObVC20H%&`~ZP(_iGTW7ID%%gNfx`1WMj(;!qZ zv!zsj>Vu*JQ5fEptLUdo$#B&g47@$*!2f|G2Ec@ zqsKT~g?9!zmTF_J2%e7M)18 z?&&Fed-V!_%u@J+&lBDCi5EM=I~lpA{-1wQS5>H+vY|%0Pg`>meHCe=QqWbV=Y*Q! z&}UY%12>Vb>>XVv*_pK+7bQAey<>DmR@-;3aS3h2+&}h74VqWRDQK=l^*V48!yKXn z%0l5jj_4KOvwlq% zcPn6r@Tz+2^QEC8eD}tD=#l!51_Rfu3^JRez}stV{LG*I+*8}v@6zrvBK}GzykG=1 z+=@*q8FjoCXO1M3?mISF{@P*z(Jh5{!n7A;ayGDE%-Fw6Y9T}N|-^eYPF26BH zT?ZenZH|9t%AI14T}v2dvQ~p1pML))7{$Clp)vkPOe21A@@&_X6|ZgJB6krnNVg>+dhUb6g@<*Xj}=bH1^V3S5|> zs(Kr)rQY;5v=9NyaHpTg%$^w+w@rY@*sEW=C*OWnjNcvqt34(gc!812>l}Z4{@&v4 z2^+Ao@)#0sJV<6pE7G(|RVyIq!avK~j=O3t^&$3aF~9WG+v&&K@kgVKYu3jyltc7g zeszsdtjz=v1$Vt;<(WO-`s9^y@WsQK>`&11!P*(=ao$>UyXaYSJa!-ck*x%vCU?bC zdk6l?Z-+&&%{W|o{Aa6F@S0^K#u~fLz`kar58l_8)0*puvX^ z_IKS^+Jz-boV&uz&#D0pO=nHc+Xwn|j@J8=4H*!N^7Fg0FD+gV)UO{tdg7>ui$;^Y z(#8pE-^2+ERxjc#!*CHSd(N7(ip-ajG3O!2mdU!2ERFk&*Ud8aquVXT$MoiMkICE(gOWLd@IO=+UZ#*!UpCC7$I-BCKgglUwh2nc7z#hIp^)5Phr)1nuLlw}3hQSqi28_wBy zQ(IAELKF+hH48wCzom(}b$rf3oB1LYbYk;X5qM-Q8*VB&@gv+6|4BI~T)tGDiy=x;>u1NtKper%C}YG z2d~Rh70^H3RT0pC^W-tgk{Sr)qUD9Z_J-F9zj4K|G0U_Cu+zF){ltW=2C>ufQ5L4g z#^rG5n%auH#+goJdtlglOP{mTPAm`glVKwsyDG4A=zRM))PHi_DTCtPO!KE@Ko!hB z6)0^=A#tw4gc8;?h7;gjG{zOJnVgEFhQ(^t&E!JnnJWo_*j@yNpyKCOf>33No zBCfu7&xxlY1ouXLJKwvR@;}gh@Po-MV}F=Cc81B#OSZatVc)NDr9t`pi;1Yo@#@OS z8H9s*_(1H>WgB^T0QZ5}I1wtbDc>V%>ip>W56Q8eZ{FtjB_;kn0 zN}Nn*dtEk*|!?*MK&zq*Gl^>jh7X4JTz;#+maVV%oDqxbpqdU z3Cti6FcTYZ$v+2*PPc?=H0yU#KY#BIm+q*&+76!U^`TDtF}D(QJq`q%o%;!PUtjng zgygvY-8d88U&ncx5V+Io)2E0mAGZV7xLl8oO>fW3bO8b8OgW_ke(giSwNAvz(Xajr z2@mt2%a0dH@%t0WvSxegGGoH+eF{7b&6lUuS7&9DA5{SI#z{S&ON23^hVD1HE}x)~ z?I8iLKSnnFBbA&rRIaQ8QG(@xs*K z+p(#LqVnP4Sl*HAuhp7u?&?&#pL^|%?1F0^1I)sX8P6YRG2v|I)jqvX*J6>d`2Fss zS4u3GCcyk#qg|ibHrH)pbZ$haKk^NuPM`C=VIB)MeV(Rh0= zsPt_h8)sbOt4N}*Heqsv=j{#&?6*l|K5|fn(9H(Y&)C`e8AkEO6|Lw<4P{1$7ms;k z>qHDfq2$i1nMA{nx6wAL0dZq6-_v2>0j{;MdBiPqFoEjwG48eZ0+J}m*9JmyO^=&r z^>jkpWT>$?K%}!P;}B1(<&Z=?H^rgkZ@S_UuC5Vc*g96|ehHFjw=D+2vN(+->O3ES ziAM^=A<`1MMk+j&5 zi}+MQj{S^g!2NQaxQmtlGERr1m{o3b7b3@KI-?cQ!=n;d`TR~OIl}j-9}HF+!MudCfbdwz`pEN^*W62bKiN<@L&W;U*k%6+O zjT+FzCn_t^2q7Ed|Ma&G6NK^tF zI+~8K)i(tP;s#jOlX?oD7g$!ij})f{*yIIS`~jE}TS7`x8e(Yly@SEOv&uc5pkF)y zt3?SXqMwq}u1^<$oh{c$E${xga!JI;ot{*6I*`O}D#g z%lK;6nrg>{$Fo*;u?MlrE;L#t*+ zk>WJRvHjKWqKCQ5xCAtdI%XJKnTk3LK9-1vPO%G=wb45xQ$6j;h`_%A>V^h~2-Y^VD+D;SHwA+0qt59s*t5-8*jnKOk z6suRrlW}dc(zg!uQT-@wugo$3Oh;|nK!P+T^aGLu=r=_V3ny=sXWLhCBhxoNw1Y>& ztQ6{Q1u-t>%GskzD8KMrNil=N=bsAY=R4SAu6xb2GzWgzWV%~o*w5nc570J>Qu4V7 zSI6u17(|kjk2U1xVcr;QsL~b9UiXJ07g76MMH&~)Ep1~BxlKfev&^pW!PXTVZ;Zii z83W0Wakm})kYg~mlw=Jw-e(6J2=6V=0hqo%&ea-HY}wxj)Boxdj-b$XaDd!(V&%k! zwst`iU*B)cL?l5$7F+!4FGQ}xRwt4eQ7B{sagACTgXs@{!=TwmC>8ir)~r`;AeF1xAUQpjUyDc3cpZ02V72YZbcP0zQGIPg4HDQo zzEZV9Xd*#~a6Dp5pcRQwGrw``;0stgZ`BUK82{LUXDZJY7r|Gm@f0j*SwFb2 z^7E$v4erS=L?7(I&Ne#SGJ6cRUI4BtoF4SZnF~i0IHp1q-`FpLB*aNNg~%U*m990r9ZK$6jSkQH50n4HV_;VQaJhDuXCyf3 z|4Rb=!}Tw4|5EoqPI%e8xKG$Guu_E0KTfLinfdSke`q&H|3i{7dS*iZ)(j4Pg<+(s zJk1+VumZe_eidv3Nr{X9r+j@4`+6%5K@+FY^JV}>{Eq!QE6BfjTz4OV_a8CoBSP5Z z2Xhiw*7pf469X_R(sVJb{2*AQV)RVJu{RsG;t_nj&r)5`irZ#6$!-J{PIUhU=`IX# z=v+j7c+^$iif|h|V_!`yjKo9wchiX+cG0o_ueO%&FE_xKRv3GB;lH(UHCq&WW)iF3 z;lU4RfL)~jnlTb=q=#KiLSTFFABzjMVA74NoqVGdme0dxTjRicyF7kg8`{yc&5$xf z>_np$|E_=I{2VXAjX$0q$qJ<@elC-VsC9m8glY*k@q6s3r#?=u+o_Zn_h{OpbfX z_{J@kKhJoN#~v$Km2kE2blD2t)_T>l9r2R#$IrW(5qW}NN%kX)_e400QYkEmBL=e8 zQ%~S>SHCVtms$k#y93WRfaxIbZ#(z@%oCGbpPvDjH|BVibaskDaK|SOL>5oBwQ5#_ zug`H3-z}wJR|Yq$zt?z>yuLC`=)OK(8KLW0opVEVUG3+01u=I9uCZs4!UhJZTBk>uG~YNsr*h!H&nXvEe?8=5bLLF8AB0PB6PU#? zq<&RdbSb78)F*izEw=E~wIbUmTCB{bCB4irFy0ou+ZD-&$NRZu?$W7_)5FzbTZO1U zYuYVsbW^8^vu7bVBVV@CkqM_s`@>ojeEIe0{%Pn;+;LVL8&(OJ*))6SG)x<<^#;1S zms0-iVSmEz;*T9oed`2Y-VBO#Z>FhdffgrQvi`*xr}cMUd*>=Ln>0vQk3SQ z&!N1LLMuS&G|q7@JRdxr(vRtn$okA*ldlW^H91U@QTI$d{0JI7-{1Zw#VS}sSmH+< z{)p@>KC!D=1y9Br;3*vQJ4|ezuvs;8m1&86*l6&S+p5JTFie`T+2ZB^(|n$8NlvW{ z3KzO}r{=*tBxr)VnQJzq=d@s)Y)QmhwpK=)U27f^*A%*kGn&uatq`x3S(tXS4D7KY zYUX11u?)OLN7&3|Gg=$mtV{o9PvlXOaIFlF1Zeh66mJQjK6RT9QD*eD3?xjHZuZ0f zLF%e?$6$IHSzq;@#+DIK zG+R7SaGwv6Pci5{CDD1z(Cn}@!qBYsai?RAYDx9YEC@op3I?3a+d)16h?1!5BFO#a$2D7b&bRZhaNL_2cc{RoHzn z8G_AZ_^rsj;3rm`y{CEUWD>3ihW=xDrmZbG#Dy=Sre^GM1bnI@3 zhtk)tQ7Z_M&1aOorE&3SHQ#Pw%@2$?d0@m_?!XZ6pxnVnZZ&pWP+uPXqSHYZIE&nt zq0fliW?Gt8(e7MJ=4vPmhY;;#gn`lEw@th zVg|+;T9=*;`trFI9iDPQ5@sV5eL3yi?eT%+r%xDXGQ`f}Tp8;)bDn5-`tjw1(3W>M z#0+!~d@(5Z($zXOWyKQ$*JN|RrL=3EW~0juJSZInOXwa@!A&J}Em)|I&UnjdUn}pM zYIk~-ZhR~Z_%U`ke^9xujd$|>Ts?R#@~5A8nYdbaL+}ZG(Hfn<4_JKDx=n7$9M? z;L0tww3IP?7*8SP6?+HU{@&B{1Di2zb&(!#Oh*2NKO#W4f@h||Zm75Mav)Aqwtdh_UrT=#!oWBEK)R*JEw)kuz8b{ThqkH5c`?h3PkVZz;<)WkQRl{>5=J z&b4Hc(Sd`GD2Tgi(NwkX`gbVO=9iyn-(OYWBP=;TJy`?HgQQmuBRs53{mbRQzYM?8 ztrNPfDB#$F$g?PuX)0HJ>_a|l(G`f>{{Ai2z+<)JH*);MDr>nAWV6h#ItZtxj#lY} z0j=jCso#E{RcSDINJnK+Mu*}f*WKJXD|uW;yvbgq)xFcg#B_uR@FbM_p8cAa(36lr zj7T$iy$I-e9@QcLd4|8fJo1(#ZP-YPp84qHs6l)A=Nn2=(Td!blq%fFV2?+n>uBPT z$>k{~lw0_CR^S)WS@|3V$TObe4~Z`y5>BSi1iW1I;h3Q3FfeLQ=zALOtXjbmM!IN- z6EE8^6z4e1f~YA>%R6>qiBeoec!sbu5-PA{)>+_rT50o>?K;B5R0cc0hytlFj-@I$ zS>sgp$a!JyNUOW_aNiOTzC!M>Tm;~9;(B;XT z$4SzY4jMzjQsgLVAMBfD^|u>f$Bj_f2X2RvVM6LRcHELLB9$DNP#$VdP*y)3N^6Dv zXgC>0^UXVIIM;#z*cshmLOG=MuOjtv;9d8iAxzwM6ii?Y}p@;43Xq>HBkAuN8P+w^d9S9V{;{prtBKvH{U?d{(m z%_;B7A4D#*9z%wh3S&X26;Vu=!s?4;X$?`Aq2T9n`<>({3ypkfR=>0<@uYV4ve#rN zo9+fSKf`<83j4O*zwe;lYk7g_P*O|+>=qvCJzX%1pVChztu%Q87_)-*;mwdZ{{S4hUmK11jCXpoVmF0&?__g_D-%*Bs zh>2`5beZ~dV&rmz#+k$UGNteSQ0`EIG<ltLV%x%q_*rOJsC$C#sIx_XiV zk4r^ES^c1LcqrAxx0XKnftH2f(=?H|EBeR2cKi*+ed`3k-8?xE>JYLMx$LK#G@K+& zmhcE(0|9K!3)@`}qYYIf_L-l&ahWhHH)q|=txCi5yyLyK|DeKsVx`yIeZq`tg7U~< zYrpWiUfz!!Y8bw+qm>bqP1Y9p;f2JOxg)f_O~G%xe3{N$5X=2~rJx5A*5la=5)S<`R>DT4D57xZ6bMGK`1UZeoQC>Gpo2vqKWo z;4!fx$frM5(7p^i-=OyvG*CV0ls~ghxHV7KHo_HtWJc2#D|KreMDx3d^QOx^67rr4 zj)8igRSFf$`(-5Wel2KM>zSOn*Cs8wP47)>1W2?6MC@OD>3X2NAE=;QPXB5Rk3IYC zl5*Uj)y2S@-(*c5m68T{dkCxXDx}#o4`Y37TpX#m^6m?#3rOGJPfc~?>Vf{sX1^Bi z9;mxy^M<~*HW7<-1{-{7*`IuBNqi@eX8BK*{S*X_MBRa%3hiT;17xqo38A^+iVs2~@bT zt1K^iF1Rv7sd5u58shuZ+6?(N`J6RxdwMB*o=8@XyxrKI+=8i($L-kH3s7yv0SiuW z?%4i5Mp3CWsf8l{z#PuTBnMhQtxr4;p&nU`ZAHcFty0hLnrYuPTn^ileq!3LZ7N%> z;uHf2za8h&^_uos`oM~fhG}~G(OaLomM@39%d#BjB0m~`oNfpU^0aU%)EO6^*f5{k z=lH93#FnDFe8gv*GP!<~E{eS_&S-#AU^EYs~<7skXRr$EZy)93Q)})!EwpV@cZt1HcslUV>r!h!v&y|_Z)F9|c zmM+RFKd)Mt8gbB0YL~DvS?;&#^hHD_1+%_eJ zq%++>8GMivJ?0)_>1!=3B%Q3`hwduiF97AI5nD(Zp$E1|nflHE(uRiIUV8gABeT2V zyP2WmHy4td6i#0&=@Tw6DJT;zP(K7G2M*5xj^}cE-81S{7!wU((>H}iTl<UvdsV8MdpGGU~Gu5kNKK)FRd4eNckP$T)5@@ex zo2X1$X-MG*1;7$zo8l>fRGYMIB44p%(u(55GE2Ln0z>$2!byj`?;zVtrg}lgqWz@Y z^h*$x<-|H__vSFUwC&mMLR5Ax>)(BQE>bX%Fx2FGiM}pj{5gt1Tg7VD)7LqnWOW`C z2fWr9D@0P(g0VD43=Gy$ODkEaG723v(JOU^%NG>DMSogzQu_^bL6*;EI~g;*)0FtWP#xJHf0c_!~ei9P`~ucHlm>!|fXia;n= z1R+@+@Dc}JcW|Aw53CvOQxt(HgbvLJ1<0>-jbwEsrw^-FBj?SOLJAY+WxNj!`UHl^ zj0`@RKfV#F*}o&bq*C94fV%f#0ImDcN@SMX55-7KJ_am-e@IV~=Lfe0@xvnF|&Wd=+P_|o`KB-6FQv36P znpL(A#q6a&O|y2A+V3Gysbcu^a4oJ!U$SvJcAI#~BTx2Xzd`X;)JM7QneMfq>y3tP zNM7(s9ptcA6s0Ja_ikCG%3M_49mvmzJcylxDGbQq|TNwIjpa+ ziz)Al%ZMhIg%3jh{urSi`L?n?=417$9o-d8QmViDXym|{d4I>n;P;a+n5-IZ{QAF& z_5%f)!>Kn{p7|o&N)h{l)K@{HDGri|65!9b^XwK6 z3-H`tA1`J;9i&9_20xZ05#LpIOmbq^|un{jaW1|=(@f<&xY4B=KZZ;89$n9 zAqWi-7+eD`_P@u;AN;Jg_T#U({yAdtQGq$Q>e7z-uyom4D%g9}J?a#4dQVoYeQBao zBX#iP$M7}Y@ZzLT>gY2)JP~2hhKCwOWb#5%QX-NpLjU|&+*yeKJ%I`d2}A!&1{D+k zZ<(Z2>TP2Qd8n|Yr0D;v6@iNW7mA3mu-Jdufx^_+#)_m+5m8}L;)f5FwN!}z2S1`! Ac>n+a diff --git a/publication_figures/src/makefig2.R b/publication_figures/src/makefig2.R old mode 100755 new mode 100644 index 1b99c30..307965a --- a/publication_figures/src/makefig2.R +++ b/publication_figures/src/makefig2.R @@ -118,7 +118,7 @@ fig3 <- plot_grid(plot.blood + xlim(0.3,1.0) + ylim(0, 1) + theme(legend.positio # plot.subcutaneous + xlim(0.3,1.0) + ylim(0, 1) + theme(legend.position="none"), plot.thyroid + xlim(0.3,1.0) + ylim(0, 1) + theme(legend.position="none"), align = 'vh', - labels = c("a", "b", "c"), + labels = c("d", "e", "f"), hjust = -1, nrow = 1 ) @@ -131,7 +131,103 @@ fig2b <- plot_grid( fig3, legend, rel_widths = c(3, .4)) fig2 <- plot_grid(fig2a, fig2b, nrow = 2) -pdf("/work-zfs/abattle4/parsana/networks_correction/publication_figures/fig2.pdf", height = 4.5, width = 7.2) -print(fig2) -dev.off() +#pdf("/work-zfs/abattle4/parsana/networks_correction/publication_figures/fig2.pdf", height = 4.5, width = 7.2) +#print(fig2) +#dev.off() + +# source("/work-zfs/abattle4/parsana/networks_correction/src/config.R") +library(cowplot) + +theme_set(theme_cowplot(font_size=9)) # reduce default font size +## categories to plot +#cat.plot <- c("uncorrected", "RIN", "PC") + +# select category to plot +cutheights <- seq(0.9,1.0, length.out = 50) +lambda <- seq(0.3,1.0, length.out = 50) + +cat.plot <- c("uncorrected", "RIN", "multi-covariate", "PC") +# select category to plot +select_category <- function(category_name, pr_table){ + pr_table <- pr_table[which(pr_table$type %in% category_name),] + pr_table$type <- factor(pr_table$type, levels = category_name) + pr_table$cutheights <- rep(cutheights,length(category_name)) + pr_table +} + +plot.thyroid <- readRDS("/work-zfs/abattle4/parsana/networks_correction/results/PR/pr_density_wgcna-signed_canonical_thyroid.Rds") +plot.thyroid <- select_category(cat.plot, plot.thyroid) +plot.thyroid <- ggplot(plot.thyroid, aes(x = cutheights, y = density, colour = type)) + geom_point(size = 0.3) + + xlab("cut-heights") + ylab("# of edges")+ggtitle("Thyroid") + +plot.lung <- readRDS("/work-zfs/abattle4/parsana/networks_correction/results/PR/pr_density_wgcna-signed_canonical_lung.Rds") +plot.lung <- select_category(cat.plot, plot.lung) +plot.lung <- ggplot(plot.lung, aes(x = cutheights, y = density, colour = type)) + geom_point(size = 0.3) + + xlab("cut-heights") + ylab("# of edges")+ggtitle("Lung") + +plot.blood <- readRDS("/work-zfs/abattle4/parsana/networks_correction/results/PR/pr_density_wgcna-signed_canonical_blood.Rds") +plot.blood <- select_category(cat.plot, plot.blood) +plot.blood <- ggplot(plot.blood, aes(x = cutheights, y = density, colour = type)) + geom_point(size = 0.3) + + xlab("cut-heights") + ylab("# of edges")+ggtitle("Whole Blood") + +fig4.wgcna <- plot_grid(plot.blood + xlim(1.0,0.9) + theme(legend.position="none"), + plot.lung + xlim(1.0,0.9) + theme(legend.position="none"), + plot.thyroid + xlim(1.0,0.9) + theme(legend.position="none"), + align = 'vh', + labels = c("g", "h", "i"), + hjust = -1, + nrow = 1 +) +legend.wgcna<- get_legend(plot.lung + + theme(legend.key = element_rect(color = "black", linetype = "solid", size = 0.5), + legend.key.size = unit(0.3, "cm"), legend.key.height=unit(1,"line")) + + guides(colour = guide_legend(override.aes = list(size= 1)), fill=guide_legend("type"))) + +fig4.a <- plot_grid( fig4.wgcna, legend.wgcna, rel_widths = c(3, .4)) + +## Glasso +select_category <- function(category_name, pr_table){ + pr_table <- pr_table[which(pr_table$type %in% category_name),] + pr_table$type <- factor(pr_table$type, levels = category_name) + pr_table$lambda <- rep(lambda,length(category_name)) + pr_table +} + + +plot.thyroid <- readRDS("/work-zfs/abattle4/parsana/networks_correction/results/PR/pr_density_glasso_canonical_thyroid.Rds") +plot.thyroid <- select_category(cat.plot, plot.thyroid) +plot.thyroid <- ggplot(plot.thyroid, aes(x = lambda, y = density, colour = type)) + geom_point(size = 0.3) + + xlab("lambda") + ylab("# of edges")+ggtitle("Thyroid") + +plot.lung <- readRDS("/work-zfs/abattle4/parsana/networks_correction/results/PR/pr_density_glasso_canonical_lung.Rds") +plot.lung <- select_category(cat.plot, plot.lung) +plot.lung <- ggplot(plot.lung, aes(x = lambda, y = density, colour = type)) + geom_point(size = 0.3) + + xlab("lambda") + ylab("# of edges")+ggtitle("Lung") + +plot.blood <- readRDS("/work-zfs/abattle4/parsana/networks_correction/results/PR/pr_density_glasso_canonical_blood.Rds") +plot.blood <- select_category(cat.plot, plot.blood) +plot.blood <- ggplot(plot.blood, aes(x = lambda, y = density, colour = type)) + geom_point(size = 0.3) + + xlab("lambda") + ylab("# of edges")+ggtitle("Whole Blood") + + +fig4.glasso <- plot_grid(plot.blood + xlim(0.3,1.0) + theme(legend.position="none"), + plot.lung + xlim(0.3,1.0) + theme(legend.position="none"), + plot.thyroid + xlim(0.3,1.0) + theme(legend.position="none"), + align = 'vh', + labels = c("j", "k", "l"), + hjust = -1, + nrow = 1 +) +legend.glasso <- get_legend(plot.lung + + theme(legend.key = element_rect(color = "black", linetype = "solid", size = 0.5), + legend.key.size = unit(0.3, "cm"), legend.key.height=unit(1,"line")) + + guides(colour = guide_legend(override.aes = list(size= 1)), fill=guide_legend("type"))) + +fig4.b <- plot_grid( fig4.glasso, legend.glasso, rel_widths = c(3, .4)) + +fig4 <- plot_grid(fig2a, fig2b, fig4.a, fig4.b, nrow = 4) + +pdf("/work-zfs/abattle4/parsana/networks_correction/publication_figures/fig2.pdf", height = 9, width = 7.2) +print(fig4) +dev.off()