From 577c52a534adc9dedac8a493f1e99d27a1c9ca39 Mon Sep 17 00:00:00 2001 From: Sebastian Weber Date: Mon, 20 Jan 2025 15:16:59 +0100 Subject: [PATCH] prep release --- DESCRIPTION | 7 ++++--- NEWS.md | 7 +++++++ R/mixess.R | 8 ++++++-- R/sysdata.rda | Bin 6088 -> 6085 bytes man/RBesT-package.Rd | 1 + tests/testthat/test-utils.R | 7 +++++++ tools/make-ds.R | 2 +- 7 files changed, 26 insertions(+), 6 deletions(-) diff --git a/DESCRIPTION b/DESCRIPTION index 4e4d6be..a11b10a 100644 --- a/DESCRIPTION +++ b/DESCRIPTION @@ -8,8 +8,8 @@ Description: Tool-set to support Bayesian evidence synthesis. This for details on applying this package while Neuenschwander et al. (2010) and Schmidli et al. (2014) explain details on the methodology. -Version: 1.8-0 -Date: 2025-01-08 +Version: 1.8-1 +Date: 2025-01-20 Authors@R: c(person("Novartis", "Pharma AG", role = "cph") ,person("Sebastian", "Weber", email="sebastian.weber@novartis.com", role=c("aut", "cre")) ,person("Beat", "Neuenschwander", email="beat.neuenschwander@novartis.com", role="ctb") @@ -17,6 +17,7 @@ Authors@R: c(person("Novartis", "Pharma AG", role = "cph") ,person("Baldur", "Magnusson", email="baldur.magnusson@novartis.com", role="ctb") ,person("Yue", "Li", email="yue-1.li@novartis.com", role="ctb") ,person("Satrajit", "Roychoudhury", email="satrajit.roychoudhury@novartis.com", role="ctb") + ,person("Lukas A.", "Widmer", email="lukas_andreas.widmer@novartis.com", role = "ctb", comment = c(ORCID = "0000-0003-1471-3493")) ,person("Trustees of", "Columbia University", role="cph", comment="R/stanmodels.R, configure, configure.win") ) Depends: @@ -71,5 +72,5 @@ Suggests: VignetteBuilder: knitr SystemRequirements: GNU make, pandoc (>= 1.12.3), pngquant, C++17 Encoding: UTF-8 -RoxygenNote: 7.3.1 +RoxygenNote: 7.3.2 Config/testthat/edition: 3 diff --git a/NEWS.md b/NEWS.md index 1e2a68e..7a17c25 100644 --- a/NEWS.md +++ b/NEWS.md @@ -1,3 +1,10 @@ +# RBesT 1.8-1 - January 20th, 2025 + +## Bugfixes + +* Fixed an issue with the `ess` function for beta and gamma mixtures + when used inside of `lapply` or `sapply`. + # RBesT 1.8-0 - January 8th, 2025 ## Enhancements diff --git a/R/mixess.R b/R/mixess.R index 78ef8fd..fbb90b2 100644 --- a/R/mixess.R +++ b/R/mixess.R @@ -177,7 +177,9 @@ weighted_lir_link <- function(mix, info, fisher_inverse, link) { ess.betaMix <- function(mix, method = c("elir", "moment", "morita"), ..., s = 100) { method <- match.arg(method) - assert_that(!("family" %in% names(match.call())), msg = "Argument family is only supported for normal mixtures.") + call_arg_names <- names(match.call()) + family_arg_is_set <- "family" %in% call_arg_names + assert_that(!family_arg_is_set, msg = "Argument family is only supported for normal mixtures.") if (method == "elir") { if (!test_numeric(mix[2, ], lower = 1, finite = TRUE, any.missing = FALSE) || @@ -270,7 +272,9 @@ binomialInfo <- function(r, theta, n) { ess.gammaMix <- function(mix, method = c("elir", "moment", "morita"), ..., s = 100, eps = 1E-4) { method <- match.arg(method) - assert_that(!("family" %in% names(match.call())), msg = "Argument family is only supported for normal mixtures.") + call_arg_names <- names(match.call()) + family_arg_is_set <- "family" %in% call_arg_names + assert_that(!family_arg_is_set, msg = "Argument family is only supported for normal mixtures.") lik <- likelihood(mix) diff --git a/R/sysdata.rda b/R/sysdata.rda index c26bf6312111ca74e8233d535c5eeba5b7269f4b..b28e66b5dd91392b658bbd5670967ac4f84b49d2 100644 GIT binary patch literal 6085 zcmV;$7dq%dT4*^jL0KkKS>uG@G5}95|NsC0|NsC0|NrzR*h;_e-~a#r|NsC0|Ns6R zfFMi)Knvg$e^3v(>*;1t*sWNq*=$;-wyfIQK{9ErX{4%+mTa`vYg=h!8Q=f_3RD^c zMvVXf3_uz*44P;F3{3#Y27u52XfiTnG|`ZHBTX6rG8mc~4FF|S|ICdvCTP$m83uqr z2+&59^%{CjHlTWd0002c0000000TybL7?>pfHJCo?FK+57@8U;CKCu66F`X5O%q0j zfB*mh00w{n0000002qK3RZssw(S*^EXbGkhCWeA!l>Moqo|DuW0000001SZ800000 z00000K+%CRdTFEvKxAkd0BFghKrse^plAW6ng9WjpwI@4MogGYgGPa%000#w1Wb&D zJrmO=l=5nOMxq*KnM2YZr=&E~Mi9xSjF>}BJvA6W4Kx4%VHiw@sKf(mZ8bJT&8a<9 zT%uJKRYgiFl|?Jk0X2QiYE}A zF))+elim~Flcy)XCs0p*PNbgzp9GyzIyX?H4K8__8`AI0o?*j zh}QOCVFip4On+(Zpm)HI5K;QgFnl2R75;Tv&>Qh~R>_=+o_dFGgZgr3!89TaAGzHD z35)E5vn?0E&XOvXI>!3Tcxb_wfYcM#s@4AAg6n=HS(P?_3k~_=-jk0nt1K>9)~QO3 zAifxWyF4;b+s&ERU2`jR*I0Cin_Pie&I1OJF&U=WLwAw)pVs!`iQNjb7bd4XJ9(j3 z36}i})mJkoW%Gwjj@jZzFv-V`W!nSqEOEf_s3`-|P_E=DvKnYoFbs`>rJ(&u0;-u741AQ<0{Ezfr>Ra9W`+o#SEX)f*P!V*=$vTfWEN_6 z8+d${LTVIYa`#ob5iam2yLR?cHuW@DMyMv8rx7MD{bl$)a|O}&wgnB=uKx(T6}13WM0ImY@|0hGk~v=PoK z07F2gN;TKk?z-Z39dX@yu;$DZC#tp-h-{}p>yGu-zQiju@iE7dxsAA#^+?5)tA04m zvmPRV(+dN<>TX`OR&7r|@xJFWwK4TKpBIU>|wlOx@uD;FGw_H&%M@-bw;ZOrpOV7-O{eqn$*LBXk|mf z+;_-Y^2CNn1WHgQ+B4qLx0$sTqktHC>==V7_t~NAJ0ZPxF^5&DHdNgwj;D^K~yr-dr(c`yNKGUMc~%1 zizbV+*G5@Xp49>bfOl!Z&kCoaI5sSxEg($}T!=`OC~ARDy^N4`7DO@!8wGupXx1-~ zR8O`%Z$A6mf(#&>E-`Dd2?aQt$kLT~gP2iV*16Qh*mJR-MOaOZVE#bbk1SURr^z4Z zz;`P$FLElRj`D4!&-Z*|XZ!P@3-iTP_PMdU5Vt%R3C@DZRUwWC)k(2L?Z7QoH2{Ik zITUW)PhXG>sstNL<1;oM>73%*k4Iu{souG{4DbQ?7TXsEn9Z@FGpm%$0vgH{^U=w#_vdoO7nQH9hvs6MBP z3~=4GmuEZDAIg(V=Q(MnK0H+R0fH48l;BE&t+=~{l!)+C!{!| ztD>z6l!dm*I^4V?vY9v3j@LSA&>I)cuyZr~ZKTr=u;S1`wi%K^sVo_QY`unndDfU6 zB`fdCsSmNRqGDZ)xziRsPeTdXzgZW8+OPseE?c(d2PJ{qZvOk7w`fb~H8ZLD9R-8X zC{>7BkWsrphTfZl&A1orXmlBw$!4tQm8(sN-B%3aMdq|UDGAotzfn$tZ}an(kw@R|mNDT-s)+REw`k;F5Gc112> zvbt`MGuhT^5x$0uJ^fA7$Yur_`!;QqN4_UowEhIA*9`SQ57aqRl z&#_V^Ho1e6@Y(@n#Ol=}Q3nLbxwLmyRBwtGuyr{~$t6jWScA?PgPu)yP|&8UHtr%j zI!g5`43wzmT>#ldjgN>!&bvaZS)!bV7NG3wu@BbrIErk*Z7PXyxq90(8X^=wa2j z&7eY>AXC4NyF5(l+*vGBb2w>AB^<31z@+dMLXD2*W?h9Si0WW+J&Oqt?RPTejbX{w zXWBJ#4+!wlYR@#;WQ&buj+cgqSejx`)gjwbu{}hlV3~T5^_w%g5%-MH$1wHCTWfZP zXI*nK&}~enH<4?xr0jDRaj|VN43ekK$mWbQ7&}Le#Ziio+>s{+y9YI0tQt;dW`?Sk zb#Ot&`XFOyR}gwv8GuHTvKieX*gb`I8SQN$Wjo5(uM7izHmP)5a|yZO=p)apWti|> z@86%H&2m|MMLbp(ebq6&@(pw|J0R_b^v(wKInuTq+2z(I>|IDTt+u(-@Au&7RoX?+ ze45!N(CZqA;a0Nfff?eckHC<4(qrIC-w zz$$)1rNOX*LQ|fhm2^8uC|#^495L^aX;o>GMPuqx{yh~06Gzp#`W?vfW}H*@FdCJ3 zCP3B6)a`o38guE%EUai+nSD1Ixv}Ma@l0xVRKd-i{X+&SCWZ3{K3G!O)ONzHcc6QW zo#O-2h1|taHv2^iw~9-a%;jL%EY|utjZ3+o6?UGT zT^KC>z(Tgex|OdbP#bS!dMF!@ro#QDquCyEsV9(x+hC%Xjp1hMTcV;HX8~j%vIn? z%*si2+ngZ>Pi!kPQP(^!ISw2N0!C7|leux8f#>W;GcZbO>=yS|ypcz%DFkK2B1p3R z6YMHFDdrU$nkV?z5FKGx=q zXNWbVgA72!1yp61ppNOQkk1U;>MJ5}(mo&deHK*AooU?!T9E?4!Bak+9mpZV8jhxut1Stdhxymy8I)uET2`Aq1&1l% z9d1xkw%SS$b8Sd9k$c5bJ2muFW)fG-JJXc5vXxS^vB|b40V#g~w6fe0&w|heS~QQE z#IB9*MhH7$8@6G0e~$2GUwb<&M{R22hInD;UVj)i+G(NOw28o%3Qzz!T71x>s0CFN zP`=eCO@O8jMDlkX`k%8Tu3@8R$ED4@*bzF9h44g*5SyJ->9Fn>2>0)KWVvR($(wqc ztPGUNk|mlZ?#8k@mbYC*M&!u$H)OXfK~&;!Lmd53>S}qWI(uQiZjlMNuwnPHIK!gV zt-lW^*tPap)cedLOWUac)2dBOcACl4IKA5#t-{4(9XX%Ru7U$5DUhY3O97Bl#ME#Y!fi9-Hy4$QScBlbdxnKB3+RGZnb@@u z7`GcOIC9X!Fb0%1&$)Qsc0RLTTEy(Gc92yh6aq;h6_pWpNq8%nrRXzpZl!uFPh_#6 zrgJUvJfJ*kVv-H>Mcz}q0!qon(l|jinU!O$ zeTQn!;IeO4D|4bGTd3*0R&pB;ppx;d?n((k7AZ7#q3t-rGI?;lI|bjHOGg27hNKek zBV01?(P&atrZCv=KUYg@6PFA`)3@(o)5(geX@P55ToX|ZK{V=aWsbxIOjRp__z|Kl zqI<~pH$=p%cd$FD3y(ALDI?rah51(5#VIZnAOdO0Ai5bm`R6NQ0Tb=LBOKV_+94LM zk(0e+KFgHf1k(!OB&*biY?!gKSh1Cf#3!GkAiDN;=-)UNfHj1c-DVnNy0YUzcu>M`_JY8{KNz zi>WFgcb81B2)VUZP?Q1}0T)u%oGd$gQN$qbb|##JEb@VG6cgxpQUWp469IXvDbUL# z7;zF>(k2Zit0j)gIKdbhfqn-{LO-ay;d6vH=Eiwn3FWwkva6Eh-rIm-b3|5*hTnig zf%)ydY3ZiHqna^oe$$aBd*3t|W6g6`Z#0~)E4!ydz9q5;Z$ls%Sm`fnhX=ux&WKD3 z&Te}X%8>_$LM<7?vn8Lk3x^%VP@8x{H_Qeg6*l6NqQEzR={`N3y*ktAEHP_aX>;mv zN4r>?f}jf;QwF8rl69U1>*)*)KGm!f5~k>HY`6^J9VMkEsLuN^k9%>aB?(T?aeaphZCh6_9F<_?jpSQ)b;R7Iy1lwee zT%Gk%xkTI=A310fX(TJ3RwkCpMB*J!b34E19wo%{+jKU~-pj!m>EYq}X5k5QL)3>; zSH4tfrv_&5F2t=>RiE4p#MXG~ScO6p$=!LWvoal5SO_N?fToR5J$$Vzr#iv{t39DM z5uh5`Y*j*A9B-v5O{6-3MGBFQ5py)E%To*x*|tlNupPt@u_ zR9`bWkmLpSJNLkD;cP1Pp%W8t@!x_cpx z%pKDQc?qU4m@9Iqlqkg~F7ht)b)f5-BZQCA^_NdqJ5jgOWw{;pm>9xd#L6Zo^B`0* z=BS~ivX4+hC6d$onw0BPv$Om*w5DEt<8ss%U5pip@!G zl5krGRW3MAYflJyuLK1yFVq}_%-YaKHqTbOF9D)x=SA*niCy#VvsR3*cNQQfeLX~L z0*!ZpWYgcu=1(EG!UL;#_+FQNEk#4m>a$KQj1AlXd~n@xqe9Y! zj<$;GjB$)9#4XD2bwr-fq#Z#>z`C}ioyxS)Cl%vedkv?u_(q2YBdB)&m(S_%2kW2- z<#=%{C+T=rNba#=wkMLOAH?QWwU?E%b1IOD5LH4*Px<6u=y`c#1tLinZCqZm;1cX^ zWu8<4ndzjGPoZi_B(eVOy{TqMTNq?wsDcyCkbx+b1w=@VItdaI%OqI^l2ItBB!XK} zAqiy(BomcH@97bxP`a=xNLxuG1-2m(tEI)M7K!$`?OQp11?0ccs_*{g#L>ZrM-SV; zO3lZ0j9bhEw_64HFDIXy{r>Jxy`9Xqs}-N088h?rFk-x#SLWY+)~*|;nYVy6T3QIo0cZxg&Ugb# z_Yebnv9QVkr%ig9Fo0R1U<(uil1U{j)!bcx0>Bf7P8GPoGe9}O3t7a^v;_70xC1kw z_d7;oF7|*-{B@(2;Ij literal 6088 zcmV;(7dPlaT4*^jL0KkKSwVv5gaA)2|NsC0|NsC0|NrzN*h;_e-~a#r|NsC0|Ns6R zfFQsIKnvg$e|-Trp88poHY-*twp%9FYih|LF{YZ0Wi}|4HFYg1PuwH1vZ|E zsiV|10000005kvq0001J02*ikRaF1<4H!)s011IIXlN!WPuiZ6o|DuK0000000000 z00000007VgD48`qG&Ma;nwS9@8f0X|zyK2hVLe7^Hl{-jPZ6LR5r&Nf!8F=v045Vm zCIVqF1xX|bGAZJlCYq-8PewAQsL8c9Q%}_f(?T@D4FDN5^#e?XfsoOri~%rAnqqo@ zG7Sx=%}>=ZJx8jG$|X@%R8*p=R8qYl6V;MPA`+q?Q#4bi0EiAgcoB{FY|Xn%Y|+66 zn9W+^Hp9-xO`n^TrBjP${m+f%c^V!qynUAvl~SBDoqJMT;cw)UNgO;Rw2T}QTh_I% z_T<{iyu5&%B2fZ+f>c0r#H0bkB}p-ph)(05CzF@tYi!)rqc0~ePTuHh5}$(}=dU+c zbOh7wYZ8gjlZ;Fx=!wyjq$j^8z$Z#if=-y7IXZ%MB=|(n_+$(3MvRhNQa2n!LQOSdL>{Q zj1(9GdqN@v<$KqsvivL}_~r@|>s9p>{dz53)R-aZ3Y7Ss{l>pxSq!HA7UW{&lWwi6 zEG|^m29UCX*kR&V;lV>rv}YV~ysWoQkyUr{h+dc&=sLFm%M-f3g`yP8>5T$Q^M*jpqe9NFGHBtv;I{qIi*jV`K~{ z_v_I91G6g@LMXYIWP*AZsby0jbVR@tPGi;vM!Jp%>4Y_n>t?N^yl*%!5sd zrVLfVUZ&r(f5lugPLj%(zMZd$R{<@%cP@n6lF?lHUY2GWM7+(YFl7z2D(1K1Ls=ms zC@^llAHgIcquB1vYu!IP(i^K!oZ6GrP{=2AcC5f}q28JIL@V%JlBgtF_ z6C1|S1AIaOkph`0)>~Ve^u**ky9^}PhXFk_g1pW)i&rU4dlAq4F14&|M14%>!NP0TF|-i@O^F*fi3N%UKF+bg z;SHUnP%y?ScHIP6GRtU&nD(i~J$Mq`i?xhbO-gW5?1V;K``l#LRBEthY=Ic8+ZwDz zl2`9U>!fGdhUkk{I1*$4L?Z$$pFPYcbj+lr`@pl7-m}=ty1fop)*F>u2lZHkIJZSr zk(cI%d{R*nNSaNO+ zyavT7D)vd}u1GrjC&wC5*pNVg4BZ#C`RQ^416qj!!UkvY<~jtWBT@=AYg~h*v7vA@ z$Sdu|pHRM0Q8vhNy!)?Z=dgNUw6(3kBoyOoBT7}G7GcG6I_C~9L!Er`D7tIw2k{2f zd1A1KK1lyA1i49(cac~CXLp=NjP8+`?kp@ui!Q6`5I zmjJa^%i{+zBnwt*cKm>45P*!WjLcShrf-W)yB%q?9@~46so*;Dn;_nz*|o4H4?+&g z3ppgCoa6n`J37}o9r*W8TT3#n6vqSw;Ts(K>t?iBIRFruY!IDS^_^uu&*xoOxQRnVKxnHFknAZp_1-lT(DBPRP+3J%0u zf_|25Bpr}yX;#nKIJ1-kp&764b4+K<+ww037Ml%tM-5~qD=~}Ifr48;sGqJAS5ks^ zX;XsvXbuPJ>8C0XQC|jIEXGOv0%gXk&HVn(+qZRZ<#gkw*)4w>6bk+xw0&kwI5k|Af zXSb}P5E(#r3ng7clln&8Mh<>_%GI^Pe9sZfH4e^IrfOs1rhPR&XP6X*DUV-o8)r@? zj=lNU6nTcasJ1-LXHkd7>DZWix|ySzWriBM>)9BMbPch_(%IHNHAiGYQ+RaIYRT$5 zooxW!$JDM}05a+fv8e^1h zmiONSfn_a`*da#x6Kvn2uNOL(7D~gpj4etax1|!k8Q&;MMf;hUaulM0#K7Ep6&fPi z=4FctgO8z|VAaMvA;UeYIKw554i%Of92y&1SxGw-2W2s0dPzi~GPI$~HPqNg-!D4c z!_qR_TXU~6s+n%1X!TbQaE)=>b-8a#QL(SB)zuy|TjmZa>(0xxA);s%WMjw4<$xm5rn9P_Qs|TVO@rEQwUNf_}b;+cy7YQ4U2AJwLCo# zI`XDfe)E0Tr|557mERFh6=k1wOYeFDEQV(!C6L~kz|MxbDSsUCWIOeoH4+Zu*~T+fZj+bdR#`R;X?kp6J5dbVXVIHb92sH+L$-7l3M+W=+mhh@?p-njwi0WX2mGY0Bk=*^_grQ z*ERFYn5w*kS(1v~mb5|+j|?gQjrZX^)u z#vx3ogh-TErd?*0ODqgQacNM%&vm#KSWM3_Xm*->jWyJxIE*RNDKK|>rm3@AnbwK` z?Z&qI@u_{Cs*oJX6AEn$^&u8ufTMjTtRqTQe%5s^5-)|&1X528?>=&*y`Z5O$C-wm zWHv;W)Q^8C2w8-%B4c7Pj3YA>0_m&eiFhCxgv)6`eUrjUI_zh%fKn-G+;brHQ3@fqKjbtPO)8aG5mRAus_$EUI&B@>)OJVbeU zH&`J+(pG?W_Fjmqz^)wqb+lJg5_QwL__ZPhL4t;PIy;ZUgf$#Cok?g+)IH9(5lj{T zy@h3?&{%Sv@@Z~@m8OyjH#E$IRTqE2W3OLLGc1CBO}<#iV(3#E`r4CPFA|4u3koT` z{THb4Q=drq)`hUS%(Miu5;8>Qf$`2`V-l^*HO6A!$e0p2k&4dZmd@p_kPK4atPNu`OTp``L z#S+z;yJls|F0e7t9!Qc{8M_+D+F09l5gU^P*wu2}r3DWYf*9xPL#dhOnCb0?{kDv4 zMS|};iO4M$D`nUC?$wU2wx?WJ@m$R(2AxrATeik6%+j{&TBeHyXKd3upDD-;m!?9N zi&qqcVS>VE46qp4D!Q&JtR?|$7nRrt-OmP}09v_quA;KgM0b*@C)phDEwyH};@J?ggZMakIpPw|D)5JniB*x5dn8Rr^<2M&ni&(|pcJ|E* zTo%y$uCcLNA2Dt=R^iJ+=)NjY+dk!^b;$bWze=!nRy!ytk^urqAr*B%vyfg5GcvUX zYAw_&JxJ@06P9NY6y}Q6wF?yBIo)1O^>YPz8iWrzJRViH=by)e{LutYZjuiMjldEW z$}U5E5?7fuyBAhwEAI;ybzrAK;60j*i-lhzF7gR28m?p!gk_3N9jJRvu8f{wPR_Y^ z=CzZ6xy4cmcrmOgcj&MvCc@Y2Z=>0+vWQw%OOxEGM;zX75x2c}fws+XnP zXd|%Zrw#5iWX05#54p>xS4Ldgw3I~v$G~OOt%nlsUc@j6J6nNRSAkcf`zFF;z}{n@aVH$rIbB-PVu9NU>_ z3oA4QVq>!F!W+)ZWy(uM6^#607RPKc5O2BVLi{ z=qxd7S7~$VY%cb(kBbTjnif+by){iCO(W1JjFVtUHNvx^(%d?qoKM*eikiAoSI4=S zQHoU7j#?ONgVD-|GdIcJ4h9q|Ia)mI1rT4N-4xwdbj|I6XmzW?u_fi)Vp6?D=|6LW zCBfyv2nn*t+POQapKL_h8QnQ&6DcF=pH3!{$^_CIPj;sK40sn4&27)yFN%|TDx-ge zqN&io2mT(V!aGXmvtb+iRsMO=z_OiWGr<5piKP6ruWVH3DB-SIOYDWOAF1EG9H; zdU1nw#U)Cpb$DsXKo&RxIzYY%loOih;YjQ3jXM;|UFE~M3(8_rVLRwrmT_sO^(taa z$(ym-*R?bijMbI=?_jN|KGLFi2CzHm7!rw16cJ?b?)5KNoJm)- z3>&=Sy$r@yV!(mEleZ|zOHevZ<%Uv-J@)K4%oz~}gkHtXGsVu7$m2?o0cx6E3G{Lq6o;pusn{Zqp_Wtw%Yj_EKSILaISe!g zG=r`|aa3B(Ub)gd_R7B>~?M)5fg4){FErDoOWo4uSTk+fp`=Y@q*dT zNgi>k&P_5ag#&#+pe3f)$E?Zw3TyHF8^@*H_hTU)$EEak4ARdOv$A;6n5eyikn3Ed z+Dl@OR*gDk@TYg*0SILzr^Znbj)Mm1KD?PYT&{qD876 z{hakt));HOy)Sq#%_PEDtj3ZMqZFLG_`77%gN^8p7#^kC?w+o8V{4|nGCRz!LI}Ny zf|%aSfl$kwP(MQ@9*~F0joI9jo;GlhoKWhnVG`5`kuvS%b$8f}c%s)TynEIqlUAn) zrG9YP7HSK0a{|}8=xL&JQfI@>SDqLvTAi%cnt#wIxlNX5qIvky*cu^ zyjFOXXzL>w6l%B&C7S+rv2g6d5E(m%zT`Wx(o?wZs|62dgC`t{(A$eE6?}D1eBIy& z#SONKGpQI>>}aZ(M+m}<0@|MD(m~A$J=7E?1+_60ty89WxGM^(*EXKFLNl~Bpn^%z zNkHdgH^SS)^576Jrzm%0boI+kVOxC^wlUW$k9zlIw*Lk?q9jJ01c?c2NVFD` ziA6~y66%QvOHh(QKkXB=q(+TG^P5nSPyJo3#JiQJZCk$| zZGf5>w)AlJ*szk@+{^Z5@9+}rcracEi1gQq+8xe_f(<7PAON171ZI{1DN8uA%K#J# zAft)^6lNSv`(bw&Zhw!&ZjT0?Dct#eSLpN zo!@L*CGCzC^qBTB`al~F(=oAl0_(j%6EBCFXD~6c0B2@t=;!jzPzHN97M`c_;ynI) zE0O2%Vd?OAlTz7NxEyP5b3jY}rt;EBB#jN#+REPPfS6{wIj#2lS0^{E<#bJG13{&p zjDQyaZ`R&`Xy0*wI2$`GfEsk))W?JZ+4cct0Hl&hOK&yrU;wZQ#U|C+U>SfK02eb! zop1@h=3oq59)n4VlVX5O?7Gp_#+k^HMIuQQm@i-eUjZ)VfHS}{I1JWNFJ@sN