plot3.R

###
### This function plots 3 line charts corresponding to the energy measured by the 3 sub meters
### between 1/2/2007 and 2/2/2007.
### The data set is downloaded and unzipped if this was not done before. 
###
plot3 <- function(){
    
    # define some constants 
    fileURL <- "https://d396qusza40orc.cloudfront.net/exdata%2Fdata%2Fhousehold_power_consumption.zip"
    zipName <- "exdata_data_household_power_consumption.zip"
    fileName <- "household_power_consumption.txt"
    dateFormat <- "%d/%m/%Y"
    timeFormat <- "%H:%M:%S"
    startDate <- "1/2/2007"
    stopDate <- "2/2/2007"
    plotFile <- "plot3.png"
    Sys.setlocale("LC_TIME", "English")
    
    # make sure the data set is available on disk
    if(!file.exists(fileName)){
        if (!file.exists(zipName)){
            download.file(fileURL, zipName, mode="wb")
        }
        unzip(zipName)
    }
    
    # create new types in order to parse the data directly during loading
    setAs("character","hpcNumeric", function(from) suppressWarnings(as.numeric(from)) ) # numeric fails for "?" values
    setClass('hpcNumeric')
    
    # loading and general formatting the data set
    colClasses <- c(
        'character', # col 1 is the date
        'character', # col 2 is the time
        'hpcNumeric', # col 3 is the Global_active_power
        'hpcNumeric', # col 4 is the Global_reactive_power
        'hpcNumeric', # col 5 is the Voltage
        'hpcNumeric', # col 6 is the Global_intensity
        'hpcNumeric', # col 7 is the Sub_metering_1
        'hpcNumeric', # col 8 is the Sub_metering_2
        'hpcNumeric'  # col 9 is the Sub_metering_3
    );    
    raw <- read.table(fileName, sep=";", colClasses=colClasses, header=TRUE)
    data <- raw[raw$Date==startDate | raw$Date==stopDate,]
    # we don't have to check for NA's because there are none in the range we're interested in, 
    # but I keep it for the sake of completeness
    data <- data[!is.na(data[,3]), ] # no need to check other columns. Columns 3 to 9 have a ? together.
    
    # plot the chart
    # date+time (col 1 & 2) vs GAP (col 3)
    merged <- strptime(paste(data$Date, data$Time), paste(dateFormat, timeFormat))
    png(plotFile, width = 480, height = 480)
    plot(merged, data$Sub_metering_1, type="l", xlab="", ylab="Energy sub metering")
    lines(merged, data$Sub_metering_2, col="red")
    lines(merged, data$Sub_metering_3, col="blue")
    legend("topright", 
           legend=c("Sub_metering_1", "Sub_metering_2", "Sub_metering_3"), 
           col=c("black", "red", "blue"), 
           lty=c(1,1,1),
           lwd=c(1,1,1)
           );
    dev.off()
    
    # return the data set as used for plotting
    data 
}