ExData_Plotting1/plot4.R at master · waynem81/ExData_Plotting1 · GitHub

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
#pull in data from txt file
data <- read.table("C:/Users/Wayne Mogel/Downloads/household_power_consumption.txt", header=TRUE, sep=";", stringsAsFactors=FALSE, dec=".")

#exploring data
head(data, 25)

sub <- data[data$Date %in% c("1/2/2007","2/2/2007") ,]

head(sub, 25)

convertDateTime <- strptime(paste(sub$Date, sub$Time, sep=" "), "%d/%m/%Y %H:%M:%S")
head(convertDateTime, 25)

gap <- as.numeric(sub$Global_active_power)
subMeter1 <- as.numeric(sub$Sub_metering_1)
subMeter2 <- as.numeric(sub$Sub_metering_2)
subMeter3 <- as.numeric(sub$Sub_metering_3)

grp <- as.numeric(sub$Global_reactive_power)
volt <- as.numeric(sub$Voltage)

head(gap, 25)
head(subMeter1, 25)
head(subMeter2, 25)
head(subMeter3, 25)

head(grp, 25)
head(volt, 25)

png("plot4.png", width=504, height=504)

par(mfrow = c(2, 2))

plot(convertDateTime, gap, type="l", xlab="", ylab="Global Active Power", cex=0.2)

plot(convertDateTime, volt, type="l", xlab="datetime", ylab="Voltage")

plot(convertDateTime, subMeter1, type="l", ylab="Energy Submetering", xlab="")
lines(convertDateTime, subMeter2, type="l", col="red")
lines(convertDateTime, subMeter3, type="l", col="blue")
legend("topright", c("Sub_metering_1", "Sub_metering_2", "Sub_metering_3"), lty=, lwd=2.5, col=c("black", "red", "blue"), bty="o")

plot(convertDateTime, grp, type="l", xlab="datetime", ylab="Global_reactive_power")

dev.off()