forked from rdpeng/ExData_Plotting1
-
Notifications
You must be signed in to change notification settings - Fork 0
/
plot3.r
24 lines (19 loc) · 1.12 KB
/
plot3.r
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
# I will use data.table instead of data.frame, they are much faster!
library(data.table)
#Fetch the data set
power_consumption <- fread("./household_power_consumption.txt",na.strings=c("?","NA","N/A","null"))
#Change class Date column
power_consumption$Date <- as.Date(power_consumption$Date,"%d/%m/%Y")
#Obtain subset of data from 2007-02-01 to 2007-02-02
power2days <- subset(power_consumption,(Date >= as.Date("2007-02-01")) & (Date < as.Date("2007-02-03")))
#Create timestamp column with the exact date and time of each measurement
power2days$Timestamp <- paste(power2days$Date,power2days$Time)
power2days$Timestamp <- as.POSIXct(power2days$Timestamp,"%Y-%m-%d %H:%M:%S")
#Plot the data and save it to plot3.png
png(file="figure/plot3.png")
with(power2days,plot(Timestamp,Sub_metering_1,type = "n",xlab="",ylab="Energy sub metering"))
with(power2days,lines(Timestamp,Sub_metering_1))
with(power2days,lines(Timestamp,Sub_metering_2,col = "red"))
with(power2days,lines(Timestamp,Sub_metering_3,col = "blue"))
legend("topright", lwd = 1,col = c("black","red", "blue"), legend = c("Sub_metering_1","Sub_metering_2","Sub_metering_3"))
dev.off()