From bd107bcb95bf5284d6e46200db3161a33e64b985 Mon Sep 17 00:00:00 2001 From: avallecam Date: Mon, 28 Sep 2020 19:06:18 -0500 Subject: [PATCH] add hallal + lewis refs --- README.Rmd | 6 +- README.md | 37 +++++++++++- docs/articles/howto-reprex.html | 54 +++++++++--------- docs/index.html | 25 +++++++- .../figures/README-unnamed-chunk-10-1.png | Bin 7776 -> 7503 bytes .../figures/README-unnamed-chunk-14-1.png | Bin 8114 -> 8250 bytes docs/reference/unite_dotwhiskers.html | 7 ++- man/figures/README-unnamed-chunk-10-1.png | Bin 7776 -> 7503 bytes man/figures/README-unnamed-chunk-14-1.png | Bin 8114 -> 8250 bytes references.bib | 24 ++++++++ 10 files changed, 120 insertions(+), 33 deletions(-) diff --git a/README.Rmd b/README.Rmd index 2995362..9fc82b0 100644 --- a/README.Rmd +++ b/README.Rmd @@ -312,7 +312,7 @@ Feel free to fill an issue or contribute with your functions or workflows in a p Here are a list of publications with interesting approaches using R: -- @Silveira2020 in the strategy 3 analysed the survey accounting for sampling design and test validity using parametric bootstraping. +- @Silveira2020 and @Hallal2020 analysed a serological survey accounting for sampling design and test validity using parametric bootstraping, following @Lewis2012. - @Flor2020 implemented a lot of frequentist and bayesian methods for test with known sensitivity and specificity. Code is available [here](https://github.com/BfRstats/bayespem-validation-code). @@ -351,12 +351,16 @@ Flor, Matthias, Michael Weiß, Thomas Selhorst, Christine Müller-Graf, and Matt Gelman, Andrew, and Bob Carpenter. 2020. “Bayesian Analysis of Tests with Unknown Specificity and Sensitivity.” Journal of the Royal Statistical Society: Series C (Applied Statistics), August. https://doi.org/10.1111/rssc.12435. +Hallal, Pedro C, Fernando P Hartwig, Bernardo L Horta, Mariângela F Silveira, Claudio J Struchiner, Luı́s P Vidaletti, Nelson A Neumann, et al. 2020. “SARS-CoV-2 Antibody Prevalence in Brazil: Results from Two Successive Nationwide Serological Household Surveys.” The Lancet Global Health, September. https://doi.org/10.1016/s2214-109x(20)30387-9. + Kritsotakis, Evangelos I. 2020. “On the Importance of Population-Based Serological Surveys of SARS-CoV-2 Without Overlooking Their Inherent Uncertainties.” Public Health in Practice 1 (November): 100013. https://doi.org/10.1016/j.puhip.2020.100013. Larremore, Daniel B., Bailey K Fosdick, Kate M Bubar, Sam Zhang, Stephen M Kissler, C. Jessica E. Metcalf, Caroline Buckee, and Yonatan Grad.2020.“Estimating SARS-CoV-2 Seroprevalence and Epidemiological Parameters with Uncertainty from Serological Surveys.” medRxiv, April. https://doi.org/10.1101/2020.04.15.20067066. Larremore, Daniel B., Bailey K. Fosdick, Sam Zhang, and Yonatan Grad.2020.“Jointly Modeling Prevalence, Sensitivity and Specificity for Optimal Sample Allocation.” bioRxiv, May. https://doi.org/10.1101/2020.05.23.112649. +Lewis, Fraser I, and Paul R Torgerson. 2012. “A Tutorial in Estimating the Prevalence of Disease in Humans and Animals in the Absence of a Gold Standard Diagnostic.” Emerging Themes in Epidemiology 9 (1). https://doi.org/10.1186/1742-7622-9-9. + Rogan, Walter J., and Beth Gladen. 1978. “Estimating Prevalence from the Results of A Screening Test.” American Journal of Epidemiology 107 (1): 71–76. https://doi.org/10.1093/oxfordjournals.aje.a112510. Silveira, Mariângela F., Aluı́sio J. D. Barros, Bernardo L. Horta, Lúcia C. Pellanda, Gabriel D. Victora, Odir A. Dellagostin, Claudio J. Struchiner, et al. 2020. “Population-Based Surveys of Antibodies Against SARS-CoV-2 in Southern Brazil.” Nature Medicine 26 (8): 1196–9. https://doi.org/10.1038/s41591-020-0992-3. diff --git a/README.md b/README.md index 1cc44e0..dbaa6ac 100644 --- a/README.md +++ b/README.md @@ -231,9 +231,10 @@ workflows in a pull request. Here are a list of publications with interesting approaches using R: - - Silveira et al. ([2020](#ref-Silveira2020)) in the strategy 3 - analysed the survey accounting for sampling design and test validity - using parametric bootstraping. + - Silveira et al. ([2020](#ref-Silveira2020)) and Hallal et al. + ([2020](#ref-Hallal2020)) analysed a serological survey accounting + for sampling design and test validity using parametric bootstraping, + following Lewis and Torgerson ([2012](#ref-Lewis2012)). - Flor et al. ([2020](#ref-Flor2020)) implemented a lot of frequentist and bayesian methods for test with known sensitivity and @@ -305,6 +306,12 @@ with Unknown Specificity and Sensitivity.” Journal of the Royal Statistical Society: Series C (Applied Statistics), August. . +Hallal, Pedro C, Fernando P Hartwig, Bernardo L Horta, Mariângela F +Silveira, Claudio J Struchiner, Luı́s P Vidaletti, Nelson A Neumann, et +al. 2020. “SARS-CoV-2 Antibody Prevalence in Brazil: Results from Two +Successive Nationwide Serological Household Surveys.” The Lancet Global +Health, September. . + Kritsotakis, Evangelos I. 2020. “On the Importance of Population-Based Serological Surveys of SARS-CoV-2 Without Overlooking Their Inherent Uncertainties.” Public Health in Practice 1 (November): 100013. @@ -321,6 +328,11 @@ Grad.2020.“Jointly Modeling Prevalence, Sensitivity and Specificity for Optimal Sample Allocation.” bioRxiv, May. . +Lewis, Fraser I, and Paul R Torgerson. 2012. “A Tutorial in Estimating +the Prevalence of Disease in Humans and Animals in the Absence of a Gold +Standard Diagnostic.” Emerging Themes in Epidemiology 9 (1). +. + Rogan, Walter J., and Beth Gladen. 1978. “Estimating Prevalence from the Results of A Screening Test.” American Journal of Epidemiology 107 (1): 71–76. . @@ -373,6 +385,16 @@ Statistical Society: Series C (Applied Statistics)*, August. +
+ +Hallal, Pedro C, Fernando P Hartwig, Bernardo L Horta, Mariângela F +Silveira, Claudio J Struchiner, Luı́s P Vidaletti, Nelson A Neumann, et +al. 2020. “SARS-CoV-2 Antibody Prevalence in Brazil: Results from Two +Successive Nationwide Serological Household Surveys.” *The Lancet Global +Health*, September. . + +
+
Kritsotakis, Evangelos I. 2020. “On the Importance of Population-Based @@ -401,6 +423,15 @@ Optimal Sample Allocation.” *bioRxiv*, May.
+
+ +Lewis, Fraser I, and Paul R Torgerson. 2012. “A Tutorial in Estimating +the Prevalence of Disease in Humans and Animals in the Absence of a Gold +Standard Diagnostic.” *Emerging Themes in Epidemiology* 9 (1). +. + +
+
Rogan, Walter J., and Beth Gladen. 1978. “Estimating Prevalence from the diff --git a/docs/articles/howto-reprex.html b/docs/articles/howto-reprex.html index 6074f99..1a45b0c 100644 --- a/docs/articles/howto-reprex.html +++ b/docs/articles/howto-reprex.html @@ -455,7 +455,7 @@

} toc() -#> 367.84 sec elapsed +#> 241 sec elapsed outcome_01_adj <- out %>% mutate(rowname=as.numeric(rowname)) %>% @@ -530,35 +530,35 @@

E outcome_one positive -75.9% (65.3 - 84.6) -78.9% ( 67.5 - 87) -83.157% (65.90 - 95.5) +75.9 (65.3 - 84.6) +78.9 ( 67.5 - 87) +83.140 (65.87 - 95.5) covariate_01 H outcome_one positive -25.0% ( 8.7 - 49.1) -14.8% ( 3.8 - 44) -0.421% ( 0.28 - 29.9) +25.0 ( 8.7 - 49.1) +14.8 ( 3.8 - 44) +0.423 ( 0.27 - 29.8) covariate_01 M outcome_one positive -65.2% (42.7 - 83.6) -44.8% ( 16.0 - 78) -31.544% ( 0.40 - 81.1) +65.2 (42.7 - 83.6) +44.8 ( 16.0 - 78) +31.450 ( 0.41 - 81.1) covariate_01 E outcome_two positive -6.0% ( 2.0 - 13.5) -18.2% ( 5.0 - 49) +6.0 ( 2.0 - 13.5) +18.2 ( 5.0 - 49) NA @@ -566,8 +566,8 @@

H outcome_two positive -5.0% ( 0.1 - 24.9) -7.7% ( 0.9 - 44) +5.0 ( 0.1 - 24.9) +7.7 ( 0.9 - 44) NA @@ -575,17 +575,17 @@

Yes outcome_one positive -93.3% (85.9 - 97.5) -96.7% ( 88.2 - 99) -99.900% (97.02 - 100.0) +93.3 (85.9 - 97.5) +96.7 ( 88.2 - 99) +99.903 (97.00 - 100.0) covariate_02 No outcome_two positive -8.1% ( 1.7 - 21.9) -21.8% ( 6.7 - 52) +8.1 ( 1.7 - 21.9) +21.8 ( 6.7 - 52) NA @@ -593,8 +593,8 @@

Yes outcome_two positive -3.4% ( 0.7 - 9.5) -9.1% ( 2.1 - 32) +3.4 ( 0.7 - 9.5) +9.1 ( 2.1 - 32) NA @@ -602,8 +602,8 @@

positive covariate_02 Yes -100.0% (95.7 - 100.0) -100.0% (100.0 - 100) +100.0 (95.7 - 100.0) +100.0 (100.0 - 100) NA @@ -611,8 +611,8 @@

positive covariate_02 No -50.0% (11.8 - 88.2) -54.9% ( 29.4 - 78) +50.0 (11.8 - 88.2) +54.9 ( 29.4 - 78) NA @@ -620,8 +620,8 @@

positive covariate_02 Yes -50.0% (11.8 - 88.2) -45.1% ( 21.9 - 71) +50.0 (11.8 - 88.2) +45.1 ( 21.9 - 71) NA diff --git a/docs/index.html b/docs/index.html index 41d2ea8..8c14605 100644 --- a/docs/index.html +++ b/docs/index.html @@ -259,7 +259,7 @@

Feel free to fill an issue or contribute with your functions or workflows in a pull request.

Here are a list of publications with interesting approaches using R:

    -

  • Silveira et al. (2020) in the strategy 3 analysed the survey accounting for sampling design and test validity using parametric bootstraping.

    • +

  • Silveira et al. (2020) and Hallal et al. (2020) analysed a serological survey accounting for sampling design and test validity using parametric bootstraping, following Lewis and Torgerson (2012).

  • Flor et al. (2020) implemented a lot of frequentist and bayesian methods for test with known sensitivity and specificity. Code is available here.

  • Gelman and Carpenter (2020) also applied Bayesian inference with hierarchical regression and post-stratification to account for test uncertainty with unknown specificity and sensitivity. Here a case-study.

@@ -304,9 +304,11 @@

Diggle, Peter J. 2011. “Estimating Prevalence Using an Imperfect Test.” Epidemiology Research International 2011: 1–5. https://doi.org/10.1155/2011/608719.

Flor, Matthias, Michael Weiß, Thomas Selhorst, Christine Müller-Graf, and Matthias Greiner. 2020. “Comparison of Bayesian and Frequentist Methods for Prevalence Estimation Under Misclassification.” BMC Public Health 20 (1). https://doi.org/10.1186/s12889-020-09177-4.

Gelman, Andrew, and Bob Carpenter. 2020. “Bayesian Analysis of Tests with Unknown Specificity and Sensitivity.” Journal of the Royal Statistical Society: Series C (Applied Statistics), August. https://doi.org/10.1111/rssc.12435.

+

Hallal, Pedro C, Fernando P Hartwig, Bernardo L Horta, Mariângela F Silveira, Claudio J Struchiner, Luı́s P Vidaletti, Nelson A Neumann, et al. 2020. “SARS-CoV-2 Antibody Prevalence in Brazil: Results from Two Successive Nationwide Serological Household Surveys.” The Lancet Global Health, September. https://doi.org/10.1016/s2214-109x(20)30387-9.

Kritsotakis, Evangelos I. 2020. “On the Importance of Population-Based Serological Surveys of SARS-CoV-2 Without Overlooking Their Inherent Uncertainties.” Public Health in Practice 1 (November): 100013. https://doi.org/10.1016/j.puhip.2020.100013.

Larremore, Daniel B., Bailey K Fosdick, Kate M Bubar, Sam Zhang, Stephen M Kissler, C. Jessica E. Metcalf, Caroline Buckee, and Yonatan Grad.2020.“Estimating SARS-CoV-2 Seroprevalence and Epidemiological Parameters with Uncertainty from Serological Surveys.” medRxiv, April. https://doi.org/10.1101/2020.04.15.20067066.

Larremore, Daniel B., Bailey K. Fosdick, Sam Zhang, and Yonatan Grad.2020.“Jointly Modeling Prevalence, Sensitivity and Specificity for Optimal Sample Allocation.” bioRxiv, May. https://doi.org/10.1101/2020.05.23.112649.

+

Lewis, Fraser I, and Paul R Torgerson. 2012. “A Tutorial in Estimating the Prevalence of Disease in Humans and Animals in the Absence of a Gold Standard Diagnostic.” Emerging Themes in Epidemiology 9 (1). https://doi.org/10.1186/1742-7622-9-9.

Rogan, Walter J., and Beth Gladen. 1978. “Estimating Prevalence from the Results of A Screening Test.” American Journal of Epidemiology 107 (1): 71–76. https://doi.org/10.1093/oxfordjournals.aje.a112510.

Silveira, Mariângela F., Aluı́sio J. D. Barros, Bernardo L. Horta, Lúcia C. Pellanda, Gabriel D. Victora, Odir A. Dellagostin, Claudio J. Struchiner, et al. 2020. “Population-Based Surveys of Antibodies Against SARS-CoV-2 in Southern Brazil.” Nature Medicine 26 (8): 1196–9. https://doi.org/10.1038/s41591-020-0992-3.

Takahashi, Saki, Bryan Greenhouse, and Isabel Rodríguez-Barraquer. 2020. “Are SARS-CoV-2 seroprevalence estimates biased?” The Journal of Infectious Diseases, August. https://doi.org/10.1093/infdis/jiaa523.

@@ -352,6 +354,17 @@

+
+ +Hallal, Pedro C, Fernando P Hartwig, Bernardo L Horta, Mariângela F +Silveira, Claudio J Struchiner, Luı́s P Vidaletti, Nelson A Neumann, et +al. 2020. “SARS-CoV-2 Antibody Prevalence in Brazil: Results from Two +Successive Nationwide Serological Household Surveys.” *The Lancet Global +Health*, September. . + + +
+
Kritsotakis, Evangelos I. 2020. “On the Importance of Population-Based @@ -383,6 +396,16 @@

