From b697c842dd79b7069c3c71dd7a29d24a3bc62d65 Mon Sep 17 00:00:00 2001
From: bcjaeger You can also compute variable importance using
@@ -237,12 +235,14 @@ A faster alternative to permutation and negation importance is
@@ -251,12 +251,12 @@ The out-of-bag estimate of 1 (the default method to evaluate
-out-of-bag predictions) is 0.796875.Oblique RFs for
#> N trees: 5
#> N predictors total: 17
#> N predictors per node: 5
-#> Average leaves per tree: 22.4
+#> Average leaves per tree: 21.4
#> Min observations in leaf: 5
#> Min events in leaf: 1
-#> OOB stat value: 0.76
+#> OOB stat value: 0.79
#> OOB stat type: Harrell's C-index
#> Variable importance: anova
#>
@@ -159,9 +159,9 @@
Oblique RFs for
#> N trees: 5
#> N predictors total: 7
#> N predictors per node: 3
-#> Average leaves per tree: 5.6
+#> Average leaves per tree: 5.2
#> Min observations in leaf: 5
-#> OOB stat value: 0.98
+#> OOB stat value: 0.99
#> OOB stat type: AUC-ROC
#> Variable importance: anova
#>
@@ -181,9 +181,9 @@
Oblique RFs for
#> N trees: 5
#> N predictors total: 10
#> N predictors per node: 4
-#> Average leaves per tree: 5.2
+#> Average leaves per tree: 5
#> Min observations in leaf: 5
-#> OOB stat value: 0.64
+#> OOB stat value: 0.73
#> OOB stat type: RSQ
#> Variable importance: anova
#>
@@ -220,14 +220,12 @@
Variable importance
+#> bili stage ascites copper edema sex
+#> 0.160185833 0.042419910 0.037004321 0.029450477 0.025523886 0.021037248
+#> ast age hepato albumin trt chol
+#> 0.020497704 0.020317763 0.020162922 0.015865138 0.007418603 0.005756730
+#> alk.phos protime spiders trig platelet
+#> 0.003328945 -0.002559019 -0.006195550 -0.007907720 -0.009315955
orsf_vi_negate(pbc_fit)
-#> bili stage copper sex ast
-#> 0.1796153443 0.0654718662 0.0377992169 0.0300985949 0.0184556866
-#> albumin age edema alk.phos chol
-#> 0.0177062224 0.0172170446 0.0130478892 0.0083611824 0.0069716249
-#> ascites spiders trt platelet trig
-#> 0.0063832024 0.0063499175 0.0015016478 0.0007794253 0.0007182015
-#> protime hepato
-#> -0.0006160457 -0.0065521643
Variable importance
+#> bili ascites copper age ast
+#> 0.0770176686 0.0286127351 0.0222179569 0.0158979928 0.0105580692
+#> alk.phos hepato trt edema chol
+#> 0.0074408002 0.0054291589 0.0052311839 0.0050824237 0.0038163615
+#> sex stage spiders trig protime
+#> 0.0027274412 0.0002958220 -0.0001345556 -0.0034778432 -0.0042680518
+#> albumin platelet
+#> -0.0130420884 -0.0152792851
orsf_vi_permute(pbc_fit)
-#> bili copper edema ascites albumin chol
-#> 0.076377458 0.027301164 0.017138302 0.012935499 0.012645786 0.010981566
-#> stage platelet age alk.phos sex spiders
-#> 0.008385515 0.006916673 0.006300806 0.006115436 0.005605453 0.002729921
-#> ast protime hepato trig trt
-#> 0.001706344 0.001490477 -0.002072969 -0.003145736 -0.003454973
Variable importance
+#> ascites bili stage edema albumin sex platelet
+#> 0.54545455 0.40909091 0.31818182 0.31746032 0.31250000 0.26666667 0.24137931
+#> age copper spiders trt trig alk.phos protime
+#> 0.22727273 0.20000000 0.20000000 0.15789474 0.15000000 0.13043478 0.12500000
+#> chol hepato ast
+#> 0.12500000 0.10714286 0.04166667
orsf_vi_anova(pbc_fit)
-#> edema ascites bili copper stage chol alk.phos
-#> 0.42139078 0.37500000 0.31818182 0.28000000 0.26923077 0.26315789 0.25000000
-#> age albumin protime platelet sex ast trig
-#> 0.25000000 0.22222222 0.21428571 0.19047619 0.16666667 0.10714286 0.09090909
-#> hepato spiders trt
-#> 0.04761905 0.00000000 0.00000000
Don’t specify a
control
# control_fast() is much faster
time_net['elapsed'] / time_fast['elapsed']
-#> elapsed
-#> 52.71429
+#> elapsed
+#> 56.2
Use
n_thread
@@ -155,7 +155,7 @@ Use
n_thread
#> N trees: 5
#> N predictors total: 17
#> N predictors per node: 5
-#> Average leaves per tree: 20.8
+#> Average leaves per tree: 20
#> Min observations in leaf: 5
#> Min events in leaf: 1
#> OOB stat value: 0.75
@@ -206,7 +206,7 @@ Do less
#> N trees: 5
#> N predictors total: 18
#> N predictors per node: 5
-#> Average leaves per tree: 6.4
+#> Average leaves per tree: 5.6
#> Min observations in leaf: 5
#> Min events in leaf: 10
#> OOB stat value: none
diff --git a/articles/oobag.html b/articles/oobag.html
index 29271892..7f8225ad 100644
--- a/articles/oobag.html
+++ b/articles/oobag.html
@@ -139,10 +139,10 @@
Out-of-bag predictions and error
# what is the output from this function?
fit$eval_oobag$stat_values
-#> [,1]
-#> [1,] 0.796875
c2hyZpi$0GuS=cpRjX(Y zTq_-op?ZJL?`HefIaVI3IZgDpI92b;hQ>PvP)RX5&>y7oW#L#0l-!7U0qY6* z)e4qb7K1`G4&(jmY*14QC4s#S^vX+(n`(cb)mUuu>8pd*!#hc{EZ6%Mk_m93)L0V+ zadgbyJ+>nA))mx?N+_f4sc)kh614j4?2u33m}-e(-ZgD+IcDWkVnM$3c<2q0;zaAQ zT{*>$)7T16&riLzNw`7#QD6?zw&pE%D1HBQeU+J96@Orct+m1NhBvZD%DXuv`*w(3 zWs?8%e*g8uBg80fZQWHAhiuMc9vW%VBpg)TPELP8)A8EciX+l=zHxxMQ)6WotKlhe z+=Yyh1^v|b-43G*ad3PF(-%Qne++&QO$IBAhh@bSo+lJW@TB- KQcvd>8WQBiP=R_K{$U*18I2cLdWI+CG>em`RtvKc1y0X#izO~! z>@P*!{tyb0LN)SU_(LQ^JW#|#j^1V;nt@V*s<&fK*V1dRb(IKlROyS(H+cc(In 9%n(wv$TNQw~% z(~!{32jf0}H&YxNFFU!3sT4PZ=F^j?Mqwm$sphW2(v2BzpS;)vgiwF*VR72^tpaQG zP9a%1|ArcGQ(jD&Om`+8mKVD%YQkV0UyslS=`4o;F0CTF$*pU`$*$F9Zzt>Ky42=J zX}rz9R3=&iCd4K``=~J*+5e#P*@i4>gk^C;7c}<%-dj)1qo2gB`*q0+eJwWKk0~Fn zj@U$~ZOq2whtjj*7Pr&A%m tf@^#?r!@t`n&Hx=bR*wl(u`4XMf3eu~q)7j&wj z650-P_Tr8Eji@EA4}-_&KiSOsILO&rq_NS_IL#h*1F-LMtxPJSuRE#o{JRjiSSPN~ zz?G|{o$}N1_3{*3(EYKm23Hvfzg>%;pKUpPT6Yjrd3;M~;9deShNYTI?mLY=>DEC{ zvp?Gl7m}rDW4~rivxjEev)^ZfOCir1!)&-0Z(GN78A_}5rH)JoI)FA1o|TJfXBR9m z74)3@3Lq!AM84h?C!VcZ#j^KZc9PU $@e8bz671it6 z7SoLj(e;P`#a}@i p*sEJLrE#P3|52|hgb)>~YU*Y;r;llq|=dQw` zL-3jIw_|!ZnRPtBf>f%Z2-4|ag`6)BJ|bMm)BDvR!nDaC*IU1+;q2sWvwiV3T)fD5 zt%y9t;8EnqaTLSp9ljyEh~DFP%5Zhv2iiUvVsHZK9%y;M0S52A{8SA+ahVy22{eUw z@*YKBBy_4FRvrbLAfrqJef>E~7$8f}0qF7M3gyMLjAv!+=nOQ)GIZw8b$8St>vOpd z%-6|Pu-EDDzL5Og9Vpn=_W%8p#T*~Fj+9%4x|mv#in6Z>7R##vM|7{A`C8*%l3S{j z?U>R~FEayY-FjCF?ZdXZ(`PtmI=f%ot4F%!u Q$rSiA(shDWF8RiRT^QtA*7h2lE z(tTdvli$rimva}CZZ9fswU1?ndOjcz^JeY0EKmu>8{5$z_OHiYF5dNmd3Bc4n8wxY zyY^*9UhHdz)${8%V`J@!Zx*IPkfzF`J}FhpXd9hhvu3X^Mak1)cJTdmOU}>~ld`h< zs1JHYdoPTsm;(BU>h15o25ncplaC&ZnL7{YgAVA!C2ugFR^=`1?fMbsi(N8=%;+)O zjtuJR6$ki^&w`ar33EPO4$B*9Vg7lDbT+mr+TC~fpMXIenoUDKm;wi%zL$jwhoBi> zin~kFWR2aPF;*j8q|JPbSkejQrBrEaDSLf|B6o$n=hG!n0|PgflPejcW}50)<)nXU zMkDEJWao2vv6y|Y>+)*-O^&>>nf>hyYFDCGQFK9)xgtt4*shI>A*=$ zTW3#NXZs^Uy7z2-fEpAGo*fjcfxAmOKJSN=_I8;DzPJo9%34Ova74)}t3kEq*P%^6 z&S!_+i)cMnLWT9%q@L!C)_2!2 T$!hhPSabhnUgWb-J5=twj|z# z8D_Wm>$E))(9sjduXABa7QRYTYVr^11xcC9tv#H|w3WaY3R8DE!0wzn*h)4`9 1{IF)-0$%5 z@qR$J0)adp`^Yv0(#Er#!cDpr5Rd%MZ=aJo-3cBB?eQp gz|te?vb?;kxM^= z_JH0t*r4?ci}`cR4l70Y#F2oA2WEgbK|JP$%jcbyd;vH(3tHokVnSo`r^liFW{fPb zyICkng26xT|DW@9QTQ!zG2q8tfd4ro{BI;xvO#_2(6eH&uH2{$0jQy)T(a&UQJ2Dc z0gnWMIzFePCe9KPQPH|Ge>-YK%Mr{_6PF3u+SAAR%^r6|p^Ty7efrq=0y%{3@ECCP z3j&?!wMI3T7<3>IDEEz>kv?|*8xesu6`!F)>1Sp?K=WTruYy~=fl`fG8&TAZk(4%8 z!0o}8yl~mY7jURLxKgq;`ktwe@stC48lf<#UHDE ?StH3?gIh%TA;HU;W70JmcSnDlyVm0yc=yaD{;TuqZBK}9v4x=% z4Tdlp)qSzl{2>3ijuq%WtA9z{FuJ?$FkPJYEzMU5cx*IRkrT1k{^1ZxuN)HNhtl<` zA(JEKFHsAWfYz$-BeJ7mlq`=z1?z0(-6d6FC~QqVy+2q+`%bU=VM_+nrG-m$J>7dA zcv(f)y=NpUvh%hGtkl-Ns55;w@1fnde5e8bTdlMV!nhl{iUQ7dM>DMzB2vLoa@sls zPn1O6?NQ)F9~nN%SMBLK8E8<*Bk^8c>5ZmuPZdI$B$yn*2G#1-KYA;Ex6yc5K)3(m z(EUH^v xc<*7!Y z=F0>H)Y9>HK(K;9C$lV{jlq3UV#T|!h+2l4n-W(cRyDDvQ7;$KE@{($fyKa^)+;Mu zNjbifI)J^Gj ;?*_m9 Ha8ZE@{!X%M+lPV++u^Hm^+mV$e0HxyhRzeY3VH(myJl_FRyUd zZX=&yDws&;@MPrcjoi`Gtzg)laoAYnh5;`KMP|*$#!xqWWVJms1!js{{uz|lCbyAX zfSN9r*KbSp1nw-LawpiM=#n3CbwA(onfJ=)_BFtRN1*cch8aW1uW1X(w9jdWGqd++ zuA~p
LS!dow}vOMxae{=w0K z-HS_T8t!3J%+vcHnQ362t_*Kz3kw#t3XR$I!=Wz;4UFol5jrEM667`wZC=CgQooCp zn)tZDWH$q@`s!D)BfPo E3kXa?v^7^1DR-uue z+W(Rq9R=Kexem27UpN&bA&zc6M{ooDJ*1s4p|Exs6i6 ~%?{C$Mqu5u1gA_v+K?zEf#Yn@5p zul&y+livf?w95K;%i@lp&3+f~YRehI^R~5d*(Sg44bAroO#aI%M{ylQ$|>X&hp2Fj z|Dfvh0l-y8ky8o%mjnP_^%MVrHgc-5z?qL&Adllm0pFBnCsUDbS^NiaBojH4+6V+{ zDsqaf{*2Nu@|r(g5cYS_U0n*d@9;zwRq)ONV|Nj%c0l<|UhwRkI59ySh`wvXieHFZok{{Z# zyVlGPhz^i&VQu5f;Anc9FLMRQv?&qrqdWj`qI_7%(gt2m39RJUd);;Xw+q*SG_z9I zzGj17(z>OiyOA!PCJ0%TIAFLoHD`C2;YGqI^zbRQ dJiC%mnYNFj+qft;~gg B} z@#{YpRlN6Zi;@N7Roa>PchD5*!l|8ED9frQtXpASOr}YNum{>c%+U=R_=PVHIAj{w z1Vpq_s`TGFqvT4nPIoraxsw$Sf2aaLl%fZIC}T+OY@fKLD8Ea?Sy1 QXt=LePHSADr%u45;~YxaFuZ3_PmZ6}_!PXc0GEWXTH+^obPj<9B|sAxy0lmqk24 z0a@CBkiTozF=rYIjY*%zfbmwL7APy(rYPM}KHA0w-Pl#9Pc&uFgbEDL6BZ+ T+v7C2QVsb(}CX=znfp7+qclqm!a{@e#UDdu;Dl RAU3|$mkyF&uD%{-Tyx2C`3Lq_Gas?FF goax8QCcm3d9NGRO{HF4w{K2{H zYghCe-Ge{y;B75zV${3$+WaBo?0yA^b^9P>r$cCHrS==InKCIo4B4rF5+=xP*gpX0 zMivHdr-2gyMizpyiJA(gagnZj8&{fO!JA&v0T& c>G`C<;JNHxviFX-=&Zmj&)q<1HDV zz+LKZ-S{%%OfB}E!$BEMV73K9a>r5rA<;|Jl<|S$L$nMniYE)8bCs9lWzp*{qpeMF zlh=Qb#|Dov(XAI}nBM`|Ge*^`^$Oy@>Zkg(i(9k&>%%@kTrs|b4VQ9acwis!W*}+~ zJlgVDPp}@W(Rkp!9RtERN+Rv(etKgR3I+yNoPnCqCs^lX@hLi}xtI=`KRoI*wH{0E zjUa>FO+qbEIemICexs3G{c(`KN^$s{wDs6_YZSFkS+2<)m4}kZpB`@aO^^f9g#rDJ zfYHvQ9QPXG(9_?6+ti!5G?sqEl0}Ou!3@2WF`Wr6u&g;th~Qt!76R@_=R^bs|IkXc zx+KWeG!zTwrCX1$Ps-P8W#qj8GDc;_v (#CB5+J#PosP$G2l5S_sMo`+ z=v$^0{V|Jxb_mX>?fXgUyq4f-Vp*sFm&;8p_87%x4ZgdHueEV=<(Z?BkG*d48Z-Z) z7H1ow0WWcDzo>R?a^*X1d0e|d703K9g0`4IeUMRRZa-ezopAaGpLX!d{wa=9|MI?l z(U0;=aTmt|T50a9Yp2H75=iB7iTgvJx5eN_D@F%6)jSlNjUYl9@4y8nWnL|a*?-$+ z`CbGyYZ8Vcr!r&^ K F#;|vj^u0`5R{{0_l?RCIH*?LK01Uy zsrkARVPDY&ao+orJ$V$!&Okh)=1;f$=`|$J9W^})6w_;*yJJ4<%baPATH#IdO&V @)oO-LRdTzq2QAC;Ri{WO4M%1Nip zQ7sr_lx^Zl&Q33QU$1Nnzx#kBV1?jONX& `n{*beDU>DL5#);0k+tw363tk zoId_Cc$)#N>Tb%sIQ~v|U5wrHA2iGeurF>SnR+{ECTIk*L*l=qxsi;^Q?1zH*2nyd z#Q~CHnuE7jfU^K#BUa#dK{;`@w|%X%)O>NEeNEZ&N6G5G8z4<|BiT@?sNM>(0=8P_ z@LWY^ygucQ@as&Ha^^k>6vm@LK~~s%7jU6*upEoat4=!{6{O`s5eFlnC?X5VqM1_- zW?x1AU+%`8fvudJY*TE*-i7>eao}b@J&p1;IPVVicn?yIuztF?>y=}{^h2ltmBZ77 z_FFNEz$+$<`m`@AUh)t7-~jumiRf;KHb*_Of4%t(4C~jK;w!hR&0iG+)Q%=mxu|Jn zjsdk4atQFB*h%Mxo_-TXVH$ZPL-*HaN0eAfR{t!xOHRA3I+q_WMhP0spSEj6p!v_H ztH3z3&<}Bl%6SW2b(tHH@J@t-t+Z8WQCfD5czkcm;72gdqy~m!Lvnvuv|QFdj-_iz zP>bxWgodvmV-S1yU4tc>p<6G=?2imVd0I!t;}Iw>dMTIDy?qU;WFC42rOez(5m)Cc zG|n#AUT7c&b_c#sf8Xf8(NGIiS1~|C+4J1VrhlWM_ToY}8e)lNp)Q@iQiCA^l&=yI zx^aC?6xyC!A*p3CWyMFmCXv?yJJs(t5+gIMb@ZJ~L!l^Y-(_V(=)zmf0SNR!oZHY1 zLH`@Rz|ct4#5qFDDO^RDjSo7H{?Qs`3`TgNcEQX!ZYZDP4OszhClxM$=G)M%+Jq~| zl&Bv-^OjJ9330p%VMt-tIcItoRLcG`Hyvt_K21e$UqkhY%nJCJPT3Q;whS&3Mf7HT z@qtZh1GM$3YfFZl*Z7yYh}uA_HJ%zF>~xPL-+uQUY&48SEfIlHqsnODpay~h!ahBI zq*h<2A$a@X@%5#gqegD5j}#t^wS%`r0cRlN6dnL(>0rFCQVw6lfIbkAyprJUr8LVR zP4zSuYgr08#dEz(EaH4LMEM63!A%zlYBs{fv16R2oAP}k{9d5r*^zs8c92uCr`{1` z_OX-Q?Jg5{ZG9|&p+^n^B=pWV3BH?jT>NAU+K6jkTY;o`! zUmLd)(9 3Q(H6~nyy3D2JW{3rO#?P;Y0XHtIej}vygaYaH|XgBW9-mgiA zi{=Zb-o)v=S2* l()qjR;GvAOap(?E=#gLl(m!hQYxZlRv_E@; zW*GsuBV*j`u+$Z@Kb|()lF!t`FOOka6snWIe#e2^7|_P~B9_Rp >I=uaTm;6; zJyY_Nk-LgtEcZ80r4Q1nW5lIZp|x14dSwUC@|M&w+^G3YzPONcnHNDwGiHH3or<;0 zxjP~A0J+aTmiLN=hxteh&RruV!00~V-~AsiV`9SY^FACK5lkad-45i$EFzg<$L)hF zmNmY9)*Wr9eJ$d!{k4!>cYB`-sqyiuH=kaJ7!SC) C_e$P}3AY z_HR?IcQr;9u$%}Su$yR(v|%>EQ`Ys5k!mMfi?-g;@a4+|Zrm0U-Ib((<=JsfJC78K zRR)pmAmOV=>Z7n6oRlD;j;Ppywp=A!jfgRgw7(DXa_{C-L#AhYFlxrBv75ATtx~B_ zUMC3>{v}NiJ02o7`NG6c+-q=4=4%>S{PQwW#P| MWsVmGC?FXG7;gPoa@e_1SX%^KkgXfyY<@8E z4@ct!q&-oCB)^)VgOB`;*+Fv?Avcsae|SSgkmRwfC~?aJf^}r}4hzkXk`7`M%8ORG zSh4ak>dm&$Z%%we#T}bSEA1#43IWe)ZMXU*?})()T Dg6+vx_)_@PMm`aY%lB zGml0~#5&tWrB?NcVcgy7iP&YxO~I)Jh;(0Qo*aIrqxE&91N?(d183+ar$rVcSae~w z5{UcmCZ`K4KGHOqaGk!-+*Uy7j5J(9TZK0{sH^RbBeZ35${0KA&5uF-nn4>q*XfHX z_mXVn-UwFyg$H;Pcj*n1vA3aR^bPmj6F&j3c~f%Gh^hj3_229u7ta=5zy}NdAw{Ez zbrn)IIb}?TS^Q(zS2sDQNSCg&Q}6Gd-y~TjZSDV21%+B1kBe#DYQ(Y}NzT|PiI*$O zv?y#ul2nv6ocUWlUt=KKp*@9mR`d?I-ini}J+e+GaCnX=ZvTY$=&)k?Ml&r!wp=jd z!fZ_t_g7{iZM3MXU;8x?@UDw%#4;)CNtIgya08L*gzl0kv?ruIOew6uz`Vf!^zXzd zTf>Gi{@)IRJVIL@q8`Ja4#e{nzF~ISrS&p?U-Uvg4+1V}%E*7cmN{or1#*(7vp# zUmoy##RXpdQTS9Gu`TV0^{V&({T@Z<-AR6%z)<2?{bNaL*ZvgHtQ*=twt5a#Who%; z!x6EPe5c8lNEi4=S~c~JQr4Tp{BChPPrw(ii7)Y{X{Si3u+)l}immB?>0>=l7~dFq zo0&~Z#6UD$PYSGiqc2Ki07 mFt hW;aZTCWK>R4)Kg Sx8G^FD455!Y&8(KhW70fH>2NH;M`R=TI5d89Uu`F z=NqxjaOM} jE?@T_t{EJ$%G85=i>HZWt^VB$R^(S Lbqs;2nbXD<=IboU>{B2aMGvUT&sh z tqpYT(fGD3Rp?|A37zcj*){sT9lK<-2aIJ4=xucyY#FN-;5&u`$T}m zB3+Naw`*KKWOxVrkGOT$^5d1#MixKMrjrg{XY)c{2`~wSSr8sDEH`NIE3j*8jC_z_ zSb>aS24W33uoT{XChuP+gs;}?LZ-H8J0vfWmva>um@ULUxO^v~9@iYjte2E%rQFaj zY!3Q|(9 $K? pO;QJOV^^Mf$5Q~8voD}}wVLHr75Wy Jg|^^A> J<7(Y9KaOd;I~;Qzp7@D- zUQ#GBhIy>WK<-!2)+%iU3$e4pE+f-+*qD54ZYr&l^jKR}KpjQ`f#Ie|c-ML&hT=3i zwI3M)qmGB%pNg@mKST`5$t(3ijgbmM9HHq@u zS2=|s(`w>vnyti{o9pp*-*|<)5@p&dxAG1weoo1E -W6Sg&clv RYkB*4LJ!AO_*JjJ>AtGt|#P54H@ cJLZ|3`nfIW*pxuhfm;}x zFBz?We)8tdyg|N`?P~6N!-eHY2txHlFkq-%MGfbTL|IR7asal|rPZHLVkL+R7ORy6 zAg!TfJu&xgZ}Ed$`z&F+Zdx)aBRC|f#SVvCZfpk9)ee}%50{SLE#M*piCBR2NV~KJ zyt $P}> o?EI*0Ug{U|a3cBWVWOLQ1|pF}*HV{y1k#RPt>DE%!&-Tr%l88n@GR_>WNIpNt;A z=Bu`2slo5aYvv!X<(`H?s_73%%yB!26oZ=@bUT3G?x2V52yO#LJ48uJm6)w|J(eW! zkVN?4CqX#|v}B$F(vM-;!vaEblqKs8CoHTU3czLL3}`a4pgd?i |$ZuUsE~+(SJz> z7>8mFSer!Rx??DlKmRGfOlOZ;<1>`vPIn?>8or^O%MPr`rw0Wyao9PSGIrUittj$3 zGmJDhry$sPru)!N_uxz)G_B zt}};76?00vJ#>%ox1k-X(LT(?FfpO3Nx?qK<5Z9v^SJxkas7>#qIoSRt}8XqziKjz zJC-fr$tR{OBij5^5xt(XN{@~RUFhvOXY#*_8R_{gpfS5Vr#$;Sm&qHhHm)#wPPHI$ zeO}vL*hQ_RVX?v3_R?ndX(X>9(wA9E%E5{^w(9WXCh8887-s)9JZv0L)aMl7eR3yu z&epxAe~0ZfnQ!fSQ@?tvNP=JB?p-X!a{66&gUS;?8{a?M>EpJ(YHYc P)Fn z@vh8ov~RRe+NS|(l1U`kSs}V@nx`aJv?4B*>#>`Nbdlf(ofBT50A1mMsi_4?i5NFW z?sX=Owo8P(Eg?;6!hX#zoU^04aLuoRoI1;7GhkL@-M<4lFJ*Rnowx#g7qlF*?>FEq zHJO1vr@P?17AM=cMn7>&HFJe|Nwh&)*yKhU3yJMg#^%@4cwt2=ygOI51pAiqE1+hk z6vU4g#H@>dL4kbxNP9)gD(aHL6Po#!^Z<0_yz;V$g$oERPfYn&<`*8Qh#=nNbX%k_ z*UtX{YKeHz{Y6(pfFlyBKN{i|HMIphMa6vJ@&8fgw*m%syX AbW)p|ufqL0<<~P{D`gZrE|a!+{i+Rrnp$5E7|m5THs6Pj9|(f!*ZHQL zY@eeo68u^GmBt;|Z{(Qjqv;ZL-;|yGK(KXqdqC*I;#zOK9j#Q$Z{1%!MHork?J1x) z*}DtNyToCQ@_MyrLEMSorUXxq959N=A2JR+dyYR~`K_nd0v+&Cws?1B>lm!9ZHr8D z6hApNMDl?__%C>f06sS7ZkJ(aXpQ_ttB!(M;SB7NE9P*4Gyr{5Q;lKx_K$3~ZR`zw zWy53qUYSW%fs2#nTcSmr3>bAV 2I`PKFxo(gNFB3!Q^_xrU=v_OhsoGGS5;_p zf((&!#$!B^Em`m0avcv89qjPc&Adz+8{_J1E1-k|f112Eywj_*#WQE>Vp*QpIq+^l z;`TfTxDF(siq}qiB2wBx4jJ24JNi3c*B$W(j2bB~idgyZpjse6X~(e~!<+iZzq0w0 z`tk>#L1ZM4HBNMdL+q(GZ9#CJ54LpLM?e)8WzEMUd11rtoqHRBGEfCc2FtOku*@%% z<#*OPyB`ppV@&5rbg~L??Uj*4^%;Ic$CImoi&G5zVtX#9tToZvNeK{lK@lcbE1c-s zao|fg?J}0d%sDwHjlq@2L=GJn7~#OLKsT0MtjZ^cg1Zif4)A3;>?`dMs3&n|{|yi0 z&yB& Nj@jFW*yxb$(vq?X0jrrN-3JX2ZyQT_peoN!{*>VC13{cix!jy7zqS zc2{hxAwLnv1$O;UM|gks1Ta#EW)-$M0#2KF9C8dO2H+|ST7Cw3`0$IRMqG3Fr-^GJ zzX3bm%MnLF9sOWnZgI^vB(c~vu?W_v4-p_WWVfim734~kKXY9Oq;-OyA{`LhQ_|GO zVs2Uyr^;QzN4@lEH&e4)4uL!pYWE-I>-m2$m3sKJe2o`pL8Kj!grwpS8?qe^YhW10 zuMO4$d9_Mjdh>5%79`fcUq|beJ0@Qhm^BC#CTMd1{+AeL3xnIj^y{tnW?^sPdzA65 z5qQp|U1JWPj>TyMoG9}>x?p>{-;lT6)Z;NDql|!BMU4&3DgjR65q{zd-)^A);7C9% z%Hv^WGY!M-GqxAa#vMI8rHRE!7u5fB7~ap ;jwd?)5!ON;zhYU05Uc`L;unfSAC#wHPHT%C2RUFgjD-UFR&lKL(l=$UL zX)lm0$}{*J>2rEz`xcGc7iGT!gd@mf-!VSV8XxR3+*!kfuTUnn=W^A;G?HBJeK-!N z9nc@h!xb ;KVmx>}w ztq7%xsK~xVr4>;DTb2M}iJwA{06~$kCQ+$StHLLsVGCeEmIMe}h!B#XRFE|&gaiTv z1QQ5^h!6q^N#>26nKNhR%$zgdnZG9g!8z}}_kHg2{GNO7bKhhCMdANBYx%zgg#WhI zRlvS_no95lQj1x{0|W$No`AL(1`GxSvM&VS4c8kTD~V)Dy`(`UVxEtAwX!iWCU*KM z -*}~s}oVekKH=jYiuMQjJ=!EiC)g}Jqq6H?^I5o)< zZT?EUV4~DKYLDG*1RxMgJ7Y1}PvI5I6z IZy5>QfwA(oiDzh_v~7Wk2gVOVBcBD*ErTC+d*=Kf;}(cN--JcHydSAW=~ sZgp(Ln}c0^v=Mn926HQtA~NL$8>7sw>EXe^}m=G=3fPb_dnLf#x RU$otR1`4aFCRYjm={sE zq-RAMo;7QO69*xptOM!%v4@KC9&xq%$c*p2 OgE5qG^@1;abn1W+!zqdJt zDZR23zPEcPfYf9_5v9`w7DTy!!%CRI>eChSs-UL=23O($sfin4)(ERIM(5QTl()3P z$q#*t?{C9t)|B_!Q|Bxn o&RLM0(As_q1%&hbA8ZB9~g z*(O8?fP5EjnIqzzYe_9kV9j8shKc_Q{o2YD=eM5ip`aGZs{v$WqdTw4`wPF_cC#V3 zG- UvQo}mv zIO)pLGpb&_h~h( f&0FmdnmW?9u)O$a79@jC7@=ZOyp z=61@%>i@U(f7r8@uyFlFjX*=!2gQhqO-r4>7X%3mqLn8?By1vy&HTkwP4@rdXWv=V zy#03MtO>E)+mN{0AFzF43*4Jhag4MIoRt^02&(9qh~Xb#aE>>C+9P6Fz2r9c+VYl2 zNtZA$7jXwa*$7Nyqd>5ZF=i15H+?@g#Pd>JTl@}??e_G0Yk~ zQqd zFIu~|y)XUi90}t5@Xtn^e=oHBvjo&33&b7h)L#F(DpE!Nbl;D78{bafJWs)5m*_i& zJK_5|`{X;KU)j+Tqhmjf-D}gj5hbw%Fp>S6afmexO)r}we`(W_b-<~oz$xxvt@GN( z=N#Gph4TM-8nvEyf5)aoKgoSz1M2{>+$+v^vbEN)bH14h FFV>liarzc30 odnU2;^1J^97K_A zW#VN?mC%d@Cm@JYFC7Rk7+-g9x&0j8Z=5QlrwT+|`!&b6B1*inIhY5Om)X{>m&MTO zRfHxo6?fb`Vt6Te(dqs7%MX!7h u%DXu$qJ`4e8RD zxI>#}m0V`zP4}F40Bvm1+W??gwfS<>gJvnNGXjR$^gZ`wfV>wEl*Jykg5|mQ`Kebe z;UGKd{DJ7Z_QVu0h)$EaunJQLZKK$38QJUo`VN5HddzX`X!O^tj4FPMTWohYJy5f} zE*T#C`rvX QpUKbxP!LdlEk~>~NW0$osg;79l z36SqJW8hBkaq;XM_Dz{}3C8gH+F|l;!M(oGesz*0*CL?C8MDIb75z9CO`L`JC0Y z?eqwSKNjqM{Y8gBd4rcZzqek|2V}-=b& +Uw0qDJOe L4>*=H;#C2E6rU#cNHC`uP8b$RFSM~tt0ML?>H^Mps#-SE(}#qMZsRrGLSak z@y|y^_%6jT(nM{5@=&j0q}3&>et$8jXBGs|l35sk$>O%o)XLsw#sUFz U>j=)3PV=1eg>3!+T6RD8xQD2HQLxR)uDb o7C!m7?``LTiW! z>o3B?E^)q&Su_tm5ze5 jQK!UPLYTVfvGi6&UyfK}tz(lG-EZ2#N(^ zBp`m%y8?O7+AJATLtm_ (CfXk3MRZX>|9ZzfK;T-Fx z$o!#HaJ`}-6->vE)iTHaNQ6n!y;Yr73y@cO?xz6b0nz@@NPEIV>MP=6KhEiulg829 z&nhhh))k)KfxXe5AH_MwAB7t yvqctURXO%k0zem~l#4K7!4aus%>np&=!+A}1k$^<*4S`P1b*`H) zqw0YO&!eO(c@E)>YoCLaej(MG803j&3zcGTaQ4S6(e*_%e +ewpRcK-_0=(xJ+}s13w@(h3;P7AeK~SVeqhN`zvBg-r+s{pE#mI_l@fqFCh62L z)8IK$6ZiW*v4wJ=1J&?^4Z_PMn#F|R#Wr&s1vTS6qGVO@T(}z#Rs)t_h4x9t92-)= z1 Hc$PHfK<_jTqz!Hv6e?41){rhF0A(f6yVEtLGv2UfjQwZA`70S zu6Vf`9oU{lfMm!RsH~WFA5J%@u}LaFGr8r)yl}hz+Tmrsuu$LERc58XVH_7Mp&=eh zf7J?b7V!)_6Q<&GuNK)-Q&hWTJ^3=7wk1pUd~c%f7loqGBNK+zaYZBJ(fC)uNyjgi zq?LT$v3aVc85{j(gICv&JrHpZ@FdRzHTS64UZxlL@;E^Q>tI8c CemQopA(X7Kh06zKVI1! z1W!K)Aca#BdC8*%i@3$f!95wMq?f76kXlQ%>(t*#kn)0Fyv6EWih!i4dZQbey4p6c z;+<4Qr}ff%cD1%noo!FbcqbkKa4zPcXa u&wR9nnBL_kK`oRtw$J@*Dr)czng zX$r&G>;<%j9j*=I0LjGh=x+JA?&gvnBE_2z8JvOe((p=W*fx( ?37BKpsUHmZ1Hi_M_g z? Id=R>Id&^;ENB|4a|hhT?QSx3RB+X>p~O^%}?_x!4uE;>y^ z QhV-{ocTtlUV^a7RAwo<)KUxT6}I@l{4D?664ev~Zp$O?8&fy7XwodfW%4JQgc=Ia2Y_-SXRt39Mv% ^p$Cz9YRj-W$pDIP=|54{M3f1Ien7vD5(%y&ECrw1;0pH558N;u zs#9zr)(*Wy>ta%|muU^MG?tADVY$erEHRC&iI@r*Qh!~rUaRGxGN^sIq+Tm`UjbP; z*q3q3{uC%sjKZE8U3^;734ZtT-~{&iSE*LQ0ncJ-Rfiw3%#hO-PbK2?uM|LmBmbhX z(j_-3Sn^ajR zu&0ZH@K!t#&&W37Y&E^}`L9&O_0M;ih98l4-8O(YV!JLIQolM)>Yw>Y#I}}&(2jHt zRe-aajlkCvMEA1y-%EQq+dk#G81?g2mrQLRhKI%Rr1 D zuQ=m#v$-$li&qM=nONfwR>E>EaxMyHmG!#8GVX8!Ix{6v!2hKKWn<%TwVRx|*a T~uCebQf=%+tQ~>8_ErcdV1`AMPqe5bibad4p%DwzPBk&!ftqznXQQPCePb z<*12USW*J3$rKbQJoZrz5h2glKeJ gQ5JStY@`F7&)K_L zADy%xN;BDk2S2xqbLp;r@5avz$a=D^w=gL1DkSsjxyF3gPjyF?3#sg`g?gWNiZu%P z!~^HqzC)vUe{!?g;mi6g3bFp=w7SMpTT>kh%FA{Ok&)Jay6~pKqh VZ>hc{o@v%xdHpGL|>c%!CN~b1w@9sxHA+n}xc_yaCj)km+)aywjsj&CC3? z1;L-G` bnK0->0i-kSQqL}hf-E)&jFEcqWq4` zVK#8@Ckg%Y=NjiS;9T)L?URb7Euht1k*aUM1i>QLk=^MWvD;0gWFs{jieu!=WVkBv zGu3{ro*PD7Ppiz_nQr0I *XS<;F@_iHjWb7^HM^j4$ zD}fg>#pEyzki%)U&R2r~ c4c)^(gO7I2HFnF~xo$kh(PF#Ke$gBpRJ<>Sn#oVCLKOmHTj!wxJ$a%Lcl~B)N z=|*g()Wp`Y>rjJy-!Ot>RlTXjWO^Of xlIGe14PTMv%dNy-^OIcLUsF;wxtlk<+kH>|H6pI(9Iwmy5; zyu07LbW4TnekR-_40~D42L$R_?aqdU9t`mrwZTCmUHXQ(6ZVWFJ#r7>k<4W5Y;%SW z8!P#Z9a%McwDeaW=($F`TXQv?hk`oiw@gr=tBM?
ywYvlKLirk5# z8)hXc{XXOM(0}*e8jS&xSLCrbYq`ap>Tl%elrGQ%dApK-x=>6(W@dSY_Oq@G96S`R zwH2oL+d{#PsI-Xz>Z+x)sL8#8Z3S}6%z#TzFT9wL&vCPNJq^m)HEVWF8hZqu2@E%F zcI;ZVSPKwj0yD794@^iKH#g)yo9`>{#VlY=LOt(}$*8senCUfTUuZE@A-x48*PF>- zo|#R}#?Jy!v+w6kvk~1#2QgG=mV-$bhPE7IN*Ps^J}diNA%bb&(!X7pJD!q%FPqC2 zyKbwOr6a>A7F~~MMQ*WsGXosTGGaliB}R~pN>4+b#tKR{=+?0=ZD FC9^T43P z?pc_alMS78p4<>UaP2GTRc2bl x4&%6if-1$^9=P2!p$gLwYVvX7s=S}(DT3z %$M_XSx^K^i7LK4p^*3_?!f}b?+-In5j2nm$y>vlhCa;CAc4Qe5z|Y8AwF8 zZYhH(lPHXI4k!4M`ZirR1eR)l3j2BQ9=d&q1L@MIRKWM%6ZlhRa&?DJThQdi5%lg- ztCp{hUa79Y{mg0YMsEW7jNC-+1L|_V2!qG K^>u zX)_4%kNw4fz!Na155rPmkePW`R+iqAYYZdo#ZL`Jyi40Y3Lupk02xp%