From e1f1bab3d2530097e0771b8b154a6493cebd04af Mon Sep 17 00:00:00 2001 From: Lan Li Date: Fri, 19 Nov 2021 18:39:29 -0500 Subject: [PATCH 1/3] Lan updates --- answer codes_Lan.R | 31 +++++++++++++++++++++++++++++++ drug_release.txt | 9 +++++++++ 2 files changed, 40 insertions(+) create mode 100644 answer codes_Lan.R create mode 100644 drug_release.txt diff --git a/answer codes_Lan.R b/answer codes_Lan.R new file mode 100644 index 0000000..bf51a0b --- /dev/null +++ b/answer codes_Lan.R @@ -0,0 +1,31 @@ +setwd("Desktop/Biocomp_tutorial12/") +library(ggplot2) +rm(list = ls()) +#1. a scatter plot of those two variables (time and the amount of cumulative drug release) +a <- read.table("drug_release.txt", header = T, sep = "", stringsAsFactors = F) +head(a) +ggplot(a,aes(x = hours, y = cumulative_drug_release)) + + geom_point(size = 3) + + stat_smooth(se=FALSE) + + xlab("Time (hour)") + + ylab("cumulative drug release (ng)") + + theme_classic() + + theme(legend.title=element_blank()) + +# 2-1 a barplot of the means of the four populations +b <- read.table("data.txt", header = T, sep = ",", stringsAsFactors = F) +head(b) +ggplot(b, aes(x = region, y = observations)) + + stat_summary(fun = mean, geom = "bar") + + xlab("region") + + ylab("population") + + theme_classic() + +#2-2 a scatter plot of all of the observations +ggplot(b, aes(x = region, y = observations)) + + geom_jitter(size = 1)+ + theme_classic() + +#2-3 Do the bar and scatter plots tell you different stories? Why? +# Yes. Population for all groups seem to be the same but they are slightly different from scatter plot. +# The reason is that bar plots mainly show the mean values, but less direct for showing the data distribution; instead, scatter plots can directly show the data distribution including outliers, but are hard to see the mean values. diff --git a/drug_release.txt b/drug_release.txt new file mode 100644 index 0000000..5afddf9 --- /dev/null +++ b/drug_release.txt @@ -0,0 +1,9 @@ +hours cumulative_drug_release +0 0 +24 565.31 +48 990.06 +72 1330.44 +96 1642.72 +120 1870.39 +144 2081.28 +168 2224.73 \ No newline at end of file From 71f9a8f33c33247c8c93f167e3e3e375c26299c8 Mon Sep 17 00:00:00 2001 From: Lan Li Date: Fri, 19 Nov 2021 18:42:18 -0500 Subject: [PATCH 2/3] update for codes --- answer codes_Lan.R | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/answer codes_Lan.R b/answer codes_Lan.R index bf51a0b..e180e1d 100644 --- a/answer codes_Lan.R +++ b/answer codes_Lan.R @@ -27,5 +27,5 @@ ggplot(b, aes(x = region, y = observations)) + theme_classic() #2-3 Do the bar and scatter plots tell you different stories? Why? -# Yes. Population for all groups seem to be the same but they are slightly different from scatter plot. +# Yes. Population for all groups seem to be the same in the bar plot but they are slightly different from the scatter plot. # The reason is that bar plots mainly show the mean values, but less direct for showing the data distribution; instead, scatter plots can directly show the data distribution including outliers, but are hard to see the mean values. From bc7d73f292be78ea2b7f57652da5acc4d0683f7d Mon Sep 17 00:00:00 2001 From: Lan Li Date: Mon, 22 Nov 2021 11:30:20 -0500 Subject: [PATCH 3/3] Lanupdates --- .RData | Bin 0 -> 32717 bytes .Rhistory | 351 +++++++++++++++++++++++++++++++++++++++++++++ answer codes_Lan.R | 2 +- 3 files changed, 352 insertions(+), 1 deletion(-) create mode 100644 .RData create mode 100644 .Rhistory diff --git a/.RData b/.RData new file mode 100644 index 0000000000000000000000000000000000000000..a3a2f9e83b1afb3f93cacb7d79e562d379011195 GIT binary patch literal 32717 zcmb??Rag{`6YU1DbSRxmBS;8H#{yC!B_Pr$A=2HjyR;yUBE14iNl15>Al|BA1$$phAObA^0$g}jBdJJ)uQ z^iUg-^Qjth8&7)8bMl~#4npP|iCnA<1VL?xN@%9J0lMo~Pl~?aK!t9OoTJ~8y;v(r zL(7*St4m?`HojB|FOum8;&O0K2OFaHdbW^^xm}|7AA$KLV~sn*V+a#t!8-*T+9+?m zfT8&vnEwl7LL!*5ryJtTlK?6)6E9S|s& z%lEjq#US023(3Tc1!VW!2$JC?PT;cK-tbcya%_xf`fDLj=adSK0h1Rrv&m3Dxd?c@ zp6|WJrw47&d_Hn125@_=*L6F{40RO_k>+tm_~d>I5WNKG<{3L2HTA)Ew@C?(8p-Ni z?iR{ga+HEQjF>Y^FCPZONm&d^^h}>c*CdRy(=!K7fQp;gH;=XzBVo!?!n(&HU=C2&+ zPshzA)1m+rJI6Y~*9l(7emzUMD+3N2yPv_Gg<;c9LFEP%c`{NN<>t%SWFN-!PMeLb=V(tJIYNt}`` zgyuJDO0)+PzFt^Nyd)n|;^9n@1g+9uE;|{Jt4oyskjtJvjYd3`olDuc!Exoog6?td^SPb%pe_0c)SXbO=phIrJ|FCavve)R^5LAw98c0$gf z&@8vNSJWe9*Zpa9HPhar&ePA*9*m$(0SrDRC|?VjOUX|#<-1z#ePt|tLbKqFDft&1 zMSbF{``X6RR0Xn`Ioe?UmzJd=e(%vk{wIvC%mbL zC$LI#znX)I;P#_9qjINSq8I9AZEsuOLDOu?GrDPj-#4oCIGKsUei1D%VJr}{ixoXb zy5JS!jGHb|Er46=SN8H-2zmG?@vO2Tz_XX4L8%H)$ggDQi^fvOx=kkXB21_eqpKD@ z8;e9Qsxq%kviEU?&2nw?s(2En2hsh88VYw`zP3LZ1m-7d#{cg1LL5|JMqAZ^_*AmF zyY6UY=IhEdN(aa_fAXfqUokd$yvLn|oG^@g_x+g*z9yoK^DcMe(TWdp1}iC?m5iahGlMV+kN7xZr3#vFEI4rS(% zuKpI|>Z$5?u1QaS<)%F_cM*t?KG%7@N)J}suo4`hOoL)K`ue|}f_FCNpRvC)gvz~^ z={#9vjEW$Hrw!jBP&cL>rz~_)d=@j=+I-;Ev_ISN_Nq|P4a(QvJm4{v{^$Nx0tx(Z zEkjX90H)RI%;+T}0`Bsd*ZttJrAc`Ld+Q!hUA9$8lMC3d$y+QL-~)E6`CQR~AzR^i zFKJaT2pE3WV;UBRQSn+2EZoNJrm6gBPRrlkQejSF#*G`6U|&M9$405$c{J_n3{!I& zDU+2Xs;^V;%m_Y)%41GlJD{dC&P|0{4#1?P9=YHGFNA!)|FGmS?ur(Tf6XhX3}H;? zL80o2J4}$btN;|1+bF-OA`XhH^DY}z21POYi_rZrMF>jxY402flb^@NLw*<_%PO{} z%mPP{zQS|WijOv^AT*JlQM?J+0Sj$@XT+{psR%C89PrC;KYe5&(de@c_A~w}*jQux zUys(%T`k{swhsaca$K!;Y89}PuU9&2YhqB0pLEf~bG%jitBz}LO{g+`F~7_!KvwF7 zNb7eQm_OQi(kkXSAA^qHJ9xuGHWDgY&bhKTyw6CZRNtndP*4z#b4{&cv!qEGPbR!9 zs*GyZpe8bPql4_XsTT+;*qL!2m6iWx0*uH0H8?;}LVf1m#q6?xbg>UqJ$Wp^BX4z~ z!!z;r1!9|6_a4xk%-O%+k#C{dNwVnBV=A;Ii}x~?m;97@$iVbHNWoDvkzECHY_BXm zB2x;?%^03{ekr`7isLi$^&@C?_NA>#CS)K-cI{=p%qZ1Dpm6%tvlR^)QgP>Fl@Sf1 zORO}BZ}~p4<1KH3@YWYQ8G0E0L~Hf4b*cmdpH1-yXWUNN+p@r=WZ1e3Vc`Cc(A|N# zb)=k8yIdSRfOUqU|3_8G?n5N7?W2ZeD<@umgRJ<{MKEyZ$a*6qYXsqYm^?d8c83fk zVaYAn!Rb!z`i6`4gnGRDy=CwdaebZQpK2wZ;1{p8{Im^}KmY15&Tt17U1-JiX9u)u zSGYUHvV<>wR$i@I3Rtze&K7u2PF_h8FICRkJk&UqQguYwS)lw-^0BZrTJwgHLV0h& zK_Pq+jxS!W%0KwQv1X8><1r#PxpHtqc4bG}IJmk#`|S z`IZC=4ihwa@pJek*r=)>e*Kcsa!=0(F29X8OSa%=e(nx-J}Z27;KfcZyF9GNK0O;Fpjx?6X8+2ki?Sk55skn+6q*6 zcz4d}MSEPjkUKB`5b66fHxJG7QGZ1FNLLz+tT0DuN#g|xvYGbFvJqM1?)&jK-NAgj zuLIYH`P)^|S$l?x-t`7}M(7OO;-=G2H58C`50=90FO5L!qv5Da9ZEy=p(Q9>0Jz$9 z^{Zx&8s^A(?q{wdakXY0r))sH)jM~3@ z?VJjYboJ2LupZNdzBiRnr=CVIybAR3r$>!e%~4nC6<%5~q4*1Wbtknj2~cnz{y#^7RG zWj?RnUl*AKI6XzC_#Ub*edrdPN5hR@a)EcOzfJwNq{Nb)60eDDA^sMus03<&9pUZh zb{=+UUj$-PrdXVydvjGAwr_W&R*MgLzkkoY z^k9zYI`3vUf=S4=eHC%VJ-OS`mcHs21j*gI`b&uqEVT9E74OCYv<`dJ_JZ}>#%+ts zm@4qr#@Hn_qjC~R+yh$eK7p6lS_wjZ$@)tz&?SGrOZ9w7rTv{w>iciIccHq!dLb8ans3$8bT8TNa%=!~dJ$|Sb;G&x@Mcug-bI#uG<=%sK0nf*^yEop z_wDcjs%F7odSs~!c3$U~!Ww!U!28PYT_^2`LGO+l8+)JW)=Wt zF$k?N#a-KIqx+ditK~T}wz+%A1t9QQ{Vfrae9bzS<-iNP{*)@bm}e>o?Nc5JXKuVO zG}Gaw)O!>|?#j4O#6IRWrzq$P>56GCS1=zZ;K9>V;j8gHcrAU(&gSrKz|`EtSB8e32qZ z5n%Ny_{c@BX&clhI}DPusY`eB=S}!<{dJPgkfJTA)bMvFXy-QZUa{;k%$Sb4wb~34 zMKHlR5|jfVNKWTwK-^Y3j{cER#v{y$rlhtYX)LcUlbS35m~$FQ-kb(k=YRzF^qLU9 z`~icPK23;o=e2lK^7)Yx=ha8LAuzF0H$ieOz&ol)aUDg zA1_)Kv|?b_GzMt#2Xfcilgtl!u0!DJjR@bIo$rReO}=3< z4?S0R6moJJSx)N-xZ_TtqfA(1U5nB8S^~wmr)O7Y0Y9d}>UjcM;*_gDNR)R+(H^Lh zk_~JF3@7%kdVFuiZ`8{@as2j#@<5fhNN-CnVT;e6z~enaq~XK809J5CIBAIV38U)0 z6}SzBy9b(rZ|an z+v-R`kHhqv{Y)a)0QoZy)ACuhpztm(GbaZ*q#Qy6B)E<-?3ndWqyl?5;Dr- zDITtEj+NA$^al@v7Z>}c5yNZV(4&29+9mU9duFG;gUFi88s436*exOi4oz_+f5Q(Z2msZ_qXN5F8H9f>KMt`EpIP zbmCABWLD1PqR7Xk!y-!`j1%L#AW`=IuAChfH7=jG^Z7O@@1l+-TMl;(aUNeF6Gc{^ zbwB8P_DHP&8&Gr`1){tFPs=-K0yE~WkCy(yVaG;q;U_a1?cxuF&&_9uaXRmdx|j-G zmOXz=o6@h2D)=!!CAiW1aGhm6DT`7IZ`JOS!CQoPaO1URfNzwXXAjXI5dCl8Rm7JG-!z2!_wj8VK$Vs%rY@0t6lbd)`EkIS}`GWMn%4XYpHe2_GJ zq}CBV&*kHVEW8kVwm?iC=SO_!&u+S&@MBO1-~$? zp5v{=^Rq|z!&`#}rfZ1p{Kmt(o2n?)#%-OpHzn%$!;u|H9 z3B!B(dy%KdZM+Q)ENqCnI5DTs3PDw8Tc8@%P`m7GW2@9b*c=6qi-{h7et36ZPWrPb zA6`PVX%S?f$la3hH);Ph-#<>hDVR{4Q|lQ)Yz|*A{6_@|zUt3BgK;aQ-S?7Q5AezM zFo{IpheEks;m1b&sXF)7(ZS{ZE%C&+*p6 zPqD*X3A#-4(hIKikg(Pez*R<6PICr^$)ET{)2D62ORNmxfaYjE6!0HuophC}dvx5& zb@AyDWqZ6kOX0u^8K-+LQm@=cl%7m5>%1TF)_MOWc#z{KwyR^xMmC zQ!^G-f?AL6*jZT#V9rOunTTETOXZ>$SzU7a5R>Jmggijg#6frYuPYe2P|=ZgJt*vE z#kxiWhze-^cI0}Evl0hiG%#u)zTV9c?k>SE=o`3BXVYzo+TMz#T?jkYup4pdi9s!l z{mqbgOMnB<=sHlFh$tNlO=C(7dEbWGf#j}jp6xp*#PZy->@ zN!KBwRv=~}fFe_l8D`G(X=#9p6vqUTck@%Ju&m{=q4PDtIyg(1TAdiWD+w2_CV`Ce ze=1-pzycp8nJ?z?Q&sC!W`GZ&Qn-739 zW`e;^B&)r?elhV>yC!cKv4k^=08{ssd`AL#OByX4tp;CU(a_pOr8WMxinjNpA=~n% zt}eB{(hG^=FbADFRYPcJ48M4;btjOBErPy1&rQ6(NFMnP0k8`C$PsEUF7a=s| z;=&C=b$fTDzZM?QriOdBcX`|+y4KV zsb#QwiPfSp3fjmusV~=cIFNNDOsgmRHB|lZd#^ENo}A4D*)qJvZaFiOt1=U2ny8z| zA|K)^%9G!83Qeg0%5lg03i*pYVCh`}(~2|0qI`YW@wh>$#3b`nPkao;i;hRo=2vf{ zMfdS4KfiwuqS`1ZWcQNr21GSE=pEN=7f+Zb-f|Q7Z%Tp7zv#5pCrvP|Vi=Bx8?K0w z_7LaA=uRguV7~zoE7jE`x%RV9b7^pV{A5`s?kQV9fW>8j5^zmnb?|w|i8@ROlh|{m zKpds#%xP{B3_7fH`jPsF0V*-{lQDz+yIf!@6D#CmM^;QQv^W~Up>HW`jC2Lwy*B;V z!9kwF%Os(vsN2#JU+zZ$`m-AL_>wbv=4Kl*|1g9>_2Bqri+FbYEyZpv{EhWI}!B zH=X{BiBo)jWS+#G*l{;z>M<1ZB>407Ye2U`O>56lL}9$=GPm>-!m!#QC*I9MZ2WM8 zLq^!;!M!n>?T86_6uqPO`~=Vop^fX=Qd zhjDPB#a`pJsiP3Hn!G;;?^9iDVo*SuC;5i%YCcqn;RKlePpAV3$=@O}AF0bbqUKdC zpe_NxllUHrP+}_^TFIuL$#XRPg_!D6Sx&-Pd*k~<`_UgsFD42P<1h5voU?~C& z#_efzbkgR`fYXi`){obr2r{9zEjz^SH0Phnt;JW+V8PVql72W96o>GP%iPxQj*be^ z#u2_|f%UKeGWsjC*=>jJHc^t6O@Tb$+E>DD?^>wvDKsZcJ-~>q`14|OZ{X_ztCpHH zdg!&}qvv8~xSjn!i@o)+C)_fsjLqvv4n`Bn>H?hOt2V>anI*i%64x*AJ~`lO9Nlwc zyc=o5{0beM3#xdHH-+U83XA*9EiO5NMvrIek*zP;mNQE|_kKwvjKKf!1jvTcG!Y%fnJSJ?4kT&)sK^p)BKLv3wT+m*hQclmdj8?C!y z$7AxT#G%HLp~AT4!iz0RMa)U=M2mmM9fHMfTg!tTI0rKt=+LBEb-4_xI!c4tANa1% z2C(2Kq|ZJ;ad}RqwWs6@Kt;ZH6Cz`U7PSG`JDnxEb~DLdk$ReoydNM3V|{t(X04m< z4=vzHY1id^Xki?m&wy!t%x+}&+Hu)*K$PHfOyeV8klL7e@w3DXgwEhTW4$6yl773x zQ5hs^gnuZ*Sqf*}WKWQ>ze+JkVMoK+3J!h6hc;PYB5eeF?p&KPwKD_`k$vU{DulB)Xrf8O0|73444?;gVm(Y{#DVQfAiYT@u}o4b;qfP~F) zQxV?ApsFUnt`U@FPH8SnKP_1t;rT%q>StoVySNPThugQuj_n3VM08QuY+g3IR2f_Y7eQ0f_J zC~-TrTkz8IcjB43ZXf%@l^ClL{~hkA4%*zBpY;}xyR;U3EB!|)8`sVZoYXuQ$b=%im=Y`h3!C`$ zD>i3=sJ~95+)auh)H{n$SvQdkpFZ1A`PB_aLq|V}az?CMOWjGFlejy0;hSwr3_WIFEGXMa#Z6bP~}mTY+N_=$3}eC5>44d7w-*|%d8SAbip zv751+WDdP`PW%uJEjLnsaIGFV|4^}#aJA8_Cd7K4DX$L{1SowgO zqW!b*GX)p{pL#SeABY$F_*CPSlIkc)uKLgm)=Z3uu1UQvj;EMi!SP0K+Vt#%l)O_!CMAtxn=$I{)K4C-Fbn-aDheWh<7&>rAE7X#oWE<(?LE57Ukks_1D$t^`}MTT zZ11>c{`=wsI2hy0gTIzh2*5r!dPz;vPxs#iDHezd-im7fb*+gq=T3Uu3_?0K+z;rY z33WatzT8;f2HexQ^EC;VgPzUa3|4E2^`B9G!{}o1^f`k|deLHqWkEAME~#MmZS}af zhs^C}+9_#&WB?H?!O2bapB0uA___tB!TPyk`|=I31dS$-Csks^z=ut6khnRan*ZyC zbVFVw=JdfPK4h>k{{C2O^Vu#KE3y4xjF?rad}FdlsM8DAfPL6jVY7Rp-N%8yLQRvA!QNeWf)x6w`t(cLxhazL zRQ(AHV}DS#qr~+6Gd$Wxg^oWYKES~W!F+Q+{BEx-veaVLD4>fbz;ya+ThWw^)$r^$>8^?8|A@rip?)!Q-pbWZrGjkBgYv zb}-JOxzb!K6AEsbbx^){fcjwTW?TlDid9!&SQMApEqPu?`>!@Mc-}MCw1I#;c+$?b zFz-S0ZhDR-kPtN-c;d>#LwJ6_m7}i_jZon2TR3spQn+6rp%F16nY{XZ{ZBPk2W|J1 z&cl;TC+z7E{97x`*R7y&yNW=kVNIhonhtTfqm*^!F9>*cpB1T<%3~KinL5^zhe!X$ z$>y{q0LJ~|xzpWdDC{y2=Ig{vTA8`5?(-uK>^AJas53_;Y8;$KI~J*6wU&41?(lqD zM)dQHXS`ydbEiT|FHVbbJYM=86^a1AzxREBLb<$~Eoc8D{rf70%qcWt&F#y`YBd#9 z+vWNkuMt_x#g<^?&obE~t{+Xb1St*P)QivCgNa|gh!}ZCp%Y8n*-$lHzAX95#G?=y zH~5osIrGWrjU3mz6N~3U>qJ`KH~xU>%?|hAEJz`#{c`h~E0L&+=6(bH72Rsrrw@{% z@KFwbNER)hf*+Ihco$N++3H9&yxRxGkHBs;sIy#HJ#MLyo@h__70@|&A!FBOUOmP+ z4k)2ZW>|Bj<*3DXck%sv`z^Dc}g}P*Gy^O0oP;bq zKb=QNh9+N7B6MwUtHq6up9+;;-OYUEjz@Ci#q@XbwT_NnQ_7+!5j}fJ`PB}p^(4$= zMquI&b^lJmTJjxB|Gk#~os)(Lt(yGi(#$sT?c?XqvY6&!vcChQ_ohkc6&1i5d1@GoJr z8!M>FxlTt8gb%j{*%!hG3@e|=7ZO@8CurM#PzE${)LA*r(RkX%+o#cFEBao>4Bv1O zdKPH7_{$dKEM1DsluB^IRM$=4Pd~F5C6=hX3)&4xU?Sbn&2mt@(06GT*C!UkSNmz( zT!+^(d+K?bb|jE<&^C|X6L)eruUWRqS;+-@#bS_|;x zp(pMZO{UPc7dzH4h2P$n;w>|}4LN;V@Vg>T+Ai$iv7d1ugzY z4-CMJxG*Z-2tVJ8mL&Q%#DAd5_Y=Wbw>224LK(RbN9L&67Vqo15k2v<0L4($^8GAK ze7l?|>+D7vl%HFAFfieU2 z-{Y1(OU%^de>=LHs4MGXj5}HWr>vyc5UA;?q}ACG6U9ArBx*1Tn7>T_9a*Y}dQ|=0 zjs6N)*iD{z{7Dz(ogkEXl#84$Yf$~|%DtUTb;Z522zFolxMM@`97;^J?pOJ1{%pV%DkYn*b`1FvjyFwrlQ4L* zd8r2xKZ_CWhqb}Kwi#CT0;E_)qI@q7I!(}b#W49y5z9v!;9c&f@~y@Y1m?Z!nbw0M ziv9${q;rC6f2&$V!%bdpU?jD8F+n04fyLKe#s|s1bIZj1pg~=--_J*}7yEohYUCfQ zQGOUIN>dm@Fr6h`n%9vLqo+-cM%{5MkF@#(TFJvAg~yf*7hp$8*Ls25B<(Lti8Dp3 z@a%}hSxmn|^y0EZ`*PSG+&reswCsTWZJxhFR7hZ($-D9?XLwQfIr*eq)8MbB{+(Mc z5;|Vj?SAhc)5>`+okrvmO-6!a(@hsf15dn$bN(s>gaxYiUCSnP^^M2eQ;^AgnM_Y@ zAUY2p$o)re3=rqmx4*N+pZd|x#Qa1M`5v}RB~=AV@1WkQl)P4G{klaphvnqS4zyZ$ z@dprMa&k*e@f#h4(_qqyW*bxhj=_TiFG_B=GjT7D6^eWxVjf;%xWMfANl-O zre&4z{G^H1p?VONEciO$&keh9T~F2>}ns0ah!B%HlQTw%asZX zpnLol2XBwXXtIR__m^W`^2vM7V)RmviDx%#s_Acf;Ae)4Eg9jNSFD95 z3<;stU}QgMxWn=>ttYobY+=D0BxxeaMltQ^@H5QBw?E)ULoT_*&O^YZv}v)XKCp%V zJnrQ842+L4Xie?`aPGK-7jAn7EAK~$a;ih3i(5)RyOD0U&9ju%KOwRsJ9Ee&^aj`z zd{G)sqk(D3(_gxYjQ=D1D!U~B`3vlU*_g`ElBa{@s3pK@Yr54!jSQ~i)H&^2T;Fuc zh&3F(ur2y$ElGa5m7N%qz905!Z@&rkoL`zj32>%cTf{@6Lj$I(=tNM)G>qtf5x1{) zC;~qB=?d*TFH-b2gNwx*{smnP0JY@Lw(bd)M_Z)LCd%$eJe?`s`gekW2fw~_b`z+! zf5@S^ilm1U8?~C2|E5^iZaYhCX~I1x79Be5)dvql^-d}XNKTraPB9Rl?-!rp_Wo_e z92Yvd#|KI%Vc1T(lbf*FMVzu z*s0>@*TxND@+tGN)+aT{d6v8)?rsK{jPvZx(l`yi)&=bpJ&TazQj^tw5oqPZ5h>3p zGZ5pn^DRkCPv`#G;xc*zhod#(z@@Xf0~&l$c@xd9BUa#c!Tb-fn|NE_{~MsQwrE*t z<0z<8CidZ|KJL0`#$v&*WC*qS)_0EEi#y?!joItU_X;x|t1VxH*$PQEUA9+l@jONL zx?L3CLTNwChglSo)qd@KH1CRn^YQPz$&&8JMmM`R3=$KAtM-k}>J?m06K2z$=FjQTT{RB806nSW{A2mAS7^^PH5RbU>pr5HKNy|xI0YXjyOO+6 zO$}-L`fgM_wlj)NS~Z}Z{DFcll2ce~Fy_?Eo>fTgNLV{imcTg1wsj_L9+v(**#1HU zkB9e3st5fx)Zc_1+aw9N`K4@n8^W?JayM_)oC{>2__o|xk8clRv?v*Q9=bWYmRIq> zCG>WA`m0VSP)e&DZ&y$T@bZAgNI;GbwQ0q~(k?IL_?arbdk`7N-BZE3{d}8r$#ElN zMMuR6{YKLJu~>9*1u@=JbYItUKk@xO12@1ID8~~%5uQAZTdTFiDs0vTT@v3 z5EBO&XtGOo`7Di)wkQp;8BZt zbpnfow5$7)`VHMgXQ0)I-tRnC+u@%rdptyw!Pa+5W#51X`s=UGwS@` zDvdjTl}A(5WqwQh1Jp{2j%svR6rW?aMy$_6f`dnGnKDT&NoG{7`QebYi`emezaVH+ zD2-j4pYPM_2w-6e&%5&v zjmSrnS(oDiG!lYX?$_qwz)-oFq`$2!sQY~EN2OC%sOm3F@7vRmuOfd8BBFQ6zKMJN z4a8}lvN%aO%(C<~IZU5R)^9&eC$WmD3o-x1v*FZDIX-F^J>Zd zc|bX`vRpIS>Wm+XVdVSkTSswh=zjHS$1bEbHevrhs1@48UemT)PHx9f^PhfmfP6_y zK64xMfOQm42c~WTvxZWMOa=ddlm$ks;_o0USC<(#aZf|F-n+VaWl`AtOD;cxxE;^->;$fVM&Ad4>oaFv!uPCTPu85#l4gRVZ?I6I+NT6~ zaF@3U^H}gBFNI;W5CacG^}1E4Ig(UcBLXiUe4lD4^>duEkhhu4Se}4AQ2J@zj2h1Q zTK8#YTng}O=+yM&dRkKU*KIjOjJp2pvj)4 zSdjOYDBr!^C#z|0sfMMn*Lg{x;FZGI-fV-D-T`XCL9?m&H{mRd!E_nztg}T+Q^tIV z^LjqU{m)1_1nH?2eL`A$$gJ9 zguYkn#V7elofe0}Hq#G)b$*ksN(+9_a%JSUa0Y2UYsA(R8+-fL{Dwe{D)5ku(ODK^ zZAA-f*Rn=2S&kLwf9?JVKhtGv!6X~_F7P=f%}Q7(b&uXP)(5#_aFSl}>~a4^j>eMm zFzC-*%30!q@R}e*taU|(bPN#4VpNX#9<%h)wRP-L2Nx0dYc!@Nm&l2U}aMkvcpP$EjGeW5-=phlXzH z)bmPC2VRQQwEu}zSuM4WYkvgI!Guzb(kOck@y#ihk`g=7ieCni;oMiL`{;)`0r!}< zU2Civ?C@+%Wyi*e{foy+CsZIXc>8?#KNG^ynR33#MjJBTh@e~Z$AhSH?op+t1waI& z@oWqpM4=$1x;YCApndg%f?AFj=q=P0=fb|5)8^Eco1d+(w6=X@b$1z_V*r6GrH&K_ain$!1Fa`j5Ui+zbYEnh&IIeLRIO; zBBOmp>0omsu>?O}OKxAWzRYGv1yd(Wdmh&qbcmktLIOuv) zJnavY{bzUOki?Ks#FuQ`)T>aZ_s}a~!3_;awHp%ur{J*ibA?ipWZQ*U%jnt$INflY zQ^OJ+U#Ma6#zD7j^=hYHTzpqBJRVP`tRLr7WH6&#TQ>5g<(J@QIZ{!%Z{;l1Y$3=l zD!HnB7U*8#+!}pcFN=X*`}`odc6#nO?UX<#%fq#yMuUmCjFGn8NlC&EJ36Yf%9U^z zk)4Cu3LVy!E-&^`L=`a1?wthES9?}GQw?-S@Jf&EQ!xXrr*O0b10Z%UCN#TA(YUw6 zy>dhX;-mjg5|Yq04CgcAZ%w66fv(kQCJPl1U5hW1ZK{i;qRlVYzBH8~-z`;r;c5zn zQ5>G)4|%BsC5@IP`I3#+U}j}@**@JOcyg^@erh)Dn#GV^mWk>Du z;GT{|-L{!M1@1u0d<9(=(VHw>21*B{fACS;t`5CUZ(;BU?|#tg)%s@FfH>05tX@r2 z3~1W-@TSe;h@5#OV(H2kub3tTYqf~Q&(Du-`X}IkIQ-kGP@E&DC&K3l}K@|0T2Uq&grD676J}jv547QqfweE?ZfZpkmg2smFcd`C^OB_Au|Im_=n;_^BKN{+uZ|` z%_#yCp5@{gnui>3&>T*$Aec0teRObFWLaU|HRdTnT=At|_s<9$jHC|Uw9`OR7|3BH zXnO*l>kis5ZXixDkz4&jl8urnJffGMM4VpDht1vT{P%zIWZ~`CP+ub(F~Lm^oosl> zh}+Ak8B?c*r$LYGTK8<_e$t2R7ktdRF&0j^y2$rvj>T&!mi4dYV?nG-EF}%u;8$0k zRJlC5hJ8Kz+9~;3hMcO(zj!1YugBqo$PaV_Ku+P;;MRoDu7gk#?iVrCkyOp`=^(Mr zrJcp~qyaEJ(s?-c(|d&I0c-HynFcDIq?+lh9MV!HO}czi2*4nvcU%m0uz4$ra~K6& zwL@9Ws-o|xL_;nMomGOy0kIrPi__4ur6p0$a4FQ^C>~CDAb!dW7c#hm0Ofbn8E#d4 zisC!^dTqi_9E(jd83+sWz0=AjyJJFZb>3@W()u~doW?;jYeIV ztwPD^DajW(rdge1Il*UE6#CN#r5E2dL#hpjW0NrisF9&8OiL_j{S@W11#1ha;m^_` zJRiVetk?I>Lcj`+V4Yh2R%q^vW>Tu+rsuia>FK{$b_I zFX__Az-z@!DmX!*SqW#*^*!-msmjJ+FN4(lN{1vpEA)$qCGIqHQr3oj;f7=qk z@Tn03EZIzorb`4B`;JXSVNKhQZ)iIDQtU^YEa26jWe9njXt>Hv_}Cf!MyJ|@$ozxq ztKz9zkWa1lw_q7yndKfMa}~l0i*$@4x)@}Go+i0(&_!jR-p6B4UZNhH|550r5S^L{GGl7mh(fo)uu0Ct4O1MW=001q1+luQE8DEx#anhL z``zAd_=z#}_+4oDJY;p4qsc z_!wlK`RA9V34HhlC6ZPj&wuPsbFTKLB!1<#MY+;pzA{q{OkH$^6%k{s5&-05diY1)ju^uzyO+OW*DhK{=OJ zdq~&O8}j$P<-bTH=6@z0`~ab@hIjKrQ1UvYCt6D~4HVxjZbJm7?hxk87hSPEtf-)W z_O&GSxa*^74+D-YD=}V7FM4T8lDSG$%(liOjA%I-@Svd zg1WEUq~0A#vF|I@=8ucN5Zt|~B7@_VRh_@)FoCF|QKM#&A`1KAqrl5ug2;3e-a3&H zq?M7jv`D4|jz6C;>H8%Plymr3bzsN>@Iq{zDeQ4fLS%Kk<`+OI1;wdva?C){-+c@8 z2G~(2q6c;r3fFlR*w4ln&#UF0rRSPblRL{^8~y!&>u>(7kI-)qz@5l3Sgi{=$ew-< z%Awel6&g$Yqz1P2WOuy#0q&1VS^h?YnYsP(YC7VOCw+ky?`!M{5nJK~FD$>VQ z9_O12%|?;gh{FDd$`|}tg@pC*IKz{_Ab%Ups3FGOxBI%jeUK#TkvLA{aKO*8 z-tWrolpKfEyKnbAkdw+Kmf6)(LRp7HN`-w0|M2we!#-XDvD9&= zANXb!`!~35*e2G1mLp{S=KMtIhBS=J6`{ScBbTt;RoD0Dy|fyojGj)$Jn8En zKu%^U+V=Iq3%?@?q9m>Hk2Asep_68vO>axj5Q)FCn~NyGuWwHdBLig=GAoQz;BLM){Fg8{=S2{0l0SXX{#hvb`p&xE`=1dzvi zbG4N$ny~o$*~Px4fZFLT(>G=z+&@H4kt(2Ag|#gMk9TnEt)}dl4P*rLUD2wDsaWjV z4saFMip7#rRO6CS0+jTbUzPLkfC3xathW`tyqXi^a_h__aX#K=xS4Fd0_G9deYOUn z+HTo?JI*$8a&)*5FS!pfMr|M@Y=<}2&GCi$qaUSwkX z;7Rb%Nl!i)`X^91C^;GOn!9DUr)}u_^#!lRkS^ZmNz*rSez+VKoXbtlwM;_kOx#+p zh7spkPQAKU?_u*VSMdfsm5qQ)N8m~ z(bP&XUiecLzBmX+-B;yz_O#F^w%&OQLUj4!DjEelR0!Aef2aK6=Fnx=9}WWud|KKL zc{nH8H}*a!Zv*-eU#s%nWvf-Y=8HEGSNgEVl|J#zM%k#LslFWMR+1!&Vb%aPX2fSc zGUV+5UAN1-r;;wI09)lZt!p{-YVWds$hYm`{;45ZWoWrk5dIjCQ8yDEuConY!sYKr z*M;={jDhVw{B&zA?H|Ui3H~o$2zA4>Sj;TbCmcVy&8UKGAjM~1j3e)$1{GEs&1kuV z%FH@WYVnYirEnUv{}ZASw{s){|{I0$D%kt`C;2_;s{77r#(x1{2xfXT(7xs%?rfN`&}BL zse#09rM}Tm3m}fU)ri0H31n5rmKh^FAd#B2Dq9!}5((|v-?E#v zojB2_=^Mk)2Z{2ag_=u>R8mt4%j?xUy5QJOmf8R)L3j>+?yxgIKrGPm!X!F9>2{86M z|88m92NGFV=>_V?LHzOBma5XEB-9raejDupaZ^a@ z4qW3*Vl|(mINQ;f>mmRWMI$mQx0OKh@|W@rgL25vcF8@r-9gf5!QmQbD@gR6j+`-h z2jcw~uDi~(0p`3^_oEv_AQAgrw?9@BFcVacbJrw6hSc{sNG%+sKgnqAk_C|P&g;2f zUIF3(>&y<-+JNL_K3|>nB4FBM=XW|Pf!I{w{t`+R%Ig=$kM~J|)Tx|x_v0f$f?6sT zyNwJI!mq@#A}m3I`7G~VwIoPA;gxvE=M9oY8a<`G`$6KUsNNS(36M#BQTCwhIbb%W zpE$tP3Q~-#JQ8^$sE)RM+WUl225^tLNWvK0o&an|Im7+Iv#oX)aHPAtY6T#feMh!xiyiThj1>C z-z4ncFi7j@YH8^T1BRhmP?Y)&FcQ447hmXuJtv9PY@@r)ai2u5g_g4_{pt93uak4oeGa=^p$)*-CEkg#+E_?2 zb`wZlF7cae?*ysp?`2=do`bmC#Lb`oHG=fg1|g-9&mbKeWh&|I4N`kgr^YWkf<%?0 z;Bs{uu=m_zuFb-rPCmFY@<+c<5gX#7oPA|3Joky}j-8;~>3ZbWWlW)s^peVoLv>08B5Znj8HBNK=^9&=!;@ zNgzN=a zIJhg#<0s+?%7>xKvf0efB#x`P}M<7ePaHnH%{puAE1Idg)l?m=B6vv$Mgrs_qQF^X^E3^wFd}Fo8l zMC}sWC4@)j(ZxPhH$XB)c_a1W7KBr`Mz&^}gS3U}XQSLD6yG5_@t2-}jId26&)5n` z-&@rQGA=+o>HVh(FAPY^=PCpOS_(WCi+ANgNfOUP6o zO% z0J)FKPgbRGg1A`bClUPzAkOtcxOxKB?}t$;p9&*EthY%iX6qJ^Nq)qq^ywW)oY7FR z{D*Lx$;>HtaRXWBM@q>%ZUUw!DNyfXrT9=P@oIUF(|lxCilX;hB_&ojXBn`J2ygY#T^~M)-OY`apc+tK!WT%^>lb zXQzER;(v|%1j-(+>0%Q36t$Xnzz&!KYOWQda7(Z(3MqXQk;r@w><*q8gGUcg4=mJAw-N(WG13+^x{C55V1sKI!rLGXkK)bKdFZQM z(2UlnaNRotw35@+PRD|Pxg$7B)Y%!BI@@pUJQNFz9MOP|&~;!$2o;tw7Y62!^UI-U z{s7ae$U|z$9~c6^>Em+)!05mB??5sc7=yB%5rS&K7?3AX20sGr7e7w*?G!M0Z5iV3 zpMaiDn0qJp2aMM}cH)d`1ezmvsr1hdFsv(DuU;1nv>A@&An9tLY40t4E!qoola6%1 z{7ztAd9>GBLJjEeTZQd}_5g!^E4t-&8PE-r?>-Lt3(Rxtle|E(yALimHiMDQK?KGwpC5 z==XOD&!6=N+8Q%4Sp7QCb_wbWG#>-H;!IwKVj(b^)>6CMzk#8LLpS#U5*Tl4IwTaW zi0ldz`g=cFfyzh8`x@Em5n%9$N3TD!84L=M^|mj9NqzpfPp`Ir(OV4e(s!JQXkOTPoCP#c@XI1CFN#M zrGTbH``{)o3bZ3%!haja0VCGPbn=2LFpKw;Klhsi`km;LJ@3%`(KFT9ttSy~F>Ti! zk_Tr02i4(>{9FT9om%>ewrwtFO%$cH9U3?-}be ze;ETqqkZSa%v->iP?7MK2B4b;S#gjPK!5(P{l3P-Krhf5EjpNka6$G~?!hC#_@d<= zU$6j-)y_4go+=bqe9w36JOGTbpW;`fRlwNhfnwGfM_`09qz9_E0{weX-cBKVF#26_ z_p-|;FixJ1x6^S0!=fGm@B z6-`N72E$lK)ff!Q?=9aqGs@3`!GDbGGLcPSc+2vyB{d0{cgPa&qercv_082lKNP!Xvc4V||4Eb`b+Z75>(5h}+3UgR%^QK>qeH;BUcD^p zVE_zD)r-~|gmZ_Y+PG6R!9=#~*PPNI7%pk{GA@MztuAkSL*ykecy`yq>97_Ue;FP( z^?Cq|+VOj;yY10+(kN3{3iL~ZwP&7eMS1v0ed>WeFwW~LjrF^u@ws%OMs}fncb=#n z#sRI=MM_(?6&Uias}kpqfpL=G{_!WMK7U*py8SW`=n~tOS6pdeSTk@9+o}!*_hSN= zRVIK&iI9*j@CRnDQawjJs;@r(4fh^J>q^nNC12YCjNap=);$40@4r}K6_9}NxyEK) zMmaG1#*egrzXy62G234yi~%EhBCiu}0!_sE(o@-CV7z~!vG3erV2WgP@3%P*hSuxG zYhDQf!;y3+Zw}=XSBLw7l~=%2(V)y+?gtuo!TDvI%fR4EZ*=vF1m>2J@bx4H(B7X2 zwR?0049jFy9TRGRF~linJhvU$^KUCJ^BfqH`(>I?|7q|eay$K^5Exn&T*ta_gYhrB z;{RfMzz76AtzHEHZS04|Re}}Jj$KgooNEH6wBx1d2^8N`53be47lYx;jDk$!OJE-1 z{;XYx_V;4XxY#S&t$i0o{{C|(Xk#jhBkqh*0j&j~1q!sz9Ka&Jg`L0XjOuyI zgFgrGDag+Mj_N<+1jfDfZIxS}p?W>4eJI2YXm-yJ7H?ApL$#>}PwS09Qx{%wZ=MF` zqmE4hJC%@Mu9{UYECc=P$;T%?#(?2-`>o=02$!qtRiF7;0rR8z?JTt!Frqi6byqea zK9IihyY3zs%o%*Ne?J69?kBC!-&n;Vm=~6)Fd^mMu z+jTIy@@Hz{(RR?k(aP)F&P; z*tu;6m?r9SG#bL62f>f49Z_A(Z(ir&_zM`18sv(+QGGm~|HXB;4;VS?6gz>&dwJ#QB%vVMSMBvk3eo4&>zQ=x*9Ek;ioloY-AC|vtE3R*9;hkd|#crV+BT81Ui%44$RNX zTC2aP!2C*$S9*wWLios79S_3uJD*Od``$%;g-x+4T@z@=lIbRT(?Iu6(g?ZghU&~e z)r&uxLI0+S^pMm97;RjYv&rBD!+X~YUmihz=dcq{``ir7)Dtq-PW}TUqeDGY!brb( z=BV;rJ}?%!<+rSi?8$cx-*z(ujLM50X1?4*^@#HZZQmmDOK8d+b4#EliseQ!AAo+- zRtrNX#J8VcB;8v24U8VAi!|$0&Pv%ZPL1@z3o zS+&bdpi!@6^K_s%uWOB&znlYx8R{_)5{$vrGOGK6`MplKcrX$= zdr)f91?W*vf0O4Mz$EYNpWTtFz;HV>>TR_S3~EJ-Z^ZS2QPgnhzUc@sKL4a&r^^lK zFE;#pe(WmH*2%5J5dHz)l@uM>=gH?Yw(ZvLA3SunbBm`kE}3=9jc1v_+) zgHhnN{UYvfz(nIW&hs$pi`FH)k@$}bm{JCJ3M5kyPxdwa!jE|Jk9}R`M%^ebvYOZE z!vMQ)TVF0k^WFCLW=P%vI>)wMgA;dwR`l-o+GGLB_XnL5MTno?KJcwSBL(GGqN{SB z0T?V@-*O~35@-o`I*V-K}XXMoaH;o@feMMX0 z&c=B#d~t58dD$5-dQponJbnX=wM>lcI~>5+zww`JSw7H&@Ba-yNdaavB% zGnQ(55uTs)pJ_(8(AMylBF;L7uWEXqyLH(`z z7NzZ#VDvio?av1hK*K-Y^+yH=2F$9mE!P}?X|K3pka-adL?hpdQEvc!__K^LzaEMg zXGv8F)Sve3Z~pZ94G8UQB2=uc0PlqtuL9Rr!AL4#<=NUIP^>opxWIw>V3Js+BnKC$ zAJe^Vw6+Ekw;c~ubMOM$M{n^V2OlV$lKv`-o|FAo%2Kj=4i3<{ zNJx!Y8wZ6uCMPR7I82Ac{>Apq zwaxk3$PU9-Vvm^|LO?I!naDyqO|kpDAqNi__Yf5n*5*MkY!C4xnxCl3BiF_u25K>i z9JIB6rvJzvXRdRI16kgK?ycH5$iI!nm*$>5`$jKK#fpftQtV&*2gLR94_wf9dR&3`UvdZn$-Jsi0li0dEBLTXvyB!wSt)coYm@aS}%_ki6biFbfmvUwS=YssBo1$9Su{HzL zU=an)wO@cSACx+U_EE!q>Tm{!1gNUrx3@XIT?~y+|5PjEfYxW6a=*D0txGdf z(uw6)vgqU%)wMsSD{I;(I5^e;DR@uZwzUz^70W6BG>`hbfd^;MxY8YEANL@;sUZq$ zkCCp)M?!I;06NlJ43^P+YESfpKeE;zmvngt+6T@0g=KH_p7@&4Y9`t*T#J>o1_!F! zW3jbG9K5FUx%HcdkX_96BhedByvXjdw0ec!S9zu%dJ6g5IDP-In;iV6Q|lCl>d|vX z^A;NcbYDJA=Sngfk8V7jI*a^HiEzjWMEB{KwUbBDd$gt2ALrJVL5DFgn2hYB5R7FG zqr9OlE^t3V?`aJhEdhr-;3E&a^q}`iN7_YHI52?UY$q~^)<;i1Fg}6xx|ynxn0 zxEF9!2gSRJe(>r6*1Gl17#E>+>eas&C8773Eo!^M&^q)R_yskQUztW~IyuPioFcET z+9EskEI5ytqInq;7jFJWS0l+O#RF?|V894oiCy~*#)`k67jOuG%xXv~iWAd;==(mb zeNa#qnRH=&|Kle<4($VGd~iE!zV%~rS_jcM5@Q>DUZZ)b`meqnLH3zc3hHH{byK+; zYb4Nka$~U^F33)qa+|U%=sW3Q-93XSA4Z9#y>Vzo@g&Q}|aQ#DeYURA}7v_)!v6K2g4oZnq5NvmCw1qu!Pn{d-Zus9C}X$uRM4P&8vDw8qjhKR)r2+Fe_ngu?SR%RabErX8oIAq6ew+t;!^V3$xWM) zJz5&kXJyuY0<~{pBNpX_;pevT*Jyn09rGg<2&arsC#X-NammqwOf2`+|jz#%0bZ+N-;j>&zVY}8vU|h#qU_PN zRnwWSAET^zVuW*fFQI)o>yU~>2S@iY{VMxAEo zl+k(w!%E#*m?Z=VbhK?y4AE0$%pSjHMLF?BrcpjaJ z@>VZdd}W2TUv4-;E81T}KmKrEv_1}*=l(Vb$26k4w>?7Z(em&6eG}O$h2Q`E6&jcN zQ={bgDVvzG_HY?v>y>Yr{hkgvFe|A);T;2 zKQ)98N9UmR8W>v_{zl^|I`v~&{6TLv`_>f{rvw`(g?dysl)uKm^hWQC_=mq}MdLHC zDN6*Q@zy&>e?a`qbnWaf?>NL~5e~^WHU;eyNMt;%zCq0>qzS9bB z^`GI81fw*gBqs77?ZMl?T9i*FisaDOC||@QV>@rMa5vFBFAeR7K=6*lXXrbFNE>1n z%ictj91RXdpvdQ5&O>-%AkXw<)e|M2(bvhWI#g};Y6^{uA3dp+itILWd_z5q_DA3A z($#XbZql|tK~^X}jIuX1mZ7-inynrcL+ijUCzPH;{6Qt$w@Vo1n^sxKx*2rk+P?Kr zJsQvG=W^5k;C03c%N$l-FF5m;qx=?#*n(AB6jrLK!r*DWK* zVv!$ZRjM0Ub%L_sRphxgZ@OTfUZ#TfQFVL64i-Pfw+m3_(0UAaU;5_3;sIYp4+|k( z@9hZnDYDOCg3PRE#fi_KhyRej3=JMrDd;(cH-CT&J*WLwajTH!rirEF zc*IZ0iub>WqC6pP9NuGuc%OFaUUOF5=}Z(WD57!k`_)eELH1xTw4N73ek3So_t z!`N09uOk*y1w>hRyt=%Pn^pHeSGB01`3!G(5#rIh#b10jibVUxWR5RAM)qn3ZG5?c z`T@Pd@+@mKj!=H(zh~%sy^1~-UjxIVAduA;Ta@ID&H_``?Z&uE_vHN9fQkUeU?UgR+J9c4VcC>V_^+)-=K zfxgq^o&6byc%wm*pLH3^3yuldFxLKPKNd}|h|a2r-jI`a@Wmoo z{Gyq^ayS9ugTWb#n;~dC5x-tPK=-K&_hr_jaWs1PuaS^Fdib0F{$uU`Ly@-6=sEp4 z{p}gF9>UkNtrZB@_4)Mm16cKc!xxVj6qnReJH!kA0b%>g>QYqy4F}W@vHDc$0UsSQ zi~rU`>^N0iW6?0!LEdD1RZ-`|u zSb4x~xpNPVV`S7_$>PsaE{WVEI^q!;3;zCibY;j*+)iiJ zFJ1nq<1E~`z9hrq<)nk#gfb9rVSAojQ$T&$=o!ZVR{z6X8|-dk^`FthE?lgBk>||C zR1}9)=RC=AG_OLx-i06*&v0t_s)*h*5nFWFfbwkpsl3TNv@ZPvCm%MW^=P}AuFkOV zY_WIlHG0oLd-^9Yve)30^ad_g-e<1}v(6ugxLyrz#IM9(gpv?XH~q}C73M(kL|FP% z63W{5KdnKVS-3ooXFNjtpqje(&H=5?+d z;xG5NbN8@t1h=o$7x_)&Z@2v+#CJ5x`|k`Q{wcwEe&8UACk?x>y|SpUXJ~mcUm|-L z_LYK*ES^f(hM7ap#pWCYme5s?cV}29x@s*N4RNn6n9fFQI1gxC!fDZt1mq|2M{hi1 z(Y(y(WA3&rKCt`O>=~rX=mcO`=Udw8>8)jGeiPi7kP+G+lf~ZPJQjbH6|#(B)zM4Y zZUl6dGI_F@g;Of4hF|#6Rr{%Od@{1nXyC~KV^)8&)%ivw%m4TJ!fny|)EWH|J!t=R zs$aR!v*LDH{%JHT&%VeE@vwOQx)%{nXk6{kUfi@YX z$BIvE@_E6FD6Uk!zdoo(_el?EnT=>1ecdg)G+252JLN<+!UqH0ibpIyWJLLI7B~b= zC%f7S&scpbPA`*nUSk|H{6iJtkN&M2MlP&A_NA^9AUllq2j$NF&v%PTGOY8ee7~DZ z=t_$fZ`gwJ%rO6;JB#;9V{!*rb(5Cp(r=0UPAZyF9$cFRjH{Lv%Wgf36MHtJdH*Ke z=N?4%X<&c)N22%D%#IfOA-sBYEc`VGnvZV)ej{8keOI1*dyVB6tpbHJXkRh+e{Evz zj{(&vas=5UEUNy6b-qab7I7^a&8O3Ntz$Fd_ZvU}^uM7xp% zNg=d9)G{B*2v%Rnk(ss^-6w4g>J??xU#H1!qKIdz4V3(`LH1J0oLd7}=gnhy)fTjl zby;Z&(oq4%*H#0@Rb{t_8HM#dbHre7MkYFTm2uvji<@lxf0*HWzWLp783m$m3VhySC@ zgJ{1fI0uj2=sVe!mwiL5e81G1yN2e)TQ1!6M*ioHo9J&v`nr8To&>SZCGR$Mpf~_5 z568njXny^Fm@qo4-g5?i+k@;Nn!nS%$>RB9&w^9dew)@_PAXxI$7Av3iYW4{+^o1v zI*V5_xac0}J>4e8?>%TdsTxc$3+Gh_tRx;GKbfRj&MBbv%7u-ovCfMq7bjviqIfsv z-O(d~bg8P9ZaS;}yJ}zEg#1Ppv74wz>(rmwG`A0}i*)yaq9YoQcEo2R>-<)2;~t6( ztIxRdqG2b}Ri+|DzM=0;gj{QOqWQE<4T9!a=Qp+kXXns98pl30mqmV&D!ed0#yYQm zdE(9r!cnn-yH}H0yzN4J8>>&0%KT-^s*gG^q+`yZ`mg;sZ(=it9OBQMTmdYev->|? z7XKv1i4^Lv&dse3Sk<6-G%z37CyU|&7qh!c5A_>VoKry(;#r@6?mQdEI>)$m-D-w~ zcWHN*SapvgP2D%k>dT}!DZA19s+r_J7tuN>A|6{u(AA(X1mDN%Q#V8=EVKGfjx9H0 zSv;eb8qYd^A`e6h^rH0}?#$0)^{bj050mm(=RRtj*FU2ENdKST84PQD=A$kAC{IOJ zu7Q`LGS~h}o;f<|?dnzUVw|N90fAFn{qz zRzBiKg?^y*s{9b+r?YtQc*qr2zlUYY-(BU90sP6GZnhkXrr*RKHv+3pPxyU{VAZ`F zRy$bnNcOgIUPAsQUUq4@hj2)Sm&N51sBd#v-*2e4*@+iu5sV6<~HOv?d^Mh;Ltc& z<*s}Z!a)PqmfI!hKHW}%mxa^D^&+{f{*O_lMrmW6m#ZyK;n6yzMU_taGwe{YuvPl78bafs?F!s4t<> zSp5%=nyx+ze~-&?X(AnyYqoLQ2N0K8B-qlMd^55+fBp zi&?aOqs9}d)u_*CJQXp?f&5IF$#^`#I?ouZ7C`k6jKA#d*G6`5eBn`N{XSq=pxwp# zT}2W5b3-1ohwq?v6{|jSXRWPPART9$m}`&r(bzD1-72d8%*sd-i*M*X96Z75)5y-w z2CXPBIksB}vVIrSpcS5F`QN}5H~#^}1uhfJM$vP!nOrFAT$ui&@xvsFD^X7F%dGQb z{q-Bajj{H*r+up~8dpg$s-=Q;9=7{VSS71&op0B<%j$DBQl%?cea4K7;t=8$g+>!= ztn*s^Qj<+TQM{6sD!W+cj|MiHo{nhUhT#JnwxM}3y)DJstlu98g-2Mtff#A^;Wgr8 z8fQN3qa#0)FP;nZW1Sxfy*~Ssg}@yu5eo_-R&uG*5aWgz$iNBJ?63vd84ZKr<`u zbb6U4GYHQ}epN4p5WY!l70$L|;i+M$8V`$ayAKa*BEOKo@pqW8`t9~~V;YMOTsjbR zkM(=1(V6Ui*u~=q>`cpAc`4_Bn6;-8HBCFnN zkwc>#P@Krl`#fUhujqQ;r8u-sorgYFtlwEzOEJrXs9z>`pUl6Xl-y?d zZdQJ6BcwO6_}wvyTpUYREP6D+%73MDF%0T27kAq&4k15icUMO$BAitTJ|o7$b;f0} z+j!RRS9|4Xto;;|*we+@52Ndf&WkKO&@=eL>Z2L7$*Bz}erQ7pKCC_-vyp3z#hLL_b{8@5c_%2re?(~2EK4x42R(5C}?#Ew~ zNhsfqx}r~%uznv39Xct3#!+znYQ*Yii4_GE7m@#sei+KL___Wb&8gnCUmP4fT&PXu z{9pfXbl$_R0iK~qKoG;h^}lZNg?M`V1RnlhnsA_7sAtF#*D!P|^#Az=TjB)5cYilzu=<@IXcv+Ze=X%&9Fn}29>FM!* zzt0uI!J$)#{qailr$;Ed~%qXc-d|UQ7or_>NWc_&(9VN|t5gIq4>=gIphQYx)@@ek14ti(6 z?B)KeDyDlnpIYh#yrItr$N2cuw*J!_Ez^HQtSH`y_cS`|sAv&%Hzbs1G1*f3Nt^eN zzz#A)_hz(tKq=kl`y0yV<2@IYza}N!KjQkhKX9%_p36<*&;5@Ad!@Q_r1Rc|>mTrV zC$zTpTBg>6M%qSE(U*bJ!oL~yp#v!i62yvpKJB&i>UZA`uiIrO_umouQs+(1+1_S< z9NS#24(7hkGRSSJZs||wQY)FtIW~RaMZ}+jzSjaW#`gxfr$14}?AX^Gp)@x8PDlEg zzy0|Q{*IgOo7>xY3NL-DRez~H-f;YwL1Ju}*-`ubZ4)Ppw>0KIJHPnWa(@YboAqI} zxtrno;yft@u(RjkhVulMacZlZLJYA1ORW9Rddss@qS@8h|FW{Ku14&0vm5_y`Y!s_ z<5P}~_`}!C95%lYIK1uf(@nRs-@aYW7XGZI^CnQSI(B~vLNnTuSItGEJvM|TNaw$ zD#sjnl5AXUb})9=l4NIY$DxVcVg=?EwVJt}Ss+twyl}j9avOalXI^2Q9wBPyvnQ`l zn3pW#&3q|SX{+vq)rnM7qFEyKv-Yd^9=CI^wDw6dh^NV?-`VS8h!tZ(M`Cpi+ z4i3RwY`nXdFtxKB=#GrdT3an)m7G$SdX(R|;#XpZpJVdHaE(zai1?-rc3eX_ab z)gR&AHzo$YzPM{Qeeh;n^XV7|?`ZL{nnmm%-B&WUrI(xV&Pq~Y6K9TJxcoNSrT@EK z>OZxgF1OZ&42Rnc}?&ufF1 z>K)PzLK|eG8dbIyFV#*IQ*ASE1q^!)&mM}h39w!gtzRwMWB%)OkPXKX)qi1urS<$V z5$9_ZeoeP^C<{y+symf&&WYNdp0vTmr}sZ=dE*h4$jP_y8=vhPD^^>Z(9+Q|Zr0}u zxnqP2YASf6u_iuDceCxymg4-~Hha6cHq(Wk;h*6Ue@F2}y`xF9@zyc_^n+Kgsu^1x z;dd)KDZHE&o+DN+BGF^gA9u#>p5%+HDxI$ZP0VY(LT$koO#dS75dDUZKNrq${pYgl zZum0^-Xpx-H3jS0<~O&sAN5@{{Ps zf#*~uA94*h6#o1^q^-60=&aRZTbNT+L2l}6ozsiW!tm3U5`ot`e z*JWA`HJFmF^ym4`6Izo$9MUE)T;SyDRWUVq-tx0gZOAersDTtzrqi<8U4}2}tL1y@ z_e=Nu5i+m(T9HHSo$B}u?fs1h)6>0lt{zov61ua7t${JN$6zLg<|&G&GJTSrpI&Um`~*+sn3_TR^!cnj5;`2#s}7pk9s za;ct96%O2VbSSpG!HsjT#m#%x#?~npyN^yTGpOs5e+9`>LM?oJFg-%e>d{4`$>gN6 z=~}&TH-}%>)6KWT)b$@rT(#=Eo7=udn79!dH;l#nFm9~@Ut(oB#2X?!n{tlJsOSaeqAHR<4364LJO&`n3sF~=5|BQd& z+Cq-0Ndzg~#cDh|*F&Ftrz@|nZ_$S3NHI1LNyxx2KZ?KO)m(3S=7zD?ch~gOWtnO_ zyGHm^8hcVn;}wdYTQV+`&=%r0@{lE$yt{SQ=J=Su^TgZ!9_sZ!fxmbEWz=nnnVsH4 zYLeS#gyrsJspWn>Bc=TQX!RqW+jqJR2bWE)eQbYyKd!k=#Q(z2yLH5pic{PQ2ZCY* ze~A1|er2RH;i$Cvw8PuzjpxRUW*tB5OQ$!EsERgy5)}?%j1^mD9q$zny(fH`@;85e z$8g8Vl8QgA#{Sx`c_u`>obrb*23}WQkQQLH%jgcPguKrYNV}YGv|GO@b=Ealo*JWl z&h%q-@9bdr;ila8s}aBb0}dDgyZlT%c78y>hAw#bwdIT1;+x2BmVz%-4=;p-as;W%i#WSIk zk^?92YG&ocl9)&61b8T@R~4D)-| z`{~-w6Hy05--|d3RSX2u$F_&{`#q}23WpQL`Ws9Ne*K*OPdGGZ59FUe!D^E(ixZ-*DXfSSF!j%&{7Y*rM z>wX*D|DgHdMbkr9(LDQS$*DgaA@()nJvT5G-J4v$kvA+ nzv&yg;@--;ygo~4-KFw79bC1p%5HirzW09t#X6sVAWH!NJMSss literal 0 HcmV?d00001 diff --git a/.Rhistory b/.Rhistory new file mode 100644 index 0000000..61b7daf --- /dev/null +++ b/.Rhistory @@ -0,0 +1,351 @@ +ls +pwd +ls +x <- 3 +print x <- 3 +setwd +install.packages("swirl") +library(swirl) +1 +swirl (1) +library("swirl") +swirl() +1 +swirl() +numeric(0.5, 55, -10, 6) +numeric("0.5, 55, -10, 6") +numeric("0.5", "55", "-10", "6") +num_vect <- numeric("0.5", "55", "-10", "6") +num_vect <- numeric("0.5","55","-10","6") +num_vect <- numeric(0.5, 55, -10, 6) +num_vect <- c(0.5, 55, -10, 6) +swirl() +swirl() +0 +tf <- num_vect < 1 +print(tf) +tf +num_vect >= 6 +my_char <- c("my", "name", "is") +my_char <- c("My", "name", "is") +my_char +paste(my_char, collapse = " ") +my_name <- c(my_char, "Lan") +my_name +paste(my_name, collapse = " ") +paste("Hello", "world!", sep = " ") +paste(num_vect(1:3), c("X", "Y", "Z"), sep = "") +paste(num_vect(1L:3L), c("X", "Y", "Z"), sep = "") +paste(num_vect(1L:3L), c("X", "Y", "Z"), sep = "") +paste(num_vect(1L, 2L, 3L), c("X", "Y", "Z"), sep = "") +paste(c(1L, 2L, 3L), c("X", "Y", "Z"), sep = "") +paste(1:3, c("X", "Y", "Z"), sep = "") +paste(LETTERS, 1:4, sep = "-") +5+7 +x <- 5+7 +x +y <- x-3 +y +0 +exit +exit() +c(1.1, 9, 3.14) +z <- c(1.1, 9, 3.14) +?c() +?c +z +c(z, 555) +c(z, 555, z) +z*2+100 +my_sqrt <- sqrt +my_sqrt <- sqrt(z-1) +my_sqrt +my_div <- my_sqrt/z +my_div <- z / my_sqrt +my_div +c(1,2,3,4 + 0, 10) +(1, 2, 3, 4) + c(0, 10) +c(1, 2, 3, 4) + c(0, 10) +c(1, 2, 3, 4) + c(0, 10, 100) +c(1, 2, 3, 4) + c(0, 10, 100) +c(1, 2, 3, 4) + c(0, 10, 1000) +z*2+1000 +my_char +my_div +sqirl +sqwirl +swirl() +swirl +swirl::main() +swirl::main() +swirl +swirl::swirl() +getwd() +ls() +x <- 9 +ls +ls() +dir() +?list.files +args(list.files) +getwd +old.dir <- getwd() +testdir <- dir.create() +testdir <- dir.create(getwd) +testdir <- dir.create(getwd()) +dir.create("testdir") +setwd +setwd("testdir") +file.create("mytest.R") +dir() +file.exists() +file.exists("mytest.R") +file.info("mytest.R") +args(file.rename) +file.rename("mytest.R", "mytest2.R") +file.copy("mytest2.R", "mytest3.R") +file.path("mytest3.R") +file.path("mytest3.R", 'folder1', 'folder2') +file.path("folder1", "folder2") +?dir.create +dir.create(file.path ("testdir2", "testdir3"), recursive = recursive = TRUE) +dir.create(file.path ("testdir2", "testdir3"), recursive = TRUE) +setwd +setwd() +setwd(old.dir) +1:20 +pi:10 +15:1 +?`:` +seq(1:20) +seq(1,20) +seq(1,20, by=0.5) +seq(1,10, by=0.5) +seq(0,10, by=0.5) +my_seq <- seq(5, 10, length=30) +length(my_seq) +seq(1:length(my_seq)) +1:length(my_seq) +seq(along.with = my_seq) +seq_along(my_seq) +rep(0, times = 40) +rep(c(0, 1, 2), times = 10) +rep(c(0, 1, 2), each = 10) +#a scatter plot of those two variables (time and the amount of cumulative drug release) +setwd("Desktop/Biocomp_tutorial12/") +a = read.table("drug_release.txt", header = T, stringsAsFactors = F, sep = " ") +a = read.table("drug_release.txt", header = T, stringsAsFactors = F, sep = "") +a = read.table("drug_release.txt", header = T, sep = "") +a <- read.table("drug_release.txt", header = T) +a <- read.table("drug_release.txt", header = T, sep = "") +a <- read.table("drug_release.txt", header = T, sep = "", fill = T) +a +a +a +library(ggplot2) +setwd("Desktop/Biocomp_tutorial12/") +a <- read.table("drug_release.txt", header = T, sep = "", fill = T) +a +library(ggplot2) +setwd("Desktop/Biocomp_tutorial12/") +a <- read.table("drug_release.txt", header = T, sep = "", fill = T) +a +setwd("Desktop/Biocomp_tutorial12/") +a <- read.table("drug_release.txt", header = T, sep = "") +a +setwd("Desktop/Biocomp_tutorial12/") +a <- read.table("drug_release.txt", header = T, sep = "", stringsAsFactors = F) +a +ggplot(a,aes(x = hours, y = cumulative_drug_release)) + +geom_point() + +stat_smooth(method="lm") + +xlab("Time (hour)") + +ylab("cumulative drug release (ng)") +library(ggplot2) +setwd("Desktop/Biocomp_tutorial12/") +a <- read.table("drug_release.txt", header = T, sep = "", stringsAsFactors = F) +ggplot(a,aes(x = hours, y = cumulative_drug_release)) + +geom_point() + +stat_smooth(method="lm") + +xlab("Time (hour)") + +ylab("cumulative drug release (ng)") +library(ggplot2) +setwd("Desktop/Biocomp_tutorial12/") +a <- read.table("drug_release.txt", header = T, sep = "", stringsAsFactors = F) +ggplot(a,aes(x = hours, y = cumulative_drug_release)) + +geom_point() + +stat_smooth(method="lm") + +xlab("Time (hour)") + +ylab("cumulative drug release (ng)") + +theme_bw() +setwd("Desktop/Biocomp_tutorial12/") +a <- read.table("drug_release.txt", header = T, sep = "", stringsAsFactors = F) +ggplot(a,aes(x = hours, y = cumulative_drug_release)) + +geom_point() + +stat_smooth(method="lm") + +xlab("Time (hour)") + +ylab("cumulative drug release (ng)") + +theme_bw()+ +theme(axis.text.x = element_text(angle=0, vjust=0.6)) +setwd("Desktop/Biocomp_tutorial12/") +a <- read.table("drug_release.txt", header = T, sep = "", stringsAsFactors = F) +ggplot(a,aes(x = hours, y = cumulative_drug_release)) + +geom_point() + +stat_smooth(method="lm") + +xlab("Time (hour)") + +ylab("cumulative drug release (ng)") + +theme_classic() + +theme(legend.title=element_blank()) +a <- read.table("drug_release.txt", header = T, sep = "", stringsAsFactors = F) +ggplot(a,aes(x = hours, y = cumulative_drug_release)) + +geom_point() + +stat_smooth(method="lm") + +xlab("Time (hour)") + +ylab("cumulative drug release (ng)") + +theme_classic() +a <- read.table("drug_release.txt", header = T, sep = "", stringsAsFactors = F) +ggplot(a,aes(x = hours, y = cumulative_drug_release)) + +geom_point() + +stat_smooth(method="lm") + +xlab("Time (hour)") + +ylab("cumulative drug release (ng)") + +theme_classic() + +theme(legend.title=element_blank()) +setwd("Desktop/Biocomp_tutorial12/") +library(ggplot2) +rm(list = ls()) +#1. a scatter plot of those two variables (time and the amount of cumulative drug release) +a <- read.table("drug_release.txt", header = T, sep = "", stringsAsFactors = F) +ggplot(a,aes(x = hours, y = cumulative_drug_release)) + +geom_point() + +stat_smooth(method="lm") + +xlab("Time (hour)") + +ylab("cumulative drug release (ng)") + +theme_classic() + +theme(legend.title=element_blank()) +# 2-1 a barplot of the means of the four populations +b <- read.table("drug_release.txt", header = T, sep = "", stringsAsFactors = F) +b +# 2-1 a barplot of the means of the four populations +b <- read.table("data.txt", header = T, sep = "", stringsAsFactors = F) +b +head(b) +b <- read.table("data.txt", header = T, sep = ",", stringsAsFactors = F) +head(b) +b <- read.table("data.txt", header = T, sep = ",", stringsAsFactors = F) +head(b) +ggplot(b, aes(x = region, y = observations)) + +stat_summary(fun = mean, geom = "bar") + +xlab("region") + +ylab("population") + +theme_classic() +?geom_scatterplot +a <- read.table("drug_release.txt", header = T, sep = "", stringsAsFactors = F) +head(a) +ggplot(a,aes(x = hours, y = cumulative_drug_release)) + +geom_point() + +stat_smooth(method="lm") + +xlab("Time (hour)") + +ylab("cumulative drug release (ng)") + +theme_classic() + +theme(legend.title=element_blank()) +a <- read.table("drug_release.txt", header = T, sep = "", stringsAsFactors = F) +head(a) +ggplot(a,aes(x = hours, y = cumulative_drug_release)) + +geom_point() + +stat_smooth(method="lm",se=FALSE) + +xlab("Time (hour)") + +ylab("cumulative drug release (ng)") + +theme_classic() + +theme(legend.title=element_blank()) +a <- read.table("drug_release.txt", header = T, sep = "", stringsAsFactors = F) +head(a) +ggplot(a,aes(x = hours, y = cumulative_drug_release)) + +geom_point() + +stat_smooth(se=FALSE) + +xlab("Time (hour)") + +ylab("cumulative drug release (ng)") + +theme_classic() + +theme(legend.title=element_blank()) +a <- read.table("drug_release.txt", header = T, sep = "", stringsAsFactors = F) +head(a) +ggplot(a,aes(x = hours, y = cumulative_drug_release)) + +geom_point(size = 2) + +stat_smooth(se=FALSE) + +xlab("Time (hour)") + +ylab("cumulative drug release (ng)") + +theme_classic() + +theme(legend.title=element_blank()) +a <- read.table("drug_release.txt", header = T, sep = "", stringsAsFactors = F) +head(a) +ggplot(a,aes(x = hours, y = cumulative_drug_release)) + +geom_point(size = 3) + +stat_smooth(se=FALSE) + +xlab("Time (hour)") + +ylab("cumulative drug release (ng)") + +theme_classic() + +theme(legend.title=element_blank()) +a <- read.table("drug_release.txt", header = T, sep = "", stringsAsFactors = F) +head(a) +ggplot(a,aes(x = hours, y = cumulative_drug_release)) + +geom_point(size = 3) + +stat_smooth(se=FALSE) + +xlab("Time (hour)", font = 4) + +ylab("cumulative drug release (ng)") + +theme_classic() + +theme(legend.title=element_blank()) +a <- read.table("drug_release.txt", header = T, sep = "", stringsAsFactors = F) +head(a) +ggplot(a,aes(x = hours, y = cumulative_drug_release)) + +geom_point(size = 3) + +stat_smooth(se=FALSE) + +xlab("Time (hour)", font = 10) + +ylab("cumulative drug release (ng)") + +theme_classic() + +theme(legend.title=element_blank()) +a <- read.table("drug_release.txt", header = T, sep = "", stringsAsFactors = F) +head(a) +ggplot(a,aes(x = hours, y = cumulative_drug_release), font = 12) + +geom_point(size = 3) + +stat_smooth(se=FALSE) + +xlab("Time (hour)") + +ylab("cumulative drug release (ng)") + +theme_classic() + +theme(legend.title=element_blank()) +a <- read.table("drug_release.txt", header = T, sep = "", stringsAsFactors = F) +head(a) +ggplot(a,aes(x = hours, y = cumulative_drug_release, font = 12)) + +geom_point(size = 3) + +stat_smooth(se=FALSE) + +xlab("Time (hour)") + +ylab("cumulative drug release (ng)") + +theme_classic() + +theme(legend.title=element_blank()) +ggplot(b, aes(x = region, y = observations)) + +geom_scatterplot() +ggplot(b, aes(x = region, y = observations)) + +geom_jitter() +#2-2 a scatter plot of all of the observations +ggplot(b, aes(x = region, y = observations)) + +geom_jitter(size = 1) +ggplot(b, aes(x = region, y = observations)) + +geom_jitter(size = 1)+ +theme_classic() +b <- read.table("data.txt", header = T, sep = ",", stringsAsFactors = F) +head(b) +ggplot(b, aes(x = region, y = observations)) + +stat_summary(fun = mean, geom = "bar") + +xlab("region") + +ylab("population") + +theme_classic() +b <- read.table("data.txt", header = T, sep = ",", stringsAsFactors = F) +head(b) +ggplot(b, aes(x = region, y = observations)) + +stat_summary(fun = mean, geom = "bar") + +xlab("region") + +ylab("population") + +theme_classic() +#2-2 a scatter plot of all of the observations +ggplot(b, aes(x = region, y = observations)) + +geom_jitter(size = 1)+ +theme_classic() diff --git a/answer codes_Lan.R b/answer codes_Lan.R index e180e1d..d1377f1 100644 --- a/answer codes_Lan.R +++ b/answer codes_Lan.R @@ -6,7 +6,7 @@ a <- read.table("drug_release.txt", header = T, sep = "", stringsAsFactors = F) head(a) ggplot(a,aes(x = hours, y = cumulative_drug_release)) + geom_point(size = 3) + - stat_smooth(se=FALSE) + + stat_smooth(method = "loess", se=FALSE) + xlab("Time (hour)") + ylab("cumulative drug release (ng)") + theme_classic() +