+
+ +Lewis, Fraser I, and Paul R Torgerson. 2012. “A Tutorial in Estimating +the Prevalence of Disease in Humans and Animals in the Absence of a Gold +Standard Diagnostic.” *Emerging Themes in Epidemiology* 9 (1). +. + + +
+
Rogan, Walter J., and Beth Gladen. 1978. “Estimating Prevalence from the diff --git a/docs/reference/figures/README-unnamed-chunk-10-1.png b/docs/reference/figures/README-unnamed-chunk-10-1.png index 39dfc0d45fbf8f232deb02f4a340cc718db1b5c8..74beebb591e3d1ea83ed4af07d0143d61beeda75 100644 GIT binary patch literal 7503 zcmb_Bc{r5o+vA-XGiEGfFMCudQciMw}?ysF<xVAS9d| zceo)CNC<%tNx&lE6_56Xaro=8>mK)=@HYZshCl=&%rX#xRfdL!@Piq`%nT7|meEV% z|4~S!G!lhQN^9}<_6`()23GOlRSJnhA*E3$@IReGA5Ig%rIFIo_^b4^7G8!KA|nGP z!=qD3bQma&L~o(fhllCIysAJ%6_3}FMrz^ldNa&=c~!kU-Y}gq%;SxXjqwl&7y{mg z_uy3?j}Hq2!XP{zZ@iC+hq)G`I(kPS5Ry&&k4SUay~7BEijC6_8+TgnL|*8j-tE`y zKRQ>{MCLr`buM-?AEOwWXKjB|ddN7^?+1@pwykV!y?>_02@#;fRQ)$6g&$arG@DG6 zYbPQ@zq3mnfmpO5CtUjX$6J@(AYzT9;c6!*4o$x4C7wx^72!0>zD7<;SXJf|k{3Ee zVyccgmTHuTB7luzHr8YjFceAfa>4=t9vh;Z0PV6R5WK9(y$F=2Okx5A0KK9b<(CO4 zlmek=I|96#tf5h!09A=#OV*-9Ks2gK!kWBp4Hm0sBDVPR%d=ATfPlKi&Hb;l!pBsF z(5$uWiOwgEu?w9vLzllEo)QvHya-9fFDyhZafe7t?7p_6|J+N0TzK5>h^YnkQp)^v zO7+Y^jV9H5%&;Yqu84xcBlWG1?%J_5*c0?G?75-VLg@6Ag)T=*+k($%H2I)6YowGD z)Pdv}|0rTcw|lVs?C$aTF4tpEcrCLkhr130z|`32HU z9H00E^p2*%LTA=nbLY?-qW8N`tR00J7zGzPkpDVC!w{q&nVi~)vJ?Y?ulP)T#QKW0 z6^^KRscDsq3RZe?wR#_yRa2j3*^MT(R{~{^SHUO2%D4UwRI8TFsxofQSrJVBj(78Qh15rAfGgQrl2H6{Xm2%aPl^9ENh%7idG8X8Z@`6#F(5o{WX#QgTYfjRojq%o0WGd{UYkZ8gimnfDypdgP<5|)>^)1G)zJI*WJ#iej~}* zudHkVImJb>sC*+al=D{p5$|MZ(dolD&>@daalEH#k|c<>qP0{HQk$<wtM}q$_a6B zY)Jpw%HrBJ?4E%wOPtA*6rkJu&OYOfi5RRNUt>7HU~Wv8&1NvGuxi3M@5K`>yPn9> z*zvszjYvTj(0Byx)6Jf-uw(Vdhk2p`1oFu{&=m)O%#Y6qRSo*yl4%ii#p|c|Dq>>F8vmP?&`n0BU04o zKj`KleKSym0yL1g31!{djCa7fZYne z7|{YZhHm{>4oq!8IcxCEv0a+Rc*HuOg6g`j(0+Ig?EVSnx!~sf1y+pdcHQe@89aX>G@m-5>K&U zS3dQVVh-?AutgPA{DqRfzWmwc!Eos}7!cju>Nr(8Cpj)hod!<%yipj8C3tx$aH->P z_Le_#D;ueaC*CYPh6f!)wWv%)iW+zJ+qXIM&D~%ga2rRgjNY{aObNw+L!wmeYoXSA zK+F)PCg8f*cuV=T28;kl**P7mN3kZgVP@r-u?LK4<74n$j0ietamNPVp98AZ?=E0l zcbcdRo$4kuBmxx{XQC&?CoNxx2$+S5l`GQ7I~Lv@MzBY};`X4_)H)2@FN+D>6?%yy znrPR%I;|JFTP7+$jiqkkCjzvCkka2>F)H9y`Bl$2r7s;>%eCL}R}Tw|9APIp6n%BF zF1SL578@|~b=Ov5Av?FG{2AcJ>3?!v%-5xBk1dNO0+&RB4f(WK5g%HX_9#$Bm;w`8 z6==2pd0+ztUcmjBsRIO>B&tiNx<3g|tVKrumVvnjq$|L2A|!Nxl@)xwvpBBe+46*n z6xMc)Tqw)bq{NzKO7JRVMoX<{pOBQ?2w-Cu87tTXE_Bo@$wHVabij=y8Rh3#qEZ)G zYNdB@Aeijpi-eN#%G=p|&Ki(vUiVR$)sm%6iW(YGYU07=uoJ>|#cEPrH_q7 zU;zB9ZjQ%jwPg)CsJczOrWd}fzgwBC21~L-N~LQc!Gap*d)eZrAYq^gCglxl(ESfX zJcCzS7+9AnlKkQn1KKSpyR#BRwS;{NG1}sKqk!jxr%>aPgdVl`(X;B0!q&zMUs6{1 z+b%MgPNK);NC?61Lp7}I45eSjOhu=$%2J`?qquiDyIz!@B+#VzrH!3&R|BiEnV&gA z&}oBVG%#kYv8PiKn6lw_PimU$zW#M^zd5!R14HaVr+p3+2Yr$zj$uF!90v(^pyP`> zS>6I&SIaPhCeUIkbz)A;}qJn!*CJJ^MkjA+9`LXAL`EVT# z1RQ=C--KR8FxQk*cdC@uAA(HB{gDAza9j%>Z5& zET|dyqL3z9k1Uh2hIQ6smH_(&WX)<>wpZwRO%^B(0UpWxo)@qaHsU>CN~&;V@{2WD z3+xT-lj$9}2jI8E4)F(WvGwBoQgK85$kuoB2$a7lD1mhLwV(e{tjKH`l%-Zd(p_$= zrTxe`Nk>r%&RD+ws{yOnn8--Ke^!CGQt0Y4YbV#w_>kXV4sSy7%RlGB3Kp3vj>g1g zc2yjAthd;>duC$)P#aOVy1!kH@xZ{Aj8BsGa&18Jmc3_gtKt_cDrWm;hSH5^nob!P z3gDz^Hn%s_pbw$dMOW44bbeg)6$OX*nh(Nzyrm5-+F;Y0%Qua;-_ODfuU3{mpOs59Mez|qy)A+L#a+Ka}=mPQ5_dEjzB4#00;JW8c&UM zY40HWfb{9z(l#v9M{?ZkhQ+D49fO|*u8gfse5PvgwE`YvOiXjx%J6Y0ni$bH%CO2l z2j?4Ou$U?wW7xPNMm;_rb^|rrhA_!nIa`)43Y;-Uj2O~XvF(!=dEiWgv{#4AZ%bXi zqWRN`pgvknb?2+Dqrsp>w}d?zY;~*>Q3LkZb;Z(jH(-evK?gb_R(syQS1k5La3Ac7 zSMGRx>J&3Q|HKj^k%OvI9NrdZs#exPLk7~$8s&c!yg2qD&##aGK1B612jo!wQMngj z!sz-?WEZZs4vsv$fapfFGc>RkHOsh`3_WMAG3o}pF(tAl%Ui{9PT1OIcOw9A*bN@3 z2edW=K9#flXJ7kA5?<@rqGR`T)@rgT#M|MwB&GSYT*5XWtU^4=b^|bFAEBi8VU4gq zvSpD1kMIAL@EyqU_U{w*S2kIj*xOq->lA6RslUBBW41#uC+4C{)2iJl&bENmm#Tn! zcEV)k!8(=#9lAL`s+#?aG)spXGxWCVFRYUo)8~)lD{LT)!5W`<<86UbE5)rQ3#xk? z;e(6N&Ci%vG1_6FHPy+LoGoTHDsV?H0>dh?Vd?KzMLOk8C{A|8-B4K8nv3i&8vQ{mJo4=EzRe#VFD&Lgf@uV~`B~kfEaNA!G zIjg~DEv|!$vDRKNTijltS!XD;6{{xt8Jms z6wcV(57=etuQ+1w1_3z2kM-D|o$lzxiVp0wUCD@T8K2Bf?}J?f?qmMIv# z_1<`=ur}U=f)*}9m4cbD6r0OvIFFD0eiu%BaQ{HwK>fCScN3K-S>;5^?#=Lpt6IW& z$f;z|V)r+IGr6YV!-=O7f6Bw9+?s57^HUZEOGLtT@e;b!t(I!EHVT}1+AkC&6Kax% zf)ZDTPfIzItJ`XiI0cX+D@1?sIoh4KwcA|4=Vk(pm~2G~CwSd}$54K&Fe>Idlbz#o zT@TeyI=fBS>yauy88=kf8wV|I{N=!5Xy+&tJHs8VwnZ&FSs3!v%^Ogj{0Jz7oYE0$ zmoGwv2ZdE0vSBU%o=#G(L;ecm7%*KcWB3k#=_8hyt1C=O)~t-5E-U)dq?)6={6Ce} z)nsNuefLL4())-Ts@FVCgzrN*Z4xS>FC8ia#9u^GamNGD!r| zdL+Z}xg`ps72){;9xtrwo|c9+C<*g(K_h}CjUBHj$P%LxxP6Me*hx@5zV;*MMUt>W zm>4xnYr=7%Zfh0GyA`)M7+UGnqM!JuJU<@_iqEcEy;ZbJ5G(?>f|h7?%r6 zP2Rt+x3i-F>gz}Uz$g%-S<7(dJTN71L8TEG>6K6{#$WOt!30QIeUjisRp5SY+xW(| z9b+`H!Em+GcUiAZoomXWNOC2@^ik`cvuDN8lxsECWYD~#1elfW;nPk7rrNh(`IylK zLiQpltywqz?8vwS)!_IB6f@~ZzHUxyv&B{G~D|@(hHH z^@k{<20+S)>l7y?LMo*$0rH3Q`eg_iFKt3hi5@68GpTPSf5IvS306BX##~tanJ_!v zc%~0xYfqdc9tg{#5*W6u9>*Ciwg_tWbI-B|CkRjTuT6Xh8vXv zN%u@8R6sRX^?T&uMSqs;ye1WO7_$c*x16K&n~PSt{UlKEQxR;^im{faMQWV{z-?O< zX@MfoUmrf&LSi8=Kx^PslUj?;_i50xTaA-3Nod5DWIloV7EKM?HfZgVVQlBY5Q#Ivp1N*#A@f3;m@7MAaZ7;#mdMp51e3PNJxuazV z&%(Ai10{J(?a6@`&q0yA6Nf*Xl`}$sE?bx%@3w8w-{e3p5P;P~#h(6~gW#hcL8=pR z54-e}&q7M~@xE=Z7ehXkxUp`dEH7=m*qQMhz647=GBx%VUJFNAT9$fV3WrTugr2pC zo(=NUrA7f$FAi3(PluGg;`A1_`B&}&FUm0*@3+BLI|fwOoaj59HDXMcV|anV+5W4* zZI|V|-ng7}DKvuD#6kI+eswgRR7gn-H*WO1-=#qLyo(BqgIxzvU5<7vih<{mJjS=I zTjGiv%U#u#t-DpC%e2a!-GJVYHFhkTSKPO|+uXy`9m@Pvl;5%LTVTiYI>xUb>!dFI zzV(U-c&h13N1Z#2pb4n+hBc{YQPQ+i^K9Jlkew_JhVTog9JDx zkIz0(>LYXno4q~oKvogHse3$hwEcRdmpKT`}M8?YbUx{vP#f=6Ke-WVUc4(ZR`S|ElA}yzj8tO|` z>LzEJdI>$NQL}|=+n-H`(T^Seaelk{Zg`Fmjw%W_ zV*MNaShK2~=YoDe+v7z4SnK@|1KPPt-5h4j%wDPA=Tk}#m;|hrY5S~VPZ@p+$^2S; zk_yWqDNXv`h}}M-UpX8% zCRLulAa3=;@7gw1;0S0{PsA}f{aHq@kSSP(R6bda8pgdbldQo^k?h&``ij}F>#UVG z3WRy{a%cWUn&4%H*KvjUtL@%L^ZxV3q5NwkD?y_?0s*w|O@qD(r-XmMuYiZ}|DRVN z0icKmdJ$!rWoXWp`va?MEuE&p-lb!Ub4}L5*2-MS0OST8eL)dO4kg;M>K~u$Sb#G! zocVP+zu))t^vRQ1AFj&{4))Rw!UBz*QcJPz1xw$ejXc~G_FSNHGo`aM$FnU>lh^>K zJ8({?ZMpBs>+N|`S(-*a#`^;gB!8*2TAZ1n{+X)+Q(+$(Yk!jHB%QC`m(MZQsfFsI zf(A}cw~t`j@KIZI6Zc%p^VsrCGd4u2qh@UIM5F0u-RFbNuhfGjvNVaD7x4=*p8=`B zv+)^SJgXb3wk{Vk0;Gt!1-!Yp>@zU^TO2LfyE|#~%?)l=$!aBSnm~Mv{<=R?oTR@m z0H!*&hTR+?BexbZ@$WU6^}CkM%ikeKx=z4Y%d_9N&?~)`G>P?b+a|@s$!^KB4%?5$ ztce0K0ae4?9}VHbW?7oU-KS?Rbf}1=X21RubdpZ9`hL(YxtD%6Vj1_+1mEYAv(s{^ z2EJR>-q})ynY=(fI@7c+f=_IGF;C&v8`jo$gA1RRrP` z+wf*8c(DlEeEIN;!v@0|n-R$XD6&=be|-!xtd@U_JhTEy182( zP(%U&ga%%7{QiRF`QD2IJXLVF2xxF_z$&JKcOH6~04cthDmY#F>-Sz?)Fuq7UW33O z61a)htnb$FB{kTJ1P!pLO;Z2xiGTXXzkgrnVODu!eyWqk8u>3qopw6!D77Ug{tt?+ BDP;fv literal 7776 zcmcgxc|6qJ_t&SH8YBi|Dmx(*CE1rKX`{$Kw&zJ?X^dUgKB+uCW9gBoY;AbTTDGx& zS}7qF27^zFEo)(H!~8z8)c5)RzJL6FfBarQujF&jJ@=mXJ@=e*&$;*NNlSA9zHNLw zJUjwN4j;7U;aLyy@T|LlS&vX`+gB!#KjkNn+n6DLd3f}AcszOZuJL$Q?%lf=Y3T9j z>G62#UF)-8{SYlIk}QY?7D>&Hj*gyeR?kW%LM2)di55vjBJ!_*SnxH8&6i}6l*FPI zBsDXy>G52@NtzZQIDR(kSOGMUXu7R^j%-!;8HW@R6f`L%%f zmB}0*A7}FLARI_H(uYu)OcpP~h_Eo3%!z&yKZ2T1W2=EuW?K| zK8$CiOZOJEwi_+uzEWW`dCc_SB-G$_raFZ40Snra!}YiqlF zX3LUe5?$e)k>Sl^@VAYinX@36!~9+{JG!#C8qXMuuUeXSg+(6EmWHg)v*xP{k!J`; zyB)#_g41re+IDF-JD_i%z1?VSw+5n&t)K$U8*D!?UA7nK`zIVr{;1efUvW8L69rS) zWaiG>y`uSkHp5=xveZ{Dw6^sHB&L&znPrp2gr#2pCuk~Kdm=%A)_675#_=ezfCFeFtiB+j&U@0c0{W5!%OUsEw|@)rKowW+lb=yzf#`RI1E;rw=9Cuqs$pbFVFa-3 z8#Oftst$+=-SQ7mMQE+Ws!9KOOs4+*`Q6Ud^)V*ySQa>Oem{*{P_+r1g;>kw+9F7C z7Qg88*aVJdx6_~`yUIpuFleFc!PPK6qW^v&hZl&lrqX~6&+aY%kgN|D=x^-zKf#}= zzaVB=cZw9b{vssDTCcVUm1_AXYD#KpTTDJo-+)3K1o3d9LTNx|DlbiUL6c$mbW#Ar z(W`|uOz)I_J}pWcxWLtDB1)_Hdh|dNZ)T!%%4Nu#jVuwOzqOyQ!iHb?saFvkr6KQ$ z$y~!fWZUB@62|9@$8+Vt0$&xZjx7t&bsRPO=@Vv_XN9&5X6I*?-Ymd^xwoqckUZ^7 zDItiZ2@o=kjYqW+yzWhe-5&gfWgeZjtiOXayCuF&zt#VzI*4z8?!tc(roI8>RY?(Q z!`;_0`AfeRXXakPT~h@*ZAD*&)K#y)1^mFa#t^Wqt^X+|a|3HzRQV;~=}XXj|5N%u z#7xxy@r&aaKQfJBHCp$9@k^Ie`{B321Ku&ebZ(E+X&&?jUE;7h0<0y6$&zM+-%Eg} zkig1X?Q_OFR>vIv$~K~_klj`p(l*wn6z?G2X^a81LoRUDEDnyF?uRugQb257_0kWu z_8iy2I>tNKb20Rel5t*%ZzUy=tpYL{iCNN(wJp<)Z;)sz#;mZ*$cdgCPeQddd6LsA z*@|ACG|f?nS^i!e+~dL?U%)FdQPK!4HG7o3=d8t;75Kg^NdL%}9uO-K)kT~_wLbX? zgwLCO^hdRt%UHV4nW6@oc@!#!K_^!&HKpHxDi@FdY@Fbm2z^?k^G&>p5I+RlP-5g~ zdT5@gN_!R#RiFTkSY7*v!_&m2O$M)Aj=)>jA~4Rn8S|DM*!Kh1kX_N=H4`$pwA8Aa zWGWcApg1#B{i&=^To`dq!$a^&M#ij{d)xZtfNm`H_N8cVu=Z3is$6>qVn2ZjX?$U^ zKv>===B)oou79WieS&fXSC}p;`GbCiB^=T2lojk(r^tQJV=`BMNWuaA5E!FQ|?ga!tZs2m<1FEFLpVUCBFF_8T0y&cTALRXc^(ZjzR}3-% z^#xa*-Qa<5`#J0Qm(>rUw6YY$F*JKAE}Td_CEtjv5WEO=ZoAT%b9SMyWjl-?yN7ErNSlrLg^e0TIB2}@YCK}Btn0&)Lx)gI)+Dh$D3+pNZVPu z@!|K2XjATVyQE+P+lg-i%wG!N77tDsymcO&9}3n`8aG?%`|XxIDJzI@eE#LbWjEcR zl^+hjAtT!)_c>c+jtGS|N{Zn?@+;?d%e@@+78gT%rB8QBqXraMTl)xTk4|um{fKIC zg|Bkd9k~M}`x;<%gjp`V53S`R9htR*FH@GJ$-<<<=Qo!`(CrYVP>-DD&D1UNa5;B@7L#D!xifRR%zF>n^`+BrI7cYM51yq|Nsb2HbZ{ zi=om2%n!NebPjt-uj<*uL#r+a;KHD|yLAL1hgU;RXOUDb-30EemIT|M`v>rvV~y5F z26^$-U2J><*n7v^(AWUl4{IO>f?eG9Wk7Z(dq%{F^-s88U#PyD8G<7AbE%Q&hKWPT zF0CL$^yKyXVaPi{la+Y?6iI#Qe2f+{V>Q3AVMH@U8xLkEVCodZV9obA#`Uo2=5=97 z5ZQoSn$d|W_IbEgSsI7rj*>#GeKdDphX*GWS|uc9B@-ZFVh*Cmxh4M`hV#VX6UAD_ zE2G_80rlCP&xm%Sl%*_Fy)8(-JJ1d_#(KJHE&#?ckwHNjBPUAyU1QFcu-NTUq3Je} zIQDE)bKVb-lVLbHHhP6ItkuYTRNu>(xbOBHU72@@pjE)zu(TScqo-)+K-i>L~zJ>LJgQC?#1~XgmV3tC9iBLFNtxDt( z={bs+om_v3P_$Ei=Ikawy^WR1UW2+Y6&EeUF`2vt-#8n}k663^K(b+u5ZEt6iUl?u za>rINW8g-{^ia~X*KBj!CRZAr8j4QW!aCG=rOTII(+}7Gq)oHj0n3qMg~GjAVQ2lH zZPnJ^==4gg193Ovob|QurpYd&pR?CEbNmUdn!vf$s~aE3s9LU1UL-?B&E&S|Yw!5f z&TTKHx+7Z9O?>do3|>oZe8{#N1aNWbs9tip+fxN3_En8^sLX%h)%|0r^#-3`A)t=m znbV3A+m?VhCDWwPw7b1a>z7$qu&=oZJf!zWUu(+)G-OqeF5P-@g#hWtx zNHnvQe@I#-7aPTW6@2gpc2sLNOdKpKFJ6?t>a8qL{T=O@yM@A>O`I826=V~ z;;ar@)CZGgIH}^GYJzjbzagdH zNA7w#IzFCP4T^qzOWz8Ud8@Qx`vKPBWtjxcE2p-8ym}BG+w93%fWE@a5MN=$nAp_= zus0Xq)E&s|n`&uhuJ1-p@ij$(I?)ns?<+*ELf3Cx;T#03r~yq(-7%K#%+UjMzCxA9 z7z0}$j?e|}DXNpH_Z3U|0pa|AExS7v+ z(4?wJ!_fvvwiD+pSW=jB>(cLpGf)0TwN2>j>;2rbji}1ufIjPV4X&IObY;Hj0MfSN z^=l=*V}}xClgQ9BnNnfasay}(pZJYZ*Q*!6OLNZd`pM8|A3T4`=8b*d!$1LGI)zmS zX5aJ8zG%}xdz-5A7jI@y)mh!2LTk;ugu)G228HZo{^&&f->LjsBDqKXfW9&bou6EN zX84zq5b>S}^|n{@FF7hTI#WEMaZoK%^7D1?syU3bL~@|>(;wP6C$1VmzAj(NE+Q;D zTvh^o4>(;t3D-?~dvT3Wt>%$xmDU3sH_fo?=rpXZkS%(#gEaVsk8571f-TI6=MJk{ z@zG|>_j2rmfOpBviRD@qWZ2g+AxvChT>kVznG_bYlg$GX`0G~pV&d8GNL#8@o{p2_ zhwD5GuiC&#NFe(x$0b_3$3YiA)>BmsS>tpP2shyQmj{BvEDEpM!FH3frGJbo%zOz8 zJHIx!JUkRn_nEzZw$V-@Jo@{|6BlsjST-*W8V^LBPd9F?Zk=xQK~9RRbYI(x-FU%{ z74a;{_J`ef)-19!wd!l)FB0V@q*giKk5xh1G6zl zbs-%cvyIzoeWmRNQoWl9aI(58cy4+`slzzIY^<~8O*`~h$qsChDu0uVC`@1rIN3R@ z)Dcs3tR(6;xf{>lb-RI z=xyXWW}rIt!)u3b|BqK)l6LL^GsqqvwcPYBDf*uZz+Jnp%g+3Gj4CEB-KJ`|_U*oEdOZkBRZ4%Ou@xv0 zBEX)NJ&$*NsSvo%8DMSqen;2mJPbi z9QTi|0EQWS<^wjoP;WdZ>$vb`#SyvE6KX5%wr>=o&F~@CV&UOJZMSo79p-k?`eMln z^+=Df2<2!!@A9;8TSPitfrCuGsBeL6cHaztMi@__@@rsdyub;CmfjEk&S^=aE~!Uc z;kUiQZPuu^As}b-^`u6}>fplBd&lnq=?^1mH(QE>%R(~kH*S)jvC&)80Ms+GL-h=& zCXRLwx|Ya%{wz0hux3qe=E*H#9aYH&V)k-h0q&rOGfL`&*K$pE;5YRe>Juf9oFPzj z%c^VVA5;p);P%v9Nc7Mp==OIJokks6om0Ndy)SYzhmW_|ChVLd8zAL9-Ch}C`=KmxX zNtnQe(@D?IqxWn@W_=1|@Wf&r_RRL`>V%nlT+VxG=B_RKe0t&@L180SM_gs7na_eu zhAb~$8C)`}I<#4S20}Rn;BIfjgX+rK^<%XQXA|!cJE&|SCsyS{v?SYum)S_a8)h1{%yLMP#=b6uX_wdMRw>u)3ls&(n&Q1OGowDUA|OO>=n6Hf3Qkfc?Tr=#7K*ft3H*Mv z`;vXQ{`bfX^ZcisVxMH;$zr`fXkU0I{5(cnNh`I_w#zgHpc|qNA!b$9%kbM#8eLi1 zn+M{qUyv{tThy4h94VQ9a=CKz!W|$zHf0jssU=W&6ZKUm1+x5>NL#8a_q$v494Sb< zP&?Dnv&g7;95pfEgxW5Lv*RQ|MN-AmBhePkH*)KAg(%0XU#_TCeB4#uNr}hm6d0jG zwDV1cudN5OsM7ar_$^f9g(+0I`!sqpA=l|Vk=0Y%LD~Z1to;&ViuZIISF93Be|9|N z^RPrtr5)YG9M=eq`8>zF!mI%+r)&;)e_kf^q9?K>`6uRz>SYtMYzt%%`QdV}G-GF1 za$Y>Jf(CRiLU}j#z3N@tF0&Dcd{FX;E;sY@xgmLC!9VO=!L)Yo1g7v-5t?Xk)im+xjGPCFo3Y~jIEOk_kW8^DXI!ov4P0_jr;qx=l+3&l&M2pahA zTWYNS#=x)t!kEGKcds-amhVjp=V$yXdyK~8bD;>@#>-mtfQle%qE{>H9!lH&*DV|S z4Z&r3?~Dg>JIHyYuqTNw+csN`zDCM4&c!|P_Yv^)5a-2$_ObZIpQ&V4;>ON~|4+XT zpcOB@xV9Udzo(s6a{Xg(hI5=ZazRJk`>;a0dc%i-J2ItS*1rCM0(tc&J;&Wmn;HM{ zYSmT)^nI93-8*F1SiSUla#STF@52C6z%|r~yPWCM6fmlq8?M|ZYaKOKzOd|8i2VO! zG)1875kx=jWR#{@SG!9oJ4rK(E;c(p(A4Kz@YL*>G;_HtdbedTYs6sJ4wFTCGzhkFUHENGPmo0*3uSOb4cx}C4^bSqwl?^=Nl6<_T99MY?uE& zxE3P;W*zw~eN1*uY1n1dGV-)E*zdU5?pDoz?^u3)SL~j=zJZZg#`yl)cWtMgj)IN# zyYv;_2ltAvm>P<4J~L;|QW94WMAp{ZmYuf9j}=thEqid^zda@j-Mn@w4!X2m0%hdBw%DFiMFgz) z0O{S0Q5T^|UO-wcHVVlUwN2plLkrAB|JkvU#`?w8L2S4>o7C>=yk-vfqd>vbLF1L- zAI8A6A|~pN4RT1uT!c3A0n_S8Q)XuEzcp9_6}Nx-4^v@3aQ+Ju!oo1d99CaD{1x8A QdMQ3)W_hsK)cxXr0Wgk^`~Uy| diff --git a/docs/reference/figures/README-unnamed-chunk-14-1.png b/docs/reference/figures/README-unnamed-chunk-14-1.png index 1541450a1ab6fb11b664887d84b51e1e7fb8c490..4e63f6c096d9ce339a019afae0d00db99f7554db 100644 GIT binary patch literal 8250 zcmc(Fc|4SD+y8aV3}as+yHSW1TM3i(j!Lq1x7e4l#9hRwFk-l*v?xSbvZYO>ED0lw z%UzZdl4NH@Au?kf`@Gjk&;2~ld;gx#`}_U=c;_?o`JC5v9_MkM=l3|i$8j!m^{|yG zKaUg-008`E2lpKT05l8$AjE}6-W+dU8b=-k4qF`Ck30f^E&zA|y77QliKL{Fq@=r~ zQKY1MmXwr~rluxRrwi!n0$#fD-AJ90q`Q%kyHRANCHv!NX&GzjmuneYZ*On!zSh!R z(%l^qayN3%>hxOE>s3f;tYvI0 zTOv2Mo)xbP#Fxa!cbCL>v$EVJvpPF-{Va2lma&$(LqoYktP(GvgvF|iG^*_EtdF&< zXR$i7+&e2HJ3A3s_kdMKC0LUbTa z*&-|!Yu}6eF$f;q#Dn%{0YKn2`v=zfKlTIw+()y0dykRQfAsQ3^wb7&5bJHyQvBMm zr=_gAw%#6p?5nUX+d=P~@eX0<Y2wR`JyB&TdhnjAHxKaQCIrZ+ z5O$B;>L#Jov5)Evxk)SuX}reBpuAEb6PQng$d%%T5&$(VU=Wnkh#r6<2S{#(nC0l5 z??5bH16aq)12#lM)j}F#s3d+cYBNCvqauysG=c57?6fg_(NLjKfiIT?ALy}}&?%%L zfvPozUu-65V6I<+PME+1E(zXs1xd5OYJS9gPSt7@J#}QCGEL5KK~Xxl-8KB4U;R0h z8t3?|I%_4dYVkQei*8r7t!wNnbKHus3=?A$8QSwNlKvaVXAK6S8e9-N_%9@eDFWuX| z9ZW{gEjDE3N^2l34#ED^-i%A1=Xe~HX=h*;)|ad{{*Lp&*mW!T+Id#`N-KLv~9V!{MHCTw@3+I%FxDPgR+nHfr3^gQm3(HBB-cCU5jG;uzJ&-eI2 zL=qWQZ78k0)g$v#Eibo(uQK6TElpy6k-;d^qzd}luFEv#z1gK;0k47^3XYj#_}ZVS zEB+iCDoORfm2#jQVUdCe2&2*_489KhzOEPXwExRe?PfwE^8)w(%(8+?|6Iaz$m|yeb&!a`2Fef2VWk;)AIw-Z{)uD&H~tzoMRxsx zpJc27E@2PCM`SNyD6=U8QUvw5JvM1kk)Q|X8F=!-moP#7J`pCo%DmfyHV+j2h@-Sg zf`S)+(nkgr>yDaOf9OHueq4iA(29g7v0M^)M_t$LfKhQ=FQK5iMR@`X;I@0JfQlab zB)&nT!wNAwhIX^*Z@4fhKzzy4XD~QgFynHYK=Z(tyI`{1gz`@f_iqDvhKc)NZt(TL zLdFHE8Fdh}{%OugKwY2xU!|LJLj@#a$ECtfXSt_4%0)( zHAAueAM9)9`ip%_%N2`fre{Mwra$|B9l{#_G&sW7%;lY-+~hTK1|?Z0Je=L%F_c}8$NLhzjd!WIW6T-LLKqZr!% zld<_?BO|gJiuZ6T5l{!%q}-{;n=uhiw~8Am`p{csw{D~dqKqG{ubZ2n_hBb^{+h

