From a95448111f2bda9ccbd219f5126e6d0671c5cf04 Mon Sep 17 00:00:00 2001 From: Phillip Alday Date: Fri, 27 Sep 2019 18:29:08 +0200 Subject: [PATCH] add test for non-zero step size correction in glmm boundary fits --- src/MixedModels.jl | 2 ++ test/lexdec.rds | Bin 0 -> 64899 bytes test/pirls.jl | 15 +++++++++++++++ 3 files changed, 17 insertions(+) create mode 100644 test/lexdec.rds diff --git a/src/MixedModels.jl b/src/MixedModels.jl index 3d2ab0494..c1cf6e242 100644 --- a/src/MixedModels.jl +++ b/src/MixedModels.jl @@ -37,6 +37,7 @@ export @formula, Gamma, GeneralizedLinearMixedModel, HelmertCoding, + IdentityLink, InverseGaussian, InverseLink, LinearMixedModel, @@ -49,6 +50,7 @@ export @formula, RaggedArray, RandomEffectsTerm, ReMat, + SqrtLink, UniformBlockDiagonal, VarCorr, aic, diff --git a/test/lexdec.rds b/test/lexdec.rds new file mode 100644 index 0000000000000000000000000000000000000000..cff9b2602993178d9f4046da3a5c380d4e23b4fc GIT binary patch literal 64899 zcmX_m2Ut_d7w@{R;);N-OD|b<6_q9`9g?-OYXMzUdSXC8MFJQ@3rSR@2mv4tCoMMfrK(p$$J~`bR&!dp~mn{Pf8?m;PuzI$Kneq%N5o8LW-@9xQD@ZSP+Ql)`IdYv@$xadJ}tDwskC2)I$JB0f7 zsR@VHJu?v4NrS@L_qfsT%vwQVZ#E9Q%WrZHIjvL`;%^ZD@_{~ zE@!@2CO9z2v7qS0Rd^_aoKvhuDO{sd>x+}HNn$h+ikBSbB}>~cKp=~9Lb9p}9n2U` zi9pq}8_QGcQ)zfC%J;kzbMDR_ghKL%srqK*fV8RRy(KLJ>;wBy*8%`(sR*d2Wp^#LB z(_+tUmKFa&J(7lzWDujZ;ii?3GS+ZRTyt_G2RiY2r0$m)(EhS$^?7|Abm`-3B#)Rdeo-hFbWvBg+;q5qFH2Sn3iIJhRnn(MY z59a_VBj&zt7WJbB3>Z~WSXEmqQ{{u}5PrzX4>oUsBEJIDs)+XV+lFbxvD{mjp>kY0 zQ39%My{oH8CFaJxr_6{qK$$hV&XLFSqnk!8v~qO?G?xXL`sEEmoUs{8j9JTUnmPXf`2-)gsi;Zjxp! zY?lQ}pGZlChhwHB^~iiJRoYP{Udun7ENp6pGSlNb8&=l#L74(vF=Ew@l1gNoHH>4e zW}Aj}quS!omCMR@vKM7~FjU2hgpemQwB&#VDLYq@(0W#Mxd$!RE)?ztLp$v0p;_3n`n@AMMDBBk3=8rllaH?9y$TzvkBt2&rv1}NfJx=gUf{8D4F5;|z zVU71rG0Q-;P(L?RnK-)G)wWk}#V?|#GfyJfhO=E$s8~DIhnN)sjCCWBnc_0OX}Y%kSf90K6JvFf=EU8 z!(zOvWrvM~1OhbN2`-={l!n}G8m;q4gHUg8pY3wF(P(dXZVS8|QD5HHakgo+4CBEs zbkauGF7=TBh*nSNBCoP6**!cxqfLys^W07P%B?V-(xqac-gh>pnxg02;kcPZ1I=hC zCP@u(%DxSV#mUHv<2I*NRuzyX4hUBx$X6E)4@HG94!&wd6e3GHzo*+tUVo>G4Iq^Is1+EiBS2<6A>};5Io+C!E|tuNW!C#gh8@_%@Vc&8jDf{^C&Bn@zqt zaBsTiJKbXcoJ>hTjw1W?v3g1utpPJkTX1e#y&hP;JVzfhQ7JE*Yf0*7#}CEN?Lu&9 z-Hsg|tIh97UO|UC7R4xyCrJKVMPT?sk{-tW_XUTlwp^Y&xsV6QEmswR1OQW9%%AB$ zLkW8bZcx#wwDkHZj_lU0nkHs22k|u<82DSqSOa%cnK_NfcaNWRG1UwHj;(_6O02JB zE@4Dt7=cJChQhi^I_>=7nZj}5cC?Ni0G#|hGFs}`VNRQ02!s@EXW86gea{0^4)m&W zZN$chLoVHc^qpE#4QT%)J6pg{Qfd(};6csvtX$fy$z#yPsjTAou_954^XQt$De47h zhYBWhi68x?Gre})zC*P-ILV~aV%Ag+WoANUD<%#|5nvssn#=a7eqCWJSr|cmBoF1} ze7HlEApBl4PWw8wsH{~z#@w#Nxlm?TsAg0$d5EofbwXUc3Y)f@DjCx8YaG*&wSHF9swnVn<+bj~GgYv1CxB+=386?*-zgYenpK5AsE@R4kS+b<@Ytd=A0BN!9JXrQ zmR_I3zC(_r{Ph@s!iZ-oDPV(HSUr4tB2=)(Ll!;MQ(u!Af_LDCy(7W?G-NRwfp=3f z#f!?ow?iNFq7?7>`5h-)J(g$Vm7kRi6^|m;vn5S{c4HTP`a;H( zRz2!P+T1(who4R3$x8>?vXiphpO{j=?Ngj{Zcp?fEh4;i%vO>vVgBaORsV zwK_E=CqPgH81Md(+^>9! zn86B=j8LcA+hYTmT1FV(+~1FV$~!oj-(ub*Z&~*Ox&%`n#PNYlZFv%3eQ7o_xP6S{ zFPY*skx8eMV|$)c60q-0?N{fIS>=*V(MnI2G(Kf*lgv4X@Y9_iGkjiQ*J<#xSf*=u#J7Qw^{Atjb0V;N{=H z_m!c|TKeXF5>BJe0GDw!$-#7$0mYiM_~>2EHBf)yv-Fi}>=A$)6jvT`niit=3 zs%XR7#nZk&-K&>5ziE{)qmg&Uq0qb=vUk8YSYlEIk6GjBU&`q0VZT2+Vu6}JK;)$n z38Su1E(IWRwRcKfJ>HNbgI_WnuU+tUibpYp0b|3_V?#5<_NHOO9P(=G6fRzK2Uq4@ zwAxy^v5Rla3iePuhC6b`r-p?b4luUmt>H<2a8!dYM>WCS@AWJ{X*@=_+m}v$#G_zA zFbBRmwM!Qsq2g0~FjChAnIB<{g-S~D8zG1HccE&x@x9blY#TQSPg!R79$9(e*-=U_ z{~H<4do-GO(oSNk8Cx}uZMJl`E7Wt%J@0K{ms~(@mtSUA4_sEmpE=zws|7BML4aBF z(0Ql!h#b1?UyT`}vPECB3->pmA4cczijM-%pdPbJlZ!c(^iqJ({mCAlV-@e)3E@y_+e zVa6~|Oy2DrZYfBo=G@hLYq`1~*dy$B?q63cl$$XhGZI+Rq8chB0Y&%i&jxOUHQM$Bm$!%+=fMQyjXe9YDZ+73P ztfhpLY*d-+4&}okKlo0DDNg{@E`{ZZsBctZ0WU`&`C8M1&Zh{7H?%Zn;((=S3793L z&Iy|fQB(Aep#uJ3P!Gle>=DEKJGZDh;dQCLq8Kefd7k1or+?lw>Su>;i7N&k1D2@%q_2@ktdLezf6s%25AGb{ z6|KeggGQ4vICHz?UGwP$nnAOz^xc91n|e(tyTiq4LG^`9AkvZp1{G8ry#Mk2^tK_A$4J~IKTc)*y+P4MdKl^{yXQnDOe;dD z8g_!vmR@K(ibV;iD7czDb0sLQoetZ{63}>{(}M9KjaNguo^b|;ap(Jeoo33!$Jb_++d`Ht!sU02}&RaPy1KMOAGz$Q2EDLgH79wBuF zD~n#Su%7PE`~Fm&DPDlL;C$<)6r{#{CC_EEEQcnjUCFAv;P6Gpb;dUgcII=WXQuZY z!TQgrbwFA2&aV;f6k&tpl+bF{=89ywezd(UPL~P(J%F(Cu#QEPbH}wg$wyw~&AV(w z?Gd@7$lM9p<0M0sJ?1Tb(yTLiK(U9x;N6_dh8r(2T(*PSgwM0Mp3-4|!Ke^5){EcB zM5q!rF@0Fw3y-)En_zc0i$meLZ&!#{y9A^ zX;GLh>=|V!Ytd_Eodl`y4+#P%IplP$*b!4sx%#ynTR7H{*i>JMJF-0=A&@JUODQK4 zQHo@7!wl=jryy z4y#YoQ$daS#u4U6-)JVGDjL0_TJbEkPO!qG>%40vP-- zvkbckjL6Xu0D<32g_f8F3;`ssjmAn>XAVie9#0osP~Dhzn~s`Ypz2*kB#hmS_B`jp zvu18YM6MAIC?8CLDXaDpnCQX;|9E^e1y*@OkEb5BTh&K;jh%M0sb4+8J}rBe+~%@J z=#%`!C0%4ceWaL<#9|h)1W=JWp`YZ3yAuB$ucv{__LR78OwOh+a5o7LcO@1?AM!nf9SH7$O(B@lqw~C;!s51m&hm9Zl;zqsK1UHF52cS9 zcW(Hyv*>I9u#&^eE_3>R=Ls)G9*DJZbOeL8LuQ*BZUrOS=#m(oh3_f&PhF+UySRa4 zXg)oOQ$l=9_~MeLAjks8F8d;NHowPz23%vrP2`R1CqBXts z&5zymSeTyPSH3q0CBAF!l{p{WXD50_=d_HbA^uL^>AAXrC!FPBeWT0~JJr;0kciP) z>WbuEw+=r{HxjNV<|RH>=fEQqV7z%C{hoJ!;aodke)NuBC3}Ty z_u(AzEcKc;%c{>Bc-o+j?K)@!2((%dHS_-4i}}3*cj)8n0)nQ?XVZy(6R8aW&v>K# zrXura$f0&Z2@x^k6Y_nDK+AJ7aLGEF?%g)w8?8F{(7*~ME=AiIxi zW1$X(HH`Y^rZu=?dWBr(z69`l45be_g|%D0)@YvZ8YPXv3Hu4%u@e28v`d47zVJ22~wxWXGm5SX~iF@-8EPp~b zhP=`j2P}+c<*0aZP*wAoSKZXwGpsj1UE&_wrg*j$i3qx#K9Xc%!~!@k-!lb$&_OoP z8Ar-dLk<-n@$*_YZyJVsBikVck3MDj@gI(4>c{G9BA**B0H^Rh{ z#QEA$zGzQ=z`hDM!GJ^b$kqBev!>a~j>LHd5OE$aF@>r=kUEX*Yft!LzJy0(5dKhe zE-9WyTx6nY{t(1S1vFkPa`}T0#_7G z#yj95Qccqe7KPtI-$eY7Sfm=TQ=EKti2mEc0HLL5<_c~?z6@Mcg*|K6d%ggF3l;>& zwk!Bcd{rT#ecY-7G~1YJAV8jdB2iZ>#0oGfhsos!&TbBNoKsYxlbOeG5K!0Dtg}7j z1gF1p?IMfQPjYrQBM$$;(~52KR5tt|bZcu+%#CbbApUh!Q0;t}xd+ql6gZ(|(uz5|`h24r}Ra=F*Qcx+GtWJ~l=y#|Tlm;XOrZ$$!xP zkwkwT=?5Vqv1p}Bpl0s}pMv&7^sDIOswmXIIOd4UJ;||(Jn1Udg2(9anDcdRgQ`Q7 z61dN}_KY%9c%por^Fk5-;tFfP?QnJc?5);74DN-{qRZxsQ_-G5vyX8XB$jdYhaBy$ zTKoxWaewH_k(=TO|3DILd5?<(t6&3!lW^PW0T5HBC&<6A1=kz*`6w&kLFR47eQT2X zg-!yy3wrjX2-sg9;fe~k>*-qfhEU=D@&3ZKfsXT0Xddqo%Z|;&tM5#OZdG;3vvi#o z30evTf*S1uTFAVDqtc&o=rA_a7e;1&fR6PL2zMA($Wtf_)$;+(^hOY2MrAG!D%6rH zQp!IyEav@%lp$TFpu4_!lnCIt^EIdQ zCCTXX4T~y6;b|njMNyca8Re6Rnn2?oSNSi&hpPA6z0pB-JhB$uDE zrQgoQoJ6IFrkeZshH|6wjmicOVHlgIZ5OoG)KrabpSAATt1|o2b)tc3gR|T81~UYf z+l`5=gG9%^)6s@&2ZFSJ#c$GPA5SgEa093@HW%wQ8{m(nM!s)ysEg6D3TijN&unF9 zrbfPLBA(#}5MxqPBZ*DKB(BNDI@ZCUNJBhAo1K&znbJhO!5uD%iMm`@Ya29TfL}?C zq%;vx+~HR-s`oLAf7iv?1o>=b!%`y&O~iQaK|D&+g!qH&S`}j+RHthfa zDu1vEk<4{{A5;8yotITm#IN`x+HC(+KD7xE%O$^#2@a^k?+@}dz?W}j52W(RO$Zd1 z{60oSj2Zj84rUb;|0~{2hZ4k{H)TXUZX&*eI{ub6a-s3oCWdvipH429c$ynf6O(Nl z6uONRn;OYzB3|SUzmAE5Nu!JzF@;SIjmw38b^G-3#ar23X;ar2u{XH~)qKMD)X3Z> zQAFd;ZCQiG7wWtoapRM?^0V6*HQMa@RQ~-Y#4Bi`KbJfjL*17)m7NCNK?fgu4nH-T z85B5J6==O54!g^BTu+f>Y(fWZqNVzYzSBWQCzDP6>rQuWVi51ADNQ5%VBvlS_M?s& zQs3tDGOci=X}H1|oW5y?B$Mj4Z$DgcvT<9*dXlqD-k;&@*Q}$~kEJvoNU-FVT! zzR=`}gPTp7d7Zuk@=vbo(F|$a9j>cw#!=+$?JM{U$VJMjuHfyA#a<&Advh38Ki2A| zIGlmcFeO)JLk%Pp^sy7;T=JtcTFL$57J4tG@lQJ(@{* zDsrwi1A38w+#J?mP#=Hcf|$1Z+47;T#@*Q1PVCMnviMg}Ee>91H%-4fC7Pu5m@}~X zdyxyspx!F3RM0hrsDkFL+}l?&U({vJL*SbNnOj zolNK^W~N_@C}7&Wem8y~1b+c%bz4wtgpTcOa=T;FhW|p@l{RwNo|>8LlWTSRmJQ>5 z#>I=c>4%-$0JjfP56=FS;NY;If26CtLFdTr;}^@dHa))5s5?68BJO|LVc+?)bY$lK z5HZcaDE_*EdMtc!VDne`Pr3#&-|5kb-zH7bzSG7g*^1}0_jiXcm%nR?`iWFx@l|X3 z0<|*39Pt%sKdQUUWPT^X?O+CV>48yPLXbPDD?8=tR=xLucx{ta{JWtL)Zf!V)YBU% zB|Lrj$%~Bs1S#0wyBU6x>UDs%|3&Cm!t~|U6G3ICwq+5wnAn8;3C@@*Lj>=j>=wM^ zPrmT5zg}RK>HII8d1#?a&xFm_S?{A3hrnjSyTH{`J3Mb&feo4Z%XVNM7t!C!2$kL| zo_(m3k|uj9&v~)rowgTY@8Z~X>gcv-TAOt06AHNg%{%Z1n3IVHs#YQr)qSX~$*n4j zZK?fJm|9;}du>nMJ;OxDE=oeqU3jdH51_5u2VlMWPx335xDQt!V%K$d6m!9f(YRZa zX#$M5b#7hn!rkuHlbhMU)+C?@Oj3t+nLfo$)uIr;dpnJ|h&4z*?Ns*P@g3)J(*t#< zv+B98Uqa12ZSFI}R%~-4e9mf60l~L2SUu}EKxc6c; zed^tK?rtrT?(Xq>mhL{C2GyyNu??+?%sN2}vL-L(-RXVt-H>KQKqo17VQI_1V9(&9 z7D`BQiE2~go}(S_{O2Or+9u+iD@s###VjPr|Ij6-Rmh*6hn*JY+Gkt;s!GW<)cT4( z7nzq3W!B~`GiZyc;=Um+kIL62tPfxaUrUG24e30dGw7HWOwGL;ajIg{y)Y2{2ET{i zTHWLNJ9wudyjlMk!zA~GQ9~ILP3mw(4}A6OXA)??Xh(je}0Ui z+`+6lwT*X$r(YdearzbM>N7*0S#sd6jD*v$KMPBn6Ac7?1sj`HUJbg`<@yASfuMa- zw50*#?%iBn>J@H#d3quh-l!|+OWZ2|J^!#3`0$@%)w0bJqUD9I$7Iy7Ff!N>=faJL zm-8nbK8nf1rRKk0V|pKD-^`7@@=p=!W>=!}3$qZ769M!>qG@LX%qcRqGE%)BBZ66v>3|3c z^}(1s*+(^UdZMa_iu=EyreEZ#So~kUo4uEnIMogD*$Zlj{6?3?qIk%l|b}4LiR@f$3@_;>S)OY|oPj)iBXPRFkbrW=RgVC^YOy zm_ivnqM+h=^cw(r{uK08#vNfdEiuZ}!+}1d?h=Wp%n>}YXf(%kL4oc13 za449nVY9~R;Y`g}sk+QwT`G`i_@cnulHBojC7*5w6S@KlMwmw(^<>2od&$n%Kl{8r zU+4@atB4}6JlSy{ugAS{gM)0I7=B)bg3Y}Nd)!JcsjqafGJS*E7nOTH5$1z+AmPpi zda}TUzvRDdwbL$d?e8r(ZbFF;b$$!&GqTNIsrr}aAg+0O_2L0eFE+n6y?*6@PTQHf zu>u3rscQzd{=;3v@X(^03w(dSYb_u?@`ii+9tmB+PSGI@cBKYathhG}Ppo4zy zHC5%8NKdamZy}wj`ZF+oZ%~(fVN{%#dfot3lfd>ncaVGHndh||E>_R7Z|vg3v!UHa z%-;VRRelY1gLy_V)LM1xR;9o5rIHnEm7@W>mF4@_LhthXX~~&MzxRfhzKPD@#s-C^ zS&45HUz2jN!PVS7y~k1joJuC8`)O#6$))R-(? zYpJxSTFL5Lox2Ce4|1;?Xn3`66Cc7u>l|8BtKDO^xV}v3k z!Q<>U;Yy^u;l+IpbzF$?u+>CW7h@cBMmQT4nr8{&d&kbOA%7iLH-uh%Xn6S}7Ghd> zX6v^^)Rll;(n#^&f4y1@3|)zwln1uXYtqYq0jms()O)I7XW$ZCb0QcW4hFKHwDQ=^ zAKzPukO|7@)b9dg9%vq`>5NoF=d@uI7Cc0@Mq44TWUD;6Yn%rnOkaQ+)1cnj&?G=` zIOV!BNFCY=Bh_rk(p)OSl*^U$dn660O(h{8bn*L30T1PrPbX$~iq_X4XOMBjWUGZT zbQ>n^QxG@cLkz+y$jSg8#J%ZT$J!sXPM|;3W>%+KLF%^Rw4dz{irV@gvybPR1lL{J zA5^)OeKU(Wv+{>vDmP0&DU<>b1_0=>D zl*bVF2W9Ew5nFxPO+@Oh(dsS1+<;FpHoxM+cvUy!Fc1X^Czqg60>ybtiucJV39wim@wit zFh1nBl%}+V>tvmBW4tY0L>w5w*q^pCc~-d%*O)KQ5U#A)}pqwVwK%aQHM@)NASzSM7_X%Yep;(A?#z#rOgJS3ALrY=E zsAwbK5qbV6m54+n)v6Fg5&r5Rf>ihE(6-IS(@;F7LponaP$`Y<9apf)@oFr>$}B=oin84BS6kK zthSzSEfNo9{xq^74MD67-OC6wW>nX=l;(0dtSSo_GO(u#7Hu?O0y{g1JsqHQHSzMW z)~(bb6mgxuEguM}r9x@1xnq4Hy^YJ@(Pit^H($!Z1|-+4ncd1oT$ z(n{3N3|cBQqc#i4kBIf_TwgX-o0^>*4H!Q3Aqz@-)bfGCB|hWQGND7QF-uXnz;^Lq z=bdQ)rEpKwUS(2kQL=mas5mCBLXE!Kt zJcRfWIyfM{7>cx>9Ak;=Kef~uX=rKeKeIgTYpwjKnbPbG=;cS=l`R!Uyf_0~h`V8p z^w-du+LrZ})Z@544AyAR@-2&#Cc$m2j}JWDBAk7P(7q&ubyQ^V5bM)FEg?O>ecK@Q zX@I&liGvNZh7P8kbsK5_kmjwO1&PItS*Jw~VK)a*0|FkMHpXiQx)=HvW!VR?CNDOJ z|CQBt(;&CS8qd%v+pS)-Nwe}BeO9BY43pf9;DbA>K2$s8endTvyh7p)m~7V8x|rSF z#`;vkb?>bLgh$TZOu&Z}TpFUe4wJ7qZe!IF#Yd^+L13B*t-;!>QYSZwfw`FV-6So_ zH0%7xZV{jOjZvu5D09YwFePeH z6_w4<;lpY*sZloCgJ|TsEm)tPd+5g=0g``uoKJwo293s4d~_g0Asf@nqv*jlG@B^; z^p{lYxE5yvm%pb^K(d$t2IUrMeP;~vtJ4thtiWt&ZAQ!epz0hH#yC*8ySiZ_-6KSS zCpugp;Z07C27ShQ3?CAN@`J{zmag*oAqqga-vzt1!|o%4XhV|plBsT}4J|4P5@|wX z*whx=pi|awYXg<_m8{mPMWd4k?XgQ% z)=_?yH;@72A^oF=Ug{c$uPj`%2e*hvdnbFxCn>cp*Ppc14#!yCe*NX?c-qWW1|h@$ zvD+z{bx3ulFt1=>P~y7zBK~GXf#{Qpa040mWjX2waz47C+8~)ufj-C=)es+jGP1Tf z=|ZVoOtfK8#+SiuVrs@BJ#R7WA1J#~-0&~Ar#A{siKK>YA04K8 zg75DK*o5+f$ExPSk!GQOq-vC-ll9kK<$ghP@OQas{>h%}WFKyjPu$iQz#tP_LN+ks zUbeh8DvwJGeyr7LZP6k%s3%r&k7qwWV%;_v!T!kXZl8{dTKm{=i*7YIoWAKIY#z!4 zJ7gqHnlK{oY#;c!-FQj+r;B4lEz71it8Yfd8ebdKtM;&tj5P_fe@R#jtuwGUwK%Er ziqD!VIo1!-ul8rhY)g43_8*V2zI_jf(aGxjGb_|IP>TqTX5{#Pn?H7`%vS%ThlTo~ zC6UDQq2$#djT}8v-Z?grd>yS_jddCK-_3sY?|4@2aS_n7=dLg5xB0xR*!f%b2X9}! z+YB@VYRC!!j( zhvg$fNt=ymswe>>t!3i3-W{k@ir-#@*+7|u!|P`2e0Yz)dq((g+JnHXmFv&~-nMr8 zTBkKg^h#u;{dJhV?c|y0i;qUOXN6u*AQ+z{PRB4lvSMyUt@t;G5dUp?`h2t`PQ$83 z`t-$p^dQx!@D?YDn^fPn){gH;>-)Qw1FfxXG$?(lLwb7u3!k-RH*DR_^=fUa^6SmZ zcC2vy>7uHCcShY2I zc8rqrN>bwX0St6|r;%dQMJB0aEHXB5r^}rkjG8Xir$^lJk&xUp9}Rvtol?GHzDsZ^Ly@4jsry%1GVD1EzWI;Oqva z1#HV4y&hhDD=PMy?Zs)Y$EH^PW~~85S9EqmET9D!EZM`nk--}IpfI1lC^NV|%ddPyX$LN74DgW#?ugyNE`h z=fc5VDN#cB!$!#uFFH5!y!2y;T;=4$G`$4!+f zoc85#;FCu+ttTzIKP?|UoVIed%MF8$dp+9`JU6Dnagj)@K{9tI41qEgp1ydZzt7D3N;+~VIyhh8y*I(nYbJ~T@ z=)dpeJGvA%k)pfh}^Mw>kDDN z*dVE5ZAGp3{e$+S<^2Qao;#vfe^;|YzJYJm&GhdsjZZ6a@UhInQMC##>#Jv>MdY0? zKVivz(+^`W_n#|1j8fjQRDD_*@zhX4D&L+Ty|F6ZS*$xi5jdlVyLuq|yo%;A&)v}x z^qD1}_+g;<)tGqjVK(ugI&Ao;WvpnauS{w_vH(O4OB*Iv-wmS(D;D>@`^04%41ZXz zUr*frBXFs{&9+ZcFne||QfX^%Ek1a?*m8+DmB2|ZTwmf0q34_L1K(c0JrQ)Y{{@PH z_bDrm`37>-qbwx%=~%@f_bteM=%>V#hN!_AMt=MtS=I){McN8LGdR-+#a9vq+%6CN zR!VIhUM#DG=uwoy=Kqe@w$~1abcZ_?f%Yt{jk-j^(8#hX#Xle?T+QXIs=m8Rmg674 z*+a8P$YkZ@5oW4K&fx9&`9RHfal%dFF6}XNyq_DN*Jk$)>7026!E0(FmRY zkUnZU3kKH06U4irBu4oabJBzslrrWB5O(~tmmvl zi2~*C3(k@AuykwcsX+(K*knXh@*jSrrC*hdx;r2~{G*J)u@9}~GyGZ^UofT*w~#$@TUw;-*n z0|Tn2YJ9j;b(bRp-;depo_BJq-&Z_N9YddE}}%MX+9dr0_#k z98m;DNm9}o;(AY*S-#uRg<6TW8v!*@b%->|Q75nUB6584DY~vp|8$Q3jEs_xw&z#Y zXS9_n9=wk%y3=Na@oIxBJls%?`V5-K+oJUtPJ3vs+KN^K73^d@LrKPzrMpH_2Ye-S zHA~o$wE7cz?-kDkQObKR& zlX~70JEfn`iZ{?EnrS_osHdI28PW+WX332g6Mx8AHUKtRO=bcgS>jFC=3I{!f3(GV z2lcqDczVv6gG)ft6QdKmnHiwwIQ~OVeP>h6%d3Ze#grl2lKATCi*p zT!%3pZ`LpBCw$&P^#jEI6fgb5-CryD0kdKxM&jIbB;nj>J&1q5twLMbvtb1eIW*d1 zxU$i3d9UuWuE$HmOfHv-nO=$|9m zZ)98K4|11{)lEq0CIHugmWn^h{Ops!r$swgTSzBkNz!-BMyJW+X3d;pZ+Ae@5JCQm z=MZQqUsmp;ij*l#6TQZbsLzJdjewrVe1=p@(`X7tVn5aG@J?6+hPdCSFYB07H2!8Q z9#U=GpkC4+PRrWhM$gGDm`en(({mu#NuA;~;i42%p~_~+vX!z?l!jjpoZ6*Wmi}vJ zgP-UWaZ+i4{Lf?J1HKWU;|Xo}v# z?6D-TCUr?uqZ{KBbnRC3r)m_rASTe;?4g$X50HV4EBCY}e$<=TlfcmC8Oz3{UFZSQ z=eD)S@5(=E?vmyhkk)KURCpIz>I44ztI#2!xMZsH-FLxtutlZUrK_>Jb3eRWi zA$GpI^wxdhpYA(S@LGXUuOK&r+`RjRl45ZI4d z(U<4B6)zrr83*bWjeBJ2`O#w#|0zb7RJ8`v#IPUP|641mRT)Ra2H+0AQtM}slI zgAi!;xlyrs%947|av7Vjl1wOj#9RIW?2K=Xmqhnfwx!fRT4Ad?*QVtpRf(1~H)nlx zJ|H#k9816)KrgvFcqDPoU~E)byz;}GG03N!%bX_P?jwmC&^tjm+@L&2l9nrpNTWz1 zrp*~C3Bu$_-K7dWvmryOIH`S7Al$@g;cv;ZTS?1bHZ;G1{yTUlM8L?l3&<}uT$_X2 zP3@s>ikNxBVk~D#OSG|BgdLE@Jw1Bl6S+%!X7%PQR1;Puk-VMmC+)xxBLY`X(G*YB zZG)UM9QK{?&CwQaSk)2{nosdhFd4r@GxzwI%1smIg8KMoQ~CsnjmR?&p%7X-+Gmgy(1~^EDt029u}SHEVk`LU4UOxB^-afP4zm%;fLu zb0QySehXY@_8_5ac74{$5O3?O=iG6{!Gsy~HEJ3!fH%2*>NUiEXsOuXxx6C$fniJ| z+nLx60{DA%=ObMqY5C3V!rC>DW>eaI9CPm@&m@Zz%I=*j-AU%-vgD!3MO4FzKW}TZ z9=P$y#urhZOPAHf`O>sG_bQsOZlV?trh=9(G zr6+^KlYTT@D;3^etZ4vMIaBcOk7Tzkn_{;|FaZV_X(_^$n!>?)7BEdg-Stf4g3uSB z7>M+q#TeqXh1VeBsSaNEvkDPFaMDiJ#O<*`w%XK`+RMIPRV%04C} znMYCklb8yV*$cO^rK&RFM$JCXpV(B=TnjI<)n{!pK{eWX=M_9+Ez_&}sqAdYcvEz@ zmb`_uY%lS69{LfT@Ex^u3I0?n*s*9QYNz4|1Ofi)QjH5N z4V(%NvA7=lK!D|#bh(uZo#t&B`oTLuYbA&RTzQ~bvhB^J-Xy13@L5$za%zP19#-%& zbS*=7M)@@{g;U3TA($I&kC%fpdH&DKl4y_N z?v)^-R|sSQWFMao3}iAU9(HAa8sFyy#(YV9C;=Si9m}2DW}axiU%1D^^8z|iW6S%h zG^0(Rk2RbxKBu{k>4b4DKQobs zK|*Pyn z_B(1XBIWxI;G|@(jd$R!)t!4_k%voV)>?xOgw%0^o7(nV-8899xC48|vzpjIYt{qG zT-;F4SNb^cJeO)Vgc{LObZOpot(2Ckn^BI{Du4y$Cwa}HUauVFPPGE0?`hdxDBDwj zNY@o|NAeG9ezKYQd&rH3tfdKK;56dx%4z035o)WljW>KaSLhzfPcu}5cL)RT>=lC8 zS=6>IT}fS9oV8J%qEC7viQyMvs}{yCH8J72;9~F}XeX?T_Zu}^VuXO@02IomwJE{b zzcjA}%p~;_mmief5*w6mlj5(#g1d(E{3YIU)7m|6DRO7W;yKTxoOJL+S=UZ#*J0Yc ztco2t_cb9GtMstVvcm20YA#NQ0kEwE@H6ITlZyrarlXeHM;oOl-9I<-f9>My7)tf9 z-y&mrBoC;zceU0Y8$hrGQ3e}o5Z%c^+gb~c;7pAGYsx>ELtGc~sw8kPrU@f4zH==p zjlTQ~=(Ey_}=z!jILTpiHF8}^QV-`h*y$@!9h&7%7r z)AC;m;!Az|;zM+2h3F|g_Za#)m-IihmWMbypVj{(shrhAceNqEcN{pEKR2k@mL@rP=EFSpo7PYaq-$*^*0JF_?g4GU zI|G#5m$Yd-so0ix3EcKgac(Iz40Rfvs`;BxNUzg-?)}PAd6o$}pQ}`0l?%tpNopU8 zSZ3#?x~>@!p!aKz$prUQTLa@cm@d&2p`Rc;OYcJUqVZL@xXyUl$W9r=;6>cUcB0~q zDroT5l4J!}qj+7pXC+8%iL+{rO+D*Hby}%^B9if52$K6e^5f^8N=GFUf>>nik(TC= zlE#$Px2x32Enm^Li|Y3ZuZMY<;AKl%g(n?5=v!#zq-k4a&;GqqG?qX`nHz(&2!|!* zD`%(jx$plERev6qRQ~?)<1=H8)25~6t~8lhxlpOOq>x!TX~v|PxiV&DW@?U_pe$Lr zWlD?t0%b0_lDTh?d%02WxgnsUpaLQS0`g_vpX>L#zTdyOxGv7?oWr^A*Lgl3cWB@I zr9!@2O9Xn$i{d^OB`8E(^@L2~YedmrdZ1m}ZRC4l7{5vG2K~GJ*p`{XE^N^>n1>Q1 z*^Ecr^7R@fR9f}U+oh8VIdVE~YJ4~^?JC;gM+Adx0?YeiA5KS@v?!roQ-@3V%II&s za_&rCL;|Y}S|kinxu#KXiZI|sLZQ`hqKc7m7w(kmK3XcwCw;av@{FYMf=vO1O1K@m z&hqbq^1dKnx|KxO_0q*G^@vN1s5{^ifrC!qXLM0H`1<&rxr|!HZLF(&!J^z~*kv2{ z(3v+y=LTE`7PqI%x30o71-XJJ^y?dM(lD|50pN@+;Xy5Dq0gU+O2T7)d5IPu!7qR} zrsE4Z1s9M!Qwkbd#H$uoW;@2=qTO=4~-*%Z*9V4n^xs?c3^V@PtAqDw^e zf_r0ThsH>W@VNf#;ogZ`2XP;3K&gadt^j9NVMQ|)Ul#K2?MqPe5VBC>FzOU5yoqj; z)(9^b5;CIDCDGOEZj0L`^);A)+t|{BxQjN0!n0kNq4h`YK#C1qzso$Ac&PGR5THbq z(_@W`kk&MwXCS|F1Hyi<5S^50ma@6YkQSRuBrekK1q-`;C(*MqgYZmf{^`M**spOQR-exfeI628U0fq6UN#u8H&raytQOpMK&Lhl%D2uIhCoDcY3>+Ao$z+9 zdhur}+`o{!+HX}n6a`+$)r)z%Zc=!y{6>8M^NJ^N3+h5r9!41(gcOfv1c30&8{UPZ zjfqG%3-l0OALS!5FEK#q96@?%=p;%8Kd?l|aKH)4?Ly`iN%SQA_N3N*AG({9S4d%4 ze@){XdUH*qp9kQTb5FqlZFp<%O#QmK0e&e=!U{k>U+mheQJLkIU^X0tUX^%xwBlZ8xUhyQD1mCXEaAR!snh@f850^zd`3(sa`3`a3PtT8SMGI`F7Li%9z$y&b4!q zeYfn#tZgcEHnxn*a0BeV5nh7A=MAi@zIgsS-oI^=wiV_cCZFeKn$x^?av&p- zUa-(n!-Z|&Q9n9`uhC$pbv&WEo$!7L8UQaEiXzJDtyW`Xcsg!bwC`6b!JcW$xyz3& zj{;U~55t(&=h2X5mb2tK3*W$yM6BVQi~TFjSHn((NjpWBO(`{^n=f)^q10EfkRye_ zAKOPODjjor^kyMfe`3~8t|Csgk@ilTpvNN;74=X)TbVad97zqYkaM12f`fRvjnRk0 zbGy)s!@Xovog#``K<+&GZh-dNDB^Bj-W<{ao!w;Ib$S8G+_4>CQQYOu{~)=qf>I~s zZ#m56^RHxmD4?djA)S$2p<#A~ zUllZiic6xtOI$tlizSF`&^*&=Q{=#NU2fU|c;r+OHs6~1g6r!5&SbYaNDM2mUrV72 z%q`X%CFC^^L=D!ljQLqm(XbR2flFHicU{5eG~sKSEc;jn$u1tPlNZTl7SaMMn_?^rVg7m@Qt2-OU^uX}lu4z^KO%kUZTx3uO-Q*LB9p@ze@J>U!>t zyt8KL`#!($Yb-AnBJFRszDk!L+RVGDU&QkdfY%PeFL-FD$2mHH9K9hop+(}a>*F;D zXaF>Q5$IH5H&TDFemC;?D!J_rr-8bCyo`LM!{E|*BFx{WxP*|k@F74@IN&q8vI}hj=#_%m876@Pda8&Q)Hvr2Y%dHKe->Z8e~e7YK#aH|eaL zTNu;EYYrlBxO0iLP>{xNOB+;ohGnyj^6ZL6ce%s8Lh&17V6whDJ=h)f0L@3<6P60> z7+%J*2PfTy>CuXzW)_VbdYg~-3=6U!<8~u~)z?T{LO3tOP-{qRz(oEs>Y@_t&TFRr zO%dKKRrDVPx(iJW07d<2Vls42KY*KM3auQ0oLnfIW1_NL`I4LP5A98-SA``J&Nr<` ziPamko%|clwDeK_o&FP_*DoDGxzknB#l|BR?})-L7%nUc4Ap%8Xa3G#N~b{YR0WBBuWKp`usVBviYcf>UBj^L&WDACf_ zL6Y?oaFbE(&g^5De`X}5!r_x5y%KJTDS4|_c&>wJFxcl#czc#k+Vb)s+-GPtwl@Z~ zV|Wr)5uxX10AP4ToUPJRvu!Y-I<%GKKudGYFjW>Ci z@5Y$rgWna0P@v)Q;lYs4;y?HF7$P_~s-djkn5odQ-AEV|aeAqO1kBIgF znBj$gt}@-XUqzsDshASf_zfMmY(NCG$#|6*?j1lk3lM%T;l+;a%*S^MqqFh%R?MqY zDsGr_J#crr%$bs>hG0k;aA;GXps~YNL?yCLh-dR3xgM(25VZ>UXCK#o#;f8Vh@t~U zCb%T&RZgIhADHCohYLI;@=2O7q49$P(8205zs|&DsUNtfl?HY(;s+fPU3Fan8Ts*( zSa_!u*>{)}u2~hUho4i&O0S}K^U2QJUS=6KwXgi2!E+|)w*Y($8bkn~@aZBB9dp#Z&G2)y;Nt;B69>f3}L$L;FVR99O}(xvpf5R*}n|#xaxdCWCr}?re{p z`z#M)Y{i7q`1DP)1w9c)>Q(m6SA0Y6pOW0s+T)PY+t9MCK2o>Yz_v#x)>At=IG0~u ze=?zF;>ct2%_>nRZ9g51OJ*5tOB#%6ff#Oy13_VN1puMU&g~~d#dO;2_cNR27V=+W zB(=GH{8YP8qXPIp;;cE9@rTxxp#wL*A*l_uuSGp*1=ioeUE`(w0VpH`dx zVwcv7@G??jPvzYepRgqR+G+3_F`MZj`*Bq-d%Pl)pnlUC)P}N~p5Oeb(ZR~mK4KGw zziXHfRE1#zu?A+UCkDf~D;n@?VrFQ+(1f(uSTWDnXbfp7N}#Kv(nYneqsPDUX;y*C zrX!BI%u(A8L0Nr|P>FQN)wUy|%dcId^jVxJCHks0sS`eY@_s|{17UQoy$0VRKCyM( z7G6fMtj{gb*{5>J}x7!CA3m?nE6M=%Ys5 zvA)TB$a>=!4NDu3s~9^Cbi@jn{3bP2p%mRMIYi3xOwPxCcfe}WV;6fZaI}f8*7LQRr0X0(SE#dM65V1U zq-88Or_QXnJ30xL@=|os`u9>COl%71@dt$()_{^l&}?!iYy2_@z}iI<32S#fH zq0$%8zSF8ywfjC8QMjMhi}5$s=LRDf)39LK3*vh%`X#TyJz!P(MQiAT!2}W7dc0eN zO8$TXR<6Pxn)-J(s*YEVF-s*d)b}Pc9YE{&^!jm9GU-pxo%{TBFay)<6uxPwPdBd= zdDZ%CsMu8+DyH&-jmT~TPZHPcg`PGUqLHsQp1uCG2-}wU5sUtT98)l`O~+vc(sO=H=c%!*=Qz3=w-UL^NRi%n;;c&s;Wq4-Sy?IcWiW%NKVPDX+lO_ zT(`xWLYAq`omG5^_*)aC5K($`%)N><;BUHLo30E|-M^7cFBzmg`d{YMPtWUt;aVEwT%b*|@7qJtZqt#} z+Pk6pf#Bo{vt+*jQAjD_b4~i89#aMXQ1B|>!SB^(OE-I9=(Zp8{)&(Lglqm=2f>Ki zJ$H?g9pl}bqnRj8QP?ChTPeOqa>s^vI+zq>lf3-yw2paESN9n7%0ms>rJGtO-64f) zZoS`XkgiR$-3{17Qis}`i8`NQuVWz^bxuE_*})pLn|^`cJltG9eH;4TGWU0)%1}*u zE!p&;ChZ@WP4bP3egl8aj>C?20nFZ^%YHDDqqrL8`Az3!m%dYA0O?!%>7+3Cz--mv zyPQU)z16U>33~PI;TxbEy+Vfuw_Ui;u-wANZt7I3nOmu-yTKIi6|$Q=O!d4nt&!xY zaW;+JbmDaHq`zKNNvdBm(l{w3ZK?dU;SlxjUERqx)u*xMvBPPg4{s!H`DNX=XP#T! zIJ7;zlF^Pkmqze$1!ZV};l?j1)+^+x&7_6Ig;tNSW@MkX!p-g`kOcDuQ#XL zqe{U>VnQRb^;6Hifxf_i*$-?{Q7rR<2zi6TfTAk>MoD8fyI991WVCwisn;JR< zY6b9ff}qf|k1_PUbhC%6S}FXR)_RPjH4rW-RO>(8a<&$hd|GGaUi7Uo>ZkV=tlLw8 z5qqUR3;2#?r6S#Z^)7~sd3&3YtNu+HoD<#kof*GA#_uu0uI{7%8i(vTlcHD7G`xAg zsCSHPQ&;3M#mqcRoqjmg*Jaxo=>xOqk!4HeCv4&;P1-iEb~FXTwYy-ckA8l~f`_P8 zIh2_npRj(>G=kBpe}>ugy`R{16=G#w{56$$)ZG?`Tgi!(;jv1qB>zz99?z>*>a(Mp z+)IR)CZH$R(arq*p5tqE<0+fMADfL$>pBjZGY2nj$Ej`_q6sw_sH%uMLz;JI$wrvm zvjk<%t>kL@#+jGnr&=ygR{Luf1sUD!d`|M|xMQ#A8x;C}pJL%KJh@iT;g4y^* z@>O#C{{96FYEv@P%M4|e)YeASx31jhSn4H+_g3yb?sx0O0(SG~vwf_TI?}-E&kshk zs-ETc9uEgJ9R$3WzOlg3pVdr5P+%(2?FA$h@L_Kfs_x#pufE93tt0Ht{($f+o_*!( zCJxNZAkxH-F{>|QSWLxoTcFv@vIo&#idU(-yRK6c#1?>~rJo}o{Jk`qRmBkGF zzD*z}<@f#Tw5xa~NJ5{O{CME|#4jC|ArsJBK6Z&FB=xhMkuFoi`$bxuAkXgB=i}!+ zr6=N8Z7Rhm+k&3cDj$BvQ7oqXNoUpuxpfN`XZIM%@}!n==Y;R>&P|^Ty=x~Dnmx=L z(Vze7@^H1amk+Fv69Q&|B5KU_I50`!A&^x2*pBbriW%egxJzx4I%GvE@9+5t-5s-H0 z=YFqiS#<7yFtFtb~Wqfcox_@-K_!c4wjlI)wE7`qCno_s;+IA$+1L)FT9=MTEkzA3Xij zHuVE%F!otOSzmgnA*HTa0vfzs1)npqP;lA*SHC@NwU%RpHhrX8;rQK~Jxvgoet?AokS?2R|S=asZWA|cbzpE2>BuH=XvM+G2H6#t&<`r5)nLLqj z4^~@;wTz4O?08CeXz%IAU!QbjdQ6Pn6V|E=AK%!&{IKJ2c4YXOt8TmVg_)Z-`k&l$ zzZ$1;k7gL$rxAa)9=RyT0#ibl`->Y`pvTe!w%>0rw6l=^fi;08@5Go?8L17*(aC{&YF`L1GDKhea;M?7OMjXLXp`fX;> zp00z77=Ob$yJN|ilnSJQi0fC!ml3jm?;mYc)x7E5hXYnLmsNz6Bkwsjk2}P<>QA9V zqEUez>s9Mfj>M00QSJ}F6w1UjL)#B-tBt7Dd1~I280y4y-7{mCkF?c0Jc zNm;h3^yphx0}YRPYUm(WOptr)&#lbHbWu-D@&D=YKS1mIV`}(Tl7BMg;k01|{6t4n zS=a7?no|$ur?Le-2?l9=T5itm=*Is@oPH4ao2**$oU(CC!n#kpTmkd&??%PaCfWoJfKibdGyO={IMHL9zLaYTXc88dzIUiameqQ3joRQKjk+gM}DYM5ik ziQ|6PU;H?0TR_e}Txu7UwJ1{axqtVC?_`r!ne$$cpOcmbnE`6~tJK+YvSHUF-6Tb% z?zGk9DTZ=~^`~1l_h9>sX+G^O4YQsrmIEtc&w?;^d!L=St|fy9hsK;;yQ^x}cExq{ z|G^(3^kW=2U&?-q4vt?3rE~2_^zObvl)ZUz=-Rig(*ioQH z^HpIb8I2(g>pII4pe(hUUD5n~M@XzdPLAXP{jC1>t3Vn4xqh`7av1;2q+ljTA zAlT2k>o5N5w@+%jU;)2&EMxP8)azTMe%I|AFPWEY!xvrd25R39)%43$tc&K=M!dy0 zKpdmv7P6JD&4a>y%`OBCn8ne~cT~l$Qf_ zrM;)P*+Z+YPVK$+=9R%S*bHoFCg%fW|3XGmZ$pOQY0u-`4ze1lP{6bh_5TsQsq|Qi zdP7pEj(BN=COnsFFNo5{V6T@o$kSifnj_}mEzcK@CLM_Sd@Eyf%BsDvv({IoN5{6a z_He@DhTs#jRXs|#=e4eC{X`Z$L-0w~H+oG_Muib`jb<_QDy6a3z4+CA@#j|5AsUxmV?=Mq_edUk*YCKV zJaVG_eU1LcFek`9+slce7zj9U5!ruPrTC`p zJQLY#k-Dk9N@_)tBXs?V?23T{1 z_f)DsB!gXTQ-|xSki)tWyN}HBj+4u>7Wi)>ZB+iFY@+)&ckpFxhLh`0pC}D}{|Cp0 zx0v5NS*JwX`;#nNSB*AnX?qr%0PgkC$W09s(ZaVjH^I7ffuP zek9LWB`^ke#f@sKlquM57@x327j?x6uc*PGo-h3q>n~xU2QNwn{T3|JZ~rA7>Hlq4WNmmy6ku!;0`l5Y>0K@r2MS=;A4riu41C_ zabWY3WB2bq+*B$o{2^q-e&dfhv!k*)X{<`aS2KyPMze|NIG9NMUnRVCGk(rezqB#*osSIFvQi|DAw7`*)VtfD38~& z`8w8^-Q{}`j?vRPH$n1L^nB6vBZf3Ton+Oy=H%3z?2QZ}cHFMri8Vc12Y$e-yld?k z*lck#eqRkezKYjgajnjcDd$iXpwc1zbmnFubOd7ZEMdxV$8`)bZM7b%j%IM2U+Pqo zXw80mtF!+VJ4}3+GbwK9*i4kKH6%Abzogh9Z8%&iE!y$pWvC1-5?-EYjZ!0ieXqY` zEn^@vHrP)D?L8#}vbG`6JGJm#227R%QxDAB<%RQ+d#HUetS_QB^d!d;1T^@!HN&^+ z^=bxepF+M4#?oJ3q0tO4yrP$E1oU&`-=`Uet!G5v1dcR^Zs*+rx*0K zyur@3nM>B5uRM#?HoKrU`|!TDsiGy-%tED{!l(BaA9BAA?OiT9%+OC*F8I<~T@$w} zhBW6uOb8;cgB;alzYrCxabLINnR|fR(P&cIZ#ZS4VGNxA?w#+x%QIlTgvHZicBvb^ zS2KND1W8Sq|J~7DR&cbZmmA`#(B)^ao;IdRv#^;DS)4lbl_sXs_ zK(lel&#eg?{RxkZPdu+Xuwj#-j3B?jl{Kga{S6SW3g-rI5d^}&G>yGpCt)hVIkm=EJSvd1tN z?U#WU8CyI8dO}`FAM}g?u|CarwsASh!Yyuw{qO4r(f|x8Foj49mWpz!GsMIFXj!2# z9oC5phL4t?7uHq7cl1`*QeEfta6QqXIG?CwiPm zCpJ77yebv(&Uu@Q+j_uUqBM+S&u^4@$ApoO10(-)m48G3jG~{zocpIr?8X-_KVd;V zC(vu$W}VdyepmgKX(<##2>S~rvMnJ=gBOi^Dp6Z)iFi_Do!|>*9XyW>ZX8IZb0bj= z4E6vx?-e+*Y?r;0_3H1^MzrGW{Uw1UY5{1xc00pw82h$C>H?iPyPmJVoi2^0P|T8A z)}Xkb=&@Pg_D00^fM^r3l78qc4^Korn5v(JZ7+#M`^Iz7&4Q6p>6-?rDjF4?%xqrG zLDNzlOM7^DCzSN1^AicHer1ad&H&#wk&4dM;T5H>)9mFR+?$^NB#D^mT%x8d_p}E! zJuWD0uw0Yn^KlJ$&KNvE^b5I@RAzkDd^?h|#~Lzi4|^pbuwOTv-R&&q23YA;E=^N$ z-(Pc-ByG46de)j!%PHjjX=R*wYlnF_jckt97u?ypHci@CjWVAtGtQ@*zbr$~z$j>W zV0cB;P0%NW1zp0?1`(eU8XF1`V7q|au({s|D)LeS!vlGV_t478x+gq)E}nA`WJSJe zpKU*cm?SEq+4daO@pGUz`4Fk_;A|PIx$-DUz)!0WuyCK2u&7~dqXQWS(ewkk6 z&bxR`dG=|D$M$7FK$M;RlHGcyREux{izp8D>Clxna|R?c)Oj-G2M!J8{X+a0Kg9UG zBEAz-&3G}w`@1*`S}2rj3YT_3^#~=T1E_{;iE@NLBW`qw{s#P2BA3h0R-D@z@vGI9 zwo#5VI?L&ePHGIuh4Cjw(fV&MS%f}N=PDOxN-q5ykVp#9_VJxju5ek`5&{cB!_Yyy z@Dm(1;k@`Q>?Fh*Q@Aw6v`6emx8a7?kF_YFH&`EJ8F_NL70?MBKy1M}qwcpVfQPq^ zgkPGy{a-Y`#<&wV;x@cx@W>f@g}bi;=e+gphQxtLL;R`GsrQ4Pg9~VV>kRRiI4u1E z?&hD_$SMcMnlI^Ag!G^dDsU8=h(p%jxAx-XAc%dq{jCKIJ5hG3)gb)U7_?y2^9^3i zWN6NP6sI8E2L^tei#mdPcP`u!J7!tP z-**Oo9qB*fx~ugKX(?%#@gV$9q$?xL8HTmW(8CJ!5G1xIMz!TIc`4Tm>BRMMu%VDu zZMMeSAatewqUJ*kE7)Gv+_Ua;1FDIDmfq6pPH?ymRjAviA5M5{*PRPNn8ELPn>JC~Z_`zs2bGn<_PFe` zkvECDpwFNM$83TrL^$3`tP>RSs@XP$;B}lGgna>C5^-cp533$#D$T{|_e-qxJupJ- zm{GR>@M4}i-Db_+Jwl<$fukz51pB{kxkhN?k3@2!1sJUJdNa1Kpu#?sEl>h77RtkZ zSwGqGS6_}h+>`bfYFMa&C zOTrWh!V8H?EZ@KB`NJ<;s4c_XLwVH)6!VzH12N8wf_}-3D4#tzgIvET!v5KFID~=$ z^7Ep+*kD2HPb>5NfI`xFh{yPxoZu+W(PPFDbcD2&dwc+D1}vJ_lY%5q2KAL$ZMQ%d zR8P=91`wS|lp%%-+imtKZr!|(b+UYqGWR(39r8DXZVDsP9eyroupIvTY?RHzsAIuz zq{9g7byrIDV3`Z|?ecZeQYz?p?VtnK7aU(Jo}>R0z;eq1xLMdcAZEFr7i*n`1Uib^J04mGgeid-jzry0TwT~h<6Z*auxjv^<_ef z!fRhwfUfqv28=@+zkm&F%-`#fqR1Dx7NOty`qo5}#CkxQLa32`8DD-Djw)(-BjW<1 zCbizRhfS47eUbbE_>PD{myzf66Zyu4pevUAtK}vOtt0aUm-SOkV_d!PCy>4<8{Bac zM}AC?)pmC+9##qr$oqfL-=XTkNfhx?ZVpwB$gf82tK_z)IgE)s z;Y;W}-9$cpagjqjRz%$NmnKYKU#2g1m1S!@a~S8~76(G$X0_PzqlzShBBClTXT6d9 zIc3`nhs&6*zzJI;wxY)!+;uE~(GsJ;ujCGyu&Eklpl~)0smH`+G_CPd+{U(1__7m{ z$(Dcb2ScSsA`)36sfOj+8|#1+^Os_V#@2CP1oh>}>?UGS`2uDs5{<_RYdIh9kch`B ziL${>DltG3({k2&rdX8o4HG!6-V&%dSz5oco?cF3)y@x zo}*j@$+{^RCdEk1srOpgy#aY~U<)hh_2QWfwW>L`*wY_ra(C|RR6i^z$yJ0vPL~PY zHYcP7RsQ8?XcToyjW2lBg7)*-Wgp3dGW~TDtHwLt1r&8!y41}a1vSm}n#|6MQ(EO~ zS}o91bP_K7B6;{XQ5yQg)ai(7eqOrt{o*R`*eZPfr5m)AA^9rG=@-xy9Cap_fs%FH z!~LkuMGyIbhUIFY^j~JQynyZ_d_aK!Y|ytSgp})3hTLUkq2llRq{5~H^hNUaqnx0W zz{gJ8RiQZ*1lk$;0BpT81Q1oeEl>v#KDYHVr7h%_k@{ z?e*0qbu^eksaxFySxLQ%iEHrjVENe5Dd^8&N-U8bL1?8sBY6{8aqzJ=l#j>lJfltU zqotQwl!RTLDF|b5;7YyAWb3EzFWprU{NAs6@{3Y4MvwWdH~+PG3W{$0QPb+C_#f*U z;-b=USlWK3dKMK#6x~CF&Cv(y1(aBI;KQhOVsP~0n0m(JX8LDZJQzf4)i=%!_Yk?8 zfXM4~y)`78V$BgSn#|)Op)jxQSo<}1VS+rG)6gQG0mhTVX}v;oAJF3wx_q!o9jv*> z20y3;qCL}w%m!?vPOVp`kn5r-$wiCcMpn%GVoA8`cHS;9@{>Jk@}(OPy7fBw+`fKd z9+?J=k4%)9^U`w=t8)Z7E;M8YF7Oj+zm`6nK7=7ewFynOt+vyYu`;K+w9PNawUymb zD&!vn9DLsFl_8Cq$N3OM!%p0g$Cja@yO^1yFG8~1Rx;@4W;ZQ+Fd_}2mU&Vp+tXtoVWtx*G~WCp5ca7d z!5X!qJpCi1OVR6#ImN$7emseikzfFdy*3KGTFjP1cDRtv$w<_!K}xfzX88&1gDxGj zJkC{Vxy)MBT~Y?M`wRB4+B8911h!FlAv}(x9)ir9>x}vijyeXSh7Fo0CG|kX+8`mf zElu(Q@EHm%=5+H<$Pie+mYddIfNoW~BBC?$dpsM$ZRBP(xLfCqkEpW$*rwVG&V}X& z+oFr2>2H^3__4uXv39->A_eBZCyb? zk$b?MIF2{Paw{BXism%Qk=xM}%a(${KXFt)YYd6Dhpnpf1i#cqlT&1y349SdLg+iK zJJ1C_`aX$U30oH!r5B~XEfK96z`2kk;HTydlxItt+xgLATZReV@uMy$i1wZyhS1rj zDGT(Yj$#EG!e1mg1SEKAk{jV(pHjB$$i4T}ZAl(?EF=><6Y^bgYd%b1(R3Q%LVTQ? zs4sqp?Ero^l4vwI>hP?SSa9Th!d0*xwj2DFHt)Tz+VISPQ=k}AytJGA9dw?d&r14U z+dUOtViktUkKX#D#H}Ux9k!d5ppOc)dby=lrpw+(wFAcxmE6Rc)=tRVQR}=9y5b#3 zwGeA^7AU`|d*fI9RlLBV9Xb$M0}PkSaB2+Fh5xjevWkm$YtmutUT@ko>((a&1RZkg zv4aP+{QF2t%B?c6`C(|nqf@vT+(rz0$1JV52sq%fR@}XZ#A0isvVygN!xdyjLcabm z5seeuw#K1fMSksIr^E6!AMFQ6+RRPOn{OPKQLqLKEyNY(EshbQi}wZ*mNLL$6u&}ei)r9^9ca2)1+NR~^}_ScZ5 z3-lc$s_mY5Oa$sO=mY66;1QT{O4fa(abHAzfoTgLcUZZ}Zie)huc(jpY(1F2MKghO z31*AFjM7!0%Fx%xsHzab2_%NAO>m|^Js0sF+p`~xlXJJO@hQI6-D5O67|4L;8O-6? zyRnpb57sGpJ;f#K=VDvHKi2EQc;R!|`FO9om^-hxu?Pwfk6Qggpy&#AqdtVW2;@m0 zaR~s>U(h#H+%6gQ3g%bAu0R*IV*KZAsgPF3>6ejf^&(m$|7!I`pK_Y_MAR_w;Y(iAZorULUQ>d&&-(Y)>-f38rQh-d!W_;O z#Hk;^FM%*TJ-3_m994|m zbWRhNRH+1ifxHe-20t3J@?atv`)I9Td7BtuORCimS$l3y@nbuS_6w4H=B@QlgP#JH zGOU27G3K-RstijUJ)NhL$94C{7~ff2iMph}^g^5;VU2nQ`WihFvym*w7(AU7_R7C4 zi=H2e)uo@&0W;ieGIq)*;qf^wdj^Nh;pjd>D<&9rtsf*Ic<4tFABGnAaNzD8@Z#}A zwxXBBI?l3wL01sM46{ynwtQ^pgHRcCf?f+k7m#B71v$Aa7&KO=1iT^sCLykv585-F zvP3ZOuYhlIQCODb!zg8GhmEuR5?t5fvoYf#5*XX$vsmH-S+>NE;1aD88!@ZkSI|U| zCiIDU_uGrDhk129yTBt{r7gug?%(LFLD@F%l|>m$`D!kMV-bq&Z({EhKNWXsy|z3h z9^(FkrPmXJx$Fx@YlxXP)Spa!W@s zTt`52*^tbuuWRr-SYM~Q_wX*+kbdwDhZu**V0hlvbBla?IOeDoKSZFDK-af~hNHP* zd-`8dpT$4HT;<`Hj^RYYdnv zx~tA2y;tnnB~{YaLu3)&SbyXiF$@{rd5$LdHcW|qmhSs=a9|hhIYeMRIV;s;Q5Btg zq7am{MPD%M_S!#)r|7vS(07rGd8H-`wS4z=?5pKo_!Gzh=z7VBF38+(I~hk|uFs~C zNO!g?2(K*s$G39D&q>az4_a-dcH4)cZI0_7#+PNfIev)vckG;gcl5Y9L5}#L{8SjbHK^EXnOz2fZEtu%vABMRg;ND2Vo6bXnR3PUAv0 z)Fz0mU14vfpkh(QaL05S+4>PRJt@Ai0taR+ml(u0I-N!8n}?py)C?cCv%^i8_U94) zZ`??~6Rc$5qzh0sIu*@KZimD*8oZDOw3~HIUy&m}fD%13r{Jxc#A!a>cda-%N$T3% zPL`T-0q3lh7gef7irfE^bI>S@{!wBRYHc~)cdcd!T~o=PrlN1K6HkpVvywuPC}gH5 zOA_REE~@cVq!S1Hy5X-((R!j&r3E$SQ{`ivtB>WQ3z<0^80ypZ>8^TjqO`j?u2I{8 z*goBXgZvEWEcg^0GbH;gy-cxTg%lm%c$X|}7-z_JOba+HW9WCCB~H;x(k?Z!=$?auW*}u_FE!Wzr4gW!8U4x(8(@X`lY8n-sa5 zx#?NQJ>v?9ZFB|9M)UzHWrLK{b+i^MPfE-Wbxdyzp)a29_a;~TaW%6!ViU4w`^-nS zupI*xJ;?N2q!yA^PwDn269SpOLPMY%xZoSzxEBLtxGpUgME{vjfK+gJHL)oAMfT9=lb~9Ebe*403S>3l0LCy zQ?Fa_68$7?z~ZLCPV4@PTU#LqFSpz#S_4{yT7pohU}tcr+;p-`vcZYhXYZPA|NqK$ zI%eR(vSL3?e(rut1V{s)Ae&5weWyA6f%#iryC&7~|8LjH-p&Pg?$gN!gv$fJ@7g-N z#_kATAv~E)4yn`>cfn+Soh+v23>Cbj*{aF%Ru!nil>h55*uQDDCH=3xV8_5dD`seh2>s=8)qyM5SdvoT<+PXbG6n>O+QSo{kjQLe>Ku}HN?QD?s!qBorJ|8neb zqQjpLQqIX7J1QaL2gJ5S!5&$IA^4wT=cj<#5)Wk`@S%AwIjp~=%fb4rI3_STXu>(p zb=W}gCE;^OCgw;mDH(7BkUKagL*0>xRj2KY#D zr?St7J+K{Unc3&qwax>Q_t9dx)(`Teq;4SWE+QP+4uHjBzJzV;B60xXzmX{BO;2%x zJ*ZglJM-{ZV-9%J7U3>0Jk$aQiGQ&lgZ!_xfOj|TgF;oUb$vQ_dbJJ1Lpaa*0wY-a zp77_g#vvFyHuQgT9l=#lHC;$g)a5*ppDT;{;y?Fm!xu%tAG7Rs z*dnYWl6}O-1Mpt4MCR6cB&$NK)!>> zx{2=wgx_#;V1k{_Q~uIEom*lAS)c`bT3^ut?~2;%-xC81fXd+}{_7_m4O{+ZLDOX) z0KFx04nT?6Z#HckV&Ps9ML%vw!nSW(@+v9j(6`iDU6gyU=j}ncHaa(WPSc z{#hkR_%moJu0wMp0sYQ?6*_kap-=itJBMwA*lCjkwj?G3ME8*807Qaf!!rO=UXeP$ zO>CyDG3D$xcjDV#3p(r`Raw9g9P{<#;x2hTuwJ+mA~Bw$cl|$?PJ;2=e<~eV;##dv z9JEu_d+u##mf!&(5%o1}=>VXu_w+cY-Q>nD1idUUX^`;zGtM5 zt=z(KU3{`6zx_p?`W!{s!?^%qk4rkklL-Ik&{6cxB+SrfRAcMAt2@Pewdb1I>f|=e z_t;HW_KQYP%l{d4Iy3Of5SHRzdv;9Z++HZ>5-pa6;VH@U%<_O8j?}JcZ;_5b0jdq| z(AHUo@8#G*R1_=5^CHZu!vF6AEwI&BHWK^BA5GJVE{Lh>7|B!X|aV;+X z|3?U^9d$BP+ac^Yj&|swR=bfTNev-Xn+gdv_uZ|n`y^T^ z<)(Jm&9$1^U2FH-@8a|QJifo*QoDA&U$6J;^?cXulH%MixEas##s9>e!qj^= zZ3bueKgoAV3%_A=Qk*UYfkUlFhoaC<@5sy)wYx~?jyJ-sKS0Tj3Nc8zAt>$giQ{KO zZ)P_j^vN*QRPsezLlSkSWk55GYhS5`qc_FKvbfU~FL49D^A(N0&T%b@jzlb^F$E4q>z1CU60h?@gRBcE`QOm-pBD^#`OLA{=vHr)zPK|p$ zybdyv5%s&EKB6lF{4y>?k;i!dRCQE(gYsd}(L6&avDasgG-xvjpz_!}e3M5XxHoW# z^$8;ABrYV+Jdx<-;MAByPWJbao7vJq-ajAlLy$hpizJ+lm%w=`Wf$Oy z$u^}>|21wsDo?RDFVa3SDg&d6AcG%wgOlGDX70Mj_Ajv{elnBpw2a8c-75pf7!Z-O z5eI69MZB!Ew_0kXGE@|jjPk?PJd;ox1~}T{o?vB$MGpS4kCA5){)@-~UGxxfO*^H- z(@$fpmEP5O`UY)PI>6c}&*Zc5>w^Xn6PYJ2;nOvLspnO7gM1uc>r{a@z~fzrem2jM z?aOT9Z_g=n!6KiV?nGd*+bD{&+Ib23AoIqFF!@MLyFnMo&lUQ zd}Gx7xW8;lSKd5u*cN$KMyZ!pEe5G66A}=Qqr5@fG_k$!cUuE=!XDBwmAa)xxm&QR z!;L4UU6|-p?i{{aW=~FsNjx#BbS`UcWToQvH#Cm&u*E2xDev(u1p+g?d0VNf&9@?! zwujf+fseF)dKSk1h~64>D3lu7;G1JQ+e300L|?YPjTgs#&`gO>@zUnawbZYS|5bwk zijC|15bcew`9I+Q2aF@}mIaMQ_IjWj(he}1-2&`Epuy`TW&6YceUn87q%=F6oQh`^ zX|U3}iq&VIvor`MC`%<5TIi;U)x) zF`p%KpF)?JZpe}@PaHSWZyEipQ;{9^7}b#m-UrfBfLsP0p5|BFh!a0zL5jT13xq|8 z__Qd!d)$)o#uU>>fEuvB#TiddW!7S)Gi(wUv$s3d)@-M<17SFgfdadUfW`4MAm*Ci zlf0;x6C1Ugx#@6pOL#neCP$pYg3EmU(Lv$O5qgPricK=rPfa{F5~#AHu9u(Cr;BP` z%+rNM1U1yWG3|ep!YEjQ%hJk{+Ku$3yk@_2lqxLWNIBn~FRrbc2=JM7z+^pYK^Kj)=Usq}u5hpEGEL$in^Q0j+h!k5mvIyWWuqjIP};|zA&5<~tGb}qK}OzmULfrpqrjfQZc$HBbtW#%Qw8QDToPCWm5TNhKM7&! ze?(8m_T~rJv*egsdsrb#{|J~Zo|7fEgOn&bq1+D5nv+Fdn^9sXo7k5gVsm14&NfOd zyC`v|sM)ioT({bpRpjvRN`z-}XoPA1E^dB=|(@zMdkRiM6Mx z#cQ4av8VAP;UAwWyF;Rmx||-_LqAbk`p7)AyB@qHxF`e3Cg9Cux%LfN%>%=e8jr6S zYFP#?A9xW}q)%SkFg|JUN_dR#`{SeWv6SPacVpA2sb9nR+nXln*>{QjQ&zEBGI3?dVm> z4q|&JmG4pYr(cDhJob-Z>?m7t6c6b^Y96JQa?+Y=vL#~qGxmVI^w&_2G5Xk7IXHzq z#j^xF{r?WLuHQJ3WK4yc)&$M_h|(>iw$<29sjCpaH%#Kxy>S20) zer>#fjN1$>Vt_t*CtN-ab#fzC<`h*C`qc_c_=5HOhb?pFJtdfD^hvVi+vCsJ?hwWu z6rD53PmvgiQ)53R+UJ!^h)l7Gvv`TK zyvB8tq)8z036KJ3Aj%%w<<}GS{d=el+9RSgrt%%FH{MC)K`>!d!@f#;ixjEcKahfv zcv8oI<;K(~J85|HnkzBCmML$L_GIAVoo&++w7!#^xG)A#X{|2JwmgJCgY>lOJ{7C? zXqYCr%9B3_FZT4?NqZouj|IQUz1{uN*{^=EogsWp^O+aaTWxV-i)Y78aQ_AQL}V}H zIqZ1#uxAwNlTm_-BU(RO^F^|4qf?Y|jNjb#epP>4%;TzM%tbJ27LWq{#)TB4v#Zib zDM-1R8_Pg=iRZACiw;NvL$7>>I!yvldRY<6%tsev?r0b5HiX~C)_&1k z8iZ>d^j)FUDFN}>=uEoU-ndhAmSfHv4%P*ls8AJDe~*Mlr{(D92vQ;Yk)-)G4R+&? z5|w`2)l0CVQ1kq#T5o~wYVq$`!l13R4ODPi#s0vbhWZZ(GKJHk&!CH+_QWXRk@E=O zXI1L*tcUg|&mpP^6Ys2@YmEwbL~BtyVRTgqWG(UuQYzObgYZW)7hsGoo!E{?e4kO9Ozs&e#+##a>Ge?Q zAFIGC;hi}>nN+xk+Fb!$(tS=)`AX;Dce7=)Ey^TFtH03JRs0^ZpIV9KZUo(Khq_M9 z)sYU(qlt*GX|(5M>&~L&O6x`5r||*vKcP!#J5Zn4g2m!wJ--?JA}?wufq_W16=$3& zu-E5wXWxcKD)eksy?hyX{dvJFq%`?N9Tc>>Cj`8g^esVtIcm%W^6e8Z%FZ&B&9q$7 zXTKMWBM52&_1<^bBUkZ^6Q0i2^Ew}u{&zS944;`~S~B*Fv9?oAtP2}}ctr($1YQb` zuEIPp1e11917P#*OWXPhmqh8j&BK`x-=Dn!cUeb04MB08iFo-BrEL@LA)5KRUQC31WADqiByl;w7(rZMAX+oV zmg4IOYq$o{yEq>enBp9|8&Br` ze2*T|jxMj%?+>qA9mc+d5q9r}&dGGd%vo z_K9gWNdR@y7&!W)^eYKv*hQVjb5@+LXY}6d{v`WDz^;JM?cmzJJ%C2o8&Y$Z!hJdD z+FK-XfUVkO|Ji<&ds$TNAdnQZ29_BQmSj^;{3RF5#}nSitR%d^t2HM|$_*4O|FXuJ z{D=q>fZ^Dldz9&1t@{SU5Uu2G!$&CZjv_k#Dmywq5w@M$agT>61b&QsuPz*1rr7m0 zK|fch&tYX31_eNMhC35!tblwAi=5&-n>GJ zE8S$?Hhh&nzLexB2elk?7^#bF&GBRZZTZ~$WV|u&hHrM$LOKXf?cK{~ab$_{c9Vnx@;~XPC3~H;uXi2 zCvMJRF`9a{76iN88wvh6rel7PRGsWz+F%GfVjJx)&>(xO{adco-EAZo@>q5z6H)EU zCTs=Cijv*?@iOx@>t0|zJxX9WuX02O@Y~#jL@q~Hpk>zzV@iYKoi>@pl{`+d;vJcu zqp=H!>D8Gq8##cTaz5j|;hgcNdTz~AC+aNrM%>yJ zynhns>U9Cf*{mvQo?WO@gjKEEwda7j(+oMTJF$8cA2GyOVIkfvn?uL1j8( zfHVAtKlt-Mp^u4x4xKU5w&CN!$#qD4!n+v*PP)afd(V>+>* z2{+96#TxsFTK=TrU9lk1)Er$ux$@ZOBh{p+SzlmW>57~`nQe@O`f$0qB#8awgNp;_ zq^*U0ey!ki&C&Eax71Pu>vTia7^KZM(c0=CS{+ws#&NE*yR0SO6Wec}^hcNfsJLbu zY4bvkMVX=(OnjDq&F`%)sRiuewh0AMjpbLqpUl-&tvfb3knez}jtdWJ0~mf{Q}2AI zFHR_IWzh@flYz1f57hIf%!AbmB>&~{mZtjCNnVEHTeh1&p8nMu`*dv#_`9c)I;-sg z#(|0Pd)68iR)ghGt^An27n`m25xA2_td^Ff1c-mYE;Dc>m1|N~Ok1N%9s5ctJ_B9ktv>JGk0mZ0`e4uWD<`Vo~HE_nF10@ldFNT7klWKNrQc zuV`3#@~J)uy+YGUaX@~u*PnGUt&8;G?D!U2=_iM}d%C-+Yvaz9gO2RH@QH)y4XCo9 zu2lJt!XN$2fW!R+z^!UEBaY?9J`i)C z(7VbSZ9HAmkQ6T=Usp|5oXyr>AD9*OSRj?D&O;hq*sC#qY}GQegUy@e%S`vD+r#r) zhDeo)0y`;4O#fyqpIl81Z~}F-*cz*^vy+aRn*p(}UmRnlra^lgO~TZpbp6lT^A-0U zdXoXSuJ7h@hU&96i+MRsT%`G=j|JW3_8@X5p(W2Zdr*|@|dQRmkfKVV)H)iFO4gMr5R;h;cu{He{PCSkI?R3LWp2B z+P}8H9|$}eX^}gDj&qxS6nz%{hB&%7K9ygKq&o0+p&CLsK&C*?3=m)mahHq^15@=E zVQ7r9)^Y{?FK>!aISoIB08~p-2_UKwbGK~U1YAEgv5qh1bO*4CYv1XAXll~SjjI`# zjc;|wtJLI%!1&<);HalQRg9~G2i~f+@!RR^MS~0Dckq_;rAS!@brXGF)MiVe!NJcH zxs#vB-uIym6K!yUNbW?sxM!b;Ocg@Its#`BHe zDVv7Z6Tj(zDq645rW>Fnd)gYhR9pOEmI1B}R1}3`6)OrJL0d{6 zr=KGXz-3`3D0zmcMv?8ehuTK#(OfpT@GkJ55y03^0PV?#5LUS(wg&=h{38Fs4ClY# zqh&@rYILy#1ig=A z8++sDeXz}s0k5TdmIQWq@P5{djc^Z%qHtG$OZe`!dIT*+huu?bJkT^b zqRm!q9Cj%4){l^GzPrK|7QqUfR%67_aip?Iy4kAnokO`;dB(s2@Lj7YXY>UQ6PaZl z&y_uL&s*$DCu-<+_&LR1k1f0o-W&N38c+2ybgy5pY3lIwNs_dcF80@x@f-3FQ>Ky9l?x3m#Ke@;0-@e2&`8eW#DD*qpcq zui{0W>`a8`J<#0s((q9?akDBvq)3O2VyPfRN3YM>tXtOVUblyS*x)d$euic zw$ESMenlCOneB{HJW`HD8CML6HmE8nx+W$82fH)ry#ZN zTqV&xqe$Q=emNAa`0}wsBklsDVCEj1U0pIQI zrkK{};CdafJTo{Nf8v-buMfE?fWNa^jrq3wh~DQT`F>jjhs6tdKLgLo!m+pl*l&F` z@Z%p*l?utD?(|ckUrp!m?bYWs=VWN6=14>2aDL~3bZh6y&Ea}v;$mb&l79Pl8vH57 zA?s4$rObtAc9rW6_)4E{JfS^tgnfu@V|LT>|}1X6#OZe`%rM#JuX7(*MjBpkLtf+ z;dL(Y>n&#qQoQVHRT{OC_S*CnsjYF)boMx!hGZ5bT^8}xbHLNX= z;iPehqb5wstp~O_irUn?AF+q5$u=@98rGu>aPCdx+B-M3O4Bj>Unl)JI8QgV&GUn! z>EQ$BUo7SRK!!?lTo|*)qai$?TxsnwZk_m<79St=3o&?Neg7Glq=+AL&@hudLEQQ% zg=!Tk1AJR@wDKf3jQw$+;!Gnc*HfIKbc-74+t?mn=)20gBj7N#kXH*TH$KpFGSvND zuMcuH{&H2jC4$x_k9d?lh^{3CMqMVZ5!mUcV#+iKPRb(VE7N7bH)AF5E9q8G%o=!5 z5C%A8{aIVaE?Ow)iz)GzWain)O+QYwU?6MwP5;p`|Ko}?cQxo1c*Fz%MBb+=bB#|! zuT(i8()u__AjIakSpJaW(!&8)j z!t#mlS(-M};o0~Wfh1Ls@fJa)?+|Uv?^3!8v{J(hN#2^4@Lz-WcBOT$J=<;}q&Af+B?f7NY1l0ymyO}d zGrolYQOG3hWmfR0)>VdFDo5k3z}@U>H6-Uf_ZO51Y8dO-O`zI1OL;jmE5MYrySsDP zelV`IT7m$l@(Ci5OrT5^sfg>^1KcKQ>6#fdq-#{=iDRsL3SYVDigG5#g`j|d%S(_S zFUZx?8g$~qN|tJKdp??_|L&ZV*@wy1iaXf$v`d6^A}(H}thp<*%+?&Hd~?sUMw1?y z;z5S|g5d>uZTsn4ebbcOsi1*w76(i-+cA_xj4^EZ$NM4FPuiaqAr*?71To%feenVR zHaGyX?s1Ms40(y>Z)P=?U70fS-UOGofD%4VtJ^RaYjm0%p%s5^0Ibb{%5(fvtRD{gNDcbbv%K-DBC2@y~I7O~Y;Z~udmwNl_oo$|5zs z6N-M-vj-QdU-I>Rp!Yz|+R6yW*fARfeb-!M_~)^0vKb=IM_WKU^^* z^8!7R#!if#u!dMO)WrIK+#FX4eX*w3&K(Ecgf{l4!XMH@4{C3&A$*=#DVKcdYg~pW zpQL^qzpyz9|8Tti$vZyE+TtCK^?2oU;}dvds_e4u8NECpbgvS5aK{(5tt{-rTL(U& zMiY7PU8%t`=SKk}!b-P(IT86#apZLH$I183;)8AxK)Ja({p#2Ku03H!s?~q^mzP5} z%iUpM{%H|o&xE4a|MP(^X~+GEn{fwo?_G4@F8}%mJ<3{LZ_7o?hRrC;mY0>(u{r}G z;|0UoW)9~9Cp#!JBvBrxkA^Kd^ih>9wO#7@=U{G6TdDT7D?4p1*6?cz4otS{)e1G<-jzU9Yu6*lx7d>2V{@)^nTuJ_{1Km6Ixwu4R5 z^$icn_$eR9Wl$U^Qu^4ei*e+A3%mgQ6%R?*!|OwtI%<#q!+NrpXO!?Et4!tF+ z=7Csnxd3ByU6Ff4eGAFVM*#&R_Wiy{-VB5Z6NjSTiPp{Ns@A--^(TJ%mc*Z(e~ahY zf7^DG7pVZFaWlvZg|P63V9PJ{ok}%>8B1IIA5n7zNo~%lwLp!-YXo^MvRf5ZCsh6e zSbfoVQPV%+YM7;XcU`6RsrBzlF%i7gPgS+wg^igH9Q4~z+zjb0S9Kd^iIkCDzSM-s zJ9g4Bzcp4dQ~zIK^Lk}_jQ)?|05iBBkEJ?92D9Q=iwONXmG?lVrb#-kn#U%JYKxU& z+5~Dc;9vF%uMM2c;4Q(fQ0$^O@RW;SXVYv6X|DqZ2?L&f8^1{#Qa}U-^1@xC>MPYA_XX01{a8Psvl}%j`CrAuIgKPZgZ^}j=K)P!3Gp~!R`xelzhj--4 z-^9;`e6m-)B$1*nQ4SL$bCV$sq4SsczIc7Dn}u5u~|wH+jWG7(%9#+Rj> z##bdiKWN`i zAFAr{VkmEhs3f-rhQcpqUln{llR;zoBr*GTI^wBCd>XEl=mb5Cy@lJ>{qRE ze}UNu-x5r#|DIk?x* zL3|jQJ@H8%FRYykm4C ztNzTVV2@GjbmaRJ=TNzhANP|XpkI6DXtlZRxBXA zt}cU&qS_#X4(1%RJzE*BX8u3&rp^tYz<;H^jDN85LeAt`9!IvriEDamZa_znXxu``gTOw~C54b9kq-uT{ZHbiS&pVZDqMu3&8sR$ht>>I zzKn)`7dXHQJ(hk1jNWh@j-OEPZlKp|9lRQo6EBb(LIg24kiTO~%K8<*!jzjwL4|W` z*G;dR^3@4NFz2YF0bi|w)DgZin3@UZZ$^kyyYg+ivIo0V(Om)b0Zm|+E6JUzwpok6 zX&mx#LnbW@w7qi@-;C}WS=q;-h92VHwbikgfWlREH5W1Osf!=xl;Z9?MQ06?e6(<4Aed^97~Y-Sq!`4qi*AlT zV5{&GX(w9k8X=K?KnqrmI!H=;gVtb)U?fLqFY(`#b_Mx@qqd#_1JOCk1M&mv`-0iB zlvOOfD^-UM#UG?N@|=-64{Hr`WC46PzRR$NyA-GYmdvrzs)s0}Q%|M^oq`QmgJFTX z{hS7WY96)ez*J1zfjWyTc-zW+nR_maYx|vXnQ{Suj2;#c?nhJ4Cf1skcR3^$wKt{J zU{+fTBjN|C&y66ySrL?(gFzMX*0CfCJqoi$xLqYU*RJ#wl;4MsJSs#VLZ=irisX@P zF)ZY-*s!1q*4uFLt#1de=pp$C(k+u6{SsX!l7tls*7b2k_ng+z8=|pkzrl6HC0+BZ zjtsQdaI_xD+W}K)>>h=8RI;mYhAQnj&okkmfSiKXT0fB$lf$o zN{XVj5e?Q^qRFb}QD3qI@R=Wkg&*D98$`M`^3enz2^1s0#dmJOHv`$-gd4pslpeAi zSo^m3F|#hs+}AK_M7p&fUn0DX0SiO>%uT?)iMC;s%+4Z>7<9uHU-n&2S5SH;1{2rS z>}ru!UION6Hv8s~YMk-X8ROz#Y1t;2>3--n;vo^DiX>c{hu8#A@^C5DZ~bXy(*<+G z+1a_aLa6##Z&P!OR7?yrdP$EN{uIP^A}AkBOU>(ztJOs|1@We+)v&9ZeHEU`i~@B| ziz$5AznkMI!1iN1UHl5#XW$hOw{~qI^@OoUw^N+izDrV^_f6)MWEWbQeh($!XXWdI zc5h*=OTJrh%^2;0rcgpj2jKY{ebSG6S_(V`5Jv#)1KcsI&EWm3G6Y+5GfO#jBGWcv zbDE)nIXH-xj|{>6aRk{ayn$M^QGBx{uE*mk>M76ZZo%v`cui_>$}Gs|KmdZg&;Bjjj-8 z|K_eD>PHlnYx7)!z`&M{c;h@N-&AhA(ObR;(Mma^17te|t*`Bbf}2?+>mcc}ejR3z z(k>k((ni-j`c~Km3j+^|X#dqZ%6+5T$wfe|5i`TGn%h9AVZ0LULvg1-*`hDRzpc=w z5=qjfEHDF9hf;ibhgBj+0zF0d@dJGq>21_@0gi&~naFBaE-nFsK+TPPeZUUcn~4Dc z@g6&3dOnEDx2uJGYYo%_mxcvu`o8*Ju%VHRQKA@xvlK5H)?hA`z6F25dMp4#_#Jk} zfieTF`j65n^ipbqG%Ex^vwca>0X)Wb!<-ovPan3VYr_-iKLKB_aCTuku(<*$gicRp z+{`jtPPQjMD0&BZJ+MEa0jI`VWr|((No@|@2xXh#f0a&t^@(2@Z<5;6b)K}SCc~Bn zorklZ7K0lNHP?w^WA(5+MUHt15DLRF2*IQ7rMM1%Q0pj#H$Ti%u8JE^A9v&B{P-^5 zZO%~5bQ<6PcJyr?>6_NE_WuGrpQK9sE2VxJk7R*$Idge=YJW`^%g)-721ah;|89WO}3lMkPCdFx0fTs_OASZlc*lJ%C6 zFd5)+gMa8Bdb57>{;WlnJ1S&rYv!YEuUzZL9r~#Hjd0Q&Vm18v`{j_DTw9wayFP`7 zUGep$LKp-;*>lbuzAJINQ>YK(fb7G<7$(q87X zh4-|-FlI+!UbY6SaFS^b{l8BpKjhRuOiS2*YACJD@eB4@|C(^vlZ_Q`EWKmi|ETj1 zDk;ph4dul{3Sa5Fz!SG=e_8=jZ-ADyZdF555RcK21093C@*Nq0ijt_>C)!6V@lS=@ z(tNUN?w{_9%~=a|sLqT4sup80Y@0f(p+53#c*F3JNCOAP9T5RtD{W)P)7X^;?4BuU zCsqGsDvZjtz;^1hhzfIDlHfLNH|AK6sVqBd zgui32ws=O>jxV;>lX{_BRJHyR%+MQt%gXHH{i;{&SxG;sk4;KFe%^zQ?_d1IKy0?u zw1ncYrrswz*nF1uU9PR25jDKQYDv!!X@y?c)D(^==DiDihk&6bX}?wn*Od1-JvK z`$6#%-3(j?r^PkjYn0_X)ZU`dOlmd6en+P0;(G1k=S?Z(^m$=Gfk9 zy$dxS(|{kS(KvttOtY%@`rtzHs6a)!&0a_lXS8^5KLBF4xW!JR3S~NwvXcIo{+e#U z82@6wiBK+MGz2Qw!28=A9I#$G&zws5uDtXx>n!$mglw@o zYyiqFJaGrM+6u<1;Ef)%dcx?Bqf5zIf^3wyxcA6#xZaME%1sy52ZE9GefDwp#w7;j z&8jVOS89V4KnLb@r-DOO5#|nkO1N?<(gjB*P*35{LuIq_w^~W?)7-b<7EX=kp2sV{ zUkAR(ytC5m0lShQ_TUCVlz*sSnMw?2dOg7Fz_;lK62|QuK^qWdUI?BCj%3T`Fp}|& za9b*%bdIWv8FkiDsL5bYK^naqsWQi{KLk$esT05L^`K^aQ8qO1D=uC zWY)z@Zl?B!J3`d6E1Yo?X`DZbrEc%-myRN+F5gQGUb|zIyGXbLuVa-h@*`bv)h5wk z%y?j)(y=qgv@))VsS2Pr%X!aOkwJs5YA5buQ30f@SZ)_p5Uz0rFD~3$y8near(*TAh9QPIin$olk?90``$>snrja= z{-RvN@(00E1`l4tw;IFt@-X-b>v^U`pbtnqV*~$O^?~FJ#`~DM1TB%A0wP?!1CH7@ z@;2rSE-a1yP<|J4ig8vnfog=(?j*Ie#BYb$ecKle&Sv5IR3W1ab)H|K-1+X2eDOxG z5UpCytOvnL>RlUq!Pt2%w4^N(DVry%O^~ee0b~D5v`plXQf3|7AYGfSX^#uw0L(1(bdZrQ6SnHk+_?ctf_f||MZ zx61u^(p(|)DM7syY%=J7rq<$k*Q~Cvc?|01%DZo&XhJm{YH=1 zmIt|=9#u_@NXGq$i6BNDIz4d_I|z)|bC`eAZ``c@BYYsYF^5_}b0y89pDPnr_Eb*9 zf?b9pi205DNP|6NYRqIBvfg2`B~eBBPnSX9(!1{S`Z!csEzQyVN$<>sPjRxo7O9dXULZmYWUB&=Y$>ecv&Y@ z2leNfH#rre%aG0)W{ux#vj?NDrDOL+#>4hLrrzCnh4*FvA`q`o=`=3#SV?~G)Jgxu z`tZ9}6+gBHmp^LR(NOhPdk&stivfPp-iXaj3bK-I$8rWSh+g z4G|Q$2Lk$LYlSi!Xqg#A~ zPQlfb3kGjhL0KA;&N@PkX8Un+w3TJIQH@^)GM$#<47>{tD!Z7{B1!cG_?9!nMJZ6_ zRdZn8VV;K}%dh}Uo1IIXS@W6_dWb$9+kvEA8{TJr59~Y-M>nPMt`WmHj~2D(_qV9_ zbzKK`nVE@p{0;E?c-<;?h8sQ}TmX#_YWK0^f-A&@z(UIQ_QM|b(h#l^d-4Wvg*xC1 zirYavG&!69z=gzG@1EOc$bRde&*`=oUo3mko@&!iqc#TxK)&=PQ`4L#xbtb!atyQB zdPU=b1=a5zcXW4cDJ#0RwDVt9(rfQF8BF$sx}+nzUkYUw-jb z!k{}&F4&Ymo!?RPw*RGIlHA~Cgyp*(so9aI4D>~&xAXsQPl0uU4@syW_C0a>5EOkw z^LMWcaeDkFdOf1kk;Gxy48?SdGKK#*r@|$Oyn*+DdKYjV#8Z%M1GQGkN@|zipF+<# zcDHW5GFZLB2u>m@ci~P*QH=Dyfld60x(KFEI+T}@Zj)E0cF!vd8Y?P+#(Sl6W;&W? zat1Eg8Iu{ZHzJ4Ul{a_@d*0jRwPk}2`Cau6QuE8Vt0dQX(8a(?A%p0(Jz#VG)GV1R z-S5O2#B!lB(-IO17{f$rO8Quk!0y4y*A9uH=jkO4U=kIp+(q?QFV66e+l*+0!q3?J z;5lZwu`?0cFR2#dZd?0+A@{w);tOh)%E=-Q(qrWJ+0&pS2_!2J_g=T4ZlA!=;lO)s zRe|cpVws_z4X+IrX2;>TC@?MFtfLl=V)baBJ9Z1mwxHf)Y~x}baPdy-VPG)r6D08z zZm5V zEo5Q^CJ3-Z`RbMmjZw~^as8`e*;1imlKyF874bI?N(>?TsRKXzp5rGRwq9(P*i)Ro zy)dWpl7+syzZtSHQO-{ce(y=1B5{4OEZr}3^l&en5(>0}(}6QQzLoyzdxC8iT{*>8 z!0M)*f$%_J7|#=!owd?O+NQj0gwg680+4p7lS^%Xz`>7oo==?{`Kva0Bm=Xl{kW&- z1@{hkO&lT;gBdn9VI8`7n>h(WotA&OxWkTQM^EF60|#9vHd&AWTXuyZa{0)4vdhZz z7x@pdo@kLV2WZEOyELNlj)y}}?3Sk%A?yc)24IOHjrD8qRK+L7bGmQzm63taL5JNS%8x$LnI5mk?=WMUbzm2zrA*K!M`Oz~yXHOu_e8;8m3Xdf#7(v)%awSsM#7L3sw=Z`Fmg zdh@W-#fjlOkQ0j06*i{&{H?E+-vog{M$9=>y7`}q6vAyt0R(U`M8qwIthD{k{!(Q& zNttDK)J)O>QJ*N?g&bnJQ?*T*YPzqSW4$}_LHn2aG~f2PfK#mc!cM^-iTWY(-f%RH zWhBmIg7{sOHd;G%=vGrx?+lhSVce}5Ve_LS;``jt&MIkm;X@&E%OJXPEKY)$9QO4F?NC@^g; z;6YAaF?P_lkIrKvc@HZDM?XShu|uR^cm=9n<80Bn@qU zuv&#|i}?dPjUQrR?}gB(VfLq8E$@i=%xkn4^UK<>I$DE?Y>v7{ieo7M@o3!KGoS-7 zpK-Otg5h!d62mn9G&ng^xGC*oifOJO1@&g&1^J$`mTf0$%{fwaScDf?FjmbVWHLP~ zYP+O#a%=_PL{}bHyJ=Gs z9ymzW8fqRchIl}w>#6A|8V0i*bMs^QL=KEFjMO`#yX`k*-fJe|wn3Qs}NE)|uugH;8s7RKTx@e#gf~HfZh%!>KNb-G0!% z-Qib*(99Ls2X$IWmU_{CajBJ?lMbGFDTmeA1^k@!O7SkDul&M!=TO%Y*~Z-V0v+%#=t+yzd{S`|p{ zpVp=3&7-$A!Rrn46R%i96xrhRo-c1A0?OZZJa5;T=Ls4u=F@%6w|>{=`M zW)0mo5FvH{GcKBagofKzldtFM|6L+F0=q^b9im_Bfz%=>g2h)5FL{e0+i2r-4s8Jt z*FALw@L@LNM&+p`$a^0=-uRa&EJWE^h1_?~MMZHFG7i=xh-DW84^LTdbhw%>6XA># zJ*J(`X^%C&rAhhaXUpZwhYl=OS!WPekyns*EU~qr7sv5d(Tk`gOGc=RUaj%7Zw4by zU%kUxMh@ICpYf^<5@lysy;c-k3weg9IedM-H4vD;f^<}>Y20koPj5d$JQOdm-YaYS zioFZ#sf<2kNp$*2N93wHujx{p^@Z?~fqwgIQ2jj04$HRDNRUVnlJUI9SB*&p0>}ts zSkwq5N=jR)jn>~$tbS;{OW&49cqod6ZB@+xGFi;}zz?@rU9XcarBi=WR91lfQe$zW z>G@8={dU|vCnYBYw3v5MwY01kx;K9YjkXOIyzqUj-YrNopcCI~PP}W-pKDInkCK#Vc!=;x$_QExozQbupdd%S@d z?{%v=!o0!g&^ynxgMx1wNJj?eIcw8<>`sMBbUD0rT8sVAjKmzB^Mqym1?5y>dp8$c zk~jlDs97gJ_pP@q-oZa&Ivsh$dLyyfP^!oeuE{$MyGCYb0I}9X{I7W6_?y6Y*?5Pt z6{Hd%fW~*GjU@d(5vG69tuyl2JV^fO6L*Mh(+SF=IU;`@vF%0G&&wcp+w*HTPu#lM zaDO7ZL{SjyQA(0Lt~ABf+UU43JUf;yc#%+lyZ9A`-ViEE%YgK0>tQV?TqlK==lSJ} z*8Y@-9+a=k7}FlHKKhHD0f~(1tk+mcT@Ol`?E)_wS0S*nY@e*vQ`q2z=r@gC;X`cJ z3t4l}gzwq|uUf+Gg}PloOmc0PTK1Pt9{E^e+f_Z+c;uWS`c`CbByq^?gi6s1=ZoHP z{uBv}WV$s|Z87^kjuj(X+J8k~p*>-ptXyX|{+T958Q|?vG{xp9>U#yhz8VSdH8r1* zEP+$aQf~oFK%)g%UImPHf5F4@ZQ_zrr$hykeV#OFz+}&_=4zrG&bMWmzI0W;h`LUD zf)gn&??QgkTuUxaE1fHHdKJYJRhIqJgB%;5G(Gb&RX*AL1jl(9y1#=VoRo+VnB#BQ z_RN!wW%6PCK4WA{hL#YQAh7dF9NaWHqI+sF@p=t|9|c(^^4^=(h2OY2lve%;`Rv{( z+OWs*_T0_XhIcc^6Pxk=Fz9K0yQ%_rF#4MBk%{}DN?%5N|CYp~f zQilvkeDTO6?4JG=P98pRS}`E1L+&87Xp)DHpEO9A6A_NtE-(s+iiyf z3LT_Qus&vm+Pz2hTsa#IilAZNX)3%ETxG>gj5}l?=nzf(H#HDdhVR*Iu5;v5iX@^+ zw!ELr^r~GcNs!tp>wD9qI!N{^IQ1g=qM*{L5wRQ38~hE#>dJ$@Obxpx_N2BJ54IC#7W)1tZH5cdNk!CVQiI9p^~th0ws)L~M~ zc~0*{4&fLU| zWH%FlRO1W`1SUnY(&(%tVSLn)8avG6lnUNr%q8q)1{o7qi|T-o)4gbsQd?tF+zcGK z)goSISWijLx==KTZ*gUWGn)A^AvH@3Ve(aY5WqrEPVJL?GaL1r_XwNATZ^!$1Gsuy zQ+XaoIg9?YK;l}phM_VSkfx9&w6d(4k%KN+15 zn5~}`!FWa74PJMzZ-m;zxSlPu%RkB*;MM1D;{9l?W;`Q|MvX&+j{~Jw9>C8fWT7ga zVwJ@NJ5;{szpxbMLsaP_W%$WKivIRnUuVm&J`zA&K$^wQrZ)_1V1+ZXDVysNFCqWr zgWQvC_{*;aqd~HHwm*c=2;{^6r>b}HXZrvD|6fI^lvUoP600Pyh*e&4nysQF*xY&Y1ONmOI=0_LCgC~-1y&d znwox?!`ur90Qg4N3(CJlKg40oXHtr|{=bA1KU^USFZx;24<(_t%X22eU1SncWkP`@ zFyr5?J;c3207Zb*+M2ZwdQBs*L)AU7PlweIlKx|z z|4i_F{t7pygbLN+o>)H>lJp0E7_CZFP#eu(i#6IINc;G`PmKq@G`G1~dl$6Cv6-Zb z!iQ!ng#3;9I^$7d9y<1sX%TO}oFb@bwHLuxwn$5mw@cDffvgM_i_6dwt~~{Waz%U6 z=E0V7N{KiNrGa}Z?PhBQhf^$5AkebV3cue(`H*c+ptTEE(H)}ggm#hd>PiKoun7=ddun@CZtH%6t!9YFj}g;ThzrgQ z<|W28OIdgtcAKon^jK*2fcbGIGSVQ2=Q<(2!qkp5*d)gp=4Rq&O;1s+D_hKED89Dj zS~mVFu0^f>s%XP%Hg_}6UVcWWJ_gpgVVx=c>&n+bW-U71S$(Q;HH8BB3?xoy%~YET z7RP}<%6lC*0G1u6gXh3D%mP=Xxr}Qf4?_AK4_)T_IQd}PYi8UN=_b-Pl^^th32csP zNFH({0eBbEXHxHF3SYxi`bETNuXCB|l-)cfYoSFl-{ zr6+&Yk~3xM&O5EwRvm~*`NC<3!t&p}mo*ia&--$m=5Z+e zsVs%HSo)lRunLb`ZL|qS0d8Ib)I`&*BP!n@_NOs4UfGj~^B4qe#b=J-vPkoF+2js*bO>Ql<`Wu@GQ$ zJqOlTcO(B1Uzp6107NvTo!D{id`XO|h;@jZi?YEGA--Svr@?(>Q*Z3sRR1q`Qvp38 znVeG?sB>2^MuVk4Ejf+TME6$1P(@6qNR=2;mTITtY04^Z+_Dcdm+uO^D6x!@NixMs zDYC@|Bl>{^obRQ@%fkj=?W8eGdSs5XQ6E?kq-?CGr7 zIwcEUGeZC)26UwY{~^rju*KfZGqPRsDJ<^lI+bE9e5W2oSb2(*4W`Fd-isQ+3)BJ* zhW>s>T_945WapT#JZ1mNq;}Q^O{>K?9{A)8tArGl>VKZf|I^C2$XWa+IGC#y7C*ZzD&K=H~$96zPjK=ja*2+!ku%Rb7klrqPiyJS~qng&zk_Car zHIrm^rqo|L9{6>>YHcWR8VV77JyRaG@EYeMSz^YI`A0~*L@=>mcTxC+sUYQ;ov79A zyF~^O6++z9sxR`9E?OIin3#3m59s<2?RB!i>NRem3Fb^NY?d6Bi619wvO^gCf)t8S z%q)~Rx8cIta}cjQvAyD6UmlpZAo1gd1AYh0m;00?;evB4*63n!x&+C-hi+XGU7??m zT^qjjTSr6#Ee#i5cp(CGASxuWm6Y|CpHiNS#o7qDvnziZNDq%b45$AB$y%l?Vdhbj zb`lYYM3Y1hyjXZ6<0mXGMf&HWrTvyEkrST*upjhg#^p)aB5eS-)ESevidgvQ#>Iyc z`E}f=$_`m*VhnI`ZL&i=tduO1X#I(IGV9W*`>|5Y|KbI<1V?LXx{avPjWtE| zdZ5NmD4$e)6wA29zgS~5mZcsT_-5pDl0_ehS-1(99P}z@R-ZrEBPK9o zibM(g!y<12E9Osv)0t}=@&3<%%gCAi!OIdW=tSH{$MXXqxiwF%<(RY|VuGBsf(GFJ zfgFd6r5t6?{XmJRh1eudmX67$f$;6|Cg0okCnydz3xCigX82tMPg#PbX}t&Fx+ zv)(r(A1n_xnS;;EI4WZ~N)aDX>EF@-d|LNmw2{PW{v_~UIc!Bf_5Oe2I!8g)tg~!~ zZssqn^QzO#TG1P%*>KKPi=DuT_uuNQIZxF%+iPC2rKUXn2_GScG#PVOq7<10sjsY9 zIt-6ar{(TdlOtFD%)a6()^mAGbH#M|q0GVh=pSkQ^DLdWgYwHrcye^}O1prGak2Ci zP{%>6?H1n=@32n9JfQ5m;QB2M#uZN^tW}NiA7R}I6a8YbU7iEn35@xJ2An8jY*=Dk z26qF=j|@0FW)*hxR3oD+o1!CfZ2R|%t=8A``}8d>)LgW5o)!KO)80y(BBd z|BniY-2~(8sznjpW@;xGe|*X@*l~;H%d1ygkP(N3q2xx=Kcfu@zF*{JMyw_SYIgn) z4{(Ag&t8tmTY{wW(fBNWi1L~w7ii3`VU!W)s8gHx`jS^n^#K5+Nok}W&$ETjQ+1XY zUBbEqOuc{5Q7Hs*(*2pjvi2Bfw=`R;6JF>qM?ag_t4-=ep2PDuqviOC;+nUSCOuD@ zokbwD9cCvc3mB4*JR{#-0%u#7w5nIw{OViMC&?7&xI!zXG*K zzXUYdFb}XaEb0*8{lSEI3VWc_wIiY1ZPXv?--sCuB8Vv-W= z%g@eeH&Z!)X!bzYQ%ZSNnW}Q}GMBA!hq9e6asN7A(<|5}tNSx$9|=vx5lt9?)8RoH zER??RF^RJs`irX(@=|lz+!VLMoG|8}AZsCsioo&cOlCT?D)ZOYNJRt+e>hsLrgz=p zV{Ay~9HuH`ayBv-_IyGovWfkH3ga9rVRuIk|LH*Xf1;=dOUZovL` z;wCn#{n{l5M4A5Kf{PVB=Aw@DnevsuDeED-QCp$slsNk{ClC;&e+aMf-~LAdzf480 zf}b;^p*5R;;$2cd;O5xEbpf zH~FPw4SigIU~cscPgUUn78_Sb&Wk&7Fwq?gKm&XEYJ^-DVE)bn@?%SXZHnw)pZZ|1 zLMj5aXAGdu$!()m+BH2@PcleHIoeqMQ4_yf#Q-y%we^{ofmoH_pGuLkz4Qm8WttX% zOHl}Ak1Pbyw-B$6W90L{&S#jYI~z;L1M5U#(1IihUhJldiCW%46%W$ zS)(f9Bd=FxFipv<%QHAl9Sq){dbfTD`Ln?q05%l+NRRv%kh3Ym3k_#&%oJ25_Q^ww zIbEI0UjGZ}cXmF!K;ekQXhBB{+)P_WMkK6x?iZeyTd-q28Tm6rk@|pj9QP(Qv8L6# zW*yr-XFc(g`AjufM{^{WI*X z0=Z*xgv)zIo0)&dE&wiT@DDMS;*enx@ueA+WkAB8!m<6*H$aX8G(4qRUABhQpNylZ z`J*;sCe0y$d%x9$)$+C{cJ3SsG@r1Zp~|)vqN+XElXS}d8y!bfV6o|^+0&jRW+p!S4D^@x|Sk8V>I*0O_< zwnZoZPzLJ?mq=(k2rm}`9gf3888w+iPvf@lu%3pW9I|5#=O#5iztz_3jl5hzp;mv zvF~JpU*>9XL(Zl4g(X!Gs1-&+$ zORAX?iF1hLZ@jixXd*7Wzd&hqq|i+|PsNE#B{A(!MhCFm^8GhS;&mkU_%?8H+Df-> zIl^3D_A5(g2Tl-e+Zf;`!wBc;c3qT{wc$mdV+WYfQru@B3P`g97rGL_&M4=;5h%$j z#(KR-{wcEYhQ!U}elkS17E?qY8&B6*q1FzfhJ5dj0o<EM* ztRoZXD&4Q`WB?j9sX$Rb#+lcKb26Iy$_>z4-4zRn631u&Iw%$s*;K%Rng`P1>idY1 zh=0??!?cmpVYRVP75^5XY|zh~GV@=s%Z@l?^V|JZKF zYtD=wWuD-SRg9kLe7dgn;J)Pnf&M%AECg-`Q2yU1Bu^5~lFkqyv~!MFck;AfY z_UYRbkF<6x?)p#c&O?ej-|XBnyD&Ptu+?QtI%D5|T%1lmEBm{wD*EBWM`h2B#+N-) z{`m4xPM=zJcs0E!oNda?j4DcX4OOk(6w7%hE>l^#TvG+*zn9iIrYm*t zK)-EHc8Ua49x$NM??}Dg~L`U%-dVUK{ugOy2Y|xwt*T)^j_IW#3*6EpS1gdn6=kqZkM4?xx z?J8};J?OZZ&^*(BeAj3>>TRSlG>+7$)CMqlZ}12ZvLPczIV3vLd)Z(46Y`R-*!|A=yYNWq)2y{x)vgrW^|6dj&7>4LP879s1tb2HTZ$9*6M^b#e}Zp{Uu2 zkBthJ_HO2Cm5Gk9+K+}`7g)nfc&kb?%)Wo(Cy?Oym1ZmMj&IX#bpKG z^-Y`|zvA8$&C66mHL0VUAQda) z<+r|aM_O~Zz3(%vP{2qXA_!0aBtGJq?bZhfibiZ?vy3G>csntQ6h?oV^QD8{up7&G zEF-`E@HcHDe=lc=W%_N`=tU?ar?c2mzTr)1l^OTAAvCquO ziM$*)PQ2)abog62Wa5&>MeAo;=qHrE9V-v9uLJpTix_Jd<^f7S(i-v$)taJ=NOgR# z&e%1w+c(I@^B{CV_)t!`tv-f4m5;Zg>$8vI`s`OXK(dWl+d|t{T*-LqBM{#!Xa|;| z4%wfP1>yrMNjq!H=&jIfg+2uxP=)w#MNlt!O%~dj+mB)ANzW>5MWFvYPeA#qB zsm^giDgO4-EmSl;XhylgP`qu~hx|-VKl_1vsI+oaToG{?!8*=U9IJP04_u~kuFd6l z0z)@p8}-CbV%+M}%%M4%<6d>;=EQ8zZ>BVM`j;9+ ziPEO_jUl)WLIs6(rh$~skwrAT0Nmddl9XI$&Rl^ zZd%>ro;57DOE+q*?f96!1Xj$fg9kkZSmKuz7(2)()t*vVC#8yF*Rg7zXGqiB<+@5& z;C9dn=@k1~_>_^^u~2zoBGna~D0YRs0p*k4h+_rzk*ob{Q1ezkex$832n<%4FOKgZ7FFFM0&E$*ItMb#yC4bn$y6%F$}2ce&^ zn=%`fBOXC6QEp0hqraH#{Ehzt*(b5|m0xJtYgyL-ibT%#Q7c?s?Qnh^J_M~rwH#Yc zk-K^8)Pqh``>#WPnfrROOW_fl#-14qk}(_dVfIk@|8a|q#lKVkeHhr5@pFzco!GPo za}3Y7s;?m665%1+-OVu17%A9J-u)%P-xHexa``)AFM1~@OZ|HSHG3FtV=r{FNr zvAH7dNX1$l3ibjNq+{|(ey5d0Oe}Sn=gKXp;t#!mEr{EYz+gPw+Bz5d8GVGxsT(Cd zA@*GrABgcn`W-6j3S*!R)SHSo)p_tfX#^WgQJ2C8(QX&0p)x0|}Dd^|gk~U2^VxgjME< zAfX)3n6RDuUQRp_Y^}!FW%d_k{8h8FzAQTP5{{Lj8qMtJOaq!TLq@8PLX<~aC=VT1 zUsM~2$_L;<8m2m4qMx$hrSq#hYB$oR03)Yj8Q5%6JZ!w36W15I0@JL2!fH(7eH=g? zU+uO?xuE$=D@3Ix!Fnza8E-WsZ-?Bua#Zvp8OTX9HU73Kg_BlO#iLJXzt7*R#L|uU z3Qne0iS}}K05<>$5Z9pDE^TkfJ)Rn*8hMqk#oTZYb2fauC<@lqQu$TtS+(Sfbe%hl z%A@>6XJwmL2CX@7#HR3v#L+q-PAmS)Hk8SZxC^v8bt+qE#-E(q3$AniCfbl}{3bA* z_|ohUuXEs(ukL2|GBaY<5W6k%n*7GNs%UObR^)^S?^sI+!P9b+5|d zSLbW~15q5ws$^L}U5wF34Vs+g{CzPVY|7jhS6=+$TeR`O1&lT<0C?HGUMJ!ZjFr(G zRbukca4@{ONA{h$&%d*?z?Y?8eFZ%#ZL@$X%SKHfHpdpz2;b{mxy;JZ6MNC?8B*Zt9289flpcg1tt+6gR5t@eWgISU!SWe< z6LnRAp+J`dn_I^>R3EOEnS;DBoNiOlu2hhl?#`WScL zN{8NK0}@YB*Si`!4eDL?wnTP;;e~;e-2us|JFr+|7&VcqjUKeQ+`Ow>N zrK2n3oEh`cR<|-S%46S6$Sj5ii;E49uIkI~Hz9mMl= z!G`>fA)ryYFuc*Q`3Y!luih2B52%9Zovxx~R|{ZE*#Wa( z$gZAI#)Tq~FC;&);^&sMcR3i;9Gp#A&zWjY_g{%j{(F&r(s;eESJ%gYiM5#{xL4Ih7hV{LxW;Xb`ZkPN$*A(}%!?t`~1A&FfT0d`sjuHt3??z@DeqM zU@fgG{0p?yaQK~{%njZ~AGh`UhSHgX)3y+6vkDhWv!!MhKaP^QBQAskKA!Rx3g2;n z;L_|pvMh*pj&Q3Z1ArsssebvXj=dYxHQ>h(`0i!BM9%fD$hLR8%oL++Z&)nGPs*jcpsfTIIbv-=cTDZoLmGdr6xzjV<_Ek?M1u0^#LS|jLRE{R+VeyQ!H=NWSs4F7r-KXs73m$zJB=ZeFEOw z%SAeS^hqVU9~(P~vmK7dd>ZXEay5g99rlIRV$x?b|JMi9`)_8U$LNwSkKHW^6v%t#Dgya@5h(Cr?DvF?%Pl3ov=jbPu;Z^*-$dXy@I%Dy}& z*K!F%={G&oDn3?rjTF2CxvW2gk&4HeFKG&>lRoSI_`PlLIiyb`TnNZP)#@#pz z5=348blHXyQO~Wtec9K%!Rr;k*XPOtnZ8xp6URE|76$T0}&Z&nUwrC95zxJlE#L)MZ%gIBQ zQR_v04ExEW@eR0nFY*3k#zrfL2?bXrA_Iv>NRk-?@o}#^uUmOHC?a$EYRn8<SQ_t zcSEwuMSwAUiSu_A6m<-m4gV^Aq_>jznd;7O z=PESQ>dja2%wqxMi-okpGck{i$-RC}ugdY4GV2wVH?s$Aoj9+;#YcpjVMzZFU+)Qp z29Xo#z^FLsfULanL%zm6U6SLiPK-6##ejFjE3R;LeOKUCK$bT^29#X%nU+*lO9ru_$ z)z7ah3^fbjTSxj}z_eelVEkl(n@!7qqT>5Fr-w53g$cj0)aNgckpg&+x2sI`3LS%yv3G{Zm*<)EN2)l^m~3?ZOx)6WzUM6@bXhxVhvx(sT!&YGd+CWV@`~hnSNs zgv~7v9=2+saYc$*I`XKgd z#68R&Vot#AQhUJYI}=awN$IP|i{fr=K?Y9_^Ot00P5>O%()z5;?NZ2F4la|$qif-- zUX1H#4fH26VzoWTfe!Aq%Xcc;8!Smy4nRRx!N@13$%9}U(X~QOZc4K9ShvD#`aR&Z z`rUYSVTQYVQiPN zTP#9}f@IjCR=;t2GoOsQmXsiai=R{I2&u#uwNbh!{6#Bne@o&VqXT8 z&s+If6PzZCXqr)1O7()rkzc8af+r;)2etd|y^pStD4p^>V4ptTfixQC6A%0C_Y^DY zu-y$1(x!HEhgv?RO{AnkACg@Boqq3_7AfHHMp?#{qsJnf4*sR_*3>TUu4CGevv|im zoL0#F!5rGmcfWh_o{I}Yeid8R>qzILVWoq}Nu|-plsCpj%nszoUTLr2#YVn)6QV~c z`^8|Q6KU_O2ibMJEyEde%n$T~ck@{g<`Q543+=eJL zP+W(54lfJcO!yFg3UQ0vB5eE9Bv=7X<}VtytHWx~q4$WBP!Fd!RVSWtvLg+hzBqG> z1D{JPp`%p@yRg?a8x#VQOf&XAkb??5Hu-X1M81UkN;=f4NTfW^H7s9PH_|n*bwH>3 z#^hyk9Z<}O>gea>pJo?d8>M4=LF;yjTCipAWp>wac8n6vRj#KEvw_^*s_x6L8;1hT64#whvwg! zB84f<8;Caa0MVOlUG8Aa=#h|fdedoX*LzcxisZlCLJqR)9hWkomrQn_*|5u}VFHo7 zIfO1;Y1oY{oLSe!O=Ee;2;10W+)9sB%yHIYOP)*_8f8Vn&02U5ztIb~W zCfU%eAMcQ?egJ_e&f;}y=6#X5cV|SBO&8CcM>vDtxI1?vkh2#8iEf`LPrdwO{yIFsSKeM6vH4fno_=dZW`ATLdEYcAX=Mk{T0uzxPT zO8e@{0ZOb22fbRwiLX+{C`KlZHln}0Z}u`nroVf~anNehLGI=x_YR=oZ9 z6<3nD3STvRV>}P2Gcxn=$KIq733p}!1xL=!T^lzwGA$EW><%1=rt6xN6b^XMv>*Bp zQ^#ZE8}E>SQWxy~l>XECH^$&NDmv@m0-M05rJm*;od^md%0cKjbL)}P3;!8(56OiN zW@@V@uBU97ZAS%!FEz|J*cV2Les!FMr1hFe_YhyoI*0tiWKIx!&T-^+bm^i6i)B``#z^zRGSl%*(SjmEB_oQTNly(|dHL5sAAR6y~b?Jlu#Z2X|z?N$C{3KF}5S6eQmqUl}-GF*PY{*B%4w zQ0uX;7f3@@XT-j^hp0?W=Uiq-y}s;Qpyv})R}F~c>Y?yAevDH^R%&I=x;%zg`Bc^I zrJKxD{wd+ULFi}631K5y*UZtOqpj6Lp5;_1LH7rzkfv*RXQV_Y3| z`#)+TXtia?fm;2-eqc@Bh}A&DQ6$a}>LZF*X> z0$Kaf8$g2kIRS3R%KXIFnYIBVVVb#PtBvj@X`V$3vNZR@>l{pV@o+EVm>ijZ$!s88Z|3vRQ)ja2)d?z8|xSPB=Q2# z(gMk3_oWmc;`YV!(mx=w;O#HY1nSQV9dAIpThzqa_$He*qx-TpaD3HqOmaun@4zyo zNAS1T^lyrN>bpHr5MxC61z0UVp_qwdxX;_kSVj$NDV?PEMJgxK^vteHKV}A3w zyw0N0#$3&}2T#*yOX+dALzIrV_R4~668&8fvZ!CnFf&yjsW^zgMQN>fE47}c{}YnW zdj1`MSm@8w2f;gQ?vMUkpt?Z6%^U#JXKGwUlVBax1ep&kAcns%sv~mpdMd9!WXj=>1&S_pue7PKueJToK-9yyiQNel#77E#5CXSN?Sa z#a|#VX-B7Y++Q$T9-m164gXA5w#Y9ni_hrQC9m>vipEVg5nfPvi{rvb4AWAl+iOPL zgCD|#hQPZF$BhT6Knb-%g_POc>UpHto`eo!ZQpUPYqS~1K zY!#xgDl6yO{9w0_qbOr3?r8T~_6#%=RSd9vlA}&NEa> zva#Ko#0~b2a!Mp&O^<|{RQ6mML1)(*jNKUX8(T*X;`Q%*$Rx|?sOl2~=$=^@N?%fs z_FUY&=80e4q<0mWz5S<|$G2To%A39C;d*PW?=BsohPqm-ul6(KoG+5WdX`yTTxuFN z7g5~d!conw+{b=Wy*}S@uCTY8>k+#4cRw#0k`gkwp`5YIj9Fc8VJVdOvEsK;c2Aj&G45ta~oa3|>^pMir+lSe^{@?s81czqWYEE$hXRj+j6h z^P7Q+G0_0?ImY0*qZ`Rt;*Bf6K{VYN)@YEX&~dl)>HWJkmqR{&%o`gGncd-W@YY(5 z{kxAyy^go4%7#rupVoo}5ZA4(W*hLV22=3-{1q7gn|kQru+Lm;ChZAuRbSBZxI`o- zp7wobe6^F`F7ceDUsq#pdtuwJ>;V~8{TTH*4=fvPI~@!9c-pcz3eaC}{cgrM#OxSK z3s&shW%RomgvnV9+!j{o*3LEu@ejdSTTHCj&!TRo&i3+VU9!EV5>T4FoEtG!{3)U#BOzaQC@jIjhV(^xQLnc&K#oLx8t> zqATY@!W2~Y>1TB~Gd*(cIhS6rkZtDq$g{w-L>D{LplYZzQtZzVVApd#83*93B zO2{MQncEjog}y&Pfr@L>gK`b!r)E)r?)BvozN{zD$NF;c(>}FKN>KV`cl$7Vv{(0L zZa}Z2dPx1zisEKxX#2s}q20irHOY{+IjBUNsl2i`HUBgO2I#o^QK?UvDK6inRn4A_QA13YK@WS*1xK;gg`#9_@c`fUlAd#ZrK*Y7HFfnJu(h#&DaX4X3N zm%Y|Zk+xT^gw^kCtK{1aJ2!^Ugm;KoDYl@(iKWG`*^KRusL$P_56JP{GZz+}7W}9s z4@>9nh1--&orndb1&B@%EAx z_5=2}*$x2byEegYp)b1V1MrqbM|3%$=u})~c#7HeuMVU8PUlk=6!x?f*qzhz`AOcx zyRqbCv^2AIF?3M>nD^|%@o0&edf=QXGi!NA)#sV%b=;UXS{Htm7C`@cVd8p^JxchN zy@%i{-9Y997Jm{0B4 zDJN03eDi{r@ZKFxFw~(wKX7orDq-h-jb?9Qb4#m|LGPVA=)Pj>HxqnhFWnO%NQl(( zn$-DCSEmDkMiz&4eb8?(Ic<`J-eYjP0}xY}(v#rhJ{S8su%EY{a?h$~d`xAfC0-At z6eJk>{8sD4erIsnjBGOw><7_`QGu+!#}TMPA4$%H4)wrJLqo_JgX5e-y%Tvqdk65i zLq11{fhYBKPhZ^F$Ey4x3A6K*t sum(abs2, x), gm4.u), 196.8695297987013, atol=0.1) #@test isapprox(sum(gm4.resp.devresid), 220.92685781326136, atol=0.1) end + +@testset "lexdec" begin + lexdec = load(joinpath(dirname(pathof(MixedModels)), "..", "test", "lexdec.rds")) + gmldγ = fit(MixedModel, @formula(RT_raw ~ 1 + Class * NativeLanguage + (1|Subject) + (1|Word)), + lexdec, Gamma(), IdentityLink(), fast=false) + @test all(gmldγ.optsum.initial_step .> 0) + @test first(gmldγ.θ) == 0 + + # don't know how much sense this model makes, but it's a boundary fit on two + # boundaries and it uses InverseGaussian() + gmldig = fit(MixedModel, @formula(RT_raw ~ 1 + Class * NativeLanguage + (1|Subject) + (1|Word)), + lexdec, InverseGaussian(), IdentityLink(), fast=false) + @test all(gmldγ.optsum.initial_step .> 0) + @test all(gmldig.θ .== 0) +end