forked from rdpeng/ExData_Plotting1
-
Notifications
You must be signed in to change notification settings - Fork 0
/
plot4.R
38 lines (27 loc) · 1.6 KB
/
plot4.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
library(data.table)
housePower <- fread("household_power_consumption.txt", na.strings="?")
#
# 1. Convert Date columns to date
# 2. Filter Date
# 3. Merge Date and Time, convert to DateTime format and create DateTime column
housePower$Date <- as.Date(housePower$Date , "%d/%m/%Y")
subsethousePower <- subset(housePower, Date==as.Date("2007-02-01") | Date==as.Date("2007-02-02"))
subsethousePower$DateTime <- as.POSIXct(paste(subsethousePower$Date, subsethousePower$Time), format="%Y-%m-%d %H:%M:%S")
print(head(subsethousePower))
if (dev.cur()!=0) {
dev.off() # To return back to screen graphic device if necessary
}
par(mfrow=c(2,2))
with(subsethousePower, plot(DateTime, Global_active_power, xlab="", ylab="Global Active Power", pch="."))
lines(subsethousePower$DateTime, subsethousePower$Global_active_power)
with(subsethousePower, plot(DateTime, Voltage, xlab="datetime", ylab="Voltage", pch=".", type="n"))
with(subsethousePower, lines(DateTime, Voltage))
with(subsethousePower, plot(DateTime,Sub_metering_1, type="n", xlab="", ylab="Energy sub metering"))
with(subsethousePower, lines(DateTime,Sub_metering_1, col="black", pch="."))
with(subsethousePower, lines(DateTime,Sub_metering_2, col="red", pch="."))
with(subsethousePower, lines(DateTime,Sub_metering_3, col="blue", pch="."))
legend("topright", lty=1, cex=0.8, col = c("black", "red", "blue"), legend = c("Sub_metering_1", "Sub_metering_2", "Sub_metering_3"))
with(subsethousePower, plot(DateTime, Global_reactive_power, xlab="datetime", pch=".", type="n"))
with(subsethousePower, lines(DateTime, Global_reactive_power))
dev.copy(png, file="plot4.png")
dev.off()