ExData_Plotting1/plot4.R at master · mackenziewildman/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
52
53
54
55
56
57
58
59
60
61
62
63
64
#input all data
data <- read.table("household_power_consumption.txt", header=TRUE, sep=";")

#convert date into date format
data$Date <- as.Date(data$Date, "%d/%m/%Y")

#subset dates
plotData <- data[(data$Date == "2007-02-01" | data$Date == "2007-02-02"),]

#create date/time column
  #combine Date and Time into single column
  plotData$DateTime <- paste(plotData$Date, plotData$Time)
  #convert from character to POSIXlt
  plotData$DateTime <- strptime(plotData$DateTime, "%Y-%m-%d %H:%M:%S")

#######################################################################################
#create plot4

#launch png graphics device
png(filename="plot4.png", width=480, height=480)

#make plot
par(mfrow=c(2,2))

  #top left plot (plot2)
  plot(plotData$DateTime, as.numeric(as.character(plotData$Global_active_power)),
       type="l",
       xlab="",
       ylab="Global Active Power")

  #top right plot
  plot(plotData$DateTime, as.numeric(as.character(plotData$Voltage)),
       type="l",
       xlab="datetime",
       ylab="Voltage")

  #bottom left plot (plot3)
    #initialize plot with no data
    plot(plotData$DateTime, as.numeric(as.character(plotData$Sub_metering_1)),
         type="n",
         xlab="",
         ylab="Energy sub metering")
    #add Sub_metering_1 data
    points(plotData$DateTime, as.numeric(as.character(plotData$Sub_metering_1)), type="l")
    #add Sub_metering_2 data
    points(plotData$DateTime, as.numeric(as.character(plotData$Sub_metering_2)), type="l",
           col="red")
    #add Sub_metering_3 data
    points(plotData$DateTime, as.numeric(as.character(plotData$Sub_metering_3)), type="l",
           col="blue")
    #add legend
    legend("topright", lty=1,
           col=c("black","red","blue"),
           legend=c("Sub_metering_1", "Sub_metering_2", "Sub_metering_3"),
           bty="n")
  #bottom right plot
  plot(plotData$DateTime, as.numeric(as.character(plotData$Global_reactive_power)),
       type="l",
       xlab="datetime",
       ylab="Global_reactive_power")

#close png device
dev.off()