From ab41f011ebf2eac4b82e6333b4fb6c58bd9fcca5 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Hugo=20Quen=C3=A9?= Date: Wed, 28 Aug 2024 15:22:19 +0200 Subject: [PATCH] resolve issue #21 --- .../kritieketwaarden-hulpfiguur-1.pdf | Bin 7047 -> 7047 bytes a3kritieketwaarden.Rmd | 9 + .../kritieketwaarden-hulpfiguur-1.png | Bin 27862 -> 30778 bytes docs/app-kritieketwaarden.html | 270 ++++++------------ docs/search_index.json | 2 +- 5 files changed, 101 insertions(+), 180 deletions(-) diff --git a/KMS-NL_files/figure-latex/kritieketwaarden-hulpfiguur-1.pdf b/KMS-NL_files/figure-latex/kritieketwaarden-hulpfiguur-1.pdf index 675a5161e587eed2e670767cb93767f9734428f5..1495e4fc194ae9723e43147f03e8a26dd7c950a0 100644 GIT binary patch delta 121 zcmZoSZ#Um?k(<@j(8SPqvYn(jgthq-w*xO!&|FduDi|esp3&LO)xy%rz}dyg#L>*v g&Dhkd&;S4c delta 121 zcmZoSZ#Um?k(<@r(8AbYvYn(jgthq-w*xO!&|FduDi|esp3&La*x10p)x_1*!pP0Y g)x^=j)X~z>$;i#r($vJj$=TUX!G@5M$vM)}09p?nvH$=8 diff --git a/a3kritieketwaarden.Rmd b/a3kritieketwaarden.Rmd index 0402304..5dad098 100644 --- a/a3kritieketwaarden.Rmd +++ b/a3kritieketwaarden.Rmd @@ -9,8 +9,17 @@ xx <- seq(1,4,by=.01) yy <- dt(xx,df=49) curve( dt(x,df=49), from=-4, to=4, xaxt="n", yaxt="n", xlab="", ylab="", bty="n", lwd=3 ) polygon( x=c(1,xx), y=c(0,yy), col="lavender", border="lavender" ) +# 20240828 added other polygon corresponding to B, three commands +xx <- seq(-1,1,by=.01) +yy <- dt(xx,df=49) +polygon( x=c(-1,xx,1), y=c(0,yy,0), col="wheat1", border="wheat1" ) +# curve( dt(x,df=49), from=-4, to=4, xaxt="n", yaxt="n", xlab="", ylab="", bty="n", lwd=3, add=TRUE ) text( 1.5, dt(1.5,df=49)*.4, "p", adj=.5, cex=1 ) +# 20240828 added other polygon corresponding to B, two more commands +text( -0.5, dt(0.5,df=49)*.4, "B/2", adj=.5, cex=1 ) +text( +0.5, dt(0.5,df=49)*.4, "B/2", adj=.5, cex=1 ) +# axis( side=1, at=c(0,1), labels=c("0","t*"), cex=1, tick=F, line=-1 ) # cex is ignored abline( h=0, lty=1, col="grey70" ) segments( x0=0, x1=0, y0=0, y1=dt(0,df=49), col="grey70" ) diff --git a/docs/KMS-NL_files/figure-html/kritieketwaarden-hulpfiguur-1.png b/docs/KMS-NL_files/figure-html/kritieketwaarden-hulpfiguur-1.png index 8c9e518b3432f0ddcef8dcc126f40228bce0afd0..a281ecf55faef30ef8be52b3e403ac01f2d9b4a7 100644 GIT binary patch literal 30778 zcmeFZRX~+tvo=hsfQpjRNQkt8bS(ktl5S8yx_i;0A_9V>_6QH7Q(acd*YgzYi1sQ1v&9sn8cV!NJzINB}5gGkkBDWNXXtdu7R&ud})e9 zLPEQ3Dk7p_E+Q^sZDDPvWb?wnNZiQE$j;P2QCt`ai90ArMaP6x1y3-xqKJ{E!?!Oj zE%f#&M|iP?V)vVo`BS4J%Q2=fi{hst&CBEQrVrTv9;cYn3*tS!^X=~ zSt|j|NyQ3DXuc2kyAH*_^ble!Rj?OTg_WDvw2mFj#;tkIP1vtaV1<-3X+=)ZqA(HX zGSBq>DQcYk`IFb-n@86zk1l4$j}Z(J50FXzeyd%dtT`V-rGy=P1hGXqb5 zj=sM>N4=|wcHMq`^dmpV*BpAP)@RL}Sut|dy*0i^UZu~*&3}&i$_E!c^1z!!E#Hk| zu#5cdcL%02J?8K}kHxg;*PSz?3R8>?Y9yVi>bG{78-hO#RNSgb=N&RKxw`8a)BB>0 z?@ao-IDX@AX;>9KE_tjMhw1347PbdQ=mU)Gf2eQS`K31{Qx`~C7+ zf>$UWtIWoxyn|glsjogfWm*mRyLneNUGS$`MnaVEi5_zqUIVSbJKeaO4o#07k9Q88 z)Z(LW{0ZeD#3Bx)WmPTN5N`5PE1^rdqgU$S75rzcNp>`>1sPTcBc4$RqX=8k7W|3- zIj9-(?Tza){v&-{(GW^(Q+F->-?%drdNl?`>@V@X!YP~Ux0s<7!tME8>^L+P583vLUuV!oP}T!1-i8ufe(CagKc%>U5vt-$vm-wmbZ=kedY z!Fa(6b>hmx4ImGUIeAJ8jO6(=*l$?yqrY_5bAKX;C-mLMmuNDe_WANgpWQ-X>+L^*rt!rk zF^eC5G*AD1iJ8a7KyY)WKJnpF<$M)sciQceD3ckad^a}k@79v1CxoM}sYdp&Yuv|l ze0!RuJlIp&s9i@Ut(@M*e{0r|vIWyzb`@Rde!iSU;qLeS@l8#hxHF4lX-CV&igp5N z$h}7Q?C1pQv4z?&1Wwc8Cyz2gLVgkg#5pMt4BpIgzCqC#X-42>5cTi zf78El2gykMzy2r^axgw6(qe(E^soQ^Y;UB8NB`xK@Q-oHMnYByw$x;J_usz(4KjTF zzdy9z0R4tH)S_4M@&CRI{DPE7ErS2;BH(M}j1Z(ggV@Lq_y7AhAS>+t7YV^1x~`6l zl1M&UtbFsof5Q#$LhS!4K5y?v%sWUsLusim|GRu4NZ}s;O)ubU{HoqaUmH1@@kRdo zH=ragqyG1YP8Fb_K?F)oGw%Iwg~G26R7&xGD*Zo|p4-*`c`4=pKQFCsH82WVP#-ln)Ae zip}ruiq$T*^h#-g%9+ZNuV24rb00l``~$g8JOH-QD^fWKT=a3cJ5=XY7f2zbtD3ET z|N3px5G+ELrXUioo|&R$#+}X0(E54^y;^SE6Vo11=*}YTOlwq`@fOK{V3?8*6PIOk zF(i(DK zBn6x-LmMmLA1J$g0Kq}tRC3V=5gbw=lq@VAH}1&AJdylJ#`oZ`j4sV$xKPp`hc-n# zm|V8cOCrc8_wOnwYFKvoI&&a=yFNMi(&@+iK;reOK~`!@~+1?;|KXXrB;t%2+$ zUqV+mV_X>N#GkoX5xau}iHyzuKuhjIuPW^dVrsSIJju2ab*-}`+Xj6p+fvK6ujGqG zBc%lrpQtt)NUWQ@kx|)xzQca{PhKtnmx!!8=q?YBgtEYpNL?ayvfW$gd2;B00fo96 zMNILZDptBY&Co~npWg$B$ z)!a#DEhq9%KN{s*hqt0EeR_(qr_CM}>()ToUVKLr`^;rY#C8`E74~p}v3XiH5KtpR zr6lG};wKVf!Q{+QEILwbLaFRUuiB}6D(t12lMYfeQ$u@(2<9z%v)@g3l6>w_p$&KD zM2g9W;NTi3u{3!_Sw<*m={p=HB-?H0NmSTuSZn3)y~0UV5@e*Lgd}+EY8)&Vbt=|A zlKy8)<#9j@QRXuZO z&y14TY}t!E+MAT@sYz*1U|lG{O#JBqa`y*(M&hSNWGV@75N1U=l>h?ec+vQf*`=T5 z=R+JEinmur)jW?@xcvPqcBEX*c9Po*COTXF46HCW z94k&ObZYgia9LQ4AErMXZEO@LcUyjZ@Y|2#=ingMvsT4>h#NpC?LqqTDpv^h5pf7^ zW}E4zyQZV1Im^REjSEJrOrM#xRFmsEc5@b60%}Giexp2fCXgJgV`oE(DW{^N?r3h> z*qZKqOA_Sy?C&4N$zYRIsR+Z>u?odZ<&2*c-=16reeVE}d^mQ$dr^eI4CE^BYy1?h z(|6&}X;FIlV~<0!`5U$kC6y8aJJESGg8@Xi*0TslSC=hZg43zJ*-iT0P1lak}AgzX}b zzK&A~3hD|{oTR&R+18Q@5gEbalO}(Czz@?XdL@~oS)y1>`u55O(b&YikxUx6x9Abz zbS4H$7grLAQ%e+-X!|g=41385Yh)+v9nM`e-g?w=4@Cf;;%Pt3Zj( zWWKZX?rfPvTgj4tQCtj*dVZozJg-U&r{$uuE9O7YFbo996=*xFbR{Z36#Ln>PryGU zm6cmyqS?x&@2iLFlDA!wGJbS6x1r>cBFQeZnKsxRY*%R%&gr}S8J3fsYmZfH2&Q0N z87Z;PuaGAryLMH^G*MF%mHtn6Rm`ex)u+7aBTMo&M6YSaPe#f@mZWYS) zkkU4?Rjx`JsqgkZ;K$WR3wY6oa|{)m&Ve4KN_D8Hs53`vb-)=bxm6N`5h4O`K@CBA z;-Rgso`=362y50JKC4}x#)bUb;+@{@HslRd6ui`2otz=IDtK6HX zW|eVuHnbZSe!T~qGhx312&MZoRmdV!5!OoS4H?0@gCy&wH1Hnrwa=KM;^NZDVnL$U zd88y;2uUd^>V|coXL?)VOtXqc<2*=aex!S?IIc&%%{Z>Rp&pS+b&eq^5sYdd^78VO zX6XMpWC$lp?qgDuk-oPjS3tnRy57d#DbW$fql9v}MgI(^_0CA)hpxR7q;4mH>cT!H+(RninhjNX(Q`Il zU^mr7f5hZhhlTf0e;ZpfRSMtFh#mpAHz-4RK_JMmH=w^V7-dRM&Vl(k`ycOarOn_N zM>0j~J@DEp9Mj7Z{xeL@vC@pJ;3kqYUM*i)w_QRiVA~+r`g^U?af$fE7wYE5$0U7) zLP$R07+u#*TB;ugAr8)sckhwD_x}F=`=W21+Ma@&f9FnS))&|KdG{CYbckS@AZd*% zjw+81M5i2mSY@wV3;7S*R`sTX^bUlTT)GnoxZ93MFs*j}1YtUdXf zXV`8Qw3eE)j%MK{Zi;rIbQB8?4$dEk^&;Dv7CUHbaA~dJk;a4lXgifXJz><*k%p0hk=YKK zTZv1xDMxr1SGS*o!ksg# z)A0*Oa1)(NCl7OuQl-yWr(In7;x!~sSLWvEfaiu^yG2YB-6?}m6B2SGFuBE{_c#cj zgoL6nYn49Hb#;EQzW}m?+OvHPp2}Uf)@wpxqQyfpx7?Qm6OUf6;(iQW)f&lpy=wli zTKDP002UTcBK=Ravu)ATJa*F|K(>6!U7mulg%M$vm4oK_wcg(-c=>OT z26IbQ<}Cte`NZ^<-w;E^a0@gSz20BUT60! zi#pzS^%{QS%>7N?uMkf{4uU#YiM;^KC_FjE$8}uXQIGat=%8T+h-*qXi%2(rp#!IT z&>q7PwltW>7jB|l7Ak5*!dk;CG@s9ee#b@(vPV0f?l5T=?m*G+p|Yoez$QxX`f-$M z6t6N2dqh19i)2vAY#J^yCS%w?m<^HCyA{$n(&wxmiBg8aE_IdV zDKkI;_szsBS3|y)zf_xc6bEkWHty%B_aE8y7|s-i!Y&%ZcIqD})Ic4FYdxw&0BXoTWClXGJ#Twy;%G8C}lPF zmE@k~c%!2p-M?euPvN^{mGbY1ox4r_cN7pj?Iu0;{7;?2Kj7bv>{_6Z^nZCkhEF?* zfvoG}j6#q8>|;DPbr4?-E1euKc`>72t$U&QAWQe;ze7r)%<2(oNxrT`mpQ0^`g!*jwpjAG2G+gp9#McBKqp9SCcFAc%>>km*&@@H~--+BY=bwMRO&!3mq`A zbk(>yyZ-T~qY)IGP^L<^zOBK(gGzKXyp=Z7BwDY6=AhT_p}MD2lyX!{jtXV^_oS>F z>zLIlnlaFe^8n{sxn2MYE4x2GG5j$11~Vz8EqrJW4;>^TGJ#>j7&y}jKluAc%<2e* zm&ESzXs4+cytG@Ndb+Wjrl;knj)$lp$_Ri~%d7-B$uCit6eC-=M$>M#ZQ;>=G_$$t z!RBO*;ZF7!G$FDtp=*UG#N~q?>U&1x>FuGO_Y)i6kyHPZS{QoU+6VrM*}hws598O> z?@D~ylm9}b!uIbwQynbACh8**@s`$_(U?eTTx93FT^<^HrsHMpq2wgc!1)4UVFdW` zs{rg1C^5<)yl@Nx!W$SEpFr5*RW_N-{MdH*Th?#pAxJ~OZ8<;EYe)mTZK1&8;?+)( zieP|)4*2h@AYC2{OuSMnprGc`%k&xYz0N#0c6J=^9k+q77f8H&^P9YqEStbn2c-($MAk$0?0CSD|ek59AC9J)wQ52rqbIJ|Ls(v7_RLd23NQ z{C$p)RI1&x^L#FWCh@3x^$5E!cPd{K%EVzkG4_H7^D5-7Qi4N~Y5H4=UWm@e5ItAB zf(o7$0vdr%ON_m6XyWu_s=uSqZuZW`<|Z4n^aGRxhQy=O%*l7ejjzdZUzslSq&5_! zHr~5Bpp7v&;OZFEf)647q5K!V)6HSxNl8gEw8!rD67zxg#kv}cDWh+wi$dydb{Pl< z-zrLJj3hUJW`DZMIh^6tb1rkE(?x*_WOkQYOri$5GYx2kG@J;_ypx16Zq@C6X8^9p z)QOF6;Fj>2R-j>9z1QI5Rf1Ta42n!=3yxAEgv5%?2j#BaCY@LKU6HQY+(D}+3Gs&g zD8z%*e?Wm+%s<{-4r3-$)eGhk*p;}#)+!9c3RCIT(ryb@L{6auvNctyC>CZw8KfZ% z*AzDACF@4*<}3jdW&eww_NL(;#Yj$pe1jJ^%!FaR@%62qgM+(vCt`7?6;#hYH2wbJ z_iJ=?w3zgYoACzjiuN_D7;O|mjHIFPIj`b*d3pJL>#Tl?TOpA{*8Xx>qs3TBl?u4)f1wx_o139tVOM{(Xz6Z#ceJb&C zXVTTxz1^w^%qWzw%RnMHH3}*a=loEduR^ccr3vKJ$eIscY*uh1}ol62X>quRd5R zm3KB{co9&02ZXACfae~d5{3kqtYFpM5!h`DC9i_s^xEzgfAy>Bsx!c*>@5gk50AyZ z{EE*s8ab%)quGt+?1meTvo(sE_9%oHJ*r9{F^i`yVRl{rx>Ooeq_z8_@7LzoXtJfA z^4Jrd6rCz3Hn&}i|CAuneURKxhacR$4ET@m9|DS!RM35>SVpdnwD9GjxFUN>F}Smi z7NFhC%Yuuv%8VLcnD#7q92_lj^gaVVu#(TJ*NY$Gf0e!aAm!DkBU5Ev#tV?;#tC{{ zh6jy86G4p0N0C#;i=Q}WZGRp#fL-~OVKI`2 z_p4%V*hM+${R^Gk>)7^SrEF3L-yuI4Rtqe+)Y=e$_#SQ|w!)q{yK$MzA_e12MpT_G zH+PRgu)qU$;|>R>lJG0du4;HA%?scbqFl<~Mm5p>VABq6=YKqFc?|olMN$quQe|2{ zEnJdWapQ}YRUajOG){zNb^TRA3Rg9dJj0Zob0j>c%N$Kxo?puTHKM zur=ErD&*x^;eJq7YR`UqZk{LQxr9XfDik7M#J9?q9yX2Q<)C2<39!{5mX0T;Z%G@;p@%f6^|IQNy zKToL9WjqFcRqu`5lcmP0gif|>a{Qify47oslUe)$%8l}LZ}L&NIax(+Yk!%$d6k5R zkfYpu9ZsH^jf+6XCst$Xs=3;xnNv%DnHgW*j2uyE1t{KIgmw?DxoD(!c*Fve@1bnM zew1kA&lr}Tug|G5Zr$Qw76Bh~5l}-EkrR+)Gw#h?=H8SV$R6<}Cf!N5ePQTZ=q;Z( zChs20sD|D|ipP3mWc+zoaT?ScTcZW1_4h2XZM<9K^;_a2x$FhLb-S7=v z5vp77SJRUO{P1}~`X{Um77lEGxcDXrW%{>}1qkPPzA-KEJnLkJe!PE`a_><=khXQV z)d>({hiCBc+(ijJ&q8CS)7H8_W^mbgRc8~q1gs|-~(JS{Q#dvfhl=}WA0zDysnVcTf z^WsoxN)U)WcJ6W8)c3v9YiCGRPbsmC+JVWWQzJ40(%`|o#%PwF3VY*dPRkh6s~$xC z2)v*hKkd7$Ae4MERkB(Dy)4s;HgdQ+%5PmB&Y)Wlhn24&o7wdZzd_1Fe!qm=Q-npk z(yUTlLgM|Kw{N5WsWF?eo^gG^(#`CGdIzaMcXu}qDd{rMVX5KM1M?q1hs{!<;HG<2 ziF(=daKm~fS@*TuJ~E-|O7ivb;9zTbzpW9*4$t_TYlQA zrnqC)%W!N;pEyS&(M-xuB8sp zk_hiPKnZ=FD|^gwQA^Y-CxRLiARkB!3!U3VLCP8mQZ|Sf#&m0xd@cE785C4s86m~G zp>BIpgzsbnlBI8CMSkU6R!mHZsPBMFw0;rslfg#;T=?<9R_Uh+tk`pSxYH!zGlqOx zXti<|G*n+NG%oVUbvk~$=ESapfk>yP27n&ND-VWTE)=Oqmv1%B+YA$`)2^`l<$xZ+ zpu+YVR$({8cA>`*A{0SE!9qjZ&;0vSKhh}j=6D6>7Ak@fN*KVUL^7-4Om8@MCj!EI z9A@ETUS}sbgoMTn8~cBQTUW=cWLwXZpz_!C(z?(gp}=Y-oB zL_+uC#fw7QzYUkZ507IB)cU78U+4P5VlFFH4WFSQ_2e}^@o(}g(rZnhsL_O&tI`#0 z(fCCYr6Pz6U+p}l52*_=eY+j^_2l&;&s7PM0M};GXw=9#`hzfic6!n<-v~ICG#o%k zab6Q;0GajW0iEyMz-(bOA0&???q7mu@`;Q;ATBQ~Qgvp0DD5oS*Mp7HIHNNyzZ1tw~_s~IFBN*wR9W<9=$lDk` zuhZs~s(PtLmzfAPi?f9_CP1S;iEvp+jiPz-r0+bGd3*C=x$iER98CGL(nv=#HLHqr zMnF6I?CB6j%P$X7;<*xfj&uJ46^-w;ThmKuP^v9LhFL4pbeP|=|Lt+iFU)M}hT0J;ql zdw9A5E_jRNNmQz|lLQ2kXj-yOn(0fuz!QNm(*%&u_py(B(mYookUYryLILY;`()_3 zzJ8+UvuBck3uk5#+qny`cz=C}2z4|MmRRVSk@;2*LyV=c~m-cpdS%Ds<#`+ zM;zC74VEq(oTUr2XcqdjqD0V{i#(2;E@>C_XH!RO>CGgl*ZHz@6NGlOr4i7z2WJFo zEgK`xQ(E9Rui9hSoR-!7FollSeFI6jIxhW271!nWTP~yT`$9ikTitDF1O^v*q7dR) zr3Y6M7!Arf=l$~Y>p}KR^YCHQy{4wXr`P#<(glyU7YtMb`GGF#dr)L!1!-C%SHPeJ zrP6gWe!0#X=xqd7>g@)2I>)@-m*$r@%75|WT5}j(aCcNX(K8#HeCK@*bI13t zfD|Or8vZym0G9jS^y!MYM&~!cTC-nh70x)7M}jyvXog*S?L0vOgpAQ?!%K3&;&|KPEOVT_PWsKx3CiWN$?@yd5or;KEYW@IH_dA)K#^I z6wPboGq3dJm1Ji-^y7d%2f}x>)qw9TmEO`i9}xru56DWSArNoSSo4*@J_ZG^|64$4QFGn=o+^_{<__In9imt911For4q%*{DdCkCkkA=+!b64`~I={cg zY6YBDDyXY2Kmy+5#CxcqK9I-9z>zt~u2Zp#l{sf2;2ef<$sKt8t{a>G>=GqRdeee# zl5%%_5gskIq~+%y`^^}0TTK+9sOsu!!+~sth2Eo7!`318?XjFW8wB3}PqhM~84%Wr zUls003whOwpnJ5KfX3(T*5~Ie9Hk}NgqK%&;NdD^b7Xs~!&~U>LHmk?-m6;@!*7{P zHiUV3ON-G`S-Y7QTF^w`i~xyPmBmQ0q^{cv2f{*9Y9O(kWJU|2k1%fDWYeigXx)j+ zopVAQ?tjknfNWaNedBH4*XJU+dUaL|U#M^CRl7Vl>(9KX5g@hljkEol2i= z{Po@4NTD-zW&yWd=l!S9kAwCC2>W0I*lIF6yw!Z60)qE}f#HdXlp3|}6(Votg--f1 zJfD)~zQ4$iO&QjI*VTnvi_19;v3^Zh+L~^rhOQ#$en?}%v=Kg3$&9lYv(Xk7HjA+eOgP#hP8>wJ z(gK)-=?W~Gy)Zu}raNS0%+CS{A~KZH0d3u#^!70^|A_~0N)nu3nPNj3VVz$V8f;OL z7nrk&-PTMTM62OMkp7e>AXiMQdkll)10`^)t}X?P8#hAM*KPc8Cf&D!m^6!<5FTHi zRxj>b7x_5ZSON@y$Z{#)yc$xbsUx#s-8On7S{4uh0?Zv$9R=9 z=$YP~`0@QaX=Jju8zw1sa5zr@SMz2G-*^#4BEu z=2m=Mt1nW>*a=Fi{R619O=4X*Xx?jb4+Lx{&o=N#DvXkQ39xkxb?rA8sOb|U@xf@aGNZf)OAJhC6Yc6zr z_Ved^qY|6pA!e?%z$WXS*YDo=Kg`WpQ;N5)afByVQ0f18Ib#}^E zw0c-|+5hRj>pMDi%d4jS+}U%rmUrbQxQi849*a{{ZR`tg>HHAzn@y+9f|$qFGL( zb`OS2v>1oWyEiJttDbeI3l|x+&i*NRxpe{KAYLUPIi7I7#5xbq+S>YnJ?sHEm|%vM=07R$H%KUO?Qr)F)6sIx|2*(Kwp-x;lA_WJ17 zy88rnpyYn&ka$gqqr^&4;Klgk#|1rK*3t{@Cv6RMYp0MFCJG9aS=5ws#q`43`EHZQ zzis&Ss}(-a^*y!D=ugzmj}b`oNw~w7@@uL&%70380)m6fRzwQy4y20Tz)t5{efQqU zM~IYC?$667;Ip!BO%WV*Tc%i@`3ar~oeZXKwl}@loQd-NwXA`6BHRup(FmB-DX|6j z6GsQ*FPl3d(Qph_f_0v2)yB;yLBOpgMKX+ZKm;}$``F-}!}3t~3k^)e>iRmp6(#}8 zytD@qPap-JOUMhSJ!)bQ_RVG5g`P7kp7Kqhv~O7ltz|@6ReA1@=@lPwzvxHxlaadN zxoEzlRO|8VwME`Aom`&4*85`5O==!o$LF3g6~g4E&qrN?1Z|1Zbfgn)eE17(JJOr2 z`XcI*XMV7zx4wVuMQqeC`dWW)kbxCY))G!c?fOD&Gq$5`QP--eZ;`twB)BiW@c($q zHw>M}O$Ldw4;in)Bo8!Jusl#0NzVUC7thmzIq!bSgb| zm3D;i>+W#5(S(iG8O&3R`J z!s4BQYV$N_LfK}+<{Ww?giZ^83N~AF;WtI}s=ccl6gfT7jCw;Hm{9%j1#bn$Zt}F@ zn3rdc=8n=nSoz&nvNy-6fQPLne0^n1$0a~@RtE0dZ1`vM^DY`*2(j3iR zrfbdl*9BDaGz<(74xD7^d6x31mC(5f^)^MA9zM+XMtSXu|hh?@_-q;cRB za=j&E)3qNYTG8tNgmQp_ySW-?GME%;c_)puVVRi^dosFx(u!eSJr;_gRXbXW9bJRc z?imi5DZmBQZYcT|moz9${XK_YMp98x1)2x5Ta$gx@q=f)APLES2omHEIdbvmKBLH- zm`Al}zBW<@)9881_^LC384%b{7{wA|F~aGeKY9^>HQgTn$YbtjW{_K8Ero|ibt>Sw zFovU8^dbm?GNp`0%{NMfTWH$k&4jCDo6hPwy9QS|?uprX`6&L{q*E3i9b}sS1wY*e zuJ1Junz~L+`qR_0y=U&mzjj8F2Pta^5KX~)S9k+y;zQ_VLck419F_BQbwJzlkHQyG ze;oKoHoP?ph9XvJ^IcTbzh9&CEM$K%d-h=Ygp;%8OtASAwVSXLqx&=aKefkxw>1@Z?GVc?K*|Qy*(xOK;Q~=?F$UHd|R+ zyOYH+Jkt_JjiUNIxyx4hUs}L+G0b)kp^7J*X~i!M(wOAbQ>=FEFnVRm^XmP>o^Lig z7vV4@j#GZ7-KAzH>ulQsDU+V31`WmOk#utmN31FB4i--h=K{wEM2#~GpM6;sVS&Ru z$kHA1{QFE<@)q!sO^f|msj4}eAb8d9pgAoISdmyB4J1C>1z1#PTHWBe*yb9| z;pHI=OZet`Y}BKMy8LVPkj(Qid+6vohnv)PE;?>IH@>azyN}iV9lBBD=bKmJ?CowZjPs*!@;fKfOW@nXxmhOms-tYCU7(@LEuQbH!Wwz zUk-cZS~s8N;WwDaZnR7SzsC3P-{-AIls{ZdTrq;wt;B9g$r)~*pe7Ovym#NadgI*! z67B9=^MP!7FqK+*23!*60chXjr)<=I0*t7}$a+`}yt2CSKFw8jBnQ%*xbU2t)5$f)P zL)$VW<}dVswR%Wrwj=(9rrt>si}ge`oBO_v>&ecbY-{L6-^&{x-35K(%5_5V-vpAy}V2U9#?^0XG6 zj}4qE6yF-oC+7MGJX)DW0nN|N?m>3L@D)>qqMx3IX$CpfXYPC&Vjg<|hMlXHnj1g} z#Dgp~?-oj~NW(MGOr94akU@|p@To}9DDk7NW>f;h#72X;I`B?fK$ff6PL;C-(~kMX zL}wavBqbO!YDOQ$jXpfw{NAmwFy|t`z!qJtAdl72%~Km$AN~sa<|!YY)sgb0!N1cj{r9?0 zh19Y~nUAnx;RasMTCmoFBJyh*r>nknj&b1dA9d>lM#8$&WW_b~*50o~?7yYeu92c4~`f_zf7wCBS9OfyDUY)Rq zNZKtK)*js`Hn~S7VCjQigjYKE)nVc6Etqsu+;*H>W%Xcwfh9h4jp1PKZOo1qcGE9u zwO_4_vHU}nKP?=~X2|%VQX~8v%D2g?b__FW1?vnKuPbJ{_>8t}_$^_g?H{#Z35U}#Jw?XyL_BsW7GvdiMJ`~NiFLh`tZIfHED(a2k)Ko#_DDp<#~o)moXF-) zz(_TkP)>RB(4uyUqNIj^`BxK)k16@Ku?{OSfS<7$>lJNy5V4)ho2>cr&h||HNQU4 zqR|s<)C&P6WpjDMLdoS&g?9ufFf=Q6GTrRe|!<=d>qrR|plya@cmiv?Z;0QjARJ^c%|=fyarCQtT@Z!(ZFq6e;zlESF}m zJ8@i!jaJ@Zt`ws87#LU!T`!=ZU4Qf`5OS~?VTnpZ7S$0ClS}(4RJ22!KFC`$dl?z> zGXZgenB`(($~Wnz0z&rcNOu$Qn6GYNe z&#$ICZw{Je=a)L=`Hw1oYFBZnH`~jXgbFp*(DIzPhibh!Ag5Aj-+fXbu#nY4$f{Rz zM!zhMfhWU#e|Ke;v5Nsfz9GTWSE!2YG2`)f*9#?aWH4QTS>4urdrvheYn&sqtS0bu;O-2})js?& zBVbC|eg_cFlz0ZyuzWARjx|>8s?;E(5w7s#=TS(9p^#MK3o3LJ=s2;!p<48#Q|!?j z>z=ox`rv2U8Iic)`|2@)*mAmB`rsuypy_}zQ$bS!KG$sv7LyfIrTkg8z}UFDdQj#h zz0Jwtt>_w)u;PWsKL?DLZ&7Q(HktG)^B?CtWj?Csz~;!Yo%>UEG7hG)OFbS>eBS$( z6`$IOD{I|4huw+Ky{neO&u$Xl_PL7lOXyjY37X4pG8(}HMf$Gd1(8C>R(jW|)3n=z ziG!jwoROdhSH8-l?R{B9l^1wNky!J^UNoJ&GwQ3>Y8ffOSipDlhF5QQB&)yqUT+24 zQO1B%#Kzy>^y6U@c&^Z*+rGaaYsYaIv(fkJ&#BN#b-F4a&MRL6wPneVx0h?Kz(B!n9_< z91^X30{`ajsC6s6K4ffgzkr}IO<)*WTofe1a3&!nydLyoQoxE5SaXX|-qqRqfJaXS zh>r-6$&N=f-ZZ~qqQBTNrEJHMr=0sWr>?^|#jxNQ5j~COu zu2FQSpPW`%EX6+yJ>lbR;D7&V`L&_pd+*z?III^4(gj;bN2Q^To(h6)a~Y#8MHJ%= z-zfGTg2SnqPE_acpKMbXy0-H6mMk2YfblZ1Jytl8nt;in<9#`$JthN=9?NkrHsG8D zIUk@L4P#Ym!#ZJAJCZp9yGZA(6xiU%rkesI+n8X&@EPVc>4NXydP+P*y4A*Z7TSL{ zpB}m`J?O{jtTl)3`P8k(G~Vna=<;}Wy`i#%0$sA5A$8^;T~pe}-8p!q-E(^QNcdG> z$D|PCRfqQdvy(L@*OeVX9I}(Lg;(>0ej}DLWme;co72r{M$MOvvUs$@q1u(I1v!$o z(u+pwZ`Vn;Ml#1U$lt%vy*8Li`N~*!d-R-~ol<~YXX+Uo2bAnm zel`J|XN2$Az=vPbTa@N002_9u{wM+p{>&&dXx0L*yJjgrRBRf=-@Asv<8_*?#F(jG zkTh?;0V*}8tGyYUXUJkc&c9S&#GaQPaX9niRb^<`uJMu`TTE0}{A_5zId+hW51VySqnS0ReOYI%GJx7eN9*GXyz*sYdgO8F(EZ%1 zF_>a4To_adh9`$Rg)SQ}!I>{t&g2R2u^FZ?=@g29<*@~2m7SQZLEV~bmo%^q7r4;s zz{~TQKm8XzfI_3!C`@r!9n%06*=)FIo^H)Fpn?UqtVT+qk zsl5$)x7?{yGc&Ys{trHg^?{IewT}(&@`n~to<4Oo+J%`LwapXZ+_^K4=cgO%X!Ep2 z6(I;CAc)*CUJDX<$baDj#`q;pj}KA-No4o8;Z{1$NU0?wXbMFkqh9-uWhI*`NS11_ z#L(^yK#bFZLU7osQ|1qAjo4C&X5KJ#<+717^#Rz6PJznM17LfAA0J6q2Q*+RSj@s9 zB3cL+EqM7WaBzICE&A`36uiAF>W|j2%W4R`v>p0fj6-#cc1CgrZ4R$HiUnTGI}m|J z&lBpB1E%w^Yk1zP0F<2=>?K&Fl_=Hz12#VGnn1^5W`E;v^=iVsccvGn-vfc@1;8GA zL?5#`ewl99s?Xq9pfB0cXfpw5{8ZfO*zctqx<#Dttsh|*| z)T)5k#OzU}RLj*aWxa3Q8aKRTkyjbzU^WBI88oleJ=!6^T)%>d!C7?^TS55j5;!l+ z>y+Nd$7gfQz7w9=#D&6f;wA~F41BIQN2|=qj0mxh2RX!4&QKIl%^9Q9^Ac!wvUFYz zkgVBV67ZhDnioL6AbEWPL1G{CyvP0ah1G#!N2vsU*A#GEyNy53(rF3}o2kG*B0!JY z^=PFmUh~}rQh;DpfU-id^3sj@G`xCeR$&x3p;peL9R6&d36gmMl*R#!{9XPHN5PW|0;x!n_kN^-m(~~MKC#{^L*?(`!Bla)hQFQgbd7AA{Bjx&q6x$= zYyNxyKJ=%q4!O?R_3bWv3p|Qd&q*^05VwVHdsbj^gAbLICr9k*C9o?&*3fIuJY%vo zJS5BpJ?j^O`rNM9_;t+?F*TcVUL1lO-b2QqxyWEYW4Z72tNodqGrnoDzBlO=ODGqf zvf?2jFw3^}Y1#(bv*w4vjXjc@%2Ub8&i59GBmwnHERawqOrt@JNU8l(@3Jh37+XBhc(p4RMjx!tKD>~_Pc;kuldn*s}~ z3dZ3NeJ`AMq}Vi+RZkEDG)h3TEJcC}P!J7X${h&;jKPA3CLbbGA6@Fer4$lGykq$L zbq<^HF0Yo}l!}tFMh#uJ8*Ytww;TRttb7H`v)$o*To9 zPItoeBAGO(z{Cb%%$M+(rC(ea&>b5&{oXW=&5kP8dqQ{b(gl^Vr=9RV)pe^Wdt$1o zu{gj;AI~fF(4fKng`4GZp*>^!8%34OzL2wy(T7_WLcH;zDSHp0c}aUr(6!684N7|W zTMFniZq8nC74lQt3%xxh=4xQbQns6M)hXjOI#vs)L{b2N?`j8cR) zRqE9Z8L$v9SgK_JYuk4b>_tUOV)@nj-Le2@8m7abwlamxA?1OeCQHIj&y=gTd!*?-+jU_1 zF8jazW&qa*-%-+<)QMnXkSw`yqb7GB)a$zv-AS}FT%+~E)^!{*nx*cdQ2w1bndm}D zg6ph`0cxJ^i53_$OKqD}rz=ZT_(HY$V`vhPO;|lugGd`&tqaqkWgixnHlV@D$&E__ z!4pRp9ms6>fJMxTTn-p7`Hn&@7I=p&?knhFWr&l7?H8nGOXNTHOYrJ$n z-kmh=nz`LICXjxZEdSPQPY}9>OG?TFWxW_$%!iLtn3r&s03Qv&1i@aEnA5_9ANdO^ z+WgmLdRgB5?!3BF1~4PIxsVP^bFgi>oRx;~Vvc!v!3Qs%hr&W=r3qX&-j?BR5-=x( zJzBBCg?HE-&9O&UB_o?W_Eh9IW?^73@E(&S471;jx<>qB{?Y}!C_|ys+Q1BKQ4`>I zL5yI69yH)p2=Ijw(CC!2olKB&2J2=-CvIY*q7op1V9G@3$Ho$nDkITW`tiKjG}C!{ z6fuuha&aA(o`hGF`oM}|U6IVAsb66VIYol_nR{Z;;FLWnXlKfv+yyxO$BRUO2H3% zF#Zs1$eZ>Su7+ovy!lAR-aAjXor=6RRC&-0`_n$?Mq6Xt<3>luExDF@;Ts9GK>0-@ zVIzQeeunOOyde%xnCBU-e6M56ZFq}0m4d>Z=me@;nUAZhZSz_S$!@GrgLTlrXP$KB^3!FIRuKZP$s7iEYk7#LL}^V~hKe=Zj7?G5oW8n1x#XOhk(vezLYEj*sOrZ z&6h7*%IO;+M((f<37k}Y#&cqGgwNWQHa-)`V#u`t9&!@nHA_vQp|7c_Hv-7KqdG(2GDT-WT8B zC<~l~WMUiNzl=0ZcYNr>8b>A;B|P!A9w*mH41FBP8KINWvK01>1j9>_*;;Pus(HL7 zScbd#K`5|F7Ub;-nYHe1UD>0X{cc<}eg;omvsP64z)sKDQ&45`72}Koe0J3YC1Ng#e=bU}dpJXxXKg1;`rlO!J<6jo2+TiQfyhVR2!^(N+}GaVLs*bEwO zR)9YROz{|1Wf}sBs4Gvw5|8B4!149D1c=7vi$}O{_IHI~O<8$6p#pAB_x4c0OLM|m z>kcInJaX-Yncul@2hAhBPZ>oAZxFdJ6nG-LyIMhs+i-T@aq-C=kAMdyl4}dw6LM=C z#ghs(bG>dtHR2cMWLOO$7&SP`a7$?459W^ZHaJXg{`|&yYxWk&ZoI>4UvKZd zo)SSbe$Ve-hLg2ArK|51xz?xSNxsUzAIR0wns8eaI5eok^{WH1&c^}lGIu$=HYx!w zUJqU=lIpxRQ7X5KMI1jZwB03wLrfg=;G4Y89Eo-HizE3C4*w4ysP>yzFJ4*)=BzcJ zSTnMxv4(xXw{~<548GZlxywG(^b2L#r!>dSKuPN~T`w%ejOC)qSoqJ6* zmc#tPTJ_dP@HQ%X1ICu~StfEgfg(Z@X?EWma)`yGSNs3kyYg_T_qUClqK++l5 zTgsMoNS2VbWDteyWXUoT)5^(GhY%XFW>1!}(}EWJZYESBWEYYeu;Z(6!i7r;coN+~T_mbFZ?)nwkzDUq9*w5aiuoGAp`n zR_;negb~)x#WPqqF76_WBueEBzw8vx*qREjQD@R+s=j48I@h{3kcL$=Z$NT!aV0>) zIP|5fYrMQB1BwSPfj280KpxE$P z<)n#^P*I5HFc7dmHSyUFo+`;pvX+@xapB<3D#MZGuIJyW)9;*dZ;NP%hV}{S2|y2d z`P@u&W@&>@h^Bdsa?(Mk=(!9bz&sFr`4$O+s_17uM(UE3-tV%-Mkw3gui^srnVktn zo&L2aki%uHISTUTZ6Hk6uGXF8aLYB#isU~T(iXu#=y&>T)`~9)LtbOsB#wG4k6+du zD}}7v5h|4jNT0o!12BE^77;2dXF|zwwvA9A{mLf@ghzv_R?;# zzg8SMTrI_dzrG9!h|AYW`P2v%kKknII0=-j4yOR#Y55Wron(Hb;l6s%r_8Uk<=N(# zqFw^QiUH|HZI}BSAf`ub(b>e*3GK=ybGVpuu*h*aZBX-WNbRV&2)KjI?Jg8tDPo0=VQQFTI{}? zP5n8~42=qRkOAW$t)Dczi+6etQ7!|;0p8N)u7R5egvX)x7vzcVwf(fO^Bwc3VC}T| zZQF#;9XxYK_e0g6OA=~NeZL8GZ<923-spR zHZiDS+NtSEufGYmKgz;B#`dT8u4CT7K52MwvF{K_FhTKpBSq2s@z#{Rbrh08wEX`4 z#;&kO7@S<-=g3r@Y1L_*X3TpCZ?zWVPpyNMWNjcElP2uP%n%QP=9NHrnCdu8F(&X) zbly?O3)9(RbJ;;bLEA1%ZWce$TG;_gjY7Z4PK@*-uZRpBBTq(nR+w}n>PTY}3zSn6 zTzayf1M;brbADHAVq&6fh~YC!OMy(?<&PJjS;)IKh0#s*I|EKc`(PPrabA?nJPGpN zc|zkbjsaTY5QJ0|hLJKYY>RdGQj{#O!(TLp{%bz27PpDx415M(XEstBl%;ufJ)>=Jqb8q?xCJbY?g-^%1#gu`KwG z_vuZ(4pZG#&3orVspnh$={u(1b-UOa*yaqpb_Ogd(t?<8o+}K|2ww==+_v-^Ve}3I zIXp;Sa%YtgP7`<_ZBuOy&PWBuZ45)SD|Ji7LAZLGv()Xvt7Y?FSaT>kEePFp=F^ zCmX0Ush1m7HQssF^Kfu6bCPmi+_7V1|Dm8y(gdb^38wYw zey21?>t6t%$Fr+1NB0KU3D{zBzzo&{jy8i%vY?UF?$Jo(`RKsd?ystQIWXpI&c?n# zGiXlJ1g!~3S}|8d|z=+ ztcE}!-b8eREzpC`R*=^dAqUWWAY0xjaXsw%`oJLPy8ByL*alh4ybPz!%u-ub7)=~i z(Ywq;2bW-4o|~8Q*h~g8_e^QnuWf74bk5j)hRfI9nY0FlSLO(H0 z*JAIYU0QHyXCP{*$dwx3x0sguoMt0bt?+7J6X$%L=AZ`ik|_7mp(R?<=b5}PRSdct z!fkI>Dh^Eb(Az!w8^|0ARL~L02(olZ)$G_EA}~tJ%*EjXG(VcR{Mm`ftxqrJvI6oh zay_4v`&W7<1hCpT6pA6}DH4U_O2%K2^buhWdHm^V1n)*3p3`3kFUNwRm#dHJ_1Yz8 zoZA4usb`2+2WU}QYbB;N^aC!k94t&S@H*i0?%7Ta4Gq(RlQLsgCqt6vyL5y=b5N#v zR&~|&!eV4(}Z|YaE|0(oa3wXR##rAg_Tes;Ky-0+GKK zMW*X)4+D7%k_a$+c6(IT0bQiD(sqYlAr4{U)fO6#)}NM6oAd)nb1s9=jzg9>y695p zJQcH6kc_MUREW`8_6kOH3r6;Kl&{-Ge(OYkp5_|-p=EGQ3XMi*Sc^*}QrcX~sO^D-S( z#^Y?@ddMg-LGflD{j)(#rbho$bxNpI$h5k5=kYbglnk%lp-|Fg(e2@Vqm+zCCq4A1 zQ4(`yYi-VB*SeB*=oR!yI!V^fjc%o$a3w6C0j{>hcX@?h;}mecV3sMW7L0785iTAu zNO4N-KcIKJ^gYRoTC$G#!I4)QKG}!Q;gwG9Rar*MN#{-}zn`j}p5I8X=rAG%ybV$^ z{PZVt-Qd(8=F(w!J|{*nC924;;OQP&Rur^K#{N!$iAT~d{E=mt9NlbWDmQT5xU8XT zxW;461Ho64EinSOVsmJJ-@ngBShF8@DHC-qt}Kj7f1~%ldcvAx5FZ<&lsg)+rY2%u zt}@u^)Z3tlobnW^7xiu+9-A{yYp4s}P@adUGpTh;Os9u>1&XnjG$iglu}kBmrU|{s zFZ|MlDd9`ErFqOVxv1=a#id=~CTU3SoZOk7`M1Qrj9{)KL!$@UJX?B?iqE#R(kX3; zB#8jedH0s05}Q=^(nAxDQHT1S3V3s71%0DWDq+3q*MD8i6CUc_V#_`zeQ|D}?n{D? z3R(uzXIK9wkU?geWN1_6qP+Q7SlPG{$S@{m7S?Db%+_0qryt+xL*JMF{Ddo+-*?xA zzrmdm5>s4kDDTx9^pX`qEeyLiW#|Qb7#s-E`0HGRQYO?ArpN(H$aRTNWB@`qxd z3vkA!fJ5$XamkZ9(jt4&olYA+Pb!Qq_!Jm-NPG8U-0qRlH^t_vS-9vK|7P3U#0L2$ zdR6$X#Pey}c-J>Mwu90DJ)w64-6^v@^!MtNf*b7I&Rm#j0xqkV|KMQg z{Io@s&J4e=y!)%cttQnn)XMSe&04xv`M0uo5@&|ym%sigbIZ)=1wTH6k1%i599r6< zBCP5dp+*psNA8DQV5S3IJt$_6u=YEs;#qVrl1W$f-W^@5hNq01aZ&r<3GN(v4a)Dg6Po#$c7E=Y^P>2$_5Y?5> zQ`GOqdcUMOpK>OO{kA>-J=d1!x5u~)3*BC+EJ_Yb??2m`Smr;CpJ_1RAJJgKN8P;$ zGSF;L2GZW_Zar^jF0li>PVtclCKm?JocReHftMwt-W+r7(dF3umxY9WpWd!&zgIEtj>O(v7l>2}6U3Z0afh9KYQFi` zUAww{9oZN5&KU{e~%mo&f&I|Fb!r$%5KYp6FpPAYDv_Sfy{}{`+qw$WH zJ8YuVocDJRyvojo+(6=P+xeTDB0#t#o0^sMk0-6X*o=o7o`yj@R?_nQ-A$i~BAD8oa`*gmbA5la^)QGn{X3I>H#JB#jQ#xatbbhhwVcIm zpvW-wBE@C9th2Gw=U2XTd~yx{_xJ4-VTB{XGn{#Ab#k%dph>TmwQr8}o4DhYO`|U7 z%#F2JECYp3Puc&&WRQ~cvnM*D47_^lJ05W|6S`}wQn3)AAhs83e-;DTd{7mp5g~8#TaG)o-rxcZuXSRVu&R)_8!!?yl>tS| z!%srtEBh2vp`Q$8u0<&MkrQhgq@_mE4PSk#Cz9MuSs_&a{y-#K;JW*x%&QBJ^vBEs zlMQ+>=H~@8hA~eKIi(g}YNu$!j2l(cD5M@h`R}L!NZ~hI_cS?JQr56PG#+Qx4Yl9AhoLKebXwAm)j9Zc zgO}0D?wZ_!xh`~35xM~9=LLgjTz4ha`vC~~ry@kLy0M~0MMaj?elFy|3*}=%teiXd zJ3PDT_jx!VXTC$dVf);+H)si{X`m=?Z%B7T2;2kJf}TJ4tj~p`QdSNcUJuEHKUOzg zvvI|f(}Xd=fwTxvT6!1TW2i9f*IBbj>=Pe`KGuCz;0bfsG$dJR28a={U=RT3fbq&U zQd`;uh|GhfTWoz}B};*QkB?i}^r;lB?2$6a;mhilM{%p0#)0n(NECXY;F*seqhebG z0<3dBSAKlF!HgR2d&MGtaSKB%-wmih3Q)i1_Gs1a*%{qb4LaM#CyOT8=8Jzv_P<`(@fG9q%1L85WUJOx(*?J{~BLg|}|O-@#x zJ6#?XQI5SxncDQm`D^^DO~;_XOz97hITwSQ4y)$YZ_|Eu zQyre=mFt*Ubar)8zIcO`p_df_Vpty(#Z;zu?}u@|341Y#-Bj+!_Zv1mb$>Z^fc$Xu`R)bD9ZpH90rUk@aPM(#`%jReEIy3f|>%~CezqI){UF~;U`{&lq zA;)Om)o*Z5{teqv-(Bs80q?wnY;CuZF60&NYf2~Is?^YicmXFFc2(GH)!ti4!=+FJw_69 zTwNe$7EAC^mRpZ&t{wv84NrIJCRqFEmDd22vhZk~P;UV4+`Ef4Q}MFuxxPx_{o7I zt0$}R^oFg>=s|O-U8^&uxDi7p4Oqm;=u~cwc>;vLYV=pCE^%Qf!`7|?tLg%NsV`#| zh<2Zo4^e7fDSW?h=nvMkA-&_etE<0yH#h=FF&06Lbe3q@YlFp^k)*@}1uQH)I2IN*k^mok z#r$J!G!_=_Ig|VMADP`3zi(x3W%I=Psh)wjfu(_siJpSEC>GWoUtc9HV``;K!cnCK zH?IEj>PSfmI(NVmQfRKw>fJYWU{GK&z!q#?crT!CZZO8=8u#d4k_n6OrF$e_(nf+u z4mG|oUop*?f6ty!_$UF_>-tv9j`+tmveUDr+y&*qpUqGW1KShPi|&&{whKdq0iW44 z!iJb&Y?RsTw}!nYzU74lf(rTpf|K(04B zVrxi>^`OQ zwN~Fq?T9oY8TN9r{~}obbK${W(Jh_mGxqKZ(Pt-!X|YEBu(IbOHW_JiWldP$H2A6* z^VY1q+4?*|Y6x>%U^g-m=x*UpHhg!FZQ=ds%6a8f;qNM+V#7uEb=XTT{bUk)r5%0N zuJ)$=-ujM%N=zidpCEoRLP{SdPUWIy(OM6cBIYC#onku=|33q@vi-sJ*l|^Hl;e+r z>4N7qg}+~p@U4S?@phiOd{g)0g8+JB6ITu0#MCpCQ^p^sl8rxZb&kvpe`VP|z$#HR zJ*b_XYP~6%a>%&UlKyqkSc$RXth=*uhJDLvLu+To1Ns1V2ZzT*9qO;FN zwFTTPRo<6BRE{g}vYocOYD8(uZdcpa5o6PQ_*G}JJ6C6tPyOt{DR&%?*?#2`p4r=t z%4giltTeuhuL{ZO3n#NHnLjN^Oi<%d-ht6IV7Zk#l@l1@YwV3})`%SEG z`t7T}wEWcNSKgAKH1J$W)^h)7Eyc{hR^JpGP=l( zBWF44%_u2C!SK?e9N|5`p5ovxE!3Ce`O|jUToCV>tB&h!;Y%W4t-UCQ-&a1KTkf!# z&TlyP$Hye5u;|0gyKi-4qtDLdaIun}9j}hRK3g_bPTiVvt|;7i94ps_>&{mz$%B2e ze&=KZ+c^9?d(3w?)r26*c8wc&kw&AaMXXHW6sHXZC+6?ZhGBO)y}o@>k*93V zpqpLSaI$0?!s>BFX&>(HBivSTI)<@6#~oo;jC z_SA)aJO%Tvc8OaS=|$Exj-?E`Ogak3HL%;JL@!4E4CbE1)mpxKNQDoKC&uN{%VkyzJv|NEo_xS3ekEyyhOGq3*p z2XNuNr~Z3UdR;sMB))m4!p;A*0|`A=z3~4w?a|S(Z@{r`>t=+#`2V;7>_;2_%dW?p zAW#Jc#MUUfboRf0fB>0K{2vzzek;8C3<(xdPkOTMe_J11soQ_w81DaFDf)*0?@9^i z|1U~0@%8@~rJeyuc8Y%>u-4vlw8p=|kLIraaHZRYix(MR6VrulZ8=FwNCZ4dm8*;7 zHuAAxRWD4tbVFH2DwLJUpdl>u%w<-IP*$a+hrW~wIod*R-@K_^w=jD7H-J4^0Fd0A z4OKV+5`KnUmrY9wKFeq}tun=pxn4$Q`PdNYaE|fXvc{LMyuBrs|NKZ+&Qfnnk)@jM zFUjf1*6NRLy!#L2hQaQs2-b0(x_JHsZ8e!m$jF#4UXuyfUY(?37x1UKdyR}!dyVCF zOUok`wfv;kXNC7)pKD8%r(;pcj+VDzr1Cn!(|}EY>#LRLdSDo!q;JOKwH&yfgWCAO zrcuncersKD_&O&aoG!;|u&k(J_1E3Dw8xB`;a|=FIVuA<>Lb%28;tp3kt7n5lDA;q z7n$WkrlyQr7d<0Z#v7QGGL%DyP%c#1+WT^%Cx~Vpj)ncmv6C-Z4E-fC!c5gX(IMAW zS+ye5&}bo-uFwlF3CY63*fr{h4-XtoI@3dM5M95GG0X@rKwis?=k@Q=2g|;&wbm{_ z?nNvd(9@$@>bh-{9VT@By2J}FFEI_RiRL(+GncLhP}+6q&VhBW^7Ho^yQTMs^=q}y z7As|{KGt?!7NfYS3(xM$<-!;S%X8o&J=O=8uFTUIe6PKXP4j90_iMVa&J2}F)s`n- zk;zh_L9;nk+I>Z4%)V57VFLCukejeI+%ajps)xwYMzO8_`2M{%TdQ3E+bg1c<5q^= zVvDG&SFe)ad(uFGicb5|LdvdQyR|$@@pMS~zTru#c2z=RrR8*;?SF~6J$(h;rI9aO zh&5w=w;j7{lZ|)EiwD>qCwQN~R54%N?^A7BWZKCU`^E5GJ62lfvvbsUba7dMYp_wy z2|fSg-5S(k8rvvCNK72g9P1?dikL1ovbPpvLIV1cMBp24h0l7}Fea#ei}dTu{!$im z*HNI(GOC99-o7v2P{2e}RR2V?`Kvhsq+#JN9<#1m4J%Y=7GrOBmPR)L>ulP}Rm=)l z*oT)0Hkz9ZGfo_M^UPHyrnkbP+PS*-Oxjaq!(PN)xpL*<>@Q-Gy_PAjCmR00KJeTS ze-X#8rAL|DY0=GP?2(|J0u1~&ATNRQJbSnD}E&&<%jqxuJI zXUcJZXFd7RM^P^46@7k>g9;VOAMeS+T^3Q@Q+iKMU1_W2)p7cbmlHPLz-0p6(=+gQTf^eP@wj&gJG(xb{7;b0W?kN0pJfi`&*QcRaA@ zkjUbSDwHdS8Dy77-ytU!&i-h_CXfnFC_MLLWv5QL>P*w*;S7q7@5f-`4=@QEmE z%j{>z=6Vab3>(So!Z~$P!hWYTnkGIpO*S-5eHAIGF0QVkZlvz`&$SL;pvRigSpU+6 zh5+PyJUgtowq$AP{$dLz4i1hd7y3%ZhYvRRKs4r?_Z3QhxNQ>jX}gZ?Rb=0%kLrfS zRgbc`v+icy&yxOk)1b0k_eJR5_`B(0*Hzz;>yJgj9+-|gDsFT&aO7l} zNUoLesy#8eSi$MfuD}yl8o8S#q@><0*A{KuldPw}F6yl5>jX|Q1QzZr>lElJyY|n* zzMy8Ucy5%1zR!_oo=Hhp6Ps&a5mD8*CW_NYd%t>BeOw8xP~Pcu`zIO#RzWn_ z4n^?%z4~IyLDup5&|qP=9r7&vm`C2PUpIK4zIgT4t=1{WwmQMFy1Ur(REo1wvof>4 zB-==n%H*diKkgEoAS;}H5|;ou5zQ|B@FXW_RyaLYFEqaL`0?Wlbaaanq8jJr;)J#2 zVg#5_E}IUwi6dRb;b1|z*%p)cj|Tbd=vv*xHYxsk6MTb^1gpvYN!T0oopX?p-_%7e zzghjjV@55Bd+Jn!Sq{pukfCCCp`51a{qCD07{DEe&G(n^(K?GbE<|T}Z{YDeE$Q;RY-A2NPsS8)*O5fP z)HlnBrz%I1JRS{}6l&lNcK^HD{{%n)R0W1jK`^{C*OpBlHP_}==unfvrxv_7*x zA+DWH%j04DVe6o^TZqBQiN?s$t9P9-*olqX1`tHn>x($XaMe=qQ-vn22`}wtyP62J z)y6wb1R|O?=KBrmf|%A;O% z{l-y2CrrhEvSY_H#?s|KeW@BiwL^^!F63cXwpIg{IEAt7hud`#Z5b*#XKg<3PDVtK zn{;JHQuEt8cHinb9Rzi}-TV^@e`=Y(Q7`EeN^W_B|4W5YCix*=> z-P*yR$bAu4oUdOi0T+MY?Iy9wG?+RUgg;?V*Xb7e_q+j>0MypBwl1FlP)mW4yh(p? z)>UTaektP?)+<-;4Y}`H^jGM^L6NW^T_!(WX%KD)mQmVN`x}APaiA|h==Ik53+KH1Rsj}J2da)J|yeV0p1g%n*6MW)UC1-G*X#-3DbbWndiK$vrt3gH@-1SqVG1;b0m%!&30vt50u*Cip;Ck1&YHcBf8xb zH=x-zE*i{08Sd7|eW`gFBK*tDL_L>)Ec5@dR93&6v<7-~g zGM?NkeE`Hz*~}~O9&cBZ<0in!Qjj&qLD3X;;px+-`5>aVK^w!qdh0y4$N3gFwM~=3 zlB^(FDumBep6E#sz(8YQ9N)z@e*Bf{aTv&uD_>60^{_s8_>hv(n9sDMmY0))WQiGj zFxB49X|QHKi1r4eF#zKnp72Y6qkB^%6?yy>P|hw4mM80ed692FtE};oOvm+cv3c*a zDF3qrTLnnSvHF>W4A**E_1>JI1a1%j33$3}{fi@e15kYEN6lhC(-F8j*{mVu?d@F) zYE6T!rQc!O5pUq}1XGpSyA=z&%98#xLa)|T8BTn_^MAIJBQ~DE+3HuDMXWO78Fk`> zT$I&|Obwtr!3{z#td=dRq$}1Ht0LHn%({Q^(otZJh3x@|z8T`8@4V|cm3yoVvIkd|3)s%^bzW?}9huR!6i}Zg>@L3k_ ze_DuljXX~Z;WM3=QbT=$c1z-*^F2^CGv_`ILk_a3zM_n?l(%b@@9o|cbl>~DOrB$5 z$VERQiajqZ>L1u|S>9Bc9Gb?^<&Jqr5Ki4HFlZoc)vINlzLi#!KpT7#9MyE7_}jV9p681}0o}%kTNf^0 zHt3(pblqM_U{vNaZmAxBc8Lyq5to6H(N|BxpD81W9ljpXp%?tu&}D#F(8;+{Ond^K zbG=V!J@%e!yafGP{qJwDnQq_ifMNqquZInJ;Y`P;CqE+(tlrT^A$)>P?p_DF`-Y-# z#m7zx?jDqk^(s1U1{q5SyNe9~+ALNBLFR}i+}j%UPpApgYSSwORE|S$#7T!S76~Fm zFwy7-8cPBKksG&E(;i*v1zp`kHJ{=8qN4hsZ06FhC24S#d4L5Dx7O39C*#kPQRuj+ zJsMFoo7iniA#+k{+DL)Wm#fRyL}T8PliVGp1 zB&q#?q}4-YhO()(+1x`k&g5m(-2CfHOu&~a6{*Y#F}~z98RO{Veu`74(wwC=UFinD z!(7BA7KQIzbbx=Ncfnm#m6nho5!%QMqAl-poa(o#vIzAx!5}3(8+4paN|wqP`kvx6 zsJUddJ@(gquB9Fvv`KT;0YrFGwF>NJBMGvJ?<`m zc@)4RG>9a2XpF zI&!!!sa81UDFVjmSIph;r?fHg?6rVEOtzbDJNEL@8}T(sv{4ROg6rhndRqR@C$&%j z8A$A_qA~Wz*aXf6yivN`Vg9>iZX!V622dSnWs5DJU$Vd!dbqAyB!$SFy_HiV{P$EA z;J!-Pf9)J^&$_f40LG2`ug>Nvkqg+Rxq@b}c)J#0l{062O*(CXgQUNk(5b(mKX8el z?2lw=sc095QChg1rd3`H8vmLL)zouoZo3=c6jO09fN5O36)HK}2WnWweF`8i;I(~v z!rb3A2neR2!E#3pj`wd*uHE5u3+#6JRkd^I%0U0R_1@O=)5OFjMDC<)s*j#Md*(f2 z^eRLKUY&}!*3+q{>Ti&dqIM~Bm(J|vLa=$ERL zrsJUfPo#u^>KR~1b5(v5I4ujqkAlchmZYwl-A}J}E z7xH=3@9r2V48xouL^@!KuPLPMSC(uhLin}V;{7!=3otZ4Q-jTMdjkxO#KuXo8m=rp zV^y*Isxelm|0wHJAA)f0(i`BIE){L2I6*UdW?XVy7G_1YYnZq)22ASS-fDB$nS1x{ zahY|o^^-meBYc`LeEKFq(6EwPvH4OozipCPclI;BgBkXc0R~WdHH6vclc3`qUVSMZ zFOB63h8osagVPyhDSisFSSo;f!t$+#gwDj5~&jjALm9SQsszTVpU3)P*4aO-vi}i_~*}eK%Ol+0wMw# zY@v;H>i9_*xU&mFsZ*sdPf%MQs%Mn(BIz)o?L*D4&{u46129p=AyeHsahckvJ}zQ_ zLDe)c1V@YQ3#B&Ocd;64)b~%)Q9e{dfEov#(i-q~)}|SoMGh;nzE&Nqb5d{{`dA2a zfq5xFvFP|o{PLtpa-r^%;J~&5Fp#h^|_tYw6nJ@bD21 zR>!~)mMym4u{SbTz`>iYu-aoEVLG2 zAbbsi*U0lN$Nr%gx<;Km$$D~-fIwhWAN?)#F31Q06)g#P05ITT4?813E-5Pd0`!2j zM;Co;kOat2KuxDI%`e@obrR(NuQ&=PoEt)Eyt{nkT8%$lNl42*2|!g@j1@K3f-DTD zCxjt{+iZPfZGp=a+|G&pc5*pzC>~9CwtO&`W7rsR#eTLcbZkr?^gpippd&ekD+2vg zWUTGlz{N|KLViD6A4s0ODHQ&DOwNBAQ@|cA+6kQjF z5sbY=#=2;`%*z6NZnj<*5Uq%9Sc`sma!rcWW#Z3$&+ z-6F|H{g-tf@*H^rl zGhO=2QkSpc_I(OaU`5JZFr0B_Z#F39j2A_0?T!M}Z zdc&p=HAJ>|=u8&84S#C8qs1oVugb*DH38F90 zfE=%su9$a}G~l2wU|^7{s}J2xex67g7L_sikmicE_&a8_&EO-Dj$)uJS^vnrILNj( z`mK^J2g{1HqAW|tN#BCh2HJ;GO`kuENI3mz<8WGTL$=`EMUS3*)jk+huP02(oNO%x z+FGWH)KPf^wgwVMPVC{@Fy|4>Ic@YB8T95GprdiCpnBCENKkxy{Mq*(C#8gp zW=c<*5l?22kd|gbE*$Okvt~S?Y*&_R==xb@B5rD;2^o)im&0U+I#kiT7Ues_|eJ_fWb;j%G2REjjkd`at@mDrRF_P>3I9ixktM zoz`FE;Pcvif4B0zRWLfLL-+V41g;0$d@FzAOU&PAOv(#&WBRpf50JfgtOoPO(SGpx z_bV!W&_WiUIXu{`G5}6}Tp5GbdSp8Fb?w%3Zl5!P3Qx@Y*qF*)V{Oz=BIW_`1c z+H_!T00<8rJRrxRP`fpyE8aJ=3bIsD{ZFWvs_ul94E0I6@0dhg;-q~#9q$E@My2@+ z2DO$r0j?h?>6l}}84yv7Oj02?cUB4$u7fFkw6L3c;jsLX36MqbaVaH(rWxZ;8hrI! zi5rQTlOw=i)Q;~J!qTlUV4~r4;$r+G81p=f2L?>n*wqUIfS$l0LY(=w)zAR>Gz!$} z$DIfP{gQ|$Py|HojtZ1!?0OyU-inb&oU2^^NhQ`?cd|`lXq%7a&6&{GDBv)6ofEYd zl=YRnS&$KM&&5sR_z-YKIw;gt7hi)AdRCTFw?xx4EiUr$%`3rKvD?bRz7b>NJg03Joaa~gulZb}F z3~Z*{{#zJ2e00G$_W~^~$MW)W*o7UsHNfOaHOC4aab?hTl)@XrEXF1c%9Hf#+QPV} z`PG{{5x+28FY^p%;U*RVUkGRXMiUFX-keKf)(&Biy3~GMQIkM{-#V<0o z*Q0J*c6_z@-N%_0Td+RxcXJoP=3MN*l44@pcJddOaSOf1OF-vHn44$LGkz#0WH_Qq z0S5(!id=d>xt!Mmjkp{2ol~0zSKgZ^knVTmUW!=v%;8Az31j;Q|?% z5lf(mhr3#dWz1RXJ1yw?1~*?P-gdIwU^$^lp~A^~{Oh`A_N1Lq38ovTUj%VA;YO2- zhr!pWME|7cqd$+L{39b#oPfO)nCqfbB`6D+9~GGrsTByQQ{FZT4q&59fBUv#Wn8J~&YE+* zgBZfE+g1pJj4N^g8Haic4AEwz7D6a8FfizYVy^fM5L(!FaHh?_fVyD<`hYB$$6>$7 zSX`W3K2&Y5gE%A++B`YsCUnY5elrNDAgh%OyWiNkv-|cUGaHlL*l(cqXaLi4jH}Tv z;+8p}E!hJ`_G?!Qm&w?$k+V>>vz>Za-k(z&ywLC)k<<9rK9ZRiI8nE z^U`L-l&=8=zaTT{P>j(X2%5Lh2#On=R#quMxiAC9x423A2#~N3 zNsgj)&!Ea~u4T}{YhC?g_t_JhLG4rYbAwN4oAs30XBfZ#1;&*9tWyMwQG8YyU;_#j zav{XKvot%AxTiD@yQOm|5c4f&C$s5%mFgYVXw>X8lcS{wa4V`aaw=yi@a{k0@Uff~fKRJ@ItbMEQHPkQa^A5;)Qab@syt?y$u2`cQwLv}%}XBfvMeROV+Xh5w9 zc*vUedZV-~QNRd}15Tz55OC0t2Xab+vw(k`tuHG&m8{0SYg?r)z4jrYn-p_ceJv0S zT|{ayG>8K(n++N$0QH&1TTmZJ%gFcu8A2a4d;*7n>z+2AJ+w+?3-iPo4vCMaZ}>XR zBbLLC@WrT*;a$)MoBT^LfeOTACf(UF^Vb*qHv7Ru%HXT_X%4c)ts|2(pK}sRaWDVg zK56t+`Ya8TcFrwC7G^0Vq&P#B=*eLMRm5GhK8*baP>K~QL5UEG+=a?rOUzlVU&Dj1 zu^Y^PAa71w2rS0yGLr>q3Rs1b|$?STvTbA$}LZ-h}qVsr*7 zpjSn`^Wr>qm?X|zz|S;$$h;Sz^hq%7Ot+b6q?iVJrX9jjm80K=n|jhVbvj6tZEZ)jy)&Bh#UaQG5YpKo9qQTmUm;lyg2N|R9H6}`0zFY$o&eda=kZ=4Im1O`J5J>}$(HzWYq4Qv2GsBcJx|Vz{(6MWk z4H!jfhd*G0W34q`+dG;x0fo(bGES)&0egLbHa6q+mo-lBH&5)W&mJ@K z(0W$5!n>1KDb{mxm^=<1HT153Nbcr9y_zoa!KBJVfHXTIgu<)M{Rs^aDB0(Ou3$p0~B8t{XeF7WOhGxnNgl>5zyB4vB#!IipK+C8zjJl zv`<+zT8@KEkn$7I1DFuRg7sb@IsaISl(UaPUU2yZD5=yq6mkp@vSwh$kEStzB##t0 zrR67+`lBIBKz22t2_J|L+OLq8pR+5pn^uUtrT4I`;SXBvV52@V9{#GnqO=NuFTrY7 z&Po$4pYt4%#dz%EBZkmw*>3gFgaULVfbWg{!c-< zc4seA(6c;Ts5*Sw-#Nu2mt&7;!Yro@Er%8y096rSi7$|nULz#qC{!XieL4^@i4A~J zJh^4hflK7V=7lm*Wr|e@_E9uO4&v~PC6N84be2H0T# z^GeV0|1t9fSzbvh%+Ap3PuF`o2Sc<4P^r(6(0lBft@4O;N zvF*PcB;t`ae{IHIJ|`5>gr*os(je;=8M2T8>yA9bCUOb`bNyH6Abo;nUTB(XUc&d9 znj=ZniQ8{^Rckp4U9Zaa552}$SU9MOQl8;?k`p6$U~F6t$81>-EbNvk+CvQ;(=+q* z8c;hqjR}D?N7Dbx54d zG16{jz$bq>HiR{aCxJ4ocAC z*1tu8>gX8Q;PHe&vYndm)JpnEb6Y@MtkW!8+hv^+EupC%1}5Fn@XbRVo8XzA-gwz; zHef`mbT~9Wk4PkdDF>9EOol2-6ZiCwFQS?q)|uWJ6c;`pM0N_rw7v4Ohg2g6g&^uO z%hL`?B;MVonoePac%ATolvJ&j{3hV~%rWfOJG) z?kT<+fVNZ#;e!VuL3EWmCTBw~(P|Hzf2|$0PMOthims9>K@#Z*H+yXXmF1djxrR-Y z5Hl~h-toj{x{XQvf>xPbAYfgA4jg%MbS%#dJOR6@Yf*eZ8WQ~lXPw!EhdVE}Z>Mw{ zqm*ZPG;#*fVBq=~(14s>XZoRg_^b-wy_5k)YtV$>k5*_#A_3&7pQ8()nX@_ZR(f6b zrJY$iC))S^=V^F`B-_l2waeEj-Vw<3>6 zsFc|vbi0V;gtWFV)pe<$j0to`Q?2*GBZsS1dn@ntpsKL}ye~*LYL4}ryaHVVgwXoh z@dgdN`IZ(wXQ{VRyuHU93$dUjwQ?{$-U?t>jkUY6W+Ap#lMNENGVrAE-{rN| zGN65|1Ike9=KBwiFcHkC%c=`-v;%#q{v^rHp7SG1>2^YO&GN*d9a?7Wq&q>^H9zMz#@tnLnD1i(5;DJuD;PbTehxg4K}^uYd+(m%es2A0 z+OEAHX?8aPH8ESHvxtqZfW{y`a)8vx?KmJR4-Z^ev@1RfZjS{+3S03DzzE?V@4tpB zExWm%R-lbMX;v=3@iAv>@EBsE^AK#cwQ$^cEs5Hg6+#XXP6=V7BYfZj=x8+BGEYHK zkyl-4(w_9-?S(d!izC1NObFv!P<>B>2^ZA$YrAh7{B>ts84B8UZN5gOvAk zSy2C|#RV8$anWg|?q0rI|MTg=@~W{0vOmW6*M{JW=hw9FEmB}?69V;6K?bGAwt*PF z-@ktYK_Vx3e^99R7AVKbai~Ef+q@c=>gqSih;frK6x`Z6w_<%WFyBP4!fwVXe5##V zE@6n)dQ=q(nS8fhhu%Vyi#XK$ z%9zzFS>iWgaUKW0d->1eZd)LNM}0e7PN816Y*Vn+-J&|m#<p0!x7W1mNsfICHn>Ip-6 zFBy_BwUyc$CHns&7T$tGFs&-8;VK7aU*x9`KE9^45u`glH2CR5Jj=3AS$k{kyv?X; zhEIU_#)l&a;)G2G^Z~dS!2!6{$DkKD01JWY>yu7ZD2hciot3{kH6S;ycwBjfdt8)e5Z(8Y3!YpbGIb60!JklODmDuz3|>g;uIn zzXb__S!xWTt~J?~n)rq#rrVO4KqC({m4Pg80eVv$K87SFkV-i6LwFFbQp+o7l_p?5DE zu*~9Y($HiUTBwq(-1@`VtT$>D&3>EMM(iPP7~N!MD1U6Y`!AU`mrE#>Ru@u~Bf7~A zh=sRYmY1@ERUNkz^7I@WavPBYhk2im0DIhMj)ZZMEqR#Y$W7!(PKCb#~zOo{b-P zcn0$yf>f7xxGe{?F<}dKIy=u7hhkQ~#=Ga6wD-waDTC-Mv>a5g+?l!2x_beg09oOL zk2Y|QA_9^b`ij>=JR^7$QfIgHWcPPprCa$Pps8?0dgGv|0mzKwIT`2CL*P!zPRRjv zDL&iJ(ks@xYTEL}%=FCs=;Q-pQ3CbF^HD3hvs@+h5Jl0r7}ij)v(YeL;&=;R#HI{) zo>_6xU@i6`oqx9rHF_tQhrHou0P7?b#cSbns3^EQR}ksmsHLGg@Vth`W?_IIdMA~- zn2q2FP3k#zjWJYptkB&gh2&ZNG9%rO;X0J8^Zky&Lx+o!k%G>#f*U=nKk;MlIw}Lr zQE$bt>T8WhfM(j>jF`^X+m&=qTEt>eOw0Yoe^{I1Uo$J@=SQ2fDTr7+p5R-j0)KVv zx|&81dL+Pf( zEYTttI@BigZf^>8JpHfM0^vWXG#o+-0*BRlq#7HNiTaX~2*8nk^|8H4I8?HcTw~}yH0Dkf8>xm}I zu-;hDxWc%WjrnrJ^teat1;uANT^Lc~Z7xhGqc#`jBgcU2LgTzfM6x2(pLp^OyS7iQ zRY7mK9OhK{LGhQuv-Y^)6s=SB)UYiO7_Tixn%pxya#d)C%&MIJtn@1wr;zp8B>=eky`9g<{`Z`0mJMc|H*XF4Q(6Cuvg#v6zIzJdL_^l#J+ep=LwbO#$%m z)}N|94p!+(7s_TI&^>oSHI1lXVO{5jJ_}$rtAtA^KClkNF;AjSGw5||0Y@o&7<=vP zBHAB5!+)IV!hVDjzpl38KQi+6V&vkPnws~VJ2m|f-DYn_T))Zf5f~+?1)I`ah?Etu z4jPjF)NLQ@asyR!EQ;N#0uiKb9Fr*h>20~ZGss{FsP?j zRjhK%cm=E9kjAu77{0=Ktn?!4IO+d80}DSR6U8e_$!l&b__gcf{Zp+q?M1{ltQ34H zd7!C|&B31Y_DTb%Z`j2lw%5*>U%!a=jF_%);0>ts>R-Mm^AQ}v$7NR;@wy1qj$@@p zI-x~o2?oFBc-qoCEwqb(8?d7$+gv#3AHo-Zti>J^G$SG)1ol(2dlZgt(jBt3>K!2j6l{y{wWhwYiL9EgSNhMnwz-PopA-5~oqeoDT41 z&?-6X=h_2j6vd|TTq~W74F-G$<-R--RD2$KBQPTyvNM zdi%N0R5${^INr#I4T|(lO9XL^@rq@C#&u3G1@wZNtTd)#AZCDRj{qrKTcS8I3m~j+ z!KR_ENZgNlxoiMot5fCf7J1h(8@$(O0GdlCF0QuO@mrYsY^%7s)!S)%Qe;@(k#Vy5 zP&KD4xd#eSTv(Ahh?;;DFh_Px58qh7Huv9sFHRU!FfK##qG5r`Nh`b{gp!d z+iX$;2KcN9NoHe7$t;qre}DZ+v<_bHa=2{FEtiSN2j_wL5hO|?_phz_{=HVuU&v4n z$i$C36-|bgGJ6+!U2Ykm(P4{OG-fTDy%eg^ugyv4vFG>PKmJR_?{8HB8kVJlObkP0 zt6ds*>q?%GHVg?F8T9PNNe`Aijujt39p)5lP^E9WA{JBE<2PR*b;7ToHtlUNeMa-DvOot?RBy?@^S4=~`H^=G^GMx~BKoVeX z54~Vh)4oP{42Iy&iarb{y zfEvxJ%b>D)eU?T*(S2u{c@&>_-a{~7{km;P>ydQAZ#S-P=1&A#*M|#xM9(^dFQ03P zioWw)#(m#D)cX}_04N%ywpT)7xrL<~S1n3=GrIMDN_SUS|0d>g|9V}m%01brg>mul za~2NE55P`=U+Ak+>ALNMGzJe%eVJW>UMIX9&%gqpFB}QwkCu1Is$k}~dL&2$>~k%6 z3c0Y!{P8*opb{AG-!lHNKL4z&HrF{xtLStfBUt58u!ElHL? z&+uBP&rW;o)8)JVcm{g52Vu(T9&mr4EP|Zn{3C4uH+DV52B#MBFQ6qVGHguk9nozr5p8FZH8D0HIL5IGx*=^Y`?74|_4X zG^p+u#|6N9EzXU%x95Y`>jZ7LT@n1(z7B37@uth7TdmVHFC9PSr1J;k5TmtCH6V56Xr}84 zE=ayEP`?Kx)Z^ugmE788EnwL32{4!gHVy16pt_~{(rT_oK?pS;dy0UXiBqPI0uSRN z$y%x>0RNy5G?0+aDbXO+FCn&_g^jJzNyY@s6Jn`GTx=2cY#65nTRDYSAHu=xrM*?Y z+}5N4S_gTNIGx&hC266}=i7y1UAsd#r_Pac&VQAB;t9ihpz+R>PZ}tj!l4JT zo^6Gz%g{uVK+%(4K8`^JR%9JQ2-Kf9;$2Zw{43&r2(K zT3dI)2MzH2RHwRa%+Q(9{cD_E*S?n}*Sghw`}ZzO{y}%>%I@Op2fOjj(@9y1K_LbT+$u>zlj!RiG?J>YQF8sL)X&qk50L zi|4R71U>img+5DTt-#OQ8I&OAZZfr>GS4X_xj7BlwaY>eDSmT~$&I76uAt4L51{b8 zCaeCgvJt$fklS{cdC+`O?!AzI@nUF!QFBCa{=OE~C2CLyJ6V7?>Oe1Hl*tG_IHJ%c z74)Z)V91mN#B(d)=`GNQ6bk>^kka)$5)>3-H7U!H0*Ukx!PMw9SSk!WeW~lp((o1U zl8?Ug6qnz2HgA2cSi@=wjC~9Rf4MgcZPq@K7h=LAB$eU3aoj21fn?qHujdja_AW=ZFic>i`Y>q5XB{bXDZ3=C`oY;tQ?;!ebo z9->wnu%*tQn#;~xzePG!S2`^}Q$rz&=K6}H0sEt0_(u>t%QaeS7;w`Ze5YyEXs5UV z2=4}HYzSt}1@*t6CK?Ny6I?pN^C#QM^iq?_dS*=$d>Tp_n}&c+v1LE~^XGZ7;*$kN zh1(g*V?W5;GAkif*K8|jlPK6FqIk@%0qsP&A^?^5OkB{jk`YZ1d$u=h7?9OGu)b51 zU`l5I6p67_0#-qgga>G4jAvzuEN7a+IQ2WeEQd>k*WbhmPp^R{uB?D^~unH~llI+Yc-E&4j9;@;eQm;L8A`t{6F{vixu z(UL<5_cS{tv!|cS2NOEMZToMgy;UA1)IyYT89408gQWYJM9gb>?rt%H&TZ|@)7U^f zuwB&qfUc;!D{1Gmo0F51>H9Vy!IuQ2MiOAY!8<6i+%T_mQ5C_D0VhMU14L~u<#4;V- zyft;G(Fp&<+tRX71SPTG8`JGNeff$wDbTwbARweVly9k7$|+wE82xmUDAqJCc_)M% z)D1A8O*Yo;x{-Va3mcU-`(|=?b1{5sYO13;8J~0g{3;+uoN2Wr*L)7P#7Um4MzM46 z^8}cvSRh6d_}D;ceR-Z!Lc_M$YIw+!H|0oag5Eh22O1(JH_&X7=a?0M(S@Y)-b%1} z?@Gqps*l8Bd8$V%SR05d7T2ad=6Pow-M-4|Og6{yQF(iAzEt#wt~q%>(I98)K~iJrMXz0;`I~ zk+D3zng@rJpE!vWXUg$ei*#gAOCLP)I*`c}g=>?f5xq(Vh~UTZaBT32K>?uprSf>$ z*%j?Lf=M(XlUVa7X% zm9f*!hi<>N$DbJX#eY8=t^L{hOFP$wPwIPQq6cgFAh98R=h>8*#9c~y9SQ zIW9OT$R}JI$(oM?lJ}WWvquF&_ArHra#}HXy z+?w!nbOw6=8D&`}S?0KCDTug`i*t#f7v-D}s>1YfDchqsN!+EWBsvdWCLXlWoAres#7Wcb%2U~OR%cYMy%kv1 zqKY-9R~}wZ{-T;CJ+IC-Z$wl!V?~+P8Sc)fvvVjw)bnB39(5HL`XC+Sb6q|MFQ>rR znCCztiCLIC-mo%?`CDOoWAg`QlSBny}7p z6i5Fa`+WJf3u4j4qrEw)nxy4@j5d57mge>xDxi$N^n$UuUXp#xUmZOJsE((n=w)K= z4qH^T-6-W1Z9RBcVzqg;y3Dcs`vbawZ&h?b-ri2~I1wT*l0Li+sYGrcF0G~ROyp4J z@&s@5qc#gY2MgPIKD*2p7FH$iXKptx?!;~s@XLC%+s&^PIhziPC4U9-r20Fg6Mlg>i(oc@r9{T>x@Aqlr7v-gJ>~~Pz+D*PPgm828J)eg z92bN_mpCWja4meYomQLpq`}C}#z;@~FQ@@@z>P*0OZfPhn^$C>2@kOT&i(K!Z}ZfW zi;x%hd~i>X7bqwff6~kLeOR081?kme5~ob;Q>ed)Ywu0;!MPa z_6Hls`fSmP;d$si4mL+Xgww-kf{Z;yE_~bJpD%ak-L!D)u1= zYcBXRDC4t38Qsz}cj{cbn+p}z%S`tSAEHFsn0fZcg(XMH4}$hm@Y1S~ANtjRtPwni zjoS|0q`q|1%X`OtzG2Q))vL`!XQI||qMYJQhogm@h^S|zJ@|&3y8KNv;4?-AcuDSE)d>eo_jM)`2PPN?x{KXAsuL$L z3OXo1A0=r(Y7rWbDBC2Y&UQt1UX<3Bma)}gC#qCRN!@f2*-m$vQd+8x(PXQgOjN03P<3Q% z*R85CMLKb;ahAA?6q)xoa%gA&gXIUW?x`VDa$Y?23P&-ucte z@`F7-UDUbp)nNNKtHC(qyFG-9)Lb&7fx21T;hbNxSyFb}Dumez%#|D~Epi2fvhh}a z<%Yad;cx!Uo05EonA&_IIv*QBh8F*utr*xT2Ljs!VVBg_*U3+XRsFkRG z$WulFTpTO9Z|yK|5>6fXMN&Fq7E>mFE@DTPqRg@A_o$USUhs>ds|Q@cTsY&T^%%STtQLP%tX&f`#ctlq>vtCP z6sRE}=y(VLi{U=w0^Hu-QryrQCiG_;#FK}smTKA^M#{=k0pp5{u`6I=rq8O4$^_Gm zJ^_Cf&LZ_l5W;&kOIuKh)-V!W0#5=}-YJ~ulS!n~O$l_A}3$Il;KRq~^Zv{s$a^dTyrcIY&V6NN&#LW>v;SK7(Mt~}m zY`CDv9|OdPMsnMi^wpDaN028Pf;yUv7ZaJ>d@q)eRIV!;1^!QY<|yFMWtHt+>}8&7 z8}Xgeq{W%#@qEGw5aEphP6*&z72HTKWOPJ@J@+}+4Q&~=rpvGBx0WdGUev- zf9}(Y_MG=1c9QKI5po9*Jum!{YZG>FA+~Pe8WHsvW$0OyRkyL&@e1t_JM!y#}wu?9b^p>x_sKlQJ>q$aw@IHXo zNfI(2F(4@Jc$K0jIU`x=HUQLzb}%Rj0I<#doU zEWJqrhY7TZ={dwI2Xqt~SMe8d70jYF#BY(D$T6~1)D_NoRf?h@mWi|=>5{6~4?&*O zix3aZ;@H+JB7|VTy*mIP3V=w_RKSR}x*VSPxX>xRq3GF>oSU|k9&n4qFE8cZy7dVJ zu`)hB8m%D_R1eN-$IuC2S`rTT_2m5}p4FkOhuQyD=l9d)QUKsrU)kyN#E8he4tMYc zem||_5g;hj6R8-)Kub31u5CnfY$3!ju=&x* zpK&OUuw9X$w>Q7)Gft>`>V%asdwB;VPQbKpA+~mR0hxV>sw521JgDu!IrflwX$&gl zc<$z^JI!6cXNpe72SF$1h3*yT3{!EYbz!c&nh8F!bKWjQCuG3sb2<$RB&L$~ zA_+aA_LLMXr6D=>rBZ|FZsva+Jv-V09Q3LV<5LuDkd8QorV=k)2T`pC$R;7#<|ph? z$t~22o|-OjQqW&I|BtzgOEz*j7i!Q3q8Wg#g+ay(kG7J+dX&P5e52VU=W#+`1GOtNdN6otL4flk zMIs=Vqr&DLE{&SaNxX@8EWjOc@K zt9wIrWj>7|OiwIhASB}w8Chv=1WbrY2@5;@M4nEB*EX%+Kn;s0%^}pouQ#2H;J}9C zt!$7ezhf*{+JhCwF(85XZuElAyE2on9Oms$*eamCbg3Ov+qUFEKJ%5IL^}hF^=Sh@ zAfBL;_I2ej-aSiE_8`82;V>PDsqseR8sceUz_q23*S8NE@XJ_Wi1cb%j*B+B$10%3 z$4xII=Aq!0n9E^p91eW0jriLxSOOiBT{F5r;=r7lI&iWs7tz6$%iSKa`q=K`HSR2{ z%#0PLEj5p`lz!;YgjZ4=4k$jBInZl;STp4g>PO{R@f1hD46puV+DfC`IXrVx({(yR z{mM^meG@sZKA3~^n=KuqpOGgcWKn2r64;kZ%5J#&hJ=Lp`1qtw#Uu?54&Goe80L1> z<>jBU+3eOo$ zm6g9J1|&QpNa0M8lu^mD@$c3bllA0&$IlCeo_8JI+Dg0>8hWFmqM|x*Vqzk&aN@P{ zALeFeXVVG_3Rd0XxYrf_XlIwUva%8u5wY|qKO&KMi9jI6KDQSXns{LTz?!|>x(zQE zs0f`Zsiv;p0{)QwEY`=2cZ~C#lFH_YQuJVhoQa7^>eA9uSy`F+cRd>` zVR3QmoRVrNCDIV2Mzwim3i0Ha{Q+UST&F`3W&%1HHJ0_-iBVBeZx>gIZOc!%jdVz- zrnWY(emNm7ZeX>Oi+fU0->6Umq3JeP`0i3L4na`{HstliKc76kh3YQkBKXexw^d%<;eC(l2ai z*#GtFOUbXqdcL0~No|h=ez$5n+IhFv_~TI6Cc6Vt*cAMS x-MRL`(zLsbe7NfFE+f0k2(JbHU$^YUJypf8RgcthLdD?kxTC8>$>+X*{}nh&t49C; diff --git a/docs/app-kritieketwaarden.html b/docs/app-kritieketwaarden.html index 04a14da..240dcb3 100644 --- a/docs/app-kritieketwaarden.html +++ b/docs/app-kritieketwaarden.html @@ -52,98 +52,11 @@ - - @@ -382,131 +295,130 @@
  • 12.9 JASP
  • 12.10 R
    • -
  • Deel III: Toetsende statistiek
    • -
  • 13 Toetsing +
  • (PART*) Deel III: Toetsende statistiek
    • -
  • 14 Power +
  • 13 Power
    • -
  • 15 Variantieanalyse +
  • 14 Variantieanalyse
    • -
  • 16 Chi-kwadraat-toetsen +
  • 15 Chi-kwadraat-toetsen
    • -
  • 17 Andere nonparametrische toetsen +
  • 16 Andere nonparametrische toetsen
  • Bijlagen
    • @@ -535,7 +447,7 @@

    -
    +

    Appendix C Kritieke waarden van \(t\)-verdeling

    De hieronder gegeven kritieke grenswaarde \(t^*\) heeft een overschrijdingskans \(p\) diff --git a/docs/search_index.json b/docs/search_index.json index 50da085..a9ef7f8 100644 --- a/docs/search_index.json +++ b/docs/search_index.json @@ -1 +1 @@ -[["ch-centrumenspreiding.html", "Hoofdstuk 9 Centrum en spreiding 9.1 Inleiding 9.2 Symbolen 9.3 Centrummaten 9.4 Kwartielen en boxplots 9.5 Spreidingsmaten 9.6 Over significante cijfers 9.7 Keuzemoment 9.8 Standaardscores 9.9 SPSS 9.10 JASP 9.11 R", " Hoofdstuk 9 Centrum en spreiding 9.1 Inleiding In het vorige hoofdstuk hebben we geleerd om observaties te tellen en te classificeren. Daarmee kunnen we de observaties van een variabele samenvatten, bijvoorbeeld in een tabel, een frequentieverdeling, of in een histogram. Vaak kunnen we de observaties nog verder samenvatten, in kenmerken die aangeven op welke wijze de observaties verdeeld zijn. In dit hoofdstuk maken we kennis met een aantal van dergelijke kenmerken. Sommige van die kenmerken zijn van toepassing op variabelen van alle meetniveau’s (bv. modus), andere alleen op variabelen van interval- of ratio-niveau (bv. gemiddelde). Na een inleiding over het gebruik van symbolen bespreken we eerst hoe we het centrum van een verdeling kunnen beschrijven, en hoe we de spreiding kunnen beschrijven. 9.2 Symbolen In de beschrijvende statistiek wordt veel gewerkt met symbolen. Die symbolen zijn verkorte aanduidingen voor een reeks van handelingen. Sommige van die symbolen zijn je reeds bekend: de exponent \\({}^2\\) in de uitdrukking \\(x^2\\) is een symbool met de betekenis “vermenigvuldig \\(x\\) met zichzelf”, ofwel \\(x^2 = x \\times x\\) (waarin ook \\(\\times\\) weer een symbool is). Vaak wordt een hoofdletter gebruikt om een variabele aan te duiden (\\(X\\)), en een kleine letter om een afzonderlijke score van die variabele aan te duiden (\\(x\\)). Als we de afzonderlijke scores willen onderscheiden, dan doen we dat met een subscript index: \\(x_1\\) is de eerste observatie, \\(x_2\\) is de tweede observatie, enz. Op deze manier geeft \\(x_i\\) de \\(i\\)’de score aan, of de score van proefpersoon nummer \\(i\\), van variabele \\(X\\). Als we willen generaliseren over alle scores, dan kunnen we de index weglaten, maar we kunnen ook een punt gebruiken als “lege” index: in de uitdrukking \\(x_.\\) staat de punt-index voor iedere willekeurige index. Het aantal observaties in een bepaalde groep geven we aan met kleine letter \\(n\\), en het totaal aantal observaties van een variabele met de hoofdletter \\(N\\). Als er maar één groep is, zoals in de voorbeelden in dit hoofdstuk, dan geldt dat \\(n=N\\). In de beschrijvende statistiek wordt veel opgeteld, en daarvoor is dan ook een apart symbool, \\(\\sum\\), de griekse hoofdletter Sigma, waarmee een sommering of optelling wordt aangeduid. We zouden kunnen zeggen “tel alle geobserveerde waarden van variabele \\(X\\) bij elkaar op”, maar dat doen we doorgaans korter: \\[\\sum\\limits_{i=1}^n x_i, \\textrm{of korter} \\sum x % \\Sigma_i^N x_i, \\textrm{of nog korter} \\Sigma X\\] Op deze wijze wordt aangegeven dat alle scores \\(x_i\\) bij elkaar moeten worden opgeteld (gesommeerd), voor alle waarden van \\(i\\) (vanaf \\(i=1\\), tenzij anders aangegeven) tot \\(i=n\\). Alle \\(n\\) scores van de variabele \\(x\\) moeten dus worden opgeteld. Als er haakjes gebruikt worden, let dan goed op: handelingen beschreven binnen een paar haakjes hebben voorrang, die moet je dus eerst uitvoeren. Ook als dat niet strikt nodig is, zullen we vaak haakjes gebruiken ter verduidelijking, zoals in \\((2\\times3)+4=10\\). 9.3 Centrummaten 9.3.1 gemiddelde De meest bekende maat voor het centrum van een verdeling is wel het gemiddelde. Het gemiddelde is eenvoudig uit te rekenen door alle scores bij elkaar op te tellen, en vervolgens die som weer te delen door het aantal observaties. In symbolen: \\[\\begin{equation} \\overline{x} = \\frac{\\sum x}{n} = \\frac{1}{n} \\sum\\limits_{i}^n x_i \\tag{9.1} \\end{equation}\\] We maken hier meteen kennis met een nieuw symbool, \\(\\overline{x}\\), vaak “x-bar” genoemd, waarmee het gemiddelde van \\(x\\) wordt aangeduid. Het gemiddelde wordt ook vaak aangeduid met het symbool \\(M\\) (Eng. mean), o.a. in artikelen in de APA-stijl. Voorbeeld 9.1: In een winkel wordt bijgehouden hoe lang klanten moeten wachten bij de kassa, voordat ze aan de beurt zijn. Voor \\(N=10\\) klanten worden de volgende wachttijden geobserveerd, in minuten: 1, 2, 5, 2, 2, 2, 3, 1, 1, 3. De gemiddelde wachttijd is \\((\\sum X)/N = 22/10 = 2.2\\) minuten. Het gemiddelde van \\(X\\) wordt meestal uitgedrukt met één decimaalcijfer meer dan de scores van \\(X\\) (zie ook §9.6.1 hieronder over het aantal decimaalcijfers waarmee het gemiddelde wordt weergegeven). Het gemiddelde is op te vatten als het “balanspunt” van een verdeling: de observaties aan weerszijden houden elkaar “in evenwicht”, zoals geïllustreerd in Figuur 9.1, waar de “blokken” van het histogram precies “in evenwicht” zijn op het “balanspunt” bij het gemiddelde van 2.2. Het gemiddelde is ook de waarde ten opzichte waarvan de \\(N\\) observaties tezamen het minste verschillen, en het vormt dus een goed kenmerk voor het centrum van een kansverdeling. Het gemiddelde is alleen bruikbaar bij variabelen van het interval- of ratio-meetniveau. Figuur 9.1: Histogram van N=10 wachttijden, met markering van het gemiddelde. 9.3.2 mediaan De mediaan (symbool \\(Md\\) of \\(\\tilde{x}\\)) is de observatie in het midden van de rangorde van observaties 1. Als we de scores van \\(X\\) rangschikken van klein naar groot, dan is de mediaan het middenpunt van die gerangschikte reeks. De helft van de observaties is kleiner dan de mediaan, en de andere helft is groter dan de mediaan. Bij een oneven aantal observaties is de middelste observatie de mediaan. Bij een even aantal observaties wordt de mediaan gevormd door het gemiddelde van de middelste twee observaties. Voorbeeld 9.2: De wachttijden uit Voorbeeld 9.1 worden als volgt gerangschikt: 1, 1, 1, 2, 2, 2, 2, 3, 3, 5. De mediaan is het gemiddelde van de middelste twee (cursieve) observaties, dus 2 minuten. De mediaan is minder gevoelig dan het gemiddelde voor extreme waarden van \\(x\\). In het bovenstaande voorbeeld is de extreme wachttijd van 5 minuten van aanzienlijke invloed op het gemiddelde. Als we die waarde zouden verwijderen, dan verandert het gemiddelde van 2.2 naar 1.9, maar de mediaan blijft nog steeds 2. Extreme waarden hebben dus een minder grote invloed op de mediaan dan op het gemiddelde. Pas als de rangschikking van de observaties verandert, kan ook de mediaan veranderen. De mediaan is bruikbaar bij variabelen van het ordinale, interval- of ratio-meetniveau. 9.3.3 modus De modus (bijv.nw. ‘modaal’) is de waarde of score van \\(X\\) die het meeste voorkomt. Voorbeeld 9.3: Bij de wachttijden uit Voorbeeld 9.1 komt de score 2 het meeste voor (\\(4\\times\\)); dit is de modus. Voorbeeld 9.4: In 2022 was het gemiddelde besteedbare inkomen per huishouden in Nederland ca €35000. Het modale inkomen (per huishouden) was tussen €22000 en €24000 2: in deze inkomenscategorie vielen de meeste huishoudens. De modus is nog minder gevoelig dan de mediaan voor extreme waarden van \\(x\\). In het bovenstaande voorbeeld 9.2 maakt het niet uit wat de waarde van de langste wachttijd is: ook al zou die observatie de waarde \\(10\\) of \\(1000\\) hebben, de modus blijft onveranderlijk \\(2\\) (ga dat zelf na). De modus is bruikbaar bij variabelen van alle meetniveau’s. 9.3.4 harmonisch gemiddelde Als de afhankelijke variabele een breuk of verhouding of ratio voorstelt, zoals de snelheid waarmee een taak verricht wordt, dan vormt het (rekenkundig) gemiddelde van formule (9.1) eigenlijk niet een goede aanduiding voor de meest kenmerkende of centrale waarde. In dat geval kan je beter het harmonisch gemiddelde gebruiken: \\[\\begin{equation} H = \\frac{1}{\\frac{1}{n} \\sum\\limits_{i}^n \\frac{1}{x_i} } = \\frac{n}{\\sum\\limits_{i}^n \\frac{1}{x_i}} \\tag{9.2} \\end{equation}\\] Voorbeeld 9.5: Een student schrijft \\(n=3\\) teksten. Over de eerste tekst (500 woorden) doet hij 2.5 uur, over de tweede tekst (1000 woorden) doet hij 4 uur, en over de derde tekst (300 woorden) doet hij 0.6 uur. Wat is de gemiddelde schrijfsnelheid van deze student? De schrijfsnelheden zijn resp. 200, 250, en 500 woorden per uur, en het “normale” (rekenkundig) gemiddelde daarvan is 317 woorden per uur. Maar, het kostte 7.1 uur om 1800 woorden te schrijven, dus het “eigenlijke” gemiddelde is \\((500+1000+300)/(2.5+4+0.6)\\) \\(=1800/7.1=254\\) woorden per uur. De hoge schrijfsnelheid van de korte tekst telt voor \\(1/n\\) deel mee in het rekenkundig gemiddelde, hoewel die tekst slechts \\(300/1800=1/6\\) van het totaal aantal woorden bevat. Omdat de afhankelijke variabele een breuk is (snelheid, woorden/uur), vormt het harmonische gemiddelde hier een betere centrummaat. Eerst rekenen we de snelheid (woorden per tijdseenheid) om naar de inverse daarvan (zie (9.2), in noemer, binnen somteken), d.w.z. naar tijd per woord: 0.005, 0.004, en 0.002 (tijdseenheden per woord, zie voetnoot3). Ten tweede middelen we deze tijden, tot gemiddeld 0.00366 uur per woord (13.2 seconden per woord), en tenslotte nemen we daarvan wederom de inverse. De harmonisch gemiddelde schrijfsnelheid is dan \\(1/0.00366=273\\) woorden per uur, dichter bij het “eigenlijke” gemiddelde van 254 woorden per uur. 9.3.5 gewinsoriseerde gemiddelde De grote gevoeligheid van het gewone (rekenkundige) gemiddelde voor uitbijters kan ingeperkt worden, door de meest extreme observaties te wijzigen naar de minder extreme, meer centrale observaties. Het gemiddelde van deze (deels gewijzigde) observaties wordt het gewinsoriseerde gemiddelde genoemd (“winsorized mean”; ). Voorbeeld 9.6: De wachttijden uit Voorbeeld 9.1 worden als volgt gerangschikt: 1, 1, 1, 2, 2, 2, 2, 3, 3, 5. Voor het 10% gewinsoriseerde gemiddelde worden (na rangordening) de 10% kleinste observaties gelijk gesteld aan de eerstvolgende grotere waarde, en de 10% grootste observaties worden gelijk gesteld aan de laatst voorafgaande kleinere waarde (gewijzigde waarden zijn hier cursief): 1, 1, 1, 2, 2, 2, 2, 3, 3, 3. Het winsorized gemiddelde over deze gewijzigde waarden is \\(\\overline{x}_w=2\\) minuten. 9.3.6 getrimde gemiddelde Een nog drastischer ingreep is om de meest extreme observaties geheel te verwijderen. Het gemiddelde van de overblijvende observaties wordt het getrimde gemiddelde genoemd (“trimmed mean”). Bij een 10% trim verwijderen we de onderste 10% en de bovenste 10% van de observaties; er resteren dan dus nog slechts \\((1- (2 \\times (10/100))\\times n\\) observaties (Wilcox 2012). Voorbeeld 9.7: De wachttijden uit Voorbeeld 9.1 worden weer als volgt gerangschikt: 1, 1, 1, 2, 2, 2, 2, 3, 3, 5. Voor het 10% getrimde gemiddelde worden (na rangordening) de 10% kleinste observaties verwijderd, en de 10% grootste observaties worden eveneens verwijderd: 1, 1, 2, 2, 2, 2, 3, 3. Het getrimde gemiddelde over deze \\(10-(.2)(10)=8\\) resterende waarden is hier \\(\\overline{x}_t=2\\) minuten. 9.3.7 vergelijking van centrummaten Figuur 9.2 illustreert de verschillen tussen de diverse centrummaten, voor asymmetrisch verdeelde observaties. Figuur 9.2: Histogram van een variabele met positief scheve (asymmetrische) frequentieverdeling, met daarin aangegeven (1) de mediaan, (2) het 10% getrimde gemiddelde, (3) het 10% gewinsoriseerde gemiddelde, en (4) het rekenkundig gemiddelde. De geobserveerde scores zijn gemarkeerd langs de horizontale as. Het rekenkundig gemiddelde is het meest gevoelig voor extreme waarden: de extreme waarden “trekken” erg hard aan het gemiddelde. Deze invloed van extreme waarden wordt getemperd in het gewinsoriseerde gemiddelde, en wordt nog meer getemperd in het getrimde gemiddelde. Naarmate de trimfactor (het percentage van de observaties dat wordt gewijzigd of verwijderd) toeneemt, gaan de gewinsoriseerde en getrimde gemiddelden meer lijken op de mediaan. Immers, bij een trimfactor van 50% resteert er van alle observaties nog maar één (ongewijzigde) observatie, en dat is de mediaan (ga dat zelf na). In §9.7 gaan we verder in op de keuze voor een passende maat voor het centrum van een verdeling. 9.4 Kwartielen en boxplots De verdeling van een variabele wordt niet alleen gekenmerkt door het centrum van die verdeling, maar ook door de mate van spreiding rondom het centrum, d.w.z. hoe groot het verschil is van de observaties t.o.v. het gemiddelde. We willen bijvoorbeeld niet alleen weten wat het gemiddelde inkomen is, maar ook hoe groot de verschillen in inkomen zijn. 9.4.1 Kwartielen Een eenvoudige en bruikbare maat daarvoor zijn de kwartielen (Tukey 1977). We delen de gerangschikte observaties in twee helften op; de grens daartussen is de mediaan. Vervolgens halveren we weer elke helft, tot kwarten. De kwartielen worden gevormd door de grenzen tussen deze kwarten; er zijn dus drie kwartielen. Het eerste kwartiel \\(Q_1\\) is de mediaan van de onderste helft, \\(Q_2\\) is de mediaan van alle \\(n\\) observaties, en het derde kwartiel \\(Q_3\\) is de mediaan van de bovenste helft. De helft van de observaties (nl. het tweede en derde kwart) ligt tussen \\(Q_1\\) en \\(Q_3\\). De afstand tussen \\(Q_1\\) en \\(Q_3\\) wordt de “interquartile range” genoemd (IQR). Deze IQR vormt een eerste bruikbare maat voor de spreiding van observaties ten opzichte van hun centrale waarde. Voor de uitleg maken we gebruik van de fictieve scores op een leestoets, gegeven in Tabel 9.1. Tabel 9.1: Scores van N=10 leerlingen op drie onderdelen van de CITO-toets, afgenomen in groep 8 van het basisonderwijs. Leerling Lezen Rekenen Wereldoriëntatie 1 18 22 55 2 32 36 55 3 45 34 38 4 25 25 40 5 27 29 48 6 23 20 44 7 29 27 49 8 26 25 42 9 20 25 57 10 25 27 47 \\(\\sum x\\) 270 270 475 \\(\\overline{x}\\) 27.0 27.0 47.5 Voorbeeld 9.8: De scores bij het onderdeel Lezen in Tabel 9.1 zijn als volgt gerangschikt: 18, 20, 23, 25, 25, 26, 27, 29, 32, 45. De mediaan is \\(Q_2=25.5\\) (tussen de 5e en 6e observatie in deze ranglijst). De mediaan van de onderste helft is \\(Q_1=23\\) en die van de bovenste helft is \\(Q_3=29\\). De interquartile range is \\(\\textrm{IQR}=29-23=6\\). 9.4.2 Uitbijters In de leesscores in Tabel 9.1 treffen we een extreme waarde aan, nl. de score 45, die opvallend veel verschilt van het gemiddelde. Zo’n opvallende waarde wordt aangeduid als een “uitbijter” (in het Engels als “outlier”). De grens voor wat we beschouwen als een uitbijter ligt doorgaans bij \\(1.5 \\times \\textrm{IQR}\\). Als een waarde meer dan \\(1.5 \\times \\textrm{IQR}\\) boven \\(Q_3\\) of onder \\(Q_1\\) ligt, dan beschouwen we die observatie als een uitbijter. Controleer deze observaties nog eens (denk aan het principe van zorgvuldigheid, zie §3.1). Voorbeeld 9.9: Voor de eerder genoemde leesscores in Tabel 9.1 vonden we \\(Q_1=23\\), \\(Q_3=29\\), en \\(\\textrm{IQR}=Q_3-Q_1=29-23=6\\). De bovenste grenswaarde voor uitbijters is \\(Q_3 + 1.5 \\times \\textrm{IQR} = 29 + 1.5 \\times 6 = 29+9 = 38\\). De observatie met score 45 ligt boven deze grenswaarde, en wordt daarom beschouwd als uitbijter. 9.4.3 Boxplots De frequentieverdeling van een variabele kunnen we nu weergeven met vijf kenmerken, de zgn. “five-number summary”, nl. de kleinste waarde, \\(Q_1\\), mediaan, \\(Q_3\\), en grootste waarde. Deze vijf kenmerken worden grafisch weergegeven in een zgn. “boxplot”, zie Figuur 9.3 voor een voorbeeld (Tukey 1977, sec. 2C). ## Loading required package: foreign Figuur 9.3: Boxplots van scores van \\(N=10\\) leerlingen op de onderdelen Lezen en Rekenen van de CITO-toets (zie Tabel 9.1), met uitbijters als open cirkels gemarkeerd. De geobserveerde scores zijn gemarkeerd langs de verticale assen. De box omspant het gebied van (bij benadering) \\(Q_1\\) tot \\(Q_3\\), en omspant dus de centrale helft van de observaties. De dikkere lijn in de box markeert de mediaan. De lijnen strekken zich uit naar de kleinste en grootste waarden die géén uitbijters zijn 4. De afzonderlijke uitbijters worden hier met een apart symbool aangeduid. 9.5 Spreidingsmaten 9.5.1 variantie Een andere manier om de spreiding van de observaties aan te geven, zou zijn om te kijken naar de afwijking van iedere observatie ten opzichte van het gemiddelde, dus \\((x_i-\\overline{x})\\). Maar als we al die afwijkingen optellen, dan is de som daarvan altijd nul! De positieve en negatieve afwijkingen heffen elkaar immers op (ga dat zelf na in Tabel 9.1). Daarom middelen we niet de afwijkingen zelf, maar de kwadraten van die afwijkingen. Zowel de positieve als de negatieve afwijkingen resulteren in positieve gekwadrateerde-afwijkingen. Van al die gekwadrateerde-afwijkingen berekenen we het gemiddelde, d.w.z. we tellen ze op en delen door \\((n-1)\\), zie voetnoot5. Het resultaat noemen we de variantie, aangeduid met symbool \\(s^2\\): \\[\\begin{equation} s^2 = \\frac{ \\sum (x_i - \\overline{x})^2 } {n-1} \\tag{9.3} \\end{equation}\\] De teller van deze breuk wordt wel aangeduid als de “sum of squared deviations” of “sum of squares” (SS) en de noemer wordt wel aangeduid als het aantal vrijheidsgraden of “degrees of freedom” van de teller (d.f.; zie §12.12.1). De variantie rekenen we tegenwoordig altijd uit met een rekenmachine of computer. 9.5.2 standaarddeviatie Om de bovenstaande variantie te berekenen, hebben we de afwijkingen van de observaties gekwadrateerd. De variantie is dus een grootheid die niet wordt uitgedrukt in de oorspronkelijke eenheden (bijv. seconden, cm, score), maar in gekwadrateerde eenheden (bijv. \\(\\textrm{s}^2\\), \\(\\textrm{cm}^2\\), \\(\\textrm{score}^2\\)). Teneinde weer terug te keren naar de oorspronkelijke eenheden, nemen we de wortel uit de variantie. Het resultaat noemen we de standaarddeviatie, aangeduid met symbool \\(s\\): \\[\\begin{equation} s = \\sqrt{s^2} = \\sqrt{ \\frac{ \\sum (x_i - \\overline{x})^2 } {n-1} } \\tag{9.4} \\end{equation}\\] Voorbeeld 9.10: Het gemiddelde van de eerder genoemde leesscores in Tabel 9.1 is \\(27.0\\), en de afwijkingen zijn als volgt: -9, 5, 18, -2, 0, -4, 2, -1, -7, -2. De gekwadrateerde afwijkingen zijn 81, 25, 324, 4, 0, 16, 4, 1, 49, 4. De som van deze gekwadrateerde afwijkingen is 508, en de variantie is dan \\(s^2=508/9=56.44\\). De standaarddeviatie is de wortel van de variantie, dus \\(s=\\sqrt{508/9}=7.5\\). De variantie en standaarddeviatie zijn alleen bruikbaar bij variabelen van het interval- of ratio-meetniveau. Ook de variantie en standaarddeviatie kunnen weer gebaseerd zijn op de gewinsoriseerde of getrimde verzameling van observaties. De standaarddeviatie hebben we nodig (a) als we de ruwe observaties willen omzetten naar standaardscores (zie §9.8 hieronder), (b) als we een variabele willen beschrijven die normaal verdeeld is (zie §10.3, en (c) als we hypotheses willen toetsen met behulp van een normaal verdeelde variabele (zie §12.12 e.v.). 9.5.3 MAD Behalve de standaarddeviatie is er ook een robuuste tegenhanger daarvan, die niet gebruik maakt van het gemiddelde. Deze maat is daarom minder gevoelig voor uitbijters (robuuster), wat soms handig is. We kijken hiervoor naar de afwijking van iedere observatie ten opzichte van de mediaan (niet t.o.v. het gemiddelde). Van deze afwijkingen nemen we de absolute waarde6 (niet het kwadraat). Van deze absolute afwijkingen bepalen we tenslotte weer de mediaan (niet het gemiddelde). Het resultaat noemen we de “median absolute deviation” (MAD): \\[\\begin{equation} \\textrm{MAD} = k ~~ Md ( |x_i - Md(x) |) \\tag{9.5} \\end{equation}\\] Hierbij gebruiken we meestal \\(k=1.4826\\) als constante; door deze schaalfactor komt de MAD ruwweg overeen met de standaarddeviatie \\(s\\) indien \\(x\\) normaal verdeeld zou zijn (§10.3). Voorbeeld 9.11: De mediaan van de eerder genoemde leesscores in Tabel 9.1 is 25.5, en de afwijkingen van die mediaan zijn als volgt: -7.5, 6.5, 19.5, -0.5, 1.5, -2.5, 3.5, 0.5, -5.5, -0.5. De gerangschikte, absolute afwijkingen zijn 0.5, 0.5, 0.5, 1.5, 2.5, 3.5, 5.5, 6.5, 7.5, 19.5. De mediaan van deze 10 absolute afwijkingen is 3, en \\(\\textrm{MAD} = 1.4826 \\times 3 = 4.4478\\). Merk op dat de MAD kleiner is dan de standaarddeviatie, o.a. omdat de MAD minder gevoelig is voor de extreme waarde \\(x_3=45\\). 9.6 Over significante cijfers 9.6.1 Gemiddelde en standaarddeviatie Een gemiddelde uitkomst wordt weergegeven in een beperkt aantal significante cijfers, d.i. een beperkt aantal cijfers, van links naar rechts geteld, zonder acht te slaan op het decimaalteken. Het aantal significante cijfers van een gemiddelde uitkomst moet gelijk zijn aan het aantal significante cijfers van het aantal observaties waarover is gemiddeld. (Overige cijfers in de gemiddelde uitkomst zijn niet nauwkeurig bepaald.) De gemiddelde uitkomst moet eerst afgerond worden tot het gepaste aantal significante cijfers, voordat de uitkomst verder geïnterpreteerd wordt, zie Tabel 9.2. Tabel 9.2: Het aantal significante cijfers in het gerapporteerde gemiddelde is gelijk aan het aantal significante cijfers van het aantal observaties. aant.obs. aant.signif.cijfers voorbeeld gemiddelde rapporteer als \\(1\\dots9\\) 1 21/8 = 2.625 3 \\(10\\dots99\\) 2 57/21 = 2.714286 2.7 \\(100\\dots999\\) 3 317/120 = 2.641667 2.64 \\(1000\\dots9999\\) 4 3179/1234 = 2.576175 2.576 Het aantal significante cijfers in de gerapporteerde standaarddeviatie is hetzelfde als in het gemiddelde, volgens Tabel 9.2. 9.6.1.1 Achtergrond Laten we aannemen dat ik de afstand van mijn huis naar mijn werk langs een vaste route een aantal keren heb gemeten. Het gemiddelde van die metingen bedraagt zogenaamd \\(2.954321\\) km. Door het gemiddelde te rapporteren met 7 cijfers suggereer ik hier dat ik precies weet dat de afstand \\(2954321\\) millimeter is, en ten hoogste \\(1\\) mm meer of minder: het laatste cijfer is geschat of afgerond. Het aantal significante cijfers (in dit voorbeeld 7) geeft de mate van nauwkeurigheid aan. In dit voorbeeld is de gesuggereerde nauwkeurigheid van 1 mm duidelijk onjuist, o.a. omdat beginpunt en eindpunt niet tot op de millimeter bepaald zijn. Het is daarom gebruikelijk om de gemiddelde gemeten afstand te rapporteren met een aantal significante cijfers dat de nauwkeurigheid van die metingen en van het gemiddelde aangeeft, bijv. \\(3.0\\) km (per auto of fiets) of \\(2.95\\) km (te voet). Dezelfde gedachtengang is van toepassing bij de meting van een kenmerk door middel van een enquête-vraag. Met \\(n=15\\) respondenten zou de gemiddelde score \\(43/15 \\approx 2.86667\\) kunnen zijn. Maar de nauwkeurigheid is in dit voorbeeld niet zo goed als dit decimale getal suggereert. In feite zorgt hier één afwijkend antwoord al voor een afwijking van \\(\\pm0.06667\\) in het gemiddelde. Bovendien is een gemiddelde score altijd het resultaat van een deling, en “(bij) delen en vermenigvuldigen geldt de regel dat de uitkomst evenveel significante cijfers bevat als de meetwaarde met het kleinste aantal significante cijfers.”7 In dit voorbeeld bestaan de teller (\\(43\\)) en de noemer (\\(15\\)) van het gemiddelde beide uit 2 significante cijfers, en dient de uitkomst dus ook uit 2 significante cijfers te bestaan. De gemiddelde score dient gerapporteerd te worden als \\(2.9\\) punten, met slechts één cijfer achter het decimaalteken. 9.6.2 Percentages Een percentage is een verhouding of breuk, vermenigvuldigd met \\(100\\). Gebruik en rapporteer een afgerond percentage (d.i. twee significante cijfers) alleen indien de noemer van de verhouding of breuk groter is dan 100. Deze noemer geeft het aantal waarnemingen of gevallen. Als de noemer kleiner is dan 100 (waarnemingen, gevallen), dan zijn percentages misleidend, zie Tabel 9.3. Tabel 9.3: Het aantal significante cijfers in de gerapporteerde proportie (of percentage) hangt samen met het aantal significante cijfers van het aantal observaties in de noemer van de breuk. aant.obs.(noemer) aant.signif.cijfers voorbeeld breuk rapporteer als \\(1\\dots9\\) 1 3/8 = 0.4 3/8 \\(10\\dots99\\) 2 21/57 = 0.36 21/57 \\(100\\dots999\\) 3 120/317 = 0.378 38% \\(1000\\dots9999\\) 4 1234/3179 = 0.3882 38.8% 9.6.2.1 Achtergrond De regels voor percentages vloeien voort uit die in §9.6.1 toegepast op delingen. Als de noemer groter is dan 100 is het percentage (met twee significante cijfers) het gevolg van een schaalverandering “naar beneden” (van een noemer groter dan 100 naar een noemer van precies 100 percentagepunten). De percentageschaal is minder nauwkeurig dan de oorspronkelijke verhouding; de percentages zijn afgerond tot op twee significante cijfers; het laatste significante cijfer van het percentage is dus geborgd. Als de noemer echter kleiner is dan 100 dan is het percentage (met twee significante cijfers) het gevolg van een “oprekking naar boven” (van een noemer kleiner dan 100 naar een noemer van precies 100 percentagepunten). De percentageschaal suggereert dan een pseudo-nauwkeurigheid die er niet was in de oorspronkelijke verhouding, en de nauwkeurigheid van de percentageschaal is vals. Als de noemer kleiner is dan 100, zijn percentages dus misleidend. Voorbeeld 9.12: In een cursus van 29 studenten zijn er 23 studenten geslaagd. We spreken dan vaak van een cursusrendement van \\(23/29=\\) 79%. Toch is zo’n weergave als percentage in dit geval misleidend. Laten we daarvoor eens kijken naar de 6 gezakten. Je kunt beredeneren dat het aantal van 6 gezakten een eigen afrondingsfout heeft van \\(1/2\\) student; bij omzetting naar de percentageschaal wordt ook deze afrondingsfout mee vergroot, zodat de percentages minder nauwkeurig zijn dan de hele percentages (2 significante cijfers) suggereren. Of anders gezegd: het aantal van 6 gezakten (d.i. een getal met één significant cijfer) noopt ons om ook de verhouding weer te geven met slechts één significant cijfer, en dus niet als percentage. Rapporteer bij voorkeur de verhouding zelf (\\(23/29\\)), of eventueel de “odds” (\\(23/6=4\\)) afgerond tot het juiste aantal significante cijfers8. Op grond van dezelfde overwegingen is een percentage met een decimaal cijfer (d.i. met drie significante cijfers, bv. “36.1%”) alleen zinnig als de noemer van de verhouding of breuk groter is dan 1000. Voorbeeld 9.13: In 2013 startten 154 studenten met een tweejarige research master. Na 2 jaar waren 69 daarvan afgestudeerd. Het nominaal rendement voor dit cohort is dan \\(69/154=\\) 0.448052, af te ronden en te rapporteren als 45% (niet als 44.81%). 9.7 Keuzemoment De verdeling van een variabele kun je op verschillende manieren beschrijven. Als variabele \\(X\\) gemeten is op het interval- of ratio-meetniveau, begin dan altijd met een histogram (§8.4) en een boxplot (§9.4.3). De centrummaten en spreidingsmaten zijn te ordenen zoals in Tabel 9.4. Tabel 9.4: Overzicht van besproken centrummaten en spreidingsmaten. Verdeling Centrummaat Spreidingsmaat alle mediaan kwartielen, IQR, MAD … getrimde of gewins.gemidd. getrimde of gewins.std.dev. (a & b & c) gemiddelde standaarddeviatie De meest robuuste maten staan bovenin (mediaan, kwartielen, IQR, MAD). Deze maten zijn robuust: ze zijn weinig gevoelig voor uitbijters of voor eventuele asymmetrie in de frequentieverdeling, zoals de voorbeelden in dit hoofdstuk laten zien. De meest efficiënte maten staan onderin in Tabel 9.4: gemiddelde en standaarddeviatie. Deze maten zijn efficiënt: ze geven het centrum en de spreiding het beste weer, ze hebben zelf de kleinste standaarddeviatie, en ze hebben daarvoor relatief het kleinste aantal observaties nodig. De andere maten nemen een tussenpositie in: de getrimde maten zijn wat robuuster, en de gewinsoriseerde maten wat efficiënter. De meest efficiënte maten vereisen echter ook de meest vèrgaande assumpties (en de meest robuuste maten vereisen de minste assumpties). Deze efficiënte maten zijn alleen zinnig, indien de verdeling van \\(X\\) voldoet aan drie assumpties: (a) de verdeling is min of meer symmetrisch, d.w.z. de linker- en rechter-helft van het histogram resp. de bovenste en onderste helft van de boxplot lijken elkaars spiegelbeeld, (b) de verdeling is unimodaal, d.w.z. de verdeling heeft één modus, en (c) de verdeling bevat geen of nauwelijks uitbijters. Inspecteer deze assumpties in het histogram en de boxplot van \\(X\\). Als aan één van deze assumpties niet is voldaan, dan doe je er beter aan om meer robuuste maten te gebruiken om de verdeling te beschrijven. 9.8 Standaardscores Soms kan het handig zijn om scores te vergelijken die gemeten zijn op verschillende schalen. Bijvoorbeeld: Jan had een 8 als eindcijfer voor wiskunde op het VWO, en zijn IQ is 136. Is de afwijking van Jan ten opzichte van het gemiddelde even groot op beide schalen? Om zo’n vraag te beantwoorden moeten we de scores van de twee variabelen uitdrukken op dezelfde meetschaal. Dat doen we door de ruwe scores om te rekenen naar standaard-scores, of z-scores. Hiervoor nemen we de afwijking van iedere score ten opzichte van het gemiddelde, en we delen die afwijking door de standaarddeviatie: \\[\\begin{equation} z_i = \\frac{(x_i-\\overline{x})}{s_x} \\tag{9.6} \\end{equation}\\] De standaardscore of z-score representeert dus de afstand van de \\(i\\)’de observatie tot het gemiddelde van \\(x\\), uitgedrukt in eenheden standaarddeviatie. Bij een standaardscore van \\(z=-1\\) is de geobserveerde score precies \\(1 \\times s\\) beneden het gemiddelde \\(\\overline{x}\\). Bij een standaardscore van \\(z=+2\\) dan is de geobserveerde score precies \\(2 \\times s\\) boven het gemiddelde9. De z-scores zijn ook handig om twee variabelen te vergelijken die weliswaar op dezelfde schaal gemeten zijn (bijvoorbeeld een schaal van \\(1 \\dots 100\\)), maar die toch verschillende gemiddelden en/of verschillende standaarddeviaties hebben, zoals de scores in Tabel 9.1. In Hoofdstuk 10 zullen we verder werken met z-scores. De standaardscore of z-score heeft twee handige eigenschappen die je moet onthouden. Ten eerste is het gemiddelde altijd gelijk aan nul: \\(\\overline{z}=0\\), en ten tweede is de standaarddeviatie gelijk aan 1: \\(s_z = 1\\). (Deze eigenschappen volgen uit de definitie in formule (9.6); het wiskundige bewijs laten we hier achterwege.) Dus de transformatie van een verzameling observaties naar standaardscores of z-scores levert altijd een verdeling op met een gemiddelde van nul en een standaarddeviatie van één. Bedenk wel dat deze transformatie naar standaardscores alleen zinnig is, indien en voor zover het gemiddelde en de standaarddeviatie ook zinnige maten zijn om de verdeling van \\(x\\) te beschrijven (zie §9.7). 9.9 SPSS Voor histogram, percentielen en boxplot: Analyze > Descriptive Statistics > Explore... Selecteer variabele (sleep naar Variable(s) paneel) Kies Plots, vink aan: Histogram, en bevestig met Continue Kies Options, vink aan: Percentiles, en bevestig met Continue en daarna met OK. De uitvoer bevat zowel beschrijvende statistiek als histogram en boxplot. Voor kenmerkende getallen: Analyze > Descriptive Statistics > Descriptives... Selecteer variabele (sleep naar Variable(s) paneel) Kies Options; vink aan: Mean, Sum, Std.deviation, Variance, Minimum, Maximum, en bevestig met Continue en daarna met OK. De uitvoer bevat de gevraagde statistische kenmerken van de verdeling van de variabele. Voor mediaan: Analyze > Compare Means > Means... Selecteer variabele (sleep naar Variable(s) paneel) Kies Options; vink aan: Mean, Number of cases, Standard deviation, Variance, Minimum, Maximum en ook Median, en bevestig met Continue en daarna met OK. De uitvoer bevat de gevraagde statistische kenmerken van de verdeling van de variabele. Standaardscores uitrekenen en bewaren in nieuwe kolom: Analyze > Descriptive Statistics > Descriptives... Selecteer variabele (sleep naar Variable(s) paneel) Vink aan: Save standardized values as variables en bevestig met OK. De nieuwe variabele(n) met z-scores worden toegevoegd als nieuwe kolom(men) aan het databestand. 9.10 JASP Voor histogram en boxplot: Klik in de bovenbalk op Descriptives. Selecteer vervolgens de gewenste variabele(n) en plaats ze in het veld “Variables”. Klik de balk Plots open en vink Distribution plots aan onder “Basic plots” voor een histogram (of staafdiagram, afhankelijk van het meetniveau). Vink Boxplots (met Boxplot element) aan onder “Customizable plots” voor een boxplot. Voor kenmerkende getallen en kwartielen/percentielen: Klik in de bovenbalk op Descriptives. Selecteer vervolgens de gewenste variabele(n) en plaats ze in het veld “Variables”. Klik de balk Statistics open en vink Quartiles aan onder “Percentile Values” om kwartielen te krijgen. Door Cut points for: aan te vinken en een getal in te vullen, kun je ook andere fractielen krijgen (bijvoorbeeld sextielen als je 6 invult). Door Percentiles: aan te vinken en een getal in te vullen, kun je ook een specifiek percentiel krijgen (bijvoorbeeld 44e percentiel als je 44 invult). Voor kenmerkende getallen, zoals centrummaten en spreidingsmaten, vink je Mean, Median, Mode en Sum aan onder “Central Tendency”, en Variance, Std.deviation, MAD Robust (bij MAD is de constante op 1 gezet), IQR, Minimum en Maximum aan onder “Dispersion”. De uitvoer is een tabel met alle opgevraagde beschrijvende statistiek van de variabele(n). Standaardscores moeten zelf worden uitgerekend als nieuwe variabele. Maak de nieuwe variabele aan door op de +-button te klikken rechts van de naam van de laatste kolom in de dataset. Er verschijnt een paneel “Create Computed Column”, waar je een naam voor de nieuwe variabele kunt invullen, bijvoorbeeld zLezen. Ook kun je kiezen uit R en een aanwijshandje. Dit zijn de twee opties in JASP om formules te definiëren waarmee de nieuwe (lege) variabele wordt gevuld; met R code of met drag-and-drop in JASP. Hieronder wordt voor allebei de opties uitgelegd hoe we hiermee de variabele met standaardscores maken. Als laatste kun je aanvinken welk meetniveau de variabele moet krijgen; dat is in dit geval “Scale”, net als de originele variabele. Klik vervolgens op Create Column om de nieuwe variabele aan te maken. De nieuwe variabele (kolom) verschijnt als meest rechtse in de data en is nog leeg. Als computed with R code is gekozen als optie om de nieuwe variabele te definiëren, dan verschijnt er boven de data een veld met de tekst “Enter your R code here”. Hier kan R code worden ingevoerd die met behulp van R functies de nieuwe variabele definieert. Om een standaardscore te berekenen kan de volgende code worden gebruikt (hier voor de variabele ‘Lezen’): ((Lezen - mean(Lezen)) / sd(Lezen)) Vul de R code in, vul de overige velden ook in, en klik onderaan op Compute column om de lege variabele te vullen met de standaardscores. Rekenkundige bewerkingen op de variabelen zijn in JASP alleen mogelijk indien het meetniveau van alle gebruikte variabelen is ingesteld op ‘Scale’. Als drag-and-drop (aanwijshandje) is gekozen als optie om de nieuwe variabele te definiëren, dan verschijnt er boven de data een veld met links daarvan de variabelen, boven wiskundige symbolen, en rechts een aantal functies. Hier kan handmatig de functie worden geselecteerd die de nieuwe variabele definieert. Om standaardscores te berekenen sleep je de variabele die je wilt omzetten in het lege veld. Klik dan bovenaan op \\(-\\) en selecteer daarna uit het rechter menu de functie mean(y). Sleep de om te zette variabele naar “values” en klik daarna bovenaan op \\(\\div\\). Selecteer in het rechter menu de functie σy voor op de puntjes onder het breukteken (als noemer) en sleep ook hier de om te zette variabele naar “values”. Voor de variabele ‘Lezen’ zou dit er zo uitzien in het werkblad: \\[\\frac{(Lezen-mean(Lezen)}{\\sigma Lezen}\\] Vul ook de overige velden in, en klik dan onderaan op Compute column om de lege variabele te vullen met de standaardscores. Rekenkundige bewerkingen op de variabelen zijn in JASP alleen mogelijk indien het meetniveau van alle gebruikte variabelen is ingesteld op ‘Scale’. Als je een fout hebt gemaakt of als je de code of formule voor de nieuwe variabele wilt aanpassen, dan kan je altijd terugkeren naar dit ‘Computed Column’ veld door te klikken op het formule-symbool \\(f_x\\) naast de variabele-naam, of door te klikken op de variabele-naam zelf. 9.11 R Voor kwartielen en boxplot zoals Figuur 9.3 gebruiken we de commando’s fivenum, quantile, en boxplot: require(foreign) # for foreign::read.spss cito <- read.spss("data/cito.sav") fivenum(cito$Lezen) # minimum, Q1, mediaan, Q3, maximum ## [1] 19 22 26 29 44 quantile(cito$Lezen, c( 1/4, 3/4 ) ) # Q1 en Q3, anders berekend ## 25% 75% ## 22.75 28.75 op <- par(mar=c(4,4,1,2)+0.1) # smaller margins with(cito, boxplot(Lezen, Rekenen, col="grey80", lwd=2, lty=1, ylab="Score", ylim=c(17,45) ) ) axis(side=1, at=c(1,2), labels=c("Lezen","Rekenen") ) plotrix::axis.break(axis=2) # break in linker Y-as rug(cito$Lezen, side=2) # markeringen linker Y-as rug(cito$Rekenen, side=4) # markeringen rechter Y-as Veel centrummaten zijn als functie in R voorgeprogrammeerd: mean(cito$Lezen) # gemiddelde ## [1] 27.2 psych::winsor.mean(cito$Lezen, trim=.1) # gewinsoriseerde gemiddelde, uit psych package ## [1] 26.3 mean(cito$Lezen, trim=.1) # getrimde gemiddelde ## [1] 26.125 median(cito$Lezen) # mediaan ## [1] 26 Ook diverse spreidingsmaten zijn voorgeprogrammeerd: var(cito$Lezen) # variantie ## [1] 50.17778 sd(cito$Lezen) # standaarddeviatie, sd(x) = sqrt(var(x)) ## [1] 7.083627 mad(cito$Lezen) # MAD ## [1] 5.1891 Daarentegen moeten we standaardscores zelf uitrekenen, en zelf bewaren als nieuwe variabele, hier zLezen genoemd (let op de haakjes in de eerste regel): zLezen <- (cito$Lezen-mean(cito$Lezen)) / sd(cito$Lezen) # z-scores head(zLezen) # eerste paar observaties van variabele zLezen ## [1] -1.1575990 0.6776189 2.3716662 -0.3105753 0.1129365 -0.7340872 Referenties Tukey, John W. 1977. Exploratory Data Analysis. Reading, MA: Addison-Wesley. Wilcox, Rand R. 2012. Introduction to Robust Estimation and Hypothesis Testing. 3rd ed. Burlington: Elsevier. In het Engels wordt de streep of middenberm in het midden van een weg aangeduid als de “median” (Am.Eng.) of “central reservation” (Br.Eng.); deze streep verdeelt de weg in twee meestal even grote helften.↩︎ https://www.cbs.nl/nl-nl/visualisaties/inkomensverdeling↩︎ Dit is vergelijkbaar met sporten als roeien, zwemmen, wielrennen, schaatsen, e.d., waar ook de tijd over een afgesproken afstand gemeten en vergeleken wordt, en niet de snelheid over een afgesproken tijd.↩︎ In een klassieke boxplot strekken de lijnen zich uit naar het minimum en maximum (Tukey 1977) en worden uitbijters niet apart aangeduid.↩︎ We delen door \\(n-1\\) en niet door \\(n\\), om een betere schatting te krijgen van de spreiding in de populatie. Op deze manier houden we rekening met het feit dat we één kenmerk van de steekproef (nl. het gemiddelde) gebruiken om de spreiding te bepalen. Als je alleen geïnteresseerd bent in de spreiding in je steekproef van observaties, en niet in de populatie, deel dan door \\(n\\).↩︎ Positieve afwijkingen blijven ongewijzigd, van negatieve afwijkingen wordt het teken omgekeerd.↩︎ https://nl.wikipedia.org/wiki/Significant_cijfer↩︎ Deze “odds” geeft aan dat er 23 geslaagden zijn op 6 gezakten, d.i. afgerond 4 geslaagden voor iedere gezakte.↩︎ Check: \\(z = +2 = \\frac{(x_i-\\overline{x})}{s_x}\\), dus \\(2 s = (x_i-\\overline{x})\\), dus \\(x_i = \\overline{x}+2s\\).↩︎ "],["404.html", "Page not found", " Page not found The page you requested cannot be found (perhaps it was moved or renamed). You may want to try searching to find the page's new location, or use the table of contents to find the page you are looking for. "]] +[["app-kritieketwaarden.html", "Appendix C Kritieke waarden van \\(t\\)-verdeling", " Appendix C Kritieke waarden van \\(t\\)-verdeling De hieronder gegeven kritieke grenswaarde \\(t^*\\) heeft een overschrijdingskans \\(p\\) onder \\(H_0\\), d.w.z. \\(P(t \\geq t^*|H_0)=p\\), en heeft kans \\(B\\) op een waarde tussen \\((-t^*, +t^*)\\). De \\(t\\)-verdeling is symmetrisch rond \\(t=0\\), dus \\(P(t < -t^*) = P(t > t^*)\\). De tabel hieronder geeft de kritieke grenswaarden \\(t^*\\) voor veel gebruikte overschrijdingskansen \\(p\\) en betrouwbaarheidsintervallen \\(B\\), voor de vrijheidsgraden aangegeven in de eerste kolom. p 0.2 0.1 0.05 0.025 0.01 0.005 0.0025 0.001 B 60% 80% 90% 95% 98% 99% 99.5% 99.8% 1 1.376 3.078 6.314 12.706 31.821 63.657 127.321 318.309 2 1.061 1.886 2.920 4.303 6.965 9.925 14.089 22.327 3 0.9785 1.638 2.353 3.182 4.541 5.841 7.453 10.215 4 0.941 1.533 2.132 2.776 3.747 4.604 5.598 7.173 5 0.9195 1.476 2.015 2.571 3.365 4.032 4.773 5.893 6 0.9057 1.440 1.943 2.447 3.143 3.707 4.317 5.208 7 0.896 1.415 1.895 2.365 2.998 3.499 4.029 4.785 8 0.8889 1.397 1.860 2.306 2.896 3.355 3.833 4.501 9 0.8834 1.383 1.833 2.262 2.821 3.250 3.690 4.297 10 0.8791 1.372 1.812 2.228 2.764 3.169 3.581 4.144 11 0.8755 1.363 1.796 2.201 2.718 3.106 3.497 4.025 12 0.8726 1.356 1.782 2.179 2.681 3.055 3.428 3.930 13 0.8702 1.350 1.771 2.160 2.650 3.012 3.372 3.852 14 0.8681 1.345 1.761 2.145 2.624 2.977 3.326 3.787 15 0.8662 1.341 1.753 2.131 2.602 2.947 3.286 3.733 16 0.8647 1.337 1.746 2.120 2.583 2.921 3.252 3.686 17 0.8633 1.333 1.740 2.110 2.567 2.898 3.222 3.646 18 0.862 1.330 1.734 2.101 2.552 2.878 3.197 3.610 19 0.861 1.328 1.729 2.093 2.539 2.861 3.174 3.579 20 0.860 1.325 1.725 2.086 2.528 2.845 3.153 3.552 21 0.8591 1.323 1.721 2.080 2.518 2.831 3.135 3.527 22 0.8583 1.321 1.717 2.074 2.508 2.819 3.119 3.505 23 0.8575 1.319 1.714 2.069 2.500 2.807 3.104 3.485 24 0.8569 1.318 1.711 2.064 2.492 2.797 3.091 3.467 25 0.8562 1.316 1.708 2.060 2.485 2.787 3.078 3.450 30 0.8538 1.310 1.697 2.042 2.457 2.750 3.030 3.385 40 0.8507 1.303 1.684 2.021 2.423 2.704 2.971 3.307 50 0.8489 1.299 1.676 2.009 2.403 2.678 2.937 3.261 100 0.8452 1.290 1.660 1.984 2.364 2.626 2.871 3.174 200 0.8434 1.286 1.653 1.972 2.345 2.601 2.839 3.131 400 0.8425 1.284 1.649 1.966 2.336 2.588 2.823 3.111 ∞ 0.8416 1.282 1.645 1.960 2.326 2.576 2.807 3.090 "],["404.html", "Page not found", " Page not found The page you requested cannot be found (perhaps it was moved or renamed). You may want to try searching to find the page's new location, or use the table of contents to find the page you are looking for. "]]