From b6f9872f02e9a1d1e75e64d76abe4ddb0b01c415 Mon Sep 17 00:00:00 2001 From: michellepistner Date: Mon, 13 May 2024 23:49:58 -0400 Subject: [PATCH] array dim bug fix --- R/bassetfit_s3.R | 2 +- R/fidofit_methods.R | 2 +- tests/testthat/Rplots.pdf | Bin 11761 -> 11776 bytes 3 files changed, 2 insertions(+), 2 deletions(-) diff --git a/R/bassetfit_s3.R b/R/bassetfit_s3.R index 611de23..7d9607d 100644 --- a/R/bassetfit_s3.R +++ b/R/bassetfit_s3.R @@ -126,7 +126,7 @@ predict.bassetfit <- function(object, newdata=NULL, response="Lambda", size=NULL if(!is.function(object$Theta[[k]])){ for(i in 1:object$iter){ - Lambda_tmp[[k]][,,i] <- object$Lambda[[k]][,,i, drop = FALSE] %*% newdata[object$linear,,drop = FALSE] + Lambda_tmp[[k]][,,i] <- drop(object$Lambda[[k]][,,i, drop = FALSE]) %*% newdata[object$linear,,drop = FALSE] } } else{ Theta <- object$Theta[[k]](cbind(object$X, newdata)) diff --git a/R/fidofit_methods.R b/R/fidofit_methods.R index 90442f8..664444f 100644 --- a/R/fidofit_methods.R +++ b/R/fidofit_methods.R @@ -550,7 +550,7 @@ predict.pibblefit <- function(object, newdata=NULL, response="LambdaX", size=NUL if (is.null(object$Lambda)) stop("pibblefit object does not contain samples of Lambda") LambdaX <- array(0, dim = c(object$D-1, nnew, iter)) for (i in 1:iter){ - LambdaX[,,i] <- object$Lambda[,,i, drop = FALSE] %*% newdata + LambdaX[,,i] <- drop(object$Lambda[,,i, drop = FALSE]) %*% newdata } if (use_names) LambdaX <- name_array(LambdaX, object, list("cat", colnames(newdata), diff --git a/tests/testthat/Rplots.pdf b/tests/testthat/Rplots.pdf index c9b1942b62bce2b5c5426b55b30097a2ff44a367..b8f2be11b316eaea2424a9e9f11eeacbbe8c4ede 100644 GIT binary patch delta 8065 zcmZX3bzD>J-~XfBiUNX2Ntcv#Zlh5YkrI%QhEdYpa8#sAx<`XZP9#SoU83|Dqd_SF z*(3(fK=0r0`+C0nYddGxIoG+a&*xod&+wJeJ26f%A#q_*5m7O)7?_7hP|ev=(cHs| zg9b0w%#DwtgR0%pMK*u%UHUS?T zwjMZ3ox&VXIiGCO?6sK!NcZ;5-L@8!f;g;iiysng2n1GSrPM?RHSbEjKWp`SHia-P zZiW8BwZAOp|X%}vsPPnSgdnf$>d{^+H6PE9HCxw8Y!H=1+UU4iWC5gp_2cvD z38%A^)bxnLf~zb8p~BW3VVrFaMk< zr)LZmjU1`Qot0jjEm=s!Dto{@b_C0&F8!S$U}u+{Y&Be`hHzqj6c7rr@$RYF)(&bO zy`~b`1To0c7q6&3DlBo*Vi<1ZPc6b<;-rP4amRqK?b)wdw-Ia~8d%G4gr{|d2UOxj zx^1v48TA2Zce%UvL`}l8S<~mepX_f}|K6&^ek(au77wAihGO6!ft58^8{X$f<(`fR*ROn3xm%&rADJR;C{KW1hpH)nl%ho4kd$@nbv&qCp;>orCqy`| z+OVxirV+lt+h7_=As`|YI9PEVfG>L@{aw`7dD~zBkK9JoTbTh--%nunDKm2@y!-EiIWO`P^ZKB{n`5mZ%Tvz6Hep>uD7rx_`(55&-5;hid4{~_OKVm%MojcIzYLEp zNP`0S^E6I(zN)@?C=4&E?r)f|bUg7fY?R+kkv+Fxjm9zwc#=872X~fbIyrip`is;f z=^+llOYO6jSe`$-*V$bg(Vb`9RPS~OGjyBb*s`k+iZb!>IZ|rt&f-daLgt*b>s8`U zrjgt`CyN5&Vggd7)n5x=7T|DR6oo{yW@eMJCTP9R$m=FzmY>E;&gBh!ANH7rQyuL- zi<_AzG`5iL#J%(=$z8C9!3=9w09RQ7eXa4~v86e%_FG3i*_WG~-D6pXZiY{{!cfQj zQ?Coc1vcy#-%@NX%2kRQ|7KGE84gsQwB!^Ep$``p5~WkRWZ?P5UA*cJuMpK94<}_9 zVfg!}5FutAz3294=JZ;7szS(BjzY2F>vu;N@({LRH<|`Tv&)M1>bPQo&#S8zJy+IQ z7~W9WU+c)r($+FypR*{dd-6Npi!}+nQ5)%QEi_)GKfX8E)%DqK_|<8nsyVGG*h3{& zzqUG`v=`J(yIY~=F!;W}1Jnx=2d&mcO_L$EU%QfQkw5kDH*=*k3OpW)J>Z@*MhK}= zGR;v@IBLfN(c-nKhEz5{FuR1;vhqC?Qo2;^vfS`vLp8wBn*(9{wQfv`5IEcN#!Glz z^8Lb1)({&S1S|W!es;APu1QJ{_fKFOSr;g@w!Soz5+!u8>=#9&akLreq4-9PJPP8? z>G++;0#g4oZ8XHa-cwsT!n@pL(6C3lC!@*1*wvyRKvPJ4)k)eh*VFp9U@LizP!#t2 zt@jT^xK43b)~7dOL~Dq8X%R0djaRRzz^)VB%jFBm9)^(1F+#xJ#C?J!PY?s~iN&l$y7b0Gw z3#5=S)iAG@R(T5Ij-8xKxV->QqBuH;3y3hGurSmuPDfy~40I?alz!^-J(-GXowXv$?ktc2CiE&!Kq1e(C6~XNuCCpL#j>1=I?Bj_(8>Uq4nw75dP*It^O!z?t!I zx$?az>v+LmnDDRyHiEOJW*lS{GfL%w4SQczIb&HA)xF$!6R818bO3#c9$QjZ>92ZTHmmkd|t z5|9EGa$A@ezASEF7)enZwR-tyX?M^tKi{S3)DM1#e_i_VOD7 z4ZAH@Hb{~azJ31!^+A0!pUs-@+ouat1C1i2v@aHh0jko~*cUxDS)- z)YsPW<{Rk)7goyvSm*CM7|MU!rq+~hWS-A=qnEBnk z_T;~%Dn4$IDQ+zzO)Kl!?^PzN&d){GO1h|j9_zH06~`~R?TF(gHv#^wsR5_)OW7JN zG%gJ+%P$?JGbrg8Q+edRZ>R0USR1Mz9Y#JKQurUE>03$-p2180V{;2h#_d6~6<~;4 zhLo|49ot2%Go~KC%GzDDdf=C0P6T)G%hBBr&9B8Hb7n#J;^7*%l!C@&ls1@R)h16wOv}n$3q^n+Vdh#pvm!Q{$cX~TRsBeO&Gjr z{VyZUt1x7_$Sdg_4ns6oengA1i4A-RWZ3Y17Qe5iW6u%C9jIDfTJYg`?N_-1bzPpF zGpgdl@dcYa!sGb}wu;iKbkD0*Q2ilkf_Mp-0KmYa8fzvwB5m4Xham~xb>9C(##~4Cdgl4WHe+dPw9@x6}KLF9X0Cnxi|PJa*u^>&Ph9#nw_OynhkAl zp{9WKUVXeEb}y-ChZ}4N;{Ne3G_3YYRuA`nG|Py#>Ak@9=+|uI(aLy;G!;7l4B>Fa!&Z_*0UQ%16_Tr=*g_4s$bToiMny< zP*q8+5(=QBscLnSSk%ii${C62kqsKCLU4w=iU8s*UI4(YlnUhI9KDVua)=mDZ`>XqAVtaQ5KIoXY*Gz=B=mz-Ek*Ddv)C}qn-2jzia+M@^1oAJw%}W z(5!hSxp8A!aHQDIG%bl^6M$Ab3Obu9W}w-OQZUeFV5Z1BB`d zB+OAQb~LcR4~#kM2ukwziJ@BWEB97>lhMCB%t*5cvRh7_nhw7$%~vfAo0NqYhfe9c z8t6(A5!)Cwc~`+R{z?}6YL_V)XH};Ov+hRMsb9I2*;9<%_X;TaBR#bk_&Agzn(^I; zZwtsf26`z#go78AYgE?rFMpWqQ?p_*i{iRjyu3N-c{WTZx!KT99UH$eboA-H7@jr; zB?8_bKk0Z$D^rM(`-CYwYvQ_>TW!3JhtH*A7TxZ==C4^d=G7b8UE5PFCF#i+RuOqX!=%PRN zWVk9y$83spH;u|?nIFZ!R5e%D#ClJZwRa4EfH(3N|EExH-M3hOF!=7!ewBpso4Ta; zTB&bvn&$Bk^dN%Csjkitq1FgX8Ls;HH zg^`K?m>!1`NarXh(@n9plX{{o!&)6GBd>2P_XN9UXE1zAeP%*GOn-Zl4v>0BY8mVQ zev(T)#e3P9DO-k8DplfxR=ydkgQ3f9yn*PCDIlieHiw?iYyXuRT6n2};_;I=eUI|& zwHfJ_#wvpENC{Mx^u}~v;jwNZ`^ei*mZTaj?&(KJX(OaWT~{Nl+2oI6_H*jLqzb2X zNc?LhzFfdb*uAI6i)k&NykO-Kv?mlrxw}{=hTQKc*VEU$(4`9&-TGeJ`XC-%hYLfL zgco;=2M6J#sd)$8S^7oDi0Na?HM(`Bp4F3=16Y@%G`?M1S8`oSPpO9uGm+u+3Kkj? zj4B|Wbu=Ev>Mn+__D;Jb(C@AtX=SAV`~dsag|J7>TQ;%vvI5#a#b=Oi>PO3Jw=N|4 zM(?%H(Rad!)64}?`Cz~5J}$`C1tJ{KPzQgs-VDfL@bg@mAbjAkvSQf7lAo3cx_*-2 zZKO|(XXr9Xq3WFdW(o-joTz7_TV0~NijztNt{Qn#1XrvJA`jW*}+7hvX6E-}HlUfOs4iyKu?1%qz$$HWtb zhP>9@bqRcVq-;*^mn*85mtr(mo5^2UF6UHB(VrJ)QfjbYFrC)uZ@YRlv2>8hK1sIV z+E#M|VYt2IzjgjD++63&vF|gv@!>`Joz6jtOr?J;PplXa~{N1#b&#xlC9~L4` z>c{7jF2xiK2PG4fohW(EFekCz;)yBz71O)5#yh!qq<^=DAw3sf{6V7>0CUkZgijck zcnEsu;SmgZL6c6)+l{`D(I`QKte(O6co^5<+nI&0b#ieYiBgW2Fe1@`!l{egLx^EtxhQDLgc%Qx-=V;@f zyu!H72!Ev0YMK=fE3NKZSMgpw5eI>-XsUC(MomJbqCCQOcWVxyQz79VpUA^Rz>lc0A4%aisgz8ArZodAR`14jJ@y~=r7lQJZf)^Wg! z1~2nRPX>Y~oPapMY-G5q`pe$L>F||Irb&%EUv!N2UjukO5s=I>(!oSwzBw3Nm%bm! z@oyw>VlEMMK7FtXRAjvt=pFh(HhIH`as_Ku=wKWqZgmt9~E1>d`mGFhJuZ=^$fF zU$e%{W@G_`?iBH?*-4-E{R`8wO!TwIqsR*CHKTAgw<@bm=b$I&W>ux{3w4R-UK?u% z)m(Y-Cq7_N6veDy5_a63MhtDuirnx3K| zKWi5U2?nJEHG#L{wNuwVKu$M|r@^r>Jm~yCH{nF(ytV((>!&!xS7J4+4*u2~R(||- zZBBLX!Nz{$O$eP9G4az28%peyG3_8kUu2ks1b>;~^L3Y7Iy!}+THW%nadU=6xTS>z z*x);+Z4^L5>6Y0kl6?#25i`hvH&iZr!a{rqEjcDc)l@C<>UIz8~1$UpCmGpbO? z;9q`Vo$ew#2G#FV5h6Z(zax7G)jSoW_rclNnNa8OMf)(atizLVGWSP#9@JJj(aQ|7grni$CMDx=ZB)Rt7X+^TC+ z3?4uxiW(t=Z|*}0O~O6nn#8ad!tKJk>fC+vDxl+}piId1H;3^WF-=O>a2rd4r>9nP zWoHhO8{(GS8EGeri|i8QdmFPV@<%rT{Ffo+C^jMP=ges@B(h6dAZnau-d(-yssvA4 zs*Nh7^A=QX_#dS#`WTa*`6p84C$O8bmlMYGk~XywiFNA&&VsZ;+wL7Xr0Nl@H>t|# z=M?9?-65TDah+keZgU0sY4y!QgkDY7ddiZ8aQLOGs~;0!X=^a!D&$c9$*N9f z^wrp=nX<1^*C%*!M&g?#2j3fXNM^k~STChY--S`VAO~6Ay{k1%0pH2m)4qL49SBRP zwL3JH?V=nNvnbIHU&oEFn5F^-?;X5q?88UyVVK+5fmOY@|73Ek4vBXsh7_$}i)W;RhbKYJ}z5X!uGwm_);bldB`^3;iQ{#r{Ut zfC0RQ0Ox9-03mcD!H%>=-iBzHkbcf^O69TEYnZ!#6hYB(i5tMu!qU(*SpMWc0`Hxf zfv}+R-fihMqu(cwj_!@2CZ=!3sMBc1R*y4IPB(k2{>g8g{z>kH5nX?unS^F*>bjR& zgm>05a=!hQ)29OdRT?Jk9i)A-GaZQR@>gDpCHbwqOZ~b={e%|AP&LmJ_>$OZ1sfCq z1$={|GI9DBSv4D`Yv+`|EdYlRJc$bOWhF0*0EB)Ls~v6MKFX<@stN*L^-;8j6Lbo8 zk5pc{mwIJd!28ziMv}hy3LkAI7B||T<~p>EB)@0K{D{y1mJrst+j_o+)Hs)5LT`ov zz&5vl(;@TM8KZ2|-u#xIVO3M60Im>sxEW>&oMbZRo}6Dn+cxO8pX88^7_10W|HPC) zYhw7%BQ4C&qkUbmGt>?^v$LE1;V)~}^5wW8zaLRDVRrkmLsZKAOQ)wLYK?tX@=St( zq!)9tKLYcjMQr!`_miiTjl@5UULU^?%5FuU zau>vuk~wGnZfblc_kMFER9))X4awj)44xwcjpfER0j?+eW!jz?I}u64vXz608ShU& zG;zW$zGcyVag;&UAN*=g(}dyF+^twKZcxZ{hhbzd&Kmw$4?n)DA#LfBKH>%_1J-)QN0?sB`` zg#=MbT-DHA;^y2&xpd_zmby(!c)WWGPk2J|m@BM&;BEeA`VJ+F6=WTdR@2P&V2Olt zNQck<_N3OAN#}^|Hzbhn41U$gwcNphAXR($Fz#zF$wK%Qs?!yX$fKX}U)M+@{F3G= zMtoZHX;=FWhQ%(Cdb);lZ|BY%^{(QuLzi-AF6||@ochg8#92B}JaTNDSHG zZCt9wh76EM+IHH$-e0h>1BU0!#20+xX5~|}OY;$K^~nhz1?+37|lrt((n&U@--7&gY=LKF_X7X7aNn_#UPBlgJ2&RJ@<-C|Cq zonC4C3Ptk@2_!PKctE+F>Wq~t>^Xy$>yvIkjLLo4T6n(}W0akF3-FjZ8tZ$mB_?89 zs_BxQ7oWh&5%VT5!5iZ=oYUID;J0(Acp%{~yyY1Wqb{r#C zm{lnEedwa67b5WUF+lr_plHfT{!!%8b?!S?sxTmvjn|`yFOKzwGNi8RTJ==FzG1Ax z6}CR}Hfjyfd~JM%nN4to> zh(M&s-mBs4o1M~oN)Ek`Xl&g($`#=oU?$Fxp8rO2JItZL?;5bhTa+~Xv36j%Nr^nm zadelyl9SgJWJ>A_Sy_FT`q1QT%6I0JwbGNWB8R|E!^~-2U2p#>L0A2r&wesyNu^A_ zB1pIK9HYk+@=vD64NDf6vc@JYC7=^-5SZ1f&tKD)r0i8Px=xzEs;k@tEs{F9!nZbi zGUIVAiekBOiVE<0O-MC=wAB}p!)V@(T}z_l`1D|$FL~uklnCTr%Z>TA#ZSYb$n##G znCCY(%w#(4315H1-KaXxt2aSB|M-ag4$|0k{izML^3-nO*P0w}KOVxLEf1G5 z|FdJ8N?*b^ti$^4?ef*+4?^FpIvacED_=OL%gh$uliJep%V{xYUlQ^6DEd+k`tuXI zEaGo;^+9vB-uGMnWYZy1Wvai5K3y9H+LNw6ni{~k2C}gW>?1CTxL1B=vh!{wWiA9?h ztJPgjS#d!D!T(*0G~y)?%2>E-QIOGz_F1qytFjAZ23X3M^! zyq@{8xaAT^{FO!wBiHjwCq5mkvw8f=MdTo{rRa4dJ`WE`hj{z=5q@**-s;F^!|BRi z?Z{zW_*u<(El!{R)Ne??w&S><`1}YrBHieJyuQ?Up51t!Ji_mP;9eI#@40?{{xcXr z0+Xx0r?^VsaAaibsL~0$*Lk=(GOaA^dD>h#k{`M+)s`Q(wb8p;1~0S(IKl`j*v4zC!yjIpi3a5yLFuKhjh%KAoT zDHRv%{M;Y@Q+o-6drZCgQf)qyS?P_y(V$0lTwVP{Q_>q)#T^?utveMBt6B7vMzu~I zj?^ggwc)ZsCFH9;w9RdsJpYwU`Q1R+L}pfdHrKBH645f3Q(G->0q2UmJ)F<=-kRjp zNWU7QW~>rhsxyj#Nvb4=)7I@b?3s3fryU>I?MP@(gQwM!=VyfFRu{#VRs4Z6A`3Yg zQR86k=~IJI@_8XQ>hS5c-^UkC`_z4zn$^?4?NCcfBD+n+9_t?uHaVi@|J-d0de0TI zEG{9gZMKR(0M`#|duL+8_>%P?qR%j`7ljSH@un^QsN`O)p0FjK#U%#MxP>nDDGs-W z;vlaLoFu#3EGF_OY#0zFe1#=z%PvX(H@UENdDEG9&1CE$=HF2xKNY^(0L=Pu*(11S z-B04w;9eGpfaIN-_mG{lgYI1*@sY1(ttOQ3rWQF9=2@_C zE;9{sjjoF4Bs|Vn(1A+z_eUu)M#^skCc1_wx=ep8x(j?JPlh*i26Xx5Gth)OlK!a9 zZe@`pP6_1zbd3BNC+~g`xr*^^oclRJE76F-XWuQ~sog5Rq!wdTZxkw5)}`^pl(^{wX&8p=LWsG$Or)VBT^V=e~oE;byE(Ri#)_2FD87BzeqizF%ofSp@k8&XG|*;2Tlt0AwwkkbZ1s?3`9 zG;9j*p49baXMUEi>eh2Yc;w`r`d1dIX(GW5OW%Je2%}%E8MJ0Uy+@qpe?*G9ZamdY zt5m>wQ3mK^DGCYEk0{q@OZ$D^>U9=;B#$U}@bCvGOH9-XMqC_$`EOFrDI1~5lwG)LliRSM9 znFa+CONuT#=(EdQ+q@JyGd>wwD4AS$Nd0PgWk3x)l+V>U33GlHK4q%1f}j{h{k4B~TVg$@o&9 zGj19c!lpNv_!M&jjX zKB`~ANya859PKs^K0F1v8=yEOoic^7CCR6aIZ)hxK-5)L)tXba?#Z?kbh9F6xd8KX zY|&fO+6168kE1-|eyOVM??OV%y|Vv7mGluD>U0n_tcEDO=pz@`FPBW6l-}iO&z-%7 zc4wIOI_e)79_=id*fp_Q(xcu~*sf-Mbi@9EO*|NWe!bCPGz{c~k+g|Fu0Xoe)=oH5 z7sw+g!(mHOVAd>KM&8{$f_!h1R3p_fwIpCa^8|4@!ax(#WP0 zya+DaiI%@026;R~d0+-VQv-O(uXc^^g95VtO&$&!n2Dvk zk#^55TRrk`#3sWzrI~8bQQKqSl<3uv`pTY^ljkQCiS$|tX}1X^RnA(({?xlNk71A` zi@*B+Xggn!aU5%ykG8>6>nUEm1bm{Jy@+B}OnNZozt3`>Q!PR<8l%?v*}l#`Bm-vHJ&W z$oD0=TWS#XTD|+S&Cj(U6;pogx1BlkC|&X09ZvYU1!9gft5E`v(fJwwnhs`4KV4V{ z#(GrW5M|r7i#K}(^DKWA^_$MBJ1)0xPculOP(jqB2ofzJT(eMcD>@a?+>^z&s9%_@ zXHVNo0B%z~+UD*4+7ePnG#AafD>pT7msJPrA!C6Qr@3srr}1#2VrBKvgNGyz2qXZk`u(#9jL7#Hi7tcm)k>ZG)|cN=>EajR&seTW zGc^IxK|%M>Fm`oWL&$bl5(P_2K^ zb|$rOkR!3gxXDSZOt_3&t-&hEI7q7~gyuB-dBD7S;euDK1nR=(-@ga+Wjlzj&8iqk z&A8EU=ZQ5q`v(dFX%k5yn3=uwUrU9(#Ve34dg&LI&Q50MjxMjVz)$$OE z$xL8|Uck_@ade%5$7Q*mDy{l{wxFR-v;GhG`D5>#R+Zb>`8sfdwv?p+hn~(>o=wjr zYs@sMLwzDsI!sGyAh>}y+Z4UwKA|puR&B=l?E+4D(eIx&8>lp;y!h}+wl7|;LIV3C zo1EJ&7D#~KHK2C)2O{9h_)s$W?Bc^MYiZRXqas zSL;AFl!?gYb4)zf5kS^M z#+KGFUtN;*??eNe*tgXCcjPYRN(5gibpn3FsfR1&QI52bD2x;m7}4ODurHR-P_G77 zpP9XoQ#`>nJtP6nF1ekyR{H#|_XgScV8=jh=-Gt`|2 zop4mwymz>U44j$^VtwN*Fe@~N2}-Z@WECzwkvWV*r2`kVpH_1Ny5`zBj)z~Q@!vw+ z>SQEUTd-zYvS`VT>_B{AHQW05`26Mj_UfqRbFNNdnjhisMBJYHhcMiK8;F{17#Za4 z=!1yGMh);m8%OK6<4_RuoaoavWN^^xYZ0mKf2kqZ9*d5|Z=734-vC4+$! zB_|8|fysD3BOj-Y+Uow9DLiAQgU$+5cX1vIEV?55)B$1ld4kv%k>ZGyBucbQvzX%H z2VG%Qy3vuqA| zy75i#@0_ZvaL#JY_Pme6iV$I_ORJoAf5L0ZP9 z=S}Uk&AG&|1yO2R3%Qg|3aV@&TRsOXZ_}F9ZP_6Yx1-Zhe3a~`r6p84k3bYmATh+Y0tN{AYcY`mT{ z3&1Wk-zV}{WGD!**X6r|qSB#?+qIZON=~I<8G|^Xa)oPU;7~cqCcU_FC0vKp-}_cL z1r<;t26ZItNXf1?CGv~C<9yT+ye+A%|KHjbu$>gJtH@&jeC%D+=W|c z=&keU6UAe2M}TB-DQeDJ1rCVYv9m7wK`+--kJ?qc)LRWMrE+afw!ae@DxK1 z@~v!3Q;~{i-Q^om_N5o0zKvHfZO2_81$h|9oW9aA8B7O3gU;Zwr@MHuw%=<+x!`ZHok2Cq$old^Qr}xWSoEQc4oX6 z>kkkOld%55+F|3XdC&3iLr~%c#}krz`7*R?v`QY&Xt!yZSK`^?P4Pn9hf%__t>iF- zAKY;(fV`o!WmnCyQk!p1TH` zF#kJ3@>MWuyVtCZV$S+DBJJ2p-xMCskCazC&_K`+|5oeIViKSt!d|0%Ru-I6ZoxfR zrLOk#8bRtx(VO`TA&O3ShOTbCsM`IyoT|{~40*TA2eyw2p^8;#YoB!xxD%|YLY%xNqE@VueC zdN1jM-%-~xD!Yqi-lR?G?fGbas~YZiot-)jhM1iS{a$XxXW?L=X~T3a@u^eiXvYn+ ze$xRqm%lxF2oY?dYSu5|$3`D-!bdrUwzYW8Kfx3jmJ zQc7xCSTOl_2mF@U;sN_+d~WhLJ++;c2h);wK-QfI7D%T?5rG)yU(YCtFQu^b%e#Vg65X`!}LD$hlRef=( zEIg^X_nIm6tZyFtp6h)PrcqgAbZs2(WKGEIwr&r->bFoyUaQ}`zhn9y~;kSn0ushDnr%Ux( zA7#*!(al(cF)uAUFI5Ly$*1%bTvuC=}4193co2Dd2}~P{~crZq=bNh z&;heVNm zl=%blfRSkXQn9Vqhr!%cPFjZE{+5Vox7e)adiBFlZ&|34>0$o8b_;gHQ5F5toMR*X ziorj2XS{)#7}il_1Oyi!PJc6fIjFvJ2c}k-HOdTS>E!l+TB;41v~|h&rj8P<#<($8 zhgCZPzfD#%g_x=;XU~TF*)=Z&Pb5zBdV^HD(tTxV&pu*PL~%hPZVPhK&P{E$L5|yB zb8C5-O@48lQ%k%@n#NNQNrcb~q!0 zTPC_uoqF>(1a9bS1USAGRZDWZrdxTT_AflCQ&XlBRp49bEN)CpsOEw107^@Wp{vUa zF_$bqC4}XGc!#jbVfM+tRBI&jula3L!*h~#vc`*JCbTX@iU{pere+z)*UFs1w6#fH zRCde)T3T$LmGEXhF&T@I?g(IZSm9%2Sk=${Vk_4CutF?ylwX^V+A}p~;+M*KM8a?1 z3*0tK3A{~I?xOjgJXrBQv*9q6bd(Gta&DQtJdh3Q*o?@_zov?OD;*{4Fc^jRcSe)` zNdMt7Z|PU zhF2!GXVb96_(nT#=)ZS_Ip1spPQ*Y5+ON~})7GjBqE5P3_nMi%Dcxkqf5Je0ezI+C#L+uYe9$Rn6YXV2%nIfwSTVfa#! zQa2&~G?9*>)BkD^MEVdI#om0#iW4o7-0!WZ?AwCsMEchHuVkDHA3X_MQwR3aWCM3Z z7*C`R#{lvW^L1|EsPc)vq|&-BMN|DnqblVWW1FnYlGEr(AcbBhO={cQ*UfFQyM2 z6DSJt#hB%p52`_*W~b=B*n1Jnn8Bogoy;+X`nu!hCbZq$K{S&8a30o)9Y^X(<^k{H zi(@pM2s3(!T&zlfv7cm5|J^hwH`niA%kRs}pS3KxM2FsQh>E%gUjHaWsU=u~wr&uE zO8UuJ5?yi6?C&zG-&0bUrhU*&k}>p9#OGbq;v?(V*VXS`B`_!^c=%qKm~b2%QU`o? zyprl3nXtf41mQpTt-M!1i>MWm8`@hJ-iR#A?T9>fMt;7- zQg$bK&rwc$;r{SH=eH0~IY~c?1yBUA65f0k5iS#?`kW8zTNd?dvo2YB?XkO;0M9#0{D1z8e#=xe~-TWxf^ijlF0e<#@kb#NyGnT*a4WOn&EQ;nCDUv2>kZvA`ZqrvDWeU73UKYjQ%l4TfV#_d`-VN