From 6f221841d5f6c88b3bfb2d99e2500b1ade4a41fc Mon Sep 17 00:00:00 2001
From: merliseclyde
BAS
uses an adaptive sampling algorithm to sample
without replacement from the space of models or MCMC sampling which is
recommended for sampling problems with a large number of predictors. See
-Clyde, Littman
-& Ghosh for more details for the sampling algorithms.
coef
.
confint(coef.ZS)
## 2.5% 97.5% beta
-## Intercept 6.666930e+00 6.78164619 6.72493620
-## M 0.000000e+00 2.19641638 1.14359433
-## So -4.707724e-02 0.31431450 0.03547522
-## Ed 6.496393e-01 3.17566810 1.85848834
-## Po1 -2.158795e-04 1.43681569 0.60067372
-## Po2 -2.108569e-01 1.42994718 0.31841766
-## LF -5.437181e-01 1.03209467 0.05933737
-## M.F -2.410463e+00 1.83011041 -0.02702786
-## Pop -1.250611e-01 0.01212007 -0.02248283
-## NW -1.528445e-05 0.16796356 0.06668437
-## U1 -5.363377e-01 0.36658398 -0.02456854
-## U2 0.000000e+00 0.66272189 0.20702927
-## GDP -6.834740e-03 1.23207087 0.20625063
-## Ineq 6.704244e-01 2.15281005 1.39012647
-## Prob -4.181555e-01 0.00000000 -0.21536203
-## Time -5.480326e-01 0.04211726 -0.08433479
+## 2.5% 97.5% beta
+## Intercept 6.666847210 6.780785e+00 6.72493620
+## M 0.000000000 2.180831e+00 1.14359433
+## So -0.058437459 3.010219e-01 0.03547522
+## Ed 0.597431542 3.173426e+00 1.85848834
+## Po1 0.000000000 1.403581e+00 0.60067372
+## Po2 -0.154818391 1.495914e+00 0.31841766
+## LF -0.553771939 9.673109e-01 0.05933737
+## M.F -2.150787989 2.045485e+00 -0.02702786
+## Pop -0.125843889 4.837513e-03 -0.02248283
+## NW 0.000000000 1.656774e-01 0.06668437
+## U1 -0.513971357 3.750977e-01 -0.02456854
+## U2 -0.005590497 6.478652e-01 0.20702927
+## GDP -0.022655463 1.204585e+00 0.20625063
+## Ineq 0.667892758 2.098649e+00 1.39012647
+## Prob -0.412865659 4.739855e-05 -0.21536203
+## Time -0.521332101 3.495832e-02 -0.08433479
## attr(,"Probability")
## [1] 0.95
## attr(,"class")
@@ -483,16 +483,15 @@ Alternative algorithmsmethod="BAS"
samples models without replacement using
estimates of the marginal inclusion probabilities using the algorithm
-described in Clyde
-et al (2011). The initial sampling probabilities provided by
+described in Clyde et
+al (2011). The initial sampling probabilities provided by
initprobs
are updated based on the sampled models, every
update
iterations. This can be more efficient in some cases
if a large fraction of the model space has been sampled, however, in
cases of high correlation and a large number of predictors, this can
-lead to biased estimates Clyde and Ghosh
-(2012), in which case MCMC is preferred. The
-method="MCMC"
is described below and is better for large
-\(p\).
+lead to biased estimates Clyde and Ghosh (2012),
+in which case MCMC is preferred. The method="MCMC"
is
+described below and is better for large \(p\).
A deterministic sampling scheme is also available for
enumeration;
@@ -506,7 +505,7 @@ Alternative algorithms }
)
## user system elapsed
-## 1.391 0.000 1.391
+## 1.379 0.009 1.387
system.time(
for (i in 1:10) {
@@ -518,7 +517,7 @@ Alternative algorithms }
)
## user system elapsed
-## 1.301 0.007 1.308
+## 1.283 0.012 1.295
which is faster for enumeration than the default method=“BAS”.
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
_=up6OeeWT{Xkj`(2H_0l%PESq;fERz*jnvJx*7lAWQ5wy-v7=YM*-0DDpYFgYcx7%6dMR(Wc<*|)^t5C& zz08)>uJmg^k%6|BQ@t}ZhbLM3)PO5B$2%xZ+fmD8^GIk7?(C)KFjW=7!fg;&YzoCS8HQsqtQkivd}erEwYt3mk-fP%g62v*1`VyJ$Ex-OaV#fxL?W6~ZwhnzB2C z31z=CSEIMxA72oh*MBNJ3_SY^@L^EBHssH?@e0G_#yLpqeez?V&NP@jg?hq$eiir- z_Cd?SA1m89DK9yKO_a3DHa`Li0Tb9&aKfd{OaOgZP?^Mbjz0pz{16@P4mzs6ORv3! znD-BxtsYuye7ksw?btsf9QMEtBNPhGr6 &zZnERmfmbdo^Z3*Ay@ zFz6VlAXZ*Yge3_!BD 4MWh^D4&*I &d zJV>>vBIKI>2VwP~oN-@ybpX{%l~jO?SF$WT-_b+RZ2l859-Y|#Miv+Zq09i0pTu3b zzw(mOp2KP=s0Xd>7no$OQhN(CqR4 @v4eF#V z%A)xUryL~vavyzF9Bo=AZ;84CH%C@ZKM2|@F%q0(FnDF P{&T=X7az@XBWI zqkj`eds~ppHnB91^S9Y+v`<7kh>xr58wp^KCx0dp1 xBV+nDH~2MI(x0DZ*rnS_g^<~`>_@CGO=hZNR`u1lR 1^bpX!14XL)xD1uJcHz%;kznbT4!z1kCAt!TK>uf?9bHdfem5%-23jaAe)Z! zZhH4PFd%EbhYT*8(9d5{^A6fSyFi`Y>3#l}j^wW1V& zx>5hwunv`-S5r%28i&eySnfhYDt(Md)vC)Nuz>R$`Rb4Lk }wW`F2~`5Rf=;mqACr~OXjT_DQjN8dN3 zPFkWYbXqsR0z)Eee-RSA)L7pyc59Su7MVD4KHq+8t-V{}>Y2==$rJ7BRwe*H5J&6F z@8^NUpBI50E>hjn4EgnKGnE#Nm!)nJMK$Zn3XJ>ygSqJ_V{<|G;=Isb87pz%-YXNX zAHqPGKY^%-Td0_W1bX~PsiFN_x5U^|zAO#=x4LD?5){s=zuKBV{S+wP5#B#$M%N)E z@4@|)mlkN9fQg#P0m?XnYiVw*aI{MqUHl(q%iHvmM1PUV1?|n!^=KMfBy)`qCgZ+b z6(DcNd4RCR8!zP8>0tBthY-7QUWJy7`A+jDJE08p%aFJ@*dkxivR5yNWB#Nm5* SI8UGlu z1H@JcL#nPF@aDq -9Q3LfP^F$xlgy|1lvUv825Y}3iXY)dxffq394+94SQUP|Y0K87^?6r0@+UN@% zNi1%0S6&(fLNv1y-5PFl_Z~y0Lcpj(a% qMxC`@C+^=7&cgcq3PKC87TLs70Ww zYxJ~4{tmSHEs)7sB3ADG!51mMyQ!%&n(rFz02}lGa37Ont& !yPr zOe{C4j`itG_M% mfc*840?2>^xa|{&s zszS$x+?=OF&d;^8X|8E>6kS~8tq?ojr1U11djJ6;)6Gev(>1tq31~Ae#2oQ#U|m63 znY#ROzy(yHu*;W
P+=19n z!9uc37uJl=GaOXtao7Ncl@?vL5G;`xUY50F#?mCwp^FtS4;wP 87@_LA%_+ADBt zlEmnz +CyIVVe0o@FnfD1*@vB0E6JI>qxj XBI-bN6l{*O0MjBt#Q=CGop>E+`~9ux`hdX1pJn3L2l zqzh9;DK^vMEtmJDMfXLHR|wZ3z+?STy@6pB9}M;`=Ml5hgb`Mp_w)JX9wlE8Bvquz z*f(5M_?xgQd*6ZydJeH*!pn~18f l3NzutciT@$5*k* T{<3nSVf*=esBLax$$|l1CfyjljExjPPb!N@gayC?ffZDHE5w zt#iUw?Go{ll+iuhCDePZgX_*UJLmK6hBvevDHUF{y9pzhBuN=1y;kt2g`uR3F6qDH zp$`zgyYBd(NLW#rYIkiF*K0gy1oR31_24;4C)@o*lGH^NIl_3 t`wx(s-Hq-I*chz`QY%UiU#z g_(=bUMQ llPY_5>&1#4*Odqq(Je 4C`Q@-#tM$`EcJ5SC?VK+mzZ}x%FFO6m2adYbdcWp+ZSF*T}%XaboE7_ zk`bxIH@i}~{pC $!Tz<&~_pQ0PO zN?%E;P#?K@y%?|HQ0Pta5a9@~`)_VrCd-M0WPx?A_)U6T9EO?gPUo(AF0xD@&b-{> zJ<|hJ#OHIYRU#d{6`3RR>SmJpz5`7fzUr3`)XjKOiX%D^tc#jB3s^~(|54-YYPf=< zSBAQed!y5EJWn}b9GRaIqKxFd9*2wn67va_l_;)~;4S@+O zC9a@WpOR%F!POU3`J|^#e9J{@PK$39Jo$(rPnxm3L~apKv3@PE{1A`cNGPD(@z{M~ zxh;0WoL67A *94d#WKM)5lVoxnQ7m2SM)&MO=?nk;WH+fs4Ep8?i6y7x=;Dx7AKdtI( zZ%}V>S1+sOQ{b^S_bn6B1Fj4%2?Y7 }wBUktc%AlOF+n zsw)hhaa%0BkN9dY%h~t-fN9N`*lH*>-T>D2BpzrY@eK}M_zRdl>Ci)CdtTpD4r}sJ zmi5U`PEI9-C-`eEcMoBbrHp6>qSNG=ghgk2Xg0WXJNW4`rPst#_zz(26hI|}>a+UF zW5fBw68fqG5* a6QT}ir)04b bPb%lab8692{?Aw}-&yOE}7 zdbmQfLf%h{-ze!JA?vazugh +Lb^{gu9)l%c60%v=I^%+r+-oDvnlQ^E4n7xkn`u0~*cKKV+mP;Cr?M|wsI3<^^ za`0?9kB`67R0<$<$Sp!Ej~Z`xtX bAdsXKkY><Yh?m;mH`SX_E1-m+IDaxAZtq(# zP9GRwxJ5WhC?O1%t_P+} XI0=4v%X^X-7!PoSc9cw(c`K&v*qD7ObKVRrolXd~B_6~0V|55=1`OAwSd|?|| z17OzYNVRivQK#Px!W5%y4%vV_pR3(5lAgL6BIr9ObbVz5nH>oJl+)|6i$TiZ{VhD& zhGG3pk;{z_tGTYu2nopQvl!%vIGulu_karwnu9W$Yc01V;GXyqd?82Lj0aZaBJcm0 zP%0Hy$H-fdqIU)K^~y>0XMJ+tCP6KRFw%9vJ>>HE+!*_zo%QQJ;nGE~cvmmYnCc>) zGN*)r>-Bd{5~^)qs0{A`7`jq?*+a8cKKO|G@2Aoo;DCvS0S?z_^5(l{mACg}BU%_3 zvH8-EFc)(T$ex!r$3mvMqs)#5Ob9CQVP$TsSL$R-8GQn@`xxKks`m9QMw!J7%<;IR z&fOW9rQ)u7M8uz|kqwRQ^MGdW_xy*u%tdhk`A*$=uVMxbf>)1St<=MJRnGppf(6uL z=;koLDR0958#g+fdG8;L6YlRXl;wtEB~r!Yv&>5()}jvYl@ChFS=)Ys$;EOEEM9V> z&lwYek<;l17YxJbY|xd{_0PYg0$OM!AcxF-dr9Y?5lQF_EK6ohvh)!!nbbwW1H5+o z3;6_j_ Wq7T>($X?5QkJqIA0}U-Ubeqt5AN07q4DCQlAW z?*GWbMerr$*R$#YHeUtyU;MbT<{CjQpQH9<{?aLGHz8zT`#Is;A&*YH0YQR*93YK@ z`h`x$rXSG01q)l$u@{&Qym|;GH1cZMuQne4vO#b<8U}3f($nM1!Pg_mIzk Ztgz_54M$F!-b0R!p%4<9IafRb9piO=@LtK z6Xq4>vwjFmlQQb!cUg8W*dLU&oTeAoCG38AH-Hy*HOa5f76otpbH%zL+a3p27lM^- zOBeNVO`>F_OkhiGS>apgnA#Ki;3-zkzd|l9FlY-Ppa-@eemI0nape!^$u@B{s^D9` z`P}iuCp=)Ci3gDgJ@89DRo>yVJ-)DgO*~Bpd(M r?=IhQ NIE#Arh z7SRLa@9CBu&OR9S?29{i>csW=+7j~}5w3 f2ymilzn$cJ_QeO2`44W$$p%=>XhSjbw;!fb-o(1!>ob2_ z&`rGy7>Usf^W)2X`|4h-K@ZGAR@ToQXDdJz697Fd`1SHY!hTcU;n;-c!#4Y}4d-f= zf!nY10HR}J`6d&22Qfig@yFQ_=m%ol`HAcc4g)p$k<$Z}t$fK(?D!;|?{ib5kgxaN zs|d=hC6=wOB;HyO<*KU@6UwX0KXLt_zuB26JPQ=cV+|61P!IT+p({6d7Gu3olV 3ya} zAS}t(ANHR?wB}3(7sKw=LHJGiCjV%W%mrtQAkBcOY=OAFM#O _5y-BZ;)a9g_P>teld#3A*hRT WwbuH`FSA?e^3 z{eRkf)%Dpj|W&s#HM*DI%JHp$h`i?Ha0v-a!PUH|Z@XiiM&f6sfP$r1#z;AiQ)% zdW`}i9qFCDp5XhP{p(!k=ee%E`GfL2tIe7 -|18P7W!3&5NC%y5+8XD9 z%`o#>fKjG{$eh%3s8us+%(iIF8iihK{~`qtwlfzA8rVwkYCi0UEVK&>sZy8NSAN~2 z9Xxe4va6T^H$q3y2ohdNG3)Gr&JD_6j$HS Z`r`ah(c{@%$ z?kwY)en}FNEESBJoQXB0N+oClmi#ZN??Z8-+V$M=PE)EH-IPM0{r;?6?OEFX+H z2=c^-L?dEf=p5Jv4Ap)GwBcIxoMilPQ)->p;k7iXWFq`RANfElJ_}L+HLXj#Ep2vT z1^Kg7=n$q# XU2z9+l6 z?v3P1azF@P;AtOrGKdu`nPTCFT??WM8X{?bYC;x*SbLH^)`Rz*0c8s`ucR kkkR`zPSz_Z!XU$f;8YX>GQ6zfwFIc>Cn!p--Wa3C7D^g{5t~w(Ng15)LWu0kC z4M%{V&7{dgjB;{EyYqJ1BAhA2o}F$SkHK05J$SF3Q6mo7*3N%U5At|nx6gT>b>=RR z5D@{M2ZUIaMZC^bE4+d)YCU1eEkm>?Al<=S b2 zg1TY#mS6tyyu3o!aD$~w@39FK%4ME(xTY8LP%%-d&qG}70z?4dRy(lLGk2H+WV(wV z4(^?B_VY+X-6XyI0FTvR<&ze=ST}HRu|U*(6+-DgR>y!ZbuIM8rm7%@ Ij!+qq*~dWc8hzV>5dW0S-6!Ea!TZ}G)7*M2w0Chjl;Lo7ihiI=3`R5wOdV|` zu}c;jP}khVj$xNsbm-K5GMXcV!DghHOEfgt+=bmtu?W1lKAK1+)IoYXlh1&LM=f4? z#^llBw*k$R=_`|A@p#sWEfo!<-d}P*7$@lJ)(c8~}erH*~d6A$WFQ&QWpkt)9 zs1f{}@5LjX#MK5!Kc|Q86ZMGsO W5$1@zA*2@5{5`_Cb&NsS4CnwRv&Q^ zUH=K&PTF)K+7bRCp*CA82rz+9ISnoy!?u{$pP{)}=0<&0T!?7eDsB(r#C8*JQG#RU zbD6028G9m6U#2in*jZ}DZ+&Csadz*hKpxN*oR0+g)e~1K2QSHPDw>}px?y9AQBYKo zl|R)o&@ %2p)v?G4N~T5sJmutAmu z#zYp0jSsRM!h=(eNXx)sw`Nb|gYXJx>u}x~#ydE6Sr)h)j#yQ_Ymn!(WsjmLl}j~w z4&FG$NpDXTN;OVp#e8FU95@`X)ru_K@sPuzu6Yot3 vMD zB%CtaaLgMDW+n#2UOn p4cQto%k0Lj=QynnPRQzzk6ps dVZMi<5K%u8K-B~}$Qv6}HcU=6AJf?Jbz6Bz-u9WTf-`ezGQ=k5{bzX4(N zzTAXK ZV&l)yz?@c{&L@2fP^gD>sWFe>%oM@cPbX|GY6L(`EV2vO8$V+Vy zWE>=yDe|dCkRxMQhspa9VT9zHJMCaAr->R(TN3S`78D CcX>Z2m9Xh+JNT6vR*^j{^0L Q4iqnveH79(|zQL{>`@ZMT$EmYsjR7X4 zfa`Vot;e Jd@J^Z4(cPo(3r!sD_Qr;_)GMEdj+k+cOBoE zLk{uvLosKB`Gc|Wy_nksts7n=JdHuiEaS`?RMin={IF*|!up^wfC<_^=qb7o0fnpL zD}^Rq(vb1@mZrLTI#$hZV)cwOM_(z1{Y1`f8VGXs#?pIf>lwmCnd@nGDaY!)NI{Z! zDbktNJle74pzS;10pTqPJYj;^SF*ku(DGbg@Vb75Zu?P;8_8EBmB4(-J-=5zZ&=25 zpCn-y^noQI`$d6bqHmw^otmr?2r{%i`$70aSZmti;~}3Rk0+rE+-tzGF{DtkGYMVX zzpKDB$e{O{^lsO1@CN)B%S2=i)jq$NEXkkrveSF@IX9CtecQPhLGmNxWz^|pDleCB zpCe=B6XF0ef*_YJLnZ~6^2{bZ>79OqAo%Znt2GRfL&u2!8jz6kGl&a#WQ^7M%4kG6 zO1lv2A`CX`PcA(dJmHoeC+tS`CsXnnSpG38O+5sWcz$K@TVM?`&qz_IQ`Iqsz6?Kr z$CB@#*ghX$1auw|)d*RVOm9H6!dQd+Q}W3UX;jG_FAD9FIIr_0K*lj~o+t =Llv*9ue5d6XjP2IQ0tqVN@Hhji8aPeuSU+|Dc|VF3Q-Tuh z!B}+v3Iwp?hRzUBgp-7~zX#odrRcw7Z 0%ciB O0ONy~0fLFSh!BS6%RX?)wDT6T~^7bBd zW@#ot5L~Mc5DA2XrY-e$?0k@hxsxy0)YFk B)c?LTL#bx>vnmacr#ld$wY?<= zA0S_5nHbrJZ_5-|QBLM+(74PXn| zff4Ufk$wMNlXq=k5l$ooglb ubB;y3PhS~K^hCo6Q2VV3dl6bB`bn2#=y;7cj;>wK-Zm-y8 zS~1^Zro<6_d~q#$+spc!^(F_{j@OP7M1O(3t$lhkO#phX|J q{Y46TkN2q~eRglS_&)`aI1reH2;Dx2>F{@|Q;@KErt*2F)vgIm!I(m{#Om zJ~^>BHQ>o&1!p^8fMlJ?0?|Hid#(=LvJie#cka&bOkWu>R 1 z*cCQ%1)uW3?;`7x820Ha=(hR)&?hk>8zz4E5ZC~!?0$WCVufmD(cp{_8WYJ0U$|c3 zc~*6sqHByGpJgW~`Sk1Xj}vAm%zr_E4q=%W)@EI_?kSj6I=k G~-zL z!=cn=EP2ip+VDtpHudsPe6!6hL%vOGU77xbbIrhQ^bN&-bt{jKX5C!^00@n-?TP!G zom7h4Y+tM$(O%su{RHs~n1w-E`Q{^7ci#_ELR1ovZ^ORl3zhQP3Z*m_l=P+nF@)Nj z<%m$Cyjt5F!3znW_UuakI3Ubag07zqmk_I}t9Ah=f0~LsyZ7+|lJJf|aRWl#p!d%Z zpOv3 k8ET4Zuy{u}m+I5{0ep=BIrl;Pd@|W+1J5 z{Y$x_b`6C4m`BqhMVoptc4dn{ywb`R_wWt+xAaP;>weEWU}u8f2lrL+Acg8S&zVPe z!^}<<@D*HK(KeeDu~>o&eRBT+7?@Wc4V-V2Gl^A+G)E$VjCyaFPbtwEQIDQUBqf_G zzqv;%YqSWCaSSC9Fj;$h`czzc%MF^p=uOa=aXhvaf3& T)4yb*DmBz9Koh97-mOTBU0f^T(8$Ll%cX| RM;3pwVPXI~0oP?cgMsio zA`>M80Fo9{Bdsv|p}PFuCY;DcOd*8_2o%ubQ@yTK%dmzPt^s}_Mji4)%6-;t>VlFN zP;eIa{VtXjxRr0;3Uz5eFChuO$4(9X{4&+mh=NOc8VB>yfF87v$a*m3LIuwE84pB> zhMuxt{U2+h;2;iC)iy!vm0w?ftaHt6mV|YuVgQ?>i^E5g0Ybg;NdW7PG=U#VGk}F1 zv9F|`?>s9ULazHn#)HmT&fFa5Va9FdX?D~k!^UC~s-{*5((=6(I>#*6oABDZ+fTlF z2S?L`=D1Yx&<|@lqjEugJ~Dx)Z*0fcG ewmmIm{`fc662^n~XrT>*2wS^eg-mZHK zPSF?$;8f_q!}&{$^(17asIT%=z`O*H;Rh _eq@_{0P~nlfaw$ka{5(so zdcMX@2_d^F1Q_fSdUz5wfy3L`Y~)re6e@iToqvJcTvpA5iWcaF32Z^tG|4d8Xy&DH zs(rZ+ryJAQkF7=o=sz%WaH2SLzV7W{y;Pb35PKa1ycyn2<_W7NBgMx^p(GE{4{Q+X zKx8s16BL5k9`cY|?74K#5%zz_g}9tvu=;Vq4%b`LPZgm0xV^W~d3YYAa(wsc^$ADkZk68uGT*Y%n}lrjjGA5(5p@{uQphYM5ZT@*^# zw6^~YQ9|g-Z+GM`X6VFoJOk;V%)&4VKA0%$5vL*(aPcNDAod8|#sd!g`CO<}BD2x~ zmH==j0FbrxfboKlACp69+s8JPH5YdkZv*=f%ktr%fr8lHM3+9!1W+W>ae{Pivw}*< z$G@P$fW=vLVs>W!P`}iYVxdCAlOIPf|NeM)dY_Pjn_W9?GkIt70-Yxd4#XQoKQoYhR27sR}<^zLc^H9&>SwAo|OnZk)b zg5nU|3v$zo6xEX-E$D&-u_(gY-?nq46c9N(ogdxRb9%zEQgm0@Drmg;kF^Y4D+crx ztY0NIcTd~Q?8$;?@L|;enZtjrsUZ?1PD8LjTI}Q?Qnrk%0!%?M^uQfhB19%Lgj?Fm z-w}ccc?1K3AMtF%!thKRmmP)k6=(5+2TvCK0bwpBRt&%=By-F~-=|`N3+Xv@?^4OH zvViPC&Cr+!k{7nr?nCGf+F}`~Q+?d|HF=)Y` WN8Td z-ZHOQ2W{EL8{mVXj3}nA1g0do7Uv(N5CnoXa-LV=-IO8}UDV_EvagZm8vwyjyyys7 zyteaLXVXz$H;fNnPjx0$i*zJVAtV9Ba0}Ep!kneWP`hs4wx^+%IM;TtXQ_VrX$|w^ z6doGIq;e+AVJV1kK<4_9`tyMy3y8_U#l3}MtDRNnL=x$OI2{Z9e-=e{5zDI`dKvCj z_;g+*|L0|r04XiOUSl>CX^r3=l)vC~z)4be_b<(`3xZg)8_pxV_9uf}Bg#7rS|=Ja zlmr_NA~Ek6G3KigD>eP!%V66Eq82$ZJ>swbs|A4Xn-a&GB=cbd kb5`^u zE35-{An$*6Aob1z@_^eDWcDo{k7MRN$phd*Y6)v%&gaF>%g)dRM>C9tl{0T&7W+h{ z6c3HxM+)iI{a=FxH(w{bvrITrAzr^B7xIe0|8cS%n;q0jNHNaaQ7Awb(5>j33E!Z5 zX^ovAa}&sc2t}Y@7Xb9?y_A%8RO~5N=KboR5>V$6z15%u3Y=eB$YDdrH_GW^L&vad z7=_9Cj{aKmZL!Xeq~&))K|$)O1LlYOVz8MP`ZX}1;6a?8N gjeFW*IpBjWu|f zROBYm=uGr7qx$2v5<;-|32 !Fy8G2n2r)_Bg>Ux@t=80oS4ICugJ{ zSNA(I@L}?1048zM0|o`@$! (Uo{-3seK-Bhu_FK^2SY&zXTOY% Rg6a=O9v75rJ$p39*I z>x`)~O+X|#Tqtf%vdEyFrtWM}a(;3_%fGMFL9d52eZD?WP$q|Lcqq?3VDxnwbb`(@ z0{IEz9W$9fYqLBw!tT2UzcD14oC8g6xLWB2lAC;7%e!YJjdIBa^NDg|Hh<7wV*z zY9=V{BaI=t2G`pQ2pAes<`CLvZ54Tx6*4Z1@4rD+0dw_z*k!mAi<0Ld-L{iiCyC(2 zApkXW-thwC>1Mq5B8*tqw2^UA=Lz6=QpX+0j3|T<0bPD)u2Qk%NK82{6&sKd762 zG m4@kKRxHM*~`yyV?kLfGorOS)}UtxTnEa0(T2~wtKnDfN+*D^2*MX`jAc!xj|Ax zw2d2;(}{%7hQ7r?H&YJE1e)Jd+tcme?K_Rec+^ph-is|9iTp+=4GyS`@mM=p#J(bY zCJ=rPS`ZLKA2YPU7m|}gWJ*|rs{-vDh&O$bDo5PyE;?3zJ+VMNn!J4|e&e4#)IvUl z!qK;<01Sm_JS`@}le_bbf|GGUga|BkP-GB8%H$fYzfU)%Rk~gZq%>&UE|K0Ob7z54 z=$|7A&&ZM0rb89Ly!^WxEebu(I{DW^q#6t-6cDDqp2!y^i+~@&_rs49|2V{L$bOfu zOPF~@StQV?YJ32w3G~aM?vx#?g3}5srCB_iB`Ns)GHE>9dtF4~XmoJYlw5*QM#~t- z)r2SG%zARW-;Xv^o44MdYWAY_5`OOT_vIJ%UtGz}wRaddPIw(juSDTnj4n?9K}>hY zp+a4FKl~s$BQ)focPFa}$-m`&D=-s$9X}Hj8}x9i4kl3pN(>1xcX**0AtVLh;}3Tj zht|~&`NfcFg^U*|YT7Xz8X)CN>sR!8-Wp~~ekNTLO8tULlU8iGURP@s3_aH -<%ED0gvdYPce zSKAgzyhJ=@Hu-yynxI9rCr @=+dJg2qZiCQcAqvY&5c?;p zV0)0e4VD|D2ejR~FeW&Zp-59O8emMIHoCigcW*8D%NZru;Rz$A(Kp@R zBmrf@{#Np=%g5((dr@8G9(>h28?yj3L=U7z;StPA(!ylvCzIMqBFe%`cxhkVgXMqo z{WS7@lh64)?|Bh%xJf4a9rzJ#jq<+yI^&9keHsZ)DuY6&6Dbr%cb3VP{D`@~yT`k# zUp(yoB}t50bzmhL-Pj8bpGiDR_!F!nPMS;_2yC9a@Fm_H{;xJCA wk5Ny$PhgVu0Lnus)H0RZN`%hOX|KjC>Y66n{nQ}^7G59-E7KjlNgtN zmBMBp02Eg(#GJ1r?gfdGi;$I(N^&94ECyO6D03wwm`>|iBb#MY+6n>ZC1~m7-Vm?m zbW$m!gDQX-%A14W$On2Ep_^oy-Bt*wA)k+b-<`1_q!AweoWW&aFVYr+>O=5$7C*^u ze)xE0SI;A~FS c` Iu)z7J-E)&u>1%x66SZtb%RKCHUh@Od1*DVr;Gv;{yK)F%nYgL9}i-WdaC9|2F zUmrF88;PQ`f+59B3F0AtgAyHdfOk>oWCylpf@rBtliY)MU#U1@+xv%!nbSYW>Yig& zehu)f{kQw3%j6m4N*C7q4z7rWiwQb$->4tIP2eQFA ){O&HJpNqmLUMKJpD9nzt z^MLF=&dvh`3I{cx7PKoyC1^|!yP5ol#fG?Vg))R40wuxlZ!my!g)(!i=3M4wOXk5t z`><~ENNhH3@!qiokOiS{>Ou*u8kDh%BEl!`(0CkO`+b` $9N zIRmRFPlt& =1U;@@q%C%5;tLq} z$AdVjJoq`=@xX0RX?N+j6ALQ8^Y> {wK4zQ1EzTB8%Wfq%h0 @^IM1nd4vY!pTx9fN;R<(SbLANH{yP_$t)Mszv!zU61G(zDSP=JtAr!tQT|b7$ z=p!gBa+f-_%^^m{&c V??G(mh?Av#igDG!3#$0%=Q^)`2vF{3u+K)PfR>cU;s!C zAr7l@VGv64LjKwmizc{VGdzw+L68Ja21de`a+8muy82d2G+Y`aJVOfWNqi|<00yMe z`RM-%^N1`&E;FnD8f!N}jrDl0LY1&*x#MXQ7>rJIkZSZG`!fe7W~+rHMuI4J4=()e z4>hd=yCYe&C8>}ED5ETJ%aXwH2LV-W`8|js3YdwTlSxAM2$%W-u-lmra>)JQS5I!= zdtxu*% $DAbx!Ph{pL66vuzGcta!Mo&oE&{Hur pE(HSSSNCej~2y1 zVMWP*J^Jr)I(R%i)E_)l6+-^(qlE^`u$U2@P#39+0ZiQpm`=R3aBy7}D~-Pg8HmjQ zl?MMdaIA??n2#~Zq2QuTAosf>yGiDZ0tEn=|7JcN4Q8H`3s0c%gR`IJ_I=Jgr%w6n ze`dA>bfvA0A7~0Ag#jp_xKLtVL8}+{*TkZqgKTAH5*2I@&Pq!r8hn1s7MzjC|Mh(u zgwClx?37$8pyt-`$$~MdE@ m{Krk5|Uf3|3}91rK39qdvkXTgD!HC$@0U{=d;Bg+TXx2HaQA zZJ#N)^JKx4;7xe`?HC^6Y(ZTSCanz7(XS9H=F3Jr3&sFO$Mf%gLA5Jw-^l&meKp;h zNU>NXEYo~m;z`N>` uOyPW@z~+fPne$mOO0Ua*$q!sy99Z^v4nqP>DD|DY zEhcct>Wp{(-7W+qGG967VKJ^-`|_6)fzhEdL;p8v6v$~}LZC`lvn|bUg_ZyTn{B=O zPhuD*pb4bhB?V-2sEz0o=xXBBNWld-dvD0{duyb~p{!-jDegULAXrl}p8u>NZWoXT zwjxNP1av^|n8HWzK5lg8ERc`yIku)+fTB7^OgET%EeQZIeN8eqwiSAYeHKJZ@c33J zMHFZV9n1L7`=GVOh`IGF1qub=PK*L=(HTX$QP?P%{0j*|?Ychf9|A35@b@4W@#*9E z#Qzkg=g4k(JmpdLoLo>~_42cJl{^lU3oN{36e-QO#RSB4Cs`-%sNUkHM!$S@{r2Ro z$d}=9|A<79`872W5l!sh*gxNTDtpN`;XM6QWyA$4#8ZkXWDH%5@9d6+u12Crm5FhL zruM;3j=N5B)?wfHp4Rf1_C~h<=6l(0bytVBj?sYfd0dTHr{-Kn9}?BF_P&Z!e3B7e zd2cxS%HG# D)WvsPb6ZdQf3nfwF<2aVabU<72iDysOjVp zpdJ7AUS+ieafDCv|HTt2FZ2Cq8Z?B786v#ibMRZ+G%W5S@>N_P;}YE_nnsK!4{zdX z_$)UaH~uVhs$|~ tm}EpToMpoAv4#}hgVzGY&PE<-CJ!&Fnv01RQr}_DUg7-T3vhwkKU+L)?a2X z$f>%w!npJe-KAP2*?G7Uh7eUR$~#`)`pu^?nJc0;d_T@zG1ITsZL4}jD4M(Qxmw%~ zC$H2>X|ox$h*o?*^{22-1X`3o7}?BWj6sj`oXOD9AFSorV7igjzyXdceyR_>=``x; zi;F &I>}-tYIbJjx?$GhaQh2i~Uwk8r~6pKr}X zE5=P_{MST$LLr#UKdmV6H`Y=ir|77Uto`9?2U6xjYX;6JYN`o8wM2<(m*9$>D!I1= zu3Lc-25T;;40IaTH@NmAD+@6FP7G=8lhzF5dj>I6u?V)bovpFBseFvcC+PFPVsCOj zw`&c)eWViis@jQ9rlMP%vBN(yP>!K6Y!dN#WF3oeQeI8Da#SXb^u9|~f1U>#A~ll8 zO>VamK5mRodN-vKHWHz{I)*q(nzBIGy>tCjY8W{hxhg#r`zRj~!td>ZOSVI;4WXB~ z+#Aone->H8*V-PrFj%o{8oA$%541P5wFRF8hQYozKR}nI$!NMEM>Y|nT%L^~_N86K zM3q%MBDuY%ezp>w)W59&jc@yqN>p>^5Ir6}wX_~H2-Y)&P?9+0V;ZwnSy_smsz6Mr zuPs+1TSFqbFc#HN=W~3Y 5=*~y8mk2L&XP`0?t*_8hl1qm7R<}^E#_}Z{Nta@RyPq1`7I9O9h K7s)6XSxI5dpK6yz7=kaW6L-Z*Q;=`)#0x~^nnk_`)m#3 zE!j=p775+NM4Mhkw#uN(`g$Rqgz2DYO@7-OE8o~^LPXe%mJ7E|q#`^-^D;@N2 wIJf(}$dA97>+&Zs?~xIl?2^)$gqmoBgHvcbsudh@A&&q-Dd> zd;1Y|Q4U63k25{lM`l_E!FONA;s|32R)+caR@)GAeS^^WSpBsDRAo#||JYj_D=N|2 ztBS3A#ZNX@lxsy=^Z8{wX)GDqjd1+wLX9gO&5MYg8*MqWYq(KEjqlLp|52bFpB2xs zWZ5khZbjY{m~Xi$V9hRt02(yF=Y4Ckt~xT|{ZubvBJ+4TUs=KXu0CT7u_t}gg(_F% z(3_7(J#V)4& #qYC5WR5O>qS zj0D5>&fpRvaf7wzCu_~Pp&8MKQHU1X&LBiJR^t%VQcEKVQ=y0q 6zC>QayZX!JSDg~{g^P1|=0n)91&{9%Md**GsUwF{2cK@{iN2k6L5rXu-gELiI zeaTojp}sncnBl9u)LMJdX>46q-qL?8MU!P3y&8xCQ{6+1{5>CH7n2gFhVkz!4;%-O z-sZ|rz2XvauJZhhA96}jv}1wlYu>}rQ *+4Z zB_+i41@nPY`9psTn~lfxIHJP&X>uDSwYAD^6Pzw`$Wo61x+l*}2y}2_?~#q`rm#)k z$&u{;&06oOt#l!W@KX{fJ7&>yM@Q taiX4w9KH mD@kY1h|tH1D; zLT4g2lp<>+TazpvLq-cW$qOui-Mk|~U}dH4FoIfWzFBIz>*>wN8Y# +znk0i^_7Hqs;dnouE*)vPIOa=zDZ6TXK?or~6DcojDYdLAuyu4(`QerdYT_r< zJFML~DC^svvgOfNvUVRm%cE>~w03z5vQoZv1_N$7OV4@WcOPfXFQbGbb^rC<*_sqj zb=-{9X7|0-c2u#?UTNN+ZQIK43(v!nM{r<2`Vkto_K^?E@$P};cZ<%xaBv6}W`ACa zvZiBjZDw{hwc}7K9~eFntFRUZz7L2Y=&VHYBav1~eSJpTgw4!qvSPEmC8`{9P7F)z zoB1}AxT!Dbyx%FI@JK(QR`altEyPf)_`LftSMpS&yEJ6nu& #)>J z7NRT_hOz=R(8B>2r+?a&rz0E((h*Ku%C26@;@dL@=l0eSGhum^7t2j?iF1y(nABBQ z qefYPNGV&_DPq!g1? zihOb98Alr$OX$@u{Dzs5 _H=>Glg4Op|oI1yEGw 2VbhW5y< H#n^rEB?GE%M9xKaF6ID(DOTLtHDUE#=H>aw|9iSR{w z8{nj+=k`r-BRG3?S8rUs>T1HBYO^@C9=jGsbe>94Y0;v>FDZgppG@5^oV-YDb);YM zgDv5(-HF7hS11d{^atD)c8Wk*Gx#d6wl1Q}5S8Ms`e>85M*}DO7*agy5fS%Rad^c5 z{pdQL5$mQ6jh?*+Sg!Z07wW&7`|7rFF6n1+-f tp7H#H~XE{c!|q;PNi)OWdzQO(sMMwSg{ui@J8cK*c!Uz4h>Me2xxL$t5b zi5yoAwPo555)z|<>hApvDLQwLmiRwcR~&q1o_p>RS27TOS9@;6U?_%Y$pID>KAYl5 zYm% mhX2io;kib7wN zC(n&BF~)B7)hNXBy&bM6W;pKR5{-xppXF?p?~<+lx^ncxPjGQSVy{to)E#wHl3odR z7ri`Do%$Yca8EymE&X$9SZ=xGT!nmV9I*Lga~!7gn*D>uLTf+unD9fGmcc61w4CR7 zN#zm!n5it~g*Q~6*!Vt6OAJ#ir`GaN>Qt8MSDC10ob_8nAEvPE#iH`Sh8gT`MG9<; zi#Pk2r|~)1bL*jw2+jU1HLBUNEd#VsMP-N~2g0fuu2-^bGaK*b!1O0Sr>Wi?>toIm zDo&dF*|#}z#l#uSyzN+olInZ^-mY4duCU2^$AA4JqYlC>a@cfian1d42y5xhqr13e zhDty7@~Y+mBV43F`O&3T2RNQOGCi+QoBQt8CUgD$p?XI!M!y=QWr$X%)E%{+R)r=W zqgTu88umRvS>mG)YjxE&8kE&~27UZfM(r@_4#tnDM89X{{n_}cEJKGrclDaS*sD3e z#i|NSf1_+uuA1v@q(=nV{PU_XP <*}scAX=W`!VJ|HglzMa$$PxbzSgfV9XedH zzRa(^sdPDA)|@#~yzCrqN+Ik3fD`L`(>D@<;3IMt0+m`dbo3Zivt&)H1B6>QmB;S# zR|ckaFrXP#D)04O58De6Z4G3+&)-pn_Eu_p3$~1o@za9o#LfNa=n*+}n=vul?_cYs z5i+HP^$M+fE%*89ZCRT8{csx6n;IIfv&?P5==%fPSzQAA_LDPJ2vM}F2~G|&0CcEV z8U8wYOucDc19{mT8>X?YKN^a0qY19V+dtz5*pS}#xRJ=3^R2$rpO|!>p}mw&6x`at zCqAoc5~-?OjUCG bQsHE9K+=yuKqd zptWS^=xf@BnigZ=`%!Awf)_=Ogk>(L<_MlBeu!7OqS_E^^%9b#mQw7dkI!sk)WP*} zy#Yv&sawmQ=?OYX2y4Eb!I$&T%rkxMW!Z7yMjY1b=v&vd7T-@bh8?gWD(#Mytj@!6 zs9wnOEM8q*b?0M`O`htrDoDP_%JpjQhm2~KdcV)2J=hBrW8K{&Y+eAUq+d{4RG;V2 zi|Grr8c8QDY-wUyb#;PsFEDv(P^i{cT(wBHrQZlwujb9 ^xqBBizk0ma}LSWw19v&sBxGJZGYoymyhS*?c z27Ml+gpOz`GM8JtqDt2`o}8bGc2!L5inLMBZ8G65_ImH<9+q1xXlbuj Z&B745tsg-a`?+0cJHWC&hUr7D9m9B~eecu&LyDBmfv@5;zIs5)F*zzt zSNny{gMq=}SaCaPl*LT=>1cRl4mc6dqDS(TD`8r4PB{z9F^|em_tQ;c1_ou*i;nsc z10FBOlBO!!7Rz%A)b$wn7K_ipabQ3F%8Mr2r8W9Cz8(rcOYAqc209zrTMvMV7VfCO z!ruzVma!IM9jgUq-iXU_S_qV~
Z`<;h48>u!)Fa96J z;1kj{3Y~#cF{uaI=O^yJs^)F<{+wS7hLPRdIw)FiT&=!X*}f3FJ`kQQUd41u%od!o z8J#vc!9>&IapPCDMK`%rak9<*Yl4 m0B *e&|C(TosMWU zAa7!jm4R`lAMNk00^X#o_n=8Jf{CwFnYE_2&au@J#+oT>_i<#?b5=h4r*_mJQ&^>A z?3Ko)iQ3Tox$;&zl#{%yF>_Yc;_wmg4dFRC%zS!nz3hIqjA~ zA6p)G5y!w{+_PtK9hN)Yl4FWY+~3k`0!+%M#I1sD`Z{pMnr9B{GE_!q(c$Vl)RJ_N zV%#`d^0@xVgDvDKT8h?!&I&qNS$O@bfpzA+yC{j1=5f83%M_DwQ?+R=+e)trzpT&B zd@O%6B8JNykLaSsuzyF7j2w5XcQ*y6u $mTpy< zx?FU)-mfiA@Fg(Y2Nzoa*JAd%V_*z)dRXrs<^V{@+0PO +l`)y3>(#) zCZo8PtK}WbPHL{MlUsqD$rCWFrPnj?HDLAqh@Cghht;^z2Ns5A{fYhC-p;s`IeQwI z*n21@`8}|PTj(1+x|YZgq)q&V2>d=680hzc6nDecAXxo&Dj{y9;mJr`PlQN*HLs+P zPMcf7+bLS^pXW?o`R=kGlUboYH{xEUo90ytt4h<3kF(aGS8MN#-cE5m1rsw)7>C0O z#OF0xv z)dE)|jm#`h`?Hrgg8xH6 YHw({Ek-zI6F6V#f@Z%^9(yj*~`KPKa+PUTQ_7MbP?k zJHv>;=%PA1>3N*Iz_ioqAH;;3t6^JVApi8T>!TNYXH&*v2jx0vJ|;V8#_Uqfy5;us zd+RGFS17a!hq*?m6@7?sJ(D=T?X8D9&{^? gNX_U-CuXy*sNE&`d*fwM#Ey>gh01sntCg}hL)9&1 zC64jKy%M2{x}|}*c{#kBn=-B}`E1Re?nBdoo20<@_djv=a!O;Yxo<`ax!V@_Kg;FF z?iFTYetk23XgK2I7_zx&;5|BwFr8%{sTzxZv?a0;c)69|M`#&k9BvnZe;-?0w%i@E zy V12 zwKZyFjHPBRYG$(9_#s4L{r(tdZ#Y_pkl~x|-LAaTA Z?r%ImvUyysZ)#T*#ZtZ(0`@ z(I$zwgai#h^L3YVec#~@Kk7x~Ys{K%X}rI#IB?0%0+)g_(nihS oUs?$g1k8pgI}`J)gJwv7$=E@(<&=dhtcphx0}2 zPj)?k{}Q-7g2AP r zM8?RuzvKpA8eWc1J-pF|yU>p83>liirLpnWAT4DgrgPL>uSXyH`=D2nv|a^jG5vw& zTJjC(N1n<)T%d3S!vby@S6g9B=`!imZW|uG>1St{Jltq89^KVitmCWmbK+atp<`=- zuFga0(pNGmENj)qZHDnyt}c~RbIM7{{Vn@9LIx}GbiTsPhb}Qi6Z)_G^_U!Hr~X70 zg)cBmw^1xJqU^{-QxTYfl2@Eo$tubhy#H) ~ zUZ#)kA$L*7^4{^h7BQ$9IfW+OktzktXnW~5wBO=)9}k3uwfLpoJSa~)z%joG$YhyV zYItHOa5>S+N5n@BStrR`gPX~iu-JWwJp4zbP&B8>aqfd~qCAvN@p#Fp{j^;6`(_&+ zd1m-5m{W9KZFOBp5x}imIAU%W^U}sgEQ}Cpr=U#_!yX>oZ9^@MB^VA(IJM2?B$XXa zD)RVx<9w@{7T&s&6{jL-q@~M*!xop?%1VoblgKA#f0`XhBYFGL8