From 1f38aaec00541757d7ba92a77992fe10698b59a2 Mon Sep 17 00:00:00 2001 From: GitHub Actions Date: Tue, 12 Nov 2024 11:04:27 +0000 Subject: [PATCH] differences for PR #176 --- clt.md | 4 +- data.md | 86 +++++++++++++++++- fig/clt-rendered-random-histogram-1.png | Bin 5920 -> 5864 bytes ...lt-rendered-repeated-means-histogram-1.png | Bin 5995 -> 5951 bytes kmeans.md | 30 +++--- md5sum.txt | 2 +- normal-distribution.md | 2 +- 7 files changed, 103 insertions(+), 21 deletions(-) diff --git a/clt.md b/clt.md index 9b4549e1..41994637 100644 --- a/clt.md +++ b/clt.md @@ -45,7 +45,7 @@ mean(random_numbers) ``` ``` output -[1] 0.5192774 +[1] 0.4994345 ``` The important point of the Central Limit Theorem is, that if we take a large number of random samples, and calculate the mean of each of these samples, @@ -59,7 +59,7 @@ mean(runif(100)) ``` ``` output -[1] 0.4992043 +[1] 0.5301663 ``` And we can use the `replicate()` function to repeat that calculation several times, in this case 1000 times: diff --git a/data.md b/data.md index 86703df6..894c3a68 100644 --- a/data.md +++ b/data.md @@ -1021,6 +1021,14 @@ _Dimensions:_ Rows: 631 Columns: 5 ### SEXRAT +Frequencies of different sex orders of the first 5 children born in families. + +Does the probability of a male birth differ from 50%? + +Are the sex distribution of successive offspring independent? Ie, does the +sex of the first born child, affect the probability of the second child? + + _Dimensions:_ Rows: 60 Columns: 8 [source](data.md#rosner_1)^1^ @@ -1043,8 +1051,6 @@ _Dimensions:_ Rows: 60 Columns: 8 | sexchldn* | Sex of all children | | num_fam** | Number of families | -:::: - + For families with 5+ children, the sex of the first 5 children are listed. The number of children is given as 5 for such families. @@ -1054,6 +1060,78 @@ such families. ** Number of families with specific gender contribution of children +Example; there are: + +* 4400 families with 2 children where both children are male, +* 4270 families with 2 children where the first child is male, and the second female and, +* 4633 families with 2 children where the first child is female and the second male. + +:::: + + +:::: spoiler + +## Example + +Compare P(child 2 is male | child 1 is female) with P(child 2 is male | child 1 is male) + +That is, the probability child 2 is male given that child 1 is female. + +```r +sexrat <- read_csv("https://raw.githubusercontent.com/KUBDatalab/R-toolbox/main/episodes/data/SEXRAT.csv") + +# Number of families with female first child: + +sexrat %>% + filter(sx_1 == "F") %>% + summarise(nF1 = sum(num_fam)) + +# A tibble: 1 x 1 +nF1 + + 1 25719 + +# Number of those families with a male second child: + +sexrat %>% + filter(sx_1 == "F", + sx_2 == "M") %>% + summarise(nF1M2 = sum(num_fam)) + +# A tibble: 1 x 1 +nF1M2 + + 1 12882 + +# Point estimate for probability of child 2 being male, given child 1 is female: + +pF1M2 <- 12882/25719 +pF1M2 + +[1] 0.5008748 + +# Standard error of mean for proportions: + +SEM_F1M2 <- sqrt(pF1M2*(1-pF1M2)/25719) +SEM_F1M2 +[1] 0.003117757 + +# That gives us a 95% confidence interval for P(Child 2 is male | Child 1 is female): + +pF1M2 + c(-1,1)*1.96*SEM_F1M2 + +[1] 0.4947640 0.5069856 +# Doing the same calculations for P(Child 2 is male | Child 1 is male) +# gives us an interval of (rounded): + +[1] 0.512 0.524 + +# Which would indicate the having a male child first, increases the probability +# of having a second male child. +``` + +:::: + ### SMOKE @@ -1507,6 +1585,10 @@ Den oprindelige kilde til det datasæt: https://www.who.int/teams/global-tubercu 1: Rosner, Bernard A. Fundamentals of Biostatistics, 7/e, International Edition, 2011 ISBN: 9780538735896. https://www.cengage.com/cgi-wadsworth/course_products_wp.pl?fid=M20b&product_isbn_issn=9780538733496&token +der er også guf her https://www.doc88.com/p-5925003681540.html + +https://statanaly.com/wp-content/uploads/2023/04/Fundamentals-of-Biostatistics-7th-Edition.pdf + 2: Hopper, J.H. & Seeman, E (1994). The bone density of female twins discordant for tobacco use. New England Journal of Medicine, 330, 387-392. diff --git a/fig/clt-rendered-random-histogram-1.png b/fig/clt-rendered-random-histogram-1.png index 9d67642a469a630da746792ce70e5e14b6ea649a..28a5321f82f2ea81c45052ba241cb2af84aba2c3 100644 GIT binary patch literal 5864 zcmeHLc{r49+n+JWl8`OLBvBC)(bzqeU6!obw}fmZ%V3NxTcWbdRwN-LJ7es7h%DJ= zNp@q#KDO~mec$^X-|;-}dpzIqe*b-c-1l)__j&%#^E!Xmd7t-vT~WH)stj~obN~Q= zL0wHr4*&p?eqJ=7qbD!77_|TZiVfY{2Fgc6AP`7FK>-4RC@CqasHmu^sZX3ZK|@1B zOG`^fM@LUjfAZu>1_p*xr%o|4GBPnSF*7rtK7E>ng@u)sm5q&!ot>S7gX7GZGn|~9 zTwGk-+}vl+p5@`;;pOE8gTd#{o#W%KR^G43l{_g1Ox>IFJ8QO>C&akmoEzm z2?+}eU%7HcL_|bXR1^Y%T)lc#OiWB%T>RR#YZ4L?l9H0wuV05kp;A&(($dm0GBUEV zvT|~A^78Tu3JQvfihum^$Bi2|Zr;3k>((tLB_(BLWfc__RaI3rH8pj0bqx&-O-;?) zw{L4{X=!U~>*(m{>gwLPb4O24PhVf(z`(%J(9p=p$k^EU?%lg4CMKq)res@vo12@v zyZgh34<9{x*wd^@9z(X!vg{W5C}wIU|>*C zP)JBfXlUq@Cr`q{!otJDBO)RqBO{}tqMkl|8XX<|?Af!Jn3&ku*yqomzj*QD<;$0G zadGkS@d*hDiHV6xNlD4c$tfu*si~>2UcE|7OG{5r&&bHg%*@Qn%6k3!^_w?uva_>u za&q3jeVdz`o0pfDpPye)P*7M{_~FBcqN1YW;^LB$l8+xhmX?-&`t+%+tgO7eyrQC_ zva+(Os;auWx~8V4wzjset`3Pr*4NjgP^gB6hQ`LmrlzLm=H`}`me$tRwzjtR_V$jB zj?T``uCA``?(UwRp3k2@_xASo_4W1l_ka2FWnf@naB%SJ*RMlEL&L+vBO@cDqod!x zeH$Aa8y_E^n3$NHoSd4PLZi_b34i3m<^5Nm(S<9PM4_1wDHp}w98W10sd@V5+MWYJZnm$c2D*KJMJ|0OK(U0Q>+0ja8iy6x z_|(SNAxu3=AD}_DEZ9_P*gIgRQ;-AK-xe=}AaUkqizEgN=hw*;aj`3(bE>OxYOEOI z@J$2ftanlN2D@qPg@>ByryPcBY`?#nZ3ig4TflBkaq#YYmKKkR1wD6?f%W_DTf@f< zmFwRBL9OB1ifz8#IX$5jT!1gGfvH7kp%Z1=J&YO-yx?c` zM&C8A&?geTSw88HZktfzz#P1IZ!sZ*_t|p_C!N?e*w`#)OMe}}Q?xAt(fOw3-wGJ# z^N2EXq`%(xVsx#4D`vE5DT1~%w69#dmQ;Mk%;G?J^&P>nX^&~?J1=4*Ma)hdPGrbY zUJ8-UddY#P*s?y~&AvApIT1QHwT{Q-hs5UJG`}#gNOzzocY{Nz)~5(dDUX!H2@k}K zM5<>(;%#a<*Gy4I$;rSDt`1ZNzps-b5m1n+Jh&=7I^r(JHgR{3^eRzbTzQU85!rsL zzO9Kvo(P8Z*UilgQH#yD@Z;rU#|av3wzZdP{YL4uzcg_R!bQZp>af}nxOF(v=wYP; zJgjrj^_TD>PclFCy+GeK_TBF*84Go(Qx*>+3- zWMmye2^MbxnPl;k$vy${^wdY*ze?JUuifNC-SB@@<(WVdlRU|OEonnd7Wbdt7iReI zCAHfAH^>@(b$SJvX=B&7Iq@UU>El#US5$U+oDtvnoHw~M2K%`_hOpek-4cZHZL?jl zsxsNz`S>VYa(|@>!&$F?l8(nnxYngpJJmpWV`jsRrvp#BZw; z9<1}fET#Ch*={5CIdH19&)yd-4rzt24e4MQpkM;^h-2e zV<6rJ{hP(?Qw@*a6`>=yP}OW+q_ODHamUqVM?$pOKty2Cq)G`;q4|H=dHy-Qc{%t$F-MV~1*QOCA$j;)r2?-EanoLP6mZRo0zb2FWXZKp+@A z=HChy-2G4``wNOx>&Gj^D=buW9-BVkMeS}hF>ZbB>WvvSms!L=^wK;=pAfLXd@fscH&Q4~ZWj&`e_BcuPTMFD_c3OS4~I+RZEVfhD9>S=h1OWv}S zaHVnRa9{luFmoLbg~E9tnh6VJgf)T`#({rt^rDdHUi5KzHOWQd1Hir0ki?xlzMxuY z`f1=J1(MGg7wwXPrV4^mlOic$w}MW{9fct|p7@)if)*>O2YFgS$@)BGE4SJnl^h@> zPX<|k00AM-V!eQtr-(045t+vnu?#><4(tp7yl|BF0q7{JKUFR_j-C~g)pI~{5oZM` z_VZ=+#rN)&H_Rd7AZ)KDCF#3rVfk~bKB1TH0Hkje0JHOSTi?}qhq^H>nT7~LuOr4y zFJS<%_0vFv{wQr3*w-4~ZuDDiP~qK#3@bs{Rr4M>r$jbfYtk%X^%r&-I6MtS;MR+) z^i`oa*$J5+!+~yHQEr0A)f2T93AS-ACA&3S**hC3p?Oy#!27Lt%9s<$0D?#p7D~Qg z*BWQ173-MzqlZ`zB+gz2>3G^tEW71op{&oYez8OK!&?0{rAP&dfupNidu32gm7&B0 zUqukwETpOA@J3zsTptU*U!2GHXW1RW72yl&m@rkm4i(b3yU`f(;mveI$s-B2V;z3$ zFlXJT_u;Q0=IgrgL$4_wtz%OMUqm>BFyknM?bW9yw{67Lxi3`r4Hz`0<(-rZK-ms- z$M?OXK+Aby&6ti&VwTGM9!S5UgRLh7gFlk_`}xlP@}%R^bc}*8O!Va)A?&dys#c!N z!OWrs&cIU%OtC-R7~=?Fu{Q1JvK};KDW9E$9A8RUp?O$Cj%Qs)yNn01B~ZG?l}o;Q z>gNpRME5DoTArSjCgD0nYErc7!exKOC1@sSCCvSrGvZ4J*M4yDv3Gkc)H^lBG;Xyv zCqD?P!wHiqwH;Wzmc-8h{6_`~7_&fCeV}`$I8n9&y_=W&8yYb6{n#Ttd4l4%KGT7O z$XE3F#T&sl7GsEW6vT??QKJNo9NJ^cJV0Yv$4pqQ{>a5?8(x`%dLFqBjLIk)Ks`je;adhb+JxE;l2S!q=^7aqQCyLyKo$lC4H9}Tck%#?pRvPprb&ax5q}j ztt2q(94hbw^;v^U=zkm&Mv1aNVxy6p0x80sO!4Zz2Isk8?hi$deC3 zmXM|zCWv_}x4V|(2gsT0LEEeidp_JK(tr?tI=M^(K)8(9CR~Tl38Qr?EJ89L#y9kDVCl3*O%Da z`&ZTe`VWKyS0d_N!o6l;?lP~-!O_#u)YkUN!n_c6`vW%PAmF3d-u=4jZv}>;qSyVi zqN8hkyU&VJ!yxG)!nI6m$tl~0*>XTqx&nxa>XN@E?}b&|bRBo3ie zHcimtP2ZB;kKSybzEuV*&7|Urj8qV?i*pz%5>;k{dXF5dM^{*zoXGAQP1&4Vpu(0= zck`wGX|k9<9Hzq^6K!Zuc5j*f+G2XUhmr7OzYd^CV8-9d`t&o@yUqJxhr`4Z(QWPS zxkN1YsOj)ZMtE*mZ%kEPY4W(|3ft^sI1{>S#|MbOfyN`A@$)5{-+kt&g}AG^z~OJh z6-JAx9&sK$ZK19(*Jr!uK~`B>ZwXi-w#Z5!iZ(X*eO!7!<)sjA4s!LH}s18`i+UOW2}9}luc(4pW#My<)U@)Pr^jRMB(vX~|M zZSUC~3)pkmlxe8Zbven&+!)s`X;rh@@Xfn?ipa8c!cmc#6)$+CZDl~#4}PM|#)OoG zk;7a*D!AsQkzE(Z@sv=TqFgf<8%s3N)uzZxR6zOd+w?~;^6H;2|G(VB#PGmP+$7$$ z1_nLvOdvsmaoeX-EpPX4_i*~BcT@MLr5p+l5d0V{Utca3K67;|$a|*6uEcjP_68%h zAW0Sf`C@-bds@#L0xe5uh`P?cmGN+-&IP}LlGcZ`t@VXjP_!iHp>8|qqh9Z= zebZIkE)lEu!6kJP@rzOA5i_r%4A4>YIamC=PYld=^DW!w9+2>nOS@XH*D`=EowGe z>&J_DViO!#db2bW8f33FL4^QdWX%=LIe*Cxc=DB5cJEP7g-62Dc8n^S*CiU6f5`;YbL@J)-y=1P4i*v>+K)4?5?*=fi4GN zK2>e?rGB3cw3gi!xGMLF?&pqnpozlWr`->H=U>*we>1=z0U++0o$L3#e&>5#=f89Qc;D;U?)!7!pZj_5`+c6+OQ!l<9KswR z5QxjrK-U5U0#kM$HV9+n@jBE51Y%mcbn%KFqYDOunV6U$5C}6fGYbpLo;`b5Sy|cG z*x1?GIXE~tIXU<4-OI(r#m&tPg+ll3+sDJhvw#17US3{4K0bbaegOdiK|#R-2M!1c z2?+}ei-?Gbii(PfiHVDg!(gz32M-=PbVx!%;_%_aM~)nkl$4Z`l9HB|K6>=%v17+% zWMpJzW##1Lq>C~Kd-m+PbLVt)bk3hYudA!8r>A$}!UcVOeFFmnLqkI&BO_yD zs^awQJX|U%&3==62)8jhi=by1To3czAeud3k$#`}p|y z`uh6$`T6_%-@bi2ARyq*ojZYnfk8n*!NI{HAt9lmp?B}z4GRl{!{G=7B0M}CiA17M zsL063d-v|$zkfd}Dk?fUIwmIO!Gi~}v9S*yK8%Zti;s^_NJvObOiW5jdi3bgo~5U!KY#u_BO~L*ix-)hnORv`+1c4IU%t%A$;r*l z&CAPs_3BlAettngL1AIx+qZ9vii(Phi%UvMN=r-8XmnXwS$TPRMMXtrWhDlKsj8}~ zuCA`Bsd@MAU2Sdc`}gnb>gww2>l+#xuvjb(hr{FXjg5^TK71e$2p>OwY-(z1ZfgF4}u6t!T zKtinJP#c;_uq3lO_#F4=oAKKsmfAh}ft}-3x~wtzdq+QVyop|O=w5ImYL=gGfx+_Q z5HjT&QiA^?8rh_F*MWKO1*KGtV4$jqXL^%5xHJTzI{O z_3p0_{S8eYXHGkJ`Y9Fiy>A7465}q2gMPV015~SCQQF zg*guMNn2*NNyc|_EpLBS>d#AjO3e;l7w_JLY8~?IeQ!k_aeCeqZt7g-M6diHf=S8? zj2PeKQuh#?)gQJw-pxsj;8LelM%F=Q50ek4ovcpyY%i6$gGwlTyTu=-?B_eyLJrD` zfvI}-y@?0Gw+;GF?L;WFLUBga>k^-DgjrYGf9R@POp3X=Je+?jD9a|!c-45FFZ0@u zT{wxcO)gSdC+4}yGUD0Z*ulnZ{y{MZ7#LO?UsWpo^TYbiP4fATbbLb;?A|cLujMn? zyG@vj{W)@M=%q{0>K%Rs@A2p*~Q`e0F^`pDsW z5N1Sp6krcGwtEk3B1Wz*D;zC0SVrKZi+WuagO_+P zqZ{_Wa39+GAWbX+BbLT$7(vnAQ{mZ4GfeX*urQ}6E3Vg|zvs%^d zIzd}0-YGgyt#A{}Ez=QEp{Dn1nifHk6Eu(5vG@SE_1 z8C9>mh~>X|X=OG;b9shso*@cOpbKh_y?$GpUJ{Vqduv-xJCx$(9rC4ot02v-&MYzAi5byP@){#n2>g^G^f(bTw+vdoe|mxYBRaq>x-iU9kMe|6(-O0wRg@p&vPKI z!5GX;`-I4TO+$PS%gI_;7fntD0S!YV3RC3?{_pD^vPN2d{K0+g2@SXdOeGdtxd%IxviS~j!QG3`f*7LP-o(a4$>#*X4jVom_ss)2?c<4acFs)eC+%56K8q@CBoO6ev;w@?A}q`^hzAK_%wh?W zJ9W=ZXd6^{~Ydwznp7X}I(luudV~lc}-=q+o11}=YMf;N90udCi#PS2ph7jkfvge8P zQag89aze_U@MifP{hI}y`-vPON5Jh{tv8!M-m}1Gd=3G2=Q;AC`Lm9y~7D z=?*nvzMtcFx7S)h&d%S#dsho9RtUUnCkd-ml=1Y7)h;L;7nr@Zpz{aCRk50BcSRhD zm-Gt=HdvwQG%hdp^eP}3QMZrynLs+(VWh3LqWi3I2AN^u^!R)4bXK&)l84D7~W3Ce{~NLdGW|d@rDs zO@%2MNJP|c^J~1e*PNgJVSx3$`7iG{=~Q@5$?nP!k}N#C&zbEvLsX~n;3u)a+se|v zII+5QboMv3#^LLooi$RLmKGMSh+mPwzBCOBSTg#piL!vj6=?U{{-%%K%AMk?V3eVN zsR{+6pxd1~P$q27vv5n?#%5F(LaLfEjF{f()~N)O9t|J|dnD>p%tXl zC=3UyW`^r75Ly;vtV;Ps#l?gbSae*cr4i^hTq(Br*Iz^>nZZHdIoWR_186WlI}_W&3(m zpzp7OV5L+?tDgsulw*1FCQ5b&P=G7O)^2i!)=%~gT+Wx_Ba3iXVmJVlGCN?t3-=f1 z15g)hNG6En_blHM#aeHOO5y-`z*$hni;-3JnM1S~m4gyZkSk#%3=<~!6;u)%0N;K6 zQR5_CJ2vd>l`jj*b_Wu*{42Vk3>nC%+|V?u28;~$eLA-SVdz6C8?SGI8DRJfNh%FB z&B#3+$%t@zP|RxNGvCI}i*c|9A((C({6E9bzoE-%nb|Fb5skrAQ6;RR!tA2bUL zSasB~1K4_Xp6EZ$lBI2%bOk ze^sOdfe>?q37)g7( zuK~~U-730}4Fq6YYp=WDeY9#_8>!k?Sm9eun~{Bxoo!V>Tt0z5(DT-3SWcdb?~z&T zbe=jvQ_qxY4_(DB31vie>WO!KuwMJCl$B?*G@TS7wH!|W`;(+98`3Q&z%gx_>~fD) zJ;Y`-kX!ueJlQrX?^~oGb>r&Fdb>YRmMSZz;oT|D3Jq93#;dUxhy;&LiDj@-Y+!pr z8bQ7*i5(DlF%(jN6_xb%OKslF3ahAw_g5KKXc#%*k5K!$5NV`J5xxFi%-24ZkipHn zjEU9hp^w}I(zj#(^2{quR$gBWAs8+sLt?h`ov9~TU#P{;K&rcPDih;P#4dgByXsZ? zSzg0?v5cs1cLiusN{0ZcKZZTy3Y4%Dco(<%9V=^%&K$ynx9Cm`3(UWNrb*Cxbm7a@ zmL6>iwMF*3uzN-dEjj&$U-KA753BMHC>w$6d*$}^UE zgz8mknGQ^_9XrdIUKWy-EJ&eC%CDTwGv*eZ^6EAddmhidimqTK`Sl|csMB*x>rP6- z+(qrQus8 zdA5Xm-Lr^$#0T3P^9O6RprKWzZ`VGalFrJ&D-$jX9S1W2uA2>^XS(wz6gx$$rKe$- z<6R40`um(FHiJR*gxM;zIk8{qJ%bDk7C)uT#uH;=+STrzifw#cX;g8B(x2M7x7z#t zy427p%-bd8@iYlcS?nR9Gc}E=%<@pj+Y@D#Bej-%UfXPGW$;m4MioX%(gXpw)|c6|NV`}Wp3;dWpVp@82MU76IJpAZQ-0E=Z_@B?Pk&Z zgLLtX^9y1v&aLbj(6W1-76;kwuXJmvS1$?wrC}&NVUZJjs;pmeLFxXXOpwQZ2PQ$y4wQHHrD&+iSWq zK4z;P0-D;fGb0XLz4E%2@?EJx5?fyxe>rvMMVN|ofUyi?i`h4?^_n<&V>@fGGwSPk zw|p&Yw+83ZJ^_xgKaq-a&$W2pU_>^<X3a#np|1HH9P+uBkWsN0&EdrAfhVy5|TwZT;kjFN&EHWpQN zeFY}?#g*K5@QKiS;i<#j>^nvPUPXsmDfI{PfhGxRi{gy&^X1pdn1;+ py-9!lPE+6E@AKKeIyY01+jrfa9+$WAFmC-ohI*#Dr8-w5{{wl*jnDu9 diff --git a/fig/clt-rendered-repeated-means-histogram-1.png b/fig/clt-rendered-repeated-means-histogram-1.png index 5469c9600be0491e2cbf2c7505b674e39e3cda51..8b48835ed61b915a40811f5d3c62aeb7ce6a19c1 100644 GIT binary patch literal 5951 zcmcgwc{r5qzn>w=GO4UtMhkgM4B4W=n?foiTC78o?2#B5GqQ$IBFR=HMIvi5V+J8v zCi^}!B4Zimu@8nh!`plP&hPxrd(L~V>zwn)eLdIbe!k1+`?;Uz-l8s<7(#f(c|jl$ z#OU1F%ODV#vDJ2QvU{E^3mAhy9MhLBTIjPO7!2m%;Naxs+_7T^7Z=ygojbX?xp(c_ z#lyqH%gf8h$H&jl4}m~-@7^sSAh2i8p1phb?%TIdP*6}vNJv;%SVTlbR8&+-{5)u-Ul9EzVQcx&VT3Y(RfddB*9z1mD(BZ>}Wn^SzWo3^XIU*+~_uFs39X)#V z*s){s^76-zA6HONP*hZe!C*>CO3KR0Cr+GDQBhGG&!4|= z;ljm>7mbaLO-xKoO-(Odx@2Z%cKPyUb8~YG3yUjPu3Wu(_1d*-*RNlgwv|=H~A1 z?&0C#>FMd^<#p%I9dB=MA0MB)cklZ8`rf;D&(F`#-`^jJLLIni{ zq0wjz1``|{91;=|8X6iF78V{J9uW}{85#NT;lrq?sOaeEn3$MHj~>Ow#>U0P#mC1# ze*E~!lP3uY2~VFsO-xKoN=kb6>{)Vh^7H4i%PT4>K7IOBSy}n{^XICnDl8U@!{P9FJb^%{si~>0t*xu8`|{;WeSLjHLqlU@ z~SBoe=W|K8r--qF#~+1c6E)z#hI-P6<4+uPgM*Vo_Q zKQJ&bI5;>oG&DRsOd^rUWHN<985tS*@#Dwn=;+wk*w3Fo$H&L1RO-aU#N_1U)YR1U z^fZk|o0*xJot>rA=>P!C&CSiv&o3-2EG{lCEiElCFR!euFc^&0)z!7Nwe|INCX>0b zv9Y6%6C&z$?xK*%@eB=2uzdNvnr)OOgOU23{Y8**_hSnPlr{hl#o}e5%H0!G+GP#*0|9= zk*g7VJ5G!*?>^wqfqetE>jV|L1ox^E+{7tr5a7j5q%r5;>^4?FK1(Xk)mpoTNO{~L z$moN(3mXNbWx~9_byNK zje@-Vm4We!HCuC^!nD%JZ~Gl}bDpiQ`%5z~1y6kMqfL|Tg?%4DEQigXB&C z>!NSF7^HPYiLUdFRZ9TWw2kWlTe9vr9kpuX?<{R6EsEWL>9w?+h85p-4}!>kGJHYc+}J&t5#* zETD7EH}pEvF%%=P(Q*I9bpgkAke|}$1M}N?)Z8Oh-F_3%xfNIRZk4jXhu-nQSd1dW z%?}ztzI?WmIc6Vuh%akvWAZG*K;&#@lccESr`V*@ zPL=c;zElHBc2x$?QEq%_>?*N4yM&ahOj|EgecFn2Rkc1}Q2&0-`Ail-EO~;LR!P%d za)P6LLr*f6&D>|B!IiOJJ`V$5PWl=8o?c98@qQBS5!z2?`rLK~YTxW-Ot5wMR@2#> z9cUYnJtr}LZ9cH)F~-bcrGxs&kStGEo^?`Jr1j0?ta+s{XJS_uv1^|7UZ; zFKg6~p0<*qyPouXi(u2CG-|KaOzNMYV7)-SX(lL%7l``nnZ^Y|_WACA(tbVT<}L9$ zX}m(!opYxPtB6`f8Y#&QKvGOWST=Ur%{Io@Axv8|FkNien45P@0I^sgMKE4-< zxRVDf9rv@N9LJZk;A=0>+1stZR|rhBbm6CaX%*$i#YKs_?!;D%98QlzfRJ8wy{3wv zN5!}>JquvcARq`GsK72NQOZK`#j$V8-}3;UYm#wJL+{_6!3>hXfh{h6iDYn5&|^{F<>Y_f&e$BYAZ6!hE_nZ*NRLxCk^WYIZAj}IQb6V zsr5myCV^QpOgY6QWEC2*wJ5Jz)0`qJb%?`$Cr^9X^K>|oZGdGeTFG2we zAXnM2p6v~WNmIau*Mqs?l%Je`GgzQh_0N~|JHOlTv#b8cFgECJsc=ovQ{B_Fr7#Wy z0!R_`PI2a)gRt!?(O=3bkqW|d0wA_LsS9ruQT8>7Aqs!sF09{Vf>R;-$Z>9L}1364-CY|u?v5*8~ZvmJI90{v1wT}g& z;2@@u*zuW^?zv_#-J8+&zpiRCl;M?DW|KeY3<6xj=$MjjgF@eRevnQq;HXfzIxRP> z!}-^PxZmtbPWLx3T^2blw`oV<)G-Pc<%CO9Y)x*khww;uB05GIL&|O;s(Ra z#n4^Zm`Qht(+6n*WJ)a{fk!C$)i+9gc-68r%6`9S9WlR}wz<#U zX(y9);3L@XjS}kFL;GI86ipl!a0_FF?X*k5!JH#aqd9ChCrgV}m8VQVeee65wee?WM*Hl9|wGPXtH1%NDZ zY7CHqdX>{;3hPE&pQVr8Lt^`jpJI2%Ygqdw^p4EkYBrhIf|K0)jqD;b@9+PE g= zBd_g>4XKObkH^S_y&PLLt?LEKY!xLGra```Nll||DYi$anNe6|G7pTpFISm8BU!e%0Z2_} zWnxWa-O}jOu`5U$Yjq>G^U|mS+s#IobmKpO_^|0lxu$$TFPO@Sy#Tg5LXg`H-1vV1 zZ6b{NI{5;1TW0dP4z`red`ujxw9d72Op(V~n;_Q)MNThxfl(dXLT^j=r|SP5&`Iwu zD3PQ%ee-2#vH-_?_V4~E(tq^}ryW1AyM-HjK@Q3;(%-&1AMmuZ5ChcH`!3?&Dp5SN zKBl0}90-HHH(=DP+7S?)fD#=U#{KIn6FXff+wFQx;i!aZR5KvFZRxzXJ9*`+ILn*K zxfMMI=s)|^m19%3b5t6H@8iY}`c>jsK4#ie$jy(sn1B2-uQQTNDiqvxfWU>jIQolu< zQ04mr4BN?0#YxH$6kvaPr<9}t#dIWliDzdF_MiM0J^t@Q`b$*1>*aK#f4g@BoO!Zh zM0VUXLF;3%sU(|orsaZgFnNG(l^>*^E=XTSaZxmo^Ixj8r2!-M7C~F}zj~sqMi5p` zWthdWDE~n8+KBb;wcXqa?2~M*Yf=V$&XZ`W!JHc`=8F_s{Jidd0$Cwov$8rFQD`&d zK9uns0b#UcWP1ME`(<4yZ{fRYSyygIi%>F1g&hXe@Z<{Y1W~x17iEdg*1g1M?M;QCMcZDsQUH>f>@WDu@otMa9cSu%k z9^I!*ZOhzG5#(l!fpyA2nW}_b=$|MtyznB~T-Ie&B`)?wV9Q6dwshSR)!I`l`=1`zFNYUO_4`|yGZr)~webgV_3^F_?2C^Qnv`&g!=PEcY4O}k z#^o@tHuC>ejbdmEWQ2#I?8RdG%OB_WqK8!-sU_ASnY#0L z0<X<`9 z5V!EPXy5`=$uok{Vw&rI@RwnCYyJxC!@1%s(f@4OEFz}C6&agKCKW-L}1lZDZKl^&GV}40rDH zv#6C>-+km`%iMm`CQB{z`N|_$-Pe#sxkm#wsMkw1-pg^k#G;wEVMfbR7B=sy z<1nKG)fH*#XiV`i))CEuz={q!bY*icypj4C@Q>mdfSmvTX#9F1vf?&@37wo-nIa#!(yLM#1G ztnV6L;M7WzvIHDv=gRxd7#p_!-)ZwAEMa5j+Im&>>FLJKwY=qWeXq&@dzY_+fcD_+ zE5hxov!X$jXCEla)vbJG+_yD%s*Mzd>xND&T>t{k3be_Im&YS8tSe{fBAv?*QkkC6 zt5);sOiNOWn++R!$w-@NZ^dFbl_;WxiJF>aq{0gQ@O+IJ@f0KtYUx0(u%-q z3TFDq^5#0x7XHC^%mc|enKMPnt(G5XBASHLahKaK;`z4HYIHKZNnki0y(42{v8yp% zJF%1+z%Z@8|9E-Gu-bEaqTcQB%2kpcHLtYi#`tJLsf7h(HlxkmH5&M@`k*I zz;Q?bG(-Igf0e_j)d&S~yXC86q$Y2tq*Htty(^VfdVCYF(1$Z5C5IJxiG|V;x|%`X zDVEf2YA%fDi}#;+9^ASpbJV&7rKu0D9O3NXg4T2W{loMO3rwkbJEqxt{yF<<5M-op La`uCseaOE7k!#1& literal 5995 zcmc&&c|4TuyPq+3FtP!KEAtYoc2?-6^27|Gb3VA7ctrL>So_(9KS3=gY z86f&wXDnZtH2XvkI_+ zKp=Lln`(DJATashVTRBsF(jxi2t>DV`<9V9tqTT&>FDSn5C}azJp%&+BO@ad6B9Et zGYbm~D=RA-8{3g1N7&ifIXF0=Q0UR4M>#n;xwyEvxw(0GczAhv`S|!?FxatU$N2gA z1q1{H1qF{EKQ1IBbmGJbI2?ZR2A z@!~~MQPE46E{TbWiHnO%NJw11eEG_iE0U6uQc_aV($X?AGP1I=a&mH4uU?gxmse0w zxOVNDqN1Xbk`e-eP*zsHe*OB58#h!`R8&<})zs9~)zvjLG&D6eZ{EDArKP2gtL_ zBHi5F+}+)uJbB{b;o<4&>E-3+?d|R3e0_cW{QUg={n2Q2KtMoXU|>*C5C($@ z2?+@e4Gjwmd;0Y0vuDr3!^59HfBxddi%VzJnF@7|@QrKP8*XJll& zfB!x+GczkID?2+oCnx9Ar%$=Lxp{ec`T6+;1qFqLg+)b0#l^*+KYuPMDJd;2Eh{T4 zFE6jCsHm*0tg5Q2uCA`Bsj024#o=&xJiflZzM-L^v9YnKsj0cSxuvD0wYBxjmoHzx ze*O0CTU%S(_wV1^+uJ)jIyyT$fBg8-)z#JA-QCmE)7#tI*VotI-%lVA1_lNO2M33S zhK7fSM@B|QM@PrT#>U6TiA3VW#Kh#}gw9s+WPwX#>U3x<|c_m+S=M8lgZoL+dDfuySux4dwUcL zWq*I4N~P|PdO6U-oyF~@nFk2O#&z(3yGH}}Kp-9wEwvj)sN|*50JM>DVCyO*`2Ep( z(AQ_J^vqn+^z=9Q=ET8V1Np4EA$d2u^2IE8bj-xgBD`#Qnt<}(`*cUv;k&XukG{pO`&&uW{CYpA%ghbBS7|dROm{}IiCZwn!9~GB&OE|(xN{UNj$2)92 zE_j|SxdTi@8Z-`u$bdenw5)>_4E(T8q z5J@8g9S35@r;_^H+wE%~ zd7^YNj75zFj6@%%35X-wV(_CL#aEB0*|b|2WDvXcHfXCvi5zj^4f;H zC8$vE<0b6c=9oWYGY&{YcTBs7?=JsEB3+4XH|$IF+h+gM)$wjOa>zr?)c`Kk@cDwP!W@j z*g7R)pPPDH{vAH}$4&GHqET_Pv=QXe#C z(i|R?F*A)H4oBXJlh2h}!L|8koYwSJDTz^&9-;Q!BN*d}gQs**BH{9MFoR4Wvrp+l~{#E7Jd`P2V#(rD;V-pP}aa9Dv-HcHEnn zNY)baGH%!CGOj!cng@pfEBhClP2IC-QH$Og>PaL(KtEqDiwq zU+ebUx$$14HjH~Mz~}ao8SVvmDoOJOQV5Kz&F@~oD)nNGO$%{#3d>% zEW@`(9>~y{q^(IQQYM;>7xc_t;f+6{*~j%NSmUB6g0?6%VDqxZM+@z(pcM$2GfA0c zk?^R=GJQUMqvXX0Y8V6s$HxywfB7@J7KJcdVDxr%nXoiDS!C&q@?nOVu5xuHt)ItR z@a(8Zc9;wYL_^Un@Uset^CS**%Py#{*cJCyxR`B6?0qL=Vm=>~rj^?Y;VPA0vVxyGd1+AVibLi{60>_OPq9_6(8?rHp+N`D);ZK;Cx%U zjPS$f*wGfs-qO4P=_>=(H6|k)Q?4&S#`}<+C8=Y7FcAIs#snM zi@ zMTFoU$Wudl4!F3;(o5J8g)EwtaAGqfiU1zFsZOTdP?8}-f`_)B0WQ5coE$Fl^zAI@?&qX9JYh9 z^16OnN!-Px20Z5VSa0SPlC|uomTmPUJCNTvm&{4DYnZcqz*@hMIjQ~>yrK|DHmV_> z{9Ewh5v+n>+BP)Igesyb$H>tsI%3M9zHz&UJ6TiGbG0hSMF%zMUS8|&fvF98+tb48EP6Za@vfFgpa%nhu^d~Stgc?|8(LYHf}=# zSJ`-vD;*)bA%a)PFb;t4uZ2pYL#tt6I}PBm)3n6+M*tXNqL~2|o2z&W`afd_$qI0S zt!LSoP!WHHoh}%08qXFmLHoo0m0%a}n)IkNiQ}|mAT0*%oDU9fgXm)bX&b!1vBY1g z`lM>y13`LaN*-BR?-c?j$fO`mFhTTMpB=f|o_%OoD#D>TRDTer=Z)&1#G{UkFK7d! znTVE{CD)+#>=H%4dsKNvd)*m+ub#{K&E1t6?EIZFw@4J~4?8mFwV{%TDB+ncreWT7 zLfNOS^^JN&s=+RTXucR!p|Ld7EVCGpZ;)M06^NcF-%KbYr|o{2cJV0@*IvCEn1dXs zxWE04faAo$dV|lporr$G@n2q>r%&ybX5MdW-C2^oU(2nABs{ZY#puZt=)pI`o>2+Mpfc!G1*0FaEj_b_HfTC_M z?&w|;7PiaBY^t(F4B_hP3MoCXvk4g&x3{n-(W<2U={0#*xSSQ%LdXjrdTW7T?_cJ< zv9r_AKq@GZaG`F5W=Ry_l7gr#4ameaC6tCCFx7fOEa#O1@y3TuC=6BV`Td{VLB|HT z0#53&G!&TT)_zsuA_HH(Xt0Y+VqgjV+QLe}L==&b{SV&n|(5zYarpV==tBMVLZ%j z8CzT&84({uK0B1@2`|B}bYG0#H0l(t*nf7HP&pXOMTX2TSF(57jo{S##FORBH zc0m7{{bm0Y*^{}g%G9|lo}m9d6wPiEw=K{|+JAGm^MiI-R2e&TR3XVdWSPR>=NQI; zPW7s?PX*ib`JQVRee|d%N8Ar>Jmd3lU2Vv+f2epP;{oCdbWILj&gWMYIeVYb7H6i` zp3*U2esdQ+dYcbJfYw}oJ-c)n4Nnh;ak5`Ks#G`sC|!STwh@e!5m|Tgjm1sua_4>f zPOIH12$4V8zwuSW_3h|SA!{7WLDtHJYmd`wZBed(oz4r^E%oREf$Q}_vbxLLCgUQj zzP>TI37PjMjtt7VH92ttX0?INI_?mP0&@Wz6<9`E8lgk@p_x->M@q63NZF{=l?mM5 z80V0`va$D|T7f03s+7c6bz3%Apu4UQijMWl*_LG^&#o%$6`XK(t(}=WbWTaJ)>v`z z88rLJK=;H`_|wNJDYs3lM7L&)-j-G`ft0#^nv|ik)BHf!lo$E498^GY>pmAY9=0uj zI9;mOEl)?QuOZSK)$1R_>XLYF@}hXS|9*|FgvxY%(9ieQ$FtnDG9pMzT~94v