`xU?N=53OXENJ4?ox-w(d#0X-R~?Y5AY|bOr4<_dksj zrVw^z193XH`dn8xw_U5c!arO~7~2!3@XO?%dGNL0_@4X?2mUhqUnO23FRSM7h@ndR z>=9iDhK5WP#>(*VNodBSW5q=$B2Vy^VL7|Yu(hX>04I@OV|M_mB1rVKt51G~=+S<3 zgxr2P7G7KHS2#CO&+a$6eh>Tls0YNXT#tA+m0*s6I2Qqn8>m>1={+Xje=)OBh*tXh zIV18mgBjyH5VzeUB(&dgf#qEGwcH<13pObypp*!yY=>WLHWs9w%R`~8`CoRj6lZilez!94bEIfcrmj`1DZJCJu^vZKc|Ix2OjGj-RK>1E zCIQ_e{2iMbCuONSc7vsB`JCJU1^dmStb(gLx5(IEkBj6cY5Yu|MW1jK(VuvbO)-n^mmqL;>2-3o`p$NSKWsmVC1IsHbLqAN_`=zRGc%mrt{A)hk9W2U8JBM8! zs3Fsg1yDFXBq+N9v5rto#CFv3>Hr!Q^g9yKLNy?ti^bq7<4bjUXuqfW1FQ**QE5d2 znk9WF!#e8#^Rf}IQ`WO9!n9>uA$D%@?f-ofloN$k!xb8@<)wcaBNl}(KN8UXNPxJw zp=2JNHgQ$A3=a3N^a}s5Q3)nuZBY=zo98$}SF0%qW`q zW#MqnHby_y3R0H404Lur2O!u`Ff4$Spg^6!xBVoT$%Wu{viuYrFzIqgmY7;UZ;-TQ z4u0OV{4H4d?T_i{p>Jc;5BAvK+yzm^n^nw|oN<&?vYWyDbN@6U`6Nu)shj}(y5eZ%Jw zs9e*1FgT~~vv)Q(UIreG?0qwnhDkr0asFbK82q+-e9O7H)YV;ybhMuMv)+dR>uGu> z>&^M(O8(_dk^*e^DtCul=A$!HPWLoa-S{N#14maKEH@uyCm|C~Uqw6-jxHzsP`km4 z#>miFg-aoWKdKP!DQTbE0_Bcv&2C1EFb;D`_}cQ~ujPPg<5=-MS8un-t{l!bpeg0fQb*DI>GA6X{T3cJq;F{D7#w&KN_%b#=Q zM2JUT@o_RdzX4W!hKUc5DZ6Pv(`m)WeYU=(FCAOH1rxd(Fl+qP8xk z@vb`@?Z|xdt&#^$%qe7>j1MFRY$T}D#^n}H$BC!yy$#)1@xw4?o(HbgVu?K%nYxX% zsiYX*Mj8(g$E!rU%VPq{)HuUS`O;Y5iBgPyNk~457#KWmZ!WaYZ4-1q0B&p;Mv^FZ z2h|mXe>d?FKQmougV8SxjH5=oU}+Q*{QemNuLc*~qDZR^se2A)&TnEJ0E7FSFD7rg zVmu4QI93ZW4hq3-#JA1G-f+Nh?|u^Zg{`0XDf2s|W;;)Vd2;Qikx6FYOL#~#9+$eD zrNp;moj%B5oca}YxG8VArew*Ya}GdUgqj`OD^?%{O}IGYT?x>R!}58~J3MRyaFyR3 zk-9$tKswJ56}@ZT=Hoj+aMS0BkCej(v;j12POR)b-*yb&>Avy-2P2y2ki|=Otm^_z zYa!+Z9W>=47gE*Ypi0|blas(DG10C{C_>Z15}iCb+fgLOzLy)|8<+VURSBz-Ig@Q;^#5G^t%X!mxG>D6XFn3z@4|#%d*3Uj+uh*RqDMCTXeG_= z{~@60b+{eZ(U(7GKpb?n;k9!y`3?S~Y61{+G2TNp+n}-A$@r`#93#4&!>Monshx1W zlB9X_CS+Zj2@3Yotk|<}gmTjk%cl*M0L_!-=V5!Uj(e`{P_|_Pls&K^-DeBrt#U%> z(lyZ0H=_lSnKp-&T3je}=c7e5%xi)Crx&ch+@vPybO1lE( zT>G}z%SG#}`NT#}osiaM-;B;KoBh%Y7{p) zoqYgqQGw0&A)>{?Z~)U?fKAqosJyZ+78@mf#(9!>98t5CM+4-3`>)W}5 zqsPUHEjkNkmifF56-uwx0-KEC12Nfan_X2%*0(jJ1Vzz8vmNp&T0gi!W*PJNzT_vc zBh|J1z6r|a=wm)cMfc0EJFY>6zMRCtpM4`sHl_eY{}3^04y|p z$W?lu<1pBf?V8)Fhk5ph9hPRn*vWNu1@~3YhMnb>Eocz`V6Tm$6sEH~&l6>g%o)Rt z^T-x*p27G>22!AH0?WI0wexor_J1C`T?6(T>Y4jdb<6b(%+#ol!?%KzpN#pnDGbYa zLJmB$w5-3r0aYzZ6e=K=7mcB@{`E!?AQyX)%4x;?6G=ML83DRxNpn&0&Ysz(^6GbF z3H_R9ZW}apx82#`WRIoRHI%tlSHId71~Yp{1J00-fqCWXU;fzXqnS;vJLUBFq~Q3u zR+&zBR6p7r5KHkP2!<`q;LYaPJJqkd%DKUY5_4+#s4rN77GsP_G8zXHXR3H*(Cy!h z*+kxm#w}ZH08100oZxhn+VVH``y9>2;+qN8LG6y;Xi6sV4-!XMK+4yB^JJPlMx}jX zT=bmxi9ODME0o*(@!LzOqHn-w{D@JS1rrP9wDfu<53b{EHiST53w zt=$uzVDj*C34lLBN0^Bv)a`iu{)4Pn6o>CDtTR~MB|uAoWJ}6f(D~c9!KP?4QR3?5 zb0W~-TP&>@;mJT3$mCX{MOMEPIH8T|AO&%eo;ttx2PxhS+#*Cnix<#fBcPoni^@Ee zjMjho^a94|dV&m1zQqP`O>Ghct&1BLZ>Zw#mbt;vIl&jOzEnYkuybE1=GwV{Q|tz? z{Mx9DFyoJM#@Mj?@od2l$N20F4-9*^gQV|!8M8HqkxeZPH4mzg_`iPoqtR*Cb7*kt zB{gRMto?{2W0#&71H2;zs_(KQM`;DW6-75$zY^@30UeHIT?z*qMQZ>Cv@gK0Sz*Xld}w+g@`h0w5|41j0d#iAkIiY ziAJ5Hgkewag(t&Rg&9M3WsbWX;YM%ym@Q*V+k3^n;^I}(qx!%^D_xg`tEu}-Ep&WF z38Rw^^L#dYY?{1W7jP6OrzP)}MYWli*jTYeqNRuGr84Xsx}!6J10?O!Ksq+>lMD1z zEwh6jf}xXoTS$SB^aW|tH<+GwwK*vAYVi;O;j-h4-+*@Sc&IL+cdEMefPF?C4|hRw z^7aM0uT|l+nQxsIMDb$SuPQVox3X*4gfo7)R^q3r>O`;p}~qb{H$8e+-11q{T5y!VqmT;YNEEAd&IqWh?P+h z@N@$lzEayJa@AYSo(N5!PMlW34G zB0fAq_MR%k>jI&3j^qz-z6*~Q-`ND(sy|Rb+*T2@w##dK1`wx{u81I&n?Rpwp9t}dA#g)K z*Li(|U>9d2TF|25Fcf$7w`@GEcObvS(==6#cPk&|6~XbiZnEcTLB`C&vyJG_rAv*P z^RMbP3FOypAxRTrpg4WUXf!)P-v_RaJiX-q%hpN{) zUzyx?bYffGsVFE|Qa9v>i$>Icb)Kg=XGI#e#Ntj-W1?#fM51rbk^X|*{!M7k$5dTtm2U9=V z++SdRZgimK#B3Y@Fds;yjF*o_rZ&kV;pY9&YQd%t<1Rv~Ac;Aw_hxf#rL83qkhVIb@$QQA`HL!uJ>b&j=ESu8=Cs#NH*D)iuS{??*X3jzr1=q+PNvd1dhq2oV}+|}cP`V% z8p22K47T%_Vh>|4b*P(&nZYW-suezZEXtQm|!(?{0h&428bh zof^J!Qgw1~GbxpSSt;6|h)%O$UWj}A?6vEx=ut?pV28k+kZ38 z)vFMjXU(~6{P-|?(YStmd4|sdMs8BpPMm2fGm|v?KYni4kuy|$h+GPUY{}P-L9O!M z$kUL!7xr=EU&j#tKcyS0gjx+54GkG}l_01))7!Fg-{rkCBadABTVwdITc$V0yoz;L z5h}QV>zx@3J&$rhSopd9ep)!$+8T~IHzca@TEpdEFOD`7qPg*65VLe-Yn4w=(QEk( z0e_qvF7y3JLAOUh)cItCz6@M=koHdDPyVkyWM#cx33{YnWoI3i_ns=5QaUwWl$IRx zT1ENujwh|c_@}+GzL$d^XZ3C!gqU%s5S51IpmomcaitU;*TxSmqdVW))TBRu+C7+? zxEx!qSwz=y9Eo{p@~ppU=D5!MjJ5R(?v~q~qlCG4Mb>Bg@&~GACb(qUGGpK2RG!b| zeRIDR055ePOX}C6J6_Fc9SGFU&G#aW&7`c>EH=`=0Ta%DU^|&<{Z!_jE`RUEshF=b z%}pPEDfH+ovg`@k5~zzZEzY^$^}V>@~HN0rT< z{$2~x<=PI1w}%CxGp#P01!y`hoeQnNX9n~2j81-=G&6jrc~rXmiw1VQRMUVIUz}wa z!Lck^u}e=@@hqbC`*dbRv9no5tgoQ6UT)pTL%~|~1NOr>gpAKD@A``kgF(~}z2U72 zcaZyv>hSVbaaYD@QL8*_vG=+Gv9o&ows+QI#qecn$&M%8!rD&_WILbtRfo^re>N$B zXEp3m&cMYaUezqgoiJwgans0EU*+rz5bSAL^WR<-;SY|T(+K_AnPaG)|Kj%P@&sMk zPHO#lUSC-q+?GCmpz|Bz?h_Zs0|LYrN#E>8Y%Uk+L|VnY?J2$a%7A)iYQy<{|Jj_4 zW60Df`5o3-a(~cv6tzYgRn3M}EvdD8m9Z60ME+x9or`4JxZs6WW+YSA#-Ca?wliSH zka#xxd7q@rICD|MXaI#^PF;4_s%(7c8>|t546P&~2zy<#kNP|i6MM?qbg9jUWyqjr zt!$k>q=I-d=Dh}Z8In|@{o6e$$$VkdpE_y4d%atTkz@?#@`C=`Fcq(a7knXj{$T)G^`Qx}SyM)owP z#Y99Ls*a3?Y%{KVFS1xmh4WO-dFK+Qr5kIeRGr+csHuRo^UTjKsw7USku$>b1?ErR z8s0la7VTvkq+Ji>BK3waCs~i8`;97?mOg^%MQT&4Zs8j5FG`?Xg8NvDMjzYj=~Ed* z55e@+T$|bWI+)#ECcQ|!KJ{~!VIx=-F~5oIMy96;svr3@+C1&yMy(X}a0W7ftIS}S0w3e-Mr$abyV-fLO%Aj@f!R8Pni9>TV4}l2{tShUDg#< SyUYF$fZ2YleR;-&3;zWbG){;B literal 8114 zcmcI}2{@GR*S{V^!-TOf*~uRDRZNnJL=s9RvQ31P7`trK6NO|6rK|}l!jx>;9#KRH zS;jV#$j(F=`}`lRzyI&|eXsX>UGMwOWv=U)`#$%%&;2>)e9rwG5yy@7*;s{G85kJY z436lUGB7Yg3=CVYA{pV4(`~EM@DI-8$4(uFe=;y=FfceWXv8u&mIw-B1qJN{v7n$` zx{#2Ns;VkHr@^41!QiM7+YirS1?{j{J1kgbwE1y0GKw~G%`=K_w6L(S`@Ygn(9RAP zvcuY?lN`V6bu1}yETO}qu12n|M$xXW@L!&5UZZQ?2#CdkAgqiO{asPCQFQdCL|$|w zJywGuwj?&Tza+Mwo^B_YP9o*G8s%-SjLwTT${QKU8=;ptGL+EiRUo#CL~4vSYNXRi z>2{sY}tK2@qm1^~_`)@U5(&>jT1bG|g$9cIuIjUeV7L?}K_hn*#wnD0M zQ1N+FoklzR&#nv-^#_p08n1;wp-Pzl?;(?C=@#-(nEWBa#IPLIQ?f}=Ksx5vSq)yk{4 z$aZ~UD`9c9yuah!O~gi?ZQCHi@FC9YQpoJBJ1U;Dg@j0Mzfa3Tse?KSkTJU7-1!+>=6#Zo-D^S05a{2jE!@)6iI0`ugsdBz40S0IQ-q${9tOmM}JUmr{+489%OQ6K23y55(vAF26hKw4M21#@Zb zxEa6B5bNRh(g^lu31ZwqoJc5iZXd#o3;8qzjZq4Njsl=J8VDOjn;VD|*>xyHVL*`u z(XEXXzIzC_1B5h?V4Vl@sWjSrYdLGB)Wekab1zfcZF6@qhK=3>?;uqFn?(|Mb9Pt^ zm=bpbWFssYWhiKK9o>#?f4{{a>Hh_L|BsZJ0J1EAsBDHEqbHRO6ipajV z+LZF`&l5KNdj!!cL-{(nu%h&Co6^WnPsbrsX14I}o8Br#jQtTEM2gVDRN%=+z)mhVRVZBb)<|amXtYQY%Bz;Ask+Ua$#=OlOVz%^1<=Box7J0!S_?a zGEc8Gqc{3D60>}my$8ac{~7oVO{^_Qo0r6g{4k`k0(zoN5&D5ztUhRT4chq=fWHZ) zHPnFTfgOMRHVDvc&{2>p>^*a9q{J@3@ZWRfI0IP0riUx_Kq@^|4n_CMP-|@w1@1pG z<8K_wqeE%6Lr$BqgV(#RqX8*m!oR1}??n&LfUoEqo!`3sj^kF<~<0j9m;?1K4%?)-jPKo~pndcZTx%QtOT-+8YT{{e0BLTC^b#BZcT_w?MW$Zr* zLhb?NJk{;<54HX=#gx+hY0;W?rB4iZuD=Yq+ME}jfaxb8$SAVn1A8!5~On)39jb0jE zUS4i6>6{81^R-Q_&d(j#iLd(WSJc1Q4XUU-$GTOQAF&PgSSgXb+YlhN3)s&Imxdcm z)mPS}SKZe~JHJkrnwSikW!X|bEeQS!%8bh_m1Hk1Y|W9vV46bJhNS1RQgcC?DCoPI zKQ(y66YV#9LurT8*F*(1kUUiJt?P^1nRKiHt%#ajaR>6Ai~NdA`JnQ~j@d`CVcTXx zbz>O`PVZK+-Smdh*Y3;(HTg)u_a`NIU)HIqa(b8g$yP;0m;00-fi zC%#0Vt%Bx61wF9=Uo?`hUDn1){3MTcIj_Xh{8HWt8b5dv!N(2bF%8g6YUF{)|6m16 z3jd>~nORc>kmL0)uV#3o0`GyV*C3sr1f)XV%7;t|%HINtz(sg-=VEAeLcJnXDzY$$ z@>Bzek^7#GjZcopA;rgR`Dxp3!KGveoE1O(TF-^hs=jrKVyF6Y?jrjSu9W1%*oGBh zAM!6%^iM?6stF%w_L%xuvLk#{zAL+s0)^f2UzBl}z(jqWZNE(5_Za!c*<1C!!ep=^ z*A|CF5R03oBm4|ya^=wyQ!s73b;{D#VS7jF36{XJVv%2wGIFYsN6{lrfV%>kFcC}Y zf9(raJmkaFhztFK(nJ1yE^qlB%r%+O`@?&;4~J?`~$lBzYwkJa+$p5bnjOg;=2G; z#5ltEsc&H9S&DCr4ziNY`zs}4O>t9;*%9WHAwx`(o1!EaMjO}jgQ>HwHyt`{21ZV; ze&R8w2jOl*5!%0))(hz8pw=r9#gR%Q@clU?Zy2S4RB|0hvV%L|D-;|^=lPee z@21(xDQ71fUXuqQZqB#pT}G;t5k`~()hy@WAVL;Epl#53$a2%$0MurPpt#Ka&3$f! znA3_**4`}Smc;YJu_gHzxL@VO(NZ~TYp?4znM<2{RD_SHBT9y* zmCfVNK6BNnX;jN7n)~X@3|@z+D54B8Tdjv_q3<&VJ9|jZ6j6dTlkxfO?MO|Q6Zzq< zrv^WE)wf(=2EE_!$+^w0jZ6G-{z_yD{T$l7?>(n5(9ePpy?Ud@fj*Ad|SaV`Wne6Cxc;YMmH}{%O`(0+C~Gb()F5p_#0Vt#RIZEqlRtG6CS)0Dt8r3%E1pnHu+VgxR*rer-T zKdHVTdObRL_P^@q>3!Ti@q&JByo*J%tQ~}Ek>$|!TBMRiYUj|DJg1*@ z7&P$&k`dv;tUAMJaAMjdnT+4*F{!=u_74GKZ0(p}fe1+E4oiw!b>CGeb-gQKcf(Pp zS<;+7d>VOthsL}Mf0?_VqbiWNpN<@;s7ER(k6BVadP+BV9xFCEi_oLkqs?tvi=bI;?#4412>?Ed9>Gk(Gc{e{mw?XcyT-p%R7}_H(7ZA3{cZRV+1}*;`f>& z2jWGqe;3l1#QKe6OsDuUCk9xiq>-BMH$z!0V(mqN|LC61yBH8`;GOZtWgJ8wp*5UwZg6p%tD)T`B5&R*?q=hu)_E zF)NLxutzwbM?RR5uRc(+1=L@Tgo4;aU?<*$BpshZi$+Hgv#fw35+UmNX1Tvur|=AA zcfT+zJJ_j%`*t|t^vyY#+WT1#ltO3HQ`2p;CYl!(7Gl!VD)Yk^=6%-ABG?a#Q_-p( zj2ecw09(Qr>8h1A1OF@Mu2^!8jlvKya=!L+E&zhq5{8>K2mvsOXYqDA2OrX zz^4M5<4~H=OFS6aW?*Z?b+8Ud+^8;JRj69-ra4Q)apcISf?qXoQbnbt(p%VVX9&Sw z81#(E_nTrWRT4ip7;|8+monHeeQwVxd}qjqt1*8R6zzL+U-U^n#NDr56E~Wgm=pF8T z9Z2Mm*2OiPsu2bmxzg5*R>V*hIMCX)no_$$iO2eP}5UewZ|eCZWpGS1B#fX6bQO7Zx8puJdmt_{BC zcsOGt?(32b$jZZcBg3+T==fclm8l^5y8reED{|K@5}q)r1*<%CIUfQTrb`c29h~ib ziD->^a7x{trgre#3DzRHcvdSLJf{M3ng~xeAzBM}sc7(m`|sW6_$oRrx6C|x5z$(@ zjP_GHQGgv{ezf;!Eg#d5w!>#0%!lreHxcaRTk`uJSVce<+khC^`RHxlcwX?e$enb? zgSZx;++L9z(>u{pki2EpBEf)Z@rsmEfWF}6_!q|mfk;u^{vSrEAvNm}LyWB= zgP9DGtoYLyyxo!h|InT94XIY1YLiyrE z4{+ewH>Bq4*1=ay#*fA`pQo3|QjL3s5G7pP`N_s6+5IWdt3w)>xc$bpUksu~q9usR zD1JXl!u2pOGZ`wyQ>lW#QN%I?YADkMiBIfs5bB~ThU?REv zGzjYQL>E>)%1KGW(Z(##W*w%q_qOd#l?q@{R?gH7FW8VSq%aPnaSiHloo{z31m-4L zi0>*=Ne%?%15VxSsJqX}kB(l|)=18(k@==R^#~W?!)zQTDwFybV`&6BW_dznN6|Y$ zzq25dPre)X%Q8{yz75uoFq%_nOm8A`X1Aj|;G{a774*Rq75}A~uL$~9%m^N?9|3ZZ zO5<4{t0>;AR{N7fAxTN-^1#+NPoB9d>rx`HBXRQ+J3l22Zh5^EpMn3P-jh9)04?ug z7Kom8zbI^PeG@vlc)ZAlZ6!CNw!N2H|I&RBn?H4pGH}ck94YhyNhV8amLPUNwT}jCrwq_V7KWZx`Of1CMpP5^y}Kwh}kwx zpeNw%6O#4&zYxYz-XA&3jE^ntxM6ndt01BXrAd(8x3_}(Y&i5}b4u@BnWIcsxio^$ zEAC{`{)}v-%leWIMd@`5eJq6#Kj}eUVv7D~-s}f5T@YPi3siq5vR@1NVr0w6UYdjfe zZ(Ly2{C4>LjEx83xyRNQ(qNvOEH5Em`5QAr4HWCxv0bE{?UWM_W?;dLx9Y==J~f1L zOnss?2#wM-=Zi+-Z)wxVw~@hfHsIJ$OXG8Iu!!{8&tmB_lX09}NqQdAb?iVUR3s+b z*x}IUgHYz^ZA|nZ1-1s>daNbtD#Cr7QXY_l1jOb(fY7qTFei)VMk{wMBD$U&J(IwO zw8~m@=s*$n^l~w&6>CLEd|?%is#!GdOM)V(uhLsb({7>^iOydQwh1PMjB#S_({-=x zY~|-I(55K7nqt;S3{X(cf_CR8gu(4iE9>jfXnW-#blEB3E3Z$jyCi7v836?bf9%UCaC^eP4=dws+-$STQJ$TMy_DLBI(4Q>yp825U zf&}fS6q-xz)#dRTUY}@bussVZcz11^J^XGj;ueSA1hRwn3q6Z7y)7&$w7S`)=^NGM ze%?TGxk{5*Q}|q!KR>3!=nnTM*f%92my_D=+dP9t&+L=5-f`&}`QYIfL@lPYCFi$I za$-k2x9qf`1ctpz5aN%ztyH?R8I2j*JAd!RCmF|1l(47Ugy^P-1x4~SBWpwqU7xfM z#E^l8N*y(O)-V$-lN?@5#FDvoN7+OBPC(37WNet1)e7@0CF_lc{)RJC^g6D#ox+B) zGP$f=VCLgIh8Wglb+NM*F)Z!s0EBo!yAOkee8aWOb+8Nxg(z?ldYf0}f(=x@#0ofLEc3&yI>H87N?Y|)_T>470IVL#Bd<=gBKp1&-Z;_*4#`gi$Qlv$x2%p6v^sK}ed*4? z(=K}me%UTwA09wwY69-q!1t^G!@Gkr9kRY>l6cvH(q`jw@e^m z0a8gq(wP}hKIcl+pn(u_YA zno*aKQRC&#UG~WyIytP&Z@-}V0p72DPc8ZRu?}fOTE30+XWfol7m}_*mv+L_=-vQs zwH>=8&UcINl0dNo+|Mb$A27(k%ES80x$bSzUq%bG_G4?PX9~*Ht7`n!$ti@d=P5>U z9;>g{mY$p}H;bBYxW~nAf8o-baM1H@28&8tkrjr+I4`!^&_pq6+#*0CfHujbrIniz z*cIqgH|{~&xvNmi_jJ~2Y@iuadTmE2RId4C-}yEoU=c!#LX`8jWa zI`58t{76WZYsW*fE_Rl{8E2zo8{g&Y$=99`K@O7o?2(#4qxXy)5}+_ zr!VkMj@wj04ciYV7AE;di&`tYRW*31$Iw0NWnMTq)9veJzCN>BE}uMJaQYC|dP(3= z8*GmQ?Y!@mp@YJE@+~w&T?V}yTW+g*SRe`%4#mY|4 zO%_E~dN$luq`3>jR*m}>MemqzA;(y~A)b~J*^5_xQmR3p(7@#Dv<4L9*SZ&H9%x#e zjdEwM4j!)->OHFtUqKBk$f;KrzI3{uw6(}w`glrm|DO4nn1&X};X73Rc?;I$FP+w> zj|^t#ADv}%QTfnb^<{=}Of0k2`HmvHqDGL35wZnp`J$55^d2X#D%lDl(p)_Wi8It6dFsO`?PQu{dbz3k|leCB+QiL%Zk zGyaM8Y8f%3xoxL|%FWo;z1)><&z;9aoy&VOve9iBw16pXasLhyrJGXL*DmT1K+{Oy z+Nvd@nYLByZeDBs6@^6qeH;T;XG>K>Wk-in#u_X)@kzj18DgS8I!E1xvF@V?`|zhS z=bS!5h7Z=TMbuq%$^3Cnrs(T%4L4fBBG8 z5q`zLCN(fmnTm=ZGhRvET{`9}86d zdHOp))27@R^A5gwp(Zn6iaa;i7ULXr;T@U&2osgi^6K@_mw^`PXlMiscellaneous function 

