forked from rdpeng/ExData_Plotting1
-
Notifications
You must be signed in to change notification settings - Fork 0
/
plot3.R
47 lines (46 loc) · 1.43 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
36
37
38
39
40
41
42
43
44
45
46
47
################################################################################
#
# Exploratory Data Analysis
# Project assignment 1
# Plot 3
#
################################################################################
#
# To make month and day printing US if another locale on Linux
# Uncomment codeline
#
#Sys.setlocale("LC_TIME","en_US.UTF-8")
#
# read data, note NA is ?
#
fileName <- "household_power_consumption.txt"
dataTable <- read.table(fileName,sep=";",header=TRUE,na.strings = "?")
str(dataTable)
#
# Want to get a timestamp using both date and time
#
datetime <- paste(dataTable$Date,dataTable$Time)
str(datetime)
#
# Add it to the data table
dataTable$timestamp <- strptime(datetime, format="%d/%m/%Y %H:%M:%S")
#
# Convert Date from factor to date
#
dataTable$Date <- as.Date(dataTable$Date,"%d/%m/%Y")
str(dataTable)
#
# Reduce dataset to 2007-02-01 and 2007-02-02
#
reducedData <- subset(dataTable,(dataTable$Date == "2007-02-01" |dataTable$Date == "2007-02-02" ))
str(reducedData)
#
# Plot it
png("plot3.png",width=480,height=480,units="px")
plot(reducedData$timestamp,reducedData$Sub_metering_1, xlab = "",
,ylab = "Energy sub metering",type="l")
lines(reducedData$timestamp,reducedData$Sub_metering_2,type="l",col="red")
lines(reducedData$timestamp,reducedData$Sub_metering_3,type="l",col="blue")
legend("topright", c("Sub_metering_1","Sub_metering_2","Sub_metering_3"),
,col=c("black","blue","red"),lty=c(1,1,1))
dev.off()