plot4.R

library("data.table")

## Change time locality to prevent label  in Russian language
Sys.setlocale("LC_TIME", "C")

## 
## Code for reading the data
##

firstRow <- rowNum <- NULL

readData <- function(sourceFilePath) {
  
  ## calculates index of rows for read
  if (is.null(firstRow)) {
    # takes only 2007-02-01 and 2007-02-02 measurements
    allOccurrencesIndexes <- grep("^(1|2){1}/2/2007",readLines(sourceFilePath))
    firstRow <<- allOccurrencesIndexes[1]
    rowNum <<- length(allOccurrencesIndexes)
  }
  
  ## reads only particular rows
  householdEnergyUsage <- read.table(sourceFilePath, header = F, sep = ';', na.strings = "?", skip = firstRow-1,  
                                     nrows=rowNum, check.names = F, stringsAsFactors = F, quote='\"')
  
  ## setups meaningful column names
  colNames <- c("Date", "Time", "Global_active_power", "Global_reactive_power", "Voltage", "Global_intensity",
                "Sub_metering_1","Sub_metering_2", "Sub_metering_3")
  setnames(householdEnergyUsage, names(householdEnergyUsage), colNames)
  
  householdEnergyUsage$Date <- as.Date(householdEnergyUsage$Date, format = "%d/%m/%Y")
  dateTime <- paste(householdEnergyUsage$Date, householdEnergyUsage$Time)
  householdEnergyUsage$DateTime <- as.POSIXct(dateTime)
  
  householdEnergyUsage
}

## Read data
householdEnergyUsage <- readData("data/household_power_consumption.txt")

## Reconstruct the plot 4
par(mfrow = c(2,2), mar = c(4,4,2,1), oma = c(0,0,2,0))
with(householdEnergyUsage, {
  plot(Global_active_power ~ DateTime, type = "l", 
       ylab = "Global Active Power", xlab = "")
  plot(Voltage ~ DateTime, type = "l", 
       ylab = "Voltage", xlab = "datetime" )
  plot(Sub_metering_1 ~ DateTime, type = "l", 
       ylab = "Energy sub metering", xlab = "")
  lines(Sub_metering_2 ~ DateTime, col = 'Red')
  lines(Sub_metering_3 ~ DateTime, col = 'Blue')
  legend("topright", col = c("black", "red", "blue"),lty = 1, lwd = 2, bty = "n",
         legend = c("Sub_metering_1", "Sub_metering_2", "Sub_metering_3"))
  plot(Global_reactive_power ~ DateTime, type = "l", 
       xlab = "datetime")
})
## Save plot to a PNG file with a width of 480 pixels and a height of 480 pixels.
dev.copy(png, file = "plot4.png", height = 480, width = 480)
dev.off()