unite_dotwhiskers(data, variable_dot, variable_low, variable_upp,
-  digits_dot = 3, digits_low = 2, digits_upp = 3)
+  digits_dot = 3, digits_low = 2, digits_upp = 3,
+  decimal_to_percent = TRUE)
 
 serosvy_extract_posterior(data, variable)
@@ -157,6 +158,10 @@

Arg digits_upp

digits of upper interval to round

+ + decimal_to_percent +

logical to multiply values by 100. TRUE by default.

+ variable

variable to extract from posterior distributions

diff --git a/man/figures/README-unnamed-chunk-10-1.png b/man/figures/README-unnamed-chunk-10-1.png index 39dfc0d45fbf8f232deb02f4a340cc718db1b5c8..74beebb591e3d1ea83ed4af07d0143d61beeda75 100644 GIT binary patch literal 7503 zcmb_Bc{r5o+vA-XGiEGfFMCudQciMw}?ysF<xVAS9d| zceo)CNC<%tNx&lE6_56Xaro=8>mK)=@HYZshCl=&%rX#xRfdL!@Piq`%nT7|meEV% z|4~S!G!lhQN^9}<_6`()23GOlRSJnhA*E3$@IReGA5Ig%rIFIo_^b4^7G8!KA|nGP z!=qD3bQma&L~o(fhllCIysAJ%6_3}FMrz^ldNa&=c~!kU-Y}gq%;SxXjqwl&7y{mg z_uy3?j}Hq2!XP{zZ@iC+hq)G`I(kPS5Ry&&k4SUay~7BEijC6_8+TgnL|*8j-tE`y zKRQ>{MCLr`buM-?AEOwWXKjB|ddN7^?+1@pwykV!y?>_02@#;fRQ)$6g&$arG@DG6 zYbPQ@zq3mnfmpO5CtUjX$6J@(AYzT9;c6!*4o$x4C7wx^72!0>zD7<;SXJf|k{3Ee zVyccgmTHuTB7luzHr8YjFceAfa>4=t9vh;Z0PV6R5WK9(y$F=2Okx5A0KK9b<(CO4 zlmek=I|96#tf5h!09A=#OV*-9Ks2gK!kWBp4Hm0sBDVPR%d=ATfPlKi&Hb;l!pBsF z(5$uWiOwgEu?w9vLzllEo)QvHya-9fFDyhZafe7t?7p_6|J+N0TzK5>h^YnkQp)^v zO7+Y^jV9H5%&;Yqu84xcBlWG1?%J_5*c0?G?75-VLg@6Ag)T=*+k($%H2I)6YowGD z)Pdv}|0rTcw|lVs?C$aTF4tpEcrCLkhr130z|`32HU z9H00E^p2*%LTA=nbLY?-qW8N`tR00J7zGzPkpDVC!w{q&nVi~)vJ?Y?ulP)T#QKW0 z6^^KRscDsq3RZe?wR#_yRa2j3*^MT(R{~{^SHUO2%D4UwRI8TFsxofQSrJVBj(78Qh15rAfGgQrl2H6{Xm2%aPl^9ENh%7idG8X8Z@`6#F(5o{WX#QgTYfjRojq%o0WGd{UYkZ8gimnfDypdgP<5|)>^)1G)zJI*WJ#iej~}* zudHkVImJb>sC*+al=D{p5$|MZ(dolD&>@daalEH#k|c<>qP0{HQk$<wtM}q$_a6B zY)Jpw%HrBJ?4E%wOPtA*6rkJu&OYOfi5RRNUt>7HU~Wv8&1NvGuxi3M@5K`>yPn9> z*zvszjYvTj(0Byx)6Jf-uw(Vdhk2p`1oFu{&=m)O%#Y6qRSo*yl4%ii#p|c|Dq>>F8vmP?&`n0BU04o zKj`KleKSym0yL1g31!{djCa7fZYne z7|{YZhHm{>4oq!8IcxCEv0a+Rc*HuOg6g`j(0+Ig?EVSnx!~sf1y+pdcHQe@89aX>G@m-5>K&U zS3dQVVh-?AutgPA{DqRfzWmwc!Eos}7!cju>Nr(8Cpj)hod!<%yipj8C3tx$aH->P z_Le_#D;ueaC*CYPh6f!)wWv%)iW+zJ+qXIM&D~%ga2rRgjNY{aObNw+L!wmeYoXSA zK+F)PCg8f*cuV=T28;kl**P7mN3kZgVP@r-u?LK4<74n$j0ietamNPVp98AZ?=E0l zcbcdRo$4kuBmxx{XQC&?CoNxx2$+S5l`GQ7I~Lv@MzBY};`X4_)H)2@FN+D>6?%yy znrPR%I;|JFTP7+$jiqkkCjzvCkka2>F)H9y`Bl$2r7s;>%eCL}R}Tw|9APIp6n%BF zF1SL578@|~b=Ov5Av?FG{2AcJ>3?!v%-5xBk1dNO0+&RB4f(WK5g%HX_9#$Bm;w`8 z6==2pd0+ztUcmjBsRIO>B&tiNx<3g|tVKrumVvnjq$|L2A|!Nxl@)xwvpBBe+46*n z6xMc)Tqw)bq{NzKO7JRVMoX<{pOBQ?2w-Cu87tTXE_Bo@$wHVabij=y8Rh3#qEZ)G zYNdB@Aeijpi-eN#%G=p|&Ki(vUiVR$)sm%6iW(YGYU07=uoJ>|#cEPrH_q7 zU;zB9ZjQ%jwPg)CsJczOrWd}fzgwBC21~L-N~LQc!Gap*d)eZrAYq^gCglxl(ESfX zJcCzS7+9AnlKkQn1KKSpyR#BRwS;{NG1}sKqk!jxr%>aPgdVl`(X;B0!q&zMUs6{1 z+b%MgPNK);NC?61Lp7}I45eSjOhu=$%2J`?qquiDyIz!@B+#VzrH!3&R|BiEnV&gA z&}oBVG%#kYv8PiKn6lw_PimU$zW#M^zd5!R14HaVr+p3+2Yr$zj$uF!90v(^pyP`> zS>6I&SIaPhCeUIkbz)A;}qJn!*CJJ^MkjA+9`LXAL`EVT# z1RQ=C--KR8FxQk*cdC@uAA(HB{gDAza9j%>Z5& zET|dyqL3z9k1Uh2hIQ6smH_(&WX)<>wpZwRO%^B(0UpWxo)@qaHsU>CN~&;V@{2WD z3+xT-lj$9}2jI8E4)F(WvGwBoQgK85$kuoB2$a7lD1mhLwV(e{tjKH`l%-Zd(p_$= zrTxe`Nk>r%&RD+ws{yOnn8--Ke^!CGQt0Y4YbV#w_>kXV4sSy7%RlGB3Kp3vj>g1g zc2yjAthd;>duC$)P#aOVy1!kH@xZ{Aj8BsGa&18Jmc3_gtKt_cDrWm;hSH5^nob!P z3gDz^Hn%s_pbw$dMOW44bbeg)6$OX*nh(Nzyrm5-+F;Y0%Qua;-_ODfuU3{mpOs59Mez|qy)A+L#a+Ka}=mPQ5_dEjzB4#00;JW8c&UM zY40HWfb{9z(l#v9M{?ZkhQ+D49fO|*u8gfse5PvgwE`YvOiXjx%J6Y0ni$bH%CO2l z2j?4Ou$U?wW7xPNMm;_rb^|rrhA_!nIa`)43Y;-Uj2O~XvF(!=dEiWgv{#4AZ%bXi zqWRN`pgvknb?2+Dqrsp>w}d?zY;~*>Q3LkZb;Z(jH(-evK?gb_R(syQS1k5La3Ac7 zSMGRx>J&3Q|HKj^k%OvI9NrdZs#exPLk7~$8s&c!yg2qD&##aGK1B612jo!wQMngj z!sz-?WEZZs4vsv$fapfFGc>RkHOsh`3_WMAG3o}pF(tAl%Ui{9PT1OIcOw9A*bN@3 z2edW=K9#flXJ7kA5?<@rqGR`T)@rgT#M|MwB&GSYT*5XWtU^4=b^|bFAEBi8VU4gq zvSpD1kMIAL@EyqU_U{w*S2kIj*xOq->lA6RslUBBW41#uC+4C{)2iJl&bENmm#Tn! zcEV)k!8(=#9lAL`s+#?aG)spXGxWCVFRYUo)8~)lD{LT)!5W`<<86UbE5)rQ3#xk? z;e(6N&Ci%vG1_6FHPy+LoGoTHDsV?H0>dh?Vd?KzMLOk8C{A|8-B4K8nv3i&8vQ{mJo4=EzRe#VFD&Lgf@uV~`B~kfEaNA!G zIjg~DEv|!$vDRKNTijltS!XD;6{{xt8Jms z6wcV(57=etuQ+1w1_3z2kM-D|o$lzxiVp0wUCD@T8K2Bf?}J?f?qmMIv# z_1<`=ur}U=f)*}9m4cbD6r0OvIFFD0eiu%BaQ{HwK>fCScN3K-S>;5^?#=Lpt6IW& z$f;z|V)r+IGr6YV!-=O7f6Bw9+?s57^HUZEOGLtT@e;b!t(I!EHVT}1+AkC&6Kax% zf)ZDTPfIzItJ`XiI0cX+D@1?sIoh4KwcA|4=Vk(pm~2G~CwSd}$54K&Fe>Idlbz#o zT@TeyI=fBS>yauy88=kf8wV|I{N=!5Xy+&tJHs8VwnZ&FSs3!v%^Ogj{0Jz7oYE0$ zmoGwv2ZdE0vSBU%o=#G(L;ecm7%*KcWB3k#=_8hyt1C=O)~t-5E-U)dq?)6={6Ce} z)nsNuefLL4())-Ts@FVCgzrN*Z4xS>FC8ia#9u^GamNGD!r| zdL+Z}xg`ps72){;9xtrwo|c9+C<*g(K_h}CjUBHj$P%LxxP6Me*hx@5zV;*MMUt>W zm>4xnYr=7%Zfh0GyA`)M7+UGnqM!JuJU<@_iqEcEy;ZbJ5G(?>f|h7?%r6 zP2Rt+x3i-F>gz}Uz$g%-S<7(dJTN71L8TEG>6K6{#$WOt!30QIeUjisRp5SY+xW(| z9b+`H!Em+GcUiAZoomXWNOC2@^ik`cvuDN8lxsECWYD~#1elfW;nPk7rrNh(`IylK zLiQpltywqz?8vwS)!_IB6f@~ZzHUxyv&B{G~D|@(hHH z^@k{<20+S)>l7y?LMo*$0rH3Q`eg_iFKt3hi5@68GpTPSf5IvS306BX##~tanJ_!v zc%~0xYfqdc9tg{#5*W6u9>*Ciwg_tWbI-B|CkRjTuT6Xh8vXv zN%u@8R6sRX^?T&uMSqs;ye1WO7_$c*x16K&n~PSt{UlKEQxR;^im{faMQWV{z-?O< zX@MfoUmrf&LSi8=Kx^PslUj?;_i50xTaA-3Nod5DWIloV7EKM?HfZgVVQlBY5Q#Ivp1N*#A@f3;m@7MAaZ7;#mdMp51e3PNJxuazV z&%(Ai10{J(?a6@`&q0yA6Nf*Xl`}$sE?bx%@3w8w-{e3p5P;P~#h(6~gW#hcL8=pR z54-e}&q7M~@xE=Z7ehXkxUp`dEH7=m*qQMhz647=GBx%VUJFNAT9$fV3WrTugr2pC zo(=NUrA7f$FAi3(PluGg;`A1_`B&}&FUm0*@3+BLI|fwOoaj59HDXMcV|anV+5W4* zZI|V|-ng7}DKvuD#6kI+eswgRR7gn-H*WO1-=#qLyo(BqgIxzvU5<7vih<{mJjS=I zTjGiv%U#u#t-DpC%e2a!-GJVYHFhkTSKPO|+uXy`9m@Pvl;5%LTVTiYI>xUb>!dFI zzV(U-c&h13N1Z#2pb4n+hBc{YQPQ+i^K9Jlkew_JhVTog9JDx zkIz0(>LYXno4q~oKvogHse3$hwEcRdmpKT`}M8?YbUx{vP#f=6Ke-WVUc4(ZR`S|ElA}yzj8tO|` z>LzEJdI>$NQL}|=+n-H`(T^Seaelk{Zg`Fmjw%W_ zV*MNaShK2~=YoDe+v7z4SnK@|1KPPt-5h4j%wDPA=Tk}#m;|hrY5S~VPZ@p+$^2S; zk_yWqDNXv`h}}M-UpX8% zCRLulAa3=;@7gw1;0S0{PsA}f{aHq@kSSP(R6bda8pgdbldQo^k?h&``ij}F>#UVG z3WRy{a%cWUn&4%H*KvjUtL@%L^ZxV3q5NwkD?y_?0s*w|O@qD(r-XmMuYiZ}|DRVN z0icKmdJ$!rWoXWp`va?MEuE&p-lb!Ub4}L5*2-MS0OST8eL)dO4kg;M>K~u$Sb#G! zocVP+zu))t^vRQ1AFj&{4))Rw!UBz*QcJPz1xw$ejXc~G_FSNHGo`aM$FnU>lh^>K zJ8({?ZMpBs>+N|`S(-*a#`^;gB!8*2TAZ1n{+X)+Q(+$(Yk!jHB%QC`m(MZQsfFsI zf(A}cw~t`j@KIZI6Zc%p^VsrCGd4u2qh@UIM5F0u-RFbNuhfGjvNVaD7x4=*p8=`B zv+)^SJgXb3wk{Vk0;Gt!1-!Yp>@zU^TO2LfyE|#~%?)l=$!aBSnm~Mv{<=R?oTR@m z0H!*&hTR+?BexbZ@$WU6^}CkM%ikeKx=z4Y%d_9N&?~)`G>P?b+a|@s$!^KB4%?5$ ztce0K0ae4?9}VHbW?7oU-KS?Rbf}1=X21RubdpZ9`hL(YxtD%6Vj1_+1mEYAv(s{^ z2EJR>-q})ynY=(fI@7c+f=_IGF;C&v8`jo$gA1RRrP` z+wf*8c(DlEeEIN;!v@0|n-R$XD6&=be|-!xtd@U_JhTEy182( zP(%U&ga%%7{QiRF`QD2IJXLVF2xxF_z$&JKcOH6~04cthDmY#F>-Sz?)Fuq7UW33O z61a)htnb$FB{kTJ1P!pLO;Z2xiGTXXzkgrnVODu!eyWqk8u>3qopw6!D77Ug{tt?+ BDP;fv literal 7776 zcmcgxc|6qJ_t&SH8YBi|Dmx(*CE1rKX`{$Kw&zJ?X^dUgKB+uCW9gBoY;AbTTDGx& zS}7qF27^zFEo)(H!~8z8)c5)RzJL6FfBarQujF&jJ@=mXJ@=e*&$;*NNlSA9zHNLw zJUjwN4j;7U;aLyy@T|LlS&vX`+gB!#KjkNn+n6DLd3f}AcszOZuJL$Q?%lf=Y3T9j z>G62#UF)-8{SYlIk}QY?7D>&Hj*gyeR?kW%LM2)di55vjBJ!_*SnxH8&6i}6l*FPI zBsDXy>G52@NtzZQIDR(kSOGMUXu7R^j%-!;8HW@R6f`L%%f zmB}0*A7}FLARI_H(uYu)OcpP~h_Eo3%!z&yKZ2T1W2=EuW?K| zK8$CiOZOJEwi_+uzEWW`dCc_SB-G$_raFZ40Snra!}YiqlF zX3LUe5?$e)k>Sl^@VAYinX@36!~9+{JG!#C8qXMuuUeXSg+(6EmWHg)v*xP{k!J`; zyB)#_g41re+IDF-JD_i%z1?VSw+5n&t)K$U8*D!?UA7nK`zIVr{;1efUvW8L69rS) zWaiG>y`uSkHp5=xveZ{Dw6^sHB&L&znPrp2gr#2pCuk~Kdm=%A)_675#_=ezfCFeFtiB+j&U@0c0{W5!%OUsEw|@)rKowW+lb=yzf#`RI1E;rw=9Cuqs$pbFVFa-3 z8#Oftst$+=-SQ7mMQE+Ws!9KOOs4+*`Q6Ud^)V*ySQa>Oem{*{P_+r1g;>kw+9F7C z7Qg88*aVJdx6_~`yUIpuFleFc!PPK6qW^v&hZl&lrqX~6&+aY%kgN|D=x^-zKf#}= zzaVB=cZw9b{vssDTCcVUm1_AXYD#KpTTDJo-+)3K1o3d9LTNx|DlbiUL6c$mbW#Ar z(W`|uOz)I_J}pWcxWLtDB1)_Hdh|dNZ)T!%%4Nu#jVuwOzqOyQ!iHb?saFvkr6KQ$ z$y~!fWZUB@62|9@$8+Vt0$&xZjx7t&bsRPO=@Vv_XN9&5X6I*?-Ymd^xwoqckUZ^7 zDItiZ2@o=kjYqW+yzWhe-5&gfWgeZjtiOXayCuF&zt#VzI*4z8?!tc(roI8>RY?(Q z!`;_0`AfeRXXakPT~h@*ZAD*&)K#y)1^mFa#t^Wqt^X+|a|3HzRQV;~=}XXj|5N%u z#7xxy@r&aaKQfJBHCp$9@k^Ie`{B321Ku&ebZ(E+X&&?jUE;7h0<0y6$&zM+-%Eg} zkig1X?Q_OFR>vIv$~K~_klj`p(l*wn6z?G2X^a81LoRUDEDnyF?uRugQb257_0kWu z_8iy2I>tNKb20Rel5t*%ZzUy=tpYL{iCNN(wJp<)Z;)sz#;mZ*$cdgCPeQddd6LsA z*@|ACG|f?nS^i!e+~dL?U%)FdQPK!4HG7o3=d8t;75Kg^NdL%}9uO-K)kT~_wLbX? zgwLCO^hdRt%UHV4nW6@oc@!#!K_^!&HKpHxDi@FdY@Fbm2z^?k^G&>p5I+RlP-5g~ zdT5@gN_!R#RiFTkSY7*v!_&m2O$M)Aj=)>jA~4Rn8S|DM*!Kh1kX_N=H4`$pwA8Aa zWGWcApg1#B{i&=^To`dq!$a^&M#ij{d)xZtfNm`H_N8cVu=Z3is$6>qVn2ZjX?$U^ zKv>===B)oou79WieS&fXSC}p;`GbCiB^=T2lojk(r^tQJV=`BMNWuaA5E!FQ|?ga!tZs2m<1FEFLpVUCBFF_8T0y&cTALRXc^(ZjzR}3-% z^#xa*-Qa<5`#J0Qm(>rUw6YY$F*JKAE}Td_CEtjv5WEO=ZoAT%b9SMyWjl-?yN7ErNSlrLg^e0TIB2}@YCK}Btn0&)Lx)gI)+Dh$D3+pNZVPu z@!|K2XjATVyQE+P+lg-i%wG!N77tDsymcO&9}3n`8aG?%`|XxIDJzI@eE#LbWjEcR zl^+hjAtT!)_c>c+jtGS|N{Zn?@+;?d%e@@+78gT%rB8QBqXraMTl)xTk4|um{fKIC zg|Bkd9k~M}`x;<%gjp`V53S`R9htR*FH@GJ$-<<<=Qo!`(CrYVP>-DD&D1UNa5;B@7L#D!xifRR%zF>n^`+BrI7cYM51yq|Nsb2HbZ{ zi=om2%n!NebPjt-uj<*uL#r+a;KHD|yLAL1hgU;RXOUDb-30EemIT|M`v>rvV~y5F z26^$-U2J><*n7v^(AWUl4{IO>f?eG9Wk7Z(dq%{F^-s88U#PyD8G<7AbE%Q&hKWPT zF0CL$^yKyXVaPi{la+Y?6iI#Qe2f+{V>Q3AVMH@U8xLkEVCodZV9obA#`Uo2=5=97 z5ZQoSn$d|W_IbEgSsI7rj*>#GeKdDphX*GWS|uc9B@-ZFVh*Cmxh4M`hV#VX6UAD_ zE2G_80rlCP&xm%Sl%*_Fy)8(-JJ1d_#(KJHE&#?ckwHNjBPUAyU1QFcu-NTUq3Je} zIQDE)bKVb-lVLbHHhP6ItkuYTRNu>(xbOBHU72@@pjE)zu(TScqo-)+K-i>L~zJ>LJgQC?#1~XgmV3tC9iBLFNtxDt( z={bs+om_v3P_$Ei=Ikawy^WR1UW2+Y6&EeUF`2vt-#8n}k663^K(b+u5ZEt6iUl?u za>rINW8g-{^ia~X*KBj!CRZAr8j4QW!aCG=rOTII(+}7Gq)oHj0n3qMg~GjAVQ2lH zZPnJ^==4gg193Ovob|QurpYd&pR?CEbNmUdn!vf$s~aE3s9LU1UL-?B&E&S|Yw!5f z&TTKHx+7Z9O?>do3|>oZe8{#N1aNWbs9tip+fxN3_En8^sLX%h)%|0r^#-3`A)t=m znbV3A+m?VhCDWwPw7b1a>z7$qu&=oZJf!zWUu(+)G-OqeF5P-@g#hWtx zNHnvQe@I#-7aPTW6@2gpc2sLNOdKpKFJ6?t>a8qL{T=O@yM@A>O`I826=V~ z;;ar@)CZGgIH}^GYJzjbzagdH zNA7w#IzFCP4T^qzOWz8Ud8@Qx`vKPBWtjxcE2p-8ym}BG+w93%fWE@a5MN=$nAp_= zus0Xq)E&s|n`&uhuJ1-p@ij$(I?)ns?<+*ELf3Cx;T#03r~yq(-7%K#%+UjMzCxA9 z7z0}$j?e|}DXNpH_Z3U|0pa|AExS7v+ z(4?wJ!_fvvwiD+pSW=jB>(cLpGf)0TwN2>j>;2rbji}1ufIjPV4X&IObY;Hj0MfSN z^=l=*V}}xClgQ9BnNnfasay}(pZJYZ*Q*!6OLNZd`pM8|A3T4`=8b*d!$1LGI)zmS zX5aJ8zG%}xdz-5A7jI@y)mh!2LTk;ugu)G228HZo{^&&f->LjsBDqKXfW9&bou6EN zX84zq5b>S}^|n{@FF7hTI#WEMaZoK%^7D1?syU3bL~@|>(;wP6C$1VmzAj(NE+Q;D zTvh^o4>(;t3D-?~dvT3Wt>%$xmDU3sH_fo?=rpXZkS%(#gEaVsk8571f-TI6=MJk{ z@zG|>_j2rmfOpBviRD@qWZ2g+AxvChT>kVznG_bYlg$GX`0G~pV&d8GNL#8@o{p2_ zhwD5GuiC&#NFe(x$0b_3$3YiA)>BmsS>tpP2shyQmj{BvEDEpM!FH3frGJbo%zOz8 zJHIx!JUkRn_nEzZw$V-@Jo@{|6BlsjST-*W8V^LBPd9F?Zk=xQK~9RRbYI(x-FU%{ z74a;{_J`ef)-19!wd!l)FB0V@q*giKk5xh1G6zl zbs-%cvyIzoeWmRNQoWl9aI(58cy4+`slzzIY^<~8O*`~h$qsChDu0uVC`@1rIN3R@ z)Dcs3tR(6;xf{>lb-RI z=xyXWW}rIt!)u3b|BqK)l6LL^GsqqvwcPYBDf*uZz+Jnp%g+3Gj4CEB-KJ`|_U*oEdOZkBRZ4%Ou@xv0 zBEX)NJ&$*NsSvo%8DMSqen;2mJPbi z9QTi|0EQWS<^wjoP;WdZ>$vb`#SyvE6KX5%wr>=o&F~@CV&UOJZMSo79p-k?`eMln z^+=Df2<2!!@A9;8TSPitfrCuGsBeL6cHaztMi@__@@rsdyub;CmfjEk&S^=aE~!Uc z;kUiQZPuu^As}b-^`u6}>fplBd&lnq=?^1mH(QE>%R(~kH*S)jvC&)80Ms+GL-h=& zCXRLwx|Ya%{wz0hux3qe=E*H#9aYH&V)k-h0q&rOGfL`&*K$pE;5YRe>Juf9oFPzj z%c^VVA5;p);P%v9Nc7Mp==OIJokks6om0Ndy)SYzhmW_|ChVLd8zAL9-Ch}C`=KmxX zNtnQe(@D?IqxWn@W_=1|@Wf&r_RRL`>V%nlT+VxG=B_RKe0t&@L180SM_gs7na_eu zhAb~$8C)`}I<#4S20}Rn;BIfjgX+rK^<%XQXA|!cJE&|SCsyS{v?SYum)S_a8)h1{%yLMP#=b6uX_wdMRw>u)3ls&(n&Q1OGowDUA|OO>=n6Hf3Qkfc?Tr=#7K*ft3H*Mv z`;vXQ{`bfX^ZcisVxMH;$zr`fXkU0I{5(cnNh`I_w#zgHpc|qNA!b$9%kbM#8eLi1 zn+M{qUyv{tThy4h94VQ9a=CKz!W|$zHf0jssU=W&6ZKUm1+x5>NL#8a_q$v494Sb< zP&?Dnv&g7;95pfEgxW5Lv*RQ|MN-AmBhePkH*)KAg(%0XU#_TCeB4#uNr}hm6d0jG zwDV1cudN5OsM7ar_$^f9g(+0I`!sqpA=l|Vk=0Y%LD~Z1to;&ViuZIISF93Be|9|N z^RPrtr5)YG9M=eq`8>zF!mI%+r)&;)e_kf^q9?K>`6uRz>SYtMYzt%%`QdV}G-GF1 za$Y>Jf(CRiLU}j#z3N@tF0&Dcd{FX;E;sY@xgmLC!9VO=!L)Yo1g7v-5t?Xk)im+xjGPCFo3Y~jIEOk_kW8^DXI!ov4P0_jr;qx=l+3&l&M2pahA zTWYNS#=x)t!kEGKcds-amhVjp=V$yXdyK~8bD;>@#>-mtfQle%qE{>H9!lH&*DV|S z4Z&r3?~Dg>JIHyYuqTNw+csN`zDCM4&c!|P_Yv^)5a-2$_ObZIpQ&V4;>ON~|4+XT zpcOB@xV9Udzo(s6a{Xg(hI5=ZazRJk`>;a0dc%i-J2ItS*1rCM0(tc&J;&Wmn;HM{ zYSmT)^nI93-8*F1SiSUla#STF@52C6z%|r~yPWCM6fmlq8?M|ZYaKOKzOd|8i2VO! zG)1875kx=jWR#{@SG!9oJ4rK(E;c(p(A4Kz@YL*>G;_HtdbedTYs6sJ4wFTCGzhkFUHENGPmo0*3uSOb4cx}C4^bSqwl?^=Nl6<_T99MY?uE& zxE3P;W*zw~eN1*uY1n1dGV-)E*zdU5?pDoz?^u3)SL~j=zJZZg#`yl)cWtMgj)IN# zyYv;_2ltAvm>P<4J~L;|QW94WMAp{ZmYuf9j}=thEqid^zda@j-Mn@w4!X2m0%hdBw%DFiMFgz) z0O{S0Q5T^|UO-wcHVVlUwN2plLkrAB|JkvU#`?w8L2S4>o7C>=yk-vfqd>vbLF1L- zAI8A6A|~pN4RT1uT!c3A0n_S8Q)XuEzcp9_6}Nx-4^v@3aQ+Ju!oo1d99CaD{1x8A QdMQ3)W_hsK)cxXr0Wgk^`~Uy| diff --git a/man/figures/README-unnamed-chunk-14-1.png b/man/figures/README-unnamed-chunk-14-1.png index 1541450a1ab6fb11b664887d84b51e1e7fb8c490..4e63f6c096d9ce339a019afae0d00db99f7554db 100644 GIT binary patch literal 8250 zcmc(Fc|4SD+y8aV3}as+yHSW1TM3i(j!Lq1x7e4l#9hRwFk-l*v?xSbvZYO>ED0lw z%UzZdl4NH@Au?kf`@Gjk&;2~ld;gx#`}_U=c;_?o`JC5v9_MkM=l3|i$8j!m^{|yG zKaUg-008`E2lpKT05l8$AjE}6-W+dU8b=-k4qF`Ck30f^E&zA|y77QliKL{Fq@=r~ zQKY1MmXwr~rluxRrwi!n0$#fD-AJ90q`Q%kyHRANCHv!NX&GzjmuneYZ*On!zSh!R z(%l^qayN3%>hxOE>s3f;tYvI0 zTOv2Mo)xbP#Fxa!cbCL>v$EVJvpPF-{Va2lma&$(LqoYktP(GvgvF|iG^*_EtdF&< zXR$i7+&e2HJ3A3s_kdMKC0LUbTa z*&-|!Yu}6eF$f;q#Dn%{0YKn2`v=zfKlTIw+()y0dykRQfAsQ3^wb7&5bJHyQvBMm zr=_gAw%#6p?5nUX+d=P~@eX0<Y2wR`JyB&TdhnjAHxKaQCIrZ+ z5O$B;>L#Jov5)Evxk)SuX}reBpuAEb6PQng$d%%T5&$(VU=Wnkh#r6<2S{#(nC0l5 z??5bH16aq)12#lM)j}F#s3d+cYBNCvqauysG=c57?6fg_(NLjKfiIT?ALy}}&?%%L zfvPozUu-65V6I<+PME+1E(zXs1xd5OYJS9gPSt7@J#}QCGEL5KK~Xxl-8KB4U;R0h z8t3?|I%_4dYVkQei*8r7t!wNnbKHus3=?A$8QSwNlKvaVXAK6S8e9-N_%9@eDFWuX| z9ZW{gEjDE3N^2l34#ED^-i%A1=Xe~HX=h*;)|ad{{*Lp&*mW!T+Id#`N-KLv~9V!{MHCTw@3+I%FxDPgR+nHfr3^gQm3(HBB-cCU5jG;uzJ&-eI2 zL=qWQZ78k0)g$v#Eibo(uQK6TElpy6k-;d^qzd}luFEv#z1gK;0k47^3XYj#_}ZVS zEB+iCDoORfm2#jQVUdCe2&2*_489KhzOEPXwExRe?PfwE^8)w(%(8+?|6Iaz$m|yeb&!a`2Fef2VWk;)AIw-Z{)uD&H~tzoMRxsx zpJc27E@2PCM`SNyD6=U8QUvw5JvM1kk)Q|X8F=!-moP#7J`pCo%DmfyHV+j2h@-Sg zf`S)+(nkgr>yDaOf9OHueq4iA(29g7v0M^)M_t$LfKhQ=FQK5iMR@`X;I@0JfQlab zB)&nT!wNAwhIX^*Z@4fhKzzy4XD~QgFynHYK=Z(tyI`{1gz`@f_iqDvhKc)NZt(TL zLdFHE8Fdh}{%OugKwY2xU!|LJLj@#a$ECtfXSt_4%0)( zHAAueAM9)9`ip%_%N2`fre{Mwra$|B9l{#_G&sW7%;lY-+~hTK1|?Z0Je=L%F_c}8$NLhzjd!WIW6T-LLKqZr!% zld<_?BO|gJiuZ6T5l{!%q}-{;n=uhiw~8Am`p{csw{D~dqKqG{ubZ2n_hBb^{+h

