From e42828dcd16e0ba03aa66fa3d5be8d9a4bc87060 Mon Sep 17 00:00:00 2001 From: John Helveston Date: Mon, 7 Dec 2020 15:02:52 -0500 Subject: [PATCH] updated simulations vignette --- NEWS.md | 1 + README.Rmd | 4 ++-- README.md | 8 ++++---- _pkgdown.yml | 1 + inst/example/simulations.R | 26 ++++++++++++++++---------- inst/extdata/sim_mnl_pref.Rds | Bin 282 -> 279 bytes inst/extdata/sim_mnl_wtp.Rds | Bin 278 -> 278 bytes inst/extdata/sim_mxl_pref.Rds | Bin 278 -> 279 bytes inst/extdata/sim_mxl_wtp.Rds | Bin 277 -> 278 bytes man/{figures => figs}/logitr-hex.ai | 0 man/{figures => figs}/logitr-hex.png | Bin next_release.md | 10 +--------- vignettes/figs/unnamed-chunk-11-1.png | Bin 0 -> 91668 bytes vignettes/mnl_models_weighted.Rmd | 4 ++-- vignettes/simulations.Rmd | 26 ++++++++++++++++++++++++++ 15 files changed, 53 insertions(+), 27 deletions(-) rename man/{figures => figs}/logitr-hex.ai (100%) rename man/{figures => figs}/logitr-hex.png (100%) create mode 100644 vignettes/figs/unnamed-chunk-11-1.png diff --git a/NEWS.md b/NEWS.md index b702ddad..97be8aaa 100644 --- a/NEWS.md +++ b/NEWS.md @@ -4,6 +4,7 @@ ## Summary of larger updates: - Added support for auto creating interactions amongst variables +- exported `getCoefTable()` function ## Summary of smaller updates: diff --git a/README.Rmd b/README.Rmd index 57b85563..2e9d1a7c 100644 --- a/README.Rmd +++ b/README.Rmd @@ -10,12 +10,12 @@ knitr::opts_chunk$set( warning = FALSE, message = FALSE, comment = "#>", - fig.path = "man/figures/", + fig.path = "man/figs/", fig.retina = 3 ) ``` -# logitr +# logitr [![Lifecycle: diff --git a/README.md b/README.md index d947e1d3..42c36da3 100644 --- a/README.md +++ b/README.md @@ -1,7 +1,7 @@ -# logitr +# logitr @@ -88,8 +88,8 @@ citation("logitr") #> #> To cite logitr in publications use: #> -#> John Paul Helveston. logitr: Random utility logit models with -#> preference and willingness to pay space parameterizations (2020) +#> John Paul Helveston (2020). logitr: Random utility logit models with +#> preference and willingness to pay space parameterizations. #> #> A BibTeX entry for LaTeX users is #> @@ -97,7 +97,7 @@ citation("logitr") #> title = {logitr: Random utility logit models with preference and willingness to pay space parameterizations}, #> author = {John Paul Helveston}, #> year = {2020}, -#> note = {R package version 0.0.4}, +#> note = {R package version 0.0.5}, #> url = {https://jhelvy.github.io/logitr/}, #> } ``` diff --git a/_pkgdown.yml b/_pkgdown.yml index 95470f18..4fcb67d3 100644 --- a/_pkgdown.yml +++ b/_pkgdown.yml @@ -54,6 +54,7 @@ reference: contents: - summary.logitr - coef.logitr + - getCoefTable - statusCodes - title: "Computing and Comparing WTP" desc: "Functions for computing and comparing WTP from estimated models." diff --git a/inst/example/simulations.R b/inst/example/simulations.R index d5387fc4..23888e7c 100644 --- a/inst/example/simulations.R +++ b/inst/example/simulations.R @@ -49,16 +49,6 @@ sim_mxl_pref sim_mxl_wtp <- simulateShares(mxl_wtp, alts, priceName = 'price') sim_mxl_wtp -# Plot simulation results from preference space MNL model: -library(ggplot2) -sim_mnl_pref$alt <- row.names(sim_mnl_pref) -ggplot(sim_mnl_pref, aes(x = alt, y = share_mean)) + - geom_bar(stat = 'identity', width = 0.7, fill = "dodgerblue") + - geom_errorbar(aes(ymin = share_low, ymax = share_high), width = 0.2) + - scale_y_continuous(limits = c(0, 1)) + - labs(x = 'Alternative', y = 'Expected Share') + - theme_bw() - # Save results saveRDS(sim_mnl_pref, here::here('inst', 'extdata', 'sim_mnl_pref.Rds')) @@ -68,3 +58,19 @@ saveRDS(sim_mxl_pref, here::here('inst', 'extdata', 'sim_mxl_pref.Rds')) saveRDS(sim_mxl_wtp, here::here('inst', 'extdata', 'sim_mxl_wtp.Rds')) + +# Plot simulation results from each model: +library(ggplot2) + +sims <- rbind(sim_mnl_pref, sim_mnl_wtp, sim_mxl_pref, sim_mxl_wtp) +sims$model <- c(rep("mnl_pref", 4), rep("mnl_wtp", 4), + rep("mxl_pref", 4), rep("mxl_wtp", 4)) +sims$alt <- rep(row.names(alts), 4) + +ggplot(sims, aes(x = alt, y = share_mean, fill = model)) + + geom_bar(stat = 'identity', width = 0.7, position = "dodge") + + geom_errorbar(aes(ymin = share_low, ymax = share_high), + width = 0.2, position = position_dodge(width = 0.7)) + + scale_y_continuous(limits = c(0, 1)) + + labs(x = 'Alternative', y = 'Expected Share') + + theme_bw() diff --git a/inst/extdata/sim_mnl_pref.Rds b/inst/extdata/sim_mnl_pref.Rds index 512c919cbc219f3ef7d56ad437841547840b7dab..5b93fd27fd83c78beb2e611163cf8340da42fab5 100644 GIT binary patch delta 252 zcmVT{pHus^Zv^S+H`A@=9Lv$B5hItWz41XRZebQmYdQF)2EsddFbh5*#HAVIF;jKreU z_}tXQJg5jKL?kD_992n1W_kviiJV3G<$B1bvOp!c9dk;o6jBoN^78XggfcR767y0} zgvwKa2A4pEc)&uH`2{(NnI-7%PtHj!hS>xYVYfec?_!w13=_DzrxJAu$!urr7q#dm zls=KMKe6oA>QgEw?a#kA**!~T5>Sl**nSoUP9V+gm{Ve%;#{i&sXl_tEw_NT>TH}2js3#fz%sE!fnFiw!8@)C1Xi-8OQsB1xj zT*VoQMXB+*sfl?|5l)CmPJTJ6l8nss3^Wrti}K6$kWFQQN=R@!=9E||q$K9$<>#RY zWn|_g=B1zrm8Sv?E`bX1fQ2ga3vv=OOVHh)oRe4#vk4~1m6BMJsFzj*bOez14*+68 JQ|^iZ008a3WF7zj delta 207 zcmV;=05JcS0+s@hPk*>yPF_ijA{Y9kf66d|A}aKvnxAUl%{$RR0{V?x{3~ z!t9l;_KU3;TWZxc>`!g?kzLN#X@6SQDJReIJ5UJ|P#q)CVVodG05AWR0+#}iPk(oRc)+y&oc)4t>+YP}k#2uT>fkI#<-PX%_w#<*eMtqb?nynP z$4m|n`=gt6T$%P?vp+ubm&{@3zxHSMy_6MOm=09J1XRZebQmYdQF)2Esl`Bs0MxZ0 zL9XJA#G=&r+|dIp+_oJINNddQ}-KuIOI9dk;o6jBoN^78Xg zgfcR767y0}gvwKa2A4pEc)&uH`2{(NnI-7%PtHj!hS>xYSr2$#mfP=twD;01-8{IuCq;8U zbYG3KKU$>vwruZl`{Q#L{?^mTw?C(MM`(JZ5l{&eP#q)CVVodGCS`1_eKwS$G z`%EUnR_NP0hItT10&F3oFGT#CFZ6U0~rEP*MbDO ziZc?6QsZ+|6Z4=VoDh+m{Bl$!8JX!BXeM$N<(KOro5})};74}MDX~&WNzBX3&qEQ) z$jnL1OFNBiQT<*AB8$rjzVF% zV;@7IP&ZtC_~1XJMj|4zCL*FD7N!=~a+cb_WGJ0Z#J}G0xsgQx7*58sFA7awJAC zg0QMA&zI)0yEHcK!LO>aBpRl)47r?OLiakF83Y+*))8H0KY&Z|+2(M}5m28&0h z3O*}^ye7mT4o_t|$2wx;9pEFt@amcOyK~3hE8Zh#P*0S^YO zNy)oIPDE)%p_4U{fa?>0u@o;~Z;$bK=QV}UjLO9govFO*!c{Df^M`f2Ds2Q^UmEF; zF^;tP9F>?n`&l=fLX7FZNY=94ad-d|2hr3su9HYpi$8zy9m+S@69ba%tWI)-O{+IedASqRhdz1Yy8r$kjxecayPxv?=ha*)w4PH=?O zm~33kjrq+-JXMxab+fpsCW({hxY{Wsh-Y^^?!P;{mE~vco~C>(>E4T%GQCU%GqcBlQQeVyY`sdD1 zmh-RhW`{L=DXj=P`1Vcn>2`i&HvL0y!6(F-WPAR}#BF-ZNuFxDXN6oDxiE!jNh#0t zJZ;io729}i#fow1kuSg0M&O0avkw<{y2VXzYw&Y*?iascyB`;Edw?gq=$+Of8@;5O zTvsw9Ual8)J>Bk6^^(1>ni*=-j9m9}M@XXz6&trTw`ccQHyHMeKmBkUQ(Vrou2K13 zlqp#`XR2i5AWvp?Y?OGZC9L{+HTkO5@|Yt29oLZBu*tMFb$(?FX1qJep+|R&w-rdG zU*1Yyd%EVG%~*29lBoLf zW?-&K9mvZ+n6?a}BUiAH;x4^$_Qhq9hcALSA8h6;uD-t%EHi%Jh(L}{?rdU4`we{W ztT_93W|+HsX4;z(uai4O!*{`czu)fH+}qvuNb*$?{U6<+)wv5&vvd$OU(g8ftTrrn=*UN4Zy@^=9^c zA%A_>8TUu~=G;cv4y&W|l;H4YJyBKHrBWk#d zU@N%wl*II5o2!PM6r6?d!eP^o|h+ezPSgHJRziG>5y&Wm&Lm}^*Tv<;#KI<^6m~~p{NzEBk8q9rn2Y22U%Uqk2NUyJ4E8VHkrd>P zG7Tn^w($DtkNojTVbp}p&)4wB%UZ!_2;5WdS57yeT`LwgcDXa&FJ4Z*kBONsZ_Hgn z`jadDd-)h|n|>yUzh3r`7Z*n5il2(-5<$OMO4KS38hIjLlP-P&yIdi0cugO@+_A8< z1dspXWd%A6jBa|4$ib877weAtK8ar5?&b6}lql6$ueJy(^ou2M-@brOH{e>7Vqu|4 z?Us+({ubX)Bn!%MJ8b2|MTjFnXnWzUDlgZ zRP$6Gq)5YRovGjwG<~#F7rC|%b{D+`Ho8BS@+CU_WW}t6Bs|s+x%ZYT*-Y=rA(x+n z#psc5{4(AFy*6xB;J(G!@H`|;U~!*}c1kVf{ltZ(KzU8~>x=Yg7mo9<@CGx z;^@!LD>e)arQhH34*~s_%TEp@7|7pp`Tt%nwyDlN9i7;Ws=18(`)h65)2*@E-Rauz zRrAbK3IyzWln&RX+bY`=3B2Q zz7_f6%|WN(`pu>mnRnZDsaH-!SX>-@G#Rxt7AWIIV#uUX%>Cj559zC$=7E9|D*W3(x3rmSO$1u@)`*-AEvzA;*$ zNLK6HE2`mE36@hOPMgx^+#^2J6;mxAD+e5g)Ey6$7|TZe`Mjb9oM`pL9pJZG3+K2zdt?b}2iIA)F`_NOTPStyAF5iR=!KtTM$gF4f zw8lNhd=j63Nm#pV+6{gnZk>-7y#VuI--D^sO?xscB$F0h7M+H}45Gyr{3XuV&vkD^ zE{ADW&zB#;XN=gE=ht)=SoucnS2Ubjw1AV5CY1-V6>rZRQ8j1?CyJAf3McF~PD**R zPc$DC5_V1drg8=>JlXpA&(AZ+%7gdxl-ia3Z&F+%tGV`|olEWS^qP2c_a(=gmAhkzz2;U;u{D|95ueSK`G-IGTlNX(z~3^6>QtWMIvJ{FEyl3; zIdyBUCoAy$?&f05oMP{EU!H~aquy+zTGxJ~*_}Fi3+roc2Yc5ZuT9Y!Lk&5;pcAn?)V3wchol=1nUK&YllC05 zn$lME1N&Aa!eJ#aTrl`lU^BTxN_WIbD!Et~KA=}0lH#mq460auPDvQ0Q1f)CyZh~$ z`Fe+<%mxEUjzKU!S*bK-mZ7t+-{j;G^Mdh^31%jTRQ7&``BCemQIxVZn9JU6Blrn+ z{Zi`G!ar~?HLn@CvwZ8@;5Wl|Vx_iv)((Rk?>7B!j*wl~UY zwT6?Wmy>Jfq8H4?+9?i@aD+VYTd6N-V_3Sz`Nn67u~uvH+8&%XM@;1lt4 z!USCpfAR+9gwf#^u%OnLf8^!Z@hin*DC|+xva3!fr?pOcR6F$*mWuCB_hcF5gX%F5 ziKR+~vC^CLia0pTwSW0oX#VGEV0kwMpGnDR&vkQVj?Q}MHOs5 z>51=euTuvdlu7LA=4q5X)@v_BlB1;p;XMt)vMa&nE<0-{o_O9ry}M2?t2SIaX;R|< z%Jnqmt8*A`(y}$w9bnWxIe=BIb&ZW`)x(Q8rUhFQ>zaQDt>PQ9fZ{vV#vA6rxfV+d z5Zt-ila?)~k~*!a0s;=+874z-Y)in{&eOzv7`XBD7)`7oWCR3VY34?|V13_QiZCxU z%A`Q7Xl5NbTGhasxdr-s*j!ijBmf|?z}b}=31+^C{^B4w$t{L{8$~@D=eu*`4ymqzYq(}LQuU-7PQ^J8 zL>X&ne$cQ{d=#;?;PzO$?R4u>1uo08voieGFIbjiXwRnP;VqVIYB=sK*Hq?PjlS8{ z*~oof3^-SPDw56IG8%JPBP$I8BxEhrP~M?hyf-Nc97 z(Hu!XeIzJzv%zvOJw#Bp*eqb|{aNN3pa>bxe>v5}lH7fIsa$Z{mUeEtDc~q`kVs z(MM(kuohJE##fbp*xG-e({HlEs1T~$xSyZ3Kd#c41$MtAgZLM3VI?GR*YioCLw7mg z&U@jk8UGf`e+luoSbiE0{1(e^vHWeLe%qFxdWqj+`7M^e{Q7SP{Zm!=EtcP6`72ES z9p(Ns&-nkBSa=P(L*f9DS*pDfldW_I19P1cNnH!5Fe<0UZtW*P zE?OVufMAk2#9MP&O?BZ=FwBT~U3Uv`xL3i#HbT@KoScBKC1Er-ik~1n6j0AK?Gf)7 z)rU02xdn1C&9x?XUZ+Ew$;Ondi_>@>wwTvZeiR61*Z#D;_|fXj`|~>xhG!w% zLqNaf@>9F?TVK!sDZlmQx4!&a0>AC<&(YD}a``Qnzg_Kbi}`PR@t@f9Z@K)I%fBV? z+hYFTU@>i_<;C9<;V`NfvP;GY3P?urauoflDrCh1KvDlp;)!s~5MmfR zCnu-VrXNcwiL}uEnk8_+JfrM;6MI;Vq?wh{DJ3H~nWS3u>#LfLZ})92k5=AkDfW&GKfe-v_;q0KOqtOZ0JLXEx0{-&U)ijEpP1%0O9dJt`mekW<3pX7VUR3Dq%`b1dnq8Zm8=WaqalKySRVi zamu?6ckn@cuXF(+s{6vQU%fGtlomac+G@S*ot!&3l2zv0XI?Z!6Z&X%!ZWLp*W|W*XcWsS!pZMpa`p`v z9>v&I&0MbC@LRI@m!f3wjU0~)jywiz4O`oUb8nl=5m1PC@o6jV=U5g zqQLk{utmv=T+@f)`bm0oBRFSzMfUS5#WX>U*5cBsRgH%Si-%sV*JgdF{qM@FP75r1 zz?cqJlGn&8D$mMv*>oXYdQ8>G|pIr3cV3^!U1q? zCbaPf{i$5~fkB|tjajQUodhIR`#u*@X6;aQ`xEEQvhR8Js-vfkWS6<47{^J1yfhk8s0Y+IV#mS8 zgfralD0qJMir4aJ&6Ul#sC^;}5`o9IZ+Y#Ew}1zRhN<2Qs zs_~PfP6gW&5tYDb?o(SHZ>%;i6L^{P@bpH@)|>_IPj=s)i1U9yt&(|SJog8R9Aw17RnqPT+lz-vnbxD`fY>A!K48-% z5V)^0XIeNNCzi8Vx@)dYt7Kh)-ju$>1BjdMc&x2{t0;v_5OSo<))fQNN;Xv93c5UQ zb+wv!VEy5x!!f)IZ`-$k%Iev=8KjnuzPTrX`j0b8!fq&gf29&%=ww5RSq$M;S=5M4 z%82%~I3(VP>$;UKYJ?Q(d9y0zQ)uQ_{ z%GB(AsONhTEfdXShw$2wWTMebaKTq%W5Dh`EGvlO-ZiNU zqVtrnnrH`~HK>-M5aowo;I`Vp73hR;`5p>2NCS}m2a}E=t@w-%Nr*mAW$R`^Ep5an zqdVyP?SzZqFspni1K!`2SnvxQ>0ngH@~0b_^7`LI_CmlTaxzB2-fvjZ8BfuQK|LhH zF|~Cz6yk+US71Ps2m2x+gW%qJnwwPJg~wCx5gPJ;I|Hy+V9StN1@dFpSu4Hq@vD^; z+4R3-3mq&6i0JB$T5}){Zg`<+kvVd)`MpR2;e`QU8`>LhVB zOaF2VVPt0L)_8#dB9^W(cnmaiS1G^CHEnTH5*u9omX*yJy0feGXQjERiuYFPhvCFy z*M~6SQr~%jy*rmB)2$b(Uh#M~MbQy|*`=D)yt+sK$Z@=rHRKLzt+~Gtp#YH9RCsvN zvK~|P{>^&)%tLArMU>ac?50$l27P(MPjH#B4PspE9f5j@6A!Nb3xQC0pMAHs_ciOW zZ_VIrCz$y}{{|}lLHJ4F^XQbJ`dpOv%EjC3)Oliaph>3B=k|G#T;&+p*<7SRLadvX zgrUNwNh$0o2IgudsevQMPGTbjNN(s%sU$N=&QK%DXgYtO>(v91J2`W~*qPxj>zmSl zb8G%kl>ZEA))}>qDFuH#E0PL|&B~~pd08@=zu?;0WW_WM?fvyG7U!x-v)CRTf?)e1 zIFoI}W*KT~jmzKx@Q#z)%5F`445)x6j`MLcF&9Uyy%BP-A^gd}JNl2G`j_6WUwy?K zn^+UwB^IO8%;C87qe3*wHc@`SC#wXw!m3+_9WQavw`a^efc6YK-fM-o(Mpn1HVR^9 zUM`2z=#z-{M0ktz@jH3w*NOERQULkM31+lM4w@Dq%j%depCbPI6#l1h|2mg;;H_r* z`+v3Gcz6phImrNeCrywb4TS2bh@xMtcjf`SCH}}q)qjclpAR~qcK6xF=G-sV8wGY` z##YPL@1NK5_cQN!!86#IFZ^n~a`2WYycg!b*n|gHfN9#;p#INg{_7b#dU#8lM+?(0 zHem(`;vEGAzgn*}xW?Z~^D`O#zpFHKD;QXR&jR?ZKtDZ%zpc*C)D$s${});v-Y%mu z^oBIyrWPbD!-?_H;tTw9>6I0<5vFf0804@)Ly#rFS~Lx28P&wYuSD>~{O1W0kbcpS z!Wnd(g`~LrYt16sYf~*^Nl24q9}Nm=*{Og?FLZ4*VrVZ8#n_TAXM!O?*{Pz0BvRWL$C zn+nJPEGqJxDJq5-QzJU~gNmOUKcEpQZwien8ji z2}Mm8J4aY5F2>+UwO?T0sMA)RV(FQUew(JHj}-e?Z5J3d?pn^Z#|o>vdVPI%(7-h$ zDL4p?SLx6o?1ePHo)`ge_8g#*X`f5)T@LR4c|ze)2yHpvC0I&=*4P@;cYyv?0a-Gs z>*4Oul}QPQk#zum*$`-}AHx1&PTg+}{Ud zG)7|eH4hIq%e>hdxj!~5WF0-QV)atQ7uvHNrrx00;5nval;lSrCbIr9@iv^r!oR0~{r)tlEKhp7-{e<27Bt5!C*qo`%iVcM;q zOy6oST>A0F(`~DQeN#z;jq}+9g0>L&?WO zHkQfD1QKZ@&OP6(;&4^z>&G4~Zw?{Nu(6*5O%RY$+g{OTSWzuE!`3gW3HsMG-g+uL zE^T#2SDZb!UtcT0Eys7iXFs9@0!0|sGc*tx1fNpW|oAF!ysY3@FdDQqAT_cLABofJTG+XE|plbZ?yAKFTLotC`3GCGnujXg zFSFuCB-6QP6Z5N;rLPO~0uYBwC9_PbYIW%=x z$nfLK3ZxGa$t7tgM;hZI%vH90ni%a*erU@of(_W z@&<@U*vfoW%Ut(y>iNE!WE_xq#SVuzsE$uFjY;*UFLaKiCVo1r?ceB;! zG$nxz+Q^GE*wf>dz41HjeqP2Ez*gzY_C-96Bagx$p<3CnnyK-v-&qBi15iD?%+YA+ zmo;Cyo?<~0f4CA_mQ8OfVV*7W?knazv3qqhgEZI6w#$C1dmkP_n``B2lMsb@@$y$> zZUZ4paEUoU9-gsh+*Tuoo*})aD#cu`9(_M|##$Z3Th%PCr$<)28G1(+{Ut)@{Kj>* z!Kmxs&LeC@ASJ!7(Odp;eK2W3(Vo7O$FTWbrcDER)vLtmtg?|-KMI$%7VB!=Zu!FM zSb~TMIt3NlPE9i7U=5SP!_}BWex$I-;{rXR-q6X)!;_kAgZLfpu-N>UB-cOW&jt$* zU}yK8@d?#*edL2*$6yZ7qAe8<^ToWiUcja+xIQ!Gc3-MW{CGXZUsbYxROF;YmxEa% za~Dk6i5FIL+XiQr@_E&&r+v$%>Yt&&D;bacF--|owPP8!ecImSbx$H%O{(`{2aEDm z9bY!J(+`Y%^51%qHj)eFE6uFOsihZW5a={o z;7#Gm&Ex+?(OF+PoGCSk6Tp{i(0A-xjmRQ%c3D}=C@+8QVWiOpG|Up3KM$Y_htAe1 zz+rpogQ#U~JnV7H#mgLZaojrO73J6UM;~--HH^df1jc@fvU^qI1OS8Wf1a^9>mjdB zc~D7A@uH`f>7e>MSLyQ0&OwE?6)R%HJ&73OTtk8`(cFIqE@Lv3f=3vb6}wwIZBJCO zXJsSZ(F$g5H4FUf?Q&Gdpb=kiltGX>i9;1Y{zWHM4*u9UjQAvl{%;$FR)Q7T@4C1lVfMwl4U%3xRRNSj6; z&m55?m&z-VmN&P0np}yYaX9cMk(KyoRI-ZHX37kv)#DYguI)2fTHS!wzVtUR|?^UdpRJ*I_1}KW!@(9 zi$%xrE+{Ezwms!modqaaQo2USRR}IOaFm@72}hK#V1v-|49>9$N=anp`A<56g=w=+)m?N{>id@SnW*E%yG#LOxD92YQQ zNaomp3M?HZ8dLNv!%)|c-G9Xo(F+c^#%$d z7_B@ldf|1XA3l>AGhx7D%RdjC@VgvLsUvAH5FjB({@IrEPAk&^UXtxBX#cdfAym}o zbl8pq<=$vwZZbZ+lOjep7otEHjd#8!X)g07^8Mq)m^o(+Lx)iKq^y^#${7T&TDE++ z?Rk6;hBPS1m(_C%WZBkDgs7S;+o_vhPXCO%w;bT6w!jqqXd-MF(0Yl+MQHt&7r4Th zNck};|DsJBM;M>o%!AIyYLm*{@<>JHzco&a21`fjL7_a`>NFTBx|B;0Ht_1zNGVF0 zsFWsKvK?m%4hSFwOF{8n)2>54Hh%8xw-VtZMvv|sldkk?V%5uzdHvQBnVJoah!gjO zu^U9@FTbd?u?!l|*y8{9&;7LkL;6>@QJh-k?wR_kIN5bLhoe8hmYF^NJm>iI6esrb z=34c9Y^VkfR1uVyn;AUhZ@yE5CD@F9>bMdtMDjfvRbl|8b^w)UzrdLc^KfS~bJX1PhAL-*Es$P&MWSFYx@l%uu>+YZl4WZxVj1NONQrH(;hY#X;U$DNrwKW`ifTq9GzQ(xlJ^hN*ZWMPVr! zE~g8#f-?yp-(KzA!H>_OucrN}aM3ZXC_V=~ z80bg0dx`l2ln6e};qs&^?|?&mmQQTd#uQh7-O3+-A?u%>?hh^ZP{FODIWfAJ@7IU8 z@1@{0XBvJ8U9};Jsc!Rpe&PB9U$#%o#Y+`q&iy>op8_jjP@$%Msw07c!9HX~OeRQS zok^1LdHvO|<)8DbdkjFyA8$eKbbe&9v(5aRrn8k#nI?<5Rdo(Tt-0x+blcuB2~ILL z+WPG??F;N95w9afHE;2mgmWJDr^+u?+&OPX7!B3GiWjL9agP+9_zt1`c_!Xgu=)(` zbmq|$*}eRS5&F|AWQPem<8>$8`&fth8ObM+309uL-$za+AxV33!9+Nfpv7#H+kuU8 z-*o~%kNdGq$BdC2c?#;G_f4V1<|#%4FoOcs1HXNrU8o63*Ri|ADOta-ZG;S%90{Ip zy7`3?Gb$0du6?8+QZPCXL=^5H?niLmEA8#?L(Z+7ds5f=pbNaHnqxrqZpj~4g0Qd%(UgC0Qv_n^;l*JfFR(Bc!?_ zkr)bvF5C5*%>0Gl>P%zlH)h4&F=M7*x(g`7%}8s~E`lg1UL&(d5PRw~8}ile-Vo&Y zklc{yRJcIW##Lmilfj(!Q)QypxeOI1JI%Cny@sV4w~j5(RDuLKOt!;X-8JR>o+n{L(f z-29rp=yovahVx!IvPYj-QR zl3+p+?`p(&(1toff}(QG)}I9wDXGtwA+&nZsXAbLb&o5%f?#eTh19I2@w2m)qf{8$P6Z4m~P~9i>#8`$8KFM|Q`w7h^ zHc5FqJBM`(d&U0HJRA%x0%6q8{_zk#G1d^FNXBJus&Vl4DQ$4JX-{kMo|Q{UoW8q} zY=e!FsM|om?4=+09`0D>sTw85oMT0{({CU_tSp}R9;W+5`@r#T>T?t(Hdc{4^Kzl( za0R=R9M`3n{{GT!1(DC-1xojS{`2gEurQ`p+1`rZ+)&xjM*EuB{!!+mutID_jDbbj zE_Z1OI?YmkjP&esD?Bt`kt|&2&{Hq`M7_wqAT`b?D6wa$$&U3Dx=a1*@-JXwoN%B1 z>d+{)!cP6-$v>C&pUV}#gN^m@MAaJ_32{lu7mS*v{rpe$pTS^Q>6wnC(9=y-93WQ@ z#w)K_GzAkMB)%~j3}w*_fpbx~V_Q&pfu~=naCNM9luYEv&xDU$t^>xV#XYurx?4-b ztr@za1t8l3godUEq{9zo{jtbuhH0Bq)%jXNvrGpHd$PhG67c z4&UUln+ZBe&imkw{{^~2{B9WO6$^v0HnLM_zW5R|8E?jZTIo>KjD1L=i*{9iyl_*N zQVSy;NEsTOrrobDPXkc*5HM3Q|6$6X@Abd-MB2lh$zv)K%kF>dGAqoSQo18%P5#UG zdAuXQ?mpuKfmkq54~`zi4E*^i`Fr!Q$vxbUmVfy8QAG6hl^Z`jQ2%->8Agck$?$+D z&;PadpB{vV&rlfYFJ#Bae){UaK6`_jAf1Nb=^gjySi$AN42n;T1&3L~wzeF*KEfz8 z`OF!E)921jJ`Q>I46BUaPoqn@{o2QdVs)9_;&r~so8%Z+6%ysE2Kj;|2a(@Gm{c8a zEM;R;5)fiWyNTA-xuu=Nk)`!OVHy!H?YW7RyxLfH_{1J)twW3Bc(BQS)}kQ=QS(=4 zuoF?;Egwb4VS4cEQkPw?c5gyyd>eoiEg$#4d(;Wz8}lbVlyNtFj;V$%RcrcCW?TEh z;GT$#XxFXQW^Tw&K4eKpa#P*AcklVm(?pR;0e}YCrKY{ko4*3K>ZfTn>6<7Ip`-p3 zLI(y<$v_tfgfZOq$m!ezDpO7iQ~Kd9IecO~!0>8WW6>-%jaG1rzwyU5+$|Dr#wv$F zjEoUn7R3AYckh$p?K~Z-n`%BJiR5#TmX7A9i{LbSV^O;Ob`noXmRV%D-b?sCCoQFe75C_uEPYge*I z)r(<2OLEVP5jB2IDim1J5jlAWg)xLV=t{f1`YrFAD1mUf8~pF$^uSiwn2=3sq!p%f zlzt?QB3#}QF}5R0@~H-L1_Oi+-s+W9?3%RjyVSm+xrgm840RB7hlVh=_=g=$AAsnT zYVh30XQ@8RP$A5Q?k30gP?nh4!D%x?%b_en7L9(EKbHtmbksy!$F(Ha`Pl9EJf}Tu zanuTJ#n_kT7z|ksGR!{cC>&&yiROPJtli z{D)g7S`I}&Q4-8ycGj^tmmJ%@QoTnYOxb-DlMJ=H)+RF!Q6vYc`WHOp>0h0IOfTVb zOs&E(gbT&Wz9nB^t!D(pqs~_bA)xBQJ~(4_!4flFX)#mN;HM{d!rY)#C+2~9RBiOz zfwfih&%fM*9S4gM#c4mM#N&2QY|`8DSw0YG+TvlXdWZ#CS?|(}7|k*6cw;}```r9p zvc%rpV2P?isv6z$^76)ZrPrxIaLm1C^nV;TkJHUD!@%tqv$auIX)Czg8evgNH{X{h z!|soV#&A%UU}E69yFkp9T2Ka1L2I6crc;?O*yC=DYXb%bP&Z6nF7#goypjd7OV0hY z;uT-$$G>gdnLJGhPP&Y zxwWb4MRUt4boW2`NQSeELC|rkEv-iz+@Kv(+BmPrv0vIyVI^#=%cxw3MU8U{?ifeZ z%*<}xxl;#nv)FWlFwlNL`O!>b{1{|h?CR$^OK1^jBL=fX#+P2XM%ybyDwr z^hz@0u$8%%WZIwK@|pIr8$hA+f&+41zxK5RrHEsOB&m^@MJfH=RK04VT0!PON{@P} zIP7dJOc9@00y_IRc-~U4)_`C7EP@eC9?WIuqd0`=77^6mi*yeX?^j>UGisyFDpQ1U z26KW2v5(hhKHLClH`T+iRz@NQ1>1{p{^-8-mFLk$!98K(XqTP1=*>!_=^~T?HGvl< zMlFf~_Cd;JZ_BE{c{_8!ULlfHD_$msxT&NvQdc6FUXP1I>CXS~6*7qU0QRFkefqG( z?&i%WrhT~$Nv=Jadevc6A?<=%SEM7;(=uFk?Te(5AnE6419<__pfjj8fil}6=A;)C zLW7?R2RXLv;Ar)@47qGZsc$B8Kh=P#pX>aR=YRc8lE^p7^}*rns{*R;8z_$UnHj`F zQf$8p0dNQaioIq@nL5=t%DJXt9kxds=A{2TR%l(BHWfiSN(rJH8dla;O3ec6+W^{9 z0V}ZxMlXntjEtD>9j9<^PNfid{D9YP=JM{weAX^3;0h#(Wdoh!)xJdPX8*%e4pYCP zmJRs4>%Iqg0191~d3k${SNw&)VMn@bVU!8fk8{kbb@}&Y-`jgsc97>hQ>b}|)1I+t zrm#QEPc(Kh=r){@M^;E}#h9O>)g=gG^xUV5YhI92gJ#si8Go2~3bLchwk6@iZ9~0x0d!pTg2SF!Eu9eTt z#@=8a!2PcJ2#Sb$V24^0?orpm#Ap*A=I2{bjD99U&1@>6kA5@2GLW4><`IQ`R?|2BYov;^!Dlp3YZ$_mL! zRLAiz27#HcOO%VhEGQ@_g)az9O!lZ6F|=pkj~A3RDSu$31DW<)R+AvZ>O?~z#GraU zS25C0YL8v|)n|o62(b-OrZ70U^}WdP!1-XB^|$Q5Hu%jauuSpOq25>zc=xY|A))R4i;*E94(^im8#1L`e8wz=(*pQEC6Rm@;1p3ty5vapDB4 zULD~%xy8W}y|Hhv@54$?#pOFI=%uG ztJ=X5r}4UA1}GWs-A0TiIKBqfJK#L#d$L3UsxhD*fLYqC=Kb8YAxseiqQA75cT>uB z0|aROh?!X9rD6mJOlyLAO`fefayTu0IZh)eI9M&4p;OIBgu?xHQ44n4_Wdc6YwkZn zC7%7#&AuwY1P>l~QQ+}q9M@?^aQsMISwSe8d zAs`ACLcnPwyH7dg%qtoK@KH27nRQS;(vg4kKm5FVIW*+V7jYUkm!}$3i$e;J(*dTS z4`L`~`zV}WDV4kM1a4Jzbq^Pl>^pb$Cf{IpMvnwppjW(ghG6Q=t5=kup`p-`dH40D z+nFDzY#^dWBOS`iSezDvj5Yq`Y)!4%Msm+zya<9+NI6~01Ax_A(9?GrZdMZVSw5)_ zNHJHUZ=BhWD$^L+0x?g`1?lt`B9Y59cuX-oKuo>Z6#V!ZF=nkTnKn!`eiP{``B8 z+Gu}eH|dlFM4W22D=icfM9#Wehx*d*=6zfN&m4(i^?$uyYG#Jhnt9x~( z9c%}Pez)s75V^e4cCy5CR>w@xmc!^@x#oYEpHISCz5an(1q&qxHd_zG{=9Jk z0(WO9skY`@>Dovh;UMR=eFXslEt%@aJEy|^s!6O|? zIC8+JSs)YCDo3@X#Il#`1w*%?A|TJoND|O7{QLcAUd1J1JI(&PK zVJmi2XnX8z;0b$}u-(*lq*o)`@b&fdCz{BB^@ie1U;R#j)$eSrMGh7PkJmqU%&)ot z^X&N(OP#g0-Jr)fZ&;G0Nw$nIwQ>bDqebVD8Ip`GX6Zj4NfN5BsBze>uVdmO1k&) z&PuZgel3&%^h2)fL1&>Pao|Nh%OCD1%1wt#)oo7*sz)zEPwurFWDR7yP(vAq$Rj@N z*|9n#`T#h7gNN?agrQc_{aA5U)9!0F7N;g){y9GFQ^~~?y~)O?xDbaNW&Gmf80id*%7dUlsru7Z5%FC9%m)9bf?5u$s1KiBO3{Ft%rQ7!}2C!jOxvyv^fWIkRu zUOj+vLoKwJ8O-c>Grvqo0|?|>)$_AbWydVv(39;Mx5d5u9>v!+h$~q$%PAhkYxf35 zq*u?~j&mWqixY=>yWkerV?y@Inz<=q*I;^86KIkva8rrt_J+)-h2AQ!FX!^L%rUF8 z-e;?~wYbhycv%nhN@^b(&4a)q6eB* zX@m~8K5~4YyQj~QY+|AsDzTRP@B;^D_t23-ES>bxYeML?p+cV0VWFSZX~R@;sXh}X zv?XyXTC^o!V$iK6(=Fj+i?UxFG_!MZ>0*qt8vR=EM6*-))Cy1sJB@vFCPM~JV-XEx z)sya7ON(V?v17;G-rNiqbSYf5{7@*k2;36G<7KwOpscr!s9ITBF5Evy=e7??Q%KiqBNdhM6%_Fs3I=@iAf@mQRye1b@o_2sv`;tQo}m%%9)L75mo0Fv z@lZZlZ!P-~BxixACwRUZ=K;D6pwo*5-Eqb*+4xFe%hV-sc_$O?IBbrcNSR}5vQkie z>MNh5kQ1J|3C`kMu34XC%))v{ybMK`Rg;&0%lMg;H~l3kR?6E0EE1H)b?<|B-UeH} zy;S*fMRw!;peJ*hoKwJVrLf8D8+kF6?l(#an; zI9C(k7UBdOfo9n+BRuHi3hZ+jrz~T$gP7HOCuJw^KNmc{zbFhoG4Q-jh-fyZ7LY%r z?49*b-L37f52DkRncUl6wFajZBSEVoLEL+Vo4yg%2^ivKn$1${abK6CliT(NKJt5K zP`KBrFCIs0L_N~aF6D#p_lF}1(QBZTi_8xku( zwif72Y!DanmDiVuA!N)^^EiWIswqWlNpGA~gnnM)o+9NX0s^vcFp~eiuXzciagWwo zM5Wg5Q@;o1@k>MYRm201ar#KKoFYM28~MW#f08RF?O(B@nEfuYxzw-SrE`;UY|G5zu{xK)B1@5t=eEd4WSo_AjDwr{5HAdU!rByzt+9UU zGG?*n1BH;a8@Hr7CJo8o0W%yK)|2W)i;m!*Zs-2@BdbFWso|WR$!g2jzAgtE{-GnQt z_7~O$8_52cfP#m>uK3dKFlrQc7qd^icL6N#noTXB$HmaKI~}bEA@3Mas58?(ZURu3 zQ4ByqrfCl=(6e7c&*gDiIa9cw&qhADPjNUm^Y#e})?2e**n7FbBXo^JzGVF5fmv_% zYryUE{gq$BT}a-RlqTEzU*LHN$$&M0EGp{OfIHPfnSTtp;>7CItH5oac)NhxnWV^K z)RP%25c126;fI9>>~1a8{b}GPQLJM8gA0m{v|9y-rgdAKwuA9lEw0mxFS*=UV0|6stV9^9WvKafMJX_c|?;fG^Ik>fzs za5!2PS;LOnfctmx7CnS)GL`d}*=hL@CTKmYlc~{7830FS?Zek&t zggnB4qIE6};*vfyi{SX^$91}##u4tFv;x5U^Js}*g(*R}-y?b0s-FUDNRFM8 zmG!$7Z|9RLEViU=WyKVx_Bq{d`WBZ|bXfhOQK~f>zV{Pf5KU$;xu<*dc1m zIVlgUCykb*LIv6x)4_}SGFK<;dQ6uvJV>Ct?Mo=OaAv~MI+w2S?yHk45D@QBu!dPYc5yMToU3V37u zYP~4T_jEr(X$Zh;>6ipFsdlDbkzZB|l-oSDdeU^8w+zzDR1uhXV3+gLwCbrhz|}1Sg9AQa9CB`cXGuWDcsf?>Ax_anP7xT9d?|F zP|oqeiR25ToOE*#>Yx6&(7>K{usk|Eby6|m?xPOp2Ko@^euNVhZ*yecEl9jG2 zMgaWrEUw(kQ&DfAxatejWSfQ{%gQu?!z;A_hg}2vTJ>_CH@rU&ilF!pHush$sXTw$ zcu!-(u4UBP3FK$v<)Hr5Q{<2q7B*R`F&I|_Ip?iiw;c0<-agqhmqNFz{JI&?GQSYxiz-*>s1$z0< zmmwT;Sf}ZD%ViNgYl&3V|2M|`BNU|s!*EQ$0F%tuz`agztE}z?DVv%)m_)0bKx4zW z9;j2yQ+eVLy#ZzD`K#cBR;Pw(5Oso1A{MwqULn&qt{QB<3Q#Dfkyursw4VeCmsN)i8w4&?-bi2Qd0xEThWd)2N0+z0VcO0-ykxyIrRZb#8gzT zsOcld>LPN_otN{f3Ps_=lRv#IN1`PO6eCz6U<*tazBj)XFLAyUZVi=b*m8-K%REx~ zvn9ZIsa7xxAQBGoG6L-#?hhRn;B_S{v&r}9rUw4%@E}m|Gk47p;Rj*hCH@b4?-_{o z|Mrcagrq@ameHW>tYlLuTUps7o2YCuQwfz2vLz##*_+DDUXhH*CL??QkGJc0UH9+4 zuTJ0V;s3$)>^q9{{CwWyI9}@r*Y5_n=k9wTc&{&wrkZ)C!8Kietl?ih-XG}88=xB| zebKRVN?>b6%6In0L0pYnal0S*vZs0;ba$>si;P=j?)+M7I!sZ#K{Wn6h?+moZn+IR z7R*3}zeU0B2Nym_i)7)JuA;_Ud=e*1zXLXtKlbXd!TPT?FbW==CiQ36%)Re$3&1r_ z+#l{@?5qKFi#8)I86CGX&e13f-_8&_c=Iauy7n3VRyJPzaw*U@A?VJ4yclMn(A6B{ zR{Qb|H44ke-(RX(Ul)45EX&Y^t6H_-iXVYr%%*3Z{B5S%e{CigyGJNQrJ zge*?Q!AiC+u+ucaJg9V&4}*M|B#9h&U49&4PLPfv3vi5*8?AykLm&QZ=1pB1O4bOA zr=@{`DWX1Pe?Py_vm$&fo`hf7zIh+jyN0zv0|gh#R)4zUeU|t+OEmgwI@FguKs|^1opRtK2lvaPzc#GH)^Z) z2ggefkzf{+uF(H|B+`>h$4@Sx8|h?_CsvN#r|%!A2XFWm#uqPy?Z4g$)6fBo+6Q5PFqvJlOMZXu>CL>>Nyqj7PSR)ajrWWk@-(N9$n*Vy%UKi12k_JCMK6B(N z9O{zZ#%hQQIdbaqO|bWIATqM+_wqKV;sXKUt$WR^e2#-d0eXLS)ekCxps09At9X9@ zN8?0@clS^gj2)H*-gy)Pey$fo{gG|(g`fzot>B~XlZOx?1w42p=FRf0GU#u)Ex)n- zLpSCrPPHdlF>7T2HZ9L=tQ(-)bAY$VF5ig5;TgB_=PaGS@*(rnq5aqB+q0-b2>yi> zf?%3%qPyg}b{{@VK)-^OG1y|cfv%I5l6nR8Y(0Pz=a`t7^1=*+(N`8HAdD(=l_?+C zHi&1{k$>j6yq)Xu@|0jsL)Yzh>dGnFxo=Yg{(e)W9Pfd*6(wL>6BVCYPtKxpfiIco z(lsOv_)+dw5(Q{8)7|z2UmZnR$qFYRE1I+6v5#dl%^6}EY*qC&B}Q_XJj3tTQglF316m9bf3vpl6L;JTJ|O3+R!qiBq# z2<+B5+JeRxguN@oOf#@cQp6eXcnNMLE6^4JVCH)j@6Iueg&1J&bbv(QNCnPG8ZdDk zw-oa?&f-I{h~*pdRtkZohx(@ zRYh|e7UKn`Ep%-RZ$eR#gp~4~x}W((Nm_nJV6dVQUU4USN}Y2~2{|-2rbF*YN==P4 zZcN_|Ve#5bw4P?xF3eI%*Hi`1kZAW@?Q#NWL)$D}q2P*7O+-EGxtDe`J-6cH<4Z+z z7GvZict6(u?k~A*cQ_H^9O@h2xo>Z@$4MNpy_ZS(wW=8eQgy(pHqu&xRt6U6qO+qL zfh8;B3s#UJ3QHOSLZBaH2=LFHfbNz%U}kgB`b8{Yv&G-}%kH_aFW2UiNjE6$4Nol> zGA*m5sx{q}4H0s*-V}VRen=f7C-xJat24~oPg`xk$Vts>_6p{kE7`?fJi3A|5nocq zMjPLa>I}}Nr`gkX$r0(4;koq}@P5rNbR;*faBj`AA{aWS7F#OD5XKa5J3GGBo@rP`IOg zMYiqvR$AYIG>(^Eo$gjF-`TJ`5a5+x_*4#P7AL+|+9>|m=o(#N zf9pkF=;m_zgU3JKG(Tu`=-Oe$cV7%^U*<)uln@kLvg$R$y@61(^*c^5iPi`ZUM8Cv z>aF&B1R?_fjc(HHyaxcOqde_;YVcR3>4CBmkrh~t0qc+Qgp@SUwe5C2ghh46MU%zL zmqW$XHKYcK_t0D}F$q=nteb zhs&jarYxa885C?|Cz0h@wqrd-!2DOui0b7p09-!!k9XqW4T^k{TT8n0CZJ>>oTTzAQ+w^#8br8smFc7gt zhP2IA?7C~fv?Yx$`7n;o->=f*y|Mwk(vZ)bb77Y!LJk6SwHw6%M0<~E!X7}>-th^l zrVEEv%)2w9%a-B@ra@#dWY1E(9h@_ve)NRq@!m_F6n+dSjNrNb!Kf`_gxO#&4D zDvMrYEwN(K1H#UK-n{<U=0XA%%Sdu&SAusAgOJaIu7-CQ7*&6Q^AAdVx#r5%dZ@zn|X8uhH#DGoTB#W(X0hFPWS~B@f z0!hNa?Ae5uZdWhXe_oDE8jQf~JOeJI{U_N{;x!P&!j<~T{G$8;gJSt0YaU+U_k(Fc z=kO^~oZL;cNu>i_)NOK{_awuf*@Xf_QwKWh3AiD3dRvJ?&7Y5em5me*BldpouPn!p zA5T);Tu!utKYgoneR9XJnaa=lBkzOZYdI(0*~f51znA5Es(Gk~PCw>Vu}?Q7IP2?f z;KxxbgP|u?I$Fg^#~MwlH|dvcqUaDQ~yiq zxf8?`{jU{tc%xM^Q^bz&{z!wXK0dH7A=LOtCr`p?kpNcuj3EnmV{-zo8UK*qqfzqE zEFk#0AgrZhP;t6kX~2)5g?A#pI(&Y}>Mj7&&rk9zEcNG1%#GS>F*Zqc-1%)A(m7Wx zk{K9irIXb9wYlAx#9lzkv~957vL*_A#zm4a{QIPGx92co{inXsi*3z&M^h)&mWovH zT}qq-U?r)=tsT(*bHE&MFAy~8qT~0FP0~`YhJiLn@}krO0|RyWKij^npDDIwzjd@K zF<~MiLai-MRnyOQm$Sc4;p{{)rKYQ^Ykb2b@o60;CH`|uYTS#|cJ4HBveIpN*V9QU zDNX-Fgw)Ph*%EBG@W6(GbbiGVcY8yzMQ})CVj-LboJQfEJFLp4Z61djJc6UofD(-% z2d{LSR(f-3oD|@(@uH>Cv~`ri{%Ggx6Q<&Ay{|YK$DjOo<2tlb zDa?sGdE)1|v>=Xvi0t9ZyAPy8-AS0B(A3Zs-xa%%zRUS~ls@Jlm^0$Ye$GirtOi~@ zO^;iAB!U7rU(|D=S;Y%5Jw|7-evhw2o|7`h<$jBfqk>E+9_Lguv<39SxV9t2-?3I{ zc(mk}n+6HjD;~5r|4hSYAqCta7_;>9t@~L34$PS|!LRO=E_l>#6`VG>I6c=Iz3Hfd zbFpblW^(S^?~~3c6ESz|Lx6GWy;{oNU`F2;dmJtSsUrF0u=n zPQS_LyT1bJd`nf$JPYHM9qt!)TY;ptAXDsdY{jFo6k;mj@xT3u8?yp>F@NibWJ?7L zCS_2Z7y~=RVe#UriHT9SAJg=?s6p)vX=h$uQY&O0(&$nt6FQ*3QN0E{O+C27Pjhi` zA-Q8teYCr`d8I7>iCjMyjlqj*dQf>3?^Ahb++ffvf1N2C62Wj=2{+|XF{6ESqu@7m zL8}!Lzfle4)U!VL+PZw*q$7RPQ=_`*NPNh4V(o0ZOmBb~$B8V$Ew$+sC0(ECuf2L?9bXbY#*ub3hp`ooFgJ{_|8a zbGSJ(q}_;(KpTvNX1>>b`6fowCy3AJTvf^W!Yh-y8Y5``;RrQ%72SWmDp#T^dI7dD z{(BCz6)wYtsp&8LWy7SZ=r^|oWL2VzgSG=iaQ{vN86}8c$Is)oQTXJL4%D`;T_pWB z?G8{63=v1(rPc|o?q=iaN8p5>y>V=D^C4^{xLe0Q+KJsWy&hf>ySy~nk)h?pp*kt} z2LGhQam2ZhcaJL9SGI*rK=F&9uuVe z*+#C{A@H#52e$Rr+xb?(9fq~%gjRFN2IP5@er6p<9qzOF#j3+ZI;WZAJUH8qdPFCJ zj5FTDJPCCg;p}mnR1->^YHc2MZ@Pcc&#BisA^GJNbDWs}l4(%RL_X=(nQrDQBI@{Q zzi+R(9Q~&F%=z?Biq_rTS*SO^XGnz|*UmQbHLGkO@P#s4y!l&d~$35e#JLd1q=`LqxRs?}(m$EmL_wzt}F8;{g>tee9iX z>fzJHTU|&_%6dVKQ-FDb%Z$XF-DKibk`B#4(ueWSL__}G=8EOAaI7n1vf_Q zVQ!O0@YOlnF<%n&2i#i{!~#mU`WToL*wEIb^EC+PG`9OpuzPb+++ zN2Fa~twVMK`j^1*af4nOxj&67f==^s^=JNcKXqi3d7U+mT91+#m1+SF&}vE_ZOSri zYinD{TjG4!EMOF3nnmxq4@fA%5huaHiX(YIKLXzbKd%(OGOczU3b|_fO^?atbYt@L zht3aoAtPeOy6Mj0Jkl)88v@f}huk|E>E4!+FhvI#X}Ah>oS1^~RXE zm4SG4So@n|%QghIMP21d1NFJ{re^VRZ5dpv;%o0zqrx@RV#NlcEgdW_wuNWEQH}XT z&!J7`6-T*two4#h|1ED>w&Zn!?O$9H{{^Z2t5Dyo;IF8P|2hG#)=L2e4}M)osLS^I z7@UhbxNVwDTi}Tb1aWZt^K_*|zcQ70w?Dfq&6l6qPBXSVXkGzXg)!U^%^bJg09D~O ztT|GGVw@dBR9JAda36r=A=`D`=Eq0)j@$;H@k*k#8UF<6t%+_-yi)=Jt@&VinYvdj zg6sec2^nhGnb76lkW9f5(#cj_?f#%$WE0fU=^3Tkp3VDi16(4tGt+NZ-GRYQYgl)6 zOk$hTeEaQ>kINN?m`~%^1YlcncxiLJkkrGXbF6d0r4@4Et!U{ESAp4rrw~lMy7qgy9F2@l;GOzeo!_C+>Yt_ z4U`Ao9xr*g`^`-=JhlUky%>e#_9BfS07z zY!Y==0VBu+!aDM5gs-uc&A3%srs#_)iw_B+!31BT8mKn}YXvd^#bcrg{oeA*4A<`e zfhPX#UrFxZ7C(;5>rTjf6lJ;+wp<5ygsC(Ln?61Q!mBrV!hYGErVxfpyy~CRXzRHG zQb=RyBmp?Ov*@xu8YVvo*LLGPpkMWH&i{6{kjgzkfq}fTvZI#)v#ejQ{iA2j4o3n- zKV>6JalsogzeRu5FwIh@h4_{q>c@3y??anX9BQ%T98~KCwL@WZq=0QwEuo`w*vl;l zq4?bl7wKSzR!E~e!E|})S!}!{1IY*6VqC3XHM8U~rTE25>mr*dbiXa}crZ&ulQB}i zgdgG75xo_ry`+Y^kz$!X7^{FZHcZoUJS`CVVs`!qlC_;1x?G_6$00qfm+}>+2MAS` z;($<9)W=DJ#&DF-1FE$p?dI{KjcP4W4BKi~n1QkE{CG>OBJa-N?oL5tjL45A(}qy_ zQvCR|V3Vc@#<{Z980X1|sOI0i68e(|#RsyRw-%19 z$nS$dk|aTW7O!e&TnHLEP5JJJkXWQct>Y3?LmI?_+I$41Grb=@mHZLsULtgdNM{-x zM@CLQN)x{~mLgX`x5UA!`{PXuJR{;Pz{I~lg~RK4u~Y;HG5kI)(S&x;;GVj&LCBhiX9go)>bMx!6xb2s}iz z?2QXTGJrDNQ~uV=n!oVZsoA@JAD@jEs2?jvRee zwzDxkmN-5oYlyRlHVHMcb+B>iG$|A)wwF~v%$#5wCJLtSmVuiwDaBJw!1+%CoEDNx zCmBSliOx*&rLm5s1+J~TW# zepb#eJ}>v3`=_Y}EeJaEEnOS-58}V4d(TH+&qU0-#rOWd5fcA?dhQ*eN{n${>~4v< z|0!c>A3Dv$!63_^mZ{r+%t;6ALavDui7M-GgfWpAzv!Mw6!W*H;fF3WMX#cK`&9Gh zq=_9QOmca>A8DtE9$)`3P5CE5I(P(V@er^IFwrEfUI5$?h9fE49Vkfvy={H?_K|OH z`&%})*hyX;xu|n58;aePvAb}i2Ll(9ysmLxG3IWvqD9SXTd^!&J>GCKzPe>lX~}47 zM@_^SMjQ%-DKTyGzm^g=9TOAG^jR4t!#GM2TLatB#9L5aRTzM9=?7^5S!mlf@i7-Y zx0Jgs;yoR|uL8#lb2x!b=9iE36t^o#p?Gi8@t8dCgG||Y4e1nitpaPl%{v;PKWebO z%0K-^Y=EydakcSzTx7FePMO^y>2dC$NwMp|fFrapOI43GMFSok85 zI>8&!5zvDeL7W+KVlFm<rFkTR?q}EE+j?hXK1Kr3U^=(T?-t_Wi$? z;s3f;Jc*FMWQ=@@23%N2mT$;)E#k8skTXq@$N(UvXL>*nB-?;7QEwVJ%=Ve%8V@`| zP*V_CDa5~qvReiu@2jC?Ac;ZlLS|?9J4PS@;WlqS$9yaP+6PhmezRE!yWwkKfS)U% zI?he?-}c$RBaRHb_z1dE)5*^QwFlEHJ}`*Dn@&N#_*5rAK>%h(9pWAC&jTL2Zp<*a zt^1?hJ}O_psvgIRi?3o9N8m+OJg4FNIt#U9?nJu}&=5hm_(W3hd7HR>0Ea24{Ri?{ zypowj1?S>;Rc;2(N5DaPfaJgXmSL=8Q5Xz1t$Wmivq8rvhZLg-tM{Z z(*QiZ@Ei;k2R0)RjD1E^-q){PtE&q<%Lw#axQm|kzP?ok5@B5T!=968!Q!P!M9C4w ztQg69n*)Li{Olk!R(ZBo{@>U3f4^%UDOcW;WtjR!!K`#`tT~!hEB~(a(h|JBFz}JS z*#QBrl%1VjJeRb+S?FyjRhy>4n-T)Oo~(+Bic<3Ze{jDZ!P}z2C|7`s=uKGiXyJoj zGngKTem=cC)wOu)%Klz)jF-ycHTy$VDlJd<5$x=7`BKTAgFKb4O>+>u^%_w8$SeC8!mYJb@WA*Yw;!OCN&E;9g;7Kn7q;9EnK;cLi*+}t_=#ngFJHm$)w824lupthlH z`sV(xoCH~Yp>iQ5BRjh>%SO}lM-^mQVb^?+_4p_a?<+`VWD>l75F0JAhrNqC9RlZI zC>Sh}1|7KTI!4_jIqLdop@P(w@9#lG0hHwKD=@mqfXw*!RRAXg;rjdv%=d56e-P^o zod-jI4+tNP0Nc!W0@iO7pf$>-`C^dM&bYV80OB_~nH(F6f0*LU^Grk^Uqya?m^hgG z;~z9KpB2$-h64C|FHNUrnM+&}aP3PJmK>*#FZjV?;c^Z3X91cdLO?0iCJh zsA1@Sep>=lU7Ml&4}vOK<#;FnsX82gWqmw0S=J3U&K+YCMs zzMQunwBfm3C`bxC$vHn#Qb0)Z`RFRJAd%N-OiQ8^8bKV?irf9Uu&_)YMSGILHd*mC4-&tJ*85 zf_iGB$oqNNz@ba$vGt$`Bu68Fi+m)WUM4s~rC;z5pA?FUm*OJXvt4|)V{+!LabBlH zU3+4l66B`EsX^tIC>PGH1T!uL9t=M=do=?$3fQ6b0}5D*Y{x^Hm6N2v^Ow#sfW6>F zEB`af)Vi3r)1bF=26P2MfbQ0T8YcAFF1*>0~JGa?6OB;eUhvq8AUUXc_W%nkco;mGQ3D0)0c=} zLGBb^e}CJQiELlUE(-;CLYUBdJqmDS_e|MK}0LrG_cBwb&K zPtz);Y^L)qfsiyM(9% zMH}080cgT{K;jeOezV!xaz>S`hxOEWDuETIv8APHbpNHXL(&7s((j8Fofjs*=Op(x ze+Y>?;<*jjOnhXQMkZtd>^Y))gDI9E1@tuVtAX#|Ur6=;hl5CkC`#0=SM(ita2-S- zl8gm&3FL6do_@&HIUY)sAIFf?9`Vu%C_@5VKa)QSGE*MLT;0E%U@-TBgZiO3NNZ+X z)H3iSF?S&%19ABEpyUWU(*VigsvwABgW&c^4IdkDexu+gS0!M2bcvZDj9?Nz-GAW5KjWYg+V}XXaMjA@RJtGG6Uoy5U?{naQMYe za04&S^=&G>n;FuPko~jtgG$O7V6T0f+g|qcs0+9(y%!k!Z^;1>E^3=+A!`J@)>6fe zbJpcvN`3_p{?vU_hvinqaiBjhY`nO(&#CX#c^F0h+2vhgr*RT+oQ=wbrUcT811Ijl zozGLx3@6K_T3-DQh$9gIyfgkIpGD^@5c4#kOey+~xK_)|GA5|>_}xJL#_bga5AKIy zR0xneH_>@tBI|wCi^x+T8X+^qdBK7jIE)X2JXB&}qzv`t`*F^WD?-!7f-wZtd1eeV zkTvrfXc!@%I>@dX0ktO=4#zZ_Ov8#4Xl57Nw2B}L;9BrA7M1s}CbC4o0Nz5jEH+ZKkVTQj?o8BuN`c?5VQmJI}=Na#7A(iJ4Z`{E1A^z;Us@2ge$NE zIkxO_EP9{5?UjueKuI#NLdys$mJpa9f=rR1ggV2jmWD)h&LACCZ@N?rs*dWS?$;}p zH`9Oz(hLC_NhW|!>plfY?`@s!X@NeEV`wIW@90$yzgcq zv^EgjPstaazQvJf`5C4w>-{XEKT&coAs|QoH4h31h!8OeLGmr50CR4o!VrKvc#X56 z;~P6!b)}nEe47cKOAs4YZ#t3>PuwIn3l{ib0=ee+bP;55hj)$arwnetpl@tzp8XY6 zL5gj1=XGfAK^89p$!xn&6Cy4L7budvnS_@^)7JpixuiE|b6m(?Fu-VpJgB&gfXi`+ z<;YnBX5@!e`otes8pYPco!p5(?hyhMY@?|E_Y%=zx%um4gg^mct7Hl7o-wLWqH)fZ zF23wqTF;<1)a+Cz^We>4gFW#Fz?&cAW~6~30+i(B#<(W%GP%OdrV<)1f^H<;q?CZA z1i)f=1bHppUtuV<@mPv&o~YvQZuaC{MLp*5%owugqb(QpnOdv?iiw|=(oBSRI2KuB zaeg&p`sEa5NJebREd^E!;X74?WDlAFwzUcSd3m2IXN{z_($`vO3D!>$?zxRkSN07d z7k0Q!F3r)rb^#lYoX&1*;_+pBINOEl5+8ayQ|Oo|W+r;^u@|3;BN&7oV%iI1jt74` zJ<1?ig^9;M@E$jQVpzM%^9^#nj!FTH(-IRn8=Ls43^tdVuFE82M`pvZsR8ri!SKxZ z?{F3+n1oIa+|jgSmeG+=Y=DizjD$idmilt5l9o_`m`~FezfRzAd%xH=dVCNf(fZ&`nkV?11l#Y#ig6`!zb^9KgK6& zYzMwk<{jJ*4Dm9;w+rCtH*E-~HHuE4);NCL<0%#IX8l*~>>Bkss3P({#LH5`}G_WS4&cY*PXz)Fzeb%^#4S^|>S{Ka>Q zrn7y!83bU*S}S7-@_A1-Po{bCna(#jb8WlwUNA4RaOejW)6^-__GPgOL7sIkdmiD` z37<@QGfn@D$Q1f=RadtI8nyc&sa^^Adl>_Y!Z$1~UJv8U^|O93(s(9&V_5`(!mOv- zpG6!W(^zuc0(dlUHM5=Vd_T)1*lS7qm8AU%`lA{!F1lxLS~cmX0@~a1BQ=D_%-Hze zx}b!&FUSLzU7VoHGcms+^8TUQ@?=Nt^yO6L+6fz9HR;)bZ-qCVh@~(ypsVML9MVsl>NP z>G|5*Uai@{GFd6f!xn*JU95#yB=}+E8h~lI1rpHF>8@(-8C8AW6{P_mw-W9=Q(kW$*8W58i!>zGM~g zoK8frsTCkbis=I#a@^zY`sU>q&SS;EK;M`|R0`7q?)(l+{D>zX-}VPsY#FfNBE}_& zR?g#J*$7cVJj${Mw1KaYwJFW?ZeH`5ZA)4p_jB|joa^72OIJ>Q=@yJMKgb*DV8MdA z3gmW9hRry50t;okR}a?-;v(1*x{V&P7=X^kI=(y3KY?aq4OD;?EOKc2!-OO|G? z;wl3r4c+pLx7+^_nGi%nvr>yoftl_|4j8Fh7{tqks3>S!{n<9ozTo1MPSYiz=c0Pq5?@@Gu869p`sh4;VVJ z@&=-)dX_uYqHwZ61llm$UlLqa1!5IM>($r+3RM4+I0OgsM60W+>cAXlb>oG*#kt!D z)S)t<=r3*3YSY&M&J1xrCV&8+x|sVQb4@8u9H;Glf+_`6l~5Lj z0;wD>-w%+Y(cEUDVXodPF6^R6dm`+zYe@VZ!3Dv7)&X5X)uZkLt99)RicnfLDVsfg zW8!CCqkb_29%BOlRHX9K!5&E5{&2DeriV4n7>Ana6A3|gWd!3O9oR-tNDj^cg1X%t zNBa8?M;7@rUsC|RC0X%KO&SX|e3fSf4-oxSQizUt&S^bT|NE|1Fb=%5<{>Da6ic49 zJU{ttnBd*o)ceA#95r&o3Ll|P;!AggT>qz_mQL&*jal!OE$qD?8VQ5qoIqd=WgppR zC`hRx173N_>`uZW8YHqeDILy0f9vEkSf-yaw65TSzvu27k%tc4uhv$#xSlC%(Qq3{ z75;p5_%m~IR8Zf`H=~l>aEjy=|o{X7s((wMS$Ji zmrG*1i(&>pE4_w-WH}Bkblu{bP#tW0weYOIP&rK`@lrodsQ1gZA5=px<`(@9XEkt? zJcbiz+?vw$H&B0Tji#bx!_?sqZm4S;eagf=Ly2wv!3h3dfInOSc2QjcqH_1Pgci5w zV01-p2k7vhvw}y5zsx$d@tDvI$*@p98pKFqNVU#9yCDd}ID$4|8jb2sI^WmY@v0^X z517rSI^}G z*H@$JLI{?j8YE^+UTD&Le@Cu)hu?=E03MMH-cS{xwQb}&?d~VPpi<~mrKd9q10PeS z)%?Ri7qM}=Jirpp!uk})g?>FWSL6$Y*4^WGNYdWgL0T-IO5(Rot~CnzVU8>+Mc*4- zHMixeO>4>4pzcz(ESY)&U0tT`>-Fg%B2U%s=Qtl>HXg^iO2v2twvWhgQwA4OjSiRm z&euvMX0|@%H@9k`IJag&Saj~C_+7JCv)D(L6Ovp|6UU^fGcRM?Ah6`PPu~o>!JrV% z%_s~Cqag^qJd9?&z==W#b>>kt=pBu)hlsNL7o22+x&+dmfB&q#*Y)9B4gslPp6^rZ z`06(dU26Ik`%+Ua&Inf@Vz&-a>N*(|0>+41%#%^}!nHPN3hirUg+ox2J9sav{%ttX z-f!y-sD#mMC*A)9$MXJo49p@8%%ADKyFK~d7{7OmunH=|CorRv7u1#Kv|$vJB!0!2 zdLWufIbf3*o$mqT#6Y%;%Xx7upPifW=B-=DBQ}GHUrh^D$d%jS_1Ya&S5;NLywYR% z1$J%(8>!<GxF5D*d&e6!&NIL)3u@KngFEEssw0tpoZiKR7 zpLd+g<8|q&-AKYQ^V8YM6L8O2cM(gr3kvV$(h9pDtU=Y(YqJz3D^zTDh_kOsU z9VPW|ka$7h2lK(d3r||s`e!@dY=1?9*lXze74_q#@ zmSM^t)b~Exo&?+5B_iw}KYgzpl_(>}wqFTOst z5Pb(QcDcZpkk&lqp7r-oEn|nnN{dZ^*rHxf%*N4y$Z5 z5*MaKvi%R9o_`gDjuN9wEuPl#9_%Z31Hjr}v}by;3~Eye>On*2$sbIRRTBcxJ~1fa_pS^}eQMz1?L!0~lYGYHP+dZx1#mAJ7?d z%E4Ps>4x;+`{2lW3LSO#_vWL2pBq+}R5|m}Fk9gMJ0|fsdAwdWj05Dg`v6(~1f~u0 zZlQr%&u&oHta|AHNjL;~rH5c|5XWwGCVvLT`_!er`(MDl_<0`WA&=L}Gtodbm|Izf zxtpH%ZEpy&goX|!R8m%ZLa%yKY2Fe@i&cHqM- zPQ$+O9~|LYdGEr4WttV*0!E5)O#?WgA}nW_upMNCBz2Eqjt7NWJ=LMTDOi8qmwxbB z&XVlT`-nXTli~@H*0D4+vTkffLGSP1p!G8@F>Tvt;tNBw6;8^Bh?`ppU%wZR1ZB|5 zZ>{)qPbIN0{#V90HyBKum!C)Pvk;z&f5)!`rbOK^h-klqDFFS>mzRg1%|qPFrO2gt zN&nAi{3V)iW-_gOL(d<&%r788r9_n@pW4IWBrjpA|I7mLJ^%Te*1di< zxAX;lJM0gbp*0EDiUQ-GgND}y>1g@TxtKi@Dx`^onF-YBRGeRL@eoA}vDsbQ<$d15 z*PC(j2lNG6P~~_R?)M>}2V2Pb)|TO_I}Br`ZTAX1Xg6rmGI9U<9`_!V_$CSZX0J|P zC3+07)QfX*%N?;hxpAfW{&1~cLK6UmH)nG23m(Eg0~KKLBW`Rp8+wbsW#gJ)&d~Pn zd%^u$DC4FYIw3j*sL8TmYd?j}lQ>$`2q1fGAI#tS^%>eg+v@u6vXx#|#J=o}z?LXF z{P!J+E_kmM#`)L@W!9FGt3JnuV5B0^!`UL!y|H1#$Fj7Gfgt{T>Bn8N;kIWB+qzuU z*>K_sRxzl%Pk}C=aM+JE$kBH+tie?IC)6qpW0hJ4og!L8)N~SYQFvRz%*kE6VKqcR z2$X!3?}Q6q+@G*&RkSwX8vJG@5;V_?`N$V4Lf#8mp!QE_hQ9h_TI$~8?PPB?5vsJg zN+V!nv^Ss&>b`JRnR1A|Eir)OhZvoFXn0|8v|g9D9@s>W0wqGW#q!El0_*xW8^pEB zKfo4S2H92y&;u9%WGwwi2TbY)HbOgk$ABVqcken^0sQFGwT!GU$Ar>qMDV{;pvi#w z=UX4W)9ZcYHXR`@K(ywp^fJ{#A@5Mt*#t^9^_%xo?+uo}1v7G7BZPL~XiVu(Dp=T8u$o`5VPEG5_I~gslZoQ_r z{`L0fanaSH`<=ILJORTv!)@*PWi;y@4Utp7SGG*_Rp+&j_>R!fSHtTj5}q)zbj`V4 z?N7{)UR+3{`LIy8w#|uA^OnTH!A)nlxoxVrsVYxkk%qpQ5<(5ZGMZe)l5QDF%yoxgfF@y43JpPm}`LYeU^th?Fk((wy0SO?WN zd(#Zjy`ZA1+C^jouFl86cnDWzl49R67_TW@#^X>*q_2(fS5TMFKjC}~SYes|b2O9j z^4J?gqOU;4ojfMG_LC$)?zyus_M^*FTp_OMfIe`e_p;wjqw$;8a`oL#3R`Hp8y!r2 zCz;!ZN$)NiCG0YVK#qy6MqIol&0n-?Mal4^oYC% zK%6fyN`el+efr_UGKPN3JMh!b>n*a@((3Qd?QV;eeEM^11pBG;TODn0t~lhB{T3Hc zON%Oh*Oi>u^aq@?7=NKp8L*RELA?HG^8tU){LmNgetmp44lIYVcnQg~O!W;j7k+c+b@9lM4xgj{HAG1-k0jWkbcjN6=p6TM&max6!P60C^h0& zmw1cyz=)f{FJEuE_jIpQLQSnP-yk$?p1su|#nf?p_ zG`+&*p~rmJdqKwg1>Hz+&=l3UO=3Mn7+SyK0y9RqaF$O#PPro@-8mI9!PA`2JBt9) z#Zx5;sl1W5S!DP50RwowcMqz|O3>T2O0ak3xno};c-#!=9OTwFe1>HM#_1eFN4)_G zJdJYvZXDJMopTmwzxEy%D8z*uBluX;;HohI#BRt^Z&MRZx%N_ohU%GH)v8+cZAvS+ zfDK@Tn*rcEhjohWnBSK`FvmO$_e?NG)g!b7#x6QqV*m+B1>iq#d6}tY(X<5plMld# z=K}g%7n`nv19HlB3j$ZE!QmIv7eIhE*@s!wHz6Wx9_aeyS%V}|6$$6Wd&FyUMdQ|n))l#whKP%_MTU?n=!C=_8AdNK_8;vp##Pk64_yu7! zomOx`9|n>|CNolqEI`B+W)ikrgRmid%oy7;@EXQ*U;9G;6rKZXiZHB)c`)lg1~B(A z6q4!J526XEpFz}j$^k~mjtoEmPf+)XGPrE`rfT%tGToO`(V*BR-6K`nFrRS+P(-{V! z(33_w4ccb0-K{Mab5X%sPqia-2&Nf=$@GMPklVS4ncQKp0vU$6#u23&PXoDi%N3}F zAA=n2oaH6Rx6XCvux*z<|AePSn5QapG)8lTunL4d0yFUWCdXkgu|+GFJDDqtJ()az z-Q8cF>rePDbAGf7mbK(5{(El=l2)a#uaRrH!bJEXR=stZY-Q%z=YBX5y^r)h&H_C1 z8GigpKuM`MJ5G4tDbi5^%alB>7#E3(NZI;W3S`3xnCi=pdC< zJHMtIy+8AL)`2Z3i&WkzWH1n5Y;9j|-&Ro>vF6*b>|>VLF*U_A zoXEA=4}}jWACoB%?g3;ltBwjyxxHUaemj8&ioNlmpDM!5^H6NdJG`h~eO&$SxCh)o zjLWtmze+xRNG_aaCa|`dgw0O8pekErH{d+}1hQROinL*Z?PusY>F$WXXIF~4Eu3*N z$thc>Z&Fpu@5a@!WN2m{160aZaiPBU-BH_5*=sB*PQ|)gG>%r1bJ|m&2U4PF_;p^l zJ9Jnf`*Qi(Y$nx(#zX|#d(}9dFyyFyM`jq8t*!j~%BY~$BDU{{=STW%Zo+Tw``6zq zU+3| z+|_P;l-7Fy# z)DwG-I@}u#VDwnM?b@&0oS2VMU;ahBxM?mU`qBzUJ*VfGYZq*Z8}BSm-+kS75?7yE z=kgb9?J7Gnn(yjT(&9nf1!T^2+DUJ%b>{Mp)O=nz`7XL|tu3huTW}Q^vdtcYYK`T| zi2x>a@&~^~;mbp#VR{@?Y8)3R+(aP($Ok??caa#F&gVf9LEhPfA2-~bx#s?nPkuq? z?$Jb8H=^lEKR!O6Yr4!}0AWJpj}Xywmx6&EDlUL-Q28T9Alu{Q5YC^<6bDIks1+jcZ?a<3s+= zNwSd->24Wpbd+y)n={qv4%}I$oes175upd7l-;EVaMtTbc)&QicF%T7C~Ol-Ggje^ zc~r|BJau2Zw>HAa?~tb|v06D8WG(;H1zE9T)Uss__5XI3Apq((UIgp=t~b%N$Tm^7aXbW{pU3qyCUmHp#}3PPzNvs4rfnl#R~E)X=k`W+aC_GA z3tY#xVOBeQ1nE!{M?ODGcNWbfKfkVxx$AxO5+QK#vR{W4iqZ8hy#sd<5kT1r1NbXd zWQymEx!zoYYyMnC(k>K?*S&Q{x_yOaLCD*!AIYXhf@v83K+S`tQoO+9Rk5?2X#fc@ zlWFeT(|SHo)ru6oEC~hnB5OoM)BMS$W)(4l2cYL3Lh;VmLw;GaNJ0-4MHUP! z7;&(r^`H`JAhsTj(b=|ik9IP;bL3HoA+YYXs(k0MJcL)Ly2bVQ)x*pLcx)M(`L23YdyTwQ zO(UKxsMxl05i0O!XTR|7IYynDJt+kXNaV#vwDt3#7$K^TxtE;I4Z@V#=ng-o)^DTX@28=+Fwbea zBpVC0+ritfWUTX(eh2!@&jyQ z_TZCp-C2xP^cz0@7`5R+nR9Dj`ew7$A3%L{n(nxUoFL4ow{X+bbe-*nQPxXT^{|SC zX7YbnSKEcx3*zdBYhUulcBFl$givH%xCFIfLbUd?TKchMKR>_wU}S#W_wqd)=ZsM< zQz`4Jy=?RWBr~l;dT&Fi68us&W;6SK{uw+F%9_fFp%@p^uk_UJN*G)#r78p-TTN3} zWDx?_pQ1BNfWrnN&`-M$a}=vz0>NW~3rMyPKYH1nrFN4p2gjR21t;%WSj@9JJ!L{pDm+$LCV@CX$+B7O~M$liVr5P-_7ok{8s0XWE z&}vBaFUAIeL!KdvQxPYzZLwQkXi0<6j+s_x`5ke`$=zmGP;?Qao?P>k_tl~}1puKV zj^<0;Wc+t-@l%;i>4eJ5U~X)CPdsh`7o-Yo#cg}^;K@M`=wDRd&2rVAA{$QbIT1kP z3X+BCv>hnHJDoEb??pNEngeer4c2E}uC8`LE^!BFju>NmF1aHp`H^EnV#sc^uTaA< z$49u|p8{UaRORU5w$TxLj0;2Qt^b=bPX;sNBerEQ0TPw6q*xF*e%da)j7mfVH5c)) zP8VlJVMQZe8GBLn_{H8)hU@+_%veynyPTcc&4QDKiov$~ndKW6&*=9BGYs)c52SHv zPLX0;LU@F)OT|-XFs#n_2hL$B0`b1e{2wY%p!)x!0wrZN4($|y zO^Xb!(GBL&kE1m`y?Ru8MwD#iO63zj<){WT+L$X=FBHBQ624smdg z!pivH5A(N&{J&2I+cEpUPv(E0%)cLm|9j40JEs4C_GIFeLkYC}40S*xg1E*vP`tlS zOFIIkt-HEuw7sU?d#E81a&7=RK!h|KynAs;CZ<20xJ%$~N95kU!f+NBS6^QHJy?=1 z?18?k6}xWgAcK+M1oGJaGBn^{fLb8>+P!oKayy0q>^%&@I7;gIcQqs&!E|&Oh^QNv z|J0ogj#K~Cq{#tU+z}T+=&d>1soM!~(ESGXf=+j3W~`qe?lCG|vZ%N)_6Yn1t4epT z3%JPk)LVPaASyR4qiasoV1A(DX$Qmwo!{?*t9}R?C)p{~^%wwrl92<{d!KC!qfp&c z`fEoXf=a{z{{4m(r+{BdQzWy{^Ql8Yi2bNLq)J{yN^Ho@WXAp)yrAAhU^1igF1Aw<%5raUADcJVESrWv4ki8a*z=au zw6wJ2TYjTl7s&&}PVrf&l*$cs*+bKOem!aH7bpQ_B7Z~UJp_0k`g7#N3c1kIf%h;L|$Pz@!c4$tSfKM!Qj<*tD9{C3tNYuUnF zb41BifpfNJARAsalo0@!t@Rz-d!LN2ns*#@Q|&mp@Z zHQ#03c%cZ;*{f;nqX%n-k6pyx2LbXlfTK2I;SNa}NGoOj8$pXr-66jAuxW$)+$R&m zKEtZky4|JOV7Lj0_IadUL)~B=ott{u5H#2w@!A*;_BBT)v5S{r8sJoc4t5P6H#-AH zZRAfG_yx?H)|%4_&!)c3g_b*66=nv;jGYXg1oo`#2K(4*9ym)Loh7sU(NBKB2;2lR z5mtpmPth^5C)|XDdAnzZkTLK6tfcNzzg(dVM;+B{7_qsQ*;SW&mn0I$&?d=C+SH1-dd zCXSb{md!l@eyg)E8QW=1!f^8goMn<3jR1XiL;#+z7LmyrtalG0*q2pxem&$dX(U9- zvECLy!D;V!rn-D3_LK*Xx=K-J5}5B}NJNV0lFOjRzLKm^?NNfbCIc9Y{XTX2(Vs1G zgPr3q@k6v9(;WdK{k3Q8`4_`V0$*@p_3>S&uh?Z}VX=74in=k`!z#PsCjx9*c5_r@ z#q@UZzDoQA_39xwrqzOz{w!yDh|GW#N8Uoq`o)J`VF+={Oqmx_DoVF=ZZ=*2Gg@B- zQB^uIWW20!g5?H@FC2fNMsqVol!?IGCgsg9V#9uLl#fI_hE_jmeXJSjoulm9r6=jx z41Z2*4Q{BV2LPO2O|)>icCBFnQNt#YR-s$mQ~?j8E*f0So8m~Fp4Tl%o2YEx5 zCYaBM?v5bILv{q}7=YhddYLM& z1S#$ngnp~sTP&Pr1*GQzB{+nB1^ps1+~yx}dI!TGqG1HX6PWEhF4+}OD+`v~N$LJ# zhYa$rTj94^IA%1*2$OEa9mjQ#;l-8hk~@q*-4Xh(Q|HNc4^qw*)GcMx2zxaWjqHpU z_2CUKEWpNZ(bxE{;7vsRKZif`b=USTtFj~^<)>Mz4}&D1}9~0oR9cnyp9}N%XAO) z&x4yfkq;EDkz>;hBi;TIAI78{`ySlb6MX#auJk(ag|Az?p+Ds9X`zEIzs87e%W520a2wNGBs3O$7p5Q*LU{2rwfBp6m3Rfq_rcUhVUEHyk^ zZh5wx*b+3p%lR+oCch?9rD+D+!I;6PpmIslfkZ8q?}sP~<|rd$U>nR=>~2;VfY7V{ z|6}jH!?|wz_+g1s%E)X{8pys1*-A)-vR8;?r?Rqr>x%jc4KsVo-Xg1zva?5svRC%z zd7tEpx?I2Kc>a2h-*NwQAKflK6NLts(&16W^0cM#zm zb29<^aN=E5S}>_E687i>qPI(`MPD0dzVogTw@Q+H5183UXj3+K^Q_b z=i;Z;C4VP&>}RLGeGt@0S@;FCZz;`X$We}SmUI3J3bxi%IPA?RQqZS z-H1IXz9XE!yESr=K_2X)Si+SrkIYfTOS{x3j_hD31cU@aOa0P@!zXsc=1s+$)LSR| zykGpbQ+r8GIL0R1uSl3_<`TT~Qfn5ac>TxK0l!ldz1| z(yiC-1U+oplgyvM)X^f8N)mv;oT#tH1OkN&2|!$e@Z7pcA4^I?vE0O<*19zwAz1)! zKoI4B?~)QI!yDGQYA-{OPqB(W_$*4b%f6mOX6wi@#qluYSpn4Fvtcys<#h$ryiFmC z14@G%FuV7OKF=i9hkW}r^N3lu1g}frxVICXJu(C{3^yRrxdA|r8$bwkl*aPx&5zHn zDl0e3xHqtEbw-4PMe(sNHEDe$DQ#AJ%=RAUJu`Rn?#a8Eia^8oQkt^Fr2O~8g=rvZ zg6%3|D`cZxn{=b_{nAWFN7Jo=pTV6W!x5YeC&yQR#(Ez41t4Z|URGW4Gjv49OE^SU z4=w;PA<%u`qP2Ay9YI7%a165pyncU}4exb=0lw2dHFA+ZArlNN7|1dBAU;miyb7Y* z1Za}0CIv*2owVI}PUgD#nm1n0{Mv#btgsQPdBJeRe;Z=_^KSboXCZG=)1E}&zH>N6 zpaRng$cgxC(f~PK{o8q5n2547mEJmqBhfnto0_tyo^!*~|8^_kK|%t@^t>$k|LB3F zp~^5}*TB|eTf)Ip_;>{>9bd?{w#j%ZS0VJqziTeJWpHw&dLAu^$;gx6o^qqYA9!fHY&cc&jo zNcgY!)`0h(st|a-@hvvLU@$&DZ@Om3zr2%P2`)5`q@rVD4AKe|^XKbPpNK0vV8_o3 zVEo}3rQR$4@urwt6iOg3@Df)OhJ|l8-0ge1c6&Vsv~L^?Z-3bQ2pc{MCWSI-dt97% zIbCTSzWn)s;B@ zWx@cv?<~-_7}oKD`6a~FL=NU`U+_bNI@xv0r$h!%#re6`>(37$`~gpjf2`;ZIAU)7 zI@BPC&_v{Ai@_S^O+^+y9cQ|L%f`QCVFsq8=mXOCAu&GB@5Jm49PI}nPEPxxGoKE^ zLO+L`N)=DGE z-1%7kjfc1NQh|MUii*_BjSIZ-I`B1;wwKTR_9IkC!B`x>18_-{vXYl06^gwO#A~Z? zp>hktsx@X^0bs4E@>*CWECF`Gb8?*e(8k-Z?FUQ`MH3Xj&M$tC0{hIp^4H`*JrgL( zB&=ccVJhc^<1v!d|1a$(3?~_at3VOasHPW+8)~V}hQI7Z@ZsHOagz_TW;1Rw`oe>| zA$a8lilsCOV&G5L`Ddxo^9BnLHf?JNRP6Y&^%>Oxo1)<4VsL5WX5LhN!b~%R)t}Cp z?~(d-QiSR5d~@Y?0(CV`l`vG6fKs2a8Qv=z&rQCE9mZr z{irG^m8g9S%05hQEH=a=!i8Y^_8xq_@khRdUjVD!xw1UEazLOA@=$S{xX%W|spoWb zbg20(j9foi4#>@xi~~?2ufdx0z^M(79IOtLNDXD*wuXZ{3O7K-F(pIm!+$`juOi)V{gl@Bg{b z&woh5=8xN@_iEG6h06%z)Uy6!R!#|*AdZ3ZbO7(%F_V`+&KRT=l!)8relkATdgqh0 z{UO1;IYn1v1+%5vnc3@~%@3;H?Qv!zeD!5%QUAt#Vif~}VWmszw9jY%w*H?6;p1aR zp}jY6hV0q2I>VA*18qijX#l|0@652o%ah2aak_a}Tb+pBA>%xC9GrUddJ!MJk#=W-jc1No1VHY&nH@27mI?OlK zzkWNG*REx%wEVvAL7rmn_&vQ(y#aY8qq_qe{i#d@-)0BDJyAyWH<^7~R&X2-TV;#3 zn5jh++lHqjIHOWPlqB|2|34P;YAN*DJ;w*A{{i(B9N2}1%{t-y-8Q48prAmAV*@Y% zzjpZQPQ=uZKpyOB&c1v7mKaJBG(|~yJFDAEY~~9zB@k*(EVK5PaamgfhP`l%y*}_# z{|p)b?fHaH?*a>V2F3O8Q=mUrR?0ThS80IyXf?#vlf1VeCHZ~$=-G8N z1Zg?^^dIa_@HPxeKg)qM`>QDcoGH>;z(5j?l+1&i__=e3e;?5k&OnwaU0v)yEK{{3 zY^Z}c_xyftmg;uAjOJT@)sSyK(9Z{g1^l<6fOn)qmEi}xL38E0yMKifm(FPxO~|Sm zQ4_ZNDp$`G`KG?79kB{VyW@7A_7zhL#p*^THYznI_l~`tA;)%?y_y(s;pWsenO{4@ zoUl_3UNNjMc$y&wMmQoqEE4m4e|86e!4cgk*+l_Ujp6>EA6?iOH!~El?s*PdgQsDY3-G^A{ zgP;B;H?xa>aU||eN4Z1bl$`Xd+?x__meU&v4dU7&TyT^yk2g&#c{m9>*HAlM+hk&Y zJFj{)+xQAzm0V-nbaq#_SUsKv&-b8}QFHn9_jIg{&B2c^#jIayfRqF$c4rC8`XK=; z*Wz68(cVrZiipBBa7ZP9mpj`PMZS8{xUEm`&R(@qF)>|hokSc5?{X>j2J?nnradWO zRK4APvWlhA+3mi)fs3SZw@s?Sq`7ev0C<$ZR*aVmH-^xk9#4YU}CvYjK;_*q7`Mo>q@n zx#KVYu)e+iG0FpE5G4p6p84kozy9IjTkwK{*10=&IqUFQw&Ikakdj}|UEtQ}DC;@3-T!}TSLR!KP?XKkYJ@SFP5vi|||7%$4a)<`jY{0(Q! zBOPKCWwyjmzfbZj?^DV)kiT!*8M?!n-|lnC4VBC}{HYys3QmG;*_0kxN3Qx6+)!O% z7!n$Ko8D%SRKAdZ;J6G1J;9|5kRV;4e9yk|%~!3539}lt#r)FJ-t)FIuXSNb^%jf^ zfySVoK~Xq3TR1;`$D`TjtA8cTgYuW(5x9RbS6lr!8Pc2SZc?*HOLaq26f=zy%*}xX zEO!m=;hlR3+Vay(B-}rL-S9jNk{F^09_FwOA>F2Zh?c1O3^IfNf`ANe#&?<;Io!Tj z>g`u&;tsX_I;mspUV;bhtZy-|$bQ>d&m=HXhk5^A_Vn~U^!D8P8M!1;K>9c$?ZAig zf873`P$B69d|$&W(c3rfv*6_{rzs)vW(Or;OGMw~OxvU9jmH4UXa)Bt-9x)vU}&7R&>GJHqp zrI=t`hnu8-hL`V9KOva}P4Vxyjq?W^Ys*(XW%Nn#e0z&@TA`1&y(1%xp(1MG zp`3n#%QepgXj4VJ|Q8xJCAl z{r(m+l>&v^f_RcE!~^St+Kh=5zvC=kIkoRpSbdBaGm$*tV1(u!mIJTlO4UOLJ5xW( za*032rZaYBk4AV{&+1?Xf0D`zBf=j zeuuPT-h2&#;p^MehcyHEELrc10UYBzJRfR@`do%;#A!b72T#&!3GDVej+x@1D4?^?QbW?z*oCns@59n8YE#HGUBUD@na%-#M31|=tJF?{MyE)3upzd z2I4+)czG?*bhT0b`#v)~Uj6N3m`FXnE3v79nTDZjU@XVlIR|60Eh#cf|4(v}YB#*> z+owYhS2W12IyFf8meYf{7ADv?9q8Y0^~4d!#{_S~QQ2U<(D}oppkp4}eP%sA-Dv#R zkIFXoj_vFJ{+CPXoPfQjP^Q0i3ve5lp1hsE-Fl;CE4>89KF+K5=70PSUV$y~@2|Jv zsXRQ94`PY>^dG-N;5bLK!Rq{OEE$+!_5TlbH4>Z&fnRp(_bVK^ssOL5lkR_Vqif#y zrtpGCaC&V1_`_e<432}er__K)@;_LTgmM10y0?BKVdY)mKI-qi`VW>Q*chYhzyB{E z{eK%{!&dzN&BoyMtT(!{zAviX@WuFrraDrq?2HAU;*VCI>Kn?iyydj;dk$;JMyclX z$i(?aHR&Ur3ahDNzgMZ;cE87J#RN^;6rR4|-EgLUKT5uko0XNz>}S||;|O!W0Vrm8 z{U7`P8A6|wGkaqGquWG&3~b8(9%jQHT(?{QTbZqF=l@pbe=D=`!2DliXIAoa7aU)P zM<3jww$pNF$#bI4k&n0!@Sf1U@Hqa0`~uzyg53abwZdYUiw^jq#!zSM=fQA;UHb93{lTEp7YD+#C8V>?U& z*!jU%mm@;{!xvN#JiQ(H>{rY;g^$2;q=1I{65!Sg9b?o%L3C$YP{77dycRVL zK$c&hnfz!ijL{1#KvUXiM10%f8&HNlTy1Z0zzSD!WnpnL@pY!0|E#gJx_t8T5R-ev zX_Ms+Y|gc?o{Mkert7i%OuxOc@C*Tgqj^<&jm$rbW4~K%u9~O6YK1qfDmOB3{66n+ zj|fv~IfsGNwxr#pzg;2tEktwKgVy>iYd*6Yg9|ty-T6vO;C#?orJn)n0_-A;H!SG` z*E3s#_6=9ghG?bHzdd!QSIn%xbr8#Q>GYIirxe-Vp2?l>Eju;-fjR%YN}o+f;>XlK zGJfja;9WzmT3%Uu;3JZO{5*Cera`gKpqI{F1CzY6&u2P`9nP{O7wx043ZOg#a?EnP z-=B$`YN2L|3b0T?N;(zX(uh|k@m_!zhhfiCj~HC_+?J-(qsmx^h%<&_RrbG^`-OwI z!;A8MrPEuRK%{-~f7asRpo0O4llPV{YM7206^7)zA%c zs^&0Ii5&U8?J7J2H4yXB^xAv>Ed37LxAFFRScl(l&YlAI7HgQYSZ_;VoNy-`dK}LT zdMY09ie+Ut2pp?dZup}bC_;AxhwJzd(fUASpBg5aDX}KYwf3*|d#*=1!eqFnXRe8n z{{Higq4RYJS{gXLFNyqdqhGt~6&x)?F^T_3x8U@#?^J43JGFPM8WM2qM-9!}j0lce z-SR&9(KM4~X>0wOBK5iW(zD+MbTeaO)Tj8D7Y|Zc<>!qQn9cO>RQ#T9rt|(0OT3iM zuN$u&8_aEr>g{o=iT)PrwWB|e6G0uM&at29HnlT1Unp1=>ety*rTX0K>$7nk)_zx5 za+|4R{mthRYhEc$^=pW?`lXdouhvsUr?q?Huv1Dr^BH_K`)7o%K3{pR<1H~7$p4zo zBw8`hF3qx4H%GlW#yiN*Gp4M9Q!rwT&@F~p&-IVnXrd<3EtjgwwVM_Q;|JQCTk|gE zl@=((XS3RyT}m9{q*RFx!__>kTYgE`&PpuNLP2+NB2>9Qh)Ui|RX^;Ser;|~9IL^f zY@M21czZBs58k)`Vu2^yX``! zZC0NTGeSWvbVF{2OTFxmja+y|7{|zw$_(b9`hiAX5{)KU0~qO)fiiG4)dKB+X9Gw* zV;7sR>rZj0UhfzR3|ji*D*vosJc&;jh3C^WUZB-t&ncBOF`-_FCgrf*6kq1ZRPg8kDvVl~byD=`A3Y!I-aIfr zZzEW`_4#XI)Fddf?`sE?U_U%$)igEGY9xMDG@=GW&9&GG++_Zqeu8$^^M$m_&sd1R zJFG4a_cPNJ;w={^H~{jprD|86{P8q^lqyUvmQGB7)BSQ<;|gLvtp%_oFe@^r&HKbB zBntJ0ZVfolWQQX+A7TqL3FHv~DsI!tYIy%w8S_h^Wb-h_@?;{WS(V}gJ?X30>fFWf ztGMy|s#4#C#XR{g3PcAYulUM%r~}a4rQ4dY`!57<&($)YP(>eJF|(8F_hj^+O@!x! zFohC}An^5;Aq$wiqF~f7Rw*$vsV~%p(hqT>X=OD(QlL$(pW?KkDg8b1{dJ)o2goHp5mR@H#io zfL0xMM0_|m5dLRMe&Gk%9chD6PIGC4#Kab7Z^wmeOM`Q;S9A{o=bm!@GZyqt)o59c zeseUA36pM&v-e8&gWHJB7fs~LZfIN_dVY+_r< z)BUnWfv|0pzYMJ%gI5B?-73d^bXn#%zQXbbdnHD8nz0_nX(>IKvy`mR5Seg;`Dy>v-W^Y*bjY%^h%Iemo#_T)t)5&?704Vd_Q(gOu9HMQY4%%$%M4U`N0 z6^}wrj2k&P(wE*5hbFBd^m)}1s|K9OcO{2qS?dP0GqUsUsj_T+3VB6^iIYwyUf^-8 zH=-K8mzX%tuP90L*V>*=oOrBl2h99Ku@{q#5w$o_GF3I~?2?nai77a;`+|Gv_%poE z-%=0Vq$!F#3+=&UYMKC*er+AyZPnvT?+PK;-(_d*aPgeZ7o4uS8qEg)t#my8WA!r6 z8(CHZ#0I&l$+mmNX)kdwQWOv}jF%P>OIW=ENJd^$udda{cNjd2pWp9O$(%hH%0YOR zVIVwb#npuP*WtMdMwp)P!jYUqI4)V+dcmAZ-A>6MJ;Htkzw43T0FKh)Or3Z44eNIS zlo%qa$61bLY~ApWdk%rr~3&tX3wL&ka)lzu;r0x zgtm+pyhDqQ_jPN}Pp1(N=45Kn67#~$LdS|=6|3Ny89h_CH9ttO2%ebbMZWCN{Tv>h z09W)U0(RLm5H|qjCU-=YZO>_xm1(FLr`OjY7mLVZ%|4y*ChsMCmk5)Q*-DA(vu&gH z-oQT;$ZQi}-XMV{wMCz=u9*CphTR*MX!EvM%^l{gNEI~frp$niIRXt93jmb2A2EO; z(olnj6rj9}LUJ!X8Ot{_{?A@2xD8^iLykYdngt;4*)(9|jR4r|WXuiVlIpJ(A(qrw znQxTcgK5I9nF9s?Fy5^{YCK8xb|XhD>33crORhnZD>OiJ7DwL%47ovGN?K2bI{H3W zd`86fl9|EsCFR6T1M}n7c3T%Yon(Tu7%Hx2Bo(?J#w1WEC>>qkR+2Pc4VyH1&Ju0C zTwabxNeo|qT*qtui)ZXZK}tfd*?|o})!&cM)yrU?y|7!80^z_@K_Q^!a#ZK$@3rX& zr9(7e@V|%om&EXYE3;m-{(Xo3Z)LV#8L`Gv4{^9v6>)9QuFRG&5zRE_FLPGFs3(P$0=$hde{IQ+p4D2>o+|}x$F`fMDiY#rRVJeJ z;1OpGR(Jtli+^t1mXh@2l;szDOL0IhWD46JW$kbeEdu@LT9`txrt(5)0xm+xX%14% z9Iv`)yhBl78?HG%-JPq0`#;TOWZzs%O)xj;A#b~3DG@-j)v5_=-h@)6>(Jon(AGI| zp%kRpP1j}Uw=O79H$%*GY+;_y=3lJN+7~HE0$kZa3-7vJTc8K}OHw_Lf73;iP` zI>J5#*4?S?g!%a`KQ!h!z=yd~67y5H)S4pL!w2r@=h8epiY8~Ug>s3CMrlyK#BlAd z%U-D~HPj}PgDDTp`QZ6w*<5w>u^G4VM)`%UGvaFt+v9fVYCkQbx*i>u*Z1C%vFA>5mw-8 zKcY~$2w?wH4J=@bHQDhY+kT-A%6~)8AZ@=y2%=#|hx61rfc&o(hPzG;*Z{)>(uhmR z0X_|8(3h~KXmu{b+AQ74dQtgM*aRaslhtmshd;_TPT@<~&h~v!J7PbCw_|}P&XzcK zO6}~SD*@YA->#1wzP(h|zLs;Ee-nb1kialJ&?Cg9@u~;(hkl-8*{@&~B(g#4Bt-;e z?D`;E6ck@b3uS=*@iMlm0SHF7*C6tfYaOW8ayz>}l6e%wLoSb6RoR@%Aoix1neV)J zB48(`9oKclG1SZ8a2SvlZ#{*nc}MBY_UeEXZ=4ZZ6#!r?RMI|UEJ{xGJuz}fWe z`{Xb%xy)rZq)mGn*gY;n);(Nz%YX{5=Li@Q&Ui5{E6`oms8=Wvf6e z7t0p(;b%hIb|eX?Rib9;VW()m+5M(!ZzI!S4q&F!9T#+BMeEnHpaH`2YYs10iv_>7 z{HEXzFIyVY;E}htRIUUY;-f(Y@)a5=Zlx9A65h-`+CH6;83BR0iE!^~brI=NBSqs< z!6->TuO67}OM@CrByba)j!$~O+kDUti3p(fjIlcxn$uRb2lo3~7(jm|=_m-ZZ*gv! zFkh_?Lv;zNibi_ocHIvJmL_Za!RkwKklynjJx|Eu;RP`JA&oRcH=cQW?~njvrp(&E zFrBRIkgE*YBf*hKG{WLaDESX(#!DI z77Ju1o_S_7iP#7PR3!s=)1!F@vcR#pOvx)==PW)QazXR>?(ZwKPfXoM4^>RxV(TdD zqS@E+IYyQ#|65$nc&R?nSY1LGoWghle2@w^p-}r;6ig%WPCHu40;P%3 zK*6ym(sj%SXA-3y=4+cXxO4j@83S!%YC)RB+a3rSqWWv26irGS6qR*6&7fc$2^KJ} zXZoYjv>??CmKI2&o6E$yOOfuCkvV8&Tn6nROH>>n(G1fZct?{l-=uooN&&Glg2HS{0eOD39~DOz#Q|+>$*j=p!DenXYW=$(Y%Kj zfm9r^gu&X`9ndS@s~Xi2|=VyUA6F^$=EGikgg zu!qXjsS1$FU=%FtkBLpVgfuv06bT0@7{^tw zbwI}MkEnmid;Riv(d%0C+j{|RveN(dSH$0ghYFu$hd_QiX#F2Kob{J|M+>zg< zW&&0`4df-mp^p$^9M%45CSRdru+=GA3->T9n8AbOEIA+(h3y`A1oZqN1op)l4nSmG zFN&QqOvs8Z>5#mW3oGtQs5zSH|2OZdN_I zrorQBbITv9uKj34Yu>(vai(v4PU$4O8eZs(?&O2=+7sStFW4HEkZy*DU5N`W@*Bm3 ze-F>Qk>zq7yAv>zp>gLOel6%kDyMO%A9#Gk;mcsxLa(q4Utr;m8lr`tiOv!kblDu# zFOeOoCdN=V&{LVJo;UBJ0mAf9jRQVdmSuVM}db>WL*_b;f! z!BsA`7`z11AgtX&!o_3SOUG#EE}BzN45q4ffWfUn=9S$3xz_y;(jYC$XrG4V`oL36 zVnll8X$V#;+8vP(H)r%sH7LGn%b1eSG{$p^6<*esUVN6)?WGzCckaM%uI+vGEdzCHUr*IIN)Wg0coI<7=wPDmQ$Aj45Wgy*KK@#2E>$9sR!5tl zP(`2~h?xg+l&l-yw#ze@oS(GE((HV$<&)nhL76*FFCFNQUxeb8p*+eG7xjaPxFfYU zp%2sEF$rBK9X5{csz;v>vIG)pHw{qJ48IEu2?QO{_e0I>f$v^Prkc9D(e{2TS$gBp zo152_H3)jEj6_1!A=xgwnMwNYVG!ZT!GtqX#~!G{cHw!ZZfD&sY-Q!)Zed?C=Kz{g zIu(zvapVF0y{<^6m>6DB=uUlGXBYAg9_^1!$VID~Ogpqp9a6Bfoe4J}9b64WHK(Ed zBY8&Oh3*xp$&=<_aU>T%l4^e}$=Qd7L2IF;r5fc_I)pq-acd>gmE;{t zDoo)W=l1yK>^ECOWyKrL-$R3D;^mL$f@UqL9WO(Y4$k zO!gIS8)~p@P}Zd=#4T=h>7+iFIXfnhXXK=j#yx(Fz$c1YieaX5B1T9z)3|<0nVww+ zPpc`B;7o~oA?a~rIYW=JWXxas`^uu|2rT@x{oYR(SL?N21c? z9c2hhT|~rG>DX2p^%Cl}#0l}7UJKjOmx_-*e$F|+`)Wy<)g1d22lJ$gVZ3~F+614} zq2nU~cWGUCO2)}MB8Xu8ZfK6B%B8b&A}e?zM=WUAnkZH*i5h{#M^AKKi8yKG&?~4D{$LJL`i-}B>vu|a8AULLAmP?LF zr8c0qUwkcvdMWJ$ZpJ`Ssn?(ym#t=HXEERsNIn+9EPXeSW6UO=T;qbaTk2xPcb}&` zJ!O4N6;iG9F^ycfg~~^qday@LPhRo4I&o$OC_U&LNkp zQay#+-Ts9mN6#GF9yRx)5)qCiXDK8x+ef?ETntA%BKI4fs84XvO9g5zCt2ra)8UX$ z%w1rcya-inLu^RM!pl(SdZ-F34DAz`7!;to4dao@nsE$YG;c(Ye%|{^_cSJf>hV{r zIh`Qe+j%PV3hfz)Ku@5!W$7?|%RYf0hEoq*jmL0hCwX1Tt60Go>1TW?UrC$biSzx& z^b_S-4l~M4Lvzi3I2c}B4ml~oJnKm{EVu2#B@V{B8jL-8Scz8c`^opoc-r@kFtLD~ z0)^k8M;tw>9`Bh>G!aXDjNE*gdt8{;;%1`HGGHK^>QAO=hO}ttIBaX zf#WiPK@OEtt>e75`>6GC$ENQs7n;PIx~t=M*OjYby2o9;e(v65gsr=t;$iOLFpuWS z3y$7X=NIZ0@Gz*@?ltoRyHr+4D@N=}5H5-hoJ8(CwcC1yjIG)T#_@}L?&XT!q9t@# zXzjqpgPig`E)S#U*R+qS`x|iWof61ywm2mo?1x{vvf>^RRAx%v+e;yEt&C0BOXH%^ zjHN&IsXTQOp%yAOq7l&iwvu&ePZSlSB^&Yvz{J%-3j-rI%NE|gv=Q0r_+ zUd>S53{*Ntww>>m=RVara6!ES}JO60# zH(H)&;z_C7UK!M(D#0a`Qe}MfIeS+tSq71yP44HS>;JhegSUkz=BB+I?0+6VLEXc{+SA&A&Nq9@nJ6N+lxj zf#et=Th5(|wp&rgoJ<}PAP4ik@aeR9jO&9%yW!Rj;{dZ2!3tVIdyaQRU~K=yNYR9ZQJGC!$ubACVR+vOc_N15RN27LFCdcgEU{sU zIsb}Wd;SPf#NIKZ7f~Qe-Roj!|6xMXUoossf2L1tI6Xj^dw0>7`$+ta>wtLL3WOG~ zsVO3mi((?> z;xfj0;zG54!q0?Np>EUiJ_cZGkx~IFN0=ewDZZycc#RcT$u1%qGNGypCy^n;@h%@? zx*)yW#Ca{%(?b4dV#0F3h%wE4HRy3sJ7fvWhteEp#05*l7n|hSaJcx4nS#P}aF0m@ z`v#Ie?EDPDw?*6F&ca#_RByv8s66}aY1|VI;wDJceONkxfvnhJbYd>QYXTC5PZ=fr znpQY*Uqr(A7jKX=?@a1){%R|g9QWenaxc(6e6;rCTu?(yxfJO=6FV0s{B*PeqI>Kb;ZF!$B^;2YwinP(SCuOQy4)b81p?}!AArFW%vWb|w|jER;H zR1J)VkL{@q!cl_gET&cOC%{K97wVjYlJyrIh24esq`!<36WP(zePHj*3!i78%i3+u;USuNao( z<}D<+QIhH93hyDKuHyP!HV#B@rOSM@0^Ma34&5S0N)`h-+!gtWJ$kpf&_++l1u8O+ zk6Wo|TS2zV>U8btltuWW_;~c&L@J-sk_4K{gEUILvcu@iw$8_fALI_1>g?_em?6#5K zL#rB?&uL;4UnONv6kD%uwN@Auj)3nR;mirez~fhG@)CW8rSRg$k1iT;9?jD$cISOg zH#_x^cJYYc1YKG5m=MQtaURd@7rGavR9RP!Oo3D8tk1cZHTy2EnaxyEJ^`w_BCOSG zF|P9=)$Igk&x9viusBMnFwIJMvkh2UnyuRxm{j}cf(YEYwiPcqI ztNyhuF@uBUYgR>dHgChp!%@DDR#{j0^NTD?Z^_ULIdj(9uTi8FT%_{R4&Q#(0p{h8 zB|*|eNm26mPtiPMez_emHS-+GwEaL2Q%%3W417W(?Dq#m&I? z;m=h)Sn}`AWn|{E*bMR|7SeCSZDafS#&1^)c(&ullRXYr5RAKvR}i{Wewm-Fwe)q4 zJ)l;j@_JUVOQk|$C13bmdo*3`LNiMso(e!Abks{JUdEJHQnB$FD%ZGL6qb)>H(gah zT(+m2^qGxQE;<&F>@E^m;;5qbS#B?zHWsAeKrGMeZ_M`H%+(1|{(5$hh1r~x_V zoK{0&hLV0#^{1*EdtcrDDPUE7ta{-S9w93I9w9jwW2+a}*9_T;Zhbsw9)CC4sYWI_ zd&09lMSV=gw);U^&p`C~;?6K^Igy%G319p3T%P@sgqf!+&xTbgd>~dpNtw%E^<+|1 z2bm#(PRIQjlo_(Ua3a%tD7;vGs0#o22+-h}!Re(Ge?78|vha}GL2lMonr%r)gj&^2 zSW!?;E7`Q_{2m+X6}k%02epIpSe4h{l63hDS(U-(?JPa>Pt<9bBN1|T3 zCPY;GS29V29xOt3d1WjjJW-COR=~cA(Sek4luCq_29&A#5voL2x6n-Prdi)5+vaA6 zl_i%3gw&HjTF%Mp>OOzzT$*ob%u~p@yMMx_ZlOlGy$Wns+llqvGgzv}f|({PVt;H^ zYGccmUMHjfT$J)tYHAhr1Lo1gv=d*xn}*~qw2dY(cqd*UcD=7&DvMVRUgQmVl}p;`#=tW#B&^>Oj$k0t@Q8#m`p9zi~8-V^%4eSMl34*chS97cmOf!F=6Vi3pyWFaTIS?^lM5P5qf^6Xa@?5P76Xs?WwQ|9t@S!oY7B`a6Kml5P!D4%h` zk&jo`nVYr_`h3|K6v3qtLR{gv+d%KlEvc)qBsr3+qg`nn_|R?J`Y0` zR5UU8HAr7>NZg-bpD32rtjjq6wnXPN)YllUn|&5le)(-@==cx1?g=HR`Xl}np^zgZ;yc#OtE@pl+ z;i9of=V!`fEjU#VMIEK(K5?NIMNhliF6x?yX|>ief1zzie1dv}J@9j>J`~!^to?@p zNdstnpmXk~Ggfw^odl|;$+!_cMVc#5)%BtAa8yyBim*sc1Dtr3OV&J{W?p}!9~Xa^ zgnv3s{kDVM2mw0JCe29#$gev4#9{VtCzbENhZhh_?qF(3xYJBuE4j_-WHbL z=m0!|Gi=rUu#oM8eu~##6nrFk+`iKa8$vjk`(s}fVYihGca!GBxEZB9llZDWyL_k> z-W1@WtEL=Gm&v(6o+}^-w5*LSUlMbW1o~lz#W7`R*Rz=ftfEP**0G{XW8ci5l3yKJ z;skdyJxINs@zA@h`w6AJNfm>$P@D~5n?El6bp@|$C9`{!mg9@T=@lX4Moa%b3G4T} z=x-WhkKQO_47>(>*il9{R%Z+E`HGc^j2-F;tm}(uT7%ZSBBQQ9?eW3kLycZSjtoQ| z)JejN&qzgN5sPng(jyKj+nnYv9WXeKIE4kk%{QxQy0GQl$8t|NP5f{`-FfShM;f+I zTr{;7tK~dtzU^mvaBZxn zLNk(jJbUq;T)6{fmC!>^t%T-2iG6H}d%r82C$;pS#dDJR8mqBnu-aHVc?iLh50Bm# z&-Xjc*KxPB#Qnxzn6_CMn!N^HlXhRNk;FI7#g>r-Iq7tJ>yhv^_Rn-LfbD8Ccg{p>?4wPs-D*)>j^oOEo*dQ(Xj!E_ zb7EuasT~Y}VbI2=#s;>=HECi1i%9mECt%{nCGCWmcSQDko1gHpuQgYz&2f8h2D2|m zzT0@wJ{k=5d5$2}K?J5i&55=T2#yQK=k}?Lk}sM78wfad2bz#kveE*ahSdlTV9{#a z$m3eO?`N-1*wc|lQ~rp;#uX6`PR9kOrHugXDVJck4oyNJ8G$xu@elpagG@bxS%7V< zMbNXZf`dj$=UQ&Ql*F17;t}l0D*&78r}gID?21A8J6mvKH%;-xE^AR2S?BWX%#U!p z__Hz zIp5s)YLK4p-2&7ML+w}unD|+`4mb%H*H^q?v;G6LYV3RPk?Lm}EKek&SO(3>>zq}0 z0MaZ4*@!XX@v2~wlf8Hn^Nzv?a45CF+USn}EJFp!dC=vJOtkH0(*N@2?!#7U@pahR zq^zb~)g7|b`7HOYjh_ZOmIb(ol+GG5mDCWDuotr!yX0cZ%h0Jla3{dxIIub5f;N4H zlb|iq(10qSCom&-1a5!r`LmN>tp#p@@?p%@uu)-6J^UW(mO?(>4K4W1PXR@T2Gq81 zklh-^3ONCQJ-Ldg3OYA!gUwA;6+|neBAxzooBw_r^aOj7p@c@Ff#bY2h=N2f9ugWu zf%^Z#G!#K)Q^a=Lp{*euXo~{B?v&g=40eCQ9|}kireO^K&i$Wz{m)-Pi}UohvAMze z=Y%ks7lVkx4_keIt`DGS(Fnu1WU&??2}OeiK(I>I4~~P$LRmk;vj#~D#Hmu#;Eq|E z7LkgR)29O4dK7HMw@R6td(nj&uua@|tt@fgm<>ui{SVV|A00*m)I9ZJ*T?IPfrqL= zJt#m~%RcWYFNG85Pv__r1sz?%0FW!0Gia+qVL9q1>w|{)Yt;`<@5uDb1OmZHrNPhu z)m{})!H^j`w&^1kd?aPZ>-w10Q0cM0wajtp1Mh1ka*4AgprAB&y!AOL#zWe+i%dUv zs^PL_yOk~wF6lgP=ig+foHUIZ$Ld-$9S1*Wy2o$YKvx-mbZC&Cnu!1P z;M6i0VAG8t;72`hnD)3WYJ+sv({jZXCTFdun{{2ZRj5@k%j)jX^W^Sc`QUnxJI&Q{ zvn(;S8JwYW^p?WwTL2@5ds?^h3Mo!mEgC_w017`_&INw3bf+5PNj~DYPeyV_fTTk5 zEY&vI4XjGAy~Vr`(kQ#Va-Ds*ngSP<(x%$*1Wc$s?%8OM;OSdAzA+VO2Pr1eUg_r?PkPYLzn^yvA2 z-sINhuMky=zJg(jYvOwxe4VP7FsM;`XY6$HHnn zt8+;p8O4_3cDCzr8mjaU))NBXzVSS8pT5A4K#@0^jG_`fDdq3eprq}ZhH%$FW6$nN z(>9yTUpa?B;MYo^x8vKs&^{9hQ=Z2UgVstvQaD1wy%_>ei3DfAuC{eijsi@ew~FmM zW!Aon3{e4r9gH;#fzcW5?LvB3pRlq1Tw)?emt!Lug=2VAM)FSv|LUAm4dlK7P~mhj2tebRgTOl{qgpUz1;hh6T?E z;s7wX8W^}8e&-ABG+$|@R59;G-bBDC`5>T06wwjIPUTV31@FT|@f+taP4}Y$t70vr z3$l0Pv_xDVQQMEfwAP(O7KG!OZ^+E^%$ut)iB!7n(E!kX=Cdzn?!_uDBe*+=yx1ns z0n4|%3MekrHmI6hCuJ`uEQf0`T0hE#YNtLp!?M&=8Z*+1F70MW#N@e#NvEWVwQM%1 zwY<9bWWv?48G#$HjvBey6N?0I0>BmPL*`bRno@;PckeT%ZM}vI_F|g5R%S)1#&43VJy2V0TR#;hN%qKxt=d$a5QB8?B;3OV42{# zyVp+5r2g_v!QHr^2_wYxK!fxLEJL~VP5GVq zT12S?ujHcovMA1F^|{suO@j7wBtmg;KGe7-JsB5MZ)ylSZ_3qizF!sV$litfJT!)u zS+WXuagzP#)FI$Dg|R^_jx1idK|1}Zr(Aa5kI_nGwppmQ&8ueJfY@0d{tr{kQ%l^i z-W#d?5TGzW&`6fX8=-rw7_egkUa4u*E3wyrtr45oaGY58$f(|jdn#%`Q_;v}?se(8 zWP-p+?)!3rVpdj=z+Of~+S8-Zkzp|@bvjj^S*$zt!{fnpr)4ME%a)C!^P04d#jFY!&v8FELIfi;G3zFZgBQtdtWW`Y4 zq;_5L``+mekwdbZM)xrmG9XtoW<+KS>oTxwTEmvW&1RtKQ zRCF#dE8M^@6#S zb0Nw|(5_ezd4#rc6wrz`e|b$$eKwZx{ef&=0S;G$!K9BXhn^KUO73W@;w-4Kcq2EJ`vgd zn*NT@UWX>K@K5j;VpWNVI$zTNY;q&?#~{{8RX6CiO4AZd_HmC_z#iQO1@>+&*>Vm_ zftCXiKb`OKRmN zKB9iW?sTcV7D?N%aV4cf~=MIc(AH;Rfer>P!K7@)48fr#-2k57l--Z)h=wa;H<_D58c`6gB&|Ol? zghEfQVphn!O>p{Ob!X@ z4JzIBatVoABK4&JwAJvez0c`%F^)U!=2ScBW9BIK5vJlSwzoqmqS8rurDE3=ubayc z_~yrjDp}oshuv$e@|tGH*wfRM(Kg`&#`Wi#Y)i(c4wnDRlSa&wZLePEN`{n@09W92 zeI%T*tgeewxP;y3Nk2a6>_uLYvKd7JvU`qFx9(sPHHoJWvQ0UwRxW8#eQqbkU)w|}^@CK#9-GzK&((i+y&qQHbE_%6ZH0=1bq@j9>mcxXlMwqUocXpi>flgpo$EP3BSES+2j0Ir z@#yS8eVdW}o+U*ow!B#_tnini@O#9e)>Ap%EGDlrUso*>OsDZ(2LKFXPiZ#jA*LQB z-T78g(V>#t^T9q0AOni}q9muW8P6!F;A8;^gx?{0nP&Xk!as=; zYNcaV2==tcCKvtMJar^z`@>&acMHPOy!LfK+L4KNASu|SZZ%vUO>{r4OR0X_I~8+Q zz*VKE9PBB6W;=WdqFAZb0YK=A`XXx-SE8u#E)n_?$y&>)7U{9qa_a?%$^>%9?T!U2 zO4Ws|$^g+dApy zgD539l~J6dHVSbk)A;<`abQkK&40VcYR;!5p7)vDZMi+oZ{R#eLLorJ8>V-9)XDpw(~a8 zUz2`P*s8{iodLP{_~_~L5M{Lk7_*a{cfi4-02o34deOKJ=-l8VO!CmJLe(hdEl7^e(XK09BzZ>YnoN^J3Y1Tlc=@ots-4Px%4e?iU z&8`VYQNB=dMH(EJ5KdG%N*m2v_z-FK4xcLAriK`l z7@8Ld<-o2MW~nAr=xPBf8mt=;4)b<2p0qw$8OeJ%)=^9qG61tF1LTCg`e&6hZve-9|GZ_s901zm1!_OQyt|br?3bbOEI`~ z-B;VKRJ=P6_cTHdiB180X@;p`LFGy3bVi#+D-t->Js0->w0EUJO0ur!nR*_U7KnO}&IvP=tn6!jN6BRXSzzs7}LDmLINC*@(6a+^LlxPiD z%PLT7WQi1pA!fiLy^k4|s4)H8nf^~+zP$Us_wGISob!pzHN>A`;&dcwH>wWpHRi`@>`$Ha)aBku?bon;i(M$ZFx0_rYzv z9B-wp$=Z2LWE|gUcFj2FOs)x{Xbw`GT6elCFT(k#zp4}XYG+Z=XQ;}$8Dvsw{qAjC z#UwN~C_W>$O?_JfLWeuAst7fQ%2mMDjgEk{gZsGo!Vt)^=+Z5RAK$@x)Zq%x{jE_n z)-&ckwcZK@59|qWn?*G7X3%X}UC=*yL4XhnR8fr&9@_%MdrcHA1u?AODz;%1V{V|M z05%CI!jSh!NNjpHrTP0Jg3e!vyo{1!?j$7|zIE>@R@nHgHv>*PVICdE`KAhJ_co){ z5@Xt3higpC4NDE2ey!Hs5*4AG{qZJwO&z*zF!-aF_Us??{WfR&CajzCaQsDNG=IKd z6OI%a7GaW-Mvj|4Th&|CEof$zbi1OC84_#fllLFo4Lv2B(R>R9yx$kl25K8e1CW^* zJB9N5JWXgeQcq)V8xCZQ{skx0U>19z)d zKe>PF0KCQE}kKEhrFL6t5)==2;J8`jJobRH@hY~wGGsXK* z40#1>v|?xp;H+i_{fSKS`o2LJL4`p#Sy<~Sfj zZ6X2DwRG1(QVFnQEl9d?#rG|kK98Vu8GziE?CicrZGTiKD%)m=wbc|@n+hNBp+fv~ z?B1c2a&tYy18E70K%{5@xqA$&5w;bXONQPFzyjz4plXt{kB_GsF-AZd&_v|kXGNw{ z`clBQVj#puLa2CC+*Y8t1u<64#cj<}FN2j<0vNJ_*jrFt<0*GI@(Av4BuM1|VqAUP zgKzoOZTLzWfH96>b`ted51{%rqvmVbVFq``Ar6`5rB0AVWh%`#F|k>oy$@H!0^~0y zi9d8s0`hktx=>uwyIVHoAn5?m=qeIl#IWkB{X%3750OOvz{X+zM*SzuCoi7&(LhTe z>_4-)HIi27@*LmdKbdHBIwN_0My1%sB6xUu-=)T!9Gt5>?Ynf5wHzuWrbR~0`F?`J zOA*#T2PH_mENfs>ZNG^MeNOI@{8Jr>wTL)dqBMXroQm^Ql@0cM0Kj%vO}%{a0GbzB z-H=56fb=$Q!%UCC3EKH`n#Tm~3ge>5nWOqcalu+UDh`}v4O;TJB};igF~luvU&LWR5dx(anJ#mjfqk`+{pqp6 zUk0u=H|tzv@~V(Kw9t-Uzp$5ro+IMauKNn)atl2%o(Q3&KUB7aE%*@~fVhFgdLniv zq1T?3>9PB*AqL{ANZew#()hyzpF69yvA!~SP}B SRtXdeodBP}&&1w)*#7`3A$oZL literal 0 HcmV?d00001 diff --git a/vignettes/mnl_models_weighted.Rmd b/vignettes/mnl_models_weighted.Rmd index 52b136e0..39f96f0e 100644 --- a/vignettes/mnl_models_weighted.Rmd +++ b/vignettes/mnl_models_weighted.Rmd @@ -158,7 +158,7 @@ summary(mnl_wtp_weighted) # Compare results -Here is a comparison of the coefficients between the weighted and unweighted models. All of the significant coefficients: +Here is a comparison of the coefficients between the weighted and unweighted models. All of the significant coefficients have the same sign, but the magnitudes shift some based on the differential weighting of each individual choice in the weighted model: ```{r} coef_compare <- data.frame( @@ -167,7 +167,7 @@ coef_compare <- data.frame( coef_compare ``` -Compare the log-likelihood between the weighted and unweighted models: +Here is a comparison of the log-likelihood for the weighted and unweighted models: ```{r} logLik_compare <- c( diff --git a/vignettes/simulations.Rmd b/vignettes/simulations.Rmd index ce8421e0..56cb24c9 100644 --- a/vignettes/simulations.Rmd +++ b/vignettes/simulations.Rmd @@ -13,6 +13,7 @@ knitr::opts_chunk$set( collapse = TRUE, warning = FALSE, message = FALSE, + fig.path = "figs/", fig.retina = 3, comment = "#>" ) @@ -60,6 +61,10 @@ You can also use WTP space models to simulate shares, but you must provide the a ```{r, eval=FALSE} sim_mnl_wtp <- simulateShares(mnl_wtp, alts, priceName = 'price') ``` +``` +#> **Using results for model 1 of 10, +#> the best model (largest log-likelihood) from the multistart** +``` ```{r} sim_mnl_wtp ``` @@ -83,3 +88,24 @@ sim_mxl_wtp <- simulateShares(mxl_wtp, alts, priceName = 'price') ```{r} sim_mxl_wtp ``` + +Here is a bar plot of the results from each model: + +```{r, fig.width=6, fig.height=4} +library(ggplot2) + +sims <- rbind(sim_mnl_pref, sim_mnl_wtp, sim_mxl_pref, sim_mxl_wtp) +sims$model <- c(rep("mnl_pref", 4), rep("mnl_wtp", 4), + rep("mxl_pref", 4), rep("mxl_wtp", 4)) +sims$alt <- rep(row.names(alts), 4) + +ggplot(sims, aes(x = alt, y = share_mean, fill = model)) + + geom_bar(stat = 'identity', width = 0.7, position = "dodge") + + geom_errorbar(aes(ymin = share_low, ymax = share_high), + width = 0.2, position = position_dodge(width = 0.7)) + + scale_y_continuous(limits = c(0, 1)) + + labs(x = 'Alternative', y = 'Expected Share') + + theme_bw() +``` + +# References