From ea4a57c6a1d579a564b18cbf5c4189c57f28da71 Mon Sep 17 00:00:00 2001 From: Markus Kalisch Date: Tue, 3 Mar 2020 13:40:07 +0000 Subject: [PATCH] version 2.6-10 --- DESCRIPTION | 8 +-- MD5 | 12 ++--- R/gacFuns.R | 2 +- R/optAdjSet.R | 108 +++++++++++++++++++------------------- build/vignette.rds | Bin 308 -> 307 bytes inst/doc/pcalgDoc.pdf | Bin 466491 -> 466491 bytes inst/doc/vignette2018.pdf | Bin 311635 -> 311375 bytes 7 files changed, 65 insertions(+), 65 deletions(-) diff --git a/DESCRIPTION b/DESCRIPTION index b66eb7b..becdb7c 100644 --- a/DESCRIPTION +++ b/DESCRIPTION @@ -1,6 +1,6 @@ Package: pcalg -Version: 2.6-9 -Date: 2020-02-11 +Version: 2.6-10 +Date: 2020-03-3 Title: Methods for Graphical Models and Causal Inference Description: Functions for causal structure learning and causal inference using graphical models. The main algorithms @@ -51,6 +51,6 @@ Encoding: UTF-8 License: GPL (>= 2) URL: http://pcalg.r-forge.r-project.org/ RoxygenNote: 6.1.1 -Packaged: 2020-02-11 08:49:45 UTC; kalischm +Packaged: 2020-03-03 12:31:36 UTC; kalischm Repository: CRAN -Date/Publication: 2020-02-11 13:00:06 UTC +Date/Publication: 2020-03-03 14:40:07 UTC diff --git a/MD5 b/MD5 index fb0bcc1..562b4bd 100644 --- a/MD5 +++ b/MD5 @@ -1,24 +1,24 @@ 1c213fc1c5a5ec68e7a5ba9293440dc3 *ChangeLog -37d84a7531ceea244c781882060b0d93 *DESCRIPTION +4b9b32644c238bb95ec960bbfc0bff0a *DESCRIPTION f813efd2d2208b900ae1d5ceb8dce2cf *NAMESPACE 66963591a11cfdb002f5da9070f1500b *R/Aaux.R 83c288d2212114e54bfe624e0edecaaa *R/AllClasses.R 9f99e28f01bd00dedf565c6e01978676 *R/adjustment.R 13a57456f8fd3ae055dd9a45ce6f933f *R/agesFuns.R 709461024e0516af39a962661e2930d0 *R/deprecated.R -4cccbecfa381ea1ad82c20e0151b239c *R/gacFuns.R +2e96c58990915d4dd7e51976b224bbd2 *R/gacFuns.R 8eff3a164a5b32a9c4818a1e840199c0 *R/genRandDAG.R 5c33d7a0849fd620631eb901d6a6540d *R/gies.R aca17dcd6f37ce1136086a8bc285046e *R/isValidGraph.R b363657f7e2c52dc0f2aa083780fa0f7 *R/jointIda.R 386ed9be1a524de90dfac185433b009a *R/lingamFuns.R -cc462e22779cb4c455536d1f28811eb7 *R/optAdjSet.R +ace998420ba6145fe7fe3c75c5892e66 *R/optAdjSet.R af4239ca324cb48eaac6b07fc0cf0061 *R/pcalg.R 67ec62ad56d7cffaba4bee8f88797136 *R/pcalg2dagitty.R f84d95dd282c42c32595a27ca98f9a3d *R/possDeAn.R 0db1ba6d57801a4d8fc6cdf996384a4d *R/zzz.R d599e23940bca0203504c0fa6b09e01c *TODO -bd68b3d89307e8260f55874d7dbcf993 *build/vignette.rds +94bed761e02266860aff243b04717fed *build/vignette.rds f832779eae511757a99906c7670fc0ad *cleanup 7cbe431a233e6c454e24f7b92d63fe75 *data/datalist 1d8d4c43b9cbb9dc9dd57acd59f565bf *data/gmB.rda @@ -32,10 +32,10 @@ f832779eae511757a99906c7670fc0ad *cleanup 4e33425acbef8bf36705d8a75622baa8 *inst/doc/mkVignettes.R cc46b1f49edd31faf5316d5cf4d6c1ce *inst/doc/pcalgDoc.R 6581d4416f174d4d6f07a2a2d4b9d3cb *inst/doc/pcalgDoc.Rnw -7c38555b6c1bbe4e36f53b8c8810e8f8 *inst/doc/pcalgDoc.pdf +277b633e6b422057891d5e6722d84add *inst/doc/pcalgDoc.pdf 01c5a9cb2024f1234d7c4aebb9762db5 *inst/doc/vignette2018.R 65ba145dfc8ae12f50c95d77a8bf8600 *inst/doc/vignette2018.Rnw -fe84a6e311f25d42e56c1f8042c8f07d *inst/doc/vignette2018.pdf +23165d2276ae3dd8fee95d261b9b2da9 *inst/doc/vignette2018.pdf cee0b4475c720bf58c466e1495451b22 *inst/external/N_6_1000.rds 0d84fdc119d76d992b72f9bdbd0e613d *inst/external/gac-pags.rds e7a2a9117d97986b2ab6afc6b4f77478 *inst/external/test_conservative_pc_data1.rda diff --git a/R/gacFuns.R b/R/gacFuns.R index e0dab9b..a9030af 100644 --- a/R/gacFuns.R +++ b/R/gacFuns.R @@ -197,7 +197,7 @@ isAmenable <- function(m,x,y, type = "pag") { ## check if there is a pdp from cand[j] to y without going through x[i] ## cand could already be in y ## pathOK <- ( length(intersect(y, possibleDeProper(m,cand[j],x[i]))) != 0 ) - pdpTemp <- possDe(m = m, x = cand[j], y = x[i], + pdpTemp <- possDe(m = m, x = cand[j], y = x, possible = TRUE, ds = FALSE, type = type) pathOK <- ( length(intersect(y, pdpTemp)) != 0 ) diff --git a/R/optAdjSet.R b/R/optAdjSet.R index 7ef9628..74c4e79 100644 --- a/R/optAdjSet.R +++ b/R/optAdjSet.R @@ -51,57 +51,57 @@ optAdjSet <- function(graphEst, x.pos, y.pos){ -isAmenable <- function(m,x,y, type = "pag") { - ## INPUT: adj.matrix m; sets of node positions x and y; type in DAG, CPDAG, - ## MAG or PAG - ## OUTPUT: TRUE if m is amenabel wrt x,y; o/w FALSE - found <- FALSE ## if found == TRUE at any time, graph is not amenable wrt x,y - ## DAG is always amenable - if (type %in% c("pdag", "cpdag","dag","pag","mag")) { ##changed added dag just in case, makes no difference - if (type == "dag") ## added the if case for dags - return(!found) - i <- 0 - p <- length(x) - - ## for all nodes in x, if amenability is still possible - while ( (imgOh`c^k%TdC z3^uD79=3|E@A z=zkY@ZOECBIQ}WVfD1 zPy$btX>Z=)yAe|N!$H0jO=ySRb;*A9)wb>aGEq_xWfNI2NJ&M)OlrOE_J61X+C*yt F005!!jJU2#Q99Uu?|gUf-e-sqa**pg$Z>V)&h+K$GwpE?1^T_MD9L$T(kLlf zgH0R33#e;YJU%`5_)zokO16sR_^n_dp`h@B*BRghSA>2M4!HNUZ=)yDg>ehl2uHwvionw@>z~udeHMmx*Xjlucy8FeNpKnACdR?fwApY(7V8 G0ssJH7m|$t diff --git a/inst/doc/pcalgDoc.pdf b/inst/doc/pcalgDoc.pdf index 614e293f06c1226e6a4b7ce99659e5828c942f9f..de0578778c0710170fe65a5c3fd2d189d5f6ff41 100644 GIT binary patch delta 4931 zcmZ9@Wmppm!^UwCK~zF|bP9-cj~F405F`YoW7GhVj!}~kq`}eB9l}5uDJ?M&=^E)l zVss3_QO|j=_uKjKyFdJ|>+bl*((#R@o}LJf$o@?Y-(5<$;~u>=*e;37PqwD>Sls^F zZp**{1&RDzn1xg0vZ{idw-Ut)2Jxp}ML^~n(bAVp#2UxMXi)C8;8@osR0!_nY_!cFK{{gm+7$%ozJOfkx{t1E>_b>V4Mw-S?_bm3N#wbJ& z9cXg$WdoQ|0f7np-^gWD*#G~ezR(^z^%c|`I^}W(r6}msK;qzoJHxC{ybCZ$^eSS@ zKuB&rmP3AkCaZvQhUQ6UwFzzmBG!?gJSj6HmfU5e_e$VXSQ$|jc7xGT+DwxAZ1J+R zQr|xacPu8?a%osJ!5#+v<6~hWHNw9gAE%DbPxh>JfV%wnTW@Og((Kh*S1={wGr=R- zJXy^#T#z{Dlnq?pz4{&Wt)`@n!5Bq(`kPl@CDw44tXDvPpZ(1RbVIT|!w_3i=abJK z_d0a=pH_dULJ4YfCu91IXWGq~Ul=H8OCeF5T|*ne-@=e{=I`c?F$bUKd!CS66FIJg z&r?$L<}Y4Tg9SrjzQeNOiuSYD?Ug>r-5MK4p+`3M?RW4$Ll;LB-788V$Ml`4pqwP< zicMUudjFH(&UX+=zTEAqsGh(``LmV`Z@N>0$0@{hAt)mt(uH*P#xmc8frlcRBf5pT zZT?#66p^1SSxyO5J)9kfjjnsnGxiDHcY81YUYG?OuxX~Ml%K6!c@CcJRsKLCD-yjI zX8r-}$^j>WM;B@Vgp*rIp(K6b@;h;IecN$*5*`nSC#o)aQW#e+>u2f*vhI%+Qxrzm zRH(zk4ETC!s7r#e*6jr#9rm*1X+G>Ujzdu=d1?d46PlMam5h~TB~sU!o^H~-*~6Lz zf#Ht3cOW_}<98WFJ&`H#8>zg~Osr(yy}r^$xNr2`n`a!VKzZG5ESD8+L7YoFe|M>| zxz&$&u3FIe)E>?UAB-=6fUX9aS@m!MQ6Hv<5R<@Ob`fq$*Synlg6t-b(U@bL7x38x z+^k=<3t&aB2m}!X2@C~K`v7=-zuUE$Dpz;nVDC!#jJ^JnYHFqY4x`DhADvUK>pxx-QW{cq$C*l6>RtItS`Oxe&X`iKn7mab0!iG zX*F)1%wn;7SL3}?HRK_^`rXlDIm;5XGg~hu>&d&CI@~!d)+xViYvZMx!=ty^;@R?8 zSii}WickDD5bkIuD=Z#jTxoRl=XZV&m%6uFc|a9g(Ht~E*Z<7y3AgK zuH_n1Pu*amdY$j(fl!$60Z`U5P>S1u<(gO7hzZf56sH~6iyLpaq71_UFS3#&FvNwrY`iFj4o>_9Iehp5mc(F1c>Sui{<}j1K zd~&j)uU#2w>PWaRn2;kyOB@N;A*wfjPY}c#4szKqmw~eUBqY0$Py17s^})~;S;?hP zPFdcvmM%J-!9~Z@sJhy~@qY!wYD7?L~ zH`=Im^0XQk=nn&$S1b(Re4AP|S=%%FJ1F%w zJrl8|h6hnxBB?-dIESmmH!VpB2Rzz152&^SX^AQw^?2SG2yO4W=%Nyx? zolki3MU21LO&=owT7c&t7xC|m{hCBiKj6Rn?I23#UF@NInk41*!XU-K?U=l@TdJwD z9LFZp9<(?Pb}PJW#_~#`ZVz6GjJ^{I2(q)4*@i$xlwik)V<&EDce{n|Rc>KOLBDQ~ zft5UQceQZ}%ByFldg;Up%DqBW?quGRmGcT;2aQ`P{+cK~PKFOT$r7MRuK_(O9*II? zs&e99m|#1St+CC)?$o~_ffm0TSi4wQfQ8`UNk%(;8TB_imB`o(Io6$GTH-L)mfjDR zhR6xzjMKxYFOGr95etvFKYF)tag2^RaW6?Wpen#y=kAXm&qg zrhMI{Ev4NogyEk2rozHAa-+sE3>k8>sej8~g~lv@PGk_=aJSWrElV2Psu)nH?%;D> z;|tXMZflxgi-%@q#FkMR?|vsppFgR$j}v=ME1}m^(rdy?yP7$|lSBfkqRvUfPfUE< z-o(x47_rmu4Gdd~#>0+_A!IFdq1ByfD{2M4bGAf=4e;aIV!fCDxO0HOlexNedyF->XTYAYQNRMvwQr7eI1>aNMm%1o5~AWiVpfYia)(-c(sz` zx4WgOy}{r1Y&;~*XC-jaQxZcJ$r!D{bgu^2&Pp!vnII@IYkTVTB-?Mbu|Jy}^QkFC zmzLxoW6)g+Y~65p{T+>~{C8ACs2lcY1+cOB%Fue>>MY zm$!Ch^fDT{7)128%g%a_;MP!vgZKb$sz2;>;wD!OF|9TwDT4@Zmzx>~FK6lUmi!e^ z*xWE2$UJfI#VXbt1Ms*S0J zvA7N_n=`yPc}Dyldh*t$CRima$;3}%Kb)!dsC1Ij8l(usFP0F`->*R_eKWFe1SeiK zzK*odZ&$o&$2uD&vi|y6`xpUteRMM^XG=4!5YI`M1LEqbJl{f#^)KU=&g%3E~UpTC^a7`oz# z9TKvuoDcoN;xf+=-NmiAP~N}ShBZAIeb|_q37<{o-!(H>vfi_^piqN1IlZ(1n^Y5W z$@!Ql$~RilK$=Vr-g|7ZTWm0F29#C&kzo$=jO=*vA*Q%T=FvE2K^RFR*33MgieISv z3WGeWLmGXum?b~-uAm{B;3x%4n?Dss{P5yU)8+SUu<~rN{>W6?NH;|-pU5#tS0`4< z{JFzea}VG_EWXFwRBQzw^iE9?6ecK-;{iLIYrzqc=+}85^P=+Sr52z4TnxUCnX|sRI{PP_iv+*Ja?ewD7y3}$ zG=05q#qXE&j*bu6DLRqn&3JgfHH=L9osPcWn_c~xM1pBYATgOaOof%jmzb|MCKB(D zcDGjh=$tXN~hc0+UQ8gS6VoU_jR&r1h#HlGar)emQ7z6dJUOV>CenBzpXmvFyKu zzK}n^KbdbuH=PW46ZeQZ)U2-s(1ha_zOJ07L8-Sie^312;NHL_SL*cn-sqBUCirYE zOL)_g$8hiPqkwt_e~*rfaj~eYWRNKw$nLUb zPR+`<^np85?`<>P(^cJCx4SiVp%!>47}!#Zu8dm#M=7!#bhcj1*LR$+3``|GZsK)) zrf-9&EK_R}L_c-+`!M6P3;f3WpiXcRXdoKPveUw=?V0Oy!$5_aWOH~}vG6z-m@d** zVLRvoznDio7~z4qkf-@StOZ{ow5Tfwa2Wz#r*Md|HA2nz%6dsrt6k&M4{tK&o7vX@ zW2k&Ehn6?ua^=(p!GpF^^#pyoKQBKtx|G_4k-%M7Fi^g_KE3^^{ZnXuGzV(K+ymda zF70^QRps-@J=XL4p5Se!mwNQ^bCLw>&a&&df_DVntY1AmB;=C9M_zaDMC2;91|F=J zw!&RTsX{wB&-{GrJw0BlBj;m^?HXMEX>0faSxS!zrPT&xT?p2s2^S%)HH&Drp=(9& zo2xz74(8&>%d}J&pde2z1a<7G&BQ`PxgvX>oE>&=BhrziTSIoaJ;1>D%~}Qa_GC5w z_hm=T--PBPzQ(^Wrf!*I`UaglV++s&OW= z$)|U)v%$3)VV_3dwi;)7$x3k3_G;MK*fSQY#E{4_MG@7w_bs_Rc@oFdbm#jr)Eyho zE~a8bVHZl-f9Zkx?sk`#QLB*5X-ahuytMRg`XTuc8S;)%PArSROei`cM4i?%8_wP0 z&RjfqTr*Ito7>uOk=g&VCT9E(d*{g#&Fw^`y_NIK-{AeJ!h@b0w9cYM*vvt6uS%-! zS!PEgbbF}8yFD@SeIQR5LpkA-2HXi(_Jm|9Cr90wjZ@rsvXbhL=CAL-yvs~ZrJB%j z1)!Os=f{uQxd4FG;=|1i_u&kBS2+2a-nM74L<>$pheT@pLC$7({yxh;Jy>r29bY_8 zD#=y_xx|id>s?0~{%m_3EfVi)B;@2AzN0r%{&d?SHs*&iV@h+$I&-vA-LD^$b>NK+ zgN~xCmkH$<%|>z;0lrr5Z;pKMYI?#dm?}-4O={Dy6I0}uLbKolZbd^(6dd%=xieR8 zN6O=9UdGzfBpmhfOxWEaJ#IC^3!@+19JXlnE61I-xNIQ6whG-w`7zZRios-{)keef zx7jBwtER;fSpQZ!P0}8%&k99(YLDy@bM%86JmSpl7g*^}A{v{eIiYuMm&c3UvAe&$ zjjVXo{Z0xTLK;^1wbP~;;DlEn1dy*YHU9g$?ZqvL+IhV~ybXYxL+wgpTw6dk zt?5TVOW%6oTrh*bF)zpRVm99=Ua@h1uwdp=FZgfh*zv08-?QWdHyG delta 4931 zcmZ9P_cPp$+l8a9-g}Q8y)UcRMX*+o=&W8NTJ(=Zj~ZenSY0G)uzE`@AwhJ4Aj)cq z-upN2yuUp&=Q?xGHFN&~r@e!%y@Tx&10Fz2wQG&l56?XSoMu(G!%fGI#wu4tudK}( z|Kj2gZ}c&@v!nn7Wdo;G=Ql@Qk!|*86dbz?uF$Kq*@M-a=&PZ-u)|IKJBrP_t-tLD zZc<(O*3tQrg7`O6duOv3;Q@j073%gp6a%f2l1^QG!^Aw1QqG7`Z9jovxZiP}GSK=~ zw4yKA*49h}7NZ@P`6sQLZuQjy=A^Pdrx1;q?XbV<@(Yqs7GCcQu#Ra`witi5xn+-~ zQ<5hkcYYk~duuoxeKUK4=3Y%#5gt3fi%sE;Q!n1aS1-;yx(i`=`=BLtD|mXjYY2Bk zyF*->b%KAL6ss%f)?%khdFEh%03c}>E3J|FM^AvEeoYo3<{x6Z+#@{7A%;*9fv$1` z!L(zsJV(p8vsG{EFm5@>=yeeTmTRGr3s&Cn3Ia=bVj^Pk)#pth`3ssdZoWDo>($-^ z<`JR%K&AdeQxRC-|5g8C!2g%HpYGO^{8a4L6Wg^RP)5AFnE^1y#NLW$s45Djbg>&% zs?gqmP&<_e;~pKo2*oXaJoPHlTEu1MNb)icvObo*{Y_7Myf%PKHbSXnh(AdOp~6}~ zNB`X@jYCCegTYRS?h->e9O*Myp~H*~hu7&Y9Z}b0k%^`#a55x$JQ&=Xf4=A7N2`WY zP7C_@&KAb6bqSDG7*3I=kSDGSSImAvAAZ*hTF>R^c2w>6N$q4)R5*S1cjYIhjuCNW zXrE17rqw9Xh}3ppgIl${vV836DEJEcalKKl_bfO(ZDw!Ea&VZ%w%z8Cg{k<9XJM7% z{^ghZKc#c+dJ@XS(=*-g%_2^OOSG3g4bFh2+wsq+%Ee;vICSFR zFFG&3!HwzV^^_U~2EGnj_zOJ{^CN&xy_Bi#Qvucq-8?xj@Z3Tl9+b?^+Q1$Ky z@toz{K_kbHUVlBl4;|>dRCvnFwhPg0iM1#91&zjHHfxj~ZuU#l5;ZOgLYO{`-rMN& zU4T`N>XHLazelKKh%s1dMx0?ad0&tX9C~hj;Jz-3!wjJ>z?Z_Xaw1upm@*43CUGWx zPtbxR;@iBz8i)Y#g@v31_s|@MkASlB z!;#HegPWytZsq9n_ZwT0g(JeDq=Bf7uRlbUt z0KuiW=aT`(lEGquevldr8~iks0uz%eFTzNw+nB8~l5$)sB(`Lqmt}b2l@J7jpDyNI z1IKH8X7_~a$7yws1D>`CJ_eh!;8p#cV49s`Dh^8jDhQnxOnz?LQlC}$+KOSXr!+_# z^LseXEoRL`M#uh6;Yn!Y3vs6H83~=tu(+3V;faJzU$TRm$(=~AUz}L$Eq#?6-T&n{ zR+3o~{T!rDn4|RFc?OANdMsZO)-*(cB?CU`d}eJ_=02D-r-W$i#wp_>{Xan{D5|2b z9A#V0)>;Jpta%XMrFXqEK%o??MJFrj{6o6^)AJ=@W(=GLW4dv*RQ08jWc^wo|4cJfe ztCKR1CZkK~&s?b@dF+3;ey@@scdFTvNvb&LZN6VR6efx1;R{s^eyj+gQT*xgPI0S$ zqO)Xm24+0}%u3mV$ZyqnDJtEnB0jF)$R&|tMfwI~{T+I< z7y-H;-#R?T%8Ixwc+n~|DWx+56*I2ywy+UPy820KGJ$b>-_P@QA4zE7*1%^b&l1QW z6U&1>(;LtU)+}k>r^Y`xi}FW`3dNB^mn6n3hL`@AVi9{&C#!5K0q~GmAZbP&TiA~c z5040ud(fqo#W5l#4*f&d{x(0Afzkb7{=}`GFvGxHA2_(0BA9G80A`J@KAtqvx8zuh zk1vfc)wr5HyL`6dKa;rs5*9mXFC{iO`ZDX_i9y0j8Sw*w^2YlY`1g*TI~9s{Z9y-) z9$5r`m*&csCDaZe*SG2e0rRvfbOfCWN|2Dtd)R_JvN(04<7&OdI({EP5-AVzt{2D6 zaST&$7FTs1i(k#Uk<%^7EtI~W?q1k(p54J3iIe?Zw^F<1936G;KY;P ztghqJ(S`CGg^3S+?r8mW4wfpQe43=xgvpNP0ZkGqKzPWu-*1wzpU94xl;hdFDgIo$ zD%XH0R6U#8zVFWr1nG zU&KLZnNu?yH|XL$2oO*Z4EDuZy0&ib_n4GJq?RT>?xlw{9K3YyASZiJH(57(jeZLa zJk#2_!lz7RB1@t+;3S`rs*@?!M==#1LO6m9Bwb!1))RwdetwVXxiWdzJEkULrjcpF z$&z!{k#}pExs5?|GQA^pk)*UWOE}j zrogXVMDa2QFfDZ5usZ`-UO^=3JK2Dpqjja^fiKOi3$QL?;4F#Eci_F1;gB@QmMYLt{m)}cak(Rm+cyV%ws zG*Wa0l{g0yv)TQtXM2>VfQENGJCl3cO3T;f>1*x`m5dTRM;AU!3x!C?E|h(e5jZ3T>*B^-3YVOS&P@gVKc#-IYe zLk56SHCmJpE5DDep<@KtaSGxCy)9-~o!Ey@U?K1G%By&GlwHw0Ky)WuaS{f(xR_eL zinm4CcC!3U>mfRc>10HLBH1^ki5DQ;X;#nW%H`Pq_&Z1EyAB?u(=-^}s@%+v$ycQe zv&^&|g%x;wqY|}x8Z`_a@?`&Pw@p}D2(Xiq%!))BLe0n?(hnAzCS!lmReGecHhgoE zv5PGax2`ys%1;*S1T~Rp7EC~L5Q`^${Cx;&w)#JVX?@}Fd>E+JIgrW2@uz!XH$U>( zAEmu>PmACV6@2c{`|0`HjqGycdt=#zp!vDoMq#)17G>mu?~?h^G|aC1lxIs0c(nww zsaO5u*kpViC?stV;$xl~Sa!p1dh4@ZXm@2-bh=tkwmeg|vSic$)NEH2^P|$&->MF* z_4;WNXYHirwu(ww6q9KFj%T1a(ZpEAZi2FL$q7p4^2=bk08}*${X56L@iac$!)yA( zGa$&rcj!@bxgGc>xuTBeS8JdHU}9LLY6(g{_y?jp&~>P}sGSxPm_qgIVtq))IFUZSp1Z%NAX+*%`G&>rax=Jq z_n?K3?4&I%i)ZJuXY_d~Xg(lV8O$ymso)@YcvyQ$t*tEVG4HADHTVKZnD%lX9MG3l z({lY$;p&R?)JUGLVBa6UcJ#b1zOFf*ymi4)`uINouy0IvuCr{xI*$$H{u;CjG2xsy zlY@!$on-0zRknSx*$qJ4V~Lp@(5dWo&QSOfQ%_7F)HHS-G>#nJRN;Ds7 zcQF&D-QVw99PG?|TELJA(>6uQ>BT!^NaRkxgmE{*MidGM#rWhVfOUNugX=Div{cF`iD zN*e2ycClZCDECx#mDXn~+X}uvh10%#O(Cw>)n|XzMs43FQXZu)EdC)ny5@K3iLddJ zSnF^S(ZSGon^9l8`A1|?nb|h!FT?kVpQ#NXda?M!bB1VSWN$3zNoaA?<}|XUF1ydp(`_{PI zE$c*CUXqr_0z15HfQI1hgKB*`9rDE2layQ(JQb4|okXO;gjO#)l;MmL=xp}O6u{0xC<3W zT7^sC57YRagEg+G1K0Q8BKEn9Y{FQ2uxm3##MTA_b065Gv=r(1AX8%+s=uDM?+8vX zdhyUqWO@R-28bJL%|g8aHl?9%Awv!4M1nn}40Bex)wS#Bq$`K@u@=)Jt)Ku6PX4r( zD8aTTP2(HEY`)gg>YDUWBjxy+wa;06iUK73ukmiCwTM&wYwoy<<>sPAmg&M?NS7b- zozV@+@9f>SovOXl5~+SME><+<21x1QHojdn;WGkktG;_?B={_OwhQQ;>1J4@I+S#W zl9Z5myHfJ%`OtJ*-Fjy5wmu~K0nPXwx)B_rm|-3kmw$J>@+Hx3amnIaD>l@EcE=O% zoDx+K(A*0SG%6CMP7Ox&-AUC9P(}7dja4?kS-PAc|J|&4y5p_WFKHaI&1jQ;1Y)dZ z&!Ypt-U`LFV5hEB3Wi0lb-S!2`x+0keo4CUuF?SPg_$?`m)hK$S79|mr*{U5D^zhW} zw(NTCBYT8(skLwNIQNXzYtnAKLe&}Kiz*T^>Ya^ZYOT2!zCoz7v17@c{Fu`MpO3Yo z5eGc{n}heO_8ge*s*%Av-4pd$8QxMG0+ zT9dctBuixD6Lcai`1IEAXUNI^>}?4+GOUM_C9vvGG6YDoVL`LEsJ?p;Hkj77AR#YU z?|$p|oMGoW*wm=Gj+!d1cdtk09ye0OmX|gsP6ofX@nX>6;Nk$9Fl#1rNv0#- zc_$$7kc5I7o=A&VBCJ>x-a$S)t^rgpUv04&1*{`&W4{%yU#K-)&t#%B+2kYeNl-=w zyyKs*(GG+tiZh_u6u7-sxNwH2RP_B7p}NtVMX!YmAEYx2;*AH&rBMA`U=-SMYY^kl z==Wumx5RAmwe8E?xxJ-@my1>sMy%gQr^VZats7D_(CBjt5m|*Z9LaG~L#uhIcWw%! zN7r~c!hx19QVSH#xWqW}?wJDUOkD@B@wuM!1T7cMGnECML+pn<(}q2OoKxir3`bd~ zRwR{v5WCGx7H?wJPBgg$dQf^~lw;7WvG`pBjRGdNzl%#zA2^V%2c(YpF>b)Rg&skC z2J&LpoQZ8Rbg!r8JARC!;*r)k@A~YIgfcVBlrmjK-kI9xXOowKyT_d$THdF|5^vZp zjFyD5U+=byD(A*TmSD@OZ&uJWAM0FyJ-z3%D^{K+A{&H>?#}4bR^D?M&aivc;E%-1 zjc@&iO2lSNcGu>)2IlJdnETh$k$v27zByxl{HV1`GJ6m*@ufEG=vDW#$$i*S*3QIM6El5&ufaYQ&sfWZ!6DG3*8 RFj!tv8Ubz|V*`2M{0}H`iOc{1 diff --git a/inst/doc/vignette2018.pdf b/inst/doc/vignette2018.pdf index 504ca76714d90f53915572e31201634f0a63f68b..e8a8041281cbb76d20a0552050966b98d14299de 100644 GIT binary patch delta 18299 zcma%?Q*&3=dx!1rj^v=Ztd418$Z;}?Aq=1RKqrbY>D_dID z1m**9_C{F`zBqoqEQZ*&fUH|@dYZ81=gB}|IHOgcEE>wHI{ICvcmlh#43S&>)C1>fwN^45(B&Z0gF=E{c$k0sR!rsZ=g zVn}Cv;KP-&F=^gn?XFW!rDk*1lCyIu+UeLIOs%;m~)lNa8s zS&6R!CBOyxI3@9?%JDRb&&Iq8j}kEZW4nhL!pr6{FjTMx&(PSzu~3(t&(9u&7CfQE zhM^Tvs2CU2xb{!7)jBLh9a+$5Sd%Eqm-KLr=}M>CX%yw^Y2w21E8(GS z=%MJrw%to9aUG|L~I zJk&cLG*6l{5<4|ov4FcE@zQbroL?AKrSCC&6P~qBh+v6V1cyk*tP|*SPF>ZUVwe%% z?!Jd1Z*3^6A8AaLmO1^pW;a5VQ0`xANhWbkg|KfZSlEO1KkWy+ix{Z3W) zjwOKOFAm&#MDWfPe&nYg_a6Tn4mc$79~ko;q&UhnaP~YNnxOiPQh1M^yPP{`$MHXL z;TOnTWy=F(+|&GI01gb^IezvCl0*_5zr%xMAIgdb#;`_8tjq`i7FKVNwC3J$4hyAP zUm1AVl&>$08X1_`xaQ#|R8w073r?Z^BwIktCe zqO-vAm-?$u>Kb=$PKF?_K+{$7l*@OU!gI@IY3nLCKBv)V1raXDH&LLt zOOfNUekIEC(knl>iqRt^nsy)Ey@#W{bK$nkHL$FR4T#OWa*MIKg+-{4tTNe! zgQYu>ZH3NVA$o*#U0$%Az6E}q#Z2Sg?!FLz&lB%Mbl)hX3zD*n-8rRdFL*wWhWFwfLb>p8M)>RGsRR#LSd3=sb9yr$?@{(zZ zM;p3$pnS2W?l`KPEoTVx0)ge>#X%~?sb>0h;LD+B!;l^)xg zgE0eW4Y#dQI0y4zRKy2M#O2=oBD_gR$Z#KX?=?BGbZ(07Y`iv4fH&8WNP0xI-mW-w z+-irpeRmx~yshye$v8yejBO&k+IiPZb_Xx+$MJJ~1>oxi!>5@RL3lPOeL!-Ie>LVS zh!+jG5}n2vdf0wSOiHXR<&3)`9p31k_QnG$S1)b4h)dopO8iLXUB$$YDx+)1=PTn8 zdfm}YwM;8Y{?hL3?t@^iOx~G#Yb(g~XC;s0mn>vW{+k#TnOykAsV(R?*>r|MBQuJT zxWJRv0GG>djX%*CdMJZfN!EU|FK%{D2{7q|raGdgnE|1)sigeoBPl7YD>M+H($WI7 zpNeNi>OP`V6;+aZv-UA#)D2Guu2y2+*H_xtJ?7Wa`Jk=UJfZZp)9WX=zaP#Os~GajH7By+wVQY%j2j zoZpIscG-!vaNbyM^{3A=7m?outO(2c`+Q#=S%fpje=?+ks;!_T`sZ)jF1Z18Wje&g z$@}6CTJk_%$iMVx*b%Pm@k)@hG#g^%gY~VnzM@@l@@jtR3Kzw6Gmj)iOBl@qgF2H_ zC#>?niq=`YxyJ?DIg2JQFVKf9&-FBe>Vt@gHD%7E_S@pY+g!=RqAdvd%aB-iSzS*_ z#C`^nVo|Ni24VEfv6M}>Ndp0W~F(Gi^-oqW0Xi>5KiJT_0|EF3@gPwC!07VV)11#_QO zVNj5NXs~OW%u?C5VSerLLf`F*>f8Em__(}f@Q)({SSt^-fiG-%3`hWkHYGxr>~(;y zJOyMUQAv!J`P`V+-s>Gj?uSkpS~6VN>1L^$2IgcDH-3w-@)!iIllE~Tdtv6<5a5@B z#iTa*F2v{B4bm~zLMLLq89t+}a zt#ol^%3t9{;e|B^yDtE^RR<7Y4&rDjk<%`Oj6x=i!-2sTPst8V#4hJFVAg&S*IF)< z#GO_@wH^(4U3+l26}k2Ibf{_1z%YR0E{(K0=2ROR#Sb{D3PhZ+wn=ctDaMU??U);5 zW3uZjv24X0$l|;UF2dV(YDx{1_uy);rdASVBI;@?aj;i zp3LGj7HlOgwrr!hI7N~qyIyzi>EK`(3f{J>+eMYumXoE++g^eRosb?FP}@Gi zT^q$?t1kR;H#6! zyL!Jjop&R~0k|jhYw(#Fd$ZeLf%QIa05`kW2nhxx_{=|L)%F%VTs&F@6Rim9+zlz$ zjVwK?F9O`7K48RWJdK?cFk4>(#5Fdzik=Rv;%cKv$e12y4J?VwaM{zj7Z31M(kf|) zROE+M+Fe5&sqT{VS}UGsaE>>c2bQ`zIl~5U{I7GUwA^pZ)9O=faJ$){ow3;4Jv`%c zy`4LlGSH^hWa3{lp16-BRPCbkwp|;Kb6JBP3jsl!t@7%#LaHeWIW=o3JIq4OE@UNY zWI;OnrgCH?C3Tb*4ME9p_u#Z-E+S&#Kq(=0f^O>~zy=3K^~%z!!*O)#;=ony zu+`o_$>)YACMW$b$D9>xpP2bEVO;acP-XQfxo~*C%hGD2Q;>r&`1-~)Cjsz&nqU+{ z0RZ*ME$O>t)t&3^hzJidI=JmpeEHnsv(o!L5f_0AFF32f(cf zp9=Ko*hG^MTu(CM15!*Ha`+ilu4__oP0&36XPVO!_T$V5dWBdG>P0|!Qab~;R!x9b zaA84WbL41YlpMnp5`WIa5L`{yYG+__Dsu7bH8ADYM9iDL|MkKv-0=)Pqgu0~o|C+w zP_pvMsNPyN1;oJX!?UQIdf41?8_;_tIM1%BMT>T`T)PN!nj=q234cc7u1iu5y~*s2 zS*@*2#r)|cL^Jr<>$M}5h0y#AX9z}bY60FHR>0suQz18)(EL4QXtpo92=Jaa7Gr{J zzv8d&1lJfahUG!Otu17T*wUl_P^EUzY$$=Rd;SHbWgpkH z2#@X^GHC{+dN6I#u@QT{1C$#w#`MUe+#r88GH-!!{oYzQbMyL?toSL`(MAc4xKjof z3s6%P^8SP;vrt@g`*7t756EeWr@kIt6dSlun}uWM|pPV56mTSHL3rbO7x(%&=j<7cW%1JVl0J^9*<45bjhdlfR#V?7@ZHhq- z119C*WM^VzHesgYVq@o^V`4IAqhtI}=uA0`O}IGN%$QhAjgmM}>;5-I7&IYx7CJ^2 zRVEg07A9^Mwj@n7`v31G2`%z}GbV&CLH@sdkAflT9);w8vf1>4&I$s;&Y5)Jj*%b( zpQM0=1kiKX;zak`GZ+}tLnVnkKEYtb%$K`#cTbeh0<$Xl$hB14)uuntwu#n}&o z!{}DzcO-(!@d*`1&1@G*fxTh%R+2uqTu`(!0;uluM>CE}ph{v+eN)mzANePvkX!m6 zhelM$2@*YCDJj}ZQ3{gYwC95*9D)FCshWVjtBM2#u~2FnJ*ggrs8z&B5E*I$S+D_6 z_nNFY-z*{mHYkng5H~E1mVjcA*h|?-b#tw&9l2LH4`eGJ2DKm@f*fN9vM6=;*^@<* z0o!58_tYzvsBb!^A#16yq`VEo<1h~IZiM%_Jfv|TKliW!!Agpvz(?jI?|+z-tH5#{ zg^Kp`VQobB3n(S&Fx3NFh&9zS-H6*&kHnxV$p6MKn-=|*8luNg(xIn37g(`Gva4hw z{!j(A_Jq>mqG99@@IU~ciQ|RYE`X)<1f1^_oQ;#n;TmI?i((C7`VnjQN9{5;{Fg@3 zwNmr6=?g$lGX5?+cFDq|Y%IGcKi9R_rl>UsynhHE@6v z`lg`F%f84Ep9o*W4D`{9o}xx>MjR)s(tAk&7Puz2MnZZ?eLFA6{r-$C;x!-E&hfQPcSRy6hD0d+u6BLW89 zUOm62clS`S!av*u`myqPv2zTttjXa;8@8^-s4oY$0N$KP+D#N*osz~Ct)BB6p?@CPWYAGa|Q8PWbDWlp2!@@CCOfUWszKksCz z4%VtbOm6MsKzdC=YVPb&*C;T* ziTZ)I%k^fruY;2+4<7>Lf=g<17yTez>r~5Nb3c`)zs#Lh>OSoBfD;3lh7X-o`wuS> zoV%VQec@-Fj>ZqQJr?)(G*Cj70Q?)$1{t_bgu7 zY3>|dJ@}K2@^tj31vrS%mFQwSZiu&5(Z*6EEHO4pCyfqz{l82Ohu3V|b zNV9M}MK*8t2z=Nyo|T%{`76M8|BK=1dJZhd>iF-`)s&3Ak7~BFP9qUZ{oJt#I-Y`M z7S^%Jy|$Jh;P$?P9UMXYNC2Zek(K}<6_~a-UA&8ySdCfmWPc$BT8q(@K-FuoMcC{z z-`O=T8eeoojZJ8P&=enCod*0`-x5^oKs`&AWX%^>5u5mdF-b;|eh!!UAx9DKVN-Hq zPYj-6Cmp=Rt3#B5d=RNlXjp9HK=cxM~&|cKvwanEyu4X$KOMArs7uvb?!Yo z7AR?uR2wyk`wj=@P%*87UO+8D{q)xqPIY8TW;Szr1JdFZzT6O;8Kzm-40>kzRai_E z>%T~y>4N3pF)(h@UJ|n2@>$*_&-r@ZOL>4&<53;Odj(EXLQc&?mX|K5F03eGB|#`K znX@+o6f?HXy3s|E*mv$H*`g`ywkmm5PO2~9R|kO@Y3+#PcK4gOPIXk&1w$nmW!ayG z!zN9`YRX?Trc4wdA{n>q7qlBhoAtM?M%;ZCoU6#Hw2FGwzg_88rEIRcbc)8}hbK!l zTte1=^(;5c z>>k>zh4~ArJX*L^%qVoybA0?J>{6Wo+(XeT@&;i_i?s0dL1U$Nv9s^ONk6#*&8K6x+=D`hn#jRt$r}G{!k=@ zf2V0#PH#t+50qj!9NZliTTdae$;K$AM@F-0E2LYd(V- zoep62(N&EKhzC)YBqzvLx}5~WI2lVsBhn0$O)bZr?-za!s5E=Ccl{t-GRO}-HeAAy!TEr^AtCw?L*4s;{_S9%yUoTzh`@THp97aGO z^Am1@tKIO9MD+{|Y*CB_C!wU&-6-1gu$dV9R4r;yrs||@mv6PUwidc=i(R#V(qFiD zoW0&%O?2w`Ck0-A#)}lo9PJ12)C*C%@L>>sH_5uPw5etL)##23h_)-|Kcu&Qr%*1b#odhRM9uaTn{l-*B2< zqP4tGqaaj)412v(6pLC5BTj3%AKQ5Uet|q7y)c^+Nta0hT!j4mFfOjnW=3`}p4m6L zGHEBRPP-nDB;FHQEHnGx^Y*A@tteX&HC+FM7lT-KTr7%jdyD)vu^qZ+=H$^}z+ebB z98NgTjxC807@z_|>w3K#0ler%yy$NRKfhg3MS!28q_NK{5~pr~j|a(8gSe;PT|HdD zJA}BW8}1n3JekM;>E~~5@|g0~LYL!pH9<=2m;-%mFV^%eDK*QeVYtAp{II2Q^djQY z2|{Y}E|oy`VTq$X=k+mVagYYz#Y(UhS#vyv#sSAxYDGE$P3%#j!qt}UujY%Xeocvgz5B2 z)czOJk2rJ@=me!k=r3V1DJ==muMwv%DKOM}1UTNGA`I;zp=_oO-P`2w2cXR?-U-L{ zVH$X-RH;mEPpqu)+JoAp-N8>k*{+3v7b8T#oY;*}*&ySsydc{c;fV=DO@QC%Z_U_k z<^b<`BL3A~Bok0!FaWiCmC-Ppzz#G}sA?zcE)SN9>R|aBS=l;lnW{Gg+}(|xirlg$ z(I}m>e-?QTxX8+wg91!(x64%{fE+OoS_PPr+em5T51%s`b#se8QfMWR3upf#^dvK& zU)~aVaSh#zV)GKo>J&rZ_ZKo(Es?}*{^NZdR5E6?Ha=r^{8P?EEjf542iG3i3df4L z0D-pM9SC1E`p>|{fnOmt4sadVB3uT`U_NmCA}I_l%KC|HhS8g90W2a29=9Vw7=in0 zVFqyf?cCo#MuI0pR&xj5D&N64*}{_mF#ju1+~{wf{^g!D#Nf4$lc^wiZG=RjEDRt$ zAQB9X^I6mRZ+asjxN>rJ7iF#@u1kJKez5y~EKt#oz;5z*%zFR3C>_ss;%qO&{=Q$5Pqw!F?pQ0r(hzGk+Vj z-5|gP$B1MvP0wZ2Omun~QlA|JF_`AxRPgg1R7Nhz#lU#?)Lbh<@f``#lA6FL8%tAO zViQ~~PB&|O$|&Pt`0AWKY{2-L8c4`HOUH~-VbfpER`TJ4s0|{Gp|usNAIZJlYEN#nXiW>$=B-XqP zSp@y#3(-zkz~n$k`wK$7G{`aN^<$fa!_P(50X$D0VgyJu8$xJ66$oQp zhw~GzHoC7$I(BZA1lsVP7+(jv@`8orSAy*LB0Cv9R2&WiZ8%so5q%#x1kA@33iFr= z2Y9;Ds2Zm;Mw0@)Bx7u{LjiFjnWL&O0>UHMgqSvLZkCV@aLDC=tlz|8-0{f(1f8ds z76LRo3rayDCiK^dIbS>xc4MOwG^1ACJ@|fc@saw22oBwn2~-3Cg%k#Omz8k{Dlm@M z-vy?uUzMSN+e-h=(>rV07G?d&kP`uxAG|w4n5z{g>X60papp!P6Moj!Q%B<#Jex1Z z!|S3vgJ%dgc}kl64*SIBS~B;#Uo^b@0f9WZCl73@V;z?hJTm#+D7~GoYfrNiVSk>l zs?JsabyP@Tfp{Nq>~*T*xtDL*P zrMH5)Tx6K^BZq|vAzYUCXH^Ju<$TY}&2w%GFL(z`q*@lR%5Y>wmWk}P0gs{=)Wv$= zjJbxGn_6C(o#VI8Nc9Js7#|b&au_$p#@85zRblG;mbaqI1W3jvJ z#omiY#PzQZN^68OmQLAol|)5oYMMh8xj~z++>u!o$!|PzgpkzM2=q%nB{uqGmsi|A z=x!A$Qac}DmhfHkkDGRFsT^fgI|RKcAc~Xt)m_~jEMfV+H4SyKg!j9P{I$AAG=X#~ z>dyjFg`ag@8Zv>rkm^ymr|N2o}5muW2k6pMh1Hi{xqOKe)K(SkFmb z>IvPTj!*AL84Gw2`|wK!YkI<%rXkXi17EW4>lJ^AvLRkSFG?ZXqANGEa;sTtE-L!B zBOh^~fxZX6-@#Sl0YimU^5%aj=oGCYg{+-1W3~!guu|yuV=~k;*aUKHz+EQUrC50= z3FrF&*R--kdR{^{g@+P+3D<@7BJxcG^GlumuV}(;F(?g}q`K8fOVTO9L}%4SwV$?v zCwq^BdXB&a#-b#<)P7kTAy^C0YGBPk7N^Wr{Xi39A(fXR$}3APbWUsWrldhGH#o%A z1N&6c6(D#L_7nj}s!;nB6l_DDLP+LIjYB^G6Zj(RIF?`jlQyUG^roV48pS-bQFa=! zq(3IoN)W%IH2vuZ_qnwPX>xHqccG$ zpaD#cX*BYO`}b)25Q34i_f#D#OT1t}#aHnVDQ`SsqO`u}$wYa#c5bV6PbLi;!%I>D zOr=b3hS!Byx`R`*g{bc?Y@W$d2MyeG)!I?Y$rgt;jHC85*KkjKH=5(QJ5s6z4|k2R zI@4Ad0!NXuu&WBZE=Mr4%u6ySD}xpQ?hFE(hqii9ERyKs)G(l@f}Xg(3IB;;fS5Z4 z>j`OMF;guYJ6A`bsUILU8!|S7!^Mv%Z-EQv-Ye|`8%`~X5k7e-e51LfvOBP3Mi_XF zT`4>57(Gxb+Rkim-c{H=IgSHf!qeuj`uDd^ErMJP`b@`i_weX8u08cCEien9Mb%24 z?L=jhM4bWqMGTh;QWg$1V_^+5@=|%`1gGlrF^1?f;Sg=zlzz`r)}MX<+_iXiy1OCd z%1eZ<_H`O9lb^{DkZ%e7!*V{o?mU@U*}z1lq{?3Wi^E41h53gMr)8NuPq^3o zrHOlLXgfr!v-Q7a4!6k!Qvy`8{%m!kujFXDy0;P*D@Iv3*r|}g`jt8^Rwm7r*-(`{ z3d8QT^f-Sw<6H~c5H;19W_R&h6(c5+ERhr@xW7Ur=-Y?Il*Ci%CZ&?<(DNZKii#7 z2e7@o&VE3~+3B&kMY=NScC@k})&LP*Iw@5du1u@7^?_3*6jNd=9?MTgSeI+MxZx#F zGnF;@uTEGCg4)HxJLLyJ5%?V{Xn-k%r}5|MaEsXh?E;N2g1wzByP0IR3Wddh&12xt zWYICR133(S!XK>GQ3S}~&&PVv->+KR$bUJJ8Ctelo+w<9PI0I8zp_F9>TTd+Pv!JH zCKrp}j1*@7$uq6d8iU@VqR^E-Uc(K921DCdM>fj*VQ^uZsI5H)7$c8MT07!x>B5=E z8Q)I{LfmexuiD>Ax6NvmyOH@DtT-RJZYE1a_-)#kACd%#bHi!_wVv&n)5#@7^&{IM zmZK?L&zrZy@Ys6xPxI9=mMjmFT3alXbISxsSNdXZlWgvG&`My>-R*~%pNu&;og zeE1vA!IF}3GCJ*{sK?`GRo$eysppDr=r37y!_06yXnPnt24BtGs8CH84G|Y(aujwt z6bIV6DAvTwp?P<9$ySe4(IZ(Bv~qRDJ*K+%L4d_nqu$^EkhE87_A}&?sOaPsG8*QG z%ov#5ebbfG%zTbIRZTbkn#jC6$#5St!DvZcmpyf zPbM+JVSw{+h>LlIT4JKQ11f!a23AyoZJ8CKJX0AluGAsX~TKPbfU=y8WIx3;~#q zbLYwk`ElBdi}T66K_wl-AqH51fa~wG$eOZHxABL_T!!VXwm@Dpt8B?;ivP^lTg(snuFdRV zSkmOKLCLxtCxgx$qnc3zIG0)bMNIlP9oxE?&uP#FzMRt^Tk*ciFNWa)L-8Bo_aPRz z00=Z#o6fxSRuq}Nzp^n3%6NJE7V1i8?dmH4^y9jUxOJbr?xbz^=$jGA?tZHE6xS9S zxe=EnTuzX_k)zOf;#p`hN0|CXezGv2F-D4!r|dKWj}n=8=RDk==eW*0Q}D&!lu`BA z#gghtR_c=?BLl+tZ7Dl(p0876o&}To4&7L}#UK!Lnvxdo5j($H>^0x%iokOlV(o7M-Kgej92ZOf>;k5x|6_(jXkaS!Pefjp9hAk4AfG2*> z&Cm<#JO2*dlKa*&($YHZk6lpRqL?5?^aysp^7@2e0|vo{f6*0o0YewE|4r^5ou0D% zL)oSJPs7%ZPedf&X+Ac0%|I))llRRN6^*m zF4~O+gLSF}?3r`_ngK*d=JlK4#?69P8a1XF=nZPtmn!$qg-a3GisegGY;d&9jxF_l z)BR3VCE&8sZ)8&>9E|L}Bu4ws)+SrWHwU;6^%+oi0tuhf*{phDjyJuRw z1C%B<6zoW_RF2ts5-n~?b`Wy>$OHm`P0A9Jb;R%9x?&#>I1w}}wy%LKKSriD8sd^R zvKoMC=Z60+w;;{7z1!F|C%aSv+H>>Z?B1*653|XG*SPD^yo78Y3fw*ttLV;iIU@-4 z+N|_+JFXPU5BJ^(AJ9sa8yK$dTN{manCyfK&}8ksOvXW(6`!025i*Zglc!|%Cj%f= z@}+}K=;!B_lb*_1n{7yi%E+TDRC_XOJ~AMRQ89C^tY?bmNMt;vZ~W;2lZz?47zcOW z7cHZW%?Rx13U<^u4m=my0OJxn2a=a%Kh(xW&<(6&p4U}JmxAI52tx<-A%^gUY@q1l z`y<2n?XAdrmc)+0E|fcr9IcwME?6ktVf~& z@%QAL*cC(Ugw%iF?vIe2G{c7sEUwHqz^yl52fhK*b{s@~sBN8nb+NB)TA!N(~j2SP21#GF}2^6LLQDMVg3sATC9U?e(V40no}^m-GI1qXfgeE6CGw0(O5eAvZoV0AyD)KpP* zpgxg-ErpsY?dDKC`5$CS!=y=ZtjFG3jDHdw>FtPW+uIBD5@X)j9Msp~w00bafe#pE zM3E_rYXWW#S8~nSw3uc+lhb>rAzI;-TR-g|OuQ!!Et9^;bbp0`p91G>e;Fzt;dlo8Fxa3Pk7Sbj% zS4dnW%lLbC0#D1%IlnP$EH0|Q8(N|-S+e&2bmH{ix6N9f!~L%P^m%l@eOLwD-uoJL z4vjd#Quevu0DQ)FIC!=}dH!M#blIN)-V1wS@W9HLO6ddE*iG78y5FSuG;5AO-+GYB1h%#U`%j zD(_jgEwPQP91O(3a7P7snVjJHiw=R5Ch}uhQ{)j07@rL513*qKS{sbnX zzdX}^#Gz=@?~ZO$1)UZS29tLcejK%c7Kg37VFBAhF~nyjF;@I8vA{FH&KoGuygYXJ z*}N2~sFnWc1wvn_e)FD}&`Q?d9A9{h<=-!Ev!8=a(hPXSvC)dG*Ca7Yb*40=LGHgTj-xuhgjw zM@`F-lm2Gs#%4ho@n&uG#UxGBtp%tIbZGhuNB(%#I#72llwWXl;0g2i&YU|89)s^< z!;Z;wm^5Dvsp#wh_JNG?I;gO065|H!IWEv}5Yh}T?v}Y$R-!+Ee>x_0aZvOFqJjhY zEX&~+D44hdOqdf2@WsNyh#u$y^w#^&I~yl~&de43@|{rJ%JFW2B!Y*WQd1(1R>}@k zvBJO+B+nh?e9PrgMI*T{D80jyYdFVt&}q98+O19a+@Uab!g7J1pN**@J@eV4$HStP z?o*r**JF(PPTetpnJFr2E{moq|s`*T@Clt&RrdJGp6VrSa;=xcaQLMxRJbbHEJk^5$tZC3sg@&IH44LWp} zLZ=lgsiJHIQ6*H^aT2#evG(Q*83`PtF^(?cy6-vi;8YHpwn2-%z(qnZtCxv9?Uu9+z zZ^?}^)EUw{s=FMmShf=3bwEm=nvT+Wjr^vU3OQ6_G(dK2gsiRGSxmnQ72DgChJ zFk!3{jf%7eALg*SZ(tb~6OqzOx_@;!h(%xA%?M$uRT&lzVX3lU_SLOrZ@EaH4*~*8 zvwiO1ly!FYsr)fXf?Htr5x%bb$zRmlhwzS{z7TByi8et-AzH6KR0ylryFxU{?tUHQmerk!x zHTKyA5HjjrzKP;w*)h(cgw~Rbxz($rWII7OZ}y7Muro53=25|ad1P+KHE@&KZ7>2J z0J|u_AI`vNX`K8XL$^V!&8D_#dPIxF_NjB*6C^iz!glBB$Jx=mG1B|cc>b?1OV2k) zOFkUHZcf4G&;~HW;`8?Sm4W1)Nf*KRJHCFJq)TLdrr8KX=|(0#?P32;|<@r{fPIfzt;- z9V>+}OoS6mm@ddIHzzxKSU5O`BU$r9`PYjDcVgv89HmGkVvZQIdqUmbc8#Pu8+YBg zS&gZM1987W&k>DB@ecJ#*5N*nX*c*2>Wg1JSLubvl4R^z*H~8tQ-84K%MPV=2CxC)P z{~r#$Y$`+UcxM}7h~lq4?k}pYDbU%P_0>;{h_Dg_zEfAwlY25+C`wOASWp+>`i!j{ z>EK?jDM~SlMTZ6raTS?8b_9LWsY6<*80($a@3fo)Lkn0%PQzDjHlKe}g>x!#qELh6D_fsuiP-&=k1F{_fBr|mKv2o(=#XgU< z)u;B~R8M?O6+~NQHO@dGp#x!Kdp^@!)c)077;i-ZF=UptHymczc56Q1E;K-+FXPs* zGh^1Omc=(Ilc4i_vQ}UWTi%HLRUxeTcxTb4x4j*`E>;}RMqby!4^P!8xkKR3#^lrk z^OO!&fas$DXCRbCO*K%jvA%taX0N1lvv3PB+RA~;flJk+6);{X&ibtIRdXm7X{{1* zJ0nr`_k;8P?VPnJ*xd}EYu|11;>USToXmEmcc6VcO9-*Bbyp9b!1msYYBImbX{}Pv zP=w|cthZ=S<V=bp4&+3}u_^kWf8iOu+xM7ENR6&nsK-XYWrAwumM>eP zt;iw+2(k}c9uT?WNRPO}k=SEuPEI6DZ-P3sFdSLsj=AJj+^{M@M@=$js>Tu2!-0_E zW}8j)7~2=Y<5jeBfs+iPI@|c0ZDW)G8s(M7J3>v}GmW=fLawlpUQdL(x4`^WI84X- zSkAPa1gR!Ld0G+}sSOs{gTOPY6+Q`Te2YRRXIobRR ziAJb3c3kkW%=xzYJatw*VAy8--QT}E^uWt4yXV&=`hyVob|iwAtgxn#MyD?H z{t3Tp^^*mUD6XattEL~9V-P_UC#qA~VgRl!dD68W&x}@(5~dwDwaZHdsxs21oojrmooq9`g5XjPjNvWJ8D=Sst#x&dHfghzwX!I8wR0n!}L5R`B-B1=1U zIM`yzB^aH7_JT!!O;{}*p`Dym#DD*i21Rg0kcqM)$#e@wlj{?^4nrA}GKrf8 z-xEaxj-iy zje$j@(pa_SBQ-n^jtLt}-W@|EoSr!da`-)Jz148Tp8&1@;_|^%SdJK?vHS0jzm~25 zsNQ%d4RKp{fX*e5QfGOM&K@vk%#N!3)r0k=7cQ1Fw(daR8|ob zU4MQLbUgd&&ACrl6BxWfTbzzps7Ws2f~>)-YWV8VPFp!8rnVL%n-O%w`@hL#G_Zs% zbpLg<*EIavu9b#y1w>P+T@s)$<~`1Umua%9t)G$uX*)yqjF^)WL!eTiXD(iZVi&=_ zDQk9~E6yu~`sL-;8ZKzJe&3}Hx9{WA)k10h60#!WiQDJMduz_2X3n@D_yn{89lC6- zHqgm{3Xm$f5U$F4eatg47u@5u9m@J#;inJ{E!Kr7ibBko;peU7iY0j2RlRM1u*qfS z$MAd4S*_!pjrt-pAnqa& zoGN+(VMlPUvU+>OSZF&6+G3yL%2@#i7R-N0zhl_0gRF7`uC4z&`y6Vh7L>QKK;HJV zgnGfV5cg;cGA4W%V{k)=uIL=F5s?5p^n91QKaYp~e*&ilSou$Urg9G$X7G3l4Y82Y zyxE<3DE~c?-B#G=^}1@?0GJRO@w+0NxQXrRi{kg-|4$RRHWFdSqjn^6$gnz<( zX5Itk*$YS7JM3KfnwFHqz&N|AZ}6uzE9m7aU7ie1)Ns+q*va9W*C%Y)-+=;6D;~C4w{(bm8C6)1ythu)lW+H+#M>15UXZZ)n_z>g zZkyX>ul8Piqhl@}8vF`=5x?}0$;cN4G$X&t%3?4SnGE>`UK1Ig9uG|1+I{Y|6BGsE4*8%yGEgHy+LfT-2`fJ1E#n-dJEXUjDLhjV9Cs1<54)^t0{;=!;c`s zS)ng4(HZr6BcB<%Se*n+(yJb{h_};bJG&Ewd}i~&Q>efb&1uYF*rK%dpUvA@kQ|-2 zl`)~OB$)pp-TIvW*T6+?5eL(bMn5JId7eVyLLb|xbty76+KhB}wtu_D)>zMB4#$Zb z`C(;=xRDM2P1g?b#46~5eb4_YrPhPI;6wKx+UwgsaN@$=X!a-Zdo0I9gEot1^ZXBf zaaXmD-f<4{P!HPKLkTOIEYxE)cS0Lvg!zG;pkrZ5&T0aJE_?~gpV;3HKh z+=Jtt-md+hSiXQR27j-|*Vru(ZOcYv!5Qd*kHt;6R{KHS#rIkcLmPCd?X&hQhj|bM zi(x@QTs9sOlbs5pE_rYneftQBmOS>8tP;=6YDyl-Cx; zg(tV@YMkzx=6_aIwWg#WQQM_QUr|>VISNI6j=RC>Y~awcF&>RgrPOMbN~gBLU2pSnqTV>(i57|$if~Zq52frm z?ZJ)mY)xi{UDI4swZmG$m6e&n42%Wx^h-0}TDO)jHBXp~P|E17dTV9RHAFu<1s;?4 z9*g52W1G}Yr94ZmOmE0`HPt)CJG4$Aq6R5;?Q}|TQWfnT);*jX<`e4h8|tX^^Xsnl z8$l5N0jL=NikC6R0TYLQ4g!aM4gW3 zVq`vDAU-|{b98cLVQmU{obA#-FGEon2JmwZdM`z<+kdw$)kCzGAhC3jh`5Qxpb63#h*2UI1}pb@zJR#%+u@g!b8p}GJ?B0~;b*i` zG5VzpHnb@I!yo1D*MhaRj1;{sN{YlVDJoTJL4KZ#k{>Zle#$=b(>Rd)h+*0H&nAI@(0ONBoH4meQ+-U~t^fA_dg%g*K6!mhxr^FIOl{mXY+VQ2RgGb zw*-5)P|f=^Z*&X zUDbO&5MwzQqoE8O*oQ=uXMdyq{>3=*D?8v}u6j`UU%Wrw-25hwkz)BHUN>Qm-%!E% zQ`6^@IjcY}<;HvRa`phR9Cf)SHOIanpVw$!uA;^L=A@e!k7c6@G0~nL(O0=(CbgH; z@M}d`H)~w)LfS$N6Z|Z4=j5KS~Rpy8i zS*bx_mCD#1F(m)+CTF*x3b6sA1n+es3m*CqeGdV@2u}Z7vn=+?<}u zv%o2oJbtLm%X-&u|6b1+N7HyV!N#0;vzy5_>Y+OTF4lF~dFsblqqOh8dcVG%E@R=b z*ovMM#ExkeqZxmJkYFFTXu{ajgY|xv;2e1;38{P&dcvN^EaJ3Plh%GIQd}VQ5La3+ zE$@fm*`~z=1Ig)^bVzF+IAvR!>}J7D?#fk6DCH%pcp4{;#{?(!B zPW)KHb72Sda?BuX!0w=6^7%v3Y`LcUvM-Cl{gjbLIofMP@7cVe(7Uzc*Dc@v9DtWE zrx0d;d_Tomr?71w8mTaux$TG~BLX5~3d< z%A1={?u~G|JkEnn6Aeeg5paCt=r@Y=6GTa_#@xLma^GSXDsNIE(^(p#3^hd|~>)1Lh2GH2w0wrVQ`JKEqL7ScDF! zlNZfq6uoEXUp8;p;2Vg(CgynvukU?kB0hvZFi36@1Wz#p@zy_*hRI136kZMerY!w5 z0Lk{v8&LlJ>uXEh>&JAc;qCkf!q-&)?bl}pVRYz~5T;qQ5D}8#91Wv9=pEqxF)U=n zvTW)V!#R66|M79*q!vO}Xxab}>8%JT_ywJtdfi!lpeeveHwfj-4unV9qVK+Bd;8zQ zjC$(#T{chvn$lKA<-UGjcaKMBj=zo<#Urj`2=)#f~7Z`At6-HX>KC$A?*4TAE*nR zdFexjPez+Ev~O4*n^jKHmsd5YQYsXzly^nkm;A-0^ZyHX3fXt-1?#Pq!=$9rez#3b z$S992fy{;1Kz{KcmY*F3yP*M?}E9P|d{nMMK_N^MPGl0pLby7K`isovLI} z>T#3*mpoj$qz9vvCfn!ws+qyrDTR!!Sp|~_6({6Kk+W5PkDYfZ zR0tR7DeBVC;@hV2jD9X#M6_ZQ4rlc`m(P^jeDOMNI z@MWWF2~<)=P{iQsLl)n2RWG+`JvSG4q{0a^nE{NRlPs=f4umff7W? zb^JNqs<0Ydg#@WeP1Q*vO?OFumKhAKjf2fd=3S$>SK#WxOr{XEwOLfVaUybdS+z!A zuUv_DFF=C^m?rCM75qu?Afp&Y$nW#u_vJYt#H+LPef>2rE_-p6>T<5HT%NJOo=Bui zx5d3aVA_41wdU#<7S}-k-zY5~-5u^Vj}FIWTlG_-6lma~^8n`I9N(^kf`DTnCT>Zn zL*z6a%AXYndU@GL2YjzY{D~%RPpwUnlUODG^OqA66!CJJ zY6T6PEx+uwQWu@oVM`c&KYqLrM@8AnkSBR*YBeO`vb;2&rnWV^`iCZY4l1G6E{I0d4IbVs+)$91dVz(wqfY1GdlaRM5 z5lK%}EskJ^hrs5V0Cf(py(Q;5w|`4|BsR-^waBgm zj&!CtSXp2F5$pe9R&X6$>ohuRd0x?^Zj40RyWu~`P9Iw4Y`@t3b3{vJJX|7_BBKjq zLBL!oZIM%eI3XF$TU-`_9va0{*HCQ7l@z-kLJ!44L>tuMRRr%3qwgqXV=?20xaEl} zxGgJWBcZ)R$#W}_=EKrBXPr(YJEDRED3V*A$$oT}n7O#e8NSwdrNG=!xe~tS9ZoEs zc;C@py^_8rY`TU{d+`t?F?mZW90uwTXf}F9hQu5#gOQl`(R40zI`W2W!*DrH_9uIY zp+eu4l&fSU-WV=Aq_)^rmz6tGtmQk6!O+soRVbwhko4CE(1l_m*^F}pEj=t|aOYeon|Rxr z}8Nj7}M(X{a~XKYPhTGDMj*cj+b>ZiqOp$rj+uk%ht}_3NeM zGXr_vTuU1qw%kW2|E}7ZSG^f3J6(@B(`+apY~0ru0z}{Oc)BIu+4E=T_;8S6R6245 z5z=i`eoaK(EoF8{-op)ez8`^ppMk&tsajyGbZ$METFJz*5!z!{5wkt_T&6P>PH0-IFRU}W@LBd9~_Fg9PBV<7u3X)%_Yopq^iQs1DD-s}9 zLT~wyr1>!9rL`vkZMrL(%r;{fF}_0$y)f)uOMAq*D2l#r`h>JmKRrK15onJZ23LL} zUuSfw>BD8ZkmU`Pl865AF8^qeFLm4_oi9Vm+glxFLEO7By_TD7vgPUyt<^i)ua!xh zZW!zx4tZEfo{pQW(dsh@?wQY7#}g{DG5t6$OQw7-K}*2_uvXgKn_Q5Ks{AhPdh33t z!4>6#Xlv=vcEG7^>Auf=!BT_#C3s@6O3VV8jA%0n=y)M3DEy$RhKb={;88+SRf@<0 z>>63&wMO-gR`1KjVL&7xm^$ilncB_qZ}#jnL1p6HVYxSE=ljJPRjdYrs2Vpn=(@Bp zcIU%n+;HGPi^`scPkG6u#hJWZ56IehpYzeFsHgUN-;Az$r`G8*`)E#Ndhlcu<2>^S z>?12ah+b`*aB zm)%niV)&C7G|z4y*4fS1O*^f}uO>rnqnqww%k?zpx(bg#^v<~~c9P^uB`>?KlcO^M{kEg`8Fg*IVmxnlWl&@ASl?p{jKfH!Nf!Dy z+HNmUo&6|{hf7T|3{HD4D!U>?qryvCR&KSW4&4Ra1hHc@HmJKGgz}>TyHd2`3r%Hh z>e#9F%@`~?UT}$#IDwP1la%kqs}BNHaVUF*tYM&$a@vnl<3<8nY6k05oaYPWbDS_6 zl?@wZfy0te(MJhUFRsFTnOYtvkWl~58FBB9uK4NjRvGDl9do=2ERWNSB7%bsvrL>KZ(>-rsp(`OED5by+M|@g9z5u=gtc zA~9+d8RO+Mdmth%{@2_Sv*}QYHaQKXN6+&AW!0|XyP02egV{eHt{=Vq_=uCe=NYVH zhK5_N#7sjjZSZ>wnBZG(0-Lb->{L5WjbF4gQ)sE88r#uz4m}^ZbOghF<-{Wl>K6%0 zS!_#l9C|nyg@pwxi-ie?83Q-FDHj6^iv>@%Dws z^LoA`&14$G_)m<|i=dy;MKqU3GfoaUV>B*>E+5)`bjIw4tp&JF977UupeSi%O;E}m z|2HnPj)SnfAQGrT*ZDXa_R^rvg*jaIl!ZlO{L(OGG)I7zJQgk#8qO36rEREV0Et?T z`C|eBV>yGAX8!dtGCWQcjLdo*i82zC>)S>dG5JI}F*+vUT0w1t_oeet75i8*+BZEI zV-7$~$~el{NG6S*nSipU6&{J9=8Q!vD$W`aF1?%>l%pV^z@6;VUh+ojy`nYt);}s_=1xaZou*4KHOk9E;#Uld^tV=l4l*sQ8#GJV3Ev+WH{U- zz@|4Sp-!WpR%B>@>BD55zV^YcV`pq|YqWcrP=VT1{Scu~){WNXr)t3EK~;sqmV;TK zo6N9OlfiTeNT9^8IT1|57^Zak|*}DDMsTTm)@Y7feQlyk(WoePfutwr#e z_lN!Gv6WyBPSk+)z&BNDy7F$E<;h89mhqaq(NE!tqr+{&##nrPL;*>=D@qG+{0vTS7rk7%0%jgKO_nhZ6S>e}+zCeL**|!HR@SNM{|U!Z z5(1IB#9^VnX4e>b3>n@n+_^ZuoLIB)zN`HTz*q3j1EwSHP+momPon6LrM#(xCVzg= zkGqhSyD;`o^?#Bd2Pn@yB#+%4@B4XihLqAa3RK}_67x=c2lmBYotSx<3cOsq0ln$# zVT1v{yQ}1MrYDl$0z#pMiSwt!_luLmx0~KUL3giKoCI<61_=sh`8TXS@u#Ds)n?~a z=67TFth_O`m0 zWUb?>VC-bLo}-nyXz=R6Y4U1#772m_xVUuQ5EMmviL;=4|GtNo-{96!_V9A6b#a`g%U8>{3pI1Or4Th`5n%2K&Hcu7t9Vm|2uXnEHM3rFB4^*zoj1I=gvgT-(s^WzYX`#C=1{(9-L0(G@$5{Dt zZ0K<1tGGeycjq-7J@(HtXyd4>p{4(>8Gk?Dx{B4Gw$&(oQ;M#tuB|IYGq@E|XMDRY zT9pkgazHQS&yrJ&u?{0uumMc})9&AL8SET#)fyg~DdP22);-xNd)}!WK%b8!t znZJknQrWiz-i7zlNWQdQb~VH*B^VOeu60TqW0Xps>g~k}U2{lZ)J*02G*|i#ueJb= z2c=Wg6l#3xTS!(E>*{~y%P5*dWQk{w?T(87OlMYy;IUWCB)(! zVx?yhL?(i6oe9 z+F!d2UvUE)J!z4rrJOm_fc+@w>it_r1TBq#(*{%)x%QZ=QPY0)Q`*|KOEQ0V~~6&*#4)ofpjp z{4EiR}k)GtK#@iKVMb6wyjcpKJ!(Xmn2w07G9`?3>3a@6tG0}FcxK*AJqW86L zHYdS>)EMpf^KKY51_(1X2<%`SQQ^^ppteFHx#L6^4I8IaXO6_`ed8#(JngeY!&$mtq? z$C-2XoW486@NH2YLTT;XaT#28@|3`N;AkIt~W}zc{{vEroj;;wV)ASP)j%% zZKJ33v#w4SzJrc-{~;Xd(zhswO@4}}i!Z`S?7RXEV{Bt9A9jRQtPE}8jw7Hx@ChDJ z!$+gcPi6quL3&IuevYZ~fTbl3H9H7lIn*IWc3SLX(m~5HUoZ=r0z+aKN2$<>{>3?L zXXUcQ9Q+&m9h)BGK+$P1^WgYfWSVLkjK|pi@^#IVn4xZia!uk?VL%`QY&G2SL}SLF zX|~I+6gL9@gx_}@Lv=+{G^+}kv8ty)?2=MH<3GmBDa-h)K{~bH?pp2SReLVjeA^FY zScZLmVG*WA}5c<#E^)!?y_jxfAwpIcJwyi~43?au%t+11i4 z#F!d;?GYt}=n%Y@_X3Xs{vQNLZ^Y0Zq!+9xlY108;jF|>gZp|Su^}i6I$5G;?a7KTZ4^NId?I?qEm)E8W^(D4i*AhW1S-=$_gd*cKSKeUs9A@_ zuXr-vd}zy(*1X8Ylb$}?lxw~k&nkOTyF|o)_66)uF*jK2;2*eK>u@%DWV?&lHPK7; zO_edlKpF?}{idU{@1nLc-3~2Hrp=r_R=^&=dFFM zx~bO;T|@s9gfN&vN_jK7&Gjl?g;fjE_sw$k-1vgffAIk|CsM3X0PLiKg7B_xF6PGe z@LoB8^kh6K#-00yKhd}~-jSv{nR<}qb(0vh(?Tstq=|jL$)QLi`b0_$FOF*dvy(D2 zGndpb=g@eYynekM$hxf3BkO-LVfb*Pm-ME;>Hqj@PLTR~k)lrUY?nIp@_bnSDc4VY z_Lb|;|M`ZQcXP=d4=mDp`n`Qd?xl>YTrGAx{jMQNZ5uBqN$4Y*Ii;Xu8#9U!a!?qt zGKpP6*}g)OqwQ4<;2!^Ru+#8sLRlK;=*HQ#C1a7dMF#o-Q3f zO6*IbLU|kZx4^tmSwtbrBc7wG>>zv4rP<6FVMCyWWkYBh6s zTRq_C;JBR|g?z~(&)uP5%mrCmvoa^Krxou5FQ@|$mWZe)(06IVCIfH1Ob!R1;ZCPP zU|%_9?sVSdSiCbCdu$X22g{3aM;&qhre{hgkMz^kf78`@aA>hWDw5&Sd zN!P+ch|&^Ik86WHF?Ec2C}M5?FO=}{JssH4N9N8c&7!tAp)SVZ&FXPt^?~9lZC}mD z8(J~t4L+7|^OXkE5V5Eu#-q_WR#JeR5I~;j7fEI|QC>k)t1Hfcbv~88j+~|jvc5ky zg8Q4HOz<$sTON6d=|UQA(5c3pNudjUr#Xj*8a5^aMT=gbMg1I6q;;g;By&N6I=Yj8 zAO`)QB2yIO*Q?n6y}*T;l+hnrixTV(OyhOCco}%{;5I+I|E?U^NHuD;H7IP94G5|y zL+^w8B3R}Rrrvt0PTr3CP=FyCjF^H%IQ%V+f$GqnUb0e&LbPECs9VyXB~V`9z7OIP zfkS)GNPcU!!M%!iAYY4;us&{d`Fk(MYL0lR(f#4rO)jJ;IMup85dTJ)$553vPRjDE zAPRh(LF2+Qnq94@a_r906GMJikF} zQ4n@d4}An=c2ty-axCPhOEU2UDw6hQ1!P7IR|gn3?L~+B$AaYa1HR}XBoHMU$b~e+ zEV@sgpg4bknQ=*}Mm?L&m$yg7whj68!b}(el@IPwq!eE}RP=(FHAoC-WSYV{b6f>IqA@vVzMa@KJV>wmBgj3O;lHaO2B)=jOhd3qgT2R+3r5D;q)GNh>Dg!wXd$ zA~OZy*ItPq;)7|o)NcS+oU1eC^B{O4KU}Q6&EXF>^A<9 zzR!rn&oK5FQchV#XK6`-Ee9VO9ti{e&dV%yl;W4%5Tz&=RA!8~tG59A4Kx$3u+wQs zql^FCv&YcQCLYnjZk$b%-smpd?~|fq<^9@cb7w!j{?KFRetfM$#zJFl=`1Qh6CECQ z=b%IW5{N=41L=@}0Z20<6wOVqdSFeFJq?8Yc0`jp-cy*0B~8Nj@Rk3rKee3C(P{aI zGMYl$Ahe=9@K6L}aT6z}9PEK9gMzk zRob)otl)^+XHfO@!gvOj4fOAxA52u(RI#CE-}}hVh5KmI6qw7{MBp@(DK}FH+&5CZ zF36v0jy@&|%$_%#kNwO1%`xM%tqy&Nlw3^(@zkI2jK9A0(g@SGSt zcL6g9y00{(`CtTCEOPu?jf?G>(lE1=!5ya3=;(u;9ssNj zvu8AVmGIb4T{tstC}TXQ2g?s4#tCheN)N8zjJBAo+(6DU#C8t>O#AC1SD&>Y{AqR_ z#R%@TUs>%mEMEZ1IxAI4?cGgF!TO^XE!Avwca1XQ{KMgl$P%lw=!VJq6c5UJ(;2)5 zh)Z8TwYm@QTL_+ftDl&FGKa=u7sePW%cj&)8na?zkl$(ra>&|whVXP`QnwOe>@m~; zn+w$jI{<4UabW(fWdgD=$xvWE;ajAef5{-#=U(1}9_8Kbo@T^9&VKA{3I2RoTwDAQ zrZH{t5LGXpX1bHa4u7vC?Sorh`5K*S6)WrMGowfc`9CG3d^(MZuHn?mVPoUdJ`ODP z`hIZ)tV~*Bw==0AnrS&9TPK%bavQGuK;KN4IlD|W*SPQbXw>OU*(=*d6(*sT1v#QKrwr<+NOJbr@Ef|~ z#*@eRLE5qQ?!OkCf_p%cA5P0oq3vrGSNFeM-K{QLt(Hnp%{fmmom>ZNL^9(ub*nhb zyMQ~MS~L)D+#6YZq6q}EnWSNIUZumT+EAm9!4QiyYRUM(PnL8ib}C^4rx~}B_!<7D z&?{>ZhC0s~qOAavflpyHj1SulXPeTC8qYpsJqyK55C^$ltmQYXxVZfa<5QFLqOwF4 z-TMzjsr6;7h7%>rD3lIdrIjgL#9?S!;DJ%rqay zipcj6x(2aryeHpW?gnqBh=u2dnmen0NgINQr|+KvG?k^62(5@(oNt%8Z`O}bBOn`{ zuaUvk7>~}8fau%Us*}}UrA3NS>gt_s`5?lHkW(M$|}ymW7%XUHdzsKLb292ijx$}Ney|14|ylU@rpmOcwaHZo)^q<0}q~7 z%PVr9&xp3PFUdL<^ECaNDA=+1V-sn;{1< zh&n}n!s@Z(9-*_{EBjV*=}!IA8Y^zY(tPebViVvV@L#g>Zg z%qDj00IlR^*YXTo>_#xp=qPN}D39R%p}`egmRXmPX*W66>RdOuAzV15I>0nLeh^Xx zF#vCTz@K3IZniM?Dkbvr>gK&F6%?^-As~Cc%qfnAgNC`Jtpv_2@Wo1N4X+2OVSVOx z$!+q^IjKHjJPKy9XT4h^nRn%3w2CvS3tkvp-DT-r5dN#F!GAAUA!}*B6PMILJWs)y z@@lM66C0@+i;vJZ@Q8W030TBs@UVBud5@~4Ha4-zFN470kjcp!~2za)!*%ondTO6)A+((yQjRbdde9C=O3>OyiRWwig6_d&0;Y0m#_-> z1SHd;Ojar!``pMh{ zWTl2&kf77{1fTlCtFsl=}fHa%J!dZJ{uxLXk3rO*^4-(q0SC#MwP-4Hb0o!&Uys^;-md$i`qea#j}#h7{tfLEo5zV$C;&t9@;p zvg#hN)2<%}_PR21ByBS}6YNq1P{#D=_HP=-lu*Cwg&L9L?6)BzhVi@$x`h5b5VXjE zV{`RTGLeW*87|Hn+yWla9dbL+f7*2Y! z*5G0S1XkW9+9-pq(!B$}r+M>u3t`BxhPD27rpKE2EHZG-_9dZL-fUhcUz|j=Aw0^s z1Ne-(Pd^_CvWCjX!Wbi*vZr-)sW0HX3+5>iYOTjWx-0Pzv}k+oPaA?d?F>4iSim8X(3?ir zJ!eGN)4Gh-mtC(j*oMMuzwM{B&L6V8Xqm3j8m?YP2rgI{FK*v5(8~EASC=SxWvwGl zSzY41#4`XIw|&Uzp=P(RhfJGKaLzt&ZbzwGf~p5L%I;&`b$@njzl6u!$hz_RPUYum z08@kN3%rlqyW1bCIk{;hYdOH4ocFM9lgtbpwP2Fr6{ zLRPMKSre+m90~!Y40)h?1a+*eUr1!6rV%asdjad7tnb9{h2>vSiNA-lsvHb6u)d;O zWNdQp6;)*Om}R<<_WHo1)mcjO!~0D*?L6ErWYhMZ0TbJuW|slL*pe8uoCDUSfHtj2 z1`32IGJ4)#^9OAWAyWfb6b;(urNpmF&j{}jgmOFYbbtSH5!bXkm#!ELEO}7=#RHu2 ziwbD{;nddG9LrA325GM3Ej~GUW_Tqw@5cWsj79gG=a(mU+s+-9~xa@+!yC z3}tZue@C50reYD4@;(>oKL7StOcj3+(=eBj)!68P`C=31H>;|K^$UzL;F-ee?C_#Z>e0MN92Zr#B4)4ii-Sy>?A7rUcu8_6;oj zIZsVOe$^E2)jm6)Q;5t8UfiUf)jvv?s@tm`Sgv^~3soQLKl^?8-#;@W{_gvgb%ktp z;~lr_|8soqcsI_rIY$m&*|%%|=K^R5h9UQ)@{Y`fBvV~Dyn8)WZ)7p0tk6;>k4*X* zrzEjx-_G-J;w3Li`eoQMAq8%fS3y^cqFbtI&|VaPuaTnRD~+IG*EXTfRJ@IoZI+TM zTr~cm4W}_k8~+~lthaZNGc{H^5%yvEe!arMcbPTXK4<#e_T0ec=z3)&hYEDpWagHw z5InsEP8kF#6iyLm&Qe*oUP*lMuOnkUcF9)9?Te+`Vb-516eA7XuEW$c{##N{=6}q3 z%PESoiffb+u~ss!s$P^m;dprC2Q?Kxwl-aVGjmZ;Jqhzr(6Q_ai|5W$>hsPT*HKOg z6Zcno#ms$y@s-=kjyxCsaR89g)Y8Y>60aR3xM~3D185WA5%*AP>Lxg5@bJettG-?|87D(Qv{tCsG%(Q>f!%^~B7iy%lmSzeQM9W) zaW2usF)IJ&Ae{X;wH1^|6)CPol(Hkd*{|c&R&t_*XbTN$-Sc45P>4IXnqqQa(&c?I z0twQycn1XRYNbfcT>Fl5$mp+y_N zPL)lY8RN732BmU0Z2?4uHQKV`7d5*GdeVZ`Knv2_ba-;CFb57bkd20R9K-v+=D^T| z-PXQ%fi^6T_+A87d~yH3>kw~KWug8P;%3RzWW1RK10s@Sb7I#@HB`2y2yhWtbrnv6 z3p({7f0M+=9!i2GL1Mq#fE;_^8RX7xoxOP`Qz*b*QW=b~1ejwZ(n8fX+gHxRO?X?S zWCmmTnmL~cY;|LM_>!2A$3mK|=!CgPIt4+G>jpE;*+H=^uMPj%WNJRN(f9OdwIy<| zdamINL_ANcrOrp2Y(Bgp9jIF2ixxWzIz+AKl!h5qys~IK^N*-!UYdpOoIqadERYwW z6roNo7bz*O2EJGlC{5F0uGJb0eR~ah3SyLZ=YU$7?tqog=pjA@q{1Gie=c zwfM}n2%L4uwej7y=hUaC2(QEwj8P8J1sb%$DHX)3`7;xVbsA{1b)`y{5hsP4MJ%PW za0_^89zOPo;=!?g3YcI87A7%ZBhZ?3AF{(3X!uZY01w36^OOfM z9Cr1L*l_D=r1272D4xN4m^z8rx1jLUR$)bg3nmCf&3%CiWN4mDBNn+TSI)t``FHJF zJY-MqBGhIumTdlroube;w2&Jdp0uF|+o*DEF4K9^fJJz=+cYlC|L__m$ zCdI=dOR%+!=C!hC#FAqX)>$uN9a0x-?1TxE=~+o!}jSJ(Qv39tPYJ$#pWREf7AV%S@}PUaayGC!4k)6(`#IvgEM1BI1^BKfHWzE1=Ue2w*bKj z;6Hp%$b(->+F6gzyQBgwSxXR>zbdp%3K>sIeu!^1Pk{Fb`6rW4G93m(WJlet%}Hk( zGsxLI=5B853LBD`^iK)-{Ey-GAEtfM zE=SxZYv95W+;^&&UgONSMkct9;fEuMfKFn;f#SY47f=k1%JF~{NKXce52aJk-==W#Mp5uDRG#JXV-52-j??end`-`>r`lJAovF+db)kj zN}SOjZWUxEVDLj>+3`kz_=?tO#g2nBxzgU-QNtY(u;Jw=G%NH)nkt63)gk$?0T5O_ zZ_8}QGS51+Q04-kNM_A*I-{E{?Xu1C?VHeGd+sWhz2-|Oai9VeL~358Bxl}@f7ID= zRiW0LzRUyO!e@FMo5unz9=h>)G&D-GZIG6+hKqTe7vCidP7F zQAmSE9|Z`PI724v3r=3#YBzOg@%z4Bc1q`zjTmk0vO-VDBiJLu>!( z-SFWx^m24=`771+LwI!{h;pC~tYS*moEmmfv%ujG3-(Snim>^x=a3Cy%1=3&Gq?8G zO$EV}{qvQ%D^4?*SXD))E&P3%;c0`>Fxih`jJ#)yP|{1Gra$?yF_1>PDuBZ-;1Nh8 zE6m`K=DHW#Lgrd@w1C<$PXI#5FOkh zv&+Eqw-^jzwUXzl9N+V8AmD{?UO>qTeTpQlsNxFU> z8u61Z%jKpyI2u=3!fggNOEc4!nGz1lXp|`5lHxsu$nmHm)@(#*Wx4Rkm7n)+s+t7q z_MBrT$v^Opbig2T1SR`pZ}U&l%A;A$s|TamQPW_Y>=nKiWC@2CHgHw?^q2MJt@*oR zOttjh2s&}dlbBuhqUXf^DpQ#*E>oFmhAmcZZ&?KAM$e^UU$|(-;%oAE`d4Pp_;@bH z*XrUj)2OKF%KCp~$mGuq*tI=25!#F$5GI3P%y1ce)z(l8g%5sJqYI0wC0H|9ER|1& z2>Aq9ifrewz2PVH_CO6gbI=H-Eaa6`XzMyM%-;_xd@?KhJs)b*=;w6bH9F(sHKZ-4 zdKdwnA4k6*Zx*ox>$ULRRy1la!8q+-vIa2xNXT<_Ppp+BbSB5)tihNKg104Fo-qWf zk00}yP&xgM!wk&9s-#rr9F`xYS+Qi8a49_C1c_6|V|}CYwg8bk{M_6o#J(gA4^c?U znO!RBtAr6Xf>i0&18g58OGi=(>pw{1W*|Rg4-f4MWgbc|eeTf?f}%JfbgIAU&nCt6 z4|L9630aZW`aMBzdJ?m9rLbIQ6IYQ=1ZWIDRiyEWiuUNGYxZuS83=Tc!HqSdqHJ7X z24yX&Dy%NW*MN5OVWE(A(?nG0wQnr-r3e^3=KcF`V87g|uV5q-)=@UF-4JFk^Rku( zn;%0!o2eJgZVY9=f5x4s`b1+$r(^L!dK(Baw{ZCD3%-pS^(y~>LzEy)#r-$(ms}QL zFRuQxG>%NDvrzums|+n*!|ZE+_V?=}O?+-LZ73z?krQBa`HfRLRLY@$Z!s=f|D+XN zubs4w&22Ghu8n`@YL5|`x~-aNRhgergkbaH7etg!69&U});Bl4O~_C&v-V^(wmpI> z2>_QuYdfx5pWS%4o|N12iedQ^MNT*gaZw~wwf463@7Yz2!G(i}DkbgwbzuUpf-l6x zTm_7;pou}M-qIJ!)ii{9#UaT_eE^~8%BQx=rMB1J5v(1ffm~akqJ<$?e#Jbp z7^CZ0Rs5&jn@-+A_2K8Lnyb!fS1g~iGU*mhuR^O$M!TRv(a&b>pZP?YmMpCMztlYB zY_t7r6Ks$UVCUrifB1Q=hP4`|CU#w0Nwpe9j4B?x)~W()<~}NFt~9ek(SW!FDMsiB zYDmilJ^bzhLN@aPy}aiy2am*B&jJE9C1X8s$dyEK%IvX3(-k{5Q@IH_bl1CR%w~R( zq@0}Npjn3&;HNdisHabzEq?*nX{Lh-OZ)0;owIv*`6bDe-7k#>&xUy z@%obyB4dIy9L92Kk{hZ)-_P&;61kqmA)bZ)o}?GKkhJKN>ZJ$KHGbwkHd#7%Mos+` z-0bWW!c9J_RZ5d%a+$0G+H~pSKUS@qgaMNKD1HMxSoqRS$#^Mq&$2$VV+RSc~se)d+!Js+? z7kslzqbXRZ;lhB!f5wn*5H$jy(PMsYqaSglI$s5rHM2ITAZYtqA4p;w+cCc>>`+|4HYt6oUz) znuf%N$vt^bPLrM#XM`Ttz{%D4e8|$PcRmH3pJsAcYD4^^zPu=7>R zoT^bK68=`mN2HSame8kj22f&u$o3qs5Imm1urC&4;4I?Y1WQoZG~A$a_C!jUtwim@)ej+Nd9{lUC5{BnWz zy&ps~CvxW~_wQ%s-=yskoaz?Z9XIOKqV@3RPM5~%53ZbI96Q$}T|YgLnBU+!Wu5tZSLFZf z=vc|C*szoTIv^8O}l?U z$yKbKZnz4r+q8!38xM2Hf<^?(FRj8d$8_A=KPPka5+BCQDQ6S4HZfvDsP2VP7QIgi zD{g6O0`#jnVI>98e~{X}$CW)${31s8037mooNY${=wJ$Qjb0vVQ(4^*hj_||FVfmj zfwkl6D>l*b3uxE)1Qysm$i`s3_wPKKb|A53C(xl$@sVTJCmVHq#zTZYYLC$c8_I zdIo38ZnTtw1`5kdH95kVFzAaWsG%l!pY>2b=7f5`2@<$z7nM)~A-K_gj5QTf0vL-h z3N1nuiWa)7e=DlX&AO%LN|O#MXhkaQn~8cGM|Hw5^v*`bLqW~nL3);=AR9{W$3}oG z6W9r+A{l}EuEnQ6bx_sy?&0b>XUZ|zd7u&X#Ux_Tl7~JQ$haL!X*y?)XTSpHrZ!jK zanBIY-Cq|AahtVLRVtU8k(}`?#6wzdLDygGPz~*Pf6Hbe$!3}0OFSa$J0$OqK4Q9F zf5Bvh8KjC}FQ~DjZks**VMBu&J7}pvSE>E2P~$6*sH36NuTjsVnI2Yf3i9`ic0Zh5oh>u7M zLEm|xiOPZ<0W)qrIx{{-=iE8~Ug+ZZVlg7igvEk*&Hciig723r5s4O^K!m+kS6xwN(CG|ZcF1P66lm!KZnTlJQYI@?kAgq8 zf4Jcrs1;V2(&@PgcS~14H{>3Mju zX58?jezAWaKZN_V8~2?12lZpi%CSu?f3t}OCQ$tZ$;Y9$8Mvmqk~*>LTDT(to`MK@ z$|H&>Z8Qi3tgH*`i#YoDo2n7`2!=;IV}0+KKZh||ND?kFdm+-2heRc>zy!R74dGJT z%k8hf*mDeyz?j6Xr$br2*!p2;|(SKE$(u@dcMRIa>r%sLofBnHoM9|uVjj>xfAyLTjQt_f&R?ok6U3ByMK?boqEhDMD3$h{4wu8?W-(Plyd7=8?1$S(^h_ml zHs|PWd7dUa%c|*aZK|uOWvi-8U;^3`{3KC}3AI`_`=Jm<`BTt^?*I8G419akI1 zYW;^*C%j+q*8WAM=zZR%$XZK^a^nRrKd;?3Kh|3ElfO}Z%6m3H)>`s2+N1oG#%zA9 zf2U=tc<8NTFG@?YB-2L59%W=k#%E=7k)4XCNgQc^ojsO`ZD|H%`~ucy z&Y<+&%fPa1{a_!X4>Fvm%dKqDBqk8lx_(GDRb|~S_1Bo~N~Zeho|SD(9$uxw9o71< zi>FL(_;;8Qb|!c;!c^C?VUC(h>Yln}y+yAw>lJ#P%dU*