plot4.R

#reads data from household_power_consumption.txt file between dates 2007-02-01 and 2007-02-02
#converts date and time fields to date and time objects
#returns the dataframe with the column names col_names defined in the function
read_data <- function(){
  col_names <- c("Date", "Time", "Global_active_power", "Global_reactive_power", "Voltage", "Global_intenstiy", "Sub_metering_1", "Sub_metering_2", "Sub_metering_3")
  
  data <- read.table("household_power_consumption.txt", header = FALSE, sep = ";", skip = 66638, nrows=2880, na.strings=c("?"), col.names = col_names)
  
  dt <- paste(data$Date, data$Time,  sep = " ")
  t <- lapply(dt, function(d) strptime(d, "%d/%m/%Y %T"))
  data$Time <- t
  dates <- lapply(data$Date, function(x){as.Date(x, format = "%d/%m/%Y")})
  data$Date <- dates
  return(data)
}

d <- read_data()



plot.new()
png(filename="plot4.png", width=480, height=480)

par(mfcol = c(2,2))
days <- c("Thu", "Fri", "Sat")


with(d, {
  plot(d$Global_active_power, type = "l", ylab = "Global Active Power (kilowatts)", xlab="", xaxt="n")
  axis(side=1, labels= days, at = c(0, 1500, 2900))
  
  plot(d$Sub_metering_1, type = "l", ylab = "Energy sub metering", xlab="", xaxt="n")
  lines(d$Sub_metering_2, col = "red")
  lines(d$Sub_metering_3, col = "blue")
  axis(side=1, labels= c("Thu", "Fri", "Sat"), at = c(0, 1500, 2900))
  legend("topright", bty = "n", legend = names(c(d[-1 : -6])), col = c("black", "red", "blue"), lty = 1, seg.len = 2)

  plot(d$Voltage, type = "l", xlab = "datetime", ylab = "Voltage", xaxt = "n")
  axis(side=1, labels= days, at = c(0, 1500, 2900))
  
  
  plot(d$Global_reactive_power, type = "l", xaxt = "n", xlab = "datetime", ylab = "Global_reactive_power")
  axis(side=1, labels= days, at = c(0, 1500, 2900))
})
dev.off()