ExData_Plotting1/plot4.R at master · rlapointe/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
46
47
48
49
50
51
# Note: for all scripts I'm assuming the data has already been downloaded,
# unzipped, and placed in a data folder locally. This is in line with the assignment.

## Read in data
setwd("~/Documents/Helpful Docs/Coursera/ExploratoryDataAnalysis")
power_data = read.csv2("./data/household_power_consumption.txt", )

# Cut down data to 2007-02-01 and 2007-02-02
power_data$Date = as.Date(power_data$Date, format = '%d/%m/%Y')
power_data_cut_down = power_data[power_data$Date >= as.Date('2007-02-01') & power_data$Date <= as.Date('2007-02-02'),]

# Make a datetime field for use in plotting
power_data_cut_down$DateTime = paste(power_data_cut_down$Date, power_data_cut_down$Time)
power_data_cut_down$DateTime = strptime(power_data_cut_down$DateTime, format = '%Y-%m-%d %H:%M:%S')

# Transform the value fields into numeric for plotting
power_data_cut_down$Sub_metering_1 = as.numeric(as.character(power_data_cut_down$Sub_metering_1))
power_data_cut_down$Sub_metering_2 = as.numeric(as.character(power_data_cut_down$Sub_metering_2))
power_data_cut_down$Sub_metering_3 = as.numeric(as.character(power_data_cut_down$Sub_metering_3))
power_data_cut_down$Voltage = as.numeric(as.character(power_data_cut_down$Voltage))
power_data_cut_down$Global_active_power = as.numeric(as.character(power_data_cut_down$Global_active_power))
power_data_cut_down$Global_reactive_power = as.numeric(as.character(power_data_cut_down$Global_reactive_power))

# Plot the data to PNG format
png(file = "plot4.png")
par(mfcol = c(2,2)) # build 2 rows and 2 cols, filling cols first

#upper left plot
plot(power_data_cut_down$DateTime, power_data_cut_down$Global_active_power,
     ylab = "Global Active Power", xlab = "",
     main = "", type="l")

#lower left plot
plot(power_data_cut_down$DateTime, power_data_cut_down$Sub_metering_1,
     ylab = "Energy sub metering", xlab = "",
     main = "", type="l")
lines(power_data_cut_down$DateTime, power_data_cut_down$Sub_metering_2, col = 'red')
lines(power_data_cut_down$DateTime, power_data_cut_down$Sub_metering_3, col = 'blue')
legend("topright", lty = c(1, 1, 1), col = c("black", "red", "blue"), bty = "n",
       legend = c("Sub_metering_1", "Sub_metering_2", "Sub_metering_3"))

#upper right plot
plot(power_data_cut_down$DateTime, power_data_cut_down$Voltage,
     ylab = "Voltage", xlab = "datetime",
     main = "", type="l")

#lower right plot
plot(power_data_cut_down$DateTime, power_data_cut_down$Global_reactive_power,
     ylab = "Global_reactive_power", xlab = "datetime",
     main = "", type="l")
dev.off()