`xU?N=53OXENJ4?ox-w(d#0X-R~?Y5AY|bOr4<_dksj zrVw^z193XH`dn8xw_U5c!arO~7~2!3@XO?%dGNL0_@4X?2mUhqUnO23FRSM7h@ndR z>=9iDhK5WP#>(*VNodBSW5q=$B2Vy^VL7|Yu(hX>04I@OV|M_mB1rVKt51G~=+S<3 zgxr2P7G7KHS2#CO&+a$6eh>Tls0YNXT#tA+m0*s6I2Qqn8>m>1={+Xje=)OBh*tXh zIV18mgBjyH5VzeUB(&dgf#qEGwcH<13pObypp*!yY=>WLHWs9w%R`~8`CoRj6lZilez!94bEIfcrmj`1DZJCJu^vZKc|Ix2OjGj-RK>1E zCIQ_e{2iMbCuONSc7vsB`JCJU1^dmStb(gLx5(IEkBj6cY5Yu|MW1jK(VuvbO)-n^mmqL;>2-3o`p$NSKWsmVC1IsHbLqAN_`=zRGc%mrt{A)hk9W2U8JBM8! zs3Fsg1yDFXBq+N9v5rto#CFv3>Hr!Q^g9yKLNy?ti^bq7<4bjUXuqfW1FQ**QE5d2 znk9WF!#e8#^Rf}IQ`WO9!n9>uA$D%@?f-ofloN$k!xb8@<)wcaBNl}(KN8UXNPxJw zp=2JNHgQ$A3=a3N^a}s5Q3)nuZBY=zo98$}SF0%qW`q zW#MqnHby_y3R0H404Lur2O!u`Ff4$Spg^6!xBVoT$%Wu{viuYrFzIqgmY7;UZ;-TQ z4u0OV{4H4d?T_i{p>Jc;5BAvK+yzm^n^nw|oN<&?vYWyDbN@6U`6Nu)shj}(y5eZ%Jw zs9e*1FgT~~vv)Q(UIreG?0qwnhDkr0asFbK82q+-e9O7H)YV;ybhMuMv)+dR>uGu> z>&^M(O8(_dk^*e^DtCul=A$!HPWLoa-S{N#14maKEH@uyCm|C~Uqw6-jxHzsP`km4 z#>miFg-aoWKdKP!DQTbE0_Bcv&2C1EFb;D`_}cQ~ujPPg<5=-MS8un-t{l!bpeg0fQb*DI>GA6X{T3cJq;F{D7#w&KN_%b#=Q zM2JUT@o_RdzX4W!hKUc5DZ6Pv(`m)WeYU=(FCAOH1rxd(Fl+qP8xk z@vb`@?Z|xdt&#^$%qe7>j1MFRY$T}D#^n}H$BC!yy$#)1@xw4?o(HbgVu?K%nYxX% zsiYX*Mj8(g$E!rU%VPq{)HuUS`O;Y5iBgPyNk~457#KWmZ!WaYZ4-1q0B&p;Mv^FZ z2h|mXe>d?FKQmougV8SxjH5=oU}+Q*{QemNuLc*~qDZR^se2A)&TnEJ0E7FSFD7rg zVmu4QI93ZW4hq3-#JA1G-f+Nh?|u^Zg{`0XDf2s|W;;)Vd2;Qikx6FYOL#~#9+$eD zrNp;moj%B5oca}YxG8VArew*Ya}GdUgqj`OD^?%{O}IGYT?x>R!}58~J3MRyaFyR3 zk-9$tKswJ56}@ZT=Hoj+aMS0BkCej(v;j12POR)b-*yb&>Avy-2P2y2ki|=Otm^_z zYa!+Z9W>=47gE*Ypi0|blas(DG10C{C_>Z15}iCb+fgLOzLy)|8<+VURSBz-Ig@Q;^#5G^t%X!mxG>D6XFn3z@4|#%d*3Uj+uh*RqDMCTXeG_= z{~@60b+{eZ(U(7GKpb?n;k9!y`3?S~Y61{+G2TNp+n}-A$@r`#93#4&!>Monshx1W zlB9X_CS+Zj2@3Yotk|<}gmTjk%cl*M0L_!-=V5!Uj(e`{P_|_Pls&K^-DeBrt#U%> z(lyZ0H=_lSnKp-&T3je}=c7e5%xi)Crx&ch+@vPybO1lE( zT>G}z%SG#}`NT#}osiaM-;B;KoBh%Y7{p) zoqYgqQGw0&A)>{?Z~)U?fKAqosJyZ+78@mf#(9!>98t5CM+4-3`>)W}5 zqsPUHEjkNkmifF56-uwx0-KEC12Nfan_X2%*0(jJ1Vzz8vmNp&T0gi!W*PJNzT_vc zBh|J1z6r|a=wm)cMfc0EJFY>6zMRCtpM4`sHl_eY{}3^04y|p z$W?lu<1pBf?V8)Fhk5ph9hPRn*vWNu1@~3YhMnb>Eocz`V6Tm$6sEH~&l6>g%o)Rt z^T-x*p27G>22!AH0?WI0wexor_J1C`T?6(T>Y4jdb<6b(%+#ol!?%KzpN#pnDGbYa zLJmB$w5-3r0aYzZ6e=K=7mcB@{`E!?AQyX)%4x;?6G=ML83DRxNpn&0&Ysz(^6GbF z3H_R9ZW}apx82#`WRIoRHI%tlSHId71~Yp{1J00-fqCWXU;fzXqnS;vJLUBFq~Q3u zR+&zBR6p7r5KHkP2!<`q;LYaPJJqkd%DKUY5_4+#s4rN77GsP_G8zXHXR3H*(Cy!h z*+kxm#w}ZH08100oZxhn+VVH``y9>2;+qN8LG6y;Xi6sV4-!XMK+4yB^JJPlMx}jX zT=bmxi9ODME0o*(@!LzOqHn-w{D@JS1rrP9wDfu<53b{EHiST53w zt=$uzVDj*C34lLBN0^Bv)a`iu{)4Pn6o>CDtTR~MB|uAoWJ}6f(D~c9!KP?4QR3?5 zb0W~-TP&>@;mJT3$mCX{MOMEPIH8T|AO&%eo;ttx2PxhS+#*Cnix<#fBcPoni^@Ee zjMjho^a94|dV&m1zQqP`O>Ghct&1BLZ>Zw#mbt;vIl&jOzEnYkuybE1=GwV{Q|tz? z{Mx9DFyoJM#@Mj?@od2l$N20F4-9*^gQV|!8M8HqkxeZPH4mzg_`iPoqtR*Cb7*kt zB{gRMto?{2W0#&71H2;zs_(KQM`;DW6-75$zY^@30UeHIT?z*qMQZ>Cv@gK0Sz*Xld}w+g@`h0w5|41j0d#iAkIiY ziAJ5Hgkewag(t&Rg&9M3WsbWX;YM%ym@Q*V+k3^n;^I}(qx!%^D_xg`tEu}-Ep&WF z38Rw^^L#dYY?{1W7jP6OrzP)}MYWli*jTYeqNRuGr84Xsx}!6J10?O!Ksq+>lMD1z zEwh6jf}xXoTS$SB^aW|tH<+GwwK*vAYVi;O;j-h4-+*@Sc&IL+cdEMefPF?C4|hRw z^7aM0uT|l+nQxsIMDb$SuPQVox3X*4gfo7)R^q3r>O`;p}~qb{H$8e+-11q{T5y!VqmT;YNEEAd&IqWh?P+h z@N@$lzEayJa@AYSo(N5!PMlW34G zB0fAq_MR%k>jI&3j^qz-z6*~Q-`ND(sy|Rb+*T2@w##dK1`wx{u81I&n?Rpwp9t}dA#g)K z*Li(|U>9d2TF|25Fcf$7w`@GEcObvS(==6#cPk&|6~XbiZnEcTLB`C&vyJG_rAv*P z^RMbP3FOypAxRTrpg4WUXf!)P-v_RaJiX-q%hpN{) zUzyx?bYffGsVFE|Qa9v>i$>Icb)Kg=XGI#e#Ntj-W1?#fM51rbk^X|*{!M7k$5dTtm2U9=V z++SdRZgimK#B3Y@Fds;yjF*o_rZ&kV;pY9&YQd%t<1Rv~Ac;Aw_hxf#rL83qkhVIb@$QQA`HL!uJ>b&j=ESu8=Cs#NH*D)iuS{??*X3jzr1=q+PNvd1dhq2oV}+|}cP`V% z8p22K47T%_Vh>|4b*P(&nZYW-suezZEXtQm|!(?{0h&428bh zof^J!Qgw1~GbxpSSt;6|h)%O$UWj}A?6vEx=ut?pV28k+kZ38 z)vFMjXU(~6{P-|?(YStmd4|sdMs8BpPMm2fGm|v?KYni4kuy|$h+GPUY{}P-L9O!M z$kUL!7xr=EU&j#tKcyS0gjx+54GkG}l_01))7!Fg-{rkCBadABTVwdITc$V0yoz;L z5h}QV>zx@3J&$rhSopd9ep)!$+8T~IHzca@TEpdEFOD`7qPg*65VLe-Yn4w=(QEk( z0e_qvF7y3JLAOUh)cItCz6@M=koHdDPyVkyWM#cx33{YnWoI3i_ns=5QaUwWl$IRx zT1ENujwh|c_@}+GzL$d^XZ3C!gqU%s5S51IpmomcaitU;*TxSmqdVW))TBRu+C7+? zxEx!qSwz=y9Eo{p@~ppU=D5!MjJ5R(?v~q~qlCG4Mb>Bg@&~GACb(qUGGpK2RG!b| zeRIDR055ePOX}C6J6_Fc9SGFU&G#aW&7`c>EH=`=0Ta%DU^|&<{Z!_jE`RUEshF=b z%}pPEDfH+ovg`@k5~zzZEzY^$^}V>@~HN0rT< z{$2~x<=PI1w}%CxGp#P01!y`hoeQnNX9n~2j81-=G&6jrc~rXmiw1VQRMUVIUz}wa z!Lck^u}e=@@hqbC`*dbRv9no5tgoQ6UT)pTL%~|~1NOr>gpAKD@A``kgF(~}z2U72 zcaZyv>hSVbaaYD@QL8*_vG=+Gv9o&ows+QI#qecn$&M%8!rD&_WILbtRfo^re>N$B zXEp3m&cMYaUezqgoiJwgans0EU*+rz5bSAL^WR<-;SY|T(+K_AnPaG)|Kj%P@&sMk zPHO#lUSC-q+?GCmpz|Bz?h_Zs0|LYrN#E>8Y%Uk+L|VnY?J2$a%7A)iYQy<{|Jj_4 zW60Df`5o3-a(~cv6tzYgRn3M}EvdD8m9Z60ME+x9or`4JxZs6WW+YSA#-Ca?wliSH zka#xxd7q@rICD|MXaI#^PF;4_s%(7c8>|t546P&~2zy<#kNP|i6MM?qbg9jUWyqjr zt!$k>q=I-d=Dh}Z8In|@{o6e$$$VkdpE_y4d%atTkz@?#@`C=`Fcq(a7knXj{$T)G^`Qx}SyM)owP z#Y99Ls*a3?Y%{KVFS1xmh4WO-dFK+Qr5kIeRGr+csHuRo^UTjKsw7USku$>b1?ErR z8s0la7VTvkq+Ji>BK3waCs~i8`;97?mOg^%MQT&4Zs8j5FG`?Xg8NvDMjzYj=~Ed* z55e@+T$|bWI+)#ECcQ|!KJ{~!VIx=-F~5oIMy96;svr3@+C1&yMy(X}a0W7ftIS}S0w3e-Mr$abyV-fLO%Aj@f!R8Pni9>TV4}l2{tShUDg#< SyUYF$fZ2YleR;-&3;zWbG){;B literal 8114 zcmcI}2{@GR*S{V^!-TOf*~uRDRZNnJL=s9RvQ31P7`trK6NO|6rK|}l!jx>;9#KRH zS;jV#$j(F=`}`lRzyI&|eXsX>UGMwOWv=U)`#$%%&;2>)e9rwG5yy@7*;s{G85kJY z436lUGB7Yg3=CVYA{pV4(`~EM@DI-8$4(uFe=;y=FfceWXv8u&mIw-B1qJN{v7n$` zx{#2Ns;VkHr@^41!QiM7+YirS1?{j{J1kgbwE1y0GKw~G%`=K_w6L(S`@Ygn(9RAP zvcuY?lN`V6bu1}yETO}qu12n|M$xXW@L!&5UZZQ?2#CdkAgqiO{asPCQFQdCL|$|w zJywGuwj?&Tza+Mwo^B_YP9o*G8s%-SjLwTT${QKU8=;ptGL+EiRUo#CL~4vSYNXRi z>2{sY}tK2@qm1^~_`)@U5(&>jT1bG|g$9cIuIjUeV7L?}K_hn*#wnD0M zQ1N+FoklzR&#nv-^#_p08n1;wp-Pzl?;(?C=@#-(nEWBa#IPLIQ?f}=Ksx5vSq)yk{4 z$aZ~UD`9c9yuah!O~gi?ZQCHi@FC9YQpoJBJ1U;Dg@j0Mzfa3Tse?KSkTJU7-1!+>=6#Zo-D^S05a{2jE!@)6iI0`ugsdBz40S0IQ-q${9tOmM}JUmr{+489%OQ6K23y55(vAF26hKw4M21#@Zb zxEa6B5bNRh(g^lu31ZwqoJc5iZXd#o3;8qzjZq4Njsl=J8VDOjn;VD|*>xyHVL*`u z(XEXXzIzC_1B5h?V4Vl@sWjSrYdLGB)Wekab1zfcZF6@qhK=3>?;uqFn?(|Mb9Pt^ zm=bpbWFssYWhiKK9o>#?f4{{a>Hh_L|BsZJ0J1EAsBDHEqbHRO6ipajV z+LZF`&l5KNdj!!cL-{(nu%h&Co6^WnPsbrsX14I}o8Br#jQtTEM2gVDRN%=+z)mhVRVZBb)<|amXtYQY%Bz;Ask+Ua$#=OlOVz%^1<=Box7J0!S_?a zGEc8Gqc{3D60>}my$8ac{~7oVO{^_Qo0r6g{4k`k0(zoN5&D5ztUhRT4chq=fWHZ) zHPnFTfgOMRHVDvc&{2>p>^*a9q{J@3@ZWRfI0IP0riUx_Kq@^|4n_CMP-|@w1@1pG z<8K_wqeE%6Lr$BqgV(#RqX8*m!oR1}??n&LfUoEqo!`3sj^kF<~<0j9m;?1K4%?)-jPKo~pndcZTx%QtOT-+8YT{{e0BLTC^b#BZcT_w?MW$Zr* zLhb?NJk{;<54HX=#gx+hY0;W?rB4iZuD=Yq+ME}jfaxb8$SAVn1A8!5~On)39jb0jE zUS4i6>6{81^R-Q_&d(j#iLd(WSJc1Q4XUU-$GTOQAF&PgSSgXb+YlhN3)s&Imxdcm z)mPS}SKZe~JHJkrnwSikW!X|bEeQS!%8bh_m1Hk1Y|W9vV46bJhNS1RQgcC?DCoPI zKQ(y66YV#9LurT8*F*(1kUUiJt?P^1nRKiHt%#ajaR>6Ai~NdA`JnQ~j@d`CVcTXx zbz>O`PVZK+-Smdh*Y3;(HTg)u_a`NIU)HIqa(b8g$yP;0m;00-fi zC%#0Vt%Bx61wF9=Uo?`hUDn1){3MTcIj_Xh{8HWt8b5dv!N(2bF%8g6YUF{)|6m16 z3jd>~nORc>kmL0)uV#3o0`GyV*C3sr1f)XV%7;t|%HINtz(sg-=VEAeLcJnXDzY$$ z@>Bzek^7#GjZcopA;rgR`Dxp3!KGveoE1O(TF-^hs=jrKVyF6Y?jrjSu9W1%*oGBh zAM!6%^iM?6stF%w_L%xuvLk#{zAL+s0)^f2UzBl}z(jqWZNE(5_Za!c*<1C!!ep=^ z*A|CF5R03oBm4|ya^=wyQ!s73b;{D#VS7jF36{XJVv%2wGIFYsN6{lrfV%>kFcC}Y zf9(raJmkaFhztFK(nJ1yE^qlB%r%+O`@?&;4~J?`~$lBzYwkJa+$p5bnjOg;=2G; z#5ltEsc&H9S&DCr4ziNY`zs}4O>t9;*%9WHAwx`(o1!EaMjO}jgQ>HwHyt`{21ZV; ze&R8w2jOl*5!%0))(hz8pw=r9#gR%Q@clU?Zy2S4RB|0hvV%L|D-;|^=lPee z@21(xDQ71fUXuqQZqB#pT}G;t5k`~()hy@WAVL;Epl#53$a2%$0MurPpt#Ka&3$f! znA3_**4`}Smc;YJu_gHzxL@VO(NZ~TYp?4znM<2{RD_SHBT9y* zmCfVNK6BNnX;jN7n)~X@3|@z+D54B8Tdjv_q3<&VJ9|jZ6j6dTlkxfO?MO|Q6Zzq< zrv^WE)wf(=2EE_!$+^w0jZ6G-{z_yD{T$l7?>(n5(9ePpy?Ud@fj*Ad|SaV`Wne6Cxc;YMmH}{%O`(0+C~Gb()F5p_#0Vt#RIZEqlRtG6CS)0Dt8r3%E1pnHu+VgxR*rer-T zKdHVTdObRL_P^@q>3!Ti@q&JByo*J%tQ~}Ek>$|!TBMRiYUj|DJg1*@ z7&P$&k`dv;tUAMJaAMjdnT+4*F{!=u_74GKZ0(p}fe1+E4oiw!b>CGeb-gQKcf(Pp zS<;+7d>VOthsL}Mf0?_VqbiWNpN<@;s7ER(k6BVadP+BV9xFCEi_oLkqs?tvi=bI;?#4412>?Ed9>Gk(Gc{e{mw?XcyT-p%R7}_H(7ZA3{cZRV+1}*;`f>& z2jWGqe;3l1#QKe6OsDuUCk9xiq>-BMH$z!0V(mqN|LC61yBH8`;GOZtWgJ8wp*5UwZg6p%tD)T`B5&R*?q=hu)_E zF)NLxutzwbM?RR5uRc(+1=L@Tgo4;aU?<*$BpshZi$+Hgv#fw35+UmNX1Tvur|=AA zcfT+zJJ_j%`*t|t^vyY#+WT1#ltO3HQ`2p;CYl!(7Gl!VD)Yk^=6%-ABG?a#Q_-p( zj2ecw09(Qr>8h1A1OF@Mu2^!8jlvKya=!L+E&zhq5{8>K2mvsOXYqDA2OrX zz^4M5<4~H=OFS6aW?*Z?b+8Ud+^8;JRj69-ra4Q)apcISf?qXoQbnbt(p%VVX9&Sw z81#(E_nTrWRT4ip7;|8+monHeeQwVxd}qjqt1*8R6zzL+U-U^n#NDr56E~Wgm=pF8T z9Z2Mm*2OiPsu2bmxzg5*R>V*hIMCX)no_$$iO2eP}5UewZ|eCZWpGS1B#fX6bQO7Zx8puJdmt_{BC zcsOGt?(32b$jZZcBg3+T==fclm8l^5y8reED{|K@5}q)r1*<%CIUfQTrb`c29h~ib ziD->^a7x{trgre#3DzRHcvdSLJf{M3ng~xeAzBM}sc7(m`|sW6_$oRrx6C|x5z$(@ zjP_GHQGgv{ezf;!Eg#d5w!>#0%!lreHxcaRTk`uJSVce<+khC^`RHxlcwX?e$enb? zgSZx;++L9z(>u{pki2EpBEf)Z@rsmEfWF}6_!q|mfk;u^{vSrEAvNm}LyWB= zgP9DGtoYLyyxo!h|InT94XIY1YLiyrE z4{+ewH>Bq4*1=ay#*fA`pQo3|QjL3s5G7pP`N_s6+5IWdt3w)>xc$bpUksu~q9usR zD1JXl!u2pOGZ`wyQ>lW#QN%I?YADkMiBIfs5bB~ThU?REv zGzjYQL>E>)%1KGW(Z(##W*w%q_qOd#l?q@{R?gH7FW8VSq%aPnaSiHloo{z31m-4L zi0>*=Ne%?%15VxSsJqX}kB(l|)=18(k@==R^#~W?!)zQTDwFybV`&6BW_dznN6|Y$ zzq25dPre)X%Q8{yz75uoFq%_nOm8A`X1Aj|;G{a774*Rq75}A~uL$~9%m^N?9|3ZZ zO5<4{t0>;AR{N7fAxTN-^1#+NPoB9d>rx`HBXRQ+J3l22Zh5^EpMn3P-jh9)04?ug z7Kom8zbI^PeG@vlc)ZAlZ6!CNw!N2H|I&RBn?H4pGH}ck94YhyNhV8amLPUNwT}jCrwq_V7KWZx`Of1CMpP5^y}Kwh}kwx zpeNw%6O#4&zYxYz-XA&3jE^ntxM6ndt01BXrAd(8x3_}(Y&i5}b4u@BnWIcsxio^$ zEAC{`{)}v-%leWIMd@`5eJq6#Kj}eUVv7D~-s}f5T@YPi3siq5vR@1NVr0w6UYdjfe zZ(Ly2{C4>LjEx83xyRNQ(qNvOEH5Em`5QAr4HWCxv0bE{?UWM_W?;dLx9Y==J~f1L zOnss?2#wM-=Zi+-Z)wxVw~@hfHsIJ$OXG8Iu!!{8&tmB_lX09}NqQdAb?iVUR3s+b z*x}IUgHYz^ZA|nZ1-1s>daNbtD#Cr7QXY_l1jOb(fY7qTFei)VMk{wMBD$U&J(IwO zw8~m@=s*$n^l~w&6>CLEd|?%is#!GdOM)V(uhLsb({7>^iOydQwh1PMjB#S_({-=x zY~|-I(55K7nqt;S3{X(cf_CR8gu(4iE9>jfXnW-#blEB3E3Z$jyCi7v836?bf9%UCaC^eP4=dws+-$STQJ$TMy_DLBI(4Q>yp825U zf&}fS6q-xz)#dRTUY}@bussVZcz11^J^XGj;ueSA1hRwn3q6Z7y)7&$w7S`)=^NGM ze%?TGxk{5*Q}|q!KR>3!=nnTM*f%92my_D=+dP9t&+L=5-f`&}`QYIfL@lPYCFi$I za$-k2x9qf`1ctpz5aN%ztyH?R8I2j*JAd!RCmF|1l(47Ugy^P-1x4~SBWpwqU7xfM z#E^l8N*y(O)-V$-lN?@5#FDvoN7+OBPC(37WNet1)e7@0CF_lc{)RJC^g6D#ox+B) zGP$f=VCLgIh8Wglb+NM*F)Z!s0EBo!yAOkee8aWOb+8Nxg(z?ldYf0}f(=x@#0ofLEc3&yI>H87N?Y|)_T>470IVL#Bd<=gBKp1&-Z;_*4#`gi$Qlv$x2%p6v^sK}ed*4? z(=K}me%UTwA09wwY69-q!1t^G!@Gkr9kRY>l6cvH(q`jw@e^m z0a8gq(wP}hKIcl+pn(u_YA zno*aKQRC&#UG~WyIytP&Z@-}V0p72DPc8ZRu?}fOTE30+XWfol7m}_*mv+L_=-vQs zwH>=8&UcINl0dNo+|Mb$A27(k%ES80x$bSzUq%bG_G4?PX9~*Ht7`n!$ti@d=P5>U z9;>g{mY$p}H;bBYxW~nAf8o-baM1H@28&8tkrjr+I4`!^&_pq6+#*0CfHujbrIniz z*cIqgH|{~&xvNmi_jJ~2Y@iuadTmE2RId4C-}yEoU=c!#LX`8jWa zI`58t{76WZYsW*fE_Rl{8E2zo8{g&Y$=99`K@O7o?2(#4qxXy)5}+_ zr!VkMj@wj04ciYV7AE;di&`tYRW*31$Iw0NWnMTq)9veJzCN>BE}uMJaQYC|dP(3= z8*GmQ?Y!@mp@YJE@+~w&T?V}yTW+g*SRe`%4#mY|4 zO%_E~dN$luq`3>jR*m}>MemqzA;(y~A)b~J*^5_xQmR3p(7@#Dv<4L9*SZ&H9%x#e zjdEwM4j!)->OHFtUqKBk$f;KrzI3{uw6(}w`glrm|DO4nn1&X};X73Rc?;I$FP+w> zj|^t#ADv}%QTfnb^<{=}Of0k2`HmvHqDGL35wZnp`J$55^d2X#D%lDl(p)_Wi8It6dFsO`?PQu{dbz3k|leCB+QiL%Zk zGyaM8Y8f%3xoxL|%FWo;z1)><&z;9aoy&VOve9iBw16pXasLhyrJGXL*DmT1K+{Oy z+Nvd@nYLByZeDBs6@^6qeH;T;XG>K>Wk-in#u_X)@kzj18DgS8I!E1xvF@V?`|zhS z=bS!5h7Z=TMbuq%$^3Cnrs(T%4L4fBBG8 z5q`zLCN(fmnTm=ZGhRvET{`9}86d zdHOp))27@R^A5gwp(Zn6iaa;i7ULXr;T@U&2osgi^6K@_mw^`PXl