From ef93880350ba2621abfe6681d449fced16cf689a Mon Sep 17 00:00:00 2001 From: pslcool <50282998+pslcool@users.noreply.github.com> Date: Sun, 12 Mar 2023 22:58:14 -0400 Subject: [PATCH 1/3] Add files via upload --- plot1.R | 15 +++++++++++++++ plot1.png | Bin 0 -> 3866 bytes plot2.R | 21 +++++++++++++++++++++ plot2.png | Bin 0 -> 4412 bytes plot3.R | 27 +++++++++++++++++++++++++++ plot3.png | Bin 0 -> 3893 bytes plot4.R | 35 +++++++++++++++++++++++++++++++++++ plot4.png | Bin 0 -> 7610 bytes 8 files changed, 98 insertions(+) create mode 100644 plot1.R create mode 100644 plot1.png create mode 100644 plot2.R create mode 100644 plot2.png create mode 100644 plot3.R create mode 100644 plot3.png create mode 100644 plot4.R create mode 100644 plot4.png diff --git a/plot1.R b/plot1.R new file mode 100644 index 00000000000..d561dc0688a --- /dev/null +++ b/plot1.R @@ -0,0 +1,15 @@ +#Course Project 1 +#loading and getting the data +data_raw <- read.table("household_power_consumption.txt", skip = 1,sep = ";", + na.strings = "?") +names(data_raw) <- c("Date","Time", "Global_active_power", + "Global_reactive_power", "Voltage", "Global_intensity", + "Sub_metering_1", "Sub_metering_2", "Sub_metering_3") +data_raw <- subset(data_raw, data_raw$Date=="1/2/2007" | data_raw$Date =="2/2/2007") + +#plotting(1) +hist(data_raw[, 3], col = "red", main = "Global Active Power", xlab = "Global Active Power (kilowatts)") + +#copy the graphic to the file +dev.copy(png, file = "plot1.png") +dev.off() \ No newline at end of file diff --git a/plot1.png b/plot1.png new file mode 100644 index 0000000000000000000000000000000000000000..07efb4ae77ab5c4ae7ed7c773a1aae3474abaa02 GIT binary patch literal 3866 zcmcInd03L^7N^og%_TE!ajV?Q6e~anx6#VR6kM__H#8MfGtEpfMH4DA3z^I`*GWyy zwa{=)qckhDa$j&E9a9k(TnJ0yet4(n&L1=PKJ(06o|p4+&iTC`=e*zfo%j9X>`=#~ zBz8!EKp-irpN=?yKMqn@q3?QipG8sVzBm=WH4M1arI6#tx5a2^5W5`qt0H6XxKt_PbYxBt% z5P1#so1aQ1QyF9imCOL7fQkW8F$gLd7#J#*0eGl@Q9z{%7*v6PK}Imh7zPy>YjD6X z00x605C}jZU^0LZ0)+-495B1--^Bvyf3-Sd?ubpF&Pa&qf4!aPW1!N4>^^VEsbLd2 z-IoJy4Cwk52cDs;ZQu(U&V_4<>u>awR?h$SIl|PBeRFgxPIVBnD}y5XQ(_YzB* zPw==XyI&gLWd6hutRO}kQ9G^lPa5A!@)wb2^vnc>_`a^?)sU&EDy|Hj)wQfYRnl;< zp*DgybmxoDjq@2wl91w?dvm&vDR-0A=JCST(sEe}zz1TmE= zkH7eME^C#|4)@aGK~~vCPdO18J|8Rt@bT z9{D;cYC2y)>m2?-v<&m)`o=uAP8}Xz{4m`7&SFcPWF{BMZ%sfqK(>A^tzS}j@AA2# z$4rN2xmp_?WS?Xxcp}estik$_V|D7~T#0%uyU1ylFBD3wYdmHWLrIOIwx?DfCqbvP zaFr?U4+>#WzsiPQr#uu|1KU4a&5iW;A21H56n-`_cvNELHM2T$;$(E+5y;^*5N5RR zO()OiZEN)N!`k4?-a;%IUJn0|M@qsy2?!Fbjzh{}H0t_(f#%7rZsM6-u|ttR*-)Nc zsFYthUvXLz$9XD<=c`qcitd&U&S=I-h9UhEnYKEk7IgtcH}0HP7{i)#>y33JH5vI$zwz9I%rty zQ=NNGv(Lto@yQPZ4=_NzUJv!+n87Qv?=vTOrJ0lc)h)A%MSkro?^hVm;v|opF4WHw zk(&mp>$_BTKUXc=S6+(FWRTHG?oRlrkMaKMm|^?mBl?Z&|4MQ4#NFL%y65v_yE$^k zs`o1wY-%D7m=%B0z;eIxX|@x3ogNogzA{;)3>+?)x}gk*S2Rz|SPjIo0IT&YZ?W*#s$T%+&zUeo9k%pm@g^0i=wW0-{ z1|xF}Yhf(X+(6s|y0u^Y=qc^XAJpE2z zw176k)4OJ=nhw65Eo=ieiH)9Q4e1W`i-1QTDMz%gFbvd4Yi5PY=UI*DzmX&Q-8BT3 zFGVx%71Zi`w*azl!MM>%u%z%O?6gvY zcljxwgZj6GD<%xt{JvozG7-x*rAROB`xlgL(J#Is3)iOoRRtLu$5UPxb%LUm_eED)$p7EcQa&2w~0(16yCN-n|!sFPRBj=vVJX(vlR_jM~+rc z&=#9gF17h}7$XnjWzzQpaj;8hqA){0?jr3>7}LiHJhN>m{z*msDIb6!GGJIVUHkg! z!22u}k+MT)%vFjAOf_~Z`svs;D9liq`#ER|y0`I48(8|~&C%yq8n4uQ@YZ2XU24z2bK^#* zg)%Mn17Xy6xl7Je*3kc_Jz*Ir1z zsls<_t2#jt62WL;9M)(NJbq5TZj647FakX7wJnmZP`8l!F25OA=xs;;#sCV!MtXokcTD`sr^Jr_|f9qKAY9rpj>BUoYl<_Z#0S=>MXdXK@a zq34Zqc%vJ$YV?9f8}qLQ@i9({s+-e6{N*%mWl1cNNuvFalK#Uw|D7=uJs0N8@2v?% z>o0xP{_`UPhu*ojC*E?-9_z!(Eo$A()2=&sxxJaC#?uxy8XwlZsiS*%+7*^z3+{{bukPmqgnIrD8o60s28V1 zRNrpy$gLwXXD+osAM-{ygv+u62cq(=HhjyOi<;k#w>@^x9ty}pyv+8vJ%)~TSMf6Y z2Id7`4u9vMGhQwmC9A8Y4$r#IBuG@IlY7;$mjdYGaZ-4LO6otS-olB+GPM*?Cynq4 z`vyX|>zQ9RCq%K{J%(M=VJ9!17+Kn!fWy2JsspR{8B^uP%a3(F!t*x|*byFj@t$ec zX4R}$&aoa86#@O}?wS4zeoxFoJIo^QN+khG*iyYE5_#avUP}_SM?Nv>?62qkt2#5i z6}9nCj*sG6*svtgn%d4%NbnZ*dDi=0?kL7LNA2e?vNlgz9or-~M{xZsQuxMp)LTCi zQPEo_FR(oD(}RKzx~ze$E!#7Z$kE3|E341-JA<_}Q<#(SPm2OP!SWH&-zxQyi;GQJ zovG&-YNA_=Ri*Kbx*KD^GW`UvmFM%d-tI82l^E@b`fP_MP<@f=y<)X2*=sM|>Xk;A z4{*9#z0_<^%Qo7Ivjs9Sq*@a(U|dz~8>>!~J5Sl{i`@O*JMOO^jIHeF)ajI=v#Z$y zVg{5>MvvYJSegp@rHkpFaKTYZlSi~s|5ooQl@V97>I@N4+xN4BN2|Yy*bZ%i8<@NZ z9Ef>JdA@@m)+w7G?K9o9yXE?_X|7gP-IK1|LCuPq9k7_u(j47cof0{d2Hdx$UbnI< z-0pig4!GtSi#-&8gPqDn@=7`Sq;p{S?!Nu(1}J*APtAJ!-%5h5KPfx>?;3 z(Pn$IdbiKImmka3!k)2bLzZ6*iX@>=56;xL`ydYzuU@g)#pplV}Mv2#jt-tC8g z$D7XjXFLz@T+z06?CEiFvJ5sYI2?4(eY!Jx{fR0tb|FVhT%y#dYTrZnD|lAKbdKD4 zd7CR;pX~C6k6*NSCvl2PGl<`JGcNK=;D=KjN+Zkg0KP$~(;RMH7UrycHr{_+d6zj& zlk>W@ub|K#{x&?vA^3vJoA#EdiT;-P;DW;5_$*h+jHK|8Zlw8EZ|+xO&#Yby7KctM zaHmQ&M)9mW`XN)3g!ZBgJ+t2Rv}QXljVfnSK*KDsS(vUPM&-+`SEMqB`AbdME)%_D zjLYQV`t;#wMv>7+F&NfylS()9C8> zQf~D+=>6>6f-7`qU-!TV-aZY^tF>3ES%lUsG29$w(LeulLgc5jHmXrd)CFTm5X$XN zcFEFidZt@)QW*bOQL`h==oZ||r)!FH@YcpVw+6TMCKaIdjlB9pR&*M+daF#GrQOA?aYOM zT>t`s5VEu|aYP`tus2`6ZLr1)uq%bDSbG~MQ@BPT@Cbx20*^=F;gX8*Wxz5l;1Rxf zJp54czP@-W6;EaP{`T`lz%p#3`fj!uz6>rEkD&hcqf%K^7MJ0RV88|zmBD2&xLg($ z&-$%qak;P!mf;Fk!_r2@kvRlH^pd5?Q720FOpZVl?WM$3zuH?(6Aa77ps&>t4gFe= z4|BO9Cw3bR@{C>R-XHyINIAt889J8e61ruJR=Uw?hTY|F#~`--ejT)+Fe={o|4}o! znrQs<(RJbo<+;$mA?x-=%^w`)|2Mb)n-MmMlinE*O8hqyScqL-pCiX=jrEG;Q{o?^ zr?2O-=jWkjS@qnwoE;_X=x_WRlt|sUxU#s5=iN_@3Weqq3FwAz*7N!?k2njv(INKf zp#IzJGfn*!;gcEQ!DnuV{{4$~fP_UCGu{=ej7~?o@^z3f5VEbO=jaWbj$8m}xH06F zJfEjMV0QbidPXskpu-cNrn< z{&8*L!YN5F)$B8WD7-OcL2Waq-~ehqMx)QelQJmu(P=i%8mQ@OX+gp|KDzRfb^h&6LsBY6d+3E9R zvTxtR$T5#sF|gMkn#;ckjsFvLzw_+MYUfwE7E-38#nq0}vu!gSnHoPgft*=Vv$wdK z{W}9}7@)7*jVtav+qvXvac)9w`Q1_OLV^9n1vh|WQj}z17!oe+`O-VPkox`|%Eb@= zY|Is%EtU+qqf@h9?s$#SGu~6yA+02Hgol7ttOc7{Yy&(g1Tl}oG?CD#M8vK@}3V55cSBVXuobRF=0S`JRq3DICYtk{t z&#LZlSx8g^EY)nu7yVkzN>y+a3C5-(s21#+@w)op5|@e%zFz|&QH3D5kf5;lUzH0{VEuhkE6Wu zO*v3RG)!$FEz2%1^ZfXGNp-00_$I(@63Fe!J5AO>9)H>YWbyopTb4j@8QBSP~i) zhy!v**4DnR56G3SEh#z7R+|Yy8LdT0?Cgt@Swk9t?E41;VCayS1=PK(<^{X)qNF3e z%Iv1wA(U>Q)6V!R!--@t=kR1TdyqgYCi>kv=M~mp+O_E;13k&piV`PbDXt11ci2FM zyZT-b8SE$isnFWF_qsN(R;t~Ah`v%q%1S_0Rh$FPi=NL8h_O0NpHl6A7GsTmc-dVb z;^*9ALBt!-*4n92%S`+_*uMDHP%$=6fCsU8dFYco% zN_VnXESnopg`R%U^ZG;?o~Qk}hnh9i{p&zaLHMWxy}q+gn zD3?iaxnHo`NJ^vxwE{W>G4)k5|$~G z1(u#^Y0U$37W;hoX&s&3XT83(v~pFLjR4-O! z`{`MXOPFWMZ9F6&jCI|gvHW-*+-MR)S*Q`T>jMF&uG0K%^y$vdmE(y*@ykxv$-i3uzYa5p4HE7-vhs? zZ3KWqsGiDx(P(`KenrH)c3N|qb4(NBd`Ie< zhaH9pAp>y+=RuY5bY}_Gv)jymQ=(6Ixk^xB)zV&l)V^d$vbOY%;TAjS(H3n{>G&ez zJ{;4L&gTfKp+pAG@m1!9Wz6rCVN7I5AX!a(x#DGqJ6>%Jn9mr|i`hf&IjxGcB_49r zFi4Zx-*TcbDYn9c+$}ECR2I6F28JU;=7Li1Vo2D!I=3)xjs+j13hW(<2!`f0j?;q( zG2R(CiVcSs^Z>eXl!n%#qjV(XtGJ9w1ZRA&S$-y>1^dBv$PwuQsD}v3@+L`=l)wfH zdGKp;oRx@m`n(|0r1tsa54TF`BKGRX+Q zyfpQuzuuVQXRILQ1IV#()R5*L((ep61X;y3?n_u^Y;;A4fT5j8 z!FMxqZS{#Mb=1^zRtlSy`j8kJ0@6}yT-9HJU7r4cH?~0=@>QzNBU(d_bJFFbTSa;v zW|c_;-l|bCRjXp8uJ=fk#vKVsQZ{T$X_4AoUAo(&XYxu*BlZ3~APY+N&c1Ty0ZQ1^ z$5Gn{8;zhnXZw)KO+SQ?8l{iWEWbhoHWM^{HOUt|UO(|i0m}nFj3D_I@$)ybm%6KEhyKB3aXSachZSjD+&;t5&ul1hOL#A#U8nK;cIS66lTz^ z$=x+XXi6+-AM=nye8>|wG15B}EupXcuWKRb%XLKE?bB)sY3$(Uets=+DA0MZiEWn< zL_6>JJR(L?zdA)0fYrXLlw$tR0^X6_655WMRR;@lm9x|l#ZM;>Z_B?mxAD5gd0L`} zk#jk8adK)z3MDEL@`!hekeDUbrTt*Qw+v)0)StUX3tL}2)?*CSn&t1-De|Z*KIEDU zG_T95R<>8Q_H;Tv@#6QEMM|J9(#(Ucq1=S@aKXk%K)7d5BufGYW>2<}IfCk`plNvk zn2OWGpmL!Rao0&Hlo3>j0vKo;NG$wHS3Zs@3COD1(~F5f@`FNAn19qz$H^u#Xlc3w z0uks3-DI0M5j`wZ1o*D8%oj0eVhTOGgvQqCyB83z>oL&NiM^a5wCZj1^Dx;ajXC5J zjiJ;_XdGo_Vmn$oEB-kn7i66suM)CZwb-*%rO(-`bHAuPJ%7L{Sed!AJ*)-iDlY_Z zwt3G!6P9}+CO$oM2A6XG)>ghmKi~ACNoET%&*{0ARb^QJ*CwEuB-&-@cqPJ}yl;fS?SlV+bNKt5VUf??>+E&z+E+#R^NeF7?i~V^<;{k)cF&nt`)RE(uk66aL z{DhOAXPVG^nKwC3OD(hGQQAiv(pusPgkFVAZ+Ljh0ht~zN}}w-$pz~aN2#_GAN9sB z+0iwP2`9rdIR{VG7t{8Pe`x8uKU0+(V(4z`#n^yKJs4svZ@-m2?+kF(L1?eCCrzc5 zxN3Z)pGPZ-RVhwafI_4Cwf*;f5RO~W&EdpPioMy(9747$ek8b^EBt1~8hg(v*ap?6 zd){DszbcWe^vS!|XzuqE1$}xgYC*2O7VPC;d%%vY3oi0#NkBm;F|hqFlg1BicZ#)c zREDY6eVU7bMh@PJ_@dRyb^{!in+mNq zY~gQH(EP@4U??{^TuLxB8NH2NK*s4c3B`J_Vg}5h6`}x-`FX*OE z9UrpgW$|J*z&63vZPjrADiL;fz@Fl=Z{1o?*1`*4ar7H+WR5y zPgW0CO}jHNOK#YEQg)mi(0wpN7ym|lZMnGqVT3X$C$iV?($1yIirR>^l?p0mt$)mM zwz9mPeeK=qJptv@v5{pN2Dw;f9s`w;^Y!kMkNQ=g@(v6ts!ALxy*4>Mm7hhD=r?sK zeDOBUib~7jtMq%D;pH~I(rIl8&b-ZR(U)kLijc1}5UzZ|?i^a154ko!Zu6#Z`w er~iB@FdrAk3ApN(E6&-x1+g@>GpRgA1pfmXe+}>e literal 0 HcmV?d00001 diff --git a/plot3.R b/plot3.R new file mode 100644 index 00000000000..94db1a29c9e --- /dev/null +++ b/plot3.R @@ -0,0 +1,27 @@ +#Course Project 1 +#loading and getting the data +data_raw <- read.table("household_power_consumption.txt", skip = 1,sep = ";", + na.strings = "?") +names(data_raw) <- c("Date","Time", "Global_active_power", + "Global_reactive_power", "Voltage", "Global_intensity", + "Sub_metering_1", "Sub_metering_2", "Sub_metering_3") +data_raw <- subset(data_raw, data_raw$Date=="1/2/2007" | data_raw$Date =="2/2/2007") + +#change format of time data +Sys.setlocale("LC_TIME", "English") +date_format <- as.Date(data_raw$Date, format = "%d/%m/%Y") +time_format <- strptime(data_raw$Time, format = "%H:%M:%S") +time_final <- as.POSIXct(paste(date_format, data_raw$Time)) + +#plotting(3) +plot(time_final, data_raw$Sub_metering_1, xlab = "", ylab = "Energy sub metering", + type = "l") +lines(time_final, data_raw$Sub_metering_2, col = "red") +lines(time_final, data_raw$Sub_metering_3, col = "blue") + +legend("topright", c("Sub_metering_1", "Sub_metering_2", "Sub_metering_3"), + col = c("black","red","blue"), lty = 1) + +#copy the graphic to the file +dev.copy(png, file = "plot3.png") +dev.off() \ No newline at end of file diff --git a/plot3.png b/plot3.png new file mode 100644 index 0000000000000000000000000000000000000000..2f5ba029536bc7ec7b86dbbd44b6862d096edd17 GIT binary patch literal 3893 zcmeHKX;f3!7LFyrMU;eD1j7uqg#_i5AwWn!G>8F;}K*$9v+$aOn7i%KiGqeAd^2CHZ?MtEy~0rGGQk&A~Tbm$PZUq4Rt#45lRPwOaQ7INiT} z-*L2o4A{7zWdQ(Jm!ltF@n*jA5IvIG&17W zg~zYO4Lm?n=?d}>nyr8#eV-FF=v_;K;`;D}JeV5)__gw_u&4Ds=$a2p7kP;{9|Jx^ z*`hkjZ?bu9mO8@6a_kcxf^~d>uF&>Z52Sl+XSFN8s5dg^DWx(9C^0GH&IpZ&uxX?x~N4WP7}oxsxV`ff3FetM60 zxMPpkSm8=!EHhkQD z%m6qyKa_Mk@TZ;1-iD2}`&jZpFNoCfWvf*imioqM5eF5&uGO3Tx6z`O#YkWDw9VC! zPE3OAv~5aAXX+W6Zfqy_xomJ^yq{N?@#$Bfit?Z}EbO=Zai7DfnN+t{_D>T+km|x% zRGwk6QiYYp0hV+SRzKe&4u6@aRE|}%`nr$hna#t`_gi>6Tx=_ujiJBG@UxclBISNJ zxB7f)D=FOkB+0u^Kiv=4AKM%J)R`<*B5zMqn@%tvNosg$fy~%`vCSmW&d?*icescvvdxIj&e3~CPwn7(%9WjPLBj2}cWYO`+O4Tvdu8 ze`cvDu*7&q;d@>vOW6=UvL+3qROr{f=WUs!tHnV%`dU$DaWYvhtseu82<+@rZbWRg z;hX0VJqS>>bn8vsJl_`NYZBHHT^A}(N+~5b=!C@3e_R*^=kD9_X|!A&8SZ8Wlc-Q$ zbcCROE=2Vd5|!MLSo{_bNS9h&}es>T=$6ZVSmzgEojEel-Hg2?$llR{j?_eJjkc^Yt3%r*2VMD*z z@QEmGVXGHC0&4Kqmxh$BrpuA`>D>K%UD$(_!PAt z2^*VR`_?#=$QSq0tNhk{GV09e$n4|A`jytp(`q}XBN?*5v{CAdZf^a=Oy4)mW%aN~ z0ec5oaWEnLqO6xYg>JiS5e-evbu==1=Iue{@d62mBfBO%ZqNV_dgk7UPoz$F3$D5@ z+)rTQgmQ7n4!F{2U}*~k5F{YLj+O0~cds{ePa$}~JL2tXD}C;j|N3xhx@vk){9Q^j z8x@0_4Jgl|x<)7sv`+b`$;Yj|Ozs`UR*mhFCd@zCILzN$+JV5lkxPt!lfjz=e5`GKUa$m?T#a5wr>c~6aj09?3l*QSvpf;DmN|bEA?~*S%{iqa5pDCjpWLJITxqm?-zff5WR!4vvYO-|j8;ol?D5XijAUy!J!v1^-gW{CFjq#efH zCNPl4=!H8}`rW0V;kYM=_1hLd`gA`k>XT^VkPh;R8atq-gHq>)cNFu^E1zAQq1;Qo zq9qJd0n>W#VJy$csv#gU{&Ntc1vMl*k^!A^LUIE>%U|~MV*V$9>&tiE^;;_6AWLUS zY2*88t1lG9;r%FoSj?UG330Ry#k8@q19f|dkO8y{lUL?9+UkO=pgOv+LNQ>Q&>3hy zgN@$%QJcM<>}QOD?1GP>F#nJP(*)V*j}xm)kDYHhz*3>7xh{No75;29+_z9&G+;9N zd2Rm(|9K%?B}-|`Fg+9```Fnk>Xnum0&I&ys{X}5C4w)&{3@Oy@cxmQjLd#buxc0o*R>D1X7uv>%kz*f5_%G-DU{$KwD9z1DQq8U|5bNXJD zNV`P44Pa%@DL}f<(zGgVE1i)!fN2U(Zp2Ps;G5X;dvRt#Q^6S))W-15aI#WU=N?uZ za5rrlphe>+nx5wGd34tUEMCd@u-I>$hLx^&imvNrq*=*9qL*dNcS}F`HhQ1^bODH} z9$-J+*zkH3Qkc?cP9rZ#bxGA#<4O@mvQVUx|7=4~lYL2! zBqZ6FqKs6?zJ#p5>GS*k`TX^Lp4aQSUuUlSzOUEVV zN~{`5WB^7+5=Zu`C_oj3RQ0cORU)Y>ky4eISXD)-VxSm_0c2qIs;c%zs`mCq5=kS8 zl#wdNlQ_aKG8!0RNJiRYpcrCgZ*Pw=nJNkc`q$&H;;$b=-K92=3jloauZltQqg%^`pYj#Pieqx2@o;kJ@-#LZk!7tH9t- z^;Wo_;QcW=9>>8Y*Bkoy0Fq&fyI+GRdMHBQ?2)%EaC9ulD}a_ff;n&X|G^-OAr@P- zA2I$q9giY-%zw{7QEw+tzlF`NDP`jqWelJ2qSyL%XM?{s7tEwc7V@FjdY58nzJ#m4 zsy*lIY2V4=H(kUn5ESC$dnm#7>l5W`7i9`?!xzTJbh2V21^OoGpFpuOfvzjgzNAKS z=6sm4w9@Ajcg5}Maj6evA-dz9l?!pZNn;z9K zSm%q5A_3;VIka=~IC3B&iF}JHf=_l)B<*oES{dIjgzM#ef`TO2Tp4jpUQXoK0bwGOe6UB|C?=&$xmYFe&L($##gb0)knTyGNG zKO(?Rymn#G9ar$pUgbNR^610n548I{|p__I^1@He%69c=jt)O>*blBQk_C)_m#dm^u%bHSF+OV4sExN27{^DWIMSa!T z2gXfB{5ma^HK~xV0RAD}2c21J`j1rcb;ABP>RgO-Yq;nW{*DIx=#M;$<9#mK$HX$F zQ20w=#3wIeDR=1Tg`Rh)4cbu%!+(fB&?Mq?Ejg0yrzjh9#h44K!iDV^g8S~2d+c$2NE zb&^0l&UqM5{J=SUu8_SZICQ|Yyu`L!L_SJ8*naY2tyva8?IA*?4BCbv1-hGo&Y5RcLmpW zkG5XZ!YVKiD+Q`()5kg8?0Z3@awF3Yg|l{W2-w-#4$zzqH;dAuu%m+gsQIEDZ@90) z@VzkY)?r;1+K-!|!n8${Tm@`5NM?94fvh*O__km`C)unr7kB}aGS z7oWDJhTRE@I)Os`ho5#Dtk#Axea8*grR3qr^XrL$w^5=)&pZ>)fMPa#-U$I;j-zR7 z5miT=rU1Fhhd8*2!>`}_I>3e@P@m{D3_ivsgZ%a&!qTem{I?N131rSTJ3FCH=n)^W zO{ewSskA-zcmVem4i}-F#;A`=$lU_g)HMF#XC?0xG2>Jb31beUihZ#wLiOU`Ss3Rb zoSAkJC~~`oqZ_xB#CyN+j-@$*HMeejFYg<%hxJ(z76ChRfaX1(ZsipE0DjY!M`VV?bj1U90#4zF5Sz*=7?3oX*4p@5GjX%ht7dHdXu zl!tP6yp>Ue`qZ&~M0H>!M&}@O4T(~tYj0X#VyuKan5_T(L67C1Pk_q=9~MtSB=%8z zUOmSsKUC3%s#zvjbTsfTRtxHAdS0Y7qVIV_3F~@c@#4}-UNSmOEaiCMi4SmgU^#ag zz<(JO-R?(M>ABD5ct-3g!PYDZiA|3>&{^=thvG;3J{X3o9!f?>jB}0=th5BKFZ}$bIwrmz{RZIYhWRva zfj-cEQJQppLF)~+yS88FRoepdKg-(aPDkM;zDgeDr$u41D9h$2ks}H*XMDfE6S^xQ z-#xEzo)>9N>PMMs{Xt_yAZ~X2i%KBv`BN@rtBqh4nSm$b!B}x%AJ5GkOXLJ4W}==4U4lJ%5`sqTkO>_Vv0xTVv7N2OOE9muq&BhEslIT$KP~1p%=m z9)?$A4aGO_P*};=^o47G=iLUX?#!?aXWO)B54=88i!N=I|1+1K#SAI5pib#X%7Zl) z!TN3?dEiXJG`6PP2|F7%S`(E(IQ6+znH@^cjRsu8%0!GmcOMXQ0SVf174f=1VYhRF z*CoW!s`}-!e)3JKbCQ}+)W@dN+!$Bu_ZNAi-@@T&q8`9sA3L$NF%)j?8~5pMeULud zz2T9oerCE-zaENhxXG+yNKa<|GCD!*jl?Dv8!JXW3e8-K$T(dA$M?*Emdz$$>?&Mj z;(QO_CnLr!8!B+g-+n`|uADv?Bzo}qTNI`~yXD^jdI0=0V4N8aQrAYJcy8|_lYxQ* zg*VM&B+t3V8F0Zf3iTo?P-1f>ISJ!?Y*@Z=oBJ^VY*s5i#6Zr&7k%v=_4FbaqFJ3B1uGt)vo z3}1BV;@?m5Jg`2@T!uBE!G_UIKUatv9&C<7icy(ramP5GdypwM68U>7Ui z5s8o>GZ$abT-B*NC~Gr~s`86M5spoLWTPg#5+JSz=ywJ#O32}~waZf)fUi{iT;=fb z>t}zvQ9iLf{jH2vtr@EKc=y5Z(y!C*gzhubCQXCvDG1r1#D9CajtTI0{~?tKYv(&H;j!u`WN(9r z#4T>^;rE<&N=G5v5_aOp91hwy9X{wq$+q4p*jXmuV=awowmdbrjfGlaC@40rhZ z6n>Ogt%wGSkM;4=RxK+YKUMZrEJL5VQ-eWF1Ps2(bp91XYcG1cMkg8@Z zrZ?%2&q58ar?Oe(oD-DXr)b)zBUxGSoC}ucRj)(eM&^<>JAQw^_mD3>moxjHw`4;s z{#5U+xupBG(x%?K$JzV1YoR}#kw)D1kDg%UzF&wR`rNSMQYO_j!2jc6#B`L^rEF=a zN4Fe3|H?wh+;jdW@J|28VAs0+cisoAh??7$&Rr&dsBwG(L1y zz3lU?s(`Pha`G5gU$ZOg^^;C|I&1?T6X+n3>t#0qB&2B^a3Mg0FqN+i3Vy0W=Au62SSnsZWnShQ9( z>dcS(JM-Z`@3yhylisT~jBC`ySZ6R55+$A+9-Yb4M`<^_*M0CrjeoJ@oO->iLx9(u zR`VlLQ~d?eJ;F_cW7|&>-ta@01L*umTV?xJ-=TthEJyWTvIJR_vX+g!L)Dk0)dEN1am<=dIH=RW7fRM9yn z^ylT4-5Jf=dRze(3Z-+I(K!Z|eyi;_k1@J3Vw2n5kOuiJk$J)CL;r87tRieT5K|}j z%C)Ul}^OeuoivX>49&Q4jg`E{v;iAF+m@3Wq*03FSdkhETuw-@_k1{IR(}(>N!xd+zRi+0GN)s=oiDOpsy!r#Ot@h-?*6CO0{U%Sr=gSTgW_I@cW_y{znH1Cgz8~S#0aja2C0ozqj(qk;q8(RyzOJa>Sskz69bZL3q#^OkI78|R}s1)Z=+Z6d%aO)FA+M%gCU5y{zKo}`8hXI9jugd zbNi=naOKwVXvTVWh=}Nj;;9yP=(XWo8J@+8^iWXPyN2oF@ zolj;5!E}+PNHww4gIM+hoha(fnd{Xl{(fh9sJGM@tfd&>5IQ9Ta`=r0Z*Tr+B%!Fx zIb;=J$amHB6Pxl?{4*wXdidVFMQ!f(_Z5Qc=Z!wEWphdRw^D!8I)6Kiw%IUn51seY z@V!PNm0^mr!s`p1wssY^XKo+6GId8}EM2&`kobUwUQ#vV&gpF1bf8gifrkzY>_&F& z{yd+0fEw3-Z96_RN2m-a4s?`uDfsd!doqLXUx=E**sgnmL|t1`=1Yg}%ewsVZ~4sI z4$F(C6_NaY?B!#d<(unks z+)Ci5#OE_HS7L^k{Jt~(UMySCYHyp0 zI5@~}#v3WO#H_Am!%f?=_(eVIRPX?4H^_2~HN@gDbQ&O!_+>xk@DaVS%@jA~w=p{t zoI&)|EfFBb-93e}sVi_uIz%`HCpmPZK;eMvVDr>f0)09+Z0Kly0m~IIP9!V>gOg(N z6K!E;l+(oNj6MD-K!V&-?=Y$dpjW2>)f94@;#CdQAZI{pc2Tt0yYBrlOmQj?#34s} zH`NmE2j@ifD_)-Uyr9#;MQcTcssIKLK5#ESoEvvt47!J5v9k4kf z5N=kn{0abi_9_a9yOp`0_WtAOa(mMWNiSws{B)zjv*hA1UxC(j@fQYO%vLWY4F6`>=7NSbstwJPK|0A}Svaomd8 z5;=5xCDOnmI?y>h4!U+3!w5i}m4y$BL*2J7OGJ8^nP`iTizX|f?m1$q0t4|N&Z#=*B|c% z6NFTSQ8n)aVj1Uq_fMz6nDA#|A`S?{&=Rctohw45VeFYz+5Dq75XH&C_w0)bw!7a( zOf-a{THV9k5N8q~=aZSe&Td`pi-E;~#Tw%X04MMG=UFsZ*zF2$ZI=KM`(a-J%Lz$L zsh8|g_HUbKb!SwUReT13v^{tC5Gmbtsj^n@lB@nhvv_J|!!O%6m>TDP5TglAp`)=@ z^ELVNwV~S8rW82-vYDK9xd)u5(a_Z{0^LV1UQp<_PNp%&#YPLFgcZ1$vf-Bj2Dwd8 z9XOF|8I6JrhFk~Xr82bW6{*{Zky>Tu^$ZD%Apy6A40a4olyMKVZgWqq%sYdPs1Bmh z?*w|%xJBNzJ8Pf^@AAIOwv<&y#1+i=nMSh-IbqdFnjt@!hjp(2>B2?%HX*#IVsod$ zlxzE^Spi)u4i0Oz=lhwFga?Gn=+e&f=#BngbPE}%^?X7W=pUs7ZJkfRBiaRzJ*iG% zGxQ~shhUAUK>V-xCpfZ{5Ihq?f4Y8#5uq2a8!=PAsKws?!nlx1=V!*X1Ca;)9$p4V zCriu$S*b{Yo*uSPXN_BBql^1`J5v%7yb`0e9i(nE2Sp~ zT?JiYd_|yrR~@e{{EQTsxD*H8ZErWTIc#wml@>qTr|GHIAowL-*hU3#5l&6x>GX*hnB(82p7t4Ah2hpGQaZ>KL=3{J!lz^ofUr%cOFW^4eOCaY3(v-A3b-u zv4n(v_yn_S3c*z%REqvBnzo_f9=rYlm`LTK^$Ye+oI-iSA_TT*ONuUje;?PKs{e%4cF#*N~GIM{^eD;%l6`3P=~Na9H0qc;~0X*QO}*?M@(H!1$)#vWb7)i(vLqOym0nspT%I z%Isbx7)mX1X`MMCo34I%DPO z>(Lx2NO8N}>_z3PS79--d-Rp1P&8WHdnWds%)eCW7K^j5Q{0WxF!=$NuYt=OlTS4_ z$JK7`uL3f^p{tA3Z2l!);A^3hIb{=0Xdep`+y)1GfsQf%lS+TT`?-|FIY8sf6`Y*do!69m$bEZ z1D2T7&K(nArWbU*;7E6lW^vB23F7##UZiOy*C|?k0e*@t@Wi!y#(-Uq_^E1sGCNNd zOYnN2)PdSq%gQT#IkU;@U_4*=_zF4s9;m#8)(hFh4YI}cu(YL!ii{&HD zKT?Nn1`?*5DTwk@CbRuLbHGnem4j9)a90tZVzpVN!zN8|5hvs6oa9>%e%OjBK=~ch z|09QA)cX9U;6mxv->xj6rB(wL+OW{zySV-TONIA1%^rkfyRtv;{Cx*|M&C%U^duDe EKeT(sqyPW_ literal 0 HcmV?d00001 From 689dbc1277c196c433eeb6f3afb2d87cb950d9ff Mon Sep 17 00:00:00 2001 From: pslcool <50282998+pslcool@users.noreply.github.com> Date: Tue, 14 Mar 2023 21:12:16 -0400 Subject: [PATCH 2/3] Create Course Project 2 --- Course Project 2 | 1 + 1 file changed, 1 insertion(+) create mode 100644 Course Project 2 diff --git a/Course Project 2 b/Course Project 2 new file mode 100644 index 00000000000..8b137891791 --- /dev/null +++ b/Course Project 2 @@ -0,0 +1 @@ + From 295e13a41de234ea82ec2208af7e9092633635f8 Mon Sep 17 00:00:00 2001 From: pslcool <50282998+pslcool@users.noreply.github.com> Date: Tue, 14 Mar 2023 21:28:55 -0400 Subject: [PATCH 3/3] Course-Project-2 --- Project2/plot1.R | 24 ++++++++++++++++++++++++ Project2/plot1.png | Bin 0 -> 4616 bytes Project2/plot2.R | 26 ++++++++++++++++++++++++++ Project2/plot2.png | Bin 0 -> 5511 bytes Project2/plot3.R | 24 ++++++++++++++++++++++++ Project2/plot3.png | Bin 0 -> 5696 bytes Project2/plot4.R | 28 ++++++++++++++++++++++++++++ Project2/plot4.png | Bin 0 -> 10180 bytes Project2/plot5.R | 27 +++++++++++++++++++++++++++ Project2/plot5.png | Bin 0 -> 4309 bytes Project2/plot6.R | 31 +++++++++++++++++++++++++++++++ Project2/plot6.png | Bin 0 -> 5839 bytes 12 files changed, 160 insertions(+) create mode 100644 Project2/plot1.R create mode 100644 Project2/plot1.png create mode 100644 Project2/plot2.R create mode 100644 Project2/plot2.png create mode 100644 Project2/plot3.R create mode 100644 Project2/plot3.png create mode 100644 Project2/plot4.R create mode 100644 Project2/plot4.png create mode 100644 Project2/plot5.R create mode 100644 Project2/plot5.png create mode 100644 Project2/plot6.R create mode 100644 Project2/plot6.png diff --git a/Project2/plot1.R b/Project2/plot1.R new file mode 100644 index 00000000000..ed4702110bc --- /dev/null +++ b/Project2/plot1.R @@ -0,0 +1,24 @@ +#Course Project 2 (plot1) +#loading packages +library("ggplot2") +library("dplyr") + +#read the file and take the data +NEI <- readRDS("summarySCC_PM25.rds") +SCC <- readRDS("Source_Classification_Code.rds") + +#calculate the total emissions +emission_year <- summarise(group_by(NEI, year), Emissions=sum(Emissions)) + +#using base plotting system show PM2.5 emission +plot_bar1 <- barplot(emission_year$Emissions/1000, main = "Total PM2.5 Emissions", + xlab = "Year", ylab = "PM2.5 Emissions in Kilotons", + names.arg = emission_year$year, col = "red", + ylim = c(0,8000)) +text(plot_bar1, round(emission_year$Emissions/1000), + label = round(emission_year$Emissions/1000), pos = 3, + cex = 1.2) + +#copy the graphic to the file +dev.copy(png, file = "plot1.png") +dev.off() diff --git a/Project2/plot1.png b/Project2/plot1.png new file mode 100644 index 0000000000000000000000000000000000000000..70fb96f37fa0fd81c39dbdebefd38257016d1fb8 GIT binary patch literal 4616 zcmeHLc{r4N`=3rT(=f?as5eUqIf#ldmQt21Az7kql_5*^b&GV0zxr1_atc=bf;f0L6^ovIrQ+))ux#zz6~%Kp-RtfdC zk$DF&!6A@H1QO_p2prI;$nR=M5I_Y^BoaqN0z@LvxHtp|_gw*p1mTbf9KggOa5ySY zgowitsW>XoxI`p~h(r*PNFokFBqE7KD!?HEbS{y|_e7 zK_EbxOT+;}Tq+>T69>rewgXo{AfXjgLj&80XLH%`jC`V`$ms5ZS23CKCoO^K3Dj3PI^7+%_gd$4qcitt*vR!31Rgg?wHTP7a=K$cAWmMw^6NQ z)zw)8WALZ1%G;nxAks|fg@Bh9dd_;UKVt_emUtV?MS_sALY-}BMUY+kD7iYzD);q* z_ALUI@rzD zzI6P&$bvga*HD|g1-aq42M!DMCa2Bj>GZfs~I+ok+Cn_H%vB{~*y6_SYK3!hzN@l`U>Yg-f%|a*bB$bmhryTQ3#jd8G z#w`9yRkZ>;a4(8$kPpb@?~5KboTRjk3oeV1;=N01=+{P{;8|#aTl=ocZB=4YoO0BA zbV!??0!cmmBdpMV>FWCM{=;?NcOHB6IpFmqS07}yRgV|DiK6y4GjC{MySVbha&r$H z2vTLyg$~t}nC9~BZlBO;MH8%&+)3x_iN30JTjS2@93sZ1guc|~zh2Mvm+Pk%79b;~@d(GU>tLW0gYS*_O zkyePtWi478eWF1gJ%_?4Gu=h7CckT>04%^jiF>bdQG%2bExRl3%o1A(;WJ3+OLQqw1P{Dbx4RNkANSl)_^Rzd%ra2E z>$E<`(*40*@N2LUJKNFkhxR z1q);%z#;P{%PpU>PYXSy8=Msovt(n2k7+@ZC8BzGH6S6+ch55Y`*7e#O2t*4>2EG> z`8&volaxm!W2TiFp<(mG> zK{c)u;Z{iW(gXC9m96>5w_XUIj=@87*L5vZD+lRwdEWOmD{_|_{4-n5hYFRS3T@U8 zG~As^O!UZ+JQ!Ru-m{S=eB1)>^z^ERT49a0bMA@j;-yNA)l|D4SlOU>$1fy}Q!NEA z4_wJ%yxwg<2)j%v&8ST-nh8@BABjee)X!&Rv0~eHCD|#M3mD;(I>)OLGQ9?jUitXY ztkk59t1*(uw6dlCo1)cCG0{^)yHj&U%b+?I{GCyQt&`Hl2Yf#d=HJ}R;#_8#E-C9Y z3=E(Y=gAo4IQ&$;@4#)2E#Ik`&9+@ujnXA9W-R~ZF^arj**#OwFGtp~T}K;E)Z4VV zRX3rc*VM72y7)fSBqThIR>rrtw0PrB4@|Spe*B*@z&TH04^zXOAlbDdS|!znZ;YY1%KblGwMatz&( zS@WD5iF|Lw6J7!3x1)~zmj&RmRnhwR;$ca zciF^N0Mka5jEO|#laSCTvpx4rJJ+O+>dDbG)-0Y@H|O2wB_c#BcE7xHZjQ)Z(;Z(; zd~Q11m`$`}#5YwA!`6xV?S0JT=j=(=m+f9|1kruJD))q$&o8D*b| zxUxY#w@At?g?kq8@NYAY!O%RXqm|R<)U(p?nNaBvms2SB!{$KlYM=YuM29JIUu|qY z|2!S48sd7YL;+g_bxpIPd=w0aGd_Gwt!K6XB!% z0v^@#kOn_WCt2W;u2&Zil&=W;?0F|E=P|yL6@^w?88NJ7 z1FTGqicK$mEGr?YaI!QW&nm1IeUs|(qga0wHuA9g7?$(nhRF1DZ)M{cij;@3s0>{c=}EzS*GoBoh>OlVkY zRFUy+&LrF8dLtZM+op3qXeY_ES9mrj$t|%8mXyXb1ik`zTCjp?IWmw+eD-nSm2YYI z9s25;rRifgJ@82Mtc?gMp%=wxKdG(Nb3AE9UjcEeY2)eueejl_IySX{&A6d=`o_s% zPrMIBi}?#_ZgEe_@R@*=5O710m7!X8>BWLev9vy6--}|TN4+T6NGIz^k$O-5`Ec>l zLgvL_`NI)=8JupTk~e9qk`JA9Y?a*3H)ReQmfQ@xDkrOD8pXZ7kSQOvo52BBpgIpF(~mGXNhCBz z1_Xr9FyiSCi-Q$R0pRLGNLi5;v1+_XHtRpJG*c}om%oRRozgGPAv1?3vy1{#R&ub? z!&EOohplp1@;O=Xh0wDjjpQy*@m5l6^x`(H1<1>M+jkxpPdaA z1Ojf^o!WgN_8tJj_B~U-KKdcRnSX{;b|_C7EBjJ*NKcX^<^fz)k_7lb1fF&08L%)b z6M~T=*u?C6nW?lN*0Moh|DMplDkh#s|8Mdyuxx_(Pdxp)TA9y78zr{2MdY$gl{ngX6}3Klbtvp~gR$_rg^h_0Jn&FHUF>(-CE~m)a#}8ZI}fqyU#PY@sH+=@``Iqwh82d}bw<%4!DZxUS`=O?AS^0W7Q?fqg|8MPRe;3lE?pB2c@zY^Wt9n7lDXZ(%x R7kPhlOpPoJ%g?%D{|5zD^6~%x literal 0 HcmV?d00001 diff --git a/Project2/plot2.R b/Project2/plot2.R new file mode 100644 index 00000000000..bc2d49d7bb0 --- /dev/null +++ b/Project2/plot2.R @@ -0,0 +1,26 @@ +#Course Project 2 (plot2) +#loading packages +library("ggplot2") +library("dplyr") + +#read the file and take the data +NEI <- readRDS("summarySCC_PM25.rds") +SCC <- readRDS("Source_Classification_Code.rds") + +#calculate emissions of Baltimore +emission_Balt <- summarise(group_by(subset(NEI, fips = "24510"), year), + Emissions = sum(Emissions)) + +#using base plotting system show PM2.5 emission of Baltimore +plot_bar2 <- barplot(emission_Balt$Emissions/1000, + main = "Total PM2.5 Emissions in Baltimore City", + xlab = "Year", ylab = "PM2.5 Emissions (Tons)", + names.arg = emission_Balt$year, col = "darkred", + ylim = c(0, 8000)) +text(plot_bar2, round(emission_Balt$Emissions/1000), + label = round(emission_Balt$Emissions/1000), pos = 3, + cex = 1.2) + +#copy the graphic to the file +dev.copy(png, file = "plot2.png") +dev.off() \ No newline at end of file diff --git a/Project2/plot2.png b/Project2/plot2.png new file mode 100644 index 0000000000000000000000000000000000000000..e8401eebfc1e0afab5e9a0ad28404efe1401d0fb GIT binary patch literal 5511 zcmbVQc|4R|`yXY(BqE`rERl5-B@K}tTa2-9*`{JJmXU1P8%dfPgODw$aqb zB83^oWM3m>XMVSyr}uq+pZAa7?|tsiHRn3_xxUwRpE>7R&V5H2qO=ddj=~@i$bl<5 z8paR^6Mgr`x)&VLg`x65wq=Mk)dV>Nf`C9U5Cj5(09z6S0}`1Jf)Ij0gdqr}Ut$DA z2oxhAm|fKF1$IaX3zMC%w^S+ZE2_1e+U=RW zw{BP0r<`&+p{Db#B4`_#+hJU)Q&Q+qR(GLBH(^oHUe&EQ=cW&?yE|{dnk0q`iLn_7 zJ+Jn8Y9`b|vsBKj!VD-p#Ajc%)bL@j;dw^qe7oA{xO6oyBjM&&^~1`wxInE_qZQHP z6LMsRL`ZOPrlrPac%+JK{oV0}^SB2>6%JhARZmzg=muw1Z`c)>8jG%N8ALE-*7Mh9 zHj;E&blq+am`MUK-df=}GV62puWWYdeWPL(16yWLF=pi`{ryiZ_S}RfI=#3Tw4M=S ztJGV&J$!q8TZe#~meMKCx=>92xiwp&*(d@3ZrLj>(Dj9NKPLr>1^oQ=e^h) zXu47#^e%RZPt$JA=A;kTd;v3e3PP54X5!M-`kDw=o6{$l2=ptCT`r3u zJ59Cfp5otQc~aw(3fAh7z#}~2;7*1wH#JR$^22!r%Wb1k&fSUtQ12R_f5i^b3{U98*AZtq@1UAJK*~=$rOQT{>HiZdQW3< z+l{z;bKf)QqKC*a8bGVKerMx~h4FvJ=x4ZP1BRP?0d=`c%0Dpc=&;o%_)YEy=0FPP zytFkp9}=KfCITBqN80I1*JO_h!1hIUob{!2{P6*aW_20#l+&XV{M)lUIs7 z8z-T(hoiv3$dg&18<|OGMZ@+@kTBd(?Ki@)ph$v071fb#)h2^5r$}2>``wnxCERY* zHGk-h70=nP7=9)A7px~+91=OS1RIn5I!F76u_UuFMiHw)|#U%L!6yx52%bKbh%iTmT6Ao{} za4WrVxUL!sf?C8Q6H)R+v-&fxDkhU3oanzc zjCYK@x8O|n&}P~8MNi}3-}&FQ!iW`x<)0ZS%fgw87#TiiCZqL?yia@9U#J) zBS*r4YZhXEube-K3yBq)1P|}7aJ%ByO;*{38p@5Kfwz3~av!1gH;n}?FW$Ny<$E`D z!!E{nPGdnjW&Pd9C@~7xj)C3L7|o|{7K875$0jrmx{J+pKCZgXCRQ`uxw1lZ^o-Df zh9*C0d%RG1Fg*1$RNBf+MDQCUfvff@^-c*jB}b>3$40o5YbRZNE?S*!%6`os1evCZGSU|ivy zGq);HJb(6(XXFoHQC;7CUF9P1I=u9iiT{a%QPHMKOLS*bq4?gB-zzg-u*d`*YlQaN9=Q`0cvvQsIuf3SO z0WgdCUbRnggb`C#dNHD|;$K!nA0GZZ#Wm+zK|7Hy5eyGr_ErbBncT#UO7&c25%s#f z(3;AxnIo=$4kSw|c~D%tew>Sw=wn)N;fB_L4kdVk?>&{ZQAhgC?GW>*^FSEo; zwZS%eD^ZXM_MNjG?dc8E3iy$Rx9#KG&7WDv6$jeOK9|YN4MN0Q&iG-zh41XiSB2%YnTuH>!#}?K;+Oc)o$s&W{=w0LP}ukA z*UNMI@=NJwQouNX)l zaPbWi(-vKGPr1;aBjxn!XxEuzGe?OL(k{?s&4z9B%#<+du`aQifOCE=9IXjfDq=x4;h)A4rC8z90$`&_Bv}UWF_gZ}Lc5 z!u>fK;A>GmnU3bQb#*3eJq!U8r>< zZ4wP+-Vvg~$oQ7MyhenOkXtZ(%P{6*vBFUxs=$RcLEcr-5x0#`z4W9tA{c&_7r<>( z-kN5k*S*iy1R1KWPJN4FTYA%wAGVC}<>5HARe*XWa+2C5jQwO4-I+F7iVu9O8N<3W zJ6am#|CA3|txFiHeQD-gvU=0G@JP_dz?$yT@OIqi0NA%A?E1-Wi8n8Qv&%qJyX6bq zXPmwu#`;6iPUp{Gzb8hq?t6VHyQ>xV&^uL&I5n<{Y`FuC(9UH^`j zXLJ{Iue_buqM-JY&^lT|30+xZpzZ+<$(k|}Lr5x@P7=Le3`L0?9@u_}k-hPG|MlZ1 zIS)pO3eNb6_ygdGA2Yv+XDVLb#g$xMdCv7nUSepk94?+Smh)kA4DqKj7YS`ffzO)9 zlN-80>H9{NOYrJQ&0!tTr&x7LZF&{c`itQ1_61cB9nIa3i>+FDc&c%0BfidST*9v! zdwS?mf7Yxf@!D98td;5rOj$M@Vq!~>wUz92({$N1E^&0TJ9L^emgWBMr>#F&#$u{t zXSa^kC6wiRMfJFTRZ>;?ZPBq@24+_gAE+2lPP%RMW~y&7=vo+Fo48U;U+MNWo3-+F zTPHrupvPntm5;ZxRJhB!x18UF2_8`gZfIM3d*+8oSGh~NyN4!h0D%t7(CFjtYTmfW z?|q9WKuDH0b${L;SMyi?ys2@Sr$iEnotc&2&&cEP+R@}*+Vc&^=ah?DT80agT`Lr8 zpNN-c?zw0gAh()gUM1;2jXy>kW#dtUJ@v2TB>wrYpsUgUSEq|@|B8QLy()L4&81EF zInL*S4?Wdf02VFfQu;Q_7cc98>UnvU@nCMhw)86+fv17DNqz9yc5$`-KttK;uQgA* z%f>V#t_(OYyAryb+cKxJ2Z!%rrNaalPjwXRY{FTtRA$;~E6MI{5^7uyz zQ|3|~zBD!;W@Vc69NM-(aqqmB$tw=OyyVXO6kJl8c-C;;h|*pO;H_0Z@@hvi@m8G$ z^Ng-f(o^p-#6TKWDa)n4XqIYlcs_9Oe9mfw0*zU>)F-#3yWVX$a;1ocmel66ChtG- z@a^VY{Emou?ffGJiKTU2K1c%n8fofTf(N~J&__59swfgN7wN)u2|vuo@eudYQCX+v}yuRa=Z#FgP+Q$ z--~8hKmoRDs`tR|DOwONZnRYrRvbNa&=$fYFpknZss=m0H1ZPK&i8To@*WZxm~pDLC2qI!P(}v|p5>l-)5;W1JQvc$7IG3A$Lh>byzi8G#r+6U zgigA8HVX3lB>Jz^_U{4#o=<;|9zm8KT~|r`_^2>t+jhdMXL+0=Qz513vZUY?oe|u` z!SDS}c^ELnD{%2$u#d9b>>He&dkww;VE)i^dPP$*F%~>niR)tu=eVsM#ILh3J=;FK zDJ5PcAz}>zrSqb#Wn1*@S`uw=92;oA^{jr>T2P zSPvRbul>s~=e@(J)YZJglr-B3z3sg6=9p0bS(MLy%K6{QDVlskHJLjWgR}D&dOtJ4 zP1B4#H0mq2YhTrM$vV%M?%gN2Nc2DL2J7ruzWk*Fj@>C83n?sdxt8`_Jup)6^d!lz zdGCtDo&WwFl83$XJsYnyJLL=1AS~{XOk&*G62Wvx-YLs$yRv}MamvFPQ@(cTmEThcOi2?yP+b6QxarE*2JPAw6<1B@icX#6XwG99r z2>|iaPR@~6AGjouB7H%~g4i80LBz)xDiTz#j(&29wNbG)J{`sL|MQIBPBOtbT+|z) UayStD_z8hr(L`w!Uq<8q2i1dI#sB~S literal 0 HcmV?d00001 diff --git a/Project2/plot3.R b/Project2/plot3.R new file mode 100644 index 00000000000..d25e8e34bb5 --- /dev/null +++ b/Project2/plot3.R @@ -0,0 +1,24 @@ +#Course Project 2 (plot3) +#loading packages +library("ggplot2") +library("dplyr") + +#read the file and take the data +NEI <- readRDS("summarySCC_PM25.rds") +SCC <- readRDS("Source_Classification_Code.rds") + +#calculate emissions of Baltimore +emission_Balt <- summarise(group_by(subset(NEI, fips = "24510"), type, year), + Emissions = sum(Emissions)) + +#using ggplot2 to plot +ggplot(emission_Balt, aes(x = factor(year), y = Emissions, fill = type, + label = round(Emissions))) + geom_bar(stat = "identity") + + facet_grid(. ~ type) + + ggtitle("Total PM2.5 Emissions in Baltimore City") + + xlab("Year")+ ylab("PM2.5 Emissions (Tons)") + + theme(plot.title = element_text(hjust = 0.7)) + +#copy the graphic to the file +dev.copy(png, file = "plot3.png") +dev.off() \ No newline at end of file diff --git a/Project2/plot3.png b/Project2/plot3.png new file mode 100644 index 0000000000000000000000000000000000000000..12e36fe20eb20733251308ca81e6f5c0afaab276 GIT binary patch literal 5696 zcmdrw2~<Ac_PWg+QPYC|m-KMxZeXw802E5Cie9qQ+G5VkrcQge!$Y zpiv|@E`>HoVbUmEFv%E5S{FJA4MSxJ15>ZMHe!cwbZs5~BxfTK}pv_URy z5M<#27T%ZzNx)BHfh-otXMtc<47REY;=`+p0J*U^E(JOWT!;-9g#KI*@so`Ce258)4}yaf=Rq1_kVY8%e#}~>rSshVoZe2-TKzq)FmF~@Au2pT?!?vY#kt}MMC zCN{I2DqR+EUrtfOVB5*VXQ&z5P3<&xapn?ABNVES6y_sFTz7p8dv+YrMftwed?XDS z^Pv>-*dNH9M5VL#wOrCZ?DZIg_rXi3nDs zi8AlGol8cV%R38yJG&qS=7sH5)RPsnXBXGORX?<(dVf*Ig1zgDHfrG+cq?| z;yk~##~}r+7TQqtB`2pUx@a8tX)y9J?Nx9?8#4~m$T!{Fk?OY_0&(>mB-!+wRuX(lu7zt4df0BH#C&J%2O)J zawuMUygNWAO2;uuEd>5UM!r_?M0N1UT>O2zH%m^fEC-k6g=kUN-mL4g;+^RU?CBXF zY*a!!HKOE4db}GsAMIEXZ{DoS$T8yZKh*7Yj&z&Ig^hSw)Vq>T*1s@0b8%pytLD}F zoHG_3-`tV*wCE7xh#Nc&o@nGecmGRln+iQHpC%opSK2N{Kkr(8T1vvMt~633NsO+*W!;*UaNT)Q)h%63kILda zfo!Zf3l=?8)lQmO#-16kSrWX`@hx3kC~&K1#N%DELjn~YIvU7WF!_r1&n+qaBojg4 z+Lfxl=+&X_DYL}}r@n+Lt&2L}>LIdsj$X2UKP}21xqn!@(fkmhw0yY0<1bV$l0G?P za20A)w9+RpU;BBfN~4;fVe}kjBG*%R*!S&szc1Aw&Ywwj&~Mwy%%J5xb3|N?76cx8 zU)AS5XC2vviyE*D;FqUfXeY!SPdR?L|2r#ZML^g6PkhG6q{NF2zuf+au3n*m3ob3J z&d0ccJJ;@@-8Sp{%|&ecEku!H?Zha^PEXAfxLG(T0>r6zK3Z*;89y&2UQ{_t3Mo6| z(iw=#Kd3Pq2{5Dza4M$>zgwkPQZe;v-VAQSa+#6 zb1+!ivp=~bxl558Y3*qJcTaZ6h%Ywz4~58`vjq?ht-u9qAob9{IaGk!T9iBtJZJwy zev8=!h6fj2@M8|FC>1#53?P5{{|=39DmyaksVGf+s3uQlZt$#^`rkq&-;h_+5H#OQ zFzoX1P1fwnYBPo3Yv$A_34=H?wIYF2L^`Y<45wpX%C-F^>W|>KOUvgTi7qL?h_-p=bI(+30JYwsVvERy8)o`(c zl({I#7@-wXJb<3qWAASd^n(5H>&$J14ccaSyx z27Q47O~85w89wwkKuvDp}qxGOV&&k1}Rl)E|p`wb! z0(XoqFjS<&iVS+(O}~PNb)<=;A{K>tH3pS|72@Af{G`HyrpoxNjJW4-xDs1BE3B~- zKT<|(hbS&H6%DUz3stfGCNBTU@zntrh$3cfkH|QXn4XhbNoQE3wur|cuk)Xw`PMl3 zN-Xf?>g~z!2=2Sz#~J9#Mzur6{`I1GhZ1iC*&@*kw7L@K;BRgyVdHWL65T6IDL7&Y zr4`)s?8QSBy|z#;ne(i#@aNobzEk%G_NiwNTvOuA35@~tM9;_kNdEO$ zT?XLh>u%urv=L#st6%q*_|yUe0R1ZR7{*h_eSb-3Yb-n=p=A&4IJrEcWSSiL1hNUS z0L*HW_1y-EMpqN-;x~9c5=hx^=fvF52PjWu(?qcu^3X*U`S~J}3k%~fu%|REUQTo} zbLu??7X}h_wqoNeTzUf$P_e9#h7-3opD|7N7rScV)~n=XvA=W?=@d zNAsz>gh**S`@s7dX2knVI{cq7t^R>Z~wIQ(JYex`T zq7CU?_R#@o*b1@J%o#7A^e(r->(}4;W*=Y-8aR05Ij;-QI0Kn>B}D~^DSzL7VqMqs(vV^2rQL@5D#&rbBkxOxBcJ54KL3%JS$Qx-Ho#C+8hG+H^Mn4~zNy5h5Nk;}8~muQN>v z47gD75?#J%Ij}7w?7a!u4EB+qadUSG3Md*8sGQ`m4WkycTF3AC0nJ#C-JUgUgVSWe z2v5?d92392hmsM$U>k7BRg&N@+SolskVivKZOjUmT41}zB0M!DuB=p&U;$wt;yYy} z&~5TWLa9==rq5lnldWvGKs%hcs>B&;#=~t|M2kwCClFK}A?$-Cu8kfeH)bA$6?}{X zi6$Q31)j3Z4^$*8O7aq*whuClHqNm%GpNg#Gs$0*Umf@f*l|xao1L-f6AT>k4@cTK zwxXtjs>S<`QfvBo0U7KOe@Uj?4HR+Xa1 z)7Do{9AsET6b2AJ|Xn+jybU<8?ApiI^R5NP((IC|Ems!lc`IO8VJ7si2JnMWR#ZRI4i87MU6-8HlgnjL<9FmJo; z(=RzxOq$OvN*-KoDk>SEwA)0Wkk5D6iTIBlE;T#NHNfhs<{!6x7|$%nRU_bGLum$p zT1%eIMLrGSwN&RNCnt0@coCpjV4o`(L5ClLO@Na2NI%dd$ z2Qcg$hNkJ$Ep07tZ||3p*9CN)XBzZi5t@3iAMWq^V>8G>Dk5aGxt?EtYym`}FKdmI z3&2d)_;Bgk<>m^K2YN1ayVsJ#oM8>xeE8e=v3yx!O7kP@*$?l$XS+HvT$|_2P9D|9 z6&1bB5ZA}ZD1JHf5L~_W?$VM2)NZTU+bRMVSL&RT2W(h(t}^VjxLC7H)`Lj1wY#rf z?TmKFblm=cvfqurc@9^2_$0VMy__-nUNF0&sS^^~!x%2E$ub%%2;3<_bF#gg=4CN& zuSnefq^6*L_&q;y29;kfP`N?UdB~KKULMf zhtn!<{xJI5POiLzk+w`QiGa4D>O}M|=O8pBFxkulLE-I9nAmq$>Me%lA0*%ByCL!=3(w=g=W z7eA3IjjLTK<&$%I&S`M#f@yQJ$cyHY1Fj~PSON~jAwQ>Z^M|;zkbItsx_qT`{9H;A zC)lFKT)Iu==jP@DZ1!4n*146LKYzss4#&KoNna%Ve9YP#bk$|wjhi@B?NEe*zuOiU**#={Wgoa;mh#`54KJitbr$8lt zY}YZORE}L-=Nu%TM%6zG>3lQ2dcc^p%Po2zBMuAc;DbE(9)4jUEv21fb4t1R{B6Un z!UB_qiaZxRJ*$rVBq6=a2qUa#52N8aUv}Pwde4JGn>do!p1z4noc6k-(dh^%=MNiL zwE9bvK3m=SE`b0Rba4ci`TQl!RS#RzdG2=o22XPk@2TJV(eqQK$YK-y1PRO(%xNg~ zt0ooY=lRv^h~CE590|qt&4k}abS2d0>>OW4@)fIvU<}R`?hjzq4JxB~@`R@@Ho+5 zom();EDhbuX*?aEj@G+swtq>qDHo>ebx22FHix>-t<(5ADip zjLUihypP{+>ss;ju;o4>Au%}OdwsfV72Ydzcg}kFb{t*)i91TE;Do@=Zd~o6w^h?k1BKP0nUWlzE-DLb$X)etLcW zjT)U7r)(*+^#we_lsg{kq-+(sg{p}{`R+#G#y!+i-Xcb1s>=_scexsj3iC6=HAS2% z!=0QJi;OO-)E2+IE{MNY=p?vw{iBheiaxm9r9R%HU8bv8M@)KAJ=RN|t#=Ff(#7&J z;^m!*^)K3fr|aHleqlT78pWhJ9cc^R&e`jz^gBqz2KCuil}BfL%ayr4p?JM;iVagMVf8R=csL8{)^8;TU>DsU8k7%l8}QYUe0+#Ua98|mx??;9w3JyH&tt_;4o#T;CFEQ4gTV0hdAn2oq9FM>9X#O&CJQ-nH zSb44DLr=+-J#bi8SASD97wS5{jc&y;C#I}>(N2$RJ34#a>De28)S$zurv9Y}y5L0L z;_N+JAw-$sJF9fxE}Ypz7`Ot;QO6GmKxhf7<^Ud8)u zJYRkGZU&C*^Ie&6Jiy-e2q!g-*t^lRMC4)Iz@nrO3tKV}Xj>$lF{(Yt!x*ZRH9QLC z)ffovQ!|}&6Q_-2j9blpDuM&(zSw;;Btlj5>)RyM*mApdR8B>+J8hnUd1Q4#qU@`N zYHQA+d+4uXllqTtd}d7zJ^F+I-HVHN7u4+xk&IM{Vw7K>G8r`R)+eObKd-gIlG+-whMhaH9sGo}HpBe$BIfc~`Sr$u)11+vhmA{tK!q^+F6^|JskH0wL z0{P=Tbh=-5`)*ZXW|G1@9}bZc#Ld2yh5)jtzy9Tj7~9%o)k(5g?TkyWAE8K)K7(}d z5d&hLq!5-=28rgmwMmni+L(_7rM|*@D-eHiJlV!6$qbwfzki-Kui<>xQ4F=G6RKE< zMd8z0ut3YnF@_s0BNX^RLF zd+^pHvRFE>3KDtkunW5tosQ+2Y#a0!ueWGxBV~CjW0<0H$y>X;5aQT~o*$IloT`&s zTeX;6NUYVqU0=ht;eVlvIDwUUz%LTk-TdNqc8!O*tk@p7XJw)rbG7TSCbIl-|N+M+;L;I(O8}3l=7S zVbIzJDYSI`_!{kL=8HX^{tF`{{D4OGab+u|J3}F~WgijEbavR-e9LEZ0!rZ}d}p4T zFrI47tv*24d9x@CwyiSRHg1V55tiJUx=ZNTBy`XHFL{M69HRA>ce$?Je+3-!401b;xPaD2ak!X^du`JpJ5p)>?Rjw}6bB$DYE+hWFyL zCE4t6%(;t4+>Lh!UnOtsW|K^?8lQ4R!D@^IsJm=U^{bxKKPZvY=j!T=jAn=sIaOTr zpI{O0fyBv>6{<}6$8RjydQg_;2w?-Rmr_W0xrloD``93P1ZZWG%iUyej7C}wqlqh7 zk_iP_lLsUSq6AMpN;n#vTlLM*^o{nChx9U@nu(jw$)}HkyN&}mvF?-E zd(F4f2Ltlm$x3~%9VQK(yX(#N7h!Acn19Pd{HCRfRr%;o5@fu0X3cKF5CQ9anyjnJ z#?ux~%9pb{X--Vgdx=on$K8zj`>MEH<|w}FI21Kll_RZS+o%g?Uu=Dfd zKHFA_<#o~=Y3NdLEqPlu>)h%@O3&yxbN!zScy|;sJ-8-uC`g$y{~=#>bN*~eWFAk2JVrgUOVX(sTWn{s&*-u* zb#E(BdQ)h<8dZ~wJ=<3Hyh%SDlQ7^4pNG43L&0xqubbc{{@o9!g@U*?tIOcs+FTMYKJs}I`q=Y{PUro; zdo{`njjj&c8%4~_<)}8Uk<&;q@|YubbW-%i$;M){&+agb48^n$JZn2usZg%}h3=fST149mXz<5~JTS0-Qsz}t1ffp2BcjU$_p>1t;={&2wz2N8X(lrX5 zRJA?LYm5+h=xnL6mb{nhecHB?w(fL!kvwvm3bTI;?ja{<;r8^Km)jd`O%&$pcQYox zdThhX$!Zp0ji`Y*o*G9r^+rx+W|{C_ub<4+o^QHEt(FIcF~~Wp?!Btovp74l_k{82 z^POcD8K{B6)gmW7b)PS1ga*I%lN!+A`tqDU^uOOIHEHMKFnTWerHl3 z3S1Y3`RgRai65Z{c8lofZ)D>y_)L23X8vwdhq80+PnCJAu6(|}Tl(tSp^mJ`_GJu`3hlN69hQI0|@}de;$PYC++&L-&ABgKD+E(`Fng^ zA^8r3@Oc&aSYez=l|V*9e7w=gmhCt7o0Y*F(?p~pk6Vps>`wzo>qLS7P6Ph?w(>vc zNjj;x;!>1EYT1vX8%#?kp|(nD!Z^TqIq{}WsebP$om!mOL#d;=X}iw0#O?^FhH}<$ zmh1VLS~nYzVVf?n$bh!aM@AmiGMoQ}gaA-mW5IXDbxE;(3sVIg>J;#3{sGSW3nyN)8}`;){# zq!_&*Zmq@LI}8R3W|g2l@ru+T1-(4X`%k8%Y6m$%0l+Jvz~%@yP-t+4c!uT{pnl&9l87c5l|VbC<{IFTP8}S1%{ZZa2TTB-Jert5(FwgC^ES(+#v;IM-9<( z)>g|)2Q)y^AajOx>wRTI-TT1V*D@`NniZw?GNq~L}A!_9(d(qBW#6TsHNInqwt@b%S zaJwV&0VwA<>fiv1F5pibJ%Sml{zZKD<%^&W!;K`VlmM|7myB07JFnm2aJX5e+={lZp_1aW5Y&jjhGZ6iqsnC7HWV6<$rn!F-IY zsI^vZs-?$dHZU+VBU}B{PnC-2C@Yl}B#o8vN*`I9kL4<^u$aUOc(YN!6fdY;nkU%# zP0&x~N+N||{u8REIDthzpl{1EyOGbc*k{JqNu%+#U#Y2lK?pu{BnngRYf7}^?w#I= zcVZB%lh|AV%pax2xTv^W?{&KoU#29T79%K^yK@r*S6t6M-0Me zSdeF=WLx~PC;;_2T!#)w`47Un_I%C~@G*MfUSSngv9sb#kmDgGhID+Y4hT;2HKWZ~58CaT16Ha+=`jmgunX+(wq<=@ zT~VvogavG%1UKni4Sg)98x^c^X!EjD4ln$~@zs3zpi~0%u!Bi@zheYmav8Fpp3{8% zBxIbQOTzp8_$3!Ql`=(_a$Wg$vZ)-LcR%9OK;Zlw_H4l|8r@qVbFOHGn#GznRY zinBQdkPuqn_k0ZI#d>6ny8$!(gcv@_{h2*JlR7xUEd0EUEF@r`dKFV+$PlYN+Ugl! z|GNE`l7(XaTt>J+@_7G(@dk07jFwT*bLa2r+mIGVDhVdU6SGF_fu=~K%)#6$S`yup ztzO~|Z#wIhVXHA5Z50Lu2O87{dzW=r5NB9=3cVVw_9kXk%khzzs84yg>byYgeINNc z#;5pDA)c67IHOZfVFMnm%r&=94TlQ&_v4M3eS%g`l0d zq}aTU2{Cq}1}>TkSaZJ?u{{|+^p{5u9o#J>YP`H69#!4#H2IL`zc?lT;68qGd1H0B zVZRk9SZ=?a3C9}}ZaQhh`XkcZq9a1*&X0Fqj_#q=qE3@@nF20G?5s@s#&73+Zk?1a zKKj+&wa4T3WKL#6NPwsNSK-c!h2#{zsoR>Hee&Q@FbArh56+kibLsHSt((2K(9PaG zmRRuJ((^+5Ces3T?R1$6tLweCb)H#0kjxEr_uyNTT7EF6h7S!{`&cXimYjv9#PGHs zuyH!swp8eOnx1~^JAbrh)#cYOpNWsH*d9=lyErCvqd3XKlz-08rYa^gzdE0UHP$E9 z-+~1ZOr*rOX9rQtAbc5$GgjF{D46|>3>2=+YWcPFRh5ovzr#H`^tPj~()|iG-3=SN zcpMP(QLi>pI%?zC%eH3NQ#LZ1&}WO(-}d*>KQu9QKzACryZkBqENNG{AMrZm6A-Ar z>9sd;Ho8V@WGu%Hrc#x|_RACo0QzuGTQ97Gvhv^7c7oBP`l?Go z*as^8p@Eh>h*5fNSL{HoSu(wvhTG}Q7U87D*U!0nLH(ev;<4_VWg=i#;WV(CmM}GS zJXp~JXX=?>r<};=4odg)wbGC59vLbJbxjk^iIL_2q8&}|4bd3(dl{}U@kNRi+Fqx@49?I6F5JHAg-HMHC3t3 zkhb8Yo;N7XuDgYk!Sj_b%Hc=lr=<9Xr}!Zo%163xz*%?a-3gZ2TqGL z7?F&EjBj?c(S_0H$&bkrp!qrXOFpCK1xjrEI{UuI*+HbJl@~aoyjlB<=>hw-g4I1* z8A|rDUIxRa(%OJRyLOiD)9Ld@(5){KlO86M%J7i<&CjJN5fpbYjhlVl$6b7x?NY|?3_F3Ob>i1C9abwD`k?p!jk}-Yw zF8D1yu#>}^CAtDe#gFKjbLn6*~9f8aP&J>Rn;S|KDIxQg_8EA32JNDR@eWI-ltUusbmJbTezh{ z>yxaNM4Jsl*JCwQQqw)WA_a>~fgy8;`Ht44cfd|z5nhteb)2chGrysDk^DQWnH6~7 zXz)o|G{(}&W?Y5q%cl3qH6Msl9e0#J62b{I+4Sv~>F?km1?+TCC82h&rIJ%Em0Wh* zp4z^QpWXqcLeK1O354Ud&}Lt9kbXA6N2hx(oVf9=C;U+5WQk1SaxS2eHbltt^Y!bjqkV%QA6)yxzX5h9AO4 z*yR_iJ3e%vL4$x#ymn7qZv2QJN2h;&j>MPrlsQ_b5BIRlz(x7v7%_)_i^TdB)O1V$lAIgN-jVt`w`yEw}`AtJAAkd*q=V3^NK@gbk-;1ORGLhjsgN^@;}?TgU9zM=b@pZC12qw1^GGc$@r>_R z24G6x?<~X4=pcV{i5Tp8+3A=dSqnketps!N@2-SmDn^rYf*N@(j=0g@pI?JD<<~ik z=?yWbl3WiOc5Ag{xI zmtw-IB)xbZpe{cesDU5yq28jdc-uGNUDYd`DQ)Pn9K0xdqg z3l?q+^7MdNQ^$vF`pEkReD|`iwu!?H>@4UN-TWXx0`@1W?O9MYF!ZD+dd2)JL!eN;Iwr#{cqjZ}3SOXJr5)2@lU&XIFYBAt zy?LGlY``88B3krKI?)N2h?>qLX`2zsKyQ02Y!<~l54#a|uORU1lR zPA!rX@Je#inaqeOQ9~{MPuhhGcWA#Pw%!gjHU0doU0EEjG*&(Oyxu6f!3Qxe5TaAH zrkhAv+9R;H^2;;!W;cy^2F_G$G4{tqy1;$i`n9-eVgtQCWtWt6LEhwlME8$+u#jY5 zdCVKX@M+JSifB>tGBa(UpI`tOP5n!dI8faide?!ZO!Ylh?k_SZ4q)m%Lf>H` z$@npwDsHJBr;MQ%ERYm!vpn|!Jip12x;W@Ap}|j*%4`A5x4V zQYpG;U;B|(14zeaY@1nk?hr8wGrBuYGczC*rz`3oEBMYq*{MctIH$g(k&;hW*(JCj z*Fl97OJLs#3KVt3Pcs<4%d?zk`W!*3s?J}Qf(Nw>U1g{NX+^S_(6HE@FLYts` zXb8wV!v81&L=bE9$`^FR$Ye$uCO4ZT5Tax+(jU`DDUng{&=4sR#@C1wU%lP(6z*;c z&jlM!03ZywngtkpDYhlWBZ3e&9wK z5g0%}Yu@Z-LSy4!akIE;X{{lJXfl?getJ!E}Uk+am>i#`IF#0He zfm?$J;A#8M1QIk-LiZMsL8;fAF9KxT`{#tfe_g8jpQG>m8EyDi!}Bz&epfM&D&H_O zQLhYp>-mW|@^{Ep+)W&WWX9n%R#uMv)YqqQc=ace3+=8jO%GWxk6Q{a=|AQr{`H~% zW~k!-%bUhG1Fv%LkDQj%%6AR$99;@bOKhlxyMv|XL01?A$P3Ha_Fa9dshR3)jb46B z8k4B~wxe}PHY9XsrmENb_+Y&^5^4iLywYtU2ob81UgSzD%-Gb($ab@!DBs@%jePJ2S{4gi=qDsUswUwV>XY@l86)>e{Q7`)4 z+Ig@&cA-`#2$f$X-hD z%y5fOY1rC^?BsyGX-AVdx>t(3=FL{m-ZgzN@pZvoUSP}`zHit0uZ}2Ice|c0bJttjw(Q+!QxEokgLBUfGFdL>z+kHES4h*g}PvGqxGE zveNLOK!S1d3v~4r_5~j%pLvf zNh+tx(h1KaaH>g`k^yhfouFZ1^-ts|_;QOlSXv$I2*s=}++Z$0T0$(f{M}!7+p9Ep z+t3+qJ?_p)n#o*3IoFQsmK4y=E(d}W2IO>?ZO`y(e8C_9g{S{pSHOrHISI~Rt&;>u zqFT)D!@mbL;zPCmo?+S}YVjn3RDQL#H2$DzCGW+(vyFc84G-zdwFd z{Z=nCn!IyUt*b7?Bq;|5V$pliF-v}+7)`#f@=AJoW(HUNr@iMhIdE&N);V~GoOnK8 z$!LqHu7r0`zZjw6Xs5jDv4DdMe-E<#GkLeLc)v#yEwok}>%r+fKt)kQp;+EL=zjrk C?pK5W literal 0 HcmV?d00001 diff --git a/Project2/plot5.R b/Project2/plot5.R new file mode 100644 index 00000000000..75a035268ce --- /dev/null +++ b/Project2/plot5.R @@ -0,0 +1,27 @@ +#Course Project 2 (plot5) +#loading packages +library("ggplot2") +library("dplyr") + +#read the file and take the data +NEI <- readRDS("summarySCC_PM25.rds") +SCC_data <- readRDS("Source_Classification_Code.rds") + +#calculate emissions from motor +data_motor <- SCC_data[grepl("Vehicle", SCC_data$SCC.Level.Two), ] +ssc_motor <- unique(data_motor$SCC) +emission_motor <- NEI[(NEI$SCC %in% ssc_motor), ] +motor_y <- emission_motor %>% filter(fips == "24510") %>% group_by(year) %>% + summarise(Emissions = sum(Emissions)) + +#using ggplot2 to plot +ggplot(motor_y, aes(factor(year), Emissions, label = round(Emissions))) + + geom_bar(stat = "identity", fill = "red") + + ggtitle("Total Motor Emissions in Baltimore City") + + xlab("Year") + ylab("PM2.5 Emissions (Tones)") + + ylim(c(0, 450)) + theme_classic()+ geom_text(size = 5, vjust = -1) + + theme(plot.title = element_text(hjust = 0.5)) + +#copy the graphic to the file +dev.copy(png, file = "plot5.png") +dev.off() \ No newline at end of file diff --git a/Project2/plot5.png b/Project2/plot5.png new file mode 100644 index 0000000000000000000000000000000000000000..9e30dac2ed811d4ee7c93f12151ba78d32d39e04 GIT binary patch literal 4309 zcmeHLXIN8Nx6TMb4FXCr2oXi;BH+*i7zbq#lpXcj z>FL;XCYA4J@Y`jd zKh-J*8MPYevuZ`I4Oq5cKg$tqFulFo6|Gm+gQ2`3$2v3sNqt z-1&)?I6@<_-Yar+_dDAvP2d{qi$8TP{%xI92IawgChny1k6$A?X3t7jVD(?WiwP>z@lPy6L$auS11l>>{hO_HzEg2FOr zHtrK~$n!qkcPKvaM*FZx2WrCVZ?~u&<224>S*h=GZiq2CZjSQ3Ye)0S{NB$^_;b5- z&QgHoqCJI1^>}!QS2w>-2QL@Ct12?|jL%i2a8+FN$l+DU|3X^B*wY5rg=ahMLr*0F za&wX^Sd^S>R#t!PGP`-d&3<@XaHTZ0p@i?pGAv4LqK_HhE-`B@+myn*y8sTB*hJ;3 zxxXM+_AGpU9oL(7GAwNL#2GCUpTvqX{}iEP)oG7dcduh~p!MbHF8BMM5(o3=gkOj8 z<*c8-wRL&Y-SE`+#tVn5FR_>p9ocWk@<4*Cl6ZGkZGP0#h{owhBece;Fn2;+zjJ6y z)G+$QBYKTHIM}RRNk1lWvEw&$3Vr6n|!`R!~DVF>0HSa|7j#lAaYg|*8Up6d81WYWq zo;{u6b!=x+`w|$f8(BKO@thWE40t;is%`Q~uAyi9`9YTV z*PxaRywc_mWQSaT>X>C3@e6h9(#+mMa$PEBd@`WXV@YNUnSe|qUfyG3 zEJ7@WNH6dW)sWA7O#GJiyQ79Y>TkM#mzY;YT$6taWcCBVy`$1Dype)kygEE~WrT3= zDl`;i0bHDWaM3ew&Qz*~GTb3LiIoaw+88)Jqo@~(%Em`?!Omn*Mt{!6+~LS(jTn|8 zoW5YbO*N16aqa;!XQnXY!X8@_`nbU1tV+^m&ELRjB17vEMzVXSTD4e5tISRg&qejG zh52`D^oQ6j0|#|gf>GBt{*i!|WjWfy*)1bAk%LbyZi-^O$;4bsHGuVmF#cQpVoM3pi!?`JT^zZGR`L~ zhYy8sgUa9(Rx?miRFSWJT89X^cO!xPHL7y39;(5pzILgPI=VPv4#5V~#MEWZt^Dqz zR4rC?tCpU!c7lS}LtXYJP?`fxt{KKWX;awpYk>q{L2*yktS_xXN_3OvNmFiqQ1wr$ z2McY;FMh;ZY@W95uyPD_X1u&Hnds`_*?pL9r%a`<;7Hs2bQMra8%Y_MjDRgY{`3I$ z1-MNqRvm3nm~8z`9mGE#L2P~7`D1ol5B-!*W;JRmAs<`yPcfqlPCryd3S& z7`yEzO={vO?+o}l0?NVq$Dyx+NP-OS;g`R=h522aTQ``M3osD}mJ!**7xxZ^eA6AN`N(JW-gi@GXMKM9 zMNHYU-+r$9%9Z&-i#}utc2WE2isJG_{5}kz+bC0PbcmJ9g*PVNp`I%E?jGi&Q%aAv zrWp81<;yJNBj$CkpXaG4l!l;g#P7VH*k*R^mM&`-(1QtaRE$XqjT~#qlQX-_n<;@X388~t=cVeJDPGQm*KeWUTlNgPLA>h2e)PWIVd z4GCYdJmE_gN|1UW-#&#X-Z)*}p&CtmcV%&$!~9$G$lW(MaA~?Mur;>V6nNI4_f26n zrs0Mx>u7jVW~)+zwMDsV`nhKB`HRhgTi%H3dgvq*b2q?A&R;Ni^*f|v zUC(qf-_S4r+3!7dW1FVmDc1XzqAvVQpHnCr_7afq?sZU$VU>>_3Jf$|EPh2ljSs|ap^Z!hYCbXXrg+FCBU zH#vfv+Sc&-wOWsZrb#H-yuta-oSGj(GIy_o6X_;v8*e6BSHmQbY<^O4W}ZF#ZXdch zoa?81%D`F)0eM!RwF==|MLU@;s*JTC1AcqRspT>|yi&D9x;KGQ@j{g}RiUSXjO#5lb-gO)86nav7|pqzJaBK*P3X&N zZI!B@)TJndcs&%Uw7A6@xGer^?|mqS%WkF6uIhG$nHHS9uGVGZz-{MLkE#*Z6}0iG zAeTA_*ypp~N9ULGr27*tdlVBo609h@VrDmU;kpW0H-Es7S0%W~Ln_CrKzr3E3tf1r zU@xNas$^1Wp6$z`;a^3{RYfWjaLSv_k4txz9*HVED`DR|l><1oF@ zbP?kJFREjco7x-BH2;Xk8L`!;b z$F^%H)dc*Fr%7yz$E(CM$4mMOIFZ~8q*-cShFY>#)wuUhu!;5L=)e>rUqp8h-Z)NYr zppNH2-R-X}?q7+{CwaC6x0RTY@SKJ5_&m=0>8r(8Ck#HmD&Tx?H#n+yg7jxctRhy5 zP4e*}y2hK6_X}#olYV+c5e)GTL@nlKzz(RI>Q+*`PdP#@HaQ|*pP{CF21j?_D3k^; z3pJx>0Z@8(k;*Ik3b#R=GMp_8mD~KGLoK-w5EH6h=h|P_dv>%fTPn%36zkeuTFvoC zt}!tI8G5F5^A&4&mxX4RD5Lo665n|PyAMbu=21UrtrT&HH&XYxqEBM#*q+&*5>Yi9 z4p?Hyc)%?vj3`b6<7?y{EQ{j>pZ2uec;m2&6J+~(UAvtpcZ$X#ijUCs?06?D2h13I zAd7tIu5)dq4bWNRAJc^^N}+Qv31YVXE?rUlB%qBOOCiS#IIf z0K=Ojgzq17r%HxD%U$e}#kG%wq% filter(fips == "24510" | fips == "06037") %>% + group_by(fips, year) %>% summarise(Emissions = sum(Emissions)) +compare_year <- mutate(compare_year, + Unit = ifelse(fips == "24510", "Baltimore City", + ifelse(fips == "06037", "Los Angeles County"))) + +#using ggplot2 to plot +ggplot(compare_year, aes(factor(year), Emissions, + fill = Unit, label = round(Emissions))) + + geom_bar(stat = "identity") + facet_grid(. ~ Unit) + + ggtitle("Total Motor Vehicle Emissions") + + xlab("Year") + ylab("PM 2.5 Emissions in Tons") + + theme(plot.title = element_text(hjust = 0.5)) + ylim(c(0, 8000)) + + theme_classic() + geom_text(size = 4, vjust = -1) + +#copy the graphic to the file +dev.copy(png, file = "plot6.png") +dev.off() \ No newline at end of file diff --git a/Project2/plot6.png b/Project2/plot6.png new file mode 100644 index 0000000000000000000000000000000000000000..9a71275df6ce830fc083fd4af2ab4bf6e39e3300 GIT binary patch literal 5839 zcmc&&c|4R|`yZLf7|Gb1u_u&$8B3#V*%?wNk`%XD+W5x|b z!ayJpw}rX!84!q#xqCs`Sv{68n=)1-{j~Kt6IK%hLW4j!5E>0aXM%7z5RL+BdW%4q zBM^9mc`gE9k3=FhG&ERDEDkgXx68O&SPz+K91e|R)f99lix?~!^WBnbj>nthS#`a6 z?l5MzgkUfkn2-=m2nmxx$L-C>nhCcD$AZA2Gk1|#B{LJ3$-q$nmOuz5BP4_rl0gq4 zkwQol_&tU2h+SyBdG6i~kFN(ZcA+yeLP#0AH+lwXB7+WOqCuH^`ekOaRLjg{0J-Lf zT$X0T04W4RngA#`5Cw;(;BXWQj>5>K02B(qpa4KU9#M}suV(?*v#{&y0f3%?p=X59 zSE=CgF3L${V~8X-_1V) z^}Z`(Wcp)eqkEv-rk`Go?x$cl9L4<+0!o7r**X3bm8CLwa0$-_>^)vaAKf=MA3<&o zO+|>Keg;<^U#fNFX)8M9F*~T6A(wg|bWAWXj9@kOy)3(O>s0c{5s#x;Efy2W56%kR zg&&9v_jFwB2hwLS3)}w0#82L@Q2S_FYJAoO-7mgllQjo8gY^pTGHq)!=*-r#RZr#S zh3uWV1!y@H|J}q|Gcan?QQseh-coW3SNxum`&^T6fYQ`7WkzxY*c4 zz400wz21wq$b2Eo+xiHBZ{Ea?GwWIR%7n z_5E(Sv_C57bfXe`h7yMmzV9MiYJCbmR4cCpB6msBlu+sjBE3* za+NZ^);8UVECS|k1ut*O@z7V|9`=@(zc;!IBkc4Q>E-c?+z#4u9{w#!Wur;r5BbH8 zs&|x|+R8=2MX}p~RXZW${TwQK#rB)?HynZK(%;7}kd{BhMg!JQo%>*qRr0Dfeqvm} zYaf=@mc#>n-D4NGrWEuUS9&Z@>8EOkBb?Rkn86pQZR4L>4!ADvXB?{0EGbCr$I5SX z7D;%N^4y#){iNr${#NLO&<#o3Ui8YCV~#y8`r#^MMX>M1s}25G1ILxsBK_$Wosla` zz2Yj}h+Mr6dUndbaGGMJTi0tyXtmNQ{Vfq`8sl21%%?E@*qo6DTm5~xKG|E}uKLjz z7EZ=icLY>GrxZ$GtO%ASU8;J5=Y^n0#vdPNOlkX@-ceLi-TWO!i*9jcqXjXBiv$Yo zp;_LqM4uWNE%s^0oiE?#SV@~r*t~Mh(b+v=R@xou9#yTMDPzSIqtE}uo~6Qy9m;K| zFr+|O!mS5f=7xMKCqm6(Z9mF{G^AZ|qXq@3PHlILq0%-TwFt#&U+`Yq+9QMGk4Iw{ zTPacCd%xy77@Rju>z_W>^?>j(m9UH^q_GwVRgQ6BcHV;3Tf0v9XKBv`$9 z9SVaiZluvjNp6R;WKn}>A_HGL2I@4nZ+>kLfC(vGbV=}!Elq7vmDgd7 z|NgxA@VaHAEaXi9%w6hlds7#$(*63<&k7I?Vi3msi*DSNS+GJPJ_1Sd2C86@q>~V0 zaN_??)lz}Q6L@wyfzhjV`sf$Qat4c{#Fn2Y-*Ubmll5Hx8Kobyk2noiv>j4I8@h<5 ztPnQ9;CE7W>pQZ6cl!7rEe#+rHC|6(IA4xK_auYE?%hk~B#L?2hb23D;E6Rmf`t1@ z(1>&lZCUgVRl~-1iH(-=>{wnbX6#jHm#yNxd*Z#)l9$D4njoUopz8cFc#y(BSLr%& zPRS}cpQB(>gARj&%kY~zPnvT}_$eBayY1iFnslDz=$ALd4BcVY$3WlkEAlK&flz5BoL&+Aax2pJk>S zhKQLCBnZtDTKeNh%zE$eGZX(o1xcref|7@b+|A1ISO_!;+f3$?bYk1H*Z&v`y-=o`W)orAhVZn#icQR_9{m>b zc|KIkKQR_U6g2H&9mFfW>0GQccLTzsPmO$)h7_&DW}}p2BVOssz-T{&f^>HkrWY=8 zIaihN3|uMr#VOw9cfajT<31Yfr(F$U&{5J!Dtw@zF~{o!Ymc8~_A|v+X@*=B@d!4s zf_IOT;C=V2QguDVAsJ_j3!GP;ua1-wh(C`aH;bDE(~+We*j3s?E{fPLx4~inbYRGi zlLnthLLJN@D`K|??GcsO4_Igt6UpUnFbIJ^0ny+=hyLX5fVMvAAVa(yd}fGhB=sO# zo;{)t+f{MEmHcDZnw~YJW$}P8I#C$DC$M|!+(0xKoT^U5BOj9?<0#8cc$e0qpY(dI z-j&F&YZNw{k_j7CDGm8s++j!=VdAF8kUz=-_z{Whi&OAsD>Q9SVUm(-p3R4D<5M&T zqX!``c&V3++fG;srL*sNa<#GnK2Y5G;~{)N1pOQ-Fo>4#8O1XH{F|1D~ybaQ4`7=XZq5zm{z_k?lg`o*v~-vGS|q zl~mn4ExWYXtVxcY($Wh8J?+s^)8G15+?K=i&!;EHKJFpLg_|o(yJ)A3nQEkDR|QqZ z86_Jxv%9>=`8Y5OPjSRM*G5$cKZ}TLl|5FNb7^_}e3|X{?3LjlLvvC;k!iv?HGcln zOsklCLOIjQ44UuL+(g`wCaA70M1vp*oi)hTy`%5`kf0N(T6M=hDYoJ!Yl?A0+jOF z=P@uLO&!ZX<%PIpj>9JXaVNr&HcLl;jS2l?K88(|_|^Au1Ho4Mh=5)!QTAbEXk_|L|^K~(fQ@`%x7xdWK==jKdtfm%FQv& zx=#`Gf9OKapF26~NvPlei+H0PAKSfd3C%vIS_yKOVp&P$^-<(Mxl$?Z)C`g8q;a^b zFKvF2UAB=!G4C5k%Q_PaS#bR|zK**u3?c?9H(6um{8z)F)bPl@WP_@d2ayE#W8Cx@!h% zPu~95uG~}?t*V0dqXiiWRZ56H$G&-2a!sKqQ85o|X7&fsS)1hzpD>*Y1r1c@l~~8# zZZ!b}Ml7YZ8JE+1Jd5^1?mNO_9(6I*;fDXhF6l&b4HNCf9rx#ZKX-jooE7U@iT=Z4oI>j-SgU~|i&yZtf!&yUIf+w~X#h=~_+LCF-nn6>~aUYPSeMkw|L?5VZOxz{F(t<-bv-|R!K%Ak5vI{=zmYKFCOa3bnP)~T)#C%b0;WzCv z($s;pj^;hygtTeN7YmTuFtEC@7>j$XAZV zJE}7EGp>xmydROY9 z*wcR>3$hBEdGOA$PYf^`ZQgseQj}d>5&O|946m7a;Vw511#QX!S|FTmw$AL7t zts)XnOsY0BBD){DoJ8Or;9`IG`v^yZi(RY$cYukLz61pJ-+X87eTDUdg?r^=BBbt) znMd_@ed#-clj^8Yi$l0{cFs_{50CUYv_jxP)4}_M4p8_1^*^i4b+KLV<;H_prP{n} zdOrue^n7lrY&0Vf^&3;ece843rSqi)yF=t=#UtR})k9$#<2uS0QNmkV+zXNC@R;L82_ zJ1O2%eAd^uR&NT8yKarCwcGQrg*lYY#I{P;DA^bEb@)`DLnS;{y89B`2V4eZgmY4S z^)D;NZM+fQ8cnyhi*tq_`^qi9JF9iaF4hUHm~((9598Zf$eMbbj`Z18J^rjWyrKQi zo#pc_hVC3_${q8EdYwbq!`%dXfRZSs9H1it9_J(;}^|Jw==QN!cLS3!S z#UeHXJRBdT>f%3fm3FI(*`4j?3<%?WA*{IHIpJQba3mW3