From 8ffd07f1ef882f16a28a5358275eb9a35644f2b4 Mon Sep 17 00:00:00 2001
From: caes <xaw5719@wmich.edu>
Date: Thu, 2 Feb 2017 23:47:24 -0500
Subject: [PATCH] finished h2 and added R log

---
 hw2/.RData    | Bin 0 -> 20119 bytes
 hw2/.Rhistory |  57 +++++++++++++++++++
 hw2/answers   | 155 ++++++++++++++++++++++++++++++++++++++++----------
 3 files changed, 182 insertions(+), 30 deletions(-)
 create mode 100644 hw2/.RData
 create mode 100644 hw2/.Rhistory

diff --git a/hw2/.RData b/hw2/.RData
new file mode 100644
index 0000000000000000000000000000000000000000..b6aee819ac1c511db7b278319b5242e806cd18f6
GIT binary patch
literal 20119
zcmce-XH*oywl*vvqM(e31SO1$5+o<djEIUzR&p9LNY0rVk(`m7vw-BBbIv)ZkuaoT
zfPo>+FvIIP_q*R(=bUxlweFv{e^ftRySjRJt*U4B?p+Ts?C!099z;kVb9dF*x#P5#
zc?8ez$==O=!0AsNldQUQFwrdVOPYD`_8sN`+3X@I6@pH}soT$MCd5cz@Q~z%k>0%<
z{E6m$G_%4`z0#Pr*Z87r%NiCk?gay`dr%_|>M+{a8QJkSJ`dJ>T3tlb!^Yx&gt1FF
zJ7wR<iVCGCPVwxRnD085QSDVuoo7o=2TN&OM)?!<x!(Ts$*@J`o^Xr5hmjb$eRv1=
zQ+xgn1%iGu#72OFu7*rU_IIm$&Z%V159#D4=&&kdF=!e3I}N_y9VGTSry{7ftk?`F
z|M-iSN?MRiN`B%d<ypE1XBz|+x)t((XNFlzXwE)(=B%;Ge9Z=)Q6s;5hC0hl<-Wv*
zA$mNFR_JbrJ7%v~0DmebRVdc9(_VZuq~FVxWK|;T*fq~%$4JowsFo&Qm2_Xe!{6rV
z`|AJwgI9<vHSZ8iwPS$hDjT->*KfM{5RaF0*};#Nj9S7@dz$>ur6p8=yyF#|+e~#2
zx;%3$$Q{BfGPZs<<Y59_JH10jruS!W+Z+2z8mT9q)|~X-vxy@*GxdlmfG-pBOn!sg
zzAm9b+$-+((B|0YH=+zTPPZl&X1p5KZk&Kw0MJ5_GVpKkx2-SyU4kpYH8jiZ%92uy
z>agt0*?w_OTd}jdxC_p2YKpxtYmb57S{Dnijd6WemhBXMx#9Z_6g>L@I8^-36g1P(
z$6-}>7nCz4_;&f$=}-Na$Ni8R;Qs9WrBMTjTr~I7?f!RAj^mI3R3Voo9q#diYbFCp
z4BaLf(&V-%Zw2God!K0ims5T&camrbt|!Y2Sun8MoA(Tk{P|y}n;F!*cO<nJ3hl?=
z@??3KMRHR2TDd>98?!wDd(6`(N7qrxOMa`K&_@$OkS<)|M{yNcgyZ~PcCSyzP`9h)
zd;)fV%a#be)QPF1b=%eWvrIl$JNN}VD98z#0$h_LP~yzgTO@w`&4K}V_i(`THf?GB
znn3(*3(9E9GfQHK5rjzO(RU@XXiCIWOi<e7(1*9&luL96vyaeHH~pR0v=+IHv{AIs
zjVlF0gOz#hyPKMsm2KG&;r3c#u0OX887Y=f?{vgL-B+K=lfN(DJ5>fnN8#Dab?uZ7
zu7E%2h}~H1*$GzUG}*@}){oqq$q=76aBjy+&dxZJs+)qE4?gc%BO?ooi(?tmc4%0l
zs<}fq)nY)OPjH|BD15&Cw3R<%V>dzITnTq_w2=K$(unJtQU}K_g4grPwJy8)o>{2m
zaPAfv3w<q}{K)$YeL3`t&ikbWa16&y4fVAc>dUsBMVtu&n%EDG0o}wlU6fm2DD4(q
z6ji-&9s<QUN%d+4yS%QG_Bv3Zic(TVHHcf77r6_dqne1&!J+-THK%vp-S9l&nW*n-
zX?q_=aw`+h(+(P;ldnNI88YrO?farfJhVkH@5CN#FQQ+MJ<|@~v|yQG--%2B?ni9~
z^Pg$sU*I|?KzFw<2z|hQt^kXR$7l})YP2keJ~)qiD9VvvA`6I%IPG<jIH@_EJf!9R
zM1$+$87iM1n$@7E!DxB!*Vz074(Y`}dm4HhC~%SZzBhgBLrfT&y)4}ZeRTJ&T^VdP
zelSG4HEe;7r?3C5u2K!gP)sY!@5DGG1alTsdt%TMf(msAKYxb}JM+8xP%uk&m8Iv|
z>29|R<x~n&0f&>-0H@g4z!cTO54XA-D)cVBaQ5Q-D#eMr8Z|6+l*B%Fq?Tg}k^L@2
zrllE1YX%v{4DXwtfokkPNty8C*+lq*{nGh;&c5E__i~9>AGjc`7N8nHu5`@_b50sf
zji>R7w(=U;-=!)xBy?+QpO|)&_9|C`kGSm<04@U!!-P1zoTHiH-uqF<iIFL>vSu(q
zyV+&kq*Rn8X3=1BU$!a((cDIHWM?MWOdeqWWVwY*nnnzW2hV$mrL+NIo*)BGA5p5M
z5z-=kk@LX`QLN=jQLLb$(fwJ5(iM?+(Cs@^#(FYs<=fw`a=m{OGB`<2JlDoigyE&d
zG*4=idSjoF)TuD0nNu$ng8H5;0Z$(wBVDNWvY|<79@^o5{m)@EYRhd&8sBIOIx*6d
zYsrewob~fYx^r)lY_<}&$=gav*0d)@Ex?ot;z^J{T;DZkLu1j}xCutuL|TB6<bxTU
z1U0hb;N|Yu7=6})!k`9nhu6D)So$S~Q`SB3K<@_PM*?U;A!wdX^Wy7wa7;<>C&0MI
z$G2X>6Cbb-0vyPuq(w}TEs|B#G$8=ZdB)>n9>6C3{Q8yPd9L{G{5p@MeEBc-36Imr
zllFciM)L~t<p%O)3l_3txjCydd18zONm|gnf_v+;ss0wJD;{o?cZt{pf}kN)Q4E0m
zfmJn^`mYHY=6OHp%$+eD{~lEB^Js}yresYDv<=wz#lM$OE|V1H)T@K-hbGdTu5b)Q
zI1HixQu@jZy~(2-=tLWB;wNt6x!m^2n+;F27m(H*;_ISGQvp1ymV%P(9)+CV%JorU
z4CgP|VO5evMK*|)$TS<bJO`!qm)#iZqU$w<DkA172S~<KMG01Ie?VUJHS9{)I8Si1
z0+yceVu+E0TO;iHKx?K|g07bLt?FDeWJo{H^)mte`jDN(q<SO0Zud?fqH<|bDV-HC
zubt<;li+jnom0`$mO`N^X^$kZGT>mGT^~GYTPRd!jDL^+;aw<n7{<!o$Bh9X?w55V
z?8BGKz9-`X&WQZqOi^y$I+F$kyD%UOzO;o-l5@aRxRY*3!xQ63HgJMH4Bz$LEONoh
z+-Sg<OnWM2EhY=x<C8GUCIwCF-@i8^-MlXO!#kQABfSEf{5yd{cM5O@-gj$w@sz&6
zF!pp;-5IhoLl3geq;@lFurlwVlbJ6ljK4#ZhIj4Z#fv2uVera0MHp6OOxow6G3YSU
zNYa&NE(}wB_Ad2Gy1<4>0u89+CD#%k>NTosH@IQv0a09%Vqc_2xvZ-Z<a+PQ9Q2DH
zk<%i(ww`O_X}ObvgvZVs#ibatm=2*lQCJ51pI9TFC_P}7{yUu>;YiOaoYbEvOt5FA
z5An|M68_P3A5yO7`uLg<@tc+AP$u0UrT@0RUi9xUMa+g+-lSB&)^~Y&?9g$c+On=h
zW#{bsW@LFmr58b{+JbpJ@fAGv!fKX`sVb)^lt$@`gD)nAL#ys%uiR`9{w4Hx=#KFX
zvF647W0KeUaSt}>R9gm#w(i~0Wx(`X?Q#9M3Jd)0zHXd#dT-Z@k!QE@<gS5cz2VD=
zaevpZ9imirjJLcC3f;*TOg$K)Z&)XYJk(prPSC_6<^++1_gz!!KlD0fL-C|k+)1Pe
zsj=&;iHvf1JuC6n%I%YRe^-f)34J0@+oJoYZ!;nc$|I29xdtl9F+v8q2O22P;DR48
z)fN+i1%f$lQUX!%TP)!S5oilc5{_Q2IVlAY-9Mdg`?WD@AcG$7tEzxJ-92L-K@cI{
zj!``<VDGdYmdb=@{$T8I>$K3kA^&sdtq}d?s4}XmJ8Rx8YlpF87bY^nMdIleboD3m
zOyjL>#khHd!Kw%qcYrIwf;rVz@04H1&6@>wAli!j)QO<{hfAAwYl44~yX0>3?aR^t
zcm3GAw3}_by8wswNm84E_6K2a%jdrXw$5G*G&g#wwN#MplAMo`cu-4ub65Jc&UJd{
zY{YxK)<=?K=!t!N$C)R*nEHC1cQwFcOo(SUa|F2B=mIOlwG+J#lIs@2>!cgyUb*LG
zWN)K>J)7q%bPK!Vcb)imwpU2gWzpR}eNLVrt#5G68Q>`Z+`PSO+^R`>W=k<pdH2Mz
zo3*0@x=LhoGo)_7T{u0Q%~gI2Pga?sJ3kcWSVGWTif33dHk3G^Ib|Vah~@^K2mD&&
z0i0Mf^dLSTI}!BYv{%3#k%-PCL`SC8^9*Xqi7A;jS}<AGVJU*#xz3m28T1YPEs2v3
z#y>DYqHN#>gy?{YZi%ps^b+ug2%RV;X~f7)BnFrE-tLyHf4IN+V#02btO79ge97bV
zY{?bEWOGxbXh2gaUECR2{G3buaiO~h5hqTqp}kYGmUIhDjzy-BPqeFAWvVI9Yy1Gs
zgKsrXjEhbHt}6!eH(ovK=NqUeJbV9>BhXXp*10ADs^K6m0Tjj(92Ad}(m&Q9I3Qs+
z1Zfg`mJxZ(z4I<_zGD*%i^#|#nTZfIU<n^5HpnDa27vrr<g#0zUY7Ya)e`n;NVRZ=
z-fHYzBLn?7|H7I6s{+JDB$4<XcVhYnAPk)2_HWT5y(aA#8_pgpe@cGsyT%s`C_%{a
zkU!TyeLj|`H&J)`yp6j4Hb$95{T*a9NYr=7=ISBKsnT~{F$0MaogX{C9c~moh)3e(
zem;(D3r<f?z$F_VSv%$xb|b??ZT4;ckRcMBW}8%t;i~3+EYJ43kwi`Q8fRdGc1l^i
z)vK9oD!R9SUFpWul5S9L+KBA7E&+Gnu5gW#NLGw<2lzM*x=zsTyAh;iax@T~NUv~w
zzcj%zZ;)k?988C0MRG3W^pa&qR-7`3M}7CnzBNXY&~?^fR9L=yW7?qf=SncmjU(XN
zPc*k1@!+gRL1v$ToH!ikLNHGd4x5^NjOndow71Cyw%4-m(UD%idG%VU_fmo2+_<^z
zt<v?KlII3~6>T9rz;=|yv+zX!4_<<g=jJK51MZ#n*p(@Bp>)`CpVePRP}WT<vH~HS
z-;aiv7_R5CUpBtKk`8d_Ig_)<U7>hK2RC=z`L-+!`jRx>Hao<w-*%_#8q?*jdGGX{
zT*raVLn*HF2ajq{Iy^Igt2YmqSV7qQdwaL2H!1hBfj=3y0cA9sG|SoVpwijEt58@|
z2&`eF#ln!jQYwiB^W^FUj=9SlMJI3|M2gGFFj@^}C8)Ne^X2N3mUZ0l_InnTR3Hr8
zB)tkbe=j;h?y=-y)D{L)4?B$}pOeT3c~?)-wncQ>Rd%f}KSI$p2J)0kh6hQf&Z=F1
zB(=W}=_Z6V2>h_qdO<-1ssTb$sqQ&<+cQQ=H$NY-B1GQ#o*Q6S{^C1?c0{H*?2Hk)
zO53hNG_squHu16jjr{Bq0508}0}%#2-nvVLmc~}Pe{57%T3@baCT#!66E(~G_tS?G
zZFL)QOeZlzxmVFdV3XI~W=Odubt97l^LMCeC)gzI`w1b^1cpLBV6$5+fRl=ZKH)-W
zq9vK|Smv&lyHbD|?y9gl&N1mml_>k?<WJYVpF=en-nDmnr_lt!-fjpPTyo|Z2yOX6
zq_74n{)9J}A8Eql-(4S_gd)MuNQV6IfBic37q^kj*8e$QgUL>td`r^4+HI$hI|-)U
zv4G~0;gW=`p|_hX>=J;5Q3g909w73E#}HZuOXl(8B+EzZ`c^CJ1j!~9(%-uhQh^~v
z=NPn6;E2%^<ex|oGA{>$P{D|W`r2KkaGesSADNRG9BaJXw|&B9URBU6%S_r10OrXM
z(OieXB)~wyy^<&*k$uBfDT&!LQ9Kgdd_!1s0}dp=8J<Ah2s>lI{$n?<c2EBLgXfm6
zzcD@@f)X^{&}8*P2|V$Mw=T%{K)2_OmbiL5V5L;wu9Z5H3yNdK$Gkts>}2k<o;c*+
z+~^cQjl||3ZPT#ihLxEt;`J88s%ur{VeM8&(5J>ZpbtLj#whS&^-J96lh(rWm&Y8k
zB)rr&v6`1+Ucr0~`=NaUH;+Gup2j%dfZ&{EE)*357e{k9?4VHB*NTZ72fDgGSIsT2
zBNLO9Hqtg64`8*~y9Ll!at@fgZ(Ms@nPPsNJl@dITdxMJz99YFesb%QhSbfl^u#p$
z<k4sGfPLpuQ14JTJKQ(#q~hCM<X6TQxVN2AjO<o%iFe{?GU5`I<4Dw!A14$#x%iL$
zev9%>a$+#ilJBIugb`^y7I9>i&pK8uz<)F*48T4%`+6DLr9#@z7M|@KpZIN8Br{Ju
zPgJ^+ySJNFJ!6Trz3Jr1L_1kcfp4n`sUWL&?$_~TP(jv2XjWDWkCWe2p-eWZggPUw
zWI>jA10>>HEyw+K?+{VJ(BDufOlF#m6%0Gb&syPBVeEL#EyudkH`@UeVykhkdR1Gj
zN@6mSdrTQ{q$VcI__w{$glI+P$;20rN~x>fCve~F6CN}zaZ!EV${wR^B^@z*{g~o-
znd48vb;{VSc*X41$4hsjZ-k=3lR^yEStIxIoQjLYNqq`iws;?3M&?$eowK*kXSMX8
z#f}$<#Hmu*HHK+AsE3GPn2sZJ(ppwl04L~|Kuh|@X5>(Dv>iBF1EnUSru6DCk4H^u
z)snPaTHfE+n)@v`-Ej)JtKGT@8NBR0`)4NIS*ZeX0etk2lNc$FKSI&Nuh<-K2wJF?
zA+l19)oUp~PHyq}b@7pU&W<(#+3p7lE9<{hW_)p!_jerBGb~sEui{BkhvdSCh*aQD
zB{Iu?AzYv;@}?7`L?S`cbz%$B9>c03?c5g4%)Linl^N5>33cdm3<{0~16V`_x#CHd
zhq~u<M9<t_J{BDsmXZ}ctAFztPs$1J^hp@1o|qq^D``Y)4obn{x8c3Ka7@k%WG*B=
zPF(auWlq$|qh2ObENe@_SLpbQC+JLBDUWZm9`<Ehzhjminz|51Y(F+?AIBxuJkA)B
z-d3I+3uw~K%M$NUFXLUy%G#<uc{f@An6Vfi>)oL!VR)c7@pUq4jw@u2b=8x!{0B(r
zZJKh*7fRC5@KB|A8YcCry4&aKj?E+irD)9%J$Pr!iSpI!%ZdBf7CJc|IbZc3Y7<;d
zq#ONWw#~B-=~{lgmgDo8!Clty!N3%iiJ0v212*O5{Wua7DF`lwkr6$^2RaIUZgl$H
zQ+$%6kyF%e%cqySHEBkg%Bq|HC#wP+UmG3t2sSkA(#p!}qT|u|sZr*XmqEi(s!{0I
zicO~6GAR%qP!jNXX}_jD9bsF&rEUG$|B~8M5m{=YC^wfn9W$b(wH<fAS#mJrC%<!&
z$b&#!z<Qj#&OIhAtwo+M>Q7|HJ}He=6ui|SF5vR{OqsV{yjQnpoW^#)R%gn3&UyDl
zy(!RiVVZM*+5cj<nvefraXT+O&DuwEt`IY3GBvUo>ODrEJ?v4GIF<QLvF=aLNlCZt
zT9;z;t3e~DFZC7t?9-O@W?vbS8bsLCeN5ip78E8WlJJ~^t1VDbQBY*>1#dG@X~`r=
za`ajbRxKzSjdx|wgikIgex~~+^!`0<+VV$ff2v;jp-Sec2^O8HYNrF^tZAj`-$Eam
zHs2z@kEWQ7)@2F#1$3rOifE1XJ*GA{)s87OOLeX6HIp5yOFg?IT_Y*N!Y)R{RQ4#`
z+%&bC6Qi=7{bX;fhG@nx_0@2iN&Vwy0j5|?e@Z<68<|X_>@F>XhCos4`ul>*Gr`%F
ziLu`%gH(89`Sz_I#0wgY@SGM>d=RxU!`>}D%lgV&H+M@ot2P5%CDLnIogDl7fyLvg
zPdy7JL|O^Dp5Gn4m{a+Y>aA1~S|T1vujVrhD^;$7*hFlqC!*`LVjXFNy6(M6q5WK#
zQp$n3dRrwjYCS&jMHS0dZepoPr(kWx$J~7CuR9b5mxbprc}j1_-A{1n?3HIlzYa&F
z_o<|I#Y7v=-&TVgZGWpN)7tQFSNWKnly1Xj*7JFDd??SEkzhS^X(lD%TU~P4xm8aI
z^Legw%6-R_xSxur&laugv+Ikf#R}__r@xL2{Q81&Qr-<sS+5j$Kxediye_@>q;Wfk
zkgav8#YsCxF1DBz@Jqvv=cJ1HzQ{LyUqHmrz^Bf4@L$eRlPoeH?=o=gw9H19^jVXh
zefN=1i0DhLYBzbAZScGVW+s~>8(pSUotI?u<`9|ng1>6y%^U540Sni@LT3q;B3aKj
zobHPI_Mi6>R03Z0XYdBjb`_1~VtAjNQzX0)vtNAd<f^XNBstYFqCwtdp+uc<bJ9}I
zi$M1&no$DRa+QYZUvou<S%|CaTAPFopEBZ<9kc~fsXM&T-6M;G-7zDBntbfS5hBU=
zLV9>!EKus?u|%F-q8L-;a!bFNb(+cQI@D-=+05xv&UXkP%oeJ3;5}hVme;AFOb(CB
zrfywSyfD{pnmiqA{ZLm~#q;^1&)H0Gt#{1dJ5B5fGSt%Ugy(5(0y<m3Vj&txSH+GC
zv?t=v6pc3Qrlf(cNuSjen3$V%eWDt@tNC`oqk!$3!c>x%N4XQVfsg{V>DQ@`?%!Qi
zyYn3#jF$%3(xM+dp?Hoy_9jTpOF}O8zW!-6T4Op9`>yoSx9#jbk#C-^n1^>)pX_aK
z-B%$2jT)$2@l8#YYJ-+Xo$F~WCr?I($Th~4N)GPx&RUM;{|3hU*L9KaP|}98y-A4B
zNoR5H{kGiq>Ar85E42iL9ED}mkcYF0yZA>)kz)+`EAk819Z8?L0eW_cS59Bn3hY;6
z;!A=pY1MV@W?1`__Y@++gQs0-7eu~!)cR$NvFbd1syO{S%E)!fIfxA^lm|f>{}w@`
zqinRhngpCg<+NM|iqeuTo73)HQoBy_%WFq*3?(=Bc}pX=oh>czKi!eUO&F?Ie`I@O
zFilI|Y@;1nlXxiYY9*S%d^;33mVx4+uQ4sK5T2?#7kDS=ks>nxPKi{;LRB_*(nt|F
zL^$|6>f`$uy<)OR9=)2?gO$3R38b^phMsWTH-+)B33;xA$z`W0vho{&t=Z>0RXPg_
z3Vd<t9Yx#2j0v$BMwKjDi?51B)r==e74;O$&gjiJSJE9!8x5b#_<K54O%qpSI%P_V
z2tDDi`C--(n5wK*E4G+2tf=FY-ef^~8q1~tq6(v1gn4|5C3<S4xDls;<wv5s*oeu(
zxL0+_HhHwXq~O2Kmu2WY`Zbp*jDsg#C*lif5C*^Zhnoec2EYe96yq`iC#uR=WbT8M
zrrn6uh@Yk&bCkuvl-^6NgTaKya_XBe<uID9&qGC2=|jR;t3HB-jlw6&RLUT;VM0Xh
z`4_X1TREk#D<MtQNnvRucQz_wCO<1KeQW$6;%z_ihS^7;`YZ9ltRgve(+5TGC#$0e
zgPL_j{<b?mdM~4DO~~ntX#~~07G)#sngrRy5?;L;+^Z!y1&`j>`K+$}F-EQN_Q;?*
zSD1yr_JJ8{O~!zia#1gSqwsM?TbO6yV<Fz4KuF53q$dJXRil^&;j`NrzXGSj!v!Fy
zuUZE*$s{5QS*{fYuXR|&900CQM{B1Xe5%qtL#a5+yS(UsqOS7_D*O)Bz3ql-w(i%O
zm`rDdKWfg=N-Y1HxNVtwA!EKrb`-U%1~edD_fHBg<TyYN*F?zXesHIiK-QXvZ#2TY
z@9c(r(`Ggv^NdFZtS+9K<IZB7EHa+vDvv&Ed&hZCOGH-n&3)ZBDeq*iO2y~WhmO37
ztwd*S_`>w><kP79vR@f;V||mwoI|NlHqO%|6Al*w^gz_HKSIOvL_9>Y{t{JH<y~&=
zqEa(g$sTO=NHY2{Z*K5?t!-<v)&)}`iPtm?lhhmguN3VglcJ5xOtoTm^+jreErx8K
zaXG%K1P^C;2?}+a)#-#LDI({%NGW9EYt(He;C+>aSuF1APJFF(LuTqHy<;X$Z{lC2
z|EgZ#&ga)D+^U?`@&}z4cfKQ&Y#!kd&8dH2zR~^1-gGm|mBaD*eh|Mkb;_BRH{I{=
zM9JNJ9DCvE?GvLOYCS(WoL;wT&#nX-NsH`-eo3*3R~g}K8Mm4<f8I;#D%T{^pRemK
z{FS}Gxw_IQ&eX@@0CTrF=S+IB<<>!sR?Oa7;o`YO^~m;1A8+yHqil&dwQt<Jh75-_
zwn<14>xPQnqqW5GG1}Pd$9k3kYNECU*;CHs3?G9_b!zO{8!OAY*kW&c!qJa(%^*#c
zWrR1BMUy@>WMYN$_sejyYr=qaYdHO{V9Oz@+QuX*fx22Dq8~MW(C+Zae71{bt<sY3
zey(!r3GsRtSyDBzl^Tl=Z2o?_dI)}2RS|UxwRfre39egBiA@lZ&$Nk}^)I!v?F%kZ
zK@W+g>yx~ErS(%qO`FyJS&aU3>F{rF0zws?a$YbAUZ(|q8PV##JPVlau$OzYoJo|1
zndbWvhZN_lqBVN@+%N3D#jt)vOkZuPTW`*EDBaNKx+VVIfO~=Irw;{LB|g+@R%t1U
zkflp7Pv`WP7I82Pj84dW{Ghe);koMvt3#R&d1=+hav!)4LaDcV-d7H%Ia(%u{BX0P
z@`nVqH5nFZGBq8MCSO1v5l^HUB&nr~L<r<->to|;<=0oe(=kOvUg#7{#r(vGrmjS~
zn(67$Z`5gr#OY<%V<j3Vv`woM>|e!WUnd%k4%H;b?#aGhe%rdYCzbGu-}o_^=!FV3
zV!}{ZKT=;$+ZrvtwC5YROdFOwnr}`cy%?)5^x>tN{Cwh%7ef-7PVPnNqXta+Xm4)%
zcolu9{b24*tM}E5V%k!@ufI!Y+c!QJG_l&<YKWeTB6oY^)7n!h^xE+Jrw_PLZbg}p
z#u=*FkM7crPS$IuJtOmAyZuWxGKk%>h0-lIX74ak2szK#my8z`<E?MSy7x~e*GqGQ
zPo1vMk>Z)nG&<2kHhL^ISEXFPuLgIyWX2mPTe%u1{rLH1krm(S1Pf$-eqX?<d0H>q
zdh)VbZ{bPk7&HLqMb6_$l77FYO=deUK8hU?TMfn==godWlPI^K8r|3_DRPg~Z839>
zpslP>%HkutlhnO)Y}=^%?d#`3){#Z_ZzcB{GPB4^uu^jUV5^|WsAZkm8e7fxM^^PX
z*?_|Juk4`G*I3xWMVL|^32GL3<Yw&TGd-O$E7OqaA@8j;4hUJ9T60+A*_a;~14MH<
z{D!?+KQX_Ej-=-)ABiUJyWN<AD|kAeb7xpBTTk%1y8mlND%_yf2D9IhdmO%JEU%Oa
z9Oh`4o<SXK%DF`f$j0bhyB$ghL8A7qljxd@O6;~4*4)wXsJay)J)9M)w~VIB^hS}6
z%U>ocURv~ENsNy-12wI&e_n#5qcc4k_Cly8@g@wqTwqCJm)urAIdA6Q9Hd^&pv_~)
zVuQ|0q+c@tA1<(uzD!)Lv&BZPrJ6?l;RB~5X!;KJ_$>cKn)^w{pMfsZFH3PWj?tDh
zBKBLapi$Cd^zH2|0klJQCMSM>lO5CLnO=CjRAYbMJHH1pXIRQA1f!MCG97f#^nd`H
zp%BcR@nqg8A$%mpPZhfE*?96gDr|MbZjcIV^Cr`dG4MPZCvJ@u(zY6iKkxAt-|tv<
z>`?iVZ+s$Cwevla1v~2_cm5Ur!WN-qazkI;h1!BQeNI@cRhYL_VBg{QpvHtPQtV+q
zDwa_#<<uU`J09i@?a^s^dnlpolf83d4@$?JT3YC!7jB!fx-vyYD5;{0yCTze-8LVE
z9$uFnVdgML>^pz+NubC?xNGa$#h)PBD95b%KXQAVqmeAU$!Gf+BWCzJH`QmI4?(k?
zM;Sp0=u%*%W*Ir;cUfO^3Ed<8!^=s&eB<<(*Qz+<loGjfyJ^z{ub0;S8^2>A*9zw;
zprfE6wVlDX!2+5dfwgBOE-Rkb@!jU8cA5OsF_Edt(CO~jQJ6I>D`s=QPC{Y$9JJeT
z?4pHhtyp*LEl=g}9#szN)#FrCnx5V!!tdUk2St+9x6_|q7F%rac`dGuvg1t0iRe<7
z?<6Hlyn6=eRE*sD4Qms^j1kuGpN5TZ&y>4`Ljo3c_LEWTno+uo*bRsAjOGn`CA=tl
zcocTZ(C%vcgb&`*iPkXiOZR{}J6L;Z+AU%0l=Cwc;YSVEjs}-P>kRVkjIyCX^)WJd
zJiD*tZhfFj;LBc{!q6mlf9NAO4{J~j(;V-jjqS2b0in#WR2E>{%f&)NWluCi;8N6X
zd62GZS-x7_{xZ|?(@uQ-FkFKd2dPX}s^8l*$*cvMww=0%9*~`*p_^-k1<(Ky+qi0<
zDyspr#_N#PupoIHHM}2Vw)8i}(Q}X~zV&eShNDcJM6j9e(z*Pg$r#hRn^@Ung_V%$
z-RoBXaIX|L|H2zhPdBLOaHp}c$f<P7tz9o{35yvTv}Y-LF+@2&np;7*OxJLuQMuPD
zY>~02Mn!sS{0V!1J2coRTi}e8>147z5~Jmi<6Vo~{!K{_g&2D|jaBL{%Kc*TN?&Rz
zwv36=-M+jy;#w0`D|0Rz&+2@D$!=jKnXa$<CSSc%4t;^mVNh#ME-Uwl5nCUi*tkMW
zzo{it<Zs`Xg!oLoTI9l(K&>o`?V>6gIx*VWLlM~W#<B#Qs)uYZm}aE80Vg@^RKi^m
z*^q>UJn*)R<%><Ian(!vOvsfdUxt|Nbw=XDDehhAg+S*{wIQGd`5lM?)&tVry=&LM
ziF08dy{c|GAVCwpg<KSalOpe9L*X*NP8s_BThe5%RYih}pDMllv<jm)7rI;#TUoE!
z%hm~(%reT^M4hW7KE;cK^dk3;+>pX?vxCmr3bxINdrw}Pm>&41mw}VZ!q+_UPEq?k
z-M>Q1>*wz*BH>Y{LO%Sf4jK_cJ$Z)>W)ARriI>k?A)aRuA7&NY5Q45PLKeio5%|%p
zSa8dSANi`sC)WB>T(OTOOCbl6)0>G-W3Zv@tGdhc$euMf6BuEpeteK(BJP>e>-y!T
zO<6~mexyGGaDLHKo~yjiSGO@}<&BqX-~pp|H{NkgL1g#cMe1~>&vOmo%RwA~W+q-@
zf7QDJg0g<m+j3z!{}kV2K9j5u_j{ULcVy#Wm$ytAE;JqCMqZ~yA(Y$@F$i0=T6?qc
z#%;XD3f4$@K|#rQc!Tm<+!@|uf+kBWRk~fQZ+*OY0p~l(0U{(TMIeIS77)XXmx-KC
zO)JqJ8CCOBaq94oT#ocfyFZU>M${Td(&KmXx8GA-1agwV&7UZM;&2v@7eA>O`nv<Y
zUg}hDJ&dlu?vDLB2jx>HFOnPuJ^z)w3GY9HmQJ?h9xOF?Wn5{QGjN+@R5P8B!JFUa
z#_O$XR~SHNtT59#I7lAAJ$CtRsbamVZ1O7V{F-mxfyLhJqL}VF#JO0lddaT-rM{mG
z1T8JahgyCNld|<p+yy{rb&jY8Yg`W8u;aGc?9B-`r*V7T4m~~{E$J5^+70)#<<%pc
zm3sf_`IkQ^DX~w_M4YkS5T{3ZeN2<Btn)m1U79&i#M|__Rb$y?U%Tzw48<RcFlEVY
zn~UZzZ`ntDC&aBDWkd(gmRs9=t$p>o43p|^Bl5?K#ftbhV~f4dZVu$Ydwm9jY!hl@
zJaBhkpFU$^bonXGyasH&K7RM8D!`kHW=Gb0Ide~M{ueJw``80lU@GycH4v&EhHMJG
zB=>qKX_cMv{u(beemAo`7`j`=o64ymTexR)gtN<JI^*Kfj@HWs-FqPHVC^vu_%vJb
zfS*4#(%gYw(sjOX-;U;CyF_}5d-`4bw5J>%VD0L30fDC0GYlpm(kYQei2UhN@=QzH
zqzl$xBX~S#T*IaQh)-|f2{Qfo>P!FjncH$Z^fYY`-S^eJSRZKDRBMS219%ECukj0&
zyI*$m=75c!)Y>7mz944*tix;{u;2|yqCzt7mo>ZyS4DVVe`%+Tm37%njJ3@NGEte&
zf4X6}xsL*g7;FzwatGh2EjM8I1QlfFAeSGP)*f!$pwl&mmn?rin^@-7#<P3$xO_Ps
zRaDEgvwh#<cJuc)f!6oEzTzr2!+y36#syZl_zUL^N|Ba>`Dqf0@?SKH8deYY7QhG$
z>dB@UMLF0nfA0$+{PVW)1Gf(&aaI|aKf1B(%?jMhibbwnUIJpU5(kb!rhZcb6xL^@
z_cB+i=fmtcI%_Ab9w5-X;hqFzDMX?p@<azTm^S3blNllT=3Zqp`+BvUlg+hCvYDD8
zssX}956(+?+cT~vP>e5cVXW$)H)pA<$UOP8zNk57dJ~)ZbuxK5vxHaK0G*X&Xluz{
zOm075(VRv&a+~o91scWl<*Y@*Ml(+6ckRfzP8cUZnQJ<J1jOfRN6Kt?<o0{5z+CXr
z!rn-T!yH|HnC~W85!&Xhw}hE36tzgQME3C8BVC(M?mxNK{bcWPJHKivWnb?t@^^#h
zX^MAXqaI0Wx*2lEG*r;t;`({V)x<<b<$yJXp|>Je)AezGBrV`PB`(Yt#lypRr?0!i
zo%U~qz<6(mdxa3G?)~tfpD2SPng-e@0H$1PlemcRs-G%v_=Y}yOkTZl3o(kzrG>LH
z`#j$4w54fy&&#M1F7tSkqSw7*{=RO`lPHTmrW91LnfSj#+CHX0RIs7A(cj&0{|-&X
zjh;kl_Aw=(g8v-~_c2AHf{n#DXrulEnu~8diE{5_%0LDGEBqP64*eF)`0+%V2qbXQ
z%?hC%1sGo%a)81&gVJ<f?obalesHw?`JeEKun^IX*`JmBX{o+W{PwnY#zs4|Xobm`
z%ni3V!qWnNR`WPge;n%?(4uX8_i=1=pz(u)txB5C;~nbW|AYcN)T50boNOo4bp9Ro
z+S~G`>D=F;?r!`K$gx8`+}P`ATbZWwAJD<pAx(#Vhq}M9_urxT4)u6ruahk*P3K==
zAK6`~nKXduD%%(M^vdutD11AJ_nIdG*oTrHi=cT(#Qp4!{@eR4-)Y?mc#jE<{t)+c
zkOlkxUkrY}PMfto9abT}^N)ziU)|A{yWf&P>(0r0tZnq?e>V7o0<wbA=>GeLV-4E9
z=R1A>+g|?n&;8%u;om+f^FIyP{3E9FS0DA|!vDUq2Pf}|w$WB?&wy3P&wm=Or&&D*
z8d3jR`LC^sdj`mYf6_c;;(qo~|801SC9TKopQGoK(q+WuPe1Al{#OmpNs-}}TK7P4
z;7^SQcwR<{trhp#|E`n^_?bw`#7*^4|8;naDy_$Fo|7+zTfe#n{@)mQ^{)<7B?jWa
zZpa4Q=}GmV!FVq0+||YqgHSKsUM?7l%uSbx2VAS~T#nj{I4?CE=2cm3Y7ffNd>S<w
zHg`Hm{S;vu<tB4*m$~*3>gGD;<s<vu+4fVcj+6X&rmJ{;^|nl;RVRO$!hU5Zs^ey-
zeA$d5ZKolbF-c75(rpn};C6?q#)>&eKoQ#Ailh*+$4&HtwyXtL7nh-jsW|Vg=niAt
zw3mt4;o+_)PmlO8s^;6#D?KKS<ntP^n%mm-VP02}Dr}X3We^V>xsS~rw*QiSx>`Sy
z%OtPlJ><Mt8XIyXmW^b}2Ri6B%o_9y8Z%G5;!0F3LkQJ8pZZDEv=SnU-%8RgRjEbx
zM=8b?ziS?H6T-mfVb4j5ZKmhK6Vm|)^T`=)lEz14xT?)X84c=i+f^M`)aUkw%VF<n
zk@*z`r8%(uK!i}>?bUI9y{qE2Er;n7lbfpaK@_~w5pMjI%l)>q;kD&0pQHUv!ZD1h
zvEwT2{5<jjDmGZ4A=zq`i%vc&XH%)<CQs<zA?A3dEdV^UHZ0MqF8?`k)#|9I@!HYJ
zOo<I72Jw(n(gFoI<2hPm=<}97ke5steNG9A+DvHY+LZuwZIQHRG3JBE>z>4Z6^Xi(
z*I(N0I2Nly_)6DP16&e^TXc2?H*V3PoLB3t29ziqqxZRWxivCgo=F^Tyw<Q5;(;$<
z)-3>w7spI6zM!(9_$II?Qz?v8<XLXcs*hsP!OsmCzYhNGKf+wo7+9fw`}Rd-e<t**
z!SR><bc9ffs^WRWcf66Si8DeGk6dV0Gi{%!My%kgDlSi4k6>mul3M419-<d^Bm&qq
zn)#ovAC}KTa0jQKG9QH7KD*Qdii-A#;zcoZxJ~!AS!^-aS0`RZ_}lu@<o-L1F5(_`
zJ|y7G_uF>IARX|%Yn=zzhg$Fv60O>T!P!pk*%^1Y7G2;|^}PXt=_n~HB&_Em-X)zq
zL-z|;q`yLVv4DA-;yJy(GW$H%S9kpjyl)t!%|WupxNe2@h&m8Ft!q07?=EcLn?!An
zODdIZHcaj1l0B%uaMAUi%>RSr0XG^gWw?N5%EmUz5{e4;v_|3<6I-{%4r38ot9{r<
z^GZ86b367#i$s@M{~g`N%ZypnI}um&V|vZ^>!juGWEz~iPAcsqcHMN&o2q=m#)bVh
z9f~iL?zv$_xu_mTQ_DE#?*xjhik7TcrMM#9Vd0&kGvg@7a>EXG!DV5t$blx~yFd#J
zvOAKCJ5%Ak2mEsu08`p|=966Z^&JJ}#g2<D0KfZei2{*y@aESS#P*`i*0Twc(^k=%
z;zDPJC^{imq`*zu0kwctSC#Fa<0sATpAN}Z@@^?U&CcCd431aekDM&T_$uy&Y1oBd
z+0`c4Nw*rgYOPhJ%H_7`xl_Q`o_(zelZPFay1H~xya+19n>)Z!i>T4OdD=xmlPRHm
z7n=|_X(0Ai{l|Qk2(SXJmr=l8d3md%OIr(g?$y4;I5T8j*x}bZ{ry$hNRS)UdBbyr
zdYZRh9h(8k+!Hojgb8jLN=^$MD@!GAwt}3*Pg3(wz3`fIqpb#duF#*Goq;Q<={joB
z5^A1R8nt&MNVn7CR@KX{3*RLumd#k@EQh?QWaS>}M3*8LCD9hcfmPf~Dz&M5ZeFl$
z{<4c#`a6%tE=q^4_uUuE4E1LjDXNs<!sCt(+?p!}K4;0RCsl*P;J?%S6%~AeLP64}
z^MoEjR>D+9v>geH&to~1UW38hE{q>p6uY?}U#aBN@u#?4IsbXP0oA#YVaj|tJw7}t
z0o8ju0gLRf-Y>PKNn^tjwBCw_wc1sunTCmLK~z)b`l^omwR@D}{Td2U?|;nTMjh&O
zj>Azpr}9y*hahy|IB9>b+3=TB(!<II3CO8buLs?xyhb2n`B{T0&+m=p{K;s+DEawm
z!qkIuS{d;%u%d%II$;%)>7{;%TZ|O5Y<bS?0x;wL^z>_(|Fu0HzlqGw#~SyQXn83x
zyXf4kFz*HJ<i*^)Iz?;$tXhU|*`B}VR8VRXzs&9nwRtbf=VFMHTC6&XKCEBWs~?0!
zvvej(O2FGT@;A*=KRdOPWut!g`O;tiDF<u{G5n>M%(e}5h_Hx-x9UeoaXp-r0bnB{
zVd#6%b~DNZ%Aixsvpw?0*=xT*dDNbS%YYaKL>MKI;u3#1Z;v{&uR0zL=81pFA6w|f
zwXMhu)vTZv1%M&$=(}vE_zi#wFw+K@{)E4)^6Y)}J&*dsA={`R+`BagH(jlYmbzYX
z@OfQ(h6iSRu1X^cKAg~&a1;&O|AH?Sd(U;+h)jRg-~TA8sQUGc$5Rtgufx>h3!HPL
z9r&Ikc1STr=9%4<rcNS%{q<JIDhWwbdt(R6sP8`F<NS72artoHFWbY&L`H%KVlx+R
z;tv*2mo&SNBKe60&sN=RcqGJKu5l?&p2JO!ZQU&iABAnf!C@B1Z<X?bfEvIVuAlDh
zj)iKlK4fsa!?4n{Si2&cTlC!H4I*Vvxf`@C>^X&FZRt&`4r^yRL}1Q=6?~H@`;Yjh
zdPQ!R>$7g;_I%skd%}7cW*?<==O3Z1!%yAjqrce~^oe}E`JD9z<CcfRegV;^S@LC@
zkW9B11mCT~>G8(7qt&0fC4mD$o`Np~RiTte@sh=el}mrso0`vEn-Lm|S6z*>W!&Zx
zu<c8JWb0QCHuG<-q%h;=?RP~p+MX^L8uO5&DdIKhG6v=~%H}?#Xk+FK0z14sPHsjC
z%iP*!v9*gCaI7~flfMbLEQk<mzfrz*Z<L88MZWYyR~OV#>~gsIrNAJ|dGd|x7Z(qw
zX*2=17JZPNn^*lHOLs6x#3O$_vzh_ca{Y2^RrlS~I?S>A$vtrYx5T8DyxjCprGGBG
z1RJ%&R4)s))%A`<@`r^KmC{#5Cj}5a`ycH*_T_iLzgzb{b-kGe4kt)FX=`U<OLezB
znbCZQ7Fr@<;T(HeYo-8D!n_Vty)EXhsM-2^eANC6a@khoLHV7Seef8xW&DKh-d&!#
zMIkSN-wQ|I!QGVnpL%BylGrs(E7ACj3$ILAWK8H@P}7mvmg0}McI{Nty5|g5N|w~<
zwYcjcsxPjy(dah6@ni@rX~tBYUmMJrp45)eCfTSq*#1pGfY`nIkeCQKyzhEYq1A3B
z;Y^1b!#60k29m%(bEy^Jd5$A>oxZL8NV-sesP|=)D^l>N@>D@$_Xy@Qmi6(g=pzj^
z&~zqN9k7C4wTc97>T#>C{d{lv-CJ;uaZ&K))I8d#czy7)b%q?NUp~8ql*@bFqr#9?
za6wybsvbZ6kBGL-t7M1Y>9IrkRNd)n2&}RW@c==G6H;V3>uT~Jx76u(c3Co>{&|ri
zM*g5J4_kHR{biJAMl|<S^$Wt;x^(MAMuls-{77AwAKRh&$klPX#k}_3y_M)SGibc`
zbUbPYY{ln8jvBrH`!dh8{t=@IK1_&hmQjz5tM^yK<HU_Gx}kCo|8+srxevH_4Tf|^
zX%+*xgB>pY`&%V1lLW+*_UiM@#6v7*&6isMJ3RCrL$01PL2G&2SwixD`76u>AoWiS
zD`y11(R^*s3WwP?>JKd~I@$4g=TvT9nM-n?{0i}IfEjwlrvPu369pfQk0d+2t5@6e
zv0qu7T8{U6S0I0L`WqjjB91_>6o$E=xe2+z=jFc?-LD%+zdePspB;Dw^w@y}_YuEK
z7n+iW&RsU)3?<<B9hsb){WL?v^wqr2ty{$~7ktc2lbG?R>$vwHHS#EFnwrPUzedF@
z_=4%(`8bnmXRmy;nkOVNnV&Qv%f>=hlFt?Lg@C$O>9AzT+0X~;B4xzi{i|bP$!*MK
z2Q{5}S7>Is;`1}MckO!}lawqg!v>Cm?Ip{ihq$J^#q@7-B&)5Tv`aXb*F@ps-&E$7
zoL$l%<iQI&fLQtD-tBD0(Wce3)??ggk4ZXOM-iDxEAdYyRhXwLL)*3K&Yx$GAUo5H
zM3%xyz{&Pq?zHm0abaqf1z_y)*R2#hLlWbz%#7cjM9auP9OS_hy%{3ht06iquHW5T
z0I!GL?8LH%L~|98a3yOYh)TI>m|K74n?91h=TlTT{hk{GaJWp~JP%r|o_%?q7n*z1
zd$xOVIi$ig<bTJ${cTHsLTeLt8vh6(gBwa5zdF~2fZx3aE<L$I8BU2e99L660!bh{
z!@a;fSUxYBzQu=+YF*o5llX!L@poxVw!9VJez;{4jjYj{50LB_UgxC02Zgy1E-NV6
zo}gSV9VE%zyw8p14ULx8j>baf4}T`C&YhRJ)pKmJ7=;I2`d`)O(qDy>@B97T?Mjgg
znuwGidt%s}ai|4;545gsB6NLkorYg-sP7g}(k-M}A~_2cOG>LcSR<SA#-<gHxq(bi
znEJe_5P-01m1&ZCr|Zu-04{A68*Kt!A8#R>y)o4UmeF@G@_`1LzLQ^zg=^1Y^Mbq{
z@3AtiNen9{OPo_%j#qJ$;3iU+vPT6ZCp_(U&xqj<&!~`ro|~fEO2={sgJKfQXUk)o
zMCi!IYtP|zNAg&Wh-J6rn~JtSP1JqHV&jLFV;gmFUn>x_(uuPkbliW&TJTAUjwIe4
z@*Sv59?nt%rRav6!Dh-w4?0fhPfs0(3oZ~4Ic;&!{Mnx|klYASI2rHSr_CL1@(1|f
znxqDIyX)1A5$-}749QL|j0v1ZUm<SITy9Pa?rID5%z>%=L$GuAk_Z`yB{`95`>U$}
znR}EBYK;9|Qe{3z&25gtA)zUvZ0mlk)Hi|63{w|}lEWLMSAQDR%iH9E3H5)>N5Imp
z)H_PZNlIyap<x1`%?zuvC3YG~0%#b&G*07fQ>K}9U+&r{TkLES?8deY&8>So0A)B1
zH-!ObQsKA;yU?p(=LN|ZY4a6l56$Hn!<7nw3+j8gu8li;#`IP0{+FEh>hIe)p&c^C
zjY?Y95B&((ULE==%)6#JRvtFr|0@F%I_Ph)b>KrBuZXIfk@V4Qy>FuhTCPd_Iq&n~
zG(0X!cVVPiHGiu}Z;mI%%Ix&dH=aH7hZ7AAZaf*^^ml*X{`iSZ>m+y=@Mj$T@UI3m
z>@2YN!##+hD=~g#%^v^V<!c<-1w#qwuZ$?2<~cC7?bPKL*sxP&Peb$XIUW+AhM}US
zPk;pcV*KA1)HqT(oQPaHM(HZ7W#~vfK?P{vgc<5%`np|qv^rdKqFbuRu!Ls&4~9+L
zh^a^e*4=DHIk;Z6`pqw=1B9)$vB-s%=&J}?eSL>x497(yv&<qBY_X~IjSmAj+o~JZ
zE{NoXo=Q(lKb+QwWsnST+oiaoq4o}clBn~?PscB98(W5IexC^$ZvP<*Cm(T=mk{2b
z<3}(B!sH}$QE9oM#%p^AB6^foB4$uvCHypCoI!_h127oN{I^|ddsA;5XJqm$mQmC7
zQd-yCv)N!+5qCheFuvtr?Ico@9nt*lB2E^kD;}8aN^a>@7hlVe%Rb@IQfaiYf-mQG
zbe=B7TpLaImPDLp+bnsPo6oXwFQwGA#q7G`qocmzP(@7T5kkDUI*Io49w$BOY2Z!-
zDZ=ga%O#~uXHmao=G>~bKfb^tBkUX=Zs1sIufX<oSZK9!>v!9B+#v%)F{RAf&=uQQ
zr4g>6X@)O|0pZ%_aiZ+o^Os=5L|ucQik-9c_nDkFK1gNf^L*PLe9c;4*ne19Ll6S5
zOMNq9O-W@%-+z4<{{dQgb%oPTzcI@b3!^+3$DXdJz$9@0WLcuKjEn>~*X>L}E&41*
zgUlBvWkWTiYfz5Gni2E<F?jY-)jxRFE={7^p>~N<OD9|-OQ+|bQo@H(;cE%Rv9EpH
zK$%x`*I*L~46d730`l6tisvuQrrt>?&(t{_j?$vfg?K^J<B#_p{EVw%XSwD(;jSE<
z?xp)~Hi@EU2UwhDEm$$vCt3WVd44j~-}5YaDQu9bPh>+a!<2EYiQjjTZu9<?!HV(q
zkMz(>trE%!SB|6=c)HXikXJH{#^&3uQS#w}s^NJQK;1-m?#e|RI~GXTr_@6W0FvC%
zMWGnScC#dRg=)ABwgELvL-ox4<|pU~kq_w3H*-`qOdA0>cHptcOzL2OXEuXMIa@q6
z0&-(gW$ck@Uw*3gn(vCjxGnUIWYc70=i*(|il$Y{K{8)zjrPXECH%5lV*b#~BjzR<
z>OCT6SLLmR-l0=)HSNn1|8UKhYg%7~Y;UmC2483&Zn$yN{C4^)s5OL8?2`>W-(_El
z*o#0+=<jBh6uXqE|9121CE32HkOuHA*`{=C^@@h#yD>iy>g1RB@_6X1sY7!ocHGL4
zgoF=E1g&Glu})ELeBTR#JVb~b+F4oECA$~>p*<DD0M63mo82@v$CpLstg%-or#{fQ
zBoa-T0Cy*A#%;{%<&*44HJy~zuolnAD`F%m-ViZyL@eE1CxEgVL`$^?V_)xwK>DvV
zqoTPm|76@_skV{7Laxa3ho@BVPLFSg?ZlGa+KU1gm1h^ArtrQ?Ax=hT`0fE?Cr8<L
zz%W3CyCjqKuj=<?zQ!8cY9myv%g1?@gY|tn+GdVj4y6c6lxST}NA(L!xBXK}$-Rfs
zv1rEE+cTtcA8rQs`q0Ba_^5E}Mq%@lqCS!AgzN|g^Nv+S;5#tBHhDq#&NqF}!(@VC
zF~{0+z2HY}`zG2o4y%LuM6L-900_0ix?gPi7|;Eaipze;2tp8^0KPfdy>CXc734?=
zrVdJS6OMG(8vosiA9en4{v?jFKN+fqO<lXGy6QRXgm3_RntU-KC{`nS8QcW-E?_Il
zGZ>R~)9dMndvzHw7<N6m<_-z`|G732JJuDp+^08NKYtIZ&EA{R2h{o$t6Gx0AzSbV
zykB%T>W#M?tW-7XGpv-ur@jikIyB|RzFvV|{Ve?V`isQB6jjIngW>;cr@B^(p54Qv
zU}tkHo%TZJrdmPyh|rc4t;?qW)5pEXGrj+D05_#EX$Py(1s%uFW@W3~bH9r)V&pnX
z<gOG$hQ6mTM50DwXr$0lE-{>?WFfcQavMg6PPEy<7{7JS<D6f=^Sk`}{(Rq`@Av!9
z_x=4o9^co~kNH4abA0`?&lJu1`RS{tBLn5!<6i}+Zt9<A>Xx&H<(J}K_{M?n&(#I0
z#S@x#hegUJx1_lGm5k?~*ni;O;Mnd>$HmGK*){mb@lSIrk}syW(_a{$;v{MJHC}dJ
zH}=!<I1LHyYwm^->@8KT6sW0CYRb=csw8vj%1u?HeRU)?R@GDZ7A+y=k&^n@Wrgy_
znoW%_OA|3Ae#89{L4J73%B~YSqLj_xm|$FH%VZ_7qOr9VH561wz8wVAJ;`^iA9?IL
zqSC=eJAL9f;8rraV@k-kYmmx1CO7w=C%0QuMZ=l=f>SO19>@1p11T^nu_idR>|tiv
z5s9^DVR_meK5u!A0NL#P@lQx#0OPLIYwv7R(sFM$pzhpxzmf~CDcyKbA~^@!)~Pam
zL4p>+4DP@5<Tuian)M6%YQr0Je~49)+|K#+=1uS(RGv?_MqO>MJ1+L~^pPwJf-U(~
z$mxA~qi^`Rm5DAWT%ua#$E)1j$Gb1K-U@3SP8T|eBx$OEZscO@g4#jpl51u=LUz?7
zr=sv4O_=n>1yL<n*wKWLO9|;ucg8)HXQ8Y@+wjB#qB29p&EJo^7**y7Yk+Lu^e-gv
zQm+^hu*6;g<6(8ph-Lz;Q#{S_v@&#`G}m044B?m6`GW@>OWd;Llt#SVmAj>9#udD?
z#aZSXH!st#M^gB&qBJ2hmhxPq=(9rR28tfN(CMg;f&_dx5G&?QDi#H$UXl70<&)i5
zPRvBJP2$@HQ9O-GX^zF?lafl$UTMMSlOsh1FNG9+pKExxPLrSm*SY5SL(?*Ny$;8x
zE>&&38_`ntsP8Mprd|ujwLRW5AHYhL9z}o`-m(?#-qwkDQ~DX>XI#*OlFsp*XD*$-
zZU2yI=M6)`blJJjBZ^vj<vxZU%wY605(s9Zi{JUS_@isRieeLWA6J?fRHR0q9lNq1
z$y;*kdz7DIb&+Y%N2tNBxAW452CgZREM(t@D|!!zZ_J?Hm2T^^c+)Q$dg4t>1Btbu
zIYpY~R0iMQV4y<wS>Df^ISO8Cof*KE%apYE4Ax^ngNNL|X3-kEqWKy=+97(`$6Ie9
znWsRrECDa>#-N4QioLLQD=lKLrCTS&E}YTnY7O~9JQj})9YBFrY!wWlv%yY<>6$mV
zT}3TOZ&I!X(arauP>2YxdQ(^@+#9UOmT=l|6Y`KI)zRy(NB>&XBKu~1*T?5c+~8xM
zd#YON4n0;JdKPPGu&nv`C1)Q7Q{(8^qbN6Hbdx%G3;d+W05rOD)>!Xsg4f*##CYN5
z01z90ND(8?I_&!01%F(%we>A|{y{Luw=j+EJBH6so=1xj!?!OUCixzEt@9&+xj8sp
zPn5nMZR+1H9{gofd|r{}Cmf>OX~H+uKjdlg!y{?SQ_U&2&mGlB@_bpSCp`KRT{ip>
z_r+eOPC@T}&Ytc?#lUc{yW1*LL*5=StyT1<A@;MyT~fgY3|}DfB|i49r!}smAIcNh
z4hQCa5B%y2((fQVxmF&1mLEPL^l=H<W5N?q4y*-GJ)CqEOwbd$+8;1qBu6W1sVP+M
zgR?9e%=CVI$(uEG0SB+3dE)cFo8I}}+lKd@p(_HkBZAu-@C!xGUmiR0!>u6!WALTp
zD`U{~7h{m|pNs+je6Yd7g`}-LJ(JS&AoC51uFQobflPBVT)-Ht6w{7%6Cj*S)lsXJ
z8W-j2akoMCFCsJdi6NRsXh~n)mBVLs%Aky--^@r`h9{3}b5e%DA1tpDF=%Yn9dWMl
zMxke1(YVPi`D1exF7X&!p7FNf$EBUYowp|p!JN`ug$)1gpD}1x>K!=OoPdCG#1|6y
z(XSR!qY~WLk)03Tc-?|><qG>5so|-=Pjy18?tg&{xCa-#{(ubL{RJ7wz<rt`g$UzC
zX@UnqkAapvKg4!mA?cdrs+04Bo&0dQ3PlX`o9`lvn49Y<v0<JDGV;Q~q;1owcESIq
zhu~b}IF|{!)?AhB2lyG%MxN``yXoDSK{JwO8|~(mL*4Mz@V06FbdJ4z+wkvmcO2kS
zwt{Fqep@{OA1tB#q+pWCz6|!{49?s0`(;#se5_L|0?L!vI9AoT!~7(s+4L=ZW(Q|a
z`$%&R^}XBJd|3WEd|9FOF4Sk-dp`Vkz<~EO&T#7MG_B_E_zD<gZPGuT`U)6;79R=#
zgOh&(23w)}vpxE{4)UyOYT#a&%Jy_LY4YthPx?gsmwV?WB~TQME%YbhlwZV;JY1TL
zp4wZ|3EXhIvz72A(Q6A=;Err#5FY`PLQ^L~aZ0h57kvF#haZ7O>>$8FM0TnD9-9(f
z`6?mdb?h>@WtY8Z$#<qjQh+z)4)VHvnQ@T{a1r*%vCPQG&M3O}OkN_pvs_h#RS=^G
zVbl<Lu<V!|2_#~*;Lb)hgJ!|H=D1UYwUr%#c4}BE8$00i`6Oy+mcjplRX~gwuXL|K
z4KW$&brF+>4Vo19tI`lgWJ^R>rTb7_QBO`{#FI*QFU0d&MG>X>?8SOjShh~tfPQIr
zB(-orq%=o`YgO)a%2P$%s3TA*_YPt3L1Khw>1@42OZKSxb>r_vSofl&L|9Gs4YiI1
z4G$gc3zDG4h3*EZA)+E;7=f4!MJ#EYvwu^9$x58Q#2P7Q2dUpMRt2yyb4M8y{1G4=
z_S!M(ej<BVFN3}l5{iw=qnkpO?4xq&F|h3VyEb?Tu!^)P#iY@}F#VEcNgxAhL&YS~
zm0^Akcg^q`fHrK-VJ(}k3*jJR&G8@rfU%KBDAjlAQV?t8kz#v8yf&~G=Ir2NhF1p!
zp>GEl3%nj60n>08HpQO+j>7b*%Ob!NSS@xfi+%;-g&j7+D+3;|gRpauM&zd=L6DIK
z`wcs6#}opRV0*AyHp~MM3q-6n^CV=_QR@^_2AF|eD-27dTS2IZO-rU4a01rtu$D~M
zh8Q?%*)ib|IKusFVDLPIaT0+;x)(2_fJIn{{aPCRG^E9`-i)aM7$C-q!gA?)kX58^
z;j%Jt115)!vSEq=RM-|}*#_rAh6KK}E(>E+D5)(ekyoX<Vv^@*rW{X7G`hDb*MGwL
z&$3lHoZ+EDUd=*_{J+Eg)=hw@BYVbj{Xbi)ttyc>q`Kmg=bB9~c~KhBy?-ymuxrRH
z_TDTuD2x4TCVM24{VJ3FG?RTllbxJt`&Y-&Y#ae0*i%&G|L^y&yLjPPsEzK;&-L%N
zR-0EMbELYWljoXDFL_e_UPfjL|FOpLpscI*#^(ArTdPefk(Z>pt|relntoHh2m9Y^
z-?XsFX(dg^vn>+-?S!Z*xRu)c9h%`cS5O+D)Fd<eq$ysq;G!ngwGC&}*TtNh4GJ{p
zZ3A^VVK06288g<%4%wW6Gl4(7&nig!(2zTxb^dh6u@A#rwp5p_k$T=>4#8djJ^6mH
zmDoG=4-rSSK$CjXlVYg6to`ZFQpsXbX|=8qM@i4!ej4#Zea<hoIcuKs%DJ!6m(hrO
QG_eg8GiV4rCJPe!CtxW-7ytkO

literal 0
HcmV?d00001

diff --git a/hw2/.Rhistory b/hw2/.Rhistory
new file mode 100644
index 0000000..0d1d0d5
--- /dev/null
+++ b/hw2/.Rhistory
@@ -0,0 +1,57 @@
+auto = read.table("auto.data",header=T,na.strings="?")
+auto
+attach(auto)
+horsepower
+plot(mpg ~ horsepower)
+fit = lm(mpg ~ horsepower)
+fit
+summary(fit)
+lines(fit)
+abline(fit)
+abline(fit,col="red")
+plot(mpg ~ horsepower,pch="x")
+abline(fit,col="red")
+abline(fit,col="red",size="2")
+abline(fit,col="red",lt="2")
+help(abline)
+abline(fit,col="red",lwd="2")
+abline(fit,col="red",lwd="4")
+summary(lm)
+summary(fit)
+help(predict)
+predict(fit,horsepower=98)
+predict(fit)
+predict(fit,98)
+predict(fit[98])
+help(predict)
+help(predict.lm)
+predict(lm(mpg ~ 98)
+predict(lm(mpg ~ horsepower)
+predict(lm(mpg ~ 98))
+pre
+predict(fit,data.frame(p=98(
+predict(fit,data.frame(p=98))
+predict(fit,data.frame(p=c(98)))
+predict(fit,data.frame(p=c(98)),interval="confidence")
+predict(fit,data.frame(horesepower=c(98)),interval="confidence")
+predict(fit,interval="confidence")
+fit
+predict(fit,data.frame(c(98))interval="confidence")
+predict(fit,data.frame(p=c(98))interval="confidence")
+predict(fit,data.frame(p=c(98)),interval="confidence")
+predict(fit,data.frame(p=c(98)),interval="confidence")
+names(fit)
+coef(fit)
+confint(fit)
+predict(fit,data.frame(horsepower=c(1,98)),interval="confidence")
+predict(fit,data.frame(horsepower=c(98)),interval="confidence")
+predict(fit,data.frame(horsepower=98),interval="confidence")
+predict(fit,data.frame(horsepower=98),interval="prediction")
+dev.print(pdf,"mpg_horsepower_regression.pdf")
+plot(fit)
+par(mfrow=c(2,2))
+plot(fit)
+dev.print(pdf,"fit_quality.pdf")
+save
+save()
+q()
diff --git a/hw2/answers b/hw2/answers
index 0228780..affeb44 100644
--- a/hw2/answers
+++ b/hw2/answers
@@ -17,19 +17,19 @@
 	which of these predictors have a strong relationship with sales
 	of this product.
 
-	TV marketing and radio marketing both have a strong relationship
-	to sales, according to their linear regression p-values, but
-	newspaper advertising does not appear to be effective, given
-	that the linear model does not account for much of the variation
-	in sales across that domain. We can conclude that cutting back
-	on newspaper advertising will likely have little effect on the
-	sales of the product, and that increasing TV and radio
-	advertising budgets likely will have an effect. Furthermore, we
-	can see that radio advertising spending has a stronger
-	relationship with sales, as the best-fit slope is significantly
-	more positive than the best fit for TV advertising spending, so
-	increasing the radio advertising budget will likely be more
-	effective.
+    TV marketing and radio marketing both have a strong relationship
+    to sales, according to their linear regression p-values, but
+    newspaper advertising does not appear to be effective, given
+    that the linear model does not account for much of the variation
+    in sales across that domain. We can conclude that cutting back
+    on newspaper advertising will likely have little effect on the
+    sales of the product, and that increasing TV and radio
+    advertising budgets likely will have an effect. Furthermore, we
+    can see that radio advertising spending has a stronger
+    relationship with sales, as the best-fit slope is significantly
+    more positive than the best fit for TV advertising spending, so
+    increasing the radio advertising budget will likely be more
+    effective.
 
 
 
@@ -40,47 +40,142 @@
    dollars). Suppose we use least squares to fit the model, and get
    β₀ = 50, β₁ = 20, β₂ = 0.07, β₃ = 35, β₄ = 0.01, β₅ = −10.
 
-	(a) Which answer is correct, and why?
-		i. For a fixed value of IQ and GPA, males earn more on
-		   average than females.
+    This is the model: ŷ = 50 + 20 X₁ + 0.07 X₂ + 35 X₃ + 0.01 X₄ +
+        -10 X₅
 
-		ii. For a fixed value of IQ and GPA, females earn more on
-		average than males.
+        For fixed IQ and GPA, we can infer that the starting salary
+        for a female sharing an IQ and GPA with her male counterpart
+        will make (35*1 - 10*(GPA*1)) more starting salary units
+        than her male counterpart. This means that at very low
+        GPAs(maybe this includes people who didn't attend school?),
+        males have a lower starting wage, and as GPA grows, males
+        make a larger starting salary from that point, overtaking
+        females at GPA=3.5. Therefore,
 
-		iii. For a fixed value of IQ and GPA, males earn more on
-		average than females provided that the GPA is high enough.
+    (a) Which answer is correct, and why? →  iii. For a fixed value
+        of IQ and GPA, males earn more on average than females
+        provided that the GPA is high enough.
+            
+            This one is correct.
 
-		iv. For a fixed value of IQ and GPA, females earn more on
-		average than males provided that the GPA is high enough.
+    (b) Predict the salary of a female with IQ of 110 and a GPA of
+    4.0.
 
-	(b) Predict the salary of a female with IQ of 110 and a GPA of
-	4.0.
+        ŷ = 50 + 20*4.0 + 0.07*110 + 35*1 + 0.01*(4.0*110) - 10*(4.0*1)
 
-	(c) True or false: Since the coefficient for the GPA/IQ
-	interaction term is very small, there is very little evidence of
-	an interaction effect. Justify your answer.
+        → ŷ = 137.1 salary units
 
+    (c) True or false: Since the coefficient for the GPA/IQ
+    interaction term is very small, there is very little evidence of
+    an interaction effect. Justify your answer.
 
+        False. There is still a noticeable effect because the
+        coefficient for IQ's effect alone is only 7 times greater
+        than the coefficient of the interaction term. So, this term
+        holds significant weight compared to the overall
+        response of the model to IQ.
 
 
 4. I collect a set of data (n = 100 observations) containing a
    single predictor and a quantitative response. I then fit a linear
    regression model to the data, as well as a separate cubic
-   regression, i.e. Y = β₀ + β₁ X + β₂ X² + β₃ X³ + .
+   regression, i.e. Y = β₀ + β₁ X + β₂ X² + β₃ X³ + ε.
 
 	(a) Suppose that the true relationship between X and Y is
-	linear, i.e. Y = β₀ + β₁ X + . Consider the training residual
+	linear, i.e. Y = β₀ + β₁ X + ε. Consider the training residual
 	sum of squares (RSS) for the linear regression, and also the
 	training RSS for the cubic regression. Would we expect one to be
 	lower than the other, would we expect them to be the same, or is
 	there not enough information to tell? Justify your answer.
 
+        For the training data, the cubic regression might return a
+        better RSS than the linear regression, but this would only
+        be because the cubic is fitting points that are varied
+        according to the ε random error. It also may not, depending
+        on how that random error expressed itself in this case.
+
 	(b) Answer (a) using test rather than training RSS.
 
+        For the test error, the RSS will almost certainly be greater
+        for the cubic model than the linear model, because the
+        random error ε will likely express itself in a way that is
+        inconsistent with the noise that the cubic model adopted
+        during its training. The linear model will be more likely to
+        have a lower RSS the more test data is used against the
+        models.
+
 	(c) Suppose that the true relationship between X and Y is not
 	linear, but we don’t know how far it is from linear. Consider
 	the training RSS for the linear regression, and also the
 	training RSS for the cubic regression. Would we expect one to be
 	lower than the other, would we expect them to be the same, or is
 	there not enough information to tell? Justify your answer. (d)
-	Answer (c) using test rather than training RSS.
\ No newline at end of file
+	Answer (c) using test rather than training RSS.
+
+        The cubic model will pick up more information because of its
+        additional degrees of freedom. If the true relationship is
+        more complex than linear, then the cubic model will likely
+        have a lower RSS over the linear model. If the model is less
+        complex than linear (E.G. perhaps it is just a constant
+        scalar relationship) then the linear model will still be
+        more likely to have a smaller RSS, because the cubic will
+        again pick up information from the ε noise that is not
+        inherent in the real relationship.
+
+
+
+
+
+8. This question involves the use of simple linear regression on the Auto
+data set.
+
+    (a) Use the lm() function to perform a simple linear regression
+    with mpg as the response and horsepower as the predictor. Use
+    the summary() function to print the results. Comment on the
+    output. For example:
+
+        There is definitely a correlation between horsepower and
+        mpg. The RSE is ~4.9, which is not insignificant and does
+        indicate that the response may not be truly linear, but it
+        is small enough relative to the mpg magnitude that it's
+        clear a relationship exists. The R² statistics corroborates
+        this by indicating (it has a small value at ~0.6) that a
+        large proportion of the mpg variability is explained by the
+        model. mpg has a negative correlation with horsepower,
+        indicated by the negative coefficient on the horsepower
+        factor.
+
+
+        For example, for a vehicle with 98 horsepower, one can
+        expect with 95% confidence that the mpg will be within 23.97
+        and 24.96, if the vehicles follow our model. However, after
+        incorporating the irreducible error, the prediction turns
+        out to be much less precise, with a 95% prediction interval
+        spanning 14.8 to 34.1. Some of this variability may also be
+        reduced by using a quadratic model, from visual inspection
+        of the plot.
+
+    (b) Plot the response and the predictor. Use the abline() function
+    to display the least squares regression line.
+
+        Attached.
+
+    (c) Use the plot() function to produce diagnostic plots of the least
+    squares regression fit. Comment on any problems you see with
+    the fit.
+
+        Attached. From these four plots it's clear there is a lot of
+        variability that remains unexplained by the linear model.
+        The standardized residuals plotted against the fitted values
+        shows clearly that the variability is strong, with values
+        consistenly lying outside 1 standardized residual unit, but
+        still within a tight range that doesn't extend past 3, which
+        is often considered an approximate threshold to indicate
+        values that aren't explained well by the model. There are
+        many points with high leverage, and these values have less
+        residual by default, of course, and in both of these graphs
+        we see a few points (323, 330) that are rearing their ugly
+        heads. These seems to be the bit of "uptick" toward the
+        higher end of the horsepower scale that would probably be
+        picked up by a quadratic fit.
+