forked from rdpeng/ExData_Plotting1
-
Notifications
You must be signed in to change notification settings - Fork 0
/
plot3.R
37 lines (25 loc) · 1.11 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
25
26
27
28
29
30
31
32
33
34
35
#Plot3
#Prepare data for plotting.
#First, work out what lines to read
line1 <- strptime("2006-12-16 17:24:00", "%Y-%m-%d %H:%M:%S")
dataStart <- strptime("2007-02-01 00:01:00", "%Y-%m-%d %H:%M:%S")
dataStartLine <- as.numeric(dataStart - line1) * 24 * 60
numLines <- 48 * 60 # total number of rows to read in
#Read in relevant lines of data - assuming working directory contains household_power_consumption.txt
houseData <- read.csv("household_power_consumption.txt", header=FALSE, sep=";", skip=dataStartLine,nrows=numLines)
head(houseData)
plot.new()
png(filename = "plot3.png")
par(mar=c(4,4,2,2), col="black", cex.axis=0.8,cex.lab=0.8)
plot(houseData$V7, type="n", xaxt = "n", xlab="", ylab="Energy sub metering")
#Column 7 is submetering 1
lines(houseData$V7)
#Column 8 is submetering 2
lines(houseData$V8, col="red")
#Column 9 is submetering 2
lines(houseData$V9, col="blue")
#Legend
legend("topright", lty=1 , col = c("black", "red", "blue"), cex=0.7, legend = c("Sub_metering_1", "Sub_metering_1", "Sub_metering_3"))
#Add x axis label
axis(1, at=c(1,numLines/2,numLines), labels=c("Thu", "Fri", "Sat"))
dev.off()