plot4.R

## Exploratory Data Analysis
## Course project: Plotting4
plot4 <- function() {
    # name of the plot
    plotName <- "plot4.png"
    
    # names of the data and zipped file
    fileNameZip <- "exdata-data-household_power_consumption.zip"
    fileNameTxt <- "household_power_consumption.txt"
    
    dateFrom <- as.POSIXct("2007-02-01")
    dateTo   <- as.POSIXct("2007-02-02")
    
    # check if data file exists. (in zip or txt format)
    if (!file.exists(fileNameTxt)) {
        if (!file.exists(fileNameZip)) {
            stop(paste("No data file present! (",fileNameZip,"or",fileNameTxt,")"))
        }
        unzip(fileNameZip, fileNameTxt)
    }
    
    # read the data from the text data file
    # first read the first row to get the first date in the data file
    exdata <- read.table(fileNameTxt, 
                         sep = ";",
                         header = TRUE, 
                         na.strings = "?", 
                         nrows = 1,
                         row.names = NULL,
                         colClasses = c("character","character","numeric","numeric","numeric",
                                        "numeric","numeric","numeric","numeric"))
    
    # get the columnnames
    exdataColNames <- colnames(exdata)
    
    # convert date and time column to one POSIXct Date column
    exdata$Date <- as.POSIXct(strptime(paste(exdata$Date, exdata$Time), "%d/%m/%Y %H:%M:%S"))
    
    # calculate the number of rows to skip and to get from the data file
    dateTo <- dateTo + 24*60*60
    dataSkip <- as.integer(difftime(dateFrom, exdata$Date[1], units = "mins")) + 1
    dataRows <- as.integer(difftime(dateTo, dateFrom, units = "mins"))
    
    # read the data from dateFrom to dateTo
    exdata <- read.table(fileNameTxt, 
                         sep = ";",
                         header = FALSE, 
                         na.strings = "?", 
                         skip = dataSkip,
                         nrows = dataRows,
                         col.names = exdataColNames,
                         row.names = NULL,
                         colClasses = c("character","character","numeric","numeric","numeric",
                                        "numeric","numeric","numeric","numeric"))

    # convert date and time column to one POSIXct Date column
    exdata$Date <- as.POSIXct(strptime(paste(exdata$Date, exdata$Time), "%d/%m/%Y %H:%M:%S"))
    exdata$Time <- NULL
    
    # make plot
    png(filename = plotName, width = 480, height = 480)
    par(mfcol = c(2, 2), mar = c(4, 4, 2, 1), oma = c(0, 0, 0, 0))
    
    # left upper graphic
    plot(exdata$Date, 
         exdata$Global_active_power, 
         type = "l",
         xlab = "",
         ylab = "Global active power (kilowatts)" )
    
    # left below graphic
    plot(exdata$Date, 
         exdata$Sub_metering_1, 
         type = "l",
         xlab = "",
         ylab = "Energy sub metering" )
    lines(exdata$Date,
          exdata$Sub_metering_2,
          col = "red")
    lines(exdata$Date,
          exdata$Sub_metering_3,
          col = "blue")
    legend("topright", 
           lty = c(1,1,1), 
           col = c("black","red", "blue"), 
           bty = "n",
           legend = exdataColNames[c(7:9)])
    
    # right upper graphic
    plot(exdata$Date, 
         exdata$Voltage, 
         type = "l",
         xlab = "",
         ylab = "Voltage" )
    
    # right below graphic
    plot(exdata$Date, 
         exdata$Global_reactive_power, 
         type = "l",
         xlab = "",
         ylab = "Global reactive power" )
    
    dev.off()
    message(paste("Plotted", plotName))
}