plot4.r

######## Exploratory Data Analysis -- Assignment 1

########### CLEAR AND INITIALISE
rm(list=ls()) #Clear all data
setwd('D:/Courses/DataScience/ExploratoryDataAnalysis/Assigments/Assignment1')
# Load relevant packages
require(dplyr)

########### LOAD DATA
data <- read.table(unz("exdata-data-household_power_consumption.zip", 
                       "household_power_consumption.txt"), header=T,sep=";",
                   colClasses=c(rep('character',9)))

str(data)

data<-mutate(data,Date=as.POSIXct(paste(Date,Time,sep=' '),
                                  format='%d/%m/%Y %H:%M:%S'))%>%dplyr::select(-Time)


### select relevant dates
startDate<-as.POSIXct('2007-02-01 00:00:00')
endDate<-as.POSIXct('2007-02-03 00:00:00') ### 3rd as we want all the data right up to end of the 2nd
data<-dplyr::filter(data,Date>=startDate&Date<endDate)

#### Convert numeric data
### Create function to transform data
convertNumericWithQuestionMark<-function(x)
{
        
        ix<-which(x!='?') ## assume anything but '?' is a number
        retVec<-rep(NA,length(x)) ## initialise all as NA
        retVec[ix]<-as.numeric(x[ix]) ## and only populate valid
        retVec
}

### transform to numeric data
dataNumeric<-select(data,-Date) %>% mutate_each(funs(convertNumericWithQuestionMark)) %>%
        mutate(Date=data$Date) %>% select(c(8,1:7))
str(dataNumeric)

######### PLOTS
#### Plot4
png(file='plot4.png',width=480,height=480,type='windows')
par(mfrow=c(2,2)) # set plotting to 2x2 grid, fill row wise
with(dataNumeric,{
        ## plot 1x1
        plot(Date,Global_active_power,
                type='l',xlab='',ylab='Global Active Power')
        
        ## plot 1x2
        plot(Date,Voltage,
             type='l',xlab='datetime',ylab='Voltage')
        
        ## plot 2x1
        plot(Date,Sub_metering_1,ylim=range(c(Sub_metering_1,Sub_metering_2,Sub_metering_3))
             ,type='l',xlab='',ylab='Energy sub metering')
        lines(Date,Sub_metering_2,col='red')
        lines(Date,Sub_metering_3,col='blue')
        legend('topright',bty='n',legend=
               c('Sub_metering_1','Sub_metering_2','Sub_metering_3'),
       col=c('black','red','blue'),lty=rep(1,3))

        ## plot 2x2
        plot(Date,Global_reactive_power,
        type='l',xlab='datetime',)
})
par(mfrow=c(1,1)) # set plotting grid back to single
dev.off()