From dba6373d152202ff0b0bdfd96c62a15e42a23100 Mon Sep 17 00:00:00 2001 From: Devan Goto Date: Thu, 29 Sep 2016 12:38:50 -0400 Subject: [PATCH] Hopefully I did this correctly. --- .gitignore | 4 ++++ Class 7 Files.RData | Bin 0 -> 10171 bytes Class 7 Instructions.Rmd | 19 +++++++++++++------ Data Frames.RData | Bin 0 -> 8518 bytes Rplot.pdf | Bin 0 -> 4617 bytes 5 files changed, 17 insertions(+), 6 deletions(-) create mode 100644 .gitignore create mode 100644 Class 7 Files.RData create mode 100644 Data Frames.RData create mode 100644 Rplot.pdf diff --git a/.gitignore b/.gitignore new file mode 100644 index 0000000..5b6a065 --- /dev/null +++ b/.gitignore @@ -0,0 +1,4 @@ +.Rproj.user +.Rhistory +.RData +.Ruserdata diff --git a/Class 7 Files.RData b/Class 7 Files.RData new file mode 100644 index 0000000000000000000000000000000000000000..8f1958a66ff8cada8f6ce9a88594b44cae686c8b GIT binary patch literal 10171 zcmV;sCq&pEiwFP!000001MOW0d`{Q%Z6rA(Oa4=Qy)@5P2>cZVrj?4G#+G6B|0z926TGZmvy|yl_$^N(hM!jf|*` z;1&RN4dS~Rq$FZ>BfaXt&08>;o43SQpcP)+yw!gQv@TPiO~)~9-Ha#+ z(~fxp?M@JA-vEIQ1`2etGNv=Gn9e;B=xPe4@7GA_<}_~JPs0Sd7c9`D7MPwY1bV>< z^eT&+&!q|Eyh|Y0r2=_HVe;=KJ^xRN2cyC$)iyy6|rk2xRG?UVI)=0-{vW zGaG)>b6)boGs|;JOBKc2+OIdI=XFeSWgtwS>Ad8)_J}+YRY%kikvG{6svv4^nc*Dg zi^zzm45C_ynj!K)REJ+gh-x5ejHoOkH21ltq<$@;_sOEGiKsOqenV-&uXfTtzelMN zwM6tDq5wqY5LH8@NAw}0Hi&%47NbDa3Q=K1^$|5eR2oqOM8y&DtFs-V3WyAd+VZO( zQF%o5$j(**kq*&Eh)N;SAPPj}j;Jo8iu_K(ZwiPiBl0I}xhf(*M2!$tLZm|U0qIu< z*eXq440pA{ma|q5+h1$_c9kh7*B&QRAH;hz>2hIf@p9u} z+mRbjUdJI#@==Yw@f{)hI3>-Ho1=!l@S?}n%tYu|1D^eg%g`(68VCBhXnm&Y}nRvT<-b*X?eYE`Fwj%ESI}4vgxwtHIKA8zUDBF;2;MIIztXO zlEX2yF93ZguNiQ~=*Tz#?aVx=7a0)}7Szvdrf8#p^pTjNQ3F1U3GI_JTbaQ_LgRy? z%+Yw1(a*#u5 z1p?$@IST~HLxZ~$EOqgaA`k2F5F`)vS?YVi+1vSdw*&J*1Pj`!029my4nv!jISc{w zkS7nt^6)SZ^YSn*4*~mf7;5DqU>-8%p;aD^<>65tidAqJ3g#hF9^U04Rvr%KVOJj3 zHQ3wvx3}X2^T8#b$SI%5A)m-GpU63XD90_MP?&kBISeeSF8Oz(%KbK2s2{p8MMiU$ zyOia+G*%WrlGCH*^i~=DPNFwMemPu?(H9}VL|=q@*mBnJ#=P27tl@mN;)b^3igICV zIHI-j*xH_4JJ!a-#@pI|wtf=Ms~sk<_`a;3&*rz*pERaj7Pm~xV|Aifl&yaQwZcX+d}0pv}f?PEUy{*u98I>7$VUp88B*p|=R z=lP2TZ+j2q@CC(!#yn$MrwvZ@rd(Vb-$izX&Um?h2lv5izmRL-?Cb z5QdN}+kzwax4rdph%fIBcoMn(#mLpsRu?;KTcIYeV}7#-Vt$X=7T-ZH7_db!p|PLk zgzOwEWGRUUU}#Qu)st-%;Ob>PX^u_kcWWy6N5h`_)~%R!t=qN zhZZ=9KP`ADX?KDJ60$%@ycZujNqEn^c~?{TapnsYWvK(zds5p4x5wAdj{((#>ceG8 zKY7nZOIiLDDl z$}ed9d_Is}hISd+WoVb7U50A-x&-Mhf*Hi}`5Dn^c(jQBdV@cMESeUFzr66^1hmPM z=v2iy{Iv$J5AhfX|Fs7HwFmzt3XhBPUyATwcJN<&@CX$DHXqTZ~hISd+WoVJ16Pzvi<(BYYru=fb_%Bm_xs&{tDZgB0{>zk4 zZb1KK;_Wi7teyYg?a8&zy=8f}eKsuD zj;-rDuky+3d}qhWa%=C?TL0wE1&+zjq|5HhoNYX>KXPGfdwKow2wQ(7f^sY8n63TK zcV|~?@oYO*`px#!){jQz-lMDoqs)AxAFH$Rav!O)>G07z`y&BtJX_yA<|BFFa^)bM z{kR{(Y&=^&jp=_D6Hd!xe3li?$}2PXbok4wlQR8r0;m@_{K!D2m3-i+3F&MY_IwGPKLk zEHw)@)I*dALCJ0IoJ zt<}rPFEq^dUrruuc{w?=>DJztTsd3kr?r0PwVbW< z5I8$dx%TAR=iahB+ddnXYsc1gomctfb-uIXWVyBXX{~>9=K{y%XVUGycDdJQYya)O zc2@GsjNk5S$Ar`J7@uW@v+~N!Jstk+zIOjYUptRsK_cZ%_>I2x&NX~D@xRSC&1>9s z;7r-FqMtv6-^R1J?GdNFG8mw=^b+^QpU44`VWEA3dYgxuqmfD( z1NvuF;Xz@2LSm_ugJNRMQBaXP2?nV{!bw!=@SxbpSftt_^)@m( z)EpBO6NMC7XgMY*Bqq#kjuJOf7)8c~41%Bs35v5oTtsL@e`wq-1`VvQxT`TieNol$ z#!+ZUv{>y%?SzH)9{_!EgBcWKi&iapSH;M?06?V(c~73?9lDYC2k)vdey6#rMx#>EfL38tLs1&5QY-Y3-(98AD*coSjX|NO=~^WXII4tby;@;} z^`$k!qqE*x z1n3D;P*9<%j3S|_j0PIeDfMa^&>LtrXQfdiNJFF2YpD>mTBD!Bpi?W9g0xg>Di*cY zsG}m%tEq5}@J^hKY8vZoM7?1pr^=J$iV=wvB~o)EQgI+t4z>nY5?6#svM?vbtctOc zb5SaUiQuj@sL;@Ldc8=R2E7WT>8w(TCYC{|69PsB9kPpBr%{81oYV&Tl5~2}JX6zW zC%pmQkdsj(5tjO`)@dOk2u~s$Cn7XwB1m^4Fh?R_H!QH+*QsclDPbBoD(NIp8fYw% zpx)r8&=^$)Is?$!7G|JPuK>}C*04^Y(2C@$rxQ=DQfX;LYK?}@UadkyC!}bmIvETq zR)Ctf5@1mx6lWs5B4kFn6M=#7V6$`K1>ebUZfFw&tSp=;> ztsk1Kdb+k%YL!8dq1K>8R$|UZR#+Ioa<}obFPv#Xbpm7MTd!6p%Zdx6m<1EtCVVxBY%Xb{P;ae z*fv!fqtYnIS*cVA0ToSgp#3wb=s81WsC(ONLETZICLZnne?=@XfdmF2D+glNf`Z%Mzvn26A4h*UD4`R($46N zdO?Ody^fWk{;dvNMhh!yOHPFm5=yIDudX0h3Aaqdr8MQjSaES0*9i1ph z*rL-ytD>`9C+x3iN2NQ0-e7=dbXd^2XGE`!KI`?8L&W&j>$k9IeZmz6t}E`GaEpe^ zpX_Z}Btmo)2y2iYogxXMJxW0(h}N&4QEgCZL~Bs5qr+7hh4qMzSfOf-g38B5qtpod zrAB8Ejv67>AZ*cS3sy5CR4d;490n1B(g7lbiwoS9JvqA6t%(o~O<1IA;m)trYDEiH zVGzy*l~L`7&H=qzSi=U<-5M=YVGq!t1ce|MIw>8fY_%F1>8#ZoMAK6z(nYlFS**^R z(&0sf27+JM+S#2|Arl2p9yt z+DR=Ox=1~vhX_5Sk=8~xX)3q8X>G86;oj^?XMQg@aHjJR^y!>6CJ$vUNwVUE#;CW08id7jK1nk-X`YL>jFCC!l#N=yM8+MmT8| z!tB&irO`>JHnIZDy9NiE1_%gF39-5hCyMx{BPT$o`gHfnO=mM`guRDfyTa1dD*X%! zw8Vvpr8R04Vo?eQi9tc{jgD&Jnk<~+PFkTEbbl94T@A$$(ZW@;;>)}828XU2C)}Gs z0WOBcWDz6WaI7i7QkgkYEtJvK3?j5eV-Ov{NaGg)3nEnn9UsCqhd+N&PNh+Fs}?Sc zE^576kgVFM7tRr_@QjwWl9l1xDQ}?cg7OAy8cq^7kqm8C-q3wTAv!y^SSn4gqNA%KouSnQ5Ex~6qHt)S}a6+*m7?M7OTL~?o8-Ut>& zr3Wg86DL3lp!NcUP>7Qr8@cR$I<*#FzdPEN6n=b<^}>c|px9jweNMRLsL`dG!kroc zC)TO;!aamuXLMQtu+|G&p$IZXlTC;=Y9!J#=1F>hHG@PzF$1g_kb00F+y-58nRmHi z%|@zO7E$VzXcnn8q7%fR)rwA8gIZ+7PzcbmMkUBwqZYtvop4ul5`l{bm6nRrAX18r zMirywB7G;%!h-7@*nt3;2f=_%$bn2lKwsF6P9iL^b|FGaMJJ~cby?Wm(1lsp_mG4H z2vDsT$x$Ubd38zww^mu^GE!1>qBSb0Wkx}Vqr$RDkH}%>Edvw0E8WGrs^}BFJanbX zbyvELbk)%(YH;aF54*1P@b5~G1@x#C9#z7lT6ok5k6PhTCp_wfM}zPvv?jnUT~$JB zDxo!5^M&<)tL+K*=6_v}uYg#UOOnb<=X^&>m6DB+06Q#9q0+DsS~emVV3LWDAaD{r zUcp94HZDm7yIK{Jdx}KkkmzudtF=TsHn0)wWm1{jLa^)wrCIL|xj*pAk;Kd7XT3_S z{jvGjdPV7$ew$LtBv|qt9EkapHMrd3&oSRU0rQmNIDJoJ96tUN=6hRVzV!j-I}$NJ za31Ghu>gm^-9qBKV*bqx%s1aA^*6%#zsttqDFaCRA7Os1XF`sCXRNmj&2)KNiiKYr zL-<(CuQtK;U0F)PJ~({M8`5$4tG*C$Kj(7vJGCy3KRXAHbJ_biybindn2u&( ze(E!veylC#Cs$&A;t}R4hcG`xn%cD=c0XZSkzKFBDmSAkan?tJPIekM`U9rb>R-Dw@2)Z}-+%Ag#Q1MbOTWtuIkr%1TK3>lR^dH6 zOe?mZj%a?}U|M`=;yuq7(Wa$;?T$ZQv5{%z>@5zpPpvmC+0(w*`77N`OV2JYy3S)0 z%0Fx7{rU>i@|ly$*J-=Ol-SzIbm4kG(}GJ++IFcJYFhO?Vxvdn8noY&$@@9zZyBg5 z`7RD8w!wT$VGBqe~S4SQty@@ zFyA`}^OSa&Cz1KOo0;Ep$bK`+fcY0>oxr}ednx8KTVOtite?rGeZ1V!zDCBeZ!hM< zr(r(#73nXT@2lU#d^2fhsUzmF-oI{UI`Yy`)eGY@euRG z!xo-0jm+Qae_+0P8!2}m^UY))ZoPo{D)PP;ll@mIQ}SokGh9DlJC=9l!i-ptdPxcsLizABEtautU!IOBRQ?#J=hd@#RC&Y>IG zI6Zw8=9ko@y|%bKJHJoI<9w$c;qvL^o^+m^oBPSUIWd;x>xFsRH6o`%m|s{(I5~Gy z$virX4;a%iKOBF00_JDPyi9i`avO>JcU+CbXUTg_C-da^^PK(ux63<&>|>Ykef0l( zzis{M!sf3ZY~|Lzs*&NgsiC;APj{O%(JOd&pMYPJyOdn>dSa_rUHtY;3>#Zie{jlOQ^f?ggt^}vMrMz^-gZ~1 z)?*e1-e~@!%4yFQ{n{@Wx^_aRE#t?Z{?UJBiLHkQ{p>u{w{!DZuliIhJ@>_|Klkr# z_*L?<@pWrO2X#E^eP`SCyHUrVubJ7?!L0NQR%eDU&K#BSZ6mkL+C3Z1`C)L8jM0@J zj;b{!vh4BnRWpi*eVW`hOm)^+x4!fA!FP9>rW7mxN!2W`#bFNuE7tfjtM>D<Eh?A#L&gU;GrcN?x5hW8s6#!{WOy9x>TxT%W6z%CoayXMuCM;U3%8#6O zR$m#i>{H{h6&rUzP=f9-nE>qp>{k-q(b*M0)$IP7rM>JLJ`0&~Ck5B(TN;&<}r#q{D zK6U!6b!YZCmU5jmYIEnF{yUsAzRLWhu&&yiUyUDJi>R)6rrg(Vx4zDI6-P`xHU4N^ z{25*Lh-&BJ_nVJanEmixj}4i{&y6kQQ>N!e&-jm9k9_+67Zr;S3p{YYM_B5j?-DPh zxQ_eFamUPa^Pk?G9#}2C=#%IrSJDr^e<)z$jrM=^3_Y;8bGtj&<0DfF9~m?9%SvrO zIKKK-i~EXYEtAeKnz-+2uVK9^cj%h6Z%$^|iRtT`4PU)Jpw-h^lQ*mycx|)axFzQ9 z=~|x-ixx%uyI*{?cGAmQcXSKd>6ZVn=0xbT^`;-6pJ>(iyDO^$k|y-(xHshT*J;0= ztFxr{$i?$7qy&Ur{`H36=?yCi^{rCJxz*{K)3^6uyU6p-)bgc&m=l{)?W6Q@<;Ra& zmKy(L+WScrBDy#nDOWXc$ogw*OD;Ik-u=_xSDXx5Ib*AK>_@da@9655_2KXyFRw(N zOk1*jRXv9e{?DT0s?Dgo$RqgnvQIp^$JWVwQ2k=>$T3~l3>dzt&FJypjtK75M5W!f zJSt9AJfQXL{#TA)*%n{s&H}~Y8xhG{p1tl|FR@+Ik2maJ+Iz!@*JTz7YJygVh}xUXS}W@QjBkEXLu}Sx1_fJ9g`$ z$}F7}mbxZ8acSMA6~{h0=aJ<+y57U!uq^+x;U3ZV`@V{Lq~2Qn>fTF+?5vLES1T^0 zO&?p+-)Cy@EJZ-r*dLr)r))m6#lPCJ$J2UUa<00yxWctj56$F&IRQS){Gw)b47;81 zmp0LRUbD--?hB7MES$J2?R3JxRX#aX|BpiQUQ<`>kPSskA38rVLHq6rb+>>j$q6whJv1 zymCp%h`Dj8X_^aJRhy+hTl3|%DwDtXt$1>J(qvWVdY4TvqSCj$r)}J#S@BD|=bRqr z^=$OKhdm=UbbdPJREfvo^_TwES)cl%%D2HcYJPg~H*Hg2*Do7yQ1?CM|7X)1(XGE8 zt?2pZ=hZK_X=<2ysaXbBWZbl$SI>L(c+<^UYa>=P9anAvS2pXmzsul(fu))`FE2bT zuEihm`#fIzjJn>XBG<7|;e>>{?G|=^_G7gyRlu#67wc3Cd@#vWZ`Z~xu?>S}%xE&X zcktpKpHD9O!Smo7J}qut_eogOvev-VORs-=HSqB0t<5%m{G@NKBZtdfp5M)3-uCz> zJk!_DU=&RVpnUiu$lErO?bFKhF<%H7I?LhFom z&ODKIc6-T(F223|uRd51?eaJ%>U*E*)i3UBXSnry&jsU-Jx?1vws`VrV1 z_d7MC-orh?D?XfXqImVs{_MW*r_@Rx-R@pgKfis??XKm%@;T#IyYDmqZz~UE4%*7PN)2N-SMt<^Y()KJxdL)JiE!!3a(c!{xR>jGF!$u-TUmqM=ie(>mBa?Ag!oZ=8L$) zMe1g@FO`^ibmPex%X;m4m`H!16@LPDS zROP?=j%Wa#5W^=avpq`&WBCyJ=jHUz&DNn6?y}zP9RtlnP1Z9zEPK;`GZ&&hZVZ zzS0@B>oVI-n(#&9tTRQ5elopZT!$y45=JQrtfHmrN!5Eb=U>^rZc;DBgAu!KwOv{uR9{x}3c4ot+idZ2Gy452_tn8TIP0*NYFM7kR$d zyzA|K#XFhL&UsjP>$5Z8YsQyuS^I4970JgRH>lrlQufy~cfJ?+ympx@>3e$Lt$gBM zb2n{RRHMiT=E(t>iRC*x27Wy^w${~?(_AAS=_b9})8J}&EZ-)QJ?Ym8rHe_-LOlihF2N4t?v1+x~4yy6hC}ygK7iIaMgTowi(>7X4TI|S3h=r z^TkTRaid)Z9xS`&#R1*v^P3N8N1ky}Dta!zbTnjIpy`D%z3a#a5r?v`zZm_j8FSZj zqyKUM^9iEZLHzH2e&}hx)LXmM+xQzkoF~=lz$yb)4zL;k3J<6@pu~W>0*VQ!Ako<< zY6U0|Xv+kp5MC|G57r>Cu7I@y(EpF2eS%H{+6(9@K_h8>fR+LJ1X?;lMw=jN!kbNCRhu>`Uci4unvK>1*{if4It<`tq;&&KvMyAB4{79Px}j2 zRj^Wm)eo#_V3h(Z4_HmW3PDhP+CHGpfZ_t$MNme`_`#YB)={uFg7prpVPIVXYYkXG zz?uN{0oo1dF`#S&-ITNsR%Ecsf|V1jhG2yQs}@*^!0G~44ET)}K-mGc26T_0svNuj zg7q4#!C+klYb99!z?ufuDX{hs>j@nXp!tSDfW0Q8oi=92coI{~cyU_A$GG+1}RS_;-j zu;zhv4E9rEy`k0+K#d7nFKHjVBEZW4toC392dgqzX~F6VRzzZzqx&ydjlc>6Xf?4C zaFY86@x}n}2JjXD>pNJp!8#1qR@kSBHBcB2u~vcg2T*ciosf(lygJL_Q zuu6lK7p$fkISq!?N(NRZK<9~7gk$$#@IC==4)BfuZv$ezr`B+=F2ntUSU+j|V4VYN z7gz;|HAXUi@Ol9+67b3(UJf)LSmD8{4OU`eb*1^hDhO6Kuv!r-5XZiM;%x)oE8q!VP({z75BH4QF~c?S$kP~*@95k4!r&Iz?+xl2(os+C|D7a+f>1dklY;$R)plH zT(BY}dO^X8kZ3^#D?+046|9I{en6=pMM!=Cu0Tad(Tjo=!QO)kRD?>3OcksM2@Y4V zA|$9|!HSSztOYAVf+!cP2nl{(up%VZf`Szxv5ge02+40|6{rZQ-K$_lu=k+%u008(~e02Z- literal 0 HcmV?d00001 diff --git a/Class 7 Instructions.Rmd b/Class 7 Instructions.Rmd index 5ae641a..41cadfc 100644 --- a/Class 7 Instructions.Rmd +++ b/Class 7 Instructions.Rmd @@ -56,9 +56,11 @@ The spread function requires the following input: instructor_data <- spread(data_long, variables, measure) ``` -##Now we have a workable instructor data set!The next step is to create a workable student data set. Upload the data "student_activity.csv". View your file once you have uploaded it and then draw on a piece of paper the structure that you want before you attempt to code it. Write the code you use in the chunk below. (Hint: you can do it in one step) +##Now we have a workable instructor data set!The next step is to create a workable student data set. Upload the data "student_activity.csv". View your file once you have uploaded it and then draw on a piece of paper the structure that you want before you attempt to code it. Write the code you use in the chunk below. (Hint: you can do it in one step) ... ```{r} +data_wide_2 <- read.table("student_activity_wide.csv", sep = ",", header = TRUE) +> student_data_1 <- spread(data_wide_2,variable, measure) ``` @@ -75,7 +77,7 @@ student_data_2 <- dplyr::filter(student_data, date == 20160204) Now subset the student_activity data frame to create a data frame that only includes students who have sat at table 4. Write your code in the following chunk: ```{r} - +student_data_3 <- dplyr::filter(student_data_1, table == 4) ``` ##Make a new variable @@ -89,7 +91,7 @@ instructor_data <- dplyr::mutate(instructor_data, total_sleep = s_deep + s_light Now, refering to the cheat sheet, create a data frame called "instructor_sleep" that contains ONLY the total_sleep variable. Write your code in the following code chunk: ```{r} - +instructor_data<- dplyr::mutate(instructor_data,instructor_sleep = total_sleep) ``` Now, we can combine several commands together to create a new variable that contains a grouping. The following code creates a weekly grouping variable called "week" in the instructor data set: @@ -100,7 +102,7 @@ instructor_data <- dplyr::mutate(instructor_data, week = dplyr::ntile(date, 3)) Create the same variables for the student data frame, write your code in the code chunk below: ```{r} - +student_data_1 <- dplyr::mutate(student_data_1, week = dplyr::ntile(date, 3)) ``` ##Sumaraizing @@ -117,7 +119,8 @@ student_data %>% dplyr::group_by(date) %>% dplyr::summarise(mean(motivation)) Create two new data sets using this method. One that sumarizes average motivation for students for each week (student_week) and another than sumarizes "m_active_time" for the instructor per week (instructor_week). Write your code in the following chunk: ```{r} - +student_data_4<- dplyr::mutate(student_data_1 %>% dplyr::group_by(week) %>% dplyr::summarise(mean(motivation))) +instructor_week<- dplyr::mutate(instructor_data_1 %>% dplyr::group_by(week) %>% dplyr::summarise(mean(motivation))) ``` ##Merging @@ -131,7 +134,11 @@ merge <- dplyr::full_join(instructor_week, student_week, "week") Visualize the relationship between these two variables (mean motivation and mean instructor activity) with the "plot" command and then run a Pearson correlation test (hint: cor.test()). Write the code for the these commands below: ```{r} - +plot(instructor_week,student_data_4) +plot(student_data_4,instructor_week) +x_motivation<-c(1.666667,1.851852,1.851852) +y_active_time<-c(6913.25,6240.286,5956.143) +cor.test(x_motivation,y_active_time) ``` Fnally save your markdown document and your plot to this folder and comit, push and pull your repo to submit. diff --git a/Data Frames.RData b/Data Frames.RData new file mode 100644 index 0000000000000000000000000000000000000000..1a1f371409a83c8ca8141082d776cce583dd8bf7 GIT binary patch literal 8518 zcmV-MA-UckiwFP!000001MOW0SQJOtUf__bfMUm3V;66C?{4oD5M-qp>6!gZ~^e-|VvI{iyjLeV!lgc_!b^ zx3e>|@6OCO-}{ZEyLWHSv9}|`FxHHXjTK{UgDw@;=s#ABf^k67H6n6kpdlpEFJMq) z(5RruZyCm>9MapOR2mhV&m?TjVHQ==GYRYdW)gNMFbS`7nS`_UOu{{XCgF4qW>F1S zCgJ>JW|3Y(&FN_A1Hhf;l%8lcn=rA8<RwK&c~2olxqGQWuoEqU3{8 zw-O9XFf75a1j7;x0Ye+GEw&};KPQ`#^q-SeN&3&px+MMQq$o-MH5H1>lrkDd8H^DO zw!F4PTgDI^8W|KC5@7TjY%my!*3qOK6=1qSG^N9aej5=qsBp20BSr>A`-K_82N^;l z1BMv<28D))8wN#^+qRJb1C0h*H*5zP10tXov~Pn021SO3qpN(TV25BezvW}5;HS5j zNp~#ieb>TzJ@bxwS-ocMF&6GulHIQ)lkSC|GrMo0eJWjLvXX`S^1Q4+&|kSU!$V#V z+?V&8k<}x|LH?W!2TS-;eGz%H_ptQ-@_bR>)765sgtw*d^seRIi+IkGzNqIc)gvE| zmhKl-FJqznyS8tsKSfD-{qKsaJkQ9arG8r&FElUjza>26&s)M-e!r;iW$C$h9iK(@ z`+c8#*Kr6^{+^cFv(!HGfzQj^m**|D^RDyyeerqU@h*Q)xh(2?7S%sXc7c(Vm)*Cx zE;H{Y_usWWOLoKTd3igQzKgtnMcA=-?Z2h|n6*>XKBb6u-&rqjr-**b`&m>y zmg=)u-opD+d1{gAx7bX7&w^fl|3RARmrfJ_*v|irFwtw)vMR!&qvaFL1u5hy#znAle981$bvZ}NvT z1-@0G!j*zTA|k^_nmjLhSkEsiD9}Jo>{Z0C0p>^T*%LNz*g$@lh z1`YBXhMAx7Qm)BlJDd~ zzM%v8e(Z>1XlW6L-kg%vC|QnuF92O60NN3d6a*-H zm{oSP zN>)UR<0&k>TBo592|NW?2j4`Y-?I!A(} z9YYB#id-2va<;5e90*RVmP4 zjUX9kJiWFNwD5&&bSg8j_z#+=r3er%1RNUznt}l7L;$lU06Sp7-abzy!;FPtV9nA& zz-p-$At-1)l`0*lr6T~%ZE*zZ1SLQ#nZp{TQZ2z%paYNRI5n+^SE=aeRV!6=KuSic zjaJLa0aSeufTam23Id)Z8BtCIFaQr0J6p!GQKA6S70;Gs1*HU5jb!!;bPOP5>03IY z8L0L|gH@n&o8vjHh@o1mhDxkOY*c8rQk*n`bS#S3TJc+=mD6kls^1$xh{pnm_7V`K zpeYF0wggy9Bha=9TGU|>DVnJ0^;?nEstWwNy^o|wnpO&L`hGMAsU~5}iGIOmZGuMt_=s>V_A^4UmY7pA6O3?@) zqbdTY3Bh-kpL*Ju4s0%v@;r=AY!Ny zG;$1u4>Gvc1TTo?!kC1KY{!UZ! zMoOKyV&T};9J#~VbqjM4&7kbUB#cLH$O{=vEDznikwQ}IJ?>+bw;)J9-6MF&< z3^drOlwrgbP`*8Vp|1ks3?mO4665ZDrHq2}l) z*NE#YT2bkWAZWF4jf@2ydph*l=(Ql284=wF&tE&q{DdtGY*(BZaTg7nKa0C%353WL zh;xvdP6>i&jZ#ts(fsw)@mfwLnS+9czAmQ|=OZ$)V$~`o#m82~s>JnDrO}F}M$~G> zB^oWkybMBK`N4J=AOxoaAcTzz?3G;@I@ODU5WSi>NqKS4&#Kju39HnKc7fCJp2!Xe zyf}xolD##Wq~aQ&LN}BmE_6^@QEb&Js#K^2tz>v=B)*7RZj1T%b2?lJXaI~2BRU=q zjOe8~GNRJ~W41^TCGMM&ZY*qA$+D=@=){ah9Fis*A?|oMUZ5}{=Y!)sb*62vbZ0}? zIR)}W=$z$Q4PEslR!PulD4aT0Zi>-=KFExK=OHB`6pXkJ1PE=JLeu{3W<)aFcyTdQ z>nMy)EP4q~w3!t+ai;6E;+$10C8k70SAT^GEP6wu(TEwX=vUkDqR~a}88sr*NF%R} zF47dY_w(9d{=&Z5h2SFH{sF-Q)&{mZMOqPSR=VAzGe)uiD#axqEr5c`Q>#&Nthja7 zh+G0PA*|G^6AiD_$xy4rO){E@;)Yox4o+mJso$rgz6u5K{X97-VPMaQjtIx!>V?XD zHLZkgQE4Rx7VnO1HAoC$^b4>OZm1K32UA4GbwB-d}M3&bHq8(9-*Jx_9atuGr zc>`w`oHv-$U`ZS#7{0UfhW0B;iTw~H=pk=Iw9UL|N?C-R$TFi1i>Dy4;x5=UC|NOs z%q~oJYYm;#8jV_Xbp@rUm$==?%P~Xl{k=EBMN#y?#js&SN&(znpb$zi=4Rt9tDlBf z!{>KG%aYR5)MG(h5Ve%ME9B3KyBr>Eswv&65_w_`FNpgPx}DLeMZ#JTy+R3PiUymg z)~U?U)4dmZz%>IR;Ftl|3`{+s2fIPrw~V{Dam_mFSvs-;i$)Qzl2{O}S}n1%T3!;w zP>R&CiWBix@gg~`5%-EV5^_4j$ zeUUG8W?)IZ3kg#yv79XGvbegT4YRoJA%sK1zSJzBOIEz( z#7kbhREd{r@lqpR3gV?!ycAm#=@ws3Y>g9JvlzZGi#OZ8VG4zp-#li_!^<93sY2mz#B!`$VMgr8X$npySE%KRw-A#ug&6~9=Ep1L3bTdFOd+4GO0#{6na5$~ z!+y5TreL{m8kQ+;`2PMDIG_4EmIpduxg!tD-HBKpzJSZGUX1fUZ6|tPEPtGf z<+j_T{^q#+&-pk%b0}%QHI^s)#}w*!_NIbY(_h({QXsEQBytj#S6kxxvR0712hLw} zhx_>XtM!R~GO16%GTQ;`uZ|)5v$!AGQCMCJ#PZ55EU$JU{Xc^9SqN0D?Oy`b*+Bc(dqYGUWV&e z{JSqYwOVuin)%zU8l2grU%tO;KKU-rCw9VeyIp}y8Iy$bOL}9u ziOi>?6|mev%I_lIef>W2JxG2Jd>?--=l8>MBB^)#4J;1~$1y;!W=ugr{pDnw&zef)Vj{_Ud6TSHw_EyS z>5regs>AZyFr3-|l-3|#&^u`Ad6>I>vm1D0uLu|BI;N`bswAL}ob$MUj( z>&;2e!O!<2`dV0@brt6?DsVlQ4q^Q@4=k?|J9IN2-_Kl&g}PE8`^24IVKl{XA#S?g1fzp~0RJ z216h&Z;fn*)jV{~#U)5I74^`JqKJGoe-q;O87iJ#+L+CuV{T;hMB1b+t{W2!? z)zoumVr~v=@WQpy=-MsHulqON@34;b?ypyckM5Rzp$i;P;}Ki?w<|4jC;s*G?#J_wm!Fl{vC7L4>Ahp8yZY}N`u zkC9UZw9Z$>WEm&HT!_ z`O(k2s>%+qjcj^5`AN<4bEC4m%v84WKAE|^=GJMmd?OP_6!@_O1b;1XRE$W|7$!u=V8A+wWiFP6Tk87e(Un~vE#S(?C-T(k^Ozn zH+GuZ3!dt}ycSYN`HVf-WuMUS=W1hTotb(fD*CJ@e{AjZ(T5Bts?L9KukYp@xAT)q zc~tDb#U=Xdj^mzuF{7H>XrIIP`x?`h{+xI*#eVX?*1PAOU-abe9G}{mr5}ec&&oXZ z#Zm9+H@p7bKj`qXo?Y%-j}A?-J3evTchx(8nY!+EyZg$O?UODnoqq7ifYAeLbo0$U zxFE-PYR;xMW7ci*?(if&Zu8n<*S2|1UT)}3@5EszZyPv(|YhO&ietK){tBGM=7qxP0KaWs5GRcpI-gz3G|0d3CA4nhg~l z((BIIHDtq5mpiknRJgGqGNpFw%*j=zj$fG;{do2lNmWC7Sskxj%V*@KYa7ZfKGoH! z-(Raw`>mO~Lp`Z=y`H;$9df@K)AvgM$vb;J zIsAU$>#&FXjyhKlT-N62cCWHdc`;+oq`F=nv;58Duf5xDN$g#v%!_6@ zi|ZYjmdj(=T5i=@R7Vzh-&B*o9GP8LEr9wc2DpTmRk8nsGB;xFu&M#c@3w zU(vq|%iQ^yx<$J-ZkP8hNT2NbY(m0={vn%tKACx@?4#f&D_-;z(q7j5$^T~Een(!Y zTh+J!uEl0P@Ql|#t!{>Q{9%H!|36deT_65eyneLxR{VyL)vYF1 zUd&X=z3pW?Vwg|)Hi}htv!mMm9evRGjmP-wy{a+Yo7=_2-0iZY=d<5x=W^b^zPi+~ zx=&uLzVY5I+asI#&7IpaZlM3NzEk2#fBD@1rboM9uY1I-Z(na%+T}ODzaDmM!j3ju zzJ457@A$FGR~GfQO4t?sxUV5;R_CB(wpVn|rWHCLxqq;EjX^Q7$qRn)_Bq+R`mYlX zuO7N$#E?#-&qTlW9K7$(`sGJ3$$1wevlFxJojc##HzfGyj4ibqJ=k`Bdi>I*jWho? zw)3CqzOvJsns;jq4{A6`k#j2d+^%vDZ0is3x|+8*-1d=Q*e@P)>Rj5>Mf>Yt{TEL@ z`8;F9B)8;6x`FXWBQ%q?hq$FI_+x@s+gUGWA8FE~R??mm_5tSw=U~-eiI>*YDC^dJ zb<8rKo9gCkrjtlm6u=>pY z@$y!6?`=GCVauGJ>aK&&%x(N&zyIp5rk!%D^X)%<4*s52z4h%rrG-UZ`|q-^{JqCn z&jx|dynd>FcuU9LnP30oK72)^VAuR6uX>#Es)>x%iS@5ehVIv1+_v{tc%6UdO*6jP z*{J)mBT>s{T(Ro5tLEq1ciwTYckA}z`~Az0s4>6giK_Nlm;O$8QE~fZn|t40Y~B7B zyF;x<&c2S`ko6+>o_uLX=+;QdK2zm z@ws$s(W6#PJ(rv;U*q4vac%Fb7cw(m*&JjK<|nLn|FP4W&ksKAzOTa2ANPkRKRNKd zv!+|wH5Y3q=Nc5L-9yXPid%nrX~R{9ZP>RXK+-v$KzwZN%Oc7TldYTDf!ARuWGKX za@c2Lmizu@U)Ff^_-XGS&$*}UUAJnhr{A8joAM^MV-`)zYcl3om&Svim+8yf3G<2O`4&@0$DJ)Nr&wsv&KGNY#FHm~eCq?&6r+tc^m^K*@D=A7S>SNrIiu-C_2 zUw#$7)aA3bzPAs$^)Q@U@W5`zv$MabrdDX*;GEm)%~0B7>b+(SZ~Z*_@|iI;C+}=@;d|fAXR*;^CN-@+tRhpp{;f_U2G_0i z?SwifuW!3l-9KuA?XV-2*1tTgNx!h|sCwL4TUObB)#Vcbvwie0b(y~7@xy+!CES|89EKo5YsP25+r z_909Fv^&t)gchddfVKo05a=vIQ_*rjqW~=dZa#6F&Dw|X1kn9Jp93AOB;|lE1^N-_ zJfO#*KA>-a4gq@M0)3yt?L$}u!VL4s`vVP5Xk|*%655l}h=i7-z5PpJi4unS_+yL}F(BVKY16@n#Pg)<)b3k_ieMIOO zXrJ~M!d?(Yg0KvPIUsBRG(6DSKob+%m9`JGAkb_;TM-(_?EN9!2H`6R2SIoR!WBUO z1Dy`^G|;_JAJCCN?*Uy#=qI!GAx+5Ox3>A82u)nc>f501XPXBG7b% z_G0AgFN7B%TnFJV2q!^!2ErW>J^(tN(A)HVfPMu!6X-!gx0$sMVND1VLf8$$SP&M1 zFbjk&APhigby^?LzCfb_ElFrTMvgzMV-TK%a36%vNH~gycWAhT-~)6%)CY7k(6>N` z5_-|BeF)=1SQNsHBy2~^L0AdGG!XWHFaneVnj2_ipkWEEY1Tf3k0Bfj;Y|n^Lii2B zSr8tAa0`SlU_B-DI;E?D{v~v(S^E&ShA=RMRUu3XVLu3?L0AgHJS1#F*I%IZfhGsq znb5e5+5Urs%OU&>;amuhLbwsacMuMPb((~0#P=cL1fb`E?k4oH+51D79>U%bMuxB~ zggGH>2w^w~Yh@RX=+H0*gav?RC$u#qUw`3v0EF8id=24X5?-a@N(leK{)2?4X!{U8 zfp7%S`-Cnxdw)2V0AYRzo0Bj!EeByz2s=U;kA%f&ISAW87zDx!gr;ZY-=7?Zfa3)a zu9t+rah`-{A>0Y!LrJ&-=SjE-!Y^>&LU;g9b4mW9LdoeP=2M(L`oIHbPZ;yS19nD@ zp&>(vGirmKQL7Hje{x1ko-KK{Vzo2+2%)FrDr7-jPikIRof8nM$Da`&tSFuuRiTD+OoHOcgJM z84p*y6lSbr@lu#^ti?-V#wZsrg&BWdycA~9g5sqxi;Wa7h1nl9Dpm^fxL5H~$oE12 z-xc`(Sl8hb9<{{!kAK7x>p%X{N__w0AE_k$`?yCbKj9HdtpE5&C$at`ADN{6`}jvC zasJ~Uk;L~u^3ljobtLi=9)-mBKmHL&QvXLk`iS!%|HvbG{^K5X{Demw@%@i{v=It_ zf+LNjUmy1<+(FQi8@}rH?KZWo<|3k5YGYt9v1B0zQ|LKze066jn ADF6Tf literal 0 HcmV?d00001 diff --git a/Rplot.pdf b/Rplot.pdf new file mode 100644 index 0000000000000000000000000000000000000000..d6c37114516e83166a305243420aa19492533420 GIT binary patch literal 4617 zcmZ`-c|4SB{}v%rkv)faC`I8hV+?~zOvt|PWEv06F=nP&vW_J}wm0p{8c{lukS*j$ zlx6HuQlyA1Av9FKXF8|OdC&R1&wuxIf4BSld9M5Oy>4Y&oQXP81C3CQ9*cex&5rK% zq#=-iCO|p%E5g740UFaFJcCFf}3y_((u&dL;@YKVKS&p2JB`_qxdoj5RJR5rJ<#Pgt70rC|pPtao_*~BKvZyK>e%+ z03Ac95C9tE8F&)K9{_Ff{tz9|f(`6}zZ`V_<-oP$TGJph<4<$Ym_la2OFE$Y-5s=t z=oBW60LMfBBOce_M>{MbvOnWEfYHGMpb3%0fDZtgkl^0nAOgh~;+9Qkz%2_z1ZUiE z^0M(jZeg)Y<_T%835%?whhFT;w>nl_E@-%kd6>7@Y)ko-X3UcsKa6CkwY`2?e!>ELFrW%-53_Kj#QejvX5G;{dh4~%KL$F!s2pGJoGE&BR)6C+&zaC!z4fjhl}!&O_)Sd1 zvZl3$MsF$j4~EFf>G-_iKRNRDUQi00&yX+P{r*Q@eUM)};wHuO%iZyyChIPwh3Hwe z4Q&fMcU-=Fa?YvllS$s>jC`oHc!OKw$XxCUS|Q^$rz-wH&$U1?7Dk_a)ZO~DXE+2f z%z9c~GlG1WcY5D9g!T1I=CQ6fYSdQLyS`Q&G^R=mQYM}CR7D#Dp5;yIM)WAow?V6x z&0^Hm`iuooHU&DGArc;}`u_8=+UoiSS5nEHXEv1wMQ(yVrPZM|FUcyi%D-7-=hHb` zZ6xbJm}Gzcuz%8x%cej6_8)ZngJqxv(U%T*09^A0GglxBsa)+I)jIU z8Kw>l?!{wZzzqkdBvEJ%R6OB30@CeGjR5T*=6_uPT!@>`13Y4EY=oynz5s@6<@MuA z;jVt1@er59WNs1=2c0AmAZ}358onmLKgo{L#$x_owsULx$@a`kL5>bXQt|9MojTid z+Y&Ey1su2MHQQQemZq^6DYVV};aMIf49<~?vr8-3gi?7Om1mq*s+4$ogI($LMj2e+ z395wbr6az@GOr$l%ynFQrLYnakX^@q=|9IBsNp%g9>f2;rK?cXmRwhf(TV1Slq7#2 z!3|sHdH6^?uc>8aMpCOLyy?8*gUO_H@m1-%8e!n`+TPxS$&EbUvM*G-DIs}muOv?_ ztPXj~NI%it?~e=%7f%{xv#lnQn%gynby9YYAqUfC+~h`S zJbvtTH;IrynR`n5gSuE<%|P|QdX>INuk~x^R!&~>j*Xm^!6~dZExSadZ{@jL096ce z;sc42JYLfqY-_*7cldxjiTCQ6!mJbWHR0a&kH*ye_uU>I_O#fL9iR$M&T+DZBaM%S zv+e6&&b`G4FP_&u$3Pl;>yZ6I3viEGJIAP>x4g@aUsVwN%9NK9UEfuo0A4tZYb&q> zEH71fw2CM2ryIo?u4Z3d6>0jKKU96`%=1;zu1l`GHxz#{jF9J1*((YNh57F5nqzH# zY|8h!BdI(6O6D5-)K`{EmF3q~sa=t&1XCQUO3beF!_{?nUNpdS*pP3^PDzvplV8{JK*vC@rC=25x| z8mB5BUZnDcZ<~#^X}bJc)D>|sx==J~o2cIPBu!DN-%X6tf_DjtD&9UTmd{Hxe129w zU!gkDEUmEBLqg#2VWh3)h0RYDBu|fBzm)5oi^(gp6SU2|d7#g!Pvuv{Cs8Ytfa@LI zos&EaH9eWkBn=anoNHY>Cd(#8mvlLj0x`LSWMu00NK3ihl4ZPQGG(&8so2vpkCX8Q zrLMo}4aiL=Zf&OBJw`&_M}%z)6ARmKzquuj;(XgjXw@WKBx2LS)NDl+@noSpv1Uel@#?wfx=eH}sw{#J?GpQxVEu@C zJ*q0s^@st^z(!8FUfIsfZx=>2dY59L;>?~*rI$)|!1*U8z9KXOH_H(t^6qTKS9==l za`w9J-Ka7z;k>U;g|XMuF1R4U1?{-x==s}~64%^d+zGRfx5T`kn4KwpWaoavF^}Q+ zKCY)6ON+l$&}la?p9;1 z0{$+0;*E03@tAT|$;8A=^K$bn^OUD2Bvp6h%69WjzF$Id@Ej&b-$~uc#;Ja4cR!|o zUr|Z_yPF|5KlYpVJ?}T^8@Z*`$Lj0ttGpF->s7BI@d`0EVEoq)-MY6wGK2$D+vH|5 z&SB#xh6iG`JJ37PF2`qk<+?HkdIlDHpY{~>=;eKKymn_MH?yeDcrA80U)ZU|$~33U zcp|@(5Qk_wpVoQPcL`a*&}M#ND%X6fiK)>FeH>~RdOCC_bY)41qxNOjmmQosj)Q=n zfT94uK(RpO){k*&O%hFBabMy}Qj$}SC1bWTW%i?-wS+Z0b;#QF0Y;b;I&OYFKC})A z!c&x$*2}=qp)Hs)9Sh`bq%11FV@EsxyV{O;^)4xEDWfETBw0(b-t*~ABq-IrM?FS; zvwE#swwsJ=vqxarR3+zlASILh@;Iv{?wI&-hvU73#oNW{>d)>!UViSm;Q6Zh_QXjx zi|kKsseC~Q>dpF9;aEds{j;wzF?lhqxgGeP?uW?k(t(cxNskePK)=i&BmEF;@kzaR zO9MlVZXd!2KTAs&SX{M8YRe0Gytumf%S1}`R5i_W+Uo{+OC_$dj%**~6_gqlv&JXF z63H?>-x6DR%fD#i5Zm8-z+0(Axa1khf7yIlmGbBXyY^97SNxdxTZb9?D;KuTP-TD6 zQ}zhu^effK^2oE_FyCZ1xNjJadKY!eu;t;!Xx5?p9R5T7heS?woI2dNUy*9&)4F`N zc5&;(=82%wr@ErLW^z}4Pj0$&nx&8`B(k|gC~`Yq+-_TULsX5=lT(-YOO)!A+LY#S z>9`};L-Qo_#IIk=JK7T5a-qeHP<+hgv3Fu`g7_n53*k-Pn`>i*(j2wnYN|owfn42D zj?$aMAxRo)1FjkMtE*zo${EV0V%=wuXVlNkVRJJCI+uE)Zd6zlR*_~kX2r7vDT0)9 zb5Fqun_bnrFYeiIe$r|r%OHI$%RX(v@u8IC=Bz?a;lPc0>&;d#vX@o&@2=a=cRB6K z2}d*00--fjh1Mb}zO210u`ItW zQ#zt6Vk?LNkt*dHYG8P^hM$*BnsfKn>!qI}rah7*6Aw4*V5&$rUnD$N8Pu&J&K49| zDOdz6OIs9M>@$Bwc}K2az1ZP$@hI#0xpL@kdoGFX-0nO&F#TeQFKDc0O~#=l908SHY3JFCUX`Qwl3D%k#$XikFYPwmS2v=}>&wf9|H1v;0@R z0tGAb^2Mh-N3i`Ddz5EVr{xd)5?o%{a-;A~Slir+LhpWOy64;AA)(*TTo=!_ewUqb z;39MC{rmA=LY+nET*jIG^kVP4~3sub&< z)PA0H<`hp<(#h;$H@Y;day+Z2E^+PZifjb>OR3isu>Zi+e!Y30#u$CVaOs|d(6^Oj~xd_J9YySA*A&44h9dY(BW}o99~l_)t(o@H$FUQeXS4=hAv_ z$+e?58wS7i%tCXF!KE)QG9`ZAmmk2x9h`qS_y;TacCU~)%Y0+B{w2Ktep zU;tzgNxm=?h$qk}WB@!ygSdvUC!Ron6*2ffDX)KjPFvP z3&Ik^8XpLk1gnGZefh)23?iL`7(@ad0siOUWb z|AlE`G5^YeMxp*@$DrYfB!~vT*8rdcF%*K&3xG$+eiT@par*(