From 75c5ba448a0a118536782195c4489a0867fe6702 Mon Sep 17 00:00:00 2001 From: Nicholas Clark Date: Wed, 3 Apr 2024 13:37:51 +1000 Subject: [PATCH] small bug fix in Xp generation for trend_formulae --- R/get_linear_predictors.R | 6 +++--- src/mvgam.dll | Bin 1084416 -> 1084416 bytes tests/testthat/Rplots.pdf | Bin 22921 -> 22897 bytes 3 files changed, 3 insertions(+), 3 deletions(-) diff --git a/R/get_linear_predictors.R b/R/get_linear_predictors.R index b8942b72..ba7e664a 100644 --- a/R/get_linear_predictors.R +++ b/R/get_linear_predictors.R @@ -82,9 +82,9 @@ trend_Xp_matrix = function(newdata, trend_map, series = 'all', trend_test <- newdata trend_indicators <- vector(length = length(trend_test$time)) - for(i in 1:length(trend_test$time)){ - trend_indicators[i] <- trend_map$trend[which(trend_map$series == - trend_test$series[i])] + for(i in 1:length(trend_test$series)){ + trend_indicators[i] <- trend_map$trend[which(as.character(trend_map$series) == + as.character(trend_test$series[i]))] } trend_indicators <- as.factor(paste0('trend', trend_indicators)) trend_test$series <- trend_indicators diff --git a/src/mvgam.dll b/src/mvgam.dll index 643fb1d5ab776668fce8c1897384fbae39529ff2..d049f54ddabc829ae1cfad03ed7510027f459463 100644 GIT binary patch delta 80 zcmV~$yA4240EW?Uy2k+n_THZB;J@>X j@Psdn2!s`(h(s(Bk%~;@q7bF1ZtA8rTW_DJtH1FFv0ozS diff --git a/tests/testthat/Rplots.pdf b/tests/testthat/Rplots.pdf index 0925395a869c8c554aba871eca8ff768a0d5e25c..7a644c0616a9b5b694fbc0e0763184e12732f585 100644 GIT binary patch delta 9605 zcmZu%c|4TuzqUj#l2(#6Ovs*?k*!@N5g|(?Th_6UWq7C;mF&w@LnLd)I@TCVmLXe| zWh`T|Gnlb6w)0T$d(Qct;}0L7XWY+yU*GTby{_x_@d``7SVSMC2)QIHFDEa5NmfCe zRYnVLb;H8LT2So9<%_ZxFUkJ+M691=Of z_V=T|kF<#)j-3?QRnyPxD9jl_A7^>WWHeqaj%t+#^eg7Z0=|*^@|S4C zVYs{1le|-s-mypnI?I3!eBUB%s%(i?va?=UHp*=Vta2*E>k=cESLrtVUw~wIi+ZG z3!B}@5aF4ahnV_|y5gg&3M6~~oxw;ra1A>LXp!BLDO10e_NZYsoAqUT+k^0cErGqY zElTWPVkZIbozb$9fIaOA)ZPY-)@RnZH*gKTgAZ%}R-UR2#Qr_XlyXB=y&4dP3rU0j@)$NA4v%FY{yNA~!k0IMLaq z8!>Yn(sk+=N9=vzM&)bCf=F`;gsp5D&j#fImX+Do$#^ydnnRu~*m*V(8L1q)FnB$|k7yJ^2vxw4 z_fr*&KkA2QlQ468%bAF0!CfcoT?YjPf%P|CvKlH#kJ;0oTustRBc^hS@CSqKU9$64 zP&%YDgxatSFK*yV8g)2{&p{IP@F*;~m>9&bV*T;El1jWv96y?m?Q zW{_$3Wf$!3rTiqhv)CKUgDqVJ_EA97vkj5pT?@J#BdXZm>n_;C5x!3r$dTEqRYu^` zgGrCb>6M3j`fzO-xwBfE5?oXl6xrj5#;-Ywl!1o4-d_P@voK4fZ{%XZuPUA0tnU2O z-Ga_M0UI7WL-J>3pln-XaKZ;GTlvTwd*{V&rdoAPt$jvy`?vZHtZMrTAKn?B^J3yA za49FFY}lMQFGuW39;WSyRn{csxdV|JWj)6?Evsl2J#SP$v}a5xtb3PTw%m17duTE) z+>WxDXxx6e`O`kUr*+f7$%ofYO%Yu)F|=G#(RFzpxlsDDORZzQ7vt}qe;9tgm$5I5 z?6-pZ@Sgjf*TU+geTqt+hK$S--YwV0o3rBHM;Sxq-HI}WXS5xxOiWzU)_GmOM zx9-wZnf1M?k2Nh&6Tk%c2K8Gf4nq5GIS7t8_R^ZG3C!t3Ezl8R11u7|c@1-P(}t&5 z2?4K=`8th-;;N%E;fj#RaWvX!>f+k7&9n>dkv0;=M|tP%b9SZyKgII(AkQdxDdFXA zEWETU2pR(~wI}hxcygfFd4jQ1$*}yU4cdR21N^SlFVPI`3?oKX6FBWQd}%hz&cFO- z(3GpT$}s+=;X>mueV)}u`#eR}r=CsrG6ul3HpMX~SCnGfJ`k*e!#dIIixhi8&Q*!y zu5>e4!T z9pj2kJpJp@-O82FyOqPe6>8&2!*AZ;_;c{A3O-W-A?wkIQVZA*vmt~-LGn^ZFfn-T zYnYw?V5}GCfT*u6xzr0nQEA9Tc=_qZAl9(y%09MX%%6}0*U2N)2n(6VQgfVk+WrAg z5Q3DUGlR|iBd^=1liI64?OAnFT7XC4YrB)%aq4g1oF7lh-K7#&Iwm`Ay=KDs+b(q! z90n`9OqBPV(wkQ*SCjVLqK5X&(R&SgpeEb8{5>*yn+-}4rE2_Rp;hjDlc-83lTr9m z$0SFaxfI0_t#uDDu#M$SD;7_R7A;xD5$@$K;#^)Q=B#?nXG0_EaK~_fyQMpbm!}Qv z)yj9@@DQJ6vB*z7m(KEPjfAB)SM)Y@jXM_Hj*oWpwq3r}|K^!ue5@$_N{1TF@O`*q zwlHTp!pzUmcgh4?Ar05*w$JL#ErC$xdcgPjAMY$4P*J=bi{O6~KAzNGxhCObi+*f- zhq(FBO51@{Sy7A%1xmgZVs}o`hFf`tlnh*iRbopnJNs%H%gg&(gxk#{;U*x#Y(f<> z5MJ>3r&Gr8rf<4Iemb4+@_xT$?p-m25qA`+Uoc4R1yM0t!d4r^rjNZn`Z z4q2EoJaSVC zZ6$7!oVCcqZM{hdGQU)hl}zzgG(u7SY@rMeJ8uR8)2exQpEji28(oDxz4bDQ*xj8r{M&0dFH^{(3jlFXKfoAual0GlR*e;`P(fJ z;09p1(^G^f-3<6So{vPt;g(_L$|*J2T`7Wq@s>sMjC2Win!Kt{j7x4~P{($+a&kL( z2;>IT$Ci8^E6}P8btXPa(s~XP%}}&a{Zy9i)N_Zy#4@~U`Oej7ge@F6^E;-H> zcPw=TEFjYerrt3^6cz3*2ZfCmORbznL-WghYdx;Kb|t5q7>vWKQFr}jy16S`ln3q0 zyX#8*YGUqp525*#yw-qzoAVBC#pQ@7P{+_56eQmHM(?-xN}5GBmuI3&mSgLpF~3sQ zQHZtJR$j6<)QV#269WNhMn0~oaNir`n1c`n`S?xnbW(0LwP~e;tMk@)l1zs=C<3Yl z?jU&6jN_SQ0sf!Z}t$F2l6+c6i&Ufv?X`;tGcNXHj-vAO(?sFq#Z|TujIsh8|^$K z=Qv=?qU5FMumdCHnQ|<39L7cW>RJ06!xzH4jSRvy+vbyo`NxM#DZ6)#DFLgA2y*pS zYo(8gx;!XoKZOh3$x zpe1ujgUeQfQ+>@u!T?lv&ava{WN3#6)$CQ?+T_)=s%5Uf!YuB|_kCK-7nXF1qiP)I z>|>qj?(raNa|PHl-i%j{Ec+?K<}wDmsB5cz{RjcZ~1? zRlK8dy-_q*O=1*59dfuopek1S`jB289{e2^Bobcg(J=E32 zrBV?g!2~v2x(tYA!ZBr8DRK4V72Pq|r(c8J?$9P|y?({gze*@JSx88zLT`H{8hDtB z#9AcGwejrg;VJ8v`DbTfodC+?(DzG5!Sq4*j4cBpu~R!}t)~*KE=T$Dg-B6%RKGB! zoj9E>)=7Q4r$>F8{OhE8z36UaooKa-wE25#GW4A~{NkQ~)5V&*Jm}}c z`n%+SJ3?GbyP_E;QbE<$GTh?f^4KPB%akde*WOBSuFHLfnLOE)-q-xb+le zx=b}ne@?}GKFQL}?8X5+uJPFlN3l6jyiz=~@()Qyv-Po8_wg@niKLv&3u(mp9}&@D zNEebg8FF|frry9f{rMGxUT#zSlWseqy>CK~ zq?f}6Vt85NUSsaT{P`brczpSCmz}}r(|xyWt!C@T*Zr$fXI$rzz%?M{nVuP$-8N$=p-jw5yL3(^BoWgR%;FCj3=32tU_~2zA>M6 zvKSuM%a8J(ZZCb<;MCleDs`9sJNK{KEA3V0t?AEM=kK>67K4apMQJ)-<2I}jiYUoE zIYUB<5R`mA+4q4WxxMs~N`~vyc+u6(QxaWUQgYefXy zRk%#|=wbicrHn#FbV4bqb6FcVP?^r@@C zUj3_*YTHdQkl-R=s2=Nwusdrm8aVLK6-Lm_pJ7Et@Dj)eu@B|7v}W=0+SOycw{Vxe zogODVzviz99CBe9FQanC?Jq{89l6peLPqJ~FDPd<^LfT<*P#Z7AkbRGzBej6P-f;=cSTAV#cvk2yGZVgr%gZ^`9wE^s#>7=&&dxFE;-T09~dxxW}t6G95x%Qs#MM;Xs@m0Lj-M;JF!oXYoKQWh1 zT+l7*t`W1QSxHCeO6h#{GbYLEWOY7xxjxCUEo`JKrtBtf4Q#CbjG14QerLVjqy72! zu@S|Ma483sM+;0&M&h=-8mIdkPwJ$3o<20oGwxiFGm7Xp{R26$PMR>8@HC$Nf*aO( zDVxiBpIj`1UL=+-jr;EwKJw;IdUaU@@045-OG#C`VKd5=a`jJ>Ly}V@jS@!y%s$Ki zt(UIxqc8Q^aq4b&T_b;$wHAcQCbR9bUMN#NRiR`I81d618C~M+Z#P< zC~n82YK%18zDd3Jn^b`k7NP2(CH5`>REpOi9A zsf?Qzop!%%-KT(m2;ICXacF-AteZp?LIjNc$4{xG*Y*htt#I#@d;nO2uKA0~Mul-4 zd1H0+<6Lul+`=_Y&d9$F2wGdc2PHbF(EWyXy3dVjKIx`M)9Ie((7_aDC4!_6la-HS zWz=N2`mab$NeU;TZ8&JWs2C-I_opo&(4e4eIGBBf^=dZ?4`b)6V80t_MS)dCI90J z{IR#3r)Dd8=00!mDH{F}oO@WS`b#om#}39;26t`kBgkTzW963b3a8A)y#A#M~15)FDj~9xZG*Le9Xq-Wsm_si|rWAz_i_t}QX2i%1eoQJDaQ-09>+2$kIUx4*h$ z!Z!^ZSxn6$r0527L-~GN z!)ofjHQ2rEM2?@KR%;Z4&8grF{_jCQ>Us%mn4F6e3J@% zbAG*U2x_~`w|zRA^_Hc<@NsUmcuq5oTawI|)GQ*&BHc0^gSgF?JRl`KPgLmAwVD-` zCz@E79kA=Sz*;xm`s~xN7_H|!6JNgMUf&8$herd(xgw$0{uswHMX zURRustf{E>>0=T3Lp0_eQ*gUALV zhQiD2FHH#A+FE^b#Y&bBUIB#RHx%M7*#qklnp*TTmrjH`wmC zoZ0c@oU$i#>oX1U-7AiggrKRHWV_|`XByz@zmmaWESI15EN4==m5@+K& z3AJm{*ZjE@HAhGKCGQwKG&e91_dbiW$}zXk|6d~MCfs}IjQZ2@cm|Ph#v8~ETK&Ls zdL|QSgfHE>1mvhO5Pi=9oG$g}(-!^+9h3iINM!?Gud~pdRdf3`lMH!=7*g#6W)Pj^ zSOz`AA3HD64vS;D%2r$lHtLd8V}%~9)%P+{)b$j$o?#+n61x%t8a&-@a`}(P{h;7C zllbt6_Qh_z4`~0RXHB-te z1N4or#J%UoUIf+{vzCe5V6innPyhW${BV>(jA{YifPf!(k@K%geyf}qJxJ49N&$97NxPR;5e#D%!EJo-Lm5;dsQC#y<`!wUDsmE!w+# zw_-=6yUwZ#THqtz=Ax68uQ03}D1ZQnF%VBXCd?_4C1g;0eVEMBnf>_rEayRH$t`Ud zcEg|C0EwxQk88wv;Pt0`i#%pk^481y&e!0Rpw4_Uw zq%$W9T2Q>i`Rt#0KA$f~3pxp|fA^%v(@?dA{IzoO++0Mdt&CnTaUPi0-||&db^+Hh zt5R^Q2d>1~DA)hcJ7 z-AX*c)n(CNq($l}BR(oH)H0^t^g(v3+UE|EE;_39i!-TzlLh9*?1H%zSmroI`k54H ziDxj`yt6R0-MZ+hshB0|azuQiG~Q}w5CUPuPmAc6ZTyS_2N4L$WwkM6j*w%b@Fgwz z%;2**Nd5h_a>+fThVYal|FO3{H5%emtwotA5P{g^3 zuc7&DYG>%D3C*=~%mqhxFSsvCsF@2HJpfC0iub!DLB9HS0NJ&K*`m>1Bb%$Zzx+K+IE*MSE%x%_zhRl1;&E(ht1u1ia?mGYCzJh1cLA+x_jG%m=J z(-KL~d!03v_b_}jk#9pyn<0V^T&1bV`#z+l zwwgM39b8wbMvsqT*v7pm12qbIe}Ea0#A?dBJ>hP)%|pIglm;>+1K1F|t9VP!KivLc zM@_wE8ypjF4h-D_en=Ag-l-hHJbS;kP!o0ATVDYxz`+k{%pnlRhgXsW`X^(|5iA>7 zf1(&0jFW4e8`cyzStCq}8{9ZPHN${dWhzf9k2c;ecZ!$HKJ4Py@Y=KUf)Uah9?w3@ zWB(lZvldQzJEojmE0>nic$D~AifYdSQ{<8*z)hwmg`b!frJZnjBhKh`zX(>OU0Zm6 zB>mrI{YnQU_`6Pp+T9Ml>UFvD5NF z#s>S1xTpNK<8j zmZFX(;G9@H5r?Fkpe!d{?gfegfGkwtS8EKUB*^k^-SPlra2N_9^*?O=bMXQF;g@DNg6Z)w> zal$WTC%}afFFVx|Ni{6Go+C<+FMn>T++xu7k)8Jkuio%Dz?B>Ipys=YQ(YP7J|9Vf zwUq9NjHxvi*sYPaNG<7r#Z|$BQI5O9e}^-YdctMyl*=2>!{)xh2P9Wkr^o*_YbTZX zfOP*GBg{b2rERKeU_d^B% z#~5b5adwCy0u17Y4I!f%Lt?+@|LDlDG(oubq2$)1ClWpTQ?;M&DqQcs2aBEk$P#+( zcx2wU=pakvvloLP-PrfT@IFt=U7{S3PethQA%T|Qo&O-91?9f?l4JvCqt!4`&+z1i zEIt021CpT)sMyS`eT8+JrRN>X3ppd#+=N4YQ#*~#`_cj*C=oR2^<%`-L`jc0bUJ(qB+q*}U@OJAJ2JGCzn z9U0bqu%+eee7crV@zx8K5=3)PAhzaGu;U+;XGWxx3s*lM5Bv2P&k$i3u`2yWZet?@ z}E>FsaU|$_Y4TFG)J1EvypgB5-v#@Kqcp z>oDj%`gMn`*lN##F6c};vwv8SJ8*jZj8*Ejp_AO(JjcCRxdf>`IN?tEeD zV}kE~x0Bv1ZjM{ZmXPVC^3SCToD$dNlQ6G5d2{tr@^O_Fc}Bc9_W~(j5vo4wb1=qt z>F~S*LPq>7(xiTHl(PKzxUcPJ{2tDRcFsdtV@laCY9>v$xw)bb z90XLIc88XuDxx4Z5TZ-Qjeuzg0RuL4Rl8R1L*0q!H}d_jpCB3Wv&<3kAcH!d_z@n0 z(|h@a=NU0;t-j@s+QL}yPe!BAq`EaS8>r~lF`N5+L7)O;y&E)pUAk>2XMaF$yw8sWjjc}4oa>x(3+27!};(B!L{aPL(=^)ArTpYWc?bzwBw*NJX_ zQXa98-oD_1gqaH%TXQvz{8J-)%^^Ws5Gd78@cyq4?1(uQn1PQz(9mA8v1le~#RaPC z@msUN748-xv)i?W6=02c#5HA_)%Klzmj%|#2(#%KJxn+oXfaQJV9ua#9YoxH`(cQM zgz11Iwy_!wj`jYC{Ufe)5d^ckI+KO`(Wf5t_|MLQeA~hG5=+7#4io`QA!5(1IFLKG z^!O_`7qWogHyNE}(WM^4*O{m;Tp3`5^(wkx+S?Dm^!=xw9Yt2doxMkRKc(vo zj~-mum`Yz)PJpidjqgcLIf}m=0PPsGZjx{_$FUNtJE~w@3$)CYA(U@u*2)ndfPwmt zfAoT?1!{c+JwP=G-{}YGO%A(!5MBI+eDn;!UdpA@H=j4!X6?RE3H~8(cwA;WAlrNz zWKy3fSkwg`Qm%aO$i3Ra^OquC)HDqt-v&Y|vkviIzq{R^ZxtEXRycRPtmB_vFqkxQ zUjg6Q2)xMq1$B4r^m@|43iHA!sRE!! zgb=v#8oPV<4!gRXlANII{+FP$AnU;oa!MCrvj6AVCHep7nWB7b$fKva3b6n9^@|rT pD#-tzXYvaFxz$DRf2aQQwcxW$O0k!XZn9odg0Vs%H??oE{ui=_j*9>Q delta 9589 zcmZvAc|4SD*tT7g7TL-&ib=9%-$ErzNyw5dQ%2Sq`4gOu!wy#V2x#muU#Nrv$Uq5fkZPMPzWt>wO#ntlhr|{X=)%HaPw@)z+-~VE8eW3m+}UdRUK_qid9m!6;!?M{Nf0H95vF}yhWBxkw}R-~ zb|EZU-OzkfJ?UTviQhBN3f9+M)UbO@!zBl@MQti0xHz(q5U{<|=u)@0pV=yLKwe)* z*Sd7qhP$nfw8; zB`9+RdEyl+0@zM}MOmMCRgG_G{|@iigRtubxek2pDvWRly0*XmxcX#8>YVfUfOUfG z@QH$6m%aT(lC!gvDj}Sbjt~0iYkfV9&slORZRvv&vK4l3e%syO%XvsnEaJ1!oO0AA zp;pI1vVJ2M@kx<|S=49qGuQCro^r_2jH{CEoArgg8sLWBZ=dyLlP8`h2+2!l*RJKj zntxL7iGM_{^}d72`Yofm*Vnm=7dn&|cM{H_v-#N9%?d>rXPMWBd$sCk-EY>l?hJ1g zs;X|JPs{8~toJii=o7*>(w=vH9bp`t=Zds8EZL^G&AcIh*?)NNRJhiO)Hq3=F6L~Z zkO|C<`W`@9vw?>EC)*Is6+t*pe=r>`S*_tDuhmetT5Ts(9E|z|*|B`uF*VU(Nza-C_SYF@cig@lf;TaxcMxft6m+clmYsD>jzVy~6foPvx zeq;1p_G3+@yq+;Cr}tjvBAOFd18WaRkzES5Yk3er%e3vx`nrH-Hf%D9i-j9ma|%|^ z^H`VtXXv>t!Q{_?i==X@E`dp@siBH`enF+Fq44Y;$u1oC4X=K%uoL0b^SRH9G*fNH zZF_d!X+w4KcSap%bzn!PhjA%%2V)Me8OXHN^j{x_K`_}1n_~?!xFA@^$b-sF%EHEe zTO}~+BLSV8@LiiY=``b8P6?=+=hNeKf=BrjNG~L9=p$n@TdY#o)fNO6S}4JI9sdG3 zl}MYYj)jVFB)TKLZ(w5x$-0Rl)s{;)MXi3TenXxKa8Qld$+! zx+PM_P*7xX?L%ql>L?P+Uf@Hxf1B-Uz9j*L(OlAUu|T)ZG;mJ@CX$K>eotpDMr!0c zU7k{GZ6CHpG7Ov8%t~Zyl)m#CR!-fbIBZ7Bbp+etyFQx9M|x<9oPZO$3R4e&@Y&N{ zog|Sgy;}{o9i9I4ZMMuVH6zQds<+=pInxc=!k<=>cjxEjg|#iNVXy5+sUPGHAeIOt zdDf+MBmTH)N6qko>;c5cPIEDyfxgtB7;npITW4MB$K8Ii(Ux1rmli7)rNI@5)pe1W zJcrV45onyY9XVc_o;{^$Jk1Am*lN-*eee3L#XO%U>${dO-Dq2?S)&Fn<-H})XzRo< z65D55i$YpBmE#$f@&*unp_}i)-!IEwSg+P_?Qx-yv6_W;2+S-;aZ>TN>N$j7WT>@I zd!KcOroSYrAVsUaGnv#lcxvxlIuf(amR{gj9c|}{(SYfW8CJk3C)D zvnH3=aQGsdZeS6>@l#rf0P*$Fuv@j_f^;u-+<7`uD>9OqRf<<%K9|n?K>l~$PG>{; zdQs_19~)q9(mZ+3i}Fg1NLjY;c6{48t&8@k)9cUM+`|qle@D4p3ova6;!5vWu!!um zI z?burj+%nX(i_&x`pWl09?^9&|>!%8&VCJ~2#*FI4oy$sutKnrF=swZiyeyB$p6~0H z9lxaJs5t}|NNd>nBT09hsEZ2T_Z`IsJUT-^9 z0jh2hUDDg%&1mnPU8+q<9RFq`_kI>Z(pt0i-01abrLP49H0*$) znY&$hlC;*$`~6!r1-6M|8Z(aj+h4uEmOYB>$31dPM?BH8^T*9>t9Mmqa3ocWS&Wae zrHk>I)|hwUIZ#&tIXlgh_Xp!U+v5gz1`)X<0|5UU7_fIQ zXxK@APj#O$cbRTo!%VA;1+t|pXJ^)UeNefgM42Ad?(+ME1;71-d)adTi%sISsou$M3geuG~agGHa&RT^m+7d*C2^^ zBX<_j3%;FsdnjMUB0Yg8{egJLY=G2LVWrHUwJ(%@1K3kxBM%0hTuWbXw#_aSw1{ZT zW9KbtaIFvPFGFJ=+pai=ta=A5`!%ms>^x%3p0?2K?6$7YPL& z=;jEApK<>5C7a)Rzf1@Jf$Tr;!k4}(8(g>bSX`*2sB{O&P{M~2yYTC6UHEBE6Oq)V zY0P4c`6eR&5-{ekE2d!wCYK_mL~tZ^cH!4h4tYPb;VOon7)MCu>bKqxyYIL4+Fdm4 zN^#PpTYbVn`>dtC#0r4hIKPaYLQkYQQat>TvUcz$wfl?VmO#$z&{hz;o_M1`6SqH1 zOv`RnNK2R8`anxjDYHzgQQ5x^ic_Cg+m&wkz!7{RH&~e|ywilyYC5UQPbR7pvocmE ztVD2<0;<3?j7M~&ocWu-!_CP|-0vp&zqi;L0J~3HEOU6GsrL#*3@XmD7 zUSC-$scTESYs4co#(idPM*BfGuEyZDG66fN((TwV!)bz^Lr1v;L9m9&o9C~Y=k!LR zZh+8u(5C3)P5CIS#JRY3MF~}_Av5J`vDQe5Q!!msbYDc+vK43z#zl;3`^&28dMpA9 z8>|I-g9YoGt^POUioj%`DuNae4>}n&W`eu4SMy6$I^;nnlqkp;D$w`z2BIGjFgw*c zsY~X8tq{|4AFB?+ZVEoLOu*f6n;Wu`x*M~>I)z>oFcDd?)Ae^*^Wl6ju@Z`{LJ2FO z2w>Udq-()a_H+1WuHn>`M}f6|xYSN}rk+NFM|v8sbYpPbz6Lp~n~zTo^$ge zTU6@Ci*x7&*>awAf`=V}G~eEZ-|@*E!skj1oe9p@kT!*=lR6@_HD)AEbL2^O1g1wE zD1E1~U;@}4yA{QY^XgXYpBPhTzWfQ1apC52ac%%Y=wllUAv6$A0=5(InaaEzTccO;U{;Bfnn#TH>X>L9hOR`PCF037RCva zGO9layJVE)kjk=17b}4it2P1>TcMcIrE-D65f3)BIDLHPxo6CCbrbc_ zTI7zaf<}sK^9R0%O!wo$0$#f63Rf@+zi1qu-%W8xEUGc zb8_rb2sHGj$;{*TH&ZLarG;{o?&)i@{l0#obj8)^La+Rzm5qv&#{f!L5JF9{xIoC| z#bnp*2x}LN!`DWyp9g5HNWihE6b|OqI<*FPPzY`}^964AMCR5xRWfiwm8{{q-TmsI zrg_&dX8XWCX4}|&fUWTJN*y-CulCZibcWJBMs^a-!tmWRV7vIG@kbS-)1w2~n1M1k zHcst;);HISeBGMb0(T&ZWPbR8l( zkqq53={&!8efe4u^?XwMoyoP<7d;Q^V&@l=jYW=KYn2PuLEaRfKR-9BFiK`}p8Sa~%Q%W`EdGSBP<(z0bIx2bsU-)YwIP73jgNS(|R)l_h>OG zEI0DeCwl^AVnMS|d|}0ol$%h&xbLh$lI4iL8yoBfOOU_jVs>o$`PBI2+V&Im7n$H!sN+KE96+hTn}lV02a2vk-k2fp0|K z0L_aTqPhOj$+IWv*6x-BDR>oD1;5S9My}*@ zJ~<3psQo7=lh3U+;};XBMxAHJIFSZlpQ*jeQn?%7h~wLeg}$nqa7dJC7jLni2<&CK zl3`1jG?Zwyt;HuP#JKCj#sUnO!~tprKHSB)nMW9Y^U$GpzpM(=IDxO#u#%j^Y{&M< znwI3%?R!28a>af;+PQROq3wrB=Fh&ciLBU@)oItOw%X4|fuE9NyOat1z8fX!b%_FLY!>Kpp) z=$Hpo=eQTHN0j=6cZ=L6+yC6S2zz2nj$$aW@e0}_z_xgteR4>(m+mGtJ3zW1eYTE@ z<)WFvM{d1>$u-`1XxJVKM9D|ge<-=btmm^8xuGTgP4W|5XpD1Am2;Qf$X@AnWwnto8oCITJng)(eVKzI5ZByiT|jeAEA= zK&?E)e@F*NG>YrDfcC0v#cFg{&WW!)Eg%tf(9waj9ay6r2C;f75FHMmC47e@61Dr{ z$6@9N+gaUvSGHOHgxP4kcfM8>+vB5~gxCyt+&i40SNJzsr4ra<%Q5PF`!PCq1Uc&) zFjfL&{#j-td!HS6nA+~oBXIn=OwI=4OKb@2`OP$2##kZDKfNeV=|02dm`b_wW}6U` zY|@FOw%#0C(sts6oLh>8md{}CSnuDFo5swIc))Az1~4N@r&}ASSYLC>~37a%U+N&jv9CB?e4-A zYNopY^s5w&jPbY3`&P`Zg9#Gbx~D77k_BAikGQII^0nFp<)x3--<~9wenOJ-``=%} zjLfyMjP|*vEC?hSnRZ0*-un6Zn)R`}%*RQ?EJ4P*Co@I<;@`5=ra@NX^68JhN`seR zPhTylFU^KQ_+1;N;2m{$T-RF-`~?8Sf-{<~w_0Fl?B-HKW#x@WVpNi3Al#ZKD@$cEBrK}NttIx; z3bBALIuAPEkMtf6y|(3u3YgQt7s^xaPY7t3Ael02?mV+wR__!^d{R`$QM3BX`-{Nv zP;#Tz$XA=&4+qa!K@Be}CW%oIio$r0+Ticzd52X!qD(R_rgaiO=^4@}(K&Cmfe2Wll(#7$+Z z`8kP?zfNAwZmj#F86%pHjJZN>3M*mv&Z?TnXk)~VI2{#uE=Hd40~ND`X0_b^hI#nL z$Qs+%cv;`x$Hxn%64Z{SlVrkXTBNEPjp~Vo%6n@5R+((cN~w!lX_jG-(7_AUmqFi9N{8Fl!Bn2I^s1N= z=#{&ZXcZke#~b)Zi$~}|O;(hp)Na5~UFpR1IHUQD4rGSuN?8~jL7sSNQId17W7_z; z%a%UJjDEqQ*S+$=f7R%R3ZPY|@BU*527P_|WYvnVvHA7*RPj=&P1=vEBfHHrP@PZq zp{7vAXp>l#FHI)HIP>5FBw1ae4$<@!F}X6Fc{4H4z^Ndfg-02<$fcCwD(<`>Y7qH@A+ngb zRz%(6^>bFvb-}V;vh19t<5_pJetMp2EL~)Az2(OuKZNXz{Z(et3*X_lyFO916T}SX zld%@y+V@nPI6+1WgQ7OIPaVYjsOYvTXX3yX9S=7(j|dA<%tBGYONcm7TyKJnQZ@lK4Ddo!&fgUL6-36(ZSHXtbGRGq!Iq|7nky`g_$rKoWA>5Lq9zvoN8^My3HIFDMCBb)29{^v=e{RDul$2GkWT9z{Kv>7kkj+e*GiiODsI|na|)Va$>u>jK}Nq zMW}JD47YlZF-`Wr=7@(TJcAp!(`1U99pv7{*Gqzv#$kwApaiq^WN>u|D<>mR*%1c7Rgd;;Vufh5G+oX`8z6~XgO>n@-sR0g7=@(OckG|rzgHQ#2)S``b{Q6K9^*NkI_0wO&OcU8i?FV8jXS+W*%^ddaB|D8rkc+9Zf#ODg6V)>&k4{PpVcuI=P>JtEhgUvNg3 zUJl=KIvx~r@1U3Uzdm?``TI^L_m@ z9memI%dOC4{3f_=a410V5~!NA3MRg$T@?UAJ(ycTPR{kjUnge*%SF?J=t{h6!>~VA z76Lbt9;V}g2Yd@{n;uV6*V-<3IXo?`LSJiQqii1?kFvb$f+Ol^o*Jo_REzTivB%o) zfuN_ao%Z*Al{R|;Qy{o+=nepwt@lY_cmE4K$$U4O zKS6He1tDN$$YW`O?Hw-5mH!|whan@g1ez$ttl;eoR^=pi?GUmwkt|!76Fl%1^O{On zxC`qgVltIM(7_{T)$Syyjw3D3=*4`w4#ppGr@lBpqx5f0xD}=;je5iZ;B{O*hwm)k zOu|%qeew5geDr`=5LMJsNkpi5)U8=dBB-kNo_uBDS;aJ-sue`7gHZyn4ohOa-=tM@ zWoCmR(qYjjr#uW&4{wRR1E%iBZUUm61hR|#e`y`5FP;ZG_&Cut>N+E8DljAvaHYjP zS_LMxc0}-or9#v=gx#+kA<3$%>0T7h|>#tnEce|XTWx(04XW>V19*K;}(PUCZ*ufDkC zuo=e%GQr6X{KRVT`~qjmbGK2TUh>6L6T+m=^QSxHSAWa!p6)>%PtDXsz?hROvjW3$ z(!&_;w-CQ=*Q&2*F+`HH zSnl;tPs?vv{*xAB{f!dY4KWZeY_6h%97hu&Aoj($=j{_P4EODDk2yyg$R#3IAPe?` zlss@OO@WXLLh-ZX>odWKH2?g?_i0|`xN0~nDZEkc$!>MsL@LFd0303T$ z5}VslBA3px7YRyuq|U8yiQDi$(l`pR3Y$6G5=`;6*yryx1tZk-r5Tfav3pkhen}RX z)_5Z*BpMK-3PlhH{)Yt6M0yp{OMx*Cgs!0Dwv@j63Y_A@&Mcjus=O4PWodUHgyA}7#vL-2_Ph&wwf z+;q8(o}~~-bYNb-DIND^0T-O7G`^Fj*t}^h}CV&<<#}JuNxvlR1r4|A9N65w_ zUWPknJzhiS6yydt!5O!Q4^fBxWl=9VnR*HUGwAC5hdWqn$VSDEAyQHZ7cw;z@a!XK zSy^}tnKo#WL{;a%&N(wT|6guRCJdid2w@GV`2EeoxpbeYhK^^xt$(Q5gAWZZH_vD;K5Y_&WXAMKLM#y=-z@ z?pKtDh(Br(2O*%(sDpty*l zR6#czV3}%Fa^uNihlL5prf=;9)AOK|{Ea~mZwOwZ$-5%bK%1rc&UA;)H_SI#Ny@eG zPPK(E*G!Z`-M@yXJ&^Xv-p~d$4~ko%k1%~?{+W10)OKqn54QN66VOq#OK!Yz^>^6M zY%Nv;cXx4plnR#v`sTa0)Jv!D+-6ael9c7YdicZd%zymo4@pJot1|!3Ga1?c&oc%2 zIA>Gq`wCM3>vbt9DS7$-=b5zp|K3VUN=E8`e=Yb-UOMignZ|KhDXHU7sD_^A@&5x% CO=S-N