forked from rdpeng/ExData_Plotting1
-
Notifications
You must be signed in to change notification settings - Fork 0
/
plot4.R
17 lines (16 loc) · 1.12 KB
/
plot4.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
header <- read.table("./data/household_power_consumption.txt", nrows = 1, header = FALSE, sep=';', stringsAsFactors = FALSE)
dat <- read.table("./data/household_power_consumption.txt", na.strings = '?', quote = "", sep = ';', header = FALSE, skip=66637, nrows = 2880)
colnames(dat) <- unlist(header)
dat$DateTime = paste(dat$Date, dat$Time)
dat$DateTime = strptime(dat$DateTime, "%d/%m/%Y %H:%M:%S")
png(filename = "plot4.png", width=480, height=480)
par(mfrow = c(2,2))
plot(dat$DateTime, dat$Global_active_power, type="l", ylab = "Global Active Power (kilowatts)", xlab = "")
plot(dat$DateTime, dat$Voltage, type = "l", ylab = "Voltage", xlab = "datetime")
plot(dat$DateTime, dat$Sub_metering_1, type="l", col = "black", ylab = "Energy sub metering", xlab = "")
lines(dat$DateTime, dat$Sub_metering_2, type="l", col = "red")
lines(dat$DateTime, dat$Sub_metering_3, type="l", col = "blue")
labels <- c("Sub_metering_1", "Sub_metering_2", "Sub_metering_3")
legend("topright", labels, lty=1,bty="n", col = c("black","red","blue"))
plot(dat$DateTime, dat$Global_reactive_power, type = "l", ylab = "Global_reactive_power", xlab = "datetime")
dev.off()