From 031b7bd587ecb01c5628cc1a282ca121b06c381c Mon Sep 17 00:00:00 2001 From: aaronshaw Date: Wed, 11 Nov 2020 13:48:26 -0600 Subject: [PATCH 1/1] initial commit of these materials --- os_exercises/ch8_exercises_solutions.html | 1695 +++++++++++++++++++++ os_exercises/ch8_exercises_solutions.pdf | Bin 0 -> 49052 bytes os_exercises/ch8_exercises_solutions.rmd | 106 ++ 3 files changed, 1801 insertions(+) create mode 100644 os_exercises/ch8_exercises_solutions.html create mode 100644 os_exercises/ch8_exercises_solutions.pdf create mode 100644 os_exercises/ch8_exercises_solutions.rmd diff --git a/os_exercises/ch8_exercises_solutions.html b/os_exercises/ch8_exercises_solutions.html new file mode 100644 index 0000000..ce2da90 --- /dev/null +++ b/os_exercises/ch8_exercises_solutions.html @@ -0,0 +1,1695 @@ + + + + + + + + + + + + + + + +Chapter 8 Textbook exercises + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +
+ + + +
+
+
+
+
+ +
+ + + + + + + +

All exercises taken from the OpenIntro Statistics textbook, \(4^{th}\) edition, Chapter 8.

+
+

8.6 British married straight couples

+
    +
  1. Husband and wife ages are positively, linearly correlated, with a few possible outliers.
  2. +
  3. Husband and wife heights also appear to be positively and linearly correlated, but very weakly.
  4. +
  5. The age plot shows a much more evident linear trend. The data points also align more tightly.
  6. +
  7. No. Correlation is not affected by rescaling. Revisit the formulas for calculating correlation for more depth on this..
  8. +
+
+
+

8.36 The heads of babies

+
    +
  1. Substituting into the regression equation:
    +\(\bar{y} = 3.91 + .78 * 28\)
  2. +
+
y_bar = 3.91 + .78 * 28
+
+y_bar
+
## [1] 25.75
+
    +
  1. Let’s write out the hypotheses first: \[H_0:~\beta_{gestational\_age} = 0\] \[H_A:~\beta_{gestational\_age} \neq 0\]
  2. +
+

The test statistic for this can be calculated as a t-score as follows:

+

\[T = \frac{0.78 - 0}{.35} = 2.23\]

+

We can look up a p-value for this in a t-table for \(df = 23\) or we can figure it out in R. The pt(q,df) function gives the proportion of the T-distribution which is less than q with df degrees of freedom.

+
(1 - pt(q = 2.23, df = 23)) * 2 # To get 2-tailed p-value, we multiply this by 2
+
## [1] 0.03579437
+
pt(q = 2.23, df = 23, lower.tail = FALSE) * 2 # Alternatively, we can get the area of the upper tail and multiply that by 2 (multiplying by 2 works because t-distributions are symmetrical)
+
## [1] 0.03579437
+

The p-value is less than 0.05, so with a traditional \(\alpha = 0.05\) we would reject \(H_0\) and conclude that there is substantial evidence of an association between gestational age and head circumference for low birth-weight babies (the true slope parameter is unlikely to be zero).

+
+
+

8.40 Cats

+
    +
  1. Here are the hypotheses:
    +\[H_0:~ \beta_{body\_weight} = 0\] \[H_A:~ \beta_{body\_weight} \neq 0\]

  2. +
  3. The estimate is large (and positive) with a small standard error, resulting in a large test-statistic and a correspondingly small p-value. This is likely lower than whatever critical value was selected for the significance level, so the data provide compelling evidence in support of the alternative hypothesis that the true slope parameter for body weight is not equal to zero, and is positively associated with heart weight in cats.

  4. +
  5. The confidence interval for the parameter can be calculated like this:

  6. +
+

\[b_{body\_weight} \pm t^* \times SE\] Plugging in values from the model estimates and 1.98 for the value of \(t^*\): \[4.034 \pm 1.98 \times 0.250\] \[(3.539, 4.529)\]
+In words, this model indicates that for each additional kilogram in a cat’s weight we can be 95% confident that the wight of it’s heart will increase by 3.539 to 4.529 grams on average.

+
    +
  1. Yes. The alternative hypothesis from parts (a) and (b) corresponds to the 95% CI above zero calculated in part (d).
  2. +
+
+
+

8.44 Rating professors

+
    +
  1. Here’s how to calculate the slope:
    +\[\bar{y} = b_0 + b_1 \bar{x}\] \[3.9983 = 4.010 + b_1 \times -0.0883\] I’ll let R handle the arithmetic:
  2. +
+
b_1 = (3.9983-4.010) / -0.0883
+b_1
+
## [1] 0.1325028
+
    +
  1. The null hypothesis is that \(H_0:~\beta_1 = 0\). With a t-score of 4.13 and a p-value of ~0.00 there is strong evidence of a positive relationship between beauty and teaching scores (we would reject the null hypothesis).

  2. +
  3. +
+
    +
  1. Nearly normal residuals: The distribution of residuals looks a little bit (left) skewed but fairly normal.
    +
  2. +
  3. Homoscedasticity (constant variance) of residuals: The residuals have approximately equal variance across the range of beauty around zero.
    +
  4. +
  5. Independent observations/residuals: The question does not really shed much light on this. Looking at the original paper, it is actually a bit unclear how the instructors were chosen and whether or not they are a random sample. The residuals do not seem to have any relationship with the order of data collection, so that’s good at least.
    +
  6. +
  7. Linear relationship between independent/dependent variables: From the scatterplot, the relationship appears to be linear (or at least not obviously non-linear).
  8. +
+
+ + + +
+
+ +
+ + + + + + + + + + + + + + + + diff --git a/os_exercises/ch8_exercises_solutions.pdf b/os_exercises/ch8_exercises_solutions.pdf new file mode 100644 index 0000000000000000000000000000000000000000..d23dc7622bdc84c6ca40bed80f5e02fd89d9470f GIT binary patch literal 49052 zcmb5VQ?O{;vb8z6rfu7{ZQHhOu4&u0ZQHhO+n%*=R6X2N|A~rS4lpPTpV$W^u86SQKbr0ZRD8M>^m()t@Fz!-HvXo<5 z`_XM_go<~?_1KH_-hA0IVdV4M#rj5%wdvxUSqu1^oZal-z0%zK=gr-ozw5r<8@_L? z>(a{lU7I;Inf0P3z?Mm`AwWKfmXEV^r3E-XC$p)M~TpH^+^W z`=Ltv!LNf>(Rj1Tlf6;}&yAs|o?z;Yp!kA;zt?IkXc&WTWBC{0X<8p$;!9|WUC(i9 zl9P5l+ED;u)$t8vMF~i}8!0=P>*Ova=3REh1Ca*_=OmgU$X^;SH(Fm(b;=zs)&CsK z+<5D#UOgb1taaPlNSxzqi|udlJo;{E5 zNfg(W&16G(6%9zPEtF>OXzVKtU$rK>v~f(*5bC@>1u)D~4_mf`VUBG=F!q#6 z05{0HmFxLSEAiezeUy*MO+D+#+pz_?`|;=Tr5~v@jroBl8)6^fPoWKJdgaQEPe7tS zjn1$;7)RQT2ylJ2g*F*gnsxY#0mEh^(~r?R<+OB)*@>Q3?8HS1PPh%@_;l^U*)AtG zH>@`$QXQ1=BU|B1(!ir2(X;T$JAtFD9@FWF&%A@&{e!GTu5bSrvmxMbD#e&!oQgstA;3M%Ws^sjlnVXaGq zkIEC|#3zJ#gc8kYL$=WH>gI&&Q+ytBsh2cMOw)#-i-*kkY^Rtl>Ikl+joi@)_o}( z*W4?1iA*V6EWq(qL;B{sNiB3O6LJsQh$x17B+Ww@qymFcufeppw>zlWCt(%o*qwtN?BkJUOHV* z!5K&oq(X&Yka^~6fmuNCw;$;>VUZ6wH%uhNrW*Dys~t;gj=h)J8_DOtdVhJ_-k)3D zs=flo;k@v!gmE@zNY=E6LbypSXR~29S5vP)dK%4~hXjH_zOTT}#F|vx&*MnQp0{2) z6>KA?^B~_A`c7lmdSZkazLeI+=LZqEigpLA_B5#5{KeWcs0^GTAr=y-? zztfnpkqRbEbQy3C9}&)P@lAw@2sKv^DK^>|Z$G<&IH<4=JweX?iVGGtO0sUFJMy94 zq9@N72lKviQt}-H+QC_?!D5q}S^hXADb*~a;g9=Kta5R~V?n?79Dq&QNQ<2x%B)ED zN&JQxENKs_1C5w~J1nu{LKkXIs!U>(a~N%>d9ik72MtTz@=ms$k!N5@%cC39f6SL@ z6FaO<_FI!vgfUYp%MSD`p8FXp`#pVBG!Bx&GVKO3?7#a_r-0CI#@EP^3wg&;9t4}o z-)opE!0MTKo3x_^o8peDo1;d)D1dKWfYX7Qc_qsi^+Zs#EXF>P6fC%xN>YK@z_aH>=!VZ*sk^2sZcC>hLFTSlo*j;@vQ+Tp$r?C`$xnZ1i;l$I zK`YkwCmCKF)LMpGSok$F`FNpn$p&ind&?ikAOYpjb}=-+bwc(&)Ao*GXldSwK97`w?3mV-k{ zg?2a#V-v_4Cvo>j*iuwr;Odws$^8jzWGxRr)A#cAy0<$SokMhBVt$bf58A~1q{Kd8 zgAsdyfY1{9A+xIYS)1yTVzFIl$)lQ(m9`Dg$r*=DY*2;`!!L3h0p@HIPUim5ug7;- z>EI*$YD&$dtYC6)?9(B@9wQbLkn#!D?}76i&sNsLd_pKU=>o!ocOR9^4^(U@p!B1I zV2P_EZ--{69n4*(9(57)?Z>`&Jh9|{?eZdep*vLS6;z>%z30qXyy^Q>c61UvY!+hx z;*mR9Q=KU@j%#rN$7}W)c}%#xj1C+8Ity|x_cAF|5zV6gCG$_u5pw7b`1rIB8eDWr zyIC!#Emipa`x&?RGMH1KYF<*#NAf`svj};>CAPFx`{HQ6*b^Y)Ko!ak;nEz_AJY0)%8cy zf*>CAW8$RA%&QL8fxftZHQq>)O~Cicw|$?vp_TfxZ`Hu9t$$YmF}>0UYB!Z6HaUV! z()KAPmx(w_nm9908I@Ei*Ju{M{Kz1s5d@!GvpZ{<{*~57Pi}1?!GN?IAGTpvZ!?Tv z`6}JgAgvAF`@%czd7_HHuLZ+=yx4JWJTI4tw)VNBG(Sw+tW_VteH@&!ed$@+Uxm<0 z|AB?Oo=t{2VJF|`sK2&H`@p|lR`V0{_c^iy6L7Ic3y;~pe)CFh*vZkC5I%$mMJ-5a zt73B+-A8z~S?SDv*k^NsG9|y;OAm6suWDHe1%=RhO7jVtI7D=818XzOEF%+2L{b^0irYJH1NcdY%_L>S!DHNl_D|PZr2p>~3D=t$t=(^kU2p$o>`zb9 zRcP<;n+@Yzu9@kNM|($c)v;rLgunO=?5Z(Zca1FmAVIZ`GyeyXWs|Bm)C_;za6SpYP9OJH2=W(%D$k=Q0*4 zao108<=GrFi#)7UWmjiAy_$KEx1+%vuB9RkRVk~nQx2Eh6cFjCr9~K3^-C{Lp61Uq z126-j(d}NZR|>k~A)0NlE3RLBh$MiJR5C>Nb|Gy!|T<$ zMhk+#(SUM&>L&GnvpZ`7!4n)eF8ZfH_To`(*V`nu;m*h%ccQYyif^vZI*;t4bqV^q z8X7~vE;zgDFxL==7zesFL1Xew{8G-?vQm-KDYPV;$=*+Ft!1v={u)(?h?-}I`bMZ8 zpj>4}6*RlrZSy7Hu*@7(K2)FH*mOJuEAWdBR4c}BA9@$5y0z=HPZVA+!8vJaN8dAb zOs(RWJYujU(N6d~1N}`{l-sKjIcU?xbQYk>#C*0NKDt1kcT5K2A;I35P{}weo_FpF z)be?!k(vMPxSu$!??V;A1_KrE;idbdEiJ7zl?f7Dg>NO}(!MkodV7w->TfCd+%;&1 zE5!ov^H8YebFg{>;-I&y!Ctm9(ME!~W;a7bl*(bi z<^M$#N)#h3dhw^H5K*w4gKZKoPhF@Ppu;rD-&#-Ft%su^K$cl}^wGjgo$H<>X8Q25B zuqD|xd4~F1#+)CUhSV zq@0aMPEcr)A6tcEB0jYPAq+6gI-gN7t%=E8EC_tV128~9z(F#cX24Bd9;t_3OT%)L zm1kn1_EmI{yx)2Po5HF?f52yUS6i4S&muq5V_;R%fei;=qpeX`0vi4VA&4|Hi3Ch1 z%=@E9xr2~ViDQut4-GJ>v1TuC;w%GU{xQ=Ik7P#v1;22Rl8DK57(}7u;X6Otv#A}@ zNZ45M#AV#B|Hxrp2|1_MuioM0BZf2PbT;;n5}3Zp$PwG3t0H*B9GLBOjTbUbXFE$7 zr8E@2eeiL_-+6F0>jpa<2gxLJ8<(Vxeo+c@{!RAD%ZhzgKf3y=`Q9x|IexbtR4)Ri zYv)VhP+X-;_>=VeAM4V<0KO~C=YX1zBMRcD=LnkXwWKlFk6xXHhrN`Mbrm zD`%H86*ec`4&ufl0D$@Aotu~_8Zmq<0H42p@vFHZ@sFYC7yvh%O=-Fm2VC$1Am_4u zzBD-B3DjI8R#^(N^0fhA=N|hybG=4_PUHScePkqJ0 z$T@8C5t|3D*CZ(vF61w>Lptz`+@jaW(zcdS0h@z-3f%`>U;*RQbeg&?y98WZBhxQ1 zmb4eNr#02k)eiDb+Do-*h6qE=hGu9-{s@dKG}|!13$w=ChKz9L1j-~5_OcmWRO?yG z6sv9ssAvivJQv{@N$P~2km(|wVbWt{CY+Enh=RfGdV>z0vCULg^_`5E9l*h_{T98p;Om9N^I(~ClZl`@cfhz zbylr(CFPVCC-5!=vRn~a_X_v!QUIQFjz9ML#djg?F;{caU&yQ;D$fLS%~)a`yZw9oKVK)T6v+?4$9y#;|l|rM;c+oF2rNYGU{>FeFPFm))8Rf8x;vm6#sC zd+yCauUsdyB_%Av1+#;*7{4Y3=&(2)+U+!?5o2@nrCK@i0G4XZjs*81T#ZUTu7~GE<4% z9n9OG32c(1wkjy>HM{c)E7o=7=c@^*v#)F%HC!Wa4bphLF|gvtf)Fs`_DW z671rCroXBfu?Dpt4^S$5bV_-RyNM+nd5{*EDczzTM2|{mkd&M#2%iJzemYKDwA^Dw*BN# zMXe;6fX?nP=KhGep%c9L$^6qUJOAX0Ux`1GJSNZ5wz24rAlAn9)=ICCt4cS7LigG?7r(%>NfE=k-s(L?gpj@U1?g1r8PWM@BlCl3IU42Mszp*WSJef#}M{RL_r zai;t~(VvOo-=IGOJw5aPp7v&|wZ{;#!B2mqd;?oPnj*m?qpz{-Sp%Ac0Kbko9_jz` zF;qsxQd<_eOr}|gTd(NU!qV)bb*6oorCwe1YViJ)m+D* ztbQ79TRN%BK8?_=Zr++16;#f;z8`*UJ{)&=eXZ$HXn9I?mGkEF-stx3@{`vmHE`alCtPh1^VvcBLX=4H72ywBTjX_Z1PEnKX^k>(SJBHAZu9?OOzOsrlQY64!|Y}{gux9<1M_X08R#&8E~iIEAJWO zx28*Di|(f)j2FmF+z6yX9ylfvLjhR?XAW=sV!xV^|J-Xd_|7^xzR^6irYJ$|OjFVX zA`jv2+EFI?(i045_+4Y?+59m6uKKt882n^-6$G>{MP)am*$sTWdWHX<^h8Ocao6vi zZ+GNt2#>nTN&|dZAgQA#P2Wm75<0aT>qh!KqeE{}fD?8-ZIZ~ieHnm!XZkpdmP(np z*-dJ9fljUa6}ih_fP(D z9|=+jJT@}+4qHRBi&lu^d;)93U8d?}&qY`XZaYHTGCJ66r>5tR(IEWK?4DCb7$WA@_L-lSc z=|WGG!xGk?+!9SI19vBW4Yd_>;FgMqm!J$#d88{IIXC~LU%7;ItBgY*-^%G0G#AN= zU}+`C&;6t%jn@@*6rEaDT-*4~GA`Ws2A<-{5S^M;*2xTcs(s5j!Lqto2D`pkXpdR@ zNxGvhv==q5jH)O~);rHWw@1<5&W|5t zus(ktCuB9RW6{D|t(c2%>^W*TSyTFEw|0F$&?c8vi2hH$S_Ho3Q$Y*$C55z*Y6-mn_v4mrTpSK7?9TYahbeuYu0n{U z6XHvb=MIz1oRfk%X0VE50A*9w(*UI4a8dAP(yN5Km!0u>#uMb%OwGp!jm~a#agGV2 zq!Ov+G5UP@SUUs)+?xE;TMc}EU^%}h7U?C~`p+06TpwG5ZPpV`j;mqYK_k!S9FPl{ zmU{Rxg0|QJK+X=Jnar*GBS}SU3#rB#Qe@YDKTaEO3Xr^0wTNu z4Mcm9b7%b@M0ewBL*es-1{WB0T}v$kvo~~t+-6o7mq6FT2AykS-XBmF@=C`~xM`Nm zga4ud>%Y=~gZ@7?kh2QphabN6jo@iX+RI&L0Yy->37BglAmw>7b#4;WHYp7C{=DcR z38f|fay8rO%EsXuUyyDgYJ{txbHdcvo{*fc&Q+GU5|p8$Sv6A8iV1n@y^82>Nc3Uw zM#W+HQ?WQK?v6~G0L4+vD0(FNK#m|7j;kq6MU zItYo91ZVWHlIa1SMmAQ~6_7_LJGI-}5XO`dY76{@4Uhwb)Sg2Pg2Zlx@s`O44loU4 z0yHDwB*;U|Ljt1ACwPON3t#WzWdf10Ro-Lb*8}km$(vh3!dC+1$z)^6rz4=W7vpR7 zSUuiP4flbz=AdASjm5R%Nt=XjPQXp%dm2nA?&X6Fp9uReQnLLUDOu=Q{zJ+n9ji@x zZO z8>Lsvhm%dmQ%{$w%!y7ZujgVVFt>_KatRfc)Nc8=l+%Jp(#h3=PmQeXj4$2a_n(&E zm&5079yA$5zV6P&=gK>|&sZ5l<>Q#CpEeiS;(^WO8`tw~k4znUK1Vh0jb+ zb>-`Ro*1KE?_=-ntffKUb||7wt5u)&Mk`Mgu}2T_RQ+0OD!229h4%uEc5_s?S?y`o zL8e%3rSDfgIS4N7$=QL}#l%ztPwdli?m5Odgbk2}OMZ)yUSv)9Q~_fKLE~}MIm$VD z^sp>axTNB^5c#njg1`WOYNp_vbP5M0amNTWsOoV-4iTIQqpa$@pmZhiEQn)xstE^E zUHq)$D!S2h1Vp+5ejUe)bR{&Rpjl|1XtYNEyM%cZnPg|^{V7Nb)u9h&UZwHy*`y2$ z?~*brZ%u_f$qZSn3gck=lJlnibQDA$Z7gJF9fqed$a*iva3LNf&g}cmud+OyG2YHwmd+10#Dm1?G@9r{3wA@tl#xN!9mH3u zf-X-!>~;I(phVHDmbVqE*E2qhpDM}@`U}q)_mE$Za z1s2%O0eA$BR*aN_P#%4=!MfX2wRVsm9-x1yU5 zb`2p@K}kD$=R6>ePK~zL*2Ay3P_I(76;mzAKl%62+JA8-`|X%^NP}mDf<_%Bsk;DG zAriVdYC=P~M^wfrk+rO~jAZx-HTW??!ec2y6eUgDgj|}|;T6bDT}HF;!M*zYYd;{L zKDHw))9H|D_+l6+s8taY7qn$0MC1Mw7;!q$=)c6A05NFrKUJ@ zGGaFboEb&L5R@p&pj84S3nvs%%9tx4RUjf#^4tM=QZKG5R|heR;YCSA#nLL_4Xt`K zB=r&XuueQBN}L?A%l}FaFKZXdYa?e+5eJji!2n{28-^3#8v}!Q(-sL1;0enU4y%bb1EUPxG==t7pnR;&F)5BWm;3SZIGs-?`(=2Wo^H1Q*nb2lYIG?d?@<$Nn!jEG+c@CTEyf|D((LU(kfzNi7nsM-w)80eQpL%yTt=v4&y^J^v3f zsXd0^vgX|xo9-M2gWvFk4i_WueDa%|*r^>3+n#rii{c~L3P+6>> zi;?hE=A?Cd$?U%H&EDdL=GO*Iu72)sRLPBQ!DtB4D8{-;L}Af3Jn+gcp%vh$V{1jN zhiZ6aEoE1`PTQB1n%qjU6mdEjR4Qp+WIvD|Cwn3yIaLwQJJQ~$!ofs0R7j{C2&!I9 zQ5iV!)(#QQe?>mROw8DPYn791%bYroLTnoo%&@#;9}$?Eho*sT5SfyLy=QQj1p$In z?)*U}Kn-XKOXn=a@Q)(B@Gfizwc)22TQbRmek!1_7vLCMMdrXSvri0v1Zso9F~%G? zN7$Mfd8*zg2$BSVBfu63xiifvADkxak}~e%>QE>llA25r@S4)_QL?n{$`R#JlzQb7 zi1kkn)UJ4X1+(tyl>0x6oaJBj#mo$h|6#0&0{=r;!T(S$Pe~HI%l0;iusAbdzJ-*` z+4@-;q(gEv$I-6ZD-Dmb{M-Kk7RHE5hFQCQw>{jz5sEUgHr147K@Pvm z`AGr)nQFdP`*3z-5UCsGnRs7wfE6eTs2Su;Qf8iA>WIx~7w_Z7uZ_<$y{V?)r`_6f z4qN^5)9WYnH@(^VE81n|s?vWEnDO6KBqQs8TzTrmZCeJ?BahztL^!Q8M&Pke@&{G0 z2X(}(^YTl_$d*Zj-4C?pN)(}S#RcuYPrZMN=&X1z?fs<^&{~>{&7#N}E+7{S9USYX z2Vzz1PV7DxE;>bniPtk0yyM%{N>;O5EkV%Mc0tc_yf-c#?S@3sU+diZq(mslzHHDj zi&mQU0<^dQeB;x))Tl&(MI+_p#@v-*P$t6LyCc`Uv_PZCWhYpLch@l848Czvrq3eD zF9XM$lM0+>6tA5AXR`H)Ca0~kx~Nwm?P}c6dYXU{;YLZ*J(@D(UkHabNkbA|I1wM2 zZN%ktYg>XfubVM1Qf>oQS*0>c&FJVIP8ZCq^6~{J8t^yizr>U2U&WL8KMO{_-*S*1 zy8DyL`3jvx!eE;OM6yFre_J$5*VW0IgoU|9Axr4Tw;stejsFi%FWET-MWvHTQ|M7> z)zMl0q9{U)D`EJYVJ*xgY+loo7K$DR_-Mide8QsL_~ps-qvobS+N)3W_Y-_m zFH7nbIAmw?))ef1-STq!Y~g4=V|FmN8NfKrMUNEi4*Jw&eovF5?;09` z8BBIuXRnoE{Gp@L_NRYMF{g`l^8@;0ujsiiy&J2SMF$a&jM;4o!D0ndCeFh|jxt*{ zoihnbUEd@l71}`773PHt&G^!V^bRF@Zd7%|{W!Pcbgw&h-j#boZ*%&>;By)-h(wyk z1dW3{av&i*qX+bcdW56Azq%ns3_pdqANyE2cwD8Hx+IXfEx)V@#6cd5IiTR{kVOY4 z`=J&*r?7o?NEk}qPy#qV1>-)6NH6#{M;3&1z+k!`X=p0)FYiV(gP1*PCnq~nTBLo( z|NZd0)iY$)BWH&v5IhW=jI%KB4E#ZGhypW`Z43|a)ag&Y$`45rlb0U)q>g8%c>9RmNI%N zI|G|PJe}6lmNcD9>gGm?PK?jNsytXogAJ5zaVbMYUV+7Xw2l?SPgbWPM-iiZ61jMT2=4+J(MEku^w@I`Ku44zr zHV?@Ud~t_k3QK4tMYp|CtgySVAy8JlD%p*-gi10o4bgYk0)7SKm3Miz08cg(<`emP zy2%M`7QK}{#e=POcH=WNEX&xZ#x+gZ5V%EhxiYb43*)QWy)IR>^2AtImroa}FVd5qwyJ>qY!ltJgQui*t+(z;ANM}-*L`$adOp>6KtDb^ zbMwlj6SRUC?S+F;|!L7b+_7?g$?5MhbZ;RbYkXTUFCSs>g zNXUMu)n9{RCFIj9AXaSL_)lG7^;Z}UJ*4O44&SiKK;IPxI(G@z!I&5pCz8DKD&u!3 zZl>}i8OPe)mh81XP{dNDVPed0XxCdGyljeh~;HYV-3_F;x;P-b`um{R&+9U8x%Kv%VI$MtfMw4v1&Km-U z8G#jp5`tlNUd8R|xta>Fa56oRaGy|fc%2Kv6~6g`$TweI0no&jca6@gFt}*R6LhKY zvf>?nu(+SyyLW$BZTu^|(2uu?i4CvnQ1gIUS0$d=m@+eVhRw|IE_8Pw!NN&ey4JyY zW}#TSubhO+bwu8=|Gl%jc$+0Ypx6qYth>R1I!|!U1p?P6l&FjB8+f3STZ%h<1+}8K=1P8Cb7Hod-d(5586?? z!2RsW;zxT{2DVr`iM=zGlXx+yCcD+ti;xr=jd4WAcTFdXliK7u?L@^w96}7 zgsZ5hlnpXfRO$r(!2B1xgEL*MMe`xjX>|wZnWP7z?y#38O0b#mjU53)=8ONmb#q*3 zNjHhgu=k|m)(nJ0B74*r{3=ThWH8nrC(!GBhG=fNU%poo>865H>YvlimV%<8#n$3- zYkLE>0Hm|c1vxfC$YC2p@L9-^HTq!Zk_rR=n zrj={*?~9JH9})Q!oVaj~OKcESq(Tk3+Vrb>07(W0Aw#m>j3^@U`vX=Hhy<@nidFEr z%v}mYp|D)->=v885%p;H>7QjW9=pd;f!!Sc9fFTHgUI~Ea}7A8M_ z)5?GeuE8-rL7vY#TJqDsNZo`qY>VKw2}|FjhZJvr-?arUI$|?I@Km5>gO4RLQX3mx zbEyGLuIqT=D^PvcSna{(5nOhFKjjwL=!Om2(@sC*#MlAeFZ>!#aoB^#+{z%%hzuGk z*zDwHsM?`mk_p?19V@V?d|glOy>sES&v654Xy-95o`;X`*aLM5dpRD~yu96G}|lNyAANd1lJJw#lGn6$>dMeZm!v$3lDrLF?@nZ%R;1eOs4blM(6S#DVnCG zzEYRFey_GapP3%G^ER%JEzd>OJ2M3o6Eg0Uf;Ti_vsgB%vRNQ6`uw|KZfKZcXSLfM zXS%1ixzeaiu32ncMWj5u6k0URbdb3CelN&e`GeoH=dSZ8$^#E0FtBAdwza!guU|X> zkIy(snXWm^dBt`ycpp>c1%kiO*G*26&!$u&a0-K*mr;#Vc+ivW`b{R{!oHdkQ+nIqO! z!jLot+u_r}+>;cmzHHZUh4IC|{b<9%oOXILv@JcQ)Y95do%}u!CXG;w%}pKl=>Y%5DgC-hk~pFkC# zaCouNp)IV#rA#M6aWl#3B=_CkDO=00lTnsc^bMxCg(q4fQv=nux0aeVVaklS&l^nn zp=>4sa|q`Tq@830$O6T=l};Qy+}M{{9@ism^anGZSx99KNOzydUcxY;L(CsGVBN|~ zg>as(t*!sy?R=t78uvt;RI6Or9LFfa%5UO@*H@x-Y^0aXVOa<8i}4cm6os7%)m9)) zEcn%xVPN?D_+bl+H;;^w=^wH460+N-ZeG81@Z@}K{EUo^>DRTdDBV80UtDy>_r9TJ zC!1mZ7ijudBTL5rm{?lr>+iwS)9Wi6-s!8``vC(WJOfSrpPy#^UsIO=DjI`R31D_yJYyoG)&1z&OrsZ`3EMKK+Q>$Ko0Q`@Qp@Jd-AtH_i!lQU3j>_V9 z=5;fsdG8k#cduNM$Fq^3381YA&~c zBI*xz##8>Vwr2+&(#b5(`7KvRl11)_()YlfZXdnuOXJRCF(1DFFZ9bUUeC+laj!$* z-*ohwK!-Yjt620t-2q&qLhRgMhRsYwlPu?G;fgt8%O?VX$CGZ%66yI$L=y;1hp{dX zs1wgwsUfO>AuFR-5(<}#B4LLrK|)BZOa+Cvh3&ARG-Kr=wF^Sgj zX3KI(nw0XS4dmx4dxt!+v$}7IZHhsNobNs=DPh zO=}!f|LRdzgLrl17NpI{VM7OXpp-BY14n8MseWa3%#zcFriM{y>yT7^J?hBA(ecx* z`}F<8y@$5%p6o6)4=eZSuHFZ@eW`mT-`@kZ%hzL&_n=OsR@nnN&J(_6zMbC_rv2F_ za3GLdzRvIU%BjDQHv(`Qm-(&d>_qWO+;5sZ+&O^g!{E7U7W61)aTuhB{Xl)KqY>m#W zCWVI{*nVVOehznXUa^t)cT|>0xdd3DRsX`Kt8BiD{pDnNw=2u3Hp?)EciY z_CJ%?llAdHhZ2Uw+oXx-*@)8$3-7r>vstE6vc?PM?hCU4H@HK7B`$Q~=s)=*-F6K{ z+Hn(fixYl@hKhp=(hMxG*IQzctX?PDkrg6wYJ>Id|4~;a6Les&pMb4KX;wB{l#s() z|4<@D0tJa8-VEo~gJ$Qz86-b}P0fiWg_wNEV;HrbVpgPTo*$n-3Mu|Jr`7y8zF*Wa zgOk7%Fy)R?TE#K9bXNAF?S0-+G!>kSlm=jW-_!^TEomYAWyD7xPCvFF{~!PhUUr&S z4ZIW7Y!!~Lk&$sQ*RrCNv>sMA$VU-0pC*T|mfCJJV>6-8kolwFU)`t5>W#B+GaEA1t&f$%Ml7D0rK&{SD+t9^reAs@w8z%b!kke97{IorkT7{A;rg*j2A^_((DSny;c4CuN>LDa+4hGOTm z1j8gcUOw&!j(9%u`S`lU5zw1ZAr5Y>6JTs%bC*AHNW;xQDhq%?AyTki=nZ5}W5I2p zqLrH(l*%=ApxQWB7?>VtfLfjzH3|8v`X_ufKnh82L+@JVr>120*97=8fIT=+8{pYwoiKJFmus&T^5$L4#4Y3Bs}`X zE^ykE<3^E%ugkuwD9^D70dr2su=31vyr62!P|frb83`9lElnI5nWeQkqON)pr;sqC zF~(zv6$EIqny;UMEuXq%YaSRGTYP>Jqgy9<^8z-VR_)UmZ>HhG!NI8GQGn}_SdOrH z&@&`d?Bm~2VO>=P16vLa)8=JlVB>!l^Qa5i-tI@P=p{?{w6OThQZL`PtU!RCNUgZ9 zbMi|w$A8pq+L;`m@UG7e)}j*v_e@2bm;$nBbgCvnPw?J{UJsXIF69<~Eu4N|HSBL8 z$`>=v-etDpmx-8yaJP~=Kq)1TQO#vQA-@2(iUA6(>`xH3rz$5Hw3>*K%Jy>2Lp*5d zH=?P205fs<*$#@^=af?)7I)EO^cN5+Kfe#NQizoyxv)&KGio#gg)V!Gg-?n+$-o@w%WjE46d! zNREsyy-@H=*1iwA3G3W>Ytgw3_%yTxUipgAIrz71(-Z0v>tB|SIw_0_Oi*cann#Gm z(9x^90)24rfYWeNigfe*kSO}I<%H0f^_StKVlbh}JF<>y!x?a`|IjeEtM}`-xXat~ z$apqn0mjnhZx|f3m$uAk)KLxS7ttR6T@9@q_VDUE`}r$638`(uN9$qvs+1KC-Cw*gc?c^NK#+xzv*3f})2O^!97`THMlB!OYEvm+-6(l!O z0_NiG4PPro$AKlu_|Kmsj7k^*&rUOa<0p&Rej2}CZRmcsP>D~F(IFA3649Y!bFH_H zKYM?l=~nH!8||sMY+y1eQoXKrq&>Ke+gy%N88E)3kp@ zrraF<{F(3Q=Yw-Rhejsgu05BrC_V5pq9i_!iuJl^sYqKy6H%Lm`Np7eWsH15RZ3Si zCnuxIn}(H#whe8uE(FA#DSWZp(^+dTafBPUPYs#c(r7(h$cx7w<*0mC_{3 zI3e~GY@Kr&fNN54Pa{-CgzYUal;_I0F(_#?6e2c%KjB}`@(|_2k@bEdSsXEITV?qc zx<;!}0D7bNIOaX@0Q$Ju?E=2N?G!qy_gTggq}g@HXQTYg)05RDmL^Sq_EtN5kED_D z9ip?YrQ&rv@j@6l-8Lo}ru#|`-jp2}2uaBxx^WI=90i#cy+9i6>r5S^o2`P@Z91GD zbGf2rvfyT!?gxEc13~s!erNqye^$;e)NoN1{;imy@Owuv8*)kovGuM83s!hK>g^=h zEz$)-F@8WaR{<8`;kESx%1RIaGQ|h-*D~KAia^~vwG-@YmAEPS_GCX4oWfb$uirtR zXJ`7hPzgSaG*qS4Q?#$h75wY*V2q9_1i*x)2WL|KS-w*m9l7e`1!7PB8&VrxURBY>CAX(hYO+(Xwso+j8UN7uC_ z9{y%*;(lG()3rrp2f=DLCr zrv{+qYvP-vLnO-K0()tSsw`lIFuoFtbHvZa%pcUlXqS$SSOoyN56B2BtP8SXN>$wo zI97g@s3g*D&Xg%79wnFrMFe2~ti|c9s5xeMg7TuHsm8!_v2FlquMN;HF9A!q&yHM@ zFx69x&y;fUtQ)yL&y0~3_a1T3r2j~55Rsg)K$BCy8%_2CmX*iUOAnozCQc2I*=osb z*4do#h&R_8Uj$QVV(}WjFQx!iLCG^YhtLwX$b_nB^x$(A=N8Kn@Jz-Nru{KB50B){GPL`R@Z1dPw)pHzgBd+zNl_GQNq21#oW08CLOAF2FaL2dQqx zQJuO*F$Had`u|1QJ9P;bZ0mw)+qP}nR;6v*wr$(CZQGfZwrzIZx4pYh_vrHnB38Vt zG3Q4w%dhX8T1%^jnt4lFD!K>1P(B`yqqxCeAg696^8aBG|Nlp>{+HU}-;N)TfS{7z zt}h?jAH*Gn13)?qy8l0%{U40k|8~h|W9IlTS!#?1xH`((X0AyFuCj=Ub`i6wsjSSE zj!}fzKyG=tKWX7|$y!JW6phJ-qiY&j4;@V)Fc6VC0;*9&j7SV)YB?h#qimA-G^S(w z@8et{<4)FABQWq>M$@VH^LOuSG(=O&hzQX_IqXzhwdO6K#mr_5qQTBSbE*6W1-wG2 zCv`Qv+FjARS^)15Gsaw``(|ib#5-f_nXK2nuh(BbGiQC4+E3j5I>U4dn{fNKfzbQ@ zAJ%<-(}ZvJO9s1NFtK%gqG?1gP@{&1b>wP!@OnkqpI(teMUh>NaixRJU%ZA z5#vB23IYZ(Hm8p#jij2_cc#;dMzNB(bDKt>OvF1NC*?-r{D^vw*Gu`4eeke5wb6papB{g1xNl2rY-r} z_z{pv6Y6c9ZDVSpK_g?UK=0%I1;>qKN+ykFZ93@ahzM;48_j4Yeam1X{mecFa!<{( zCdq}R>fNlyL6e-n$*fgrbb6m=4pY@k`8GsC6^b_Pxj9_!v$HdnaTjIAI}gV_xx>rmFDtx#U@T0u1|peK{ki5+k@*?F_AmgUtJQLj*VWjkXAF*Qr{yZQ zBVnUhs??t-zyb4y!_A}3Pd2~X(Eb&h)7x+0s-|kF?JcaVrKqy#3b5JPR#}&D%X!F} z_(0Zbo@%igyh;V3Wb;xgwj79yeB=##byH<3ZnAr^y_M_yu%0-j@%&6rW|7gwYCiPI zJ1N3<-U@hh!UXTd?!;HX8AtHNoxu6w;31Hu8}r%OMQIIrfu*m5GM*X|p7OC!Dh;2M z_|tfl*xt;cV%XSqD#27nf@pVpvEvlIvm1}GHciH6AMt777*Rsk*FT-uSQ%J3_<+Q1 z0H5-14bEc}`hy!Q@1?pMFv@LAt_LVw5}$nsA&KVUK|;cK#}qJHd8csT=zH>Mz1xOX zc#MPhXS3n?4;Gx~1A^3s2UtblEU7m*G8Q$wYnkLPE?`T|g0&+ngNnj@SUE`OAk&HQ zU-0*qr-D=p4_d8Lqs=p6&ZiIRjvdPEvm+xu$GCPkD4A3=|85jig|c(ZC(lr8KBR@q z`P7~}3?4Z4TV80<*Rmq^s&!-H`P6xne&v}#mO`AP3Kp#R%mHhcL9MSgoE#BibQbU# zA>DEG-VoD%7yBFGaUt&20KK}#{!;K7w1pjR`$u3#erUKvv|>P;PPw90#y>;5vd?Ag zGh6_t{T*%Goo+d6;4Dpinl9>4s49qRdTJ`Yr-vgpxA_K32%({c`*EYf+!CxppFpF4 zj(e2R?o^$x$}rpR;1pvzkZ)loG(A#;DdzMWYo>;}` zg@1j3k%@{5CyBMBq#dE5Fb@aAU2^eUEn9T!3YeT&I7Qj*gaXfB^CP(@2H-*P9Sd+{ z1Lab9i?i!+2N0*2A3D{4AO%x1GRxl zS|J6>q>RT~Z49`Rb)RDWJAD?cmAhV)xmPT|>h}?bRBo-f1lirLwgwtPA_UR|c$+NL zS)JLUE^hmua!^WZEh{9E?)m`sSuDD5(%1O1(ds+0EPB&Xkf=LmG=9UDTD>li>o*9k&~feCOQR1u5AWl9N}QYKYON+>EYnmBY} zOI2BB_)-gJg(cne*OX^0N@7pDdFi%*XBjL42qb)t zmIXwUpt&GI5)_flnNQA`PC{wPd&}XsS<8=`O?}1qYE=mAAzjbBiGGUPl;yCy=1HRD z+eYjdFQPBX9;a?@ib&f#>Yuo#R^xG=E{Xt7U;DlTCaUs+uq)!!MI#T~P#f}r%%(XW zzS)SNj>}4|X$}DStfL2q4Nh*Tl5;XP6d?WibeI=(a|X8ggf8-tq&X~>pbpUV^>2|Q z$ir&Ce=5jq#DcI?vclP~wi0L$jZz^gI6!{8Uwl{=>*J7biDMTgQb}rSla!R#X*j@@ zZ`wkdD!r0{fGZY?I@ZQHOvE@imZtSN0>j9at+tJIg6?SJ9$TxkA~{!BRueL>D$F zl~4w&0)C@HUy0d}$oqU_uPk#rbW{miT2}$}cy0N%tPkYSP~CiqUE-qG`ZAbZ4d;@t zyY4wZf|Ixd+-BgPqJxlcY@jDwBHYFOtG9j43ItuY+lUoLPr4Zcx{kReTLk+IeQ0N6PpR}2*!0{jiW|2gs_h9vY= z%wBO==-LT71ImucMfDjS>f)33O06%p-NLiBSG6TRqVdA@&wS0bX`ihSW&@e0IXDWh z%hy30C~q++kYLkFW;R_c5k~sK)n0^Qor17Q9E; z9D4DqTA(dwlteqy8g`STtny0=JoH1x;}}=;3K)0nd}x_J1ebozLi}_nQpyLA3KVKu zc4opvKI1sxsS!lu`>hT+lVBElt(T1U8$cgl=@4r`r2O)*R-O@AcY$D?6~X}~`T-|h zY~DLgmqAoW_xNsbh*HQAtV*LoNb15P>5((MjxO@lE#XUG;%x?*oKWl?kMZV^J~I)pJ*Q{LN2H$FMpq|K#*=m zGVweBEtWnKn;#4h(MInGi05S5zUaB9Fi*@5q)mLX>iKhtbYUcp} zGLHVX3(eWq%)_yWT&iGR%?bJIHwTkTnqZ!Int!)c>MfCL(CpZBQ?}oPobgbf{2RJy z#}j{M$^LGIHElR8{LN~1n?JaHs7c;P!{Oyc-Tr*tL|$hI+EFV?CgT0QOxr86t2-kf zw_$ghg88}{a{l$kOY9P>3`BsE?rtyP!gf1*B4_%p!mMZ-Rnp+%w`$JaJCGFU5k9Q& zx&?gkV<6eAAqb27CfsaTBV^5 z^;902WdLKShut*|f>0QeikZN+L#w=irOgEl*>I$(Uf}oB5B=vtwOc~G)91yW$85tp zW4oCrj1GQM-}B~!<@wOPUV$3s)01e|lir0Om$rz2SGdm)GPi3mI0BaafZ;%R0Bhop z?RhWAERfMon<{v7L{QSq7X{tM!P^rOO_wm7kDN6bpHOR5H*boA4E((~V|!xRi}5Fr z-rM#n$UT3N5n)4=#N@gC?-l4SW&jv{-sKT>on@RcNUBkJp+1Ws9RKS4ihXpirz_eZ zD5h9J@cORXXarsk+E5q8A4{OoVMj37pjbkC;A|nmQ=fTExPYPb10sDEMGvebE}Ky) zky^7>f-Gf{TxnV=ksy8nj-XI)_uFkCF8064uZx2WlIk76nZ@kN_Vw7qA#^Jyvm-2; zH4Eh}0%{rkd^mnm1cHuGv@hr%0R`g;l9HfYv?a5vM`+&Pj3jdmWBIa(D7&a1oANp{ z2pIvGnbUfixbnH&Rvv}{(6a+EXb9`Bmtc5E#? zd@U>?dR;bjSb-pZ4g2ehiqH905a&@mvu}F58Os8G?opG|MlheePkz(Ab4mhah04IE zKp-Pm1HPdw&lUF@x`NN5n1J(}7ZzLO`c zp{}PGuX8}+=1@Prt1_dxcJ)`iigbmF+_CXtUn+H1$?3I1*+qVW&*a2-`LWAC-m`tO z#9eD-g^1?sasLpnscUeoX&TLUrD_P?BU%dyp8rFo$Bd5tm)#m{!&*}OYQA<_Yr+*HMV|<+> zEKSsUjWKr7#Y!H@-vqUp%)gS-6|@1pQ72&8A4Iv7s~1g6xzb#XdW$$g1&xV}SOE=% zi?g45gg&gW5kvSshNN?If`=r(k?C$8xziJw%`CJ*y`TJMJk+4ZkRBb>J@q7A`TfKF zx4v$pP9TH<&*@^}fi4a9$mEs{{a>F9yG-8Bb;y09Sau~{3P^Fu?V>=Z$;oU3lFp*I z^QiSJSox6!!4j}Q6LpNes*cal95VR6|9}$jh&6!;U8^u{P${ZLOrw2U1<4`Fk^ytM z2-HJmv^Q6r>Ak&7jxq8ykVsjL@!frN`|_vc!5Dj}q|Z8_-DKwYpgN~B7c0=bn)tKR zun98rt7MP~p8zZp^YOL{TD6NhZFe!NIstBuRxr4)_iRy+GCqJo)@%f+H^cL5R(9Pw zQ`=iUpLa4Buhj`jJ33bl9P3!Sy~=4=ao&Ep(x^*H+k(IR3bPU`HsKG(L#?0|sGi(D zV$F0(v-4&bOB5=26lQQzwr#c9{$h0yuTPpvzGU$>vj@|k<*DP9cW3d&bLYFYp+skQ zGwgj&$jg`fx?$A5-x*FkAD2UV1MEk)Kn<>FnhVwSGIPy2u9z8B0f}N&B=Y!&_FTq% zJom_;q|CO2;KAIoS?8W8>62Meo0muzK+MhpbtZJH+2*{meop!W-l^ShFRLN8R+vYj zfetB9XN8ox@$s{BanZ5%cm2}S=i5JGmB0NELV-N~qu&{e)ZrU{Rk`L?obk*gyxu?s zTlu{_fz8bt&50TyWWti<@v8aP`lkuTNr7t7h98zV6iM@FPMqEvfC?`R2cH?AtgXV< zE9s*n$`)O)Vj?|&kBPO`!>ghM+43(o1we_P5{YQ>&aD_iAxZYUTmNMQg%}|YCppDv zj5s^I)CO9k^#a*lRZ{ZL8(|#Zc{T8*`i)bzqs=xPyqY^s-Cq2pDx;%S!WS(NT;}4R zR|OtYwy+78na&UM;j?7P4?J&lZsAzYStyWOSHh}O)-xaJ=g7`-tH1(9sy{jozi<|C zc0;XYJ$aP&<^1)b<(}#b55e>4aqpZ+R+UM^c#H%^JbmP}7u37qhlSePo7^3dSHYn$ zFZyRty&^e~;#3r^OKM+x&L~OV_ax2!_S$wrW?2%w`5$%u80L5C@!T*P)N-?qt@ zia&Ir#Y6j{Qg2_Od8TnK~~Vd{T5?Co3WpMTkpQxqSGOSYF(mOvBpbx=@&R z{{*BzT_Ppzrhn1Z-6(jwIKd~R1dpEoxXxiY)cV}1UNFP5>JZBH z;bo;KMyHs`6eUgRB)qgCpFTIb|6X}O!kUuX%d9jPziF0j=a!}07WSRgxYa{`6XS!1 zQ_%Y%frY4=)SNalR*RNC7@1ixJlc! zwkI7Y5~%t_Vd5jm`5tZr zs^Lrz@l0VWH-jDliSAwHlj>?szkCLM>m2=KsZMrL=|?>wbt77jEbT<@LT|c;{hp;Z zy3aw=#VveeOq!$ZR`Cg{WSW0UQ=! z1~iFZkBk-TAR~TkpKGAEKc~dXfx^HAxDLM(rv-hyLZTJ0izGJMQXnO46s~2=n<(+Q zQiMRDg4u8Pu45y$AMSU7GeQr(Q7v78b`XI{_)mFT{FeeFAdWI)Nl#q7`A=N95@_DQ z4v-tKXFESmUYd)Bl^FOrTtpF2wbB*U6_kISMr(#FsIZzrBhr6M4`mMQbtB=UEPu`k zgzUer%xe3)pd6hrgBHvI`qQfmgqYH5HKvMRkfAQQBL=zi z1_Gjkp11E5a$Ef#RbYlr)OTNO?Ugg%OO6@38P;eFz&1NY*Wm$Pg4zSpz%pRPE(1Np ziq5NhdP=sFt)^KsWH>uEUV?t{>6c`V{|7?lKU^&TAA#@xMV$N}vIeD|wfD`_0FnU^pn(^H;eWW`KNO$;TgHa*{}wu$qY7z{yoRYuC=FPRWvD=YHFyISZUk-} zZGZtG3tbSZe-{FKGL9jgN;AcjIBzf&W<8Net0AB2qW)k*WBne_>1}YeK1#P2j2234 zv4jF9w3rB)VVWr^wc9z>tg~Od_(VOG*`EvdD_4oJH@AuOW4GAr@UcFkav5{#J z9(UcI;XVE!S)iAb&$UlXfno;MN4iw?j}-uXnSC(PKV@M!(Z+I>!JC<+_5Q9$)ENby zS1|euPZplm4CUJ;j~Q%Wy=7|1$4^?`{VG$?2x#|nfgArV>rsg)1hfa{J7Xa6WxO3? z=E#S9Mwcj5fqS0>&6)%&DWt>wz1Ay5GAUIBiJmao@R;2O#)VtF-B^NTpms)~I&D%S zP)E2-ByrN=IEpAEGc{=us-qIcEIP_bLL?8j0}2}q!Fl; zXyaTehXVC{v?z%hZ4%XlOmD`HTYE9G4mHBHc<_=A!wyZE(b2?WO&YRCR;=T!EH;bg zM2--tM#?yiq8;V-dWjy&C%vnk<>cx0pwCy5Qe%=)3ui1Q`2r$Q!x#b0^4QO^0ZOP9 zN{ZCs8PTu}99K9m z3K^uAVCKP?laIW*NQSoxb^gmsZtF?-wfi12{o3Y@<`wr)T{nfu95CJ)+xd$d(wlx; zIUnY7r^1{PxuR*BnudKNE{!UfqKzjdPh#8_8>0Kes5I^}vGBv*fkNnYeKWWE5OaLX zFM#|fK3qLhDPLw}ueB%NUuOh5i$0O9v)pvN5jr025t`P9 zl%m?&n||k0u&k-6*MRhVepHkTsF-|$Q&k3)g1OOTXLOuk8_(9;<83h*DHYw?Dn)bQk044Z0zG{K!4^ zU1r=Ul+oH!Zl_Pc;V&K!9bKrZZPhnWiuav`CY8iQhD+Vp!s6xyoJ)5weMJG2C5&(S z5qHaxP+8SqHof?RBcoK!sUmjfNUa3WnV+(KJAB*XCIi|45RW}s%_WW55S4qptoqlo zx_sb|Z)wBK_Wr2&o@&|d%K1~BW#>N>(I=`abBzCo*0qDKi{hOGw|0({s(7`M=*hUw z_+jdiVKhSAy9RdJ?$&hRj6C^`(W^I3>Tu|~?mHp27}S~0rg-WcR$j6ls~7q)A5Hb} zFWr5c#F>itfkS3aa(pXv)SNrS5P2q{2fP@9yjkDd^Xd96pHzwB!RI=ZL&rT{mg4}y zZ&h=bmLI{%haM?^;1nT#ZAjf=#LyScwWkJEE+v1a4-WD%D1RoLz)L@UICfCd_t3iX z?F7r{lz5>1Ay9Wk+eR`r+~aQ|SBfW(I8n->Q_4n6YUuE;gNb9Of zs({(6Y6oOqpV*y%$5nz8OHV<`Z^gZb&aHqC)dQmFgXGSx5M~VH0ecKsBM48-Z#X@b zj3jpyX?DIBj;mNM_xf_hi}g3X`u4 z61@T42}T=CJ}U$llrKHJTenw_i#x6Qki+^o;UBZ~CiJs8gmGOVykGmv(9^cI_yH|^> zjlvNdPdN$nyrCV{xVe7q*q)u6O7VahtLlag7!`L{u+!GF(-*^Nz@U*WdXx9KX-den z6}@f4jKw~tKQ|4nDKszVKAMnNzpIy!U9LJ`H`PsJ`$N}sv2{?&YOA4xL{Vr?mNRTs zs5P#d+O}RfYeAn_z`#4DBVpinsm$f#LTn-iBta%g1AgCt(Ot|l2-G)N(cQD9tNk-y zbGw%4YYxUz4M#um>j|HQ zO95JukbMhD+|~$XAUn0tnMUmPSWX$+X2#;np^_1~JAxQK^pwIT)rZq*3ft^k9i(9e z)HVP#Bfz(KpCHCe8GZc`ZE9G8?`-9{yzr= z_W$s2{C`X$|79WnH;L@&BPi%g=yML-^#cMRSOSCo4=4YJAmx8Mg0QhN|JUr3;{ooC ztk&8a!t)ZsQ3zo|7?98iBGgzAYE5AvyMYjXjv>^*Sld1++e*_>u6bRqm3nO3G0`3b zTN(sODj;#jT2Dz+Y&RIRR z2xuc5!B$>kua7;>3xmn_EyzLMPc8-?Vh;yfZYOLSgVWYna~yv{+|zFy`iGi6L_e1Hjl&ZD@qc5v}A0UJ}L@bI+3Ii_FiK{^4xGCy242}tBHeGI5nl2 zSCJ@_(`R&=^3V7ujiQie>e;+6vyvM+Jumzmt!3z8Ff&B^* zSd+0N?W-yE7Lh(ddxtjAG_)53y28V0Zc(2i#XQ8Bf?Q0&hB#op@X_ANLTcM^5UbQKIz`+ zF7r0quRj>z^B+S(!fre|Big~*>AO#bgA#|u(BM2{@yFld(qetb=r`;&nH~JpyG|c5 zHjPEe*z z#i4wN;|#!c7_lS)JLePo>DOb>jQVBy?eH&Fj}gyk4vNX7+p%E5>H{$4;L_+bYE8S3 z|I47S5<-NZ^5TG58N!w5R}Ij*R$9E`>!u(7m8?+xoLl!9_>q@N+oT;nq8b04xOVKY zkhy_*@k4r8x4H`tO7rO8n55xIYoVso$V|yq%-WDk3m8+qXw;-JC7_RjgfcKVs5~+* zNu=5FZ+7k3_PON3lj6rbOcK||@l~Er+$KOXGm9F5-Jqp~g{5V&Q;K!dk4|g4grn;i zGHKnP0%|YmN^vKBsNzgnsE0x`8I_UMD9nO?)DN%`vpdV}NGkycd>(4)?4|!Z``4~t zmZeJu0D6zsM16qTYN=HSr`k5zr-TReb)O_V#EJBq^{t|+R!rd9UQU&{!x=Yt&~Z2m zq(;V#8r)2KpQ(<2DHDljXzvD%i98UypQu1BxCe^TUn$6Hf)+3H3?AQijq74w7{`UaD&aR?BMu$1G0@p!lR_mxwF$-%ph0?n0XTAQ z?>wDv1^!on!+exBbxaAmY)0D+qFN=#q(`=~R7CS|Y)<*t*Y6dw9o*|=5GTFuJmwBH z9M&wUV(iY8n2aNw;Iu2W<1K^?=WF|XRI$rpFRyRx@g|IKk6DTgDGX}wh{M%>8x_z*iKj4s;2E&0R&~f27Jz=Gao(nxF1FVW6fBs#72_IID@6s zx(JPhtOin9;yjqDlRlikynFk5Rr;!(^nhZYO43_FnX5GuHZbd-xigAECB5Ef@HsA_i=Xj4)ovs*>d4sTZ<77w|9Po zWn$2`3eI1B;4|8KT0bwM!CR0uv^J~10Q>vsHy$`(U8)d2*^Ck^StZhhXhb=PE?8B$ z!3;_PWP0OxFsT8h*{*>jmQL-|U7!4c#Pb*$L{&2N!vkaz!73XI8xD{e3I1pUn?(;v zs!1zHV8g!+J(Udwf;;Kj>E$ur*I!c{3gHW{dIo%egUM6gjyfau#jGL~SU+Dz2J5ySm8pDzPDN+0$N zI8Aue=WQ!#ev58c`=@uBlPY2o`Ec(Im>6<}`|-u_?9Y;L_i;(;?!o{%lRfucj~5Tp zL!Jn~A0;=%;kqQbE&=E^fhE3ppzk;$K}>!*w-56KC(YX?(?u&&N6V z!C#ktBfHT0He0DG5#hKFQIT5(xBgU5JD*+ji)nMcP3q%)TI>`EG}?Ew*=uNdy~0Bg ztxfp1KZi6r}02nFt-BsY9O9V&DhA^J7&M3)_XtDq2%Y)>DLu3dies zJS}0u=@0nw=?`W3`FiON1QchC$c@-QzsGLfC})OPT^#H0h-s2urPO?HqIb>E8AAgz zVoqGRt=^rKVDqj*)Tevt@UUD#Nd;%3T)tiix`svWluEk*4-Ts#6Z&;e?DNV0Db`-Rlio3**VYjBW;>l-yZyY5J~LaFUB zYgKTnx;$EI*X9z%6RXygYuV{THH-@mCHZnKN3LtK1J1~9{B8C;0$Mb{HT=P71(i@G zQs~OhUjk$xJQ*&G!*Q9{MJf)WGuZ(qQsT=9cA+yS^F0}|a#GItmYsG0^~NQ>DtuU5{4x1MMf@P%zt5t{SoOEFza$6TiW z_1x_mx2fbGwud{OU1@Ie7YZg|Qzgj&vw+I&8&0lRv_2|_Cr-Ci@h7S1jAkAPpaEEXXHaX zkrfp{)EL5EL2^W(CRa^qdO|Avfn5H;NBd!i$A$@B%*DNWBR6JZ?Wn^M2!CH~BAbNP zHNwqLWMaz1dO%t%(1d`C+P)&lXE{W!54I^#rt~O9^@sWyi#s??w>2y;d3rz=4(V>p zUJW?!RTQs9S-t5Bho^P zv&XQOUx|#E4xdN$dact3els`K=Y4ajB@gC&Edl-+s}O;KG1FVml2ipUj1;2exFOQK ztvy86fSCo@Q^R{~gk=D(7=FGI|2)<>iYellsphb13@X(G3>;FnJI!HHHeXiJ^hINH z62u$a%Y^|S%Gu#l23e2dzW5-1*f?w|k~FC0XcR@D7NvT`j#mB?)1Ph*19r1t%IfZ6 z6Wx)VKrigsxgughvFy~y+CAS$<3EVB0iu~QMWe+m0Wb;ERUG;YmFDx6NzNVSa9&AD z){CdYM4rcIk<=jj`){OQKt1j9098pF`WoI3i{sZ0@gNg3!6=sbr4KLL;k$^{-bHcyOnfT^%?L(xJ`ev${} z1c3$#Nxd|Be{ct%H4@Z$=B=!NqL%BbrlwKOM(l2;ZA;?^sJbLEkkn{XA3?>%&?O`86lkgs4O*aRidBZlOx=^)#qm))#$9>* zGb#49_tPmxl!%p%Fp#ivEct{bv6UZhoTH^j1;%1POFd4c0WLJyj}}#RxX>{Qhs2s= zsSaH6Py$&^;U{Y7j4~kydE|-Um%ngR1B`DuOGzm)H7OIpLvmd-C%6V$RYX%jeS?w@ zKGu-GaO*}$=+L8>MdL}$9b*3Y#dZczH*|R;jO)&GY+l*t7Z5K7yo{RAnOmD=is+Md zJ8Muk82W{J`s2!qy|_M_+{u4pjalk zG8ONv^lH!>U)?&J@oOtv#ll6flw*8lxuQp9FsVc{knS?^d?!13TXWEnO zg!7xhZQm}ZP~h)^sJgB zwhoADMH3*GRc@)6Rfiovq`8OKK&eKrFY!mdqm5K6dSDd7KkbQ2KXrA?fWHR`C<{-X z7fKx4gc(WgURj8T{Fkd;Q&~K<5!~%pCe<|qr`A6aNZJPq|z{oRD#I{Q%m@~ zQ4&-kU(#&oW-q~=CLulP+@lK=A&HsMM)u`IR*P-WRAN`ybJ$i5+@ZM$jD~e5$H%Kj zhs2m`*bdTTq;LJ1IWY}v2jT6s?04l$M66;#lP&IJw;#;+v>*92f7iWi^q|ob0gXxA zPnXy5w1NK%A3?X~VIKrPcmL4a9^36|pD7a@C~UK)AZoyT+5AM_C@e6V<+_3Q4l!BS zW`of=ba^)`^(T{B_843roC|HIH+VdTU*1TVJSQV0?BJcfkG5=gRCQqL%6nbR4(?(E znXxT;d9U{-+nta;a{FY+@0k&3x13i`ba!w?@2<~uGs6t@XOD=LiC#XB+)i%u9!_If z;sySy=_(kVuDY9c3*6^9r{;Z1y6bT$d0{j9t++hcCis%LV4SLn@Mmb0$qYETD3m9T zMTaIr?>XXJ2W!qB-ZlkaJbj3!#2Pg?C*nd|(30>BR_ZDEPd`N(nV8)~G&$|=KNaIE z4x*4AGo1CsUmpWMn=t%J&|-rVpSP%`I^B`?v>CNH7e{cQA7Loe%alY9MMj6Qj=7Rqu`#eBs1EK4|HfLCT@BC2PV+BKe{O+TAqTfC?f*D)V zBn~m#I%Y2IomVd9KXfKF6x-*&;3u1{twCY2y97y^#1YB5EiuO@;?zsemhnZfJ`-WT z3t{I_O2cft8gYA%il5wCbhqu@AJi`7zvA9L{z%GH-RGDGG7Vcve`VTOJ~KnU@-o4k z{Q8tA&v^DiJk4UJ{_Quo!`_pcmZX~S*KymJe4L56RX-)ue2efmc~uwuku zE6IW;N=Q0Lt~}2G6$0E+nJ=p@X?b6)Y0)FanJ6}s(0|3O=h+gywpOvKKG{2lz21V! zUT-ZsqoaIZo{E33&TSYQepAZ&exz5x1ey%CN1N5J-&b}p*1-~&i6rHO_BDy>zwP5i zn6z?c$`6cZ26Td-ayW-md0c=d9dbxNnaVO4nu9+)9dQ` zbGifh0>%+qI?1I0G@>RI?L_PKP?;bvSY~bRkv;vZo4tI?T_VqXNi6l!1ELY>My1-l zp}z~x+U`|E@1x_0I?dC0NBl)hiuOb3?pyCl%*mYvD_e|LNYkk`PgK2}qn0?-E#t2u zFq%xcU(j&i%n}4_s9zDeYP23f z8}>`yOI~?Dtm)j)+R*C1MJ7ywI zUjsQ7Ego|R?YuXhv~F9M6i`M_m4@^|8E?thq?rq#a7E0-NLB;k3^6H2|3p%4lB(p8 z?#Y2TUI}Z(AKPM9Kk~o@XYD0?2BGtJ!b7bK=>_m%!m8sA6T(rvoZA)BZGZVifL*WypRk-QAU|b{IQ|a zt}HCrP?b!d3InG60WE3DxC_9we4IKOUSkq%*veNFwQ7=3z=U2U{P8tWK@EAF@(+^yi&(}sTYvCWhR>4i`ij_~qS}dIXODflNOjiwBq2ZLo09eI z9`CE}^Y-_6cHxq@{XGx+*`V)giT%^wyN(J&nq-^Kv=$3aXM; zK_o>ErOLk!=nAUJWJF5EWV?hSDpxY@59rEz9EFUk;{Yvw8ZV zE_bS0pbNDvT7Qo>Is=cUl(#iM*ya2MOjcjc7yoaRa=mYR%+-cJ8T~ntA)e1gyYep@ONeMrAQC1yf?j^J;iRHdrM+;kN?sR{VS9uxv&_{j2n+9 z3J|Ln&h>eWE(Sub++AU~1AKfhD&}v`ne8inHobCTuPeL&acy{2&uLoGpj-4AFq#`@ z?u%c_fF#WqH3CYSub_FV;ozUG`s`b^?k@Cx3ZqG-Xl>}TLxN;E@~`iHroLHo?q~r% z!PDiH_ppF?=@WF)qlUWQc6)L?t?TedNOI^r(|(Da5;`dSRZXI%vlh^FZ~?6$WX#6M z8x9e;UE)hv*m{Cm!4VIdGi`phX)7&fde`CiIl9%Hj86k zRYg?{Lbh}+2*f4lVsU$Jdd4prk^#@FeYMqS{k@B{$<@ZD4*rZ_)o50&N=lQwSg*q^ zENkYBrxYqD`r|HAsoo9l~{~Kj%SCr55FCScT=dtFph;! z*EB09j=rO;D%!X&!_xE%OL$#SS5J%70#Osn55t+QK|nc!f-Q0)3j8C*qMH8}pHee- z!iJyzN_v7kYL6G+0|Ej|KSb(e#1ZiSjl-lPDtL|5K=4!!KhsZ_9OQ!Db9y1RMq{Zh?Id%-!9l zy(T@ginK2ZHOh7o3ndtag(??s_A65GD*j=XDRiRV*fp8eQEC>=iWln;DRe@p#juK` zn9ZI+R7!QyQVcTp3ZSB**vgb#KG9O$mn{vNErQs!Qjzd$7b+2onk$Vi0fJO|I^Kno z{GL&4JNYWIS*1#t&6+C8hK7UM@^iz*67V^Ei_Ah_qkEX27sNqjE;}r{D^-D|?{#+? zomBP};)a%mR6=EHAy&vM3k_g#zroTGefP^CmPD*6j%7R*n9DMgS((B!MQ939701c= zr~*CWpo)acrc}}Y$RMi7$kJ9N2j$GPr zxx)2@kTH)>&xl-rJ|tJ^`L6pM5X=RO?A7Bd=2h3itUZ!oTy0qFSTOo6rptfhmamT8 z5|3e72gFU!>nHeQM+e^OcRz_Dfz-pB>5&poA+~2`-o8+L_Shw9n48rp<@H4T&RA0e@l%*&N_bL~Az@ z-2Sx;b@7o43Q_=s5+n5n21Cm(`;OrtKAFvgtzg8$^R#F*Z_~l>X+Q~k6y|qbVr^k< zSph%xyW&oB$p(J9!fU@5Q5?4&w1P#XU~v_8g=zc8Y5Ri1m#`dedObY_*use!4OZ{r z!kcLlkb%NX#2Y4PgP<4GYkd%5D{W+M4MtyCXPbx|)#wDn=FZk+d+aR<@Kq9@;1C=i zR$}v>F4RcNzdSp(=@)0xa8BP+QH`Z&1JhX^5$4+Vk4Tyq#Oc`WJ&cOtY{>G2V@D8R zr4x}A19%63iQ6L(+s)aa>)NF``vav01Sy9jHIz>eZk&r5nc_*!NUGmVAfm(A<#a1R z{DU1z(4ijc8}j?dPfu^B=c!l3TiI7!yyua%hX$^ejRHaD`S!+=KwILM%JKB@eL$DaV|ML=%b+TDFMbV3s$z)0@8z1!ly&LVw z)pR1qmtls^`^!|SauJ`CfP)zm@Lmn>r@T9mXV;^!*zb}f@z)w{$a35!tRKyN?P{57 zGiy2SGlM>agt4^}COgY^dN(6V)mAG3XA+T}*#6XmgWAZ~l43wM%*oc~YQJkG0Ka<- zv94`999fa~8OY6^|Ghf&HJDw!QCf9bWnE*9&aR4A({4XE4SD24XovB!{Vrc;(FvQc zmll<(KM!;Q-%}`aK|E|0|5NIH`4J=`o;{T&y>qtt!<2)__41ETP)e{b62zhY!)u#h z?A9K&Y7$zjD>YGN-L=yHRoOQNNA|7n#&#yQZ5tC!Y-eJlV`JjZ#L2`?Cbn(cwr%_7 z-cP@C&-ve5=i9E`U2pBGUcL7Fyz2q2Goso4`Jp)C7gM^K)*TvwC$0MC+#yZH0M;5eTLpDHOBp`Z&7Uq{JiqBi>3a&k zg(X0Yh>&8)8Z)SfJPfg_X^L^-Jf6(LeKG7vqeIkJgS%Ns$h(3uj~ zp%M*n44!%>Vr^^!ao!tO47}1VfGVf!^^*qIwQ;Qx>1ZvG8VxMvas&lGJpQ7fZJrQ6{CL!P=l%zvo<8l=Wkk9EYYNhmHRgVhCIHkGRr&NGlc;$iQrWR$EHbe;=G0OGvV-47fxVlCz<* z)x?aC>X6e{-ydv-h7O{6*#+ul+ zA@oHUd3*%l3IZAt@!^Hn4aZeupI!d8QU9U@E@tr?Y=5j}W|wXwOYCeZ=3|3B=RvNC z<878itfjjGUChJ1w_S}QZYS%Q;rRIpNn;+wualW#GOW47op2sri>YA$4 z8ngQ>>r?|7m6Ef*15wfSOQ0#aE4bL_deOJR%&Ow9chdsNz zc7<*fhtD+dIGopKQFHRW_6G3kF+=HIB{HC8LP+{Bg+ykurN?g`34GTNT@@4awZEbn z5;xT+p0d%mcRJ#m?y$0^k~TH+pZov}`Sr%rS(?=}e=`dZScy%j6_6`t0RXI67lGUH zAG&PR$bq-pqoEW88pQ9vRg=N^S6j!^9YHLk#FENZsUAx;^`=YP-ar=p>r9F}V9x33 z(eGQ3u(+-^r|X;8$r?c}M})m@Y8`KG_UBuGaZ-#?vi&A!5ZD=Wse^|pW6fU8yRA{D z^CNEU_KrZmdolvvWNZ4h35wk}vCJuF+KFP?FUw2l=^X9WqImOb&!yZV{_q~h|lZEU{b@D$JE zJkwb(fOY22lIxJJ)6OnKmMn2)?DOsL$9}S^2&d5p*e;h@RKKiv@Mn#ub~VM4tS*Q*lA3^sr`_ekWeCHCZ6su->`xm9}QnYL`;PYPuOf z=VpE+=zttAqKn#xv7`lTwg2QqtT5f^E6^3YuU;~(g9KAa|A0kv;fgpjBZNg_P8G;Q zD$wZbAFLS2;v(k$(%S66`CF}U^Elt^c|WUHdtAn-r;M)mopLFpcZT)NIfH=}<~`#>hyBUb@**D@n_q8xPY8o+Et3RwD^3WYP%qMaV78MSX z&qRy%F#07$J!Cx8P7O^pyu+bbc;M%+udODjPW#F)(~2A)T0Y-faK0ZLUQglycK9$=xitppoyXYW!OezY zV{2VCd8~zW&HcJ!WNM^1DS1@ZS%!n?hiHQ$Kx1oZwc--0$rzFH|B{>FHDUKcyU?_C z(6r;SCM;g%*UVZNeBRG^R-7A!JoA5C~obd8lv6jl)!C2Yy%zRVx(+2ugh@?3QyF;hrLoj$sIFtiGPVKkU+QaA{I z5j)Z$4r1bfk)8Ynoq9|r5j{lJGP}iwR-^iPtMiuENWw?I4_E_3w^lxd%pVgxY)%s% zo0ndj08(NbTUdYLt()zYei^J+dA;q+$|#l2Lly=_oY*3`Y0bezhux+|P1qfjZWb_k zyIz{0wbcelG|)^YRBiHx%0;75nfp?(X&nfXCBYfb3m33r6aG)@Im=-&>p>!lc$jKL z65@0q3j`@8N+(JsL|>;uBeg6+$sMA~z;Ui~2+jo|VJ&9qE@EmCP%6VHBy>=g>(ry7 zXH{<-tA<0Od~raGMs%m=o9_jjpPTqE?a4<8nu4O;Z#g4FC~8fj1QjsK6)NOb3Ap8? zqIvQWiygZaFux;aXz#4#A_6`asMSVL%2{P&v9Su~sGVb_`!6*r;gp|>cAbA?DCs0P zXy+TPum_8ls?w6Wv^;st7M`M!t<9$bkGJ1AR*p;j!|P_tALVYfa8TH7gD-dTK}!=` zPUBpD+t*8=7R&{U71vz=IPtR=`a+6eK4Q>!-PsFJE>UKIUvI!Z1S3k8@MO6$68^B} zAgpbwW*JCPn1v3kD zdbnIl;%&Q(xx5n2jeR>hI$9#VkG}{JS6+DIVtZ`(Sw*BS7?L15u_S#tz8tR<0T&r& zlZwyhNI4xlb_bTg&zPWQ#=3Ov(oU^6cz}~~F0mVBTeD==03|?WjkmnODS+2~_L>G` z7XL~KB)u!6ZHO&5&5Y+5x1kHf>f*(1jqwY@BQ2iBHwhuEy*0;_s#hG=b-`ONF`N~3 z{f1>mq8R3-?iWsFm(u4sxKPNcUjN7CTlP96WPx5$wj+}c<>yJ86>MPA!A|UH;YiWB zpvHTqn{z@#hQyzwl5N=!R1X>O)q%+;_;Ta;9t7Y@8d`POc-(6!p`hmj69IEEa61X4 z-3%oKQh{F42-VRtJJ)b8B5MO_!$eV>MmBM0W&T5?v zZI2xj)h|4)`M%!%Mw~Ml33KxS=1*Nl^@KC-h^|Jz2oEPfi#Uva&=gY=a#=_xzg+J? z3GT8L)Y=CRXe*-iYh0-GC=#nOW;(b<)zbGVi8d{8TIb`eCQ-xC8;8 z=+;U%eZE#&m}?Q3Fta5a*MnkgYaYwzCDHbso9J1e9o_~K+BmrxG9sG*F>E>{Qx%Ul zypyTC?&>+kpurxqI9?md`)YkNK&z>ejS{_!cQi8ffY0-H{MU59+MPwf?Lh_Td& z83(X$`pe`*>p@oT6*LHMa6{bXnULBs|+pp;m`|j1NhLs&=fIyLrLJp^K-Zs%k4MPO*x1k|IH;eo7@bd1a-A?ckyIMH!1dWkwGk_Qu+@71M)4QV_%qHaEM{Qt;!yX`xGm3~7e<0bVgOQE< zTBNh*X0K9oBB@naKV`CDLchrYQ=3g|O-KY?sC%+$)CH(+L0;+}4c3ia= zmq~fKVW2ZTQB&nQ+i`+vlAkw@aG3`A$PhDmm=6bboQ3rfbWX-)P8C856J&Ee^^OQW z9Cq3VrS4YMP~<*o=WPFZy&`ktMb-X(z@)_^+3x*_WmVAR?pFN6RahhA4`v$9!SlIt z6=Y!{vhY$TUC&s1K9(@RT%#~%G++yqLKB2(==*YOUQFuT!QN>j3pY+DnyH=f?On}L z%^4>|fASTTN_}8!i0OB{Tj-sOU-AZ*oVeCfQ(aE6-*7#8$hKM62`wq=L&t!7KWrP( zHomLZA;+zlUicg!C{4C8{7Heg%P$|wxwZ!oH(oy(PW`sWa>XAwOYtqj83B(z#63#v zgR19cQ7X5Svaz(Zu&B06lsSAx*X*)TQ?-UlTk_1b_d(0*mSAjG(QWpNq2u?m01#Bp z_J(Sw??uKS@-wFRsf(Rv`$l460ROpsZewjwJ={gu3+=}p)wsiBW6c2v`2)%r~+-EXXK6~7;F>wU`LB7TR z`1mE#z~6*hlWO~0j?>I!k9hyAel`B~O!+qcW{y?PR$guWI*u$Wck!0L@4M1Z_96U* zWW$#E=zy|@t76D+m?Qg}a;cu4#Ffg=m4xdx9-8=Lvc4TB f zZ|`yj1mmYO7MLW*AEYdv&Vd8GDfmH>~xe+k`sX zOO!TrZC1v^FNxeYeY?#`xox~ro^e|02TE1w7P0TOEcciDn-ScNL@<0e6*zhD*)%Ww zXIxG8{|?tk-^oE6z!RC79z{m-!EW1Ot)s5u=gk$VZkNMCN3+-0=Wv1Cr$P>vb&&2M z?Pi0{aoJ-Hfo8}l<$_=5jQcy<@P%=i^PhyebsroQNBP_%L@g5&(h?jXmkcA=ArLef zUaX1xt|Jq=>U4DYIObetv)=M|vBSeKoVophP8#8vDu__vUu57agAonwpUd$7TZdC3 z`4Q22fLJO|8)||J-yQas^27d_g`3n?B(mejhC#^&Ku*K1arybpMQGF#o*8rBMY^Bv zj{ilwj|#_2>8Sc?2D1LvN>7yb&v~Kg;Y8*!=xd2(3ECe`O*{7Ca?YM#;Y+v2Mg3mUefB2T`$@%)) zMyzSuF)xXyk63BrNz)=H<}4W>%vn2Dz25i|k*&w~>9N8K*6Z$}E;_hDsU73htXOQp zk5O~7q`;ygCu{4$2jHSH_0Z+o^Ebri*G`{vaO^`XvN(6KUvcU;u+3krsL|9?Zp!^U za(6esiq=D%aAwdU!pX?`6YzEm?sl{XzkDl)U=nc1>~zy6`JqiXp21-f{pi5Rb8psH zz9E!#bjsDHw`c&a?{`;zwz>-akXVjMThR%2{QOh{pgVc{KH3jo^Yao|I5-*XIjLEh;Ecqo8e81Sv^Gz7qe?H%p+2Mvy7%h&P9S@`s zslnZrjl*?$UfW%zidt5Ru5it4d8kJquV}3Pw(7_1=ITEBfP30Q>-b=nn5&{zxY`l3 z*;?yxFoQ_+48_6?3Jqkfjo>+ERYUKjFTo~^`?LNlcI#Q)xR3F(7u(Gj@^CD`3g!be zAXt}$E6i@b1}dokugGZVm0&{@k4(kjq1f}@#DQ%gYt*V9$6li3h})Gez}&I6(Wp)Z zkfcx%?37yzUiqDuGsVNnPasQA!#+z&Lp9k|*0-!`Om~!hovDc4rizX@e9?6aRb0!J z2BY34xGk@$B?s6H=99bN*IN9aWT?y(@s zIH0uJ6B%MNe|a;U-ZJ)b#_KlZYKg0k8Hs(f!ihGpo0OYNV8cjf%;D5viIPBo^({6j z5ip3Ebe^E}Q!Yfi<5tZ3Hwnq!)k5v0e1MDxD)?-(fzoLu=jdViFTWJNkx&4>ZZ3%d zT0}qRk}GSXp%DI6dgZij*8|Q_7xMF;S@+*W4f<7G$Y+LdP1>0c8z;Wd3KGhvLg;zi z&3*R_9|>-m_gNu$6RV@m{rgqfZAIx~4-Fz#!!YVQ1`Une_{EGgRfPQoxoBa3+kG{Z zEJgIXGvPxYr0|LFnSrNFXdQQ0D)k z_CHP(|65$j!~QRcLV7AbHt43dq*U z4eh7UrEd2ZFC6NmrE|}-Gul6&X_~*- zcC7wGRWc-t|529OgD6DD2mN4+a-L`~iyS@AMr%dU@U~ zA-+#@EF>a!(tW-fbM-?Vw|^UlVIgT9-@mSNsS2aV3KxT{8sA;sD0$BK!X+)>IU=0n z$OFQ^a)cY(53?QK}Q0S=&k>;r-9qljP} zj?U>wSL~2+2;xq{z?rVr&Q0VNO$iWzYLxLSc6(mXElXT6NRsN&bIdqYtWohbLzfV@ zK_aAUyYr0N(#pllv*p$4<<-q&G*&Q}TKsGa5|r@U3du*(63{)0zG!^m9x2)AZ0-J_ zWM?ZgVTo7|XS!4gp&0=vw=|N}xOU5Lj*Xq2bS)tHb1X9cBkI1MW@uD&^%H^`ce03>=8?#B;G=)k$4WgQ*%UU1ju{ql;qQ#%@U1G5;GJF2xFvc@36lNZIzyuUX*os+I{KG`d;Iy9NrFs z$TO6uZIozt6o%#tgvvib7W!iTgR43BKTi1n=Y58S`=3$0U-DD~3lL)@Xi8`h6jm_! z|FFaSW2*Mw?lA1^|3VwkQxdU87x|X?V*zuuT04idA#&16e0m*G4YhqfMp_1!L`$U2 zwl#j8Vl{=CZFD7YOUYuf2k~w?J_cxbLDr?ZQ;creElVhBpL(O9o7ttHrnKs- zGV=AFwF6$Oq)23tGV>>qry-YJ0?-R@IZ#=URR>)y;BuK5^Vk*kPK%ZiX~{FPO%vIY4; zBPM-+(NZD(Jv4mk{2T5ur6lQbDv#H`$ZRcM>DLqJxy_+N5-<{Hk3^#%zwtE<_a|m$ zgO%wfD`y8uW9zWxO{#Q0AY+9rU5#B$`tl2!3ffCg#*GHY9wj>5t{&ScEN(PrcsKiQ=d&utp?!%0YmKH2JQhuf}9-G_6r0SkmsgADs?!Q!#i87ZDwm<~)cvz1zMe8Z<+Az_@cHm@jsPD96xO6)V7q?Z zQHhUNZ;Ov1rpA0*JxN&YPZ57Ro61|C1G6H*fu#PW2oc`XOKrzpPu=sa9}L+w$oyS^ zHev%CgD~Q^dr?lHkl$T6tR9e<{+>wmF#LAwtf+r%eoaV{0d}1mz5ZKDx6$}Mz3$ul zq~GVdU78bwR@*sGn-PnCnI?Lp=l!S-Bn2Fa5R?9xUajZ73{%L&X|v8w6CZ0WJS3|U zNmiOxdha^wpd4_Q63%2yPw>%0goP_SDNF0mD8`;+yhpKa7?A3l!)*}QkFX>C-KODZ zZ&VFR6imky07)eHF|Iy zD_n-NOq(KcuOBwZqtr2lC09*cB~YwX^j6(&H6gj}CgWUZ7yWTyFl^t*u0nWdd^?Qo zSvVpz4NRl-i<)AATa>!{2&V-@{PAjlXoRhSY3_z?CV>`4uB>&~n5=>4qECpAC z%oAZk@(e>nBNjzW0Yr%nAuF4LgAfg3(Le=4JclI=zYrt9)&@@^03tQf92j?G>XO@& z%QLHO1`?$7_7YM;%^t?yT-6`SF3rp*K2%gM@SUSrJ@xKg!?-c+1+`ie(FEDp>=JgCQ|_{&uHs}rNv5|n&q!KW3r&PKnFA0HcxfB;!J z>6>8tyX~c+U%y1|Z`+i<71F98@;#A7!&VK;%vrxw1GyVolIacF>_{n>PmdH7ZAoZU zjAg)oG#ZkN;Z(v@jgRFlTOU=~7LOo@ki0-fG`X~eHsbMLyBMNidSl4nukAPO+wNn` z^X7q4n?Qo`A|#cQ`301H{a3Q*Z54UKtbFp-gA(dW1$;6uj4z4!n-dPESlH#!e`KD0 zkjD5J1)TLlPo|_w!Y17(E=G&b*8Yw8b*otg!Uz|7g)XOkc+ac4$}pG;O)ol;IUL6^ECFSP#&nc;_ShC6Z}_-) zWTBs+<@Zy=9!}K<81_Lvepd?Sl-K?fHOi!^iO@yTV!@x0wQ+pmjE zN}VoMaH?3E0vJDH0kt?$STJo7qeNAYQ;$lXA!U2PflU$sR%_Vsaoxp$`!L**-}(9d z$9Vy#ZVDThwjdQMh{jRXrB_MMS=DKtDP@LJWSrIIg0HYLDkjn6CxTHrZS_z+22-Uh z4%S-c3Q6W5K{f$+;8$kE2WBw{10TruBa}!8g%9$hu))gt%)erL=v!e5&8Mix^Sf_G zTQ8qNQ(K^5KqGT;zi(T!!!uP9kDCS+PpBRmP&UBJ^u^_7WPNR*%lQDR)L-ZyJz8Hm zR(6opljc%`6Cw-X*sfk{^;~IB31lo+w+6P%3omqX(g;?~gl8&=-=yv3BJ+x4+4)|9 zOO>as3g0Kh&*WDOM5PlI)dtAttQANmu+ySEq;jD!w2Zq<80%NfPtO`apWd#tl3|vg zvgJsnWePvDX^4%=TO8Y^qX@1Jvpx>U-;!KW5L@lOj0EOR>`%zs`f$dm=s8G?!eAX92FzKF(tmoFF8hPnPOAyrRC5#9%f=}iLtZ*6t?4R$K zRNG~~D5K-qKX>kibjT^tax^#SuRYqzeMlTtw4nx~1NZ{XKlv8e5Tz1Ck^ew>cr*T_ zsAk6YJ!@F-thO%hcpvXtwm%Q)cl>CZJmVBkoL|6=6xLHKRyN3OxQTKY&A(td?#j6X z?>UFvQY_JAoL&pWIuI)NeMh=jpfnFfjOa>I@*6i_VW(eEsi?#-gx!6s-9!@J^0r}v zCp1)6x?9rC`I*NGQY??@hrLKpgjW?uDZ>W3Jm-I`RPgz<+M!?5jLel6;;3BOy8}p) z^%jC)!K!to3x?2C9rfE(ku?d_Onq7gMU;KsOrXUBnt$lhn9%vSV&DmeojQMQf+3N&( z$Du_l(i{#lpH}o|(TWBO26XN0AxK6;!7X9iI{6%)b+Q!mfq-hRM`tWng+U{;ZA4{I z@yhWTo4t^>@yVgw7Robf2S~G+)az$)B89Y{>!@Z`T(y}BTw(K!9ebei`?v9n!d(e&1v8>R~9Us8^lgyIJ*)vl8sfe@wC1bbKRhSR_>RgyR|CbC; z7L&D{)30{K=7*nU8ylx?;2N9E%fGlD<42HkTUIoc!mhpU5zjXg(a|jfnPy~e)y;pO z%CswQB{WAi&c{l+B++_krj*1~g50#x8?#i%QwmUs>FLFfMmDqSN&b0Od)4TItoPl- zi2nrR-ai-cfiu?Dlg2Q3AIsgtpQI`pHfqo*_&whXY2yH?^_Ri~Ptn+PXjsw1qOy@? ze)^sS?ywd|y=ZE^^H^x@EiKW!%8)%@beqkl{V%Jvm7R57{^-NdYZ_?s&=m2ZI=9+= zlP3@%0ggc`-ps^oOS*X6oNPQ1a|hQ(tTS_GMKyFQ6x+%Y$0&7PpJeJm$NgY_P)PX^h9un+zw!xgr= zUdDK5F(b2W&8LO~ZGsm*3sxBa)$rYxH{xgQ?}LI8{E|!wtmA5MZsVuWDAczBlZOS4 zG*-*1P${&p*fT0I*0|bo{L(t=D9(l$Lg*vF8TsjT-k@s2&)l~shkIE3iFPC`L=*+; zD`k7n@M5&nU$T{J6oYCKNltM7cvYuLPR89WrbzYK+KPJMNV925)r`?UU!W!yJzrn& z$qKh*18B&jzNj-5Oni*t2U6eb?zK{NG3}Sv_iwjc!{zf>R8K}=Hga$Nd|zq_Ya90W zFX+ecY|^fc`{;tupU7sattYve1a2ijv#kQlvXU{ zD0WUJO(WD!FS{ecvmnI#KO>49 zU?lWsOaDzc+Ka6wEZ6N7iePA`;_)b(1hnyu}#18U|`6r5YA}>oPuy7&W>a>#RGJC2M&cxfd0pAVK-0zi42g>eqHNj&<`-nMX<_VKk7eNeTe zaga)E!$r&fTHVXW3aW-{ilCUv!HWf292AG&K$(q`kJ6jO5a&mM1bEea-r2{}8M)MWw%-!Ui zrj^}flB8B<3>pWMkpeSJxepAZ>Q*+AYzW-i%&KhhI=eQl!6BxKJJxY=ahcahBK&KZ zKHL=)EBk{o!*qltAIw+QWAcxgHg*BREG~JXWjM#Av%Rd7Sycg~Qyd$;7g4f0J;WLK zS4RyO$z?DWG?#oJOHI;+P(CZmnp1BY0DFkfmVHB5mJGwaYLHbGr`7%UAWDm}d}<7$ zs<$Rs6W=JHlr>!o#^$dXc~-Abz;}7w(YWjx=0f7!!5?UDT@(a=ufP~4_K_U@Yx`m+ znDg>Yi8(My;lj~*F0hmLv3EJuVdK7)_iKuvQA4>kB7qn6wEXgb!MAr;R79?zr*Qf^ z0-_ET&OID|(uFMvl|J~xd`zZDaZ4L6&VK_hJ?OTrp91oz3emIMBl*u1Cs%-g-%tEV zNj58DP8&~c{Oz_mP`ML}9U`v~d7yxdd0%FpiqYgOPc_!j-=#oa)3usbopjlOB=h+) zdU>R;9I=ztGVYFZeGTa1WrcAEfZL-GVl$>J7W8-uXV4nCHp4)4W=4wJ`iI76Lu#Uq zw~lyNLm%gwKd+;X_;AA+nTMjpVO?4ma+&xc-d{7^Ee0H zd2J%a@iUu&V7W+3uN7}vTMy8WiPQ3-sA?|0FTr;ADiyek` z*1|c6P^^U&FL9Jk#zI))LqA9CqTb3SYb1p)_o`u5-nHq@_MBZc^z+2;q}^Kw2bjn> znwNMh_B<((q&N76Y{i=^)43f{iNfCFU+(S~lH1N`+!Mx~<6brl-n#+TJz*}}@|5!$ zaBlD-`<*awp!U8qZ6m4kcu4WKp3uR09qD{NBp{R*UyI zh?ogN$Cn7PcYkOi zQ3M!dtd-vESJ~RcdM&0)nk;N-5HES?-Dht}?e2IhAa!W-yt?yX1-f1KzN?5XNo}U( zMZ~E)xPE$Iit~KgAu(JTbqba5oD{b2x?w?gcKSN;cAQ5Ct@4ytTmoLuNrheoGt?Rf z^GY4W-tX;s)4VrV@FIjXR{K|nvr`cy7j*`~4ZOKuJ0jyRF6>HlZa|PUs6K}1lj64t z?4>!di}$YrW)vJmC{{=@=Vus?!rPB29B$SWeSg0bpvpQxr3hx51=wAF6ooA6qm6GY zo`!E5JOoQu4OwNanL$^X@Xh3bC^Dy_%-sDpUNUiS=4$CE#%7ALds$R5>we_p-{Ohj z##JaiZ2MsIJ_c}YM*Nzp5Xj}m-B5)~(EkmoXwG3ga417YmJ}tu(?v#JcVVe!W$L+N zAooxg9uI*ei<^x3hq5ok0PGKqC7G)5xrvqYv68%(YHD9PHi{=$^dU+BGdjetJ<^yW zRtkv?338r41Igbv`{g*0ix})g2;|u0TWo;7RF7N2K>H_d?5e{gZ<@Q2UR^J zVu47$#Us`U6pR-|H?s^E$!2lKCQ&v_E)e|@nGnVlQ*l%-8BJeR1ylsu!5Bdu&q|}4 zXV3@$PY~;fV8-!ND1_07h<2i-EU`q?b6Ku%%Ait;z_N!;zmfJ}C1vZ?d%HvR>xeyx zLQQ%8>136VP6aCes*Jhy>iLc_K&hX#HsaR%*osB;$?YN3h~mfVty3QWxK(*VDt}tY z2?L0ua!j!$ZB%q(Zsv zRUWpi271f#V0isGGrJEDuXQV|#C4=f;P;GqmWt?MC!|ZTl!LlE@nqyCkA0CCGH3;H%jbIiQiN*hZw^7I{|$^mcTf# zz9;QA`@DUed=mSwza+lgJ;KE9x$-s4tVI^2P64MR=w!aysR<}b)6i@yWVXjj*xBHI z_FI?U)DCjt?apk@q8G$uLrY46li28a#bSsA=IpeSGc$)g)~!5!-KDe0l@Gptk+|Ms zYIi1g1wyYk&fQaCrk(aD=5M06oSjwbTyJ2O(947#NEn!eghjomPP{?APc>neuj3n6 zOxn8K!SOQkIcXwFt(l^RQqLH^;Wn4BbeG4xmVG+NPT2@JZoJ2Se|M`@6#fjj6B4-PwIdN+!$OtloBi zasUtD@ug|F)O^Hd-6e1H!Gwt?`rYzD+=JI6L?4AGkjp`AX4k<2e2?XJ4*|H+6YyJ4 zYe6=4%j{lUP96MSUUhzd#61N96&29c-lP#Pi8OloDsdOSG3ub#B3c{GS;@(1jSgZ7 zN~d25%c3Ip{c(+NII1+9&Zpg&EOmF3u52y%r2uVSc<2J*)f27hoB#44q~_E5fKMzw zt1Ap!0tmzZ#RXRac&E$KwH>iuvMSy|dh>_K9=YEFC?|7R&X8Z$S#>lqEkAfm>DmxYlUzCaW_Ol5=TmFzgeSj4S^CI%Q;!SVAm{h$6;Q`yvvNf`ieCT98i zM49P_of&}me;@vQU5e@O^RO8)^BD4QvzQsNa&nm%8=4v#8S!u#adH~7u(O%53;zGV zp~s|R>1hhb!^xzf!}3*{Rqy}VKC}I^duAAKmIrhA*1zpFK!CM2U}`TW7K4g0=g^yY z1iloXjIEv3RsJP?Y?unT0-;>8$}R5*3Jz8MMMe=GAsVbF%-xX&@|NT`ZL99o`KceSz{wUUN=DX() zqGQr`lVpS#2psb7CI}+K?@|IhJ@Aq*%ds?kHHA^>mqVZ2aG_C0x=e+NWkr;zS;c0k zouvgz0_x*AFGX;QY-+gm=xV^1ax2t(rg^XmK~)}IqT(70u@W#`t`h!>O;vIJpc0tz z&vqAfCbd$F^NL#1$mq F{ueAo7r6ic literal 0 HcmV?d00001 diff --git a/os_exercises/ch8_exercises_solutions.rmd b/os_exercises/ch8_exercises_solutions.rmd new file mode 100644 index 0000000..8c9dad6 --- /dev/null +++ b/os_exercises/ch8_exercises_solutions.rmd @@ -0,0 +1,106 @@ +--- +title: "Chapter 8 Textbook exercises" +subtitle: "Solutions to even-numbered questions \nStatistics and statistical programming \nNorthwestern University \nMTS + 525" +author: "Aaron Shaw" +date: "November 11, 2020" +output: + pdf_document: + toc: no + toc_depth: '3' + latex_engine: xelatex + html_document: + toc: yes + toc_depth: 3 + toc_float: + collapsed: false + smooth_scroll: true + theme: readable +header-includes: + - \newcommand{\lt}{<} + - \newcommand{\gt}{>} + - \renewcommand{\leq}{≤} + - \usepackage{lmodern} +--- + +```{r setup, include=FALSE} +knitr::opts_chunk$set(echo = TRUE) + +``` + + +All exercises taken from the *OpenIntro Statistics* textbook, $4^{th}$ edition, Chapter 8. + +# 8.6 British married straight couples + +a) Husband and wife ages are positively, linearly correlated, with a few possible outliers. +b) Husband and wife heights also appear to be positively and linearly correlated, but very weakly. +c) The age plot shows a much more evident linear trend. The data points also align more tightly. +d) No. Correlation is not affected by rescaling. Revisit the formulas for calculating correlation for more depth on this.. + + +# 8.36 The heads of babies + +a) Substituting into the regression equation: +$\bar{y} = 3.91 + .78 * 28$ +```{r} +y_bar = 3.91 + .78 * 28 + +y_bar +``` + +b) Let's write out the hypotheses first: +$$H_0:~\beta_{gestational\_age} = 0$$ +$$H_A:~\beta_{gestational\_age} \neq 0$$ + +The test statistic for this can be calculated as a t-score as follows: + +$$T = \frac{0.78 - 0}{.35} = 2.23$$ + +We can look up a p-value for this in a t-table for $df = 23$ or we can figure it out in R. The `pt(q,df)` function gives the proportion of the T-distribution which is less than `q` with `df` degrees of freedom. + +```{r} +(1 - pt(q = 2.23, df = 23)) * 2 # To get 2-tailed p-value, we multiply this by 2 + +pt(q = 2.23, df = 23, lower.tail = FALSE) * 2 # Alternatively, we can get the area of the upper tail and multiply that by 2 (multiplying by 2 works because t-distributions are symmetrical) +``` + +The p-value is less than 0.05, so with a traditional $\alpha = 0.05$ we would reject $H_0$ and conclude that there is substantial evidence of an association between gestational age and head circumference for low birth-weight babies (the true slope parameter is unlikely to be zero). + +# 8.40 Cats + +(a) Here are the hypotheses: +$$H_0:~ \beta_{body\_weight} = 0$$ +$$H_A:~ \beta_{body\_weight} \neq 0$$ +(b) The estimate is large (and positive) with a small standard error, resulting in a large test-statistic and a correspondingly small p-value. This is likely lower than whatever critical value was selected for the significance level, so the data provide compelling evidence in support of the alternative hypothesis that the true slope parameter for body weight is not equal to zero, and is positively associated with heart weight in cats. + +(c) The confidence interval for the parameter can be calculated like this: + +$$b_{body\_weight} \pm t^* \times SE$$ +Plugging in values from the model estimates and 1.98 for the value of $t^*$: +$$4.034 \pm 1.98 \times 0.250$$ +$$(3.539, 4.529)$$ +In words, this model indicates that for each additional kilogram in a cat's weight we can be 95% confident that the wight of it's heart will increase by 3.539 to 4.529 grams on average. + +(d) Yes. The alternative hypothesis from parts (a) and (b) corresponds to the 95% CI above zero calculated in part (d). + +# 8.44 Rating professors + +a) Here's how to calculate the slope: +$$\bar{y} = b_0 + b_1 \bar{x}$$ +$$3.9983 = 4.010 + b_1 \times -0.0883$$ +I'll let R handle the arithmetic: +```{r} +b_1 = (3.9983-4.010) / -0.0883 +b_1 +``` + +b) The null hypothesis is that $H_0:~\beta_1 = 0$. With a t-score of 4.13 and a p-value of \~0.00 there is strong evidence of a positive relationship between beauty and teaching scores (we would reject the null hypothesis). + +c) + 1. Nearly normal residuals: The distribution of residuals looks a little bit (left) skewed but fairly normal. + 2. Homoscedasticity (constant variance) of residuals: The residuals have approximately equal variance across the range of beauty around zero. + 3. Independent observations/residuals: The question does not really shed much light on this. Looking at the original paper, it is actually a bit unclear how the instructors were chosen and whether or not they are a random sample. The residuals do not seem to have any relationship with the order of data collection, so that's good at least. + 4. Linear relationship between independent/dependent variables: From the scatterplot, the relationship appears to be linear (or at least not obviously non-linear). + + -- 2.39.5