plot4.R

library(dplyr)
dat <- read.csv('household_power_consumption.txt', sep=';',
                na.strings = '?')

fil_dat <- filter(dat, dat$Date == '1/2/2007' | dat$Date == '2/2/2007')
fil_dat$datetime <- strptime(paste(fil_dat$Date, fil_dat$Time), '%d/%m/%Y %T')

png(filename='plot4.png', width=480, height=480)
par(mfrow=c(2,2))
with(fil_dat,{
     plot(datetime, Global_active_power, type='n',
          xlab='',
          ylab='Global Active Power')
     lines(datetime, Global_active_power)
     plot(datetime, Voltage, type='n')
     lines(datetime, Voltage)
     plot(datetime, Sub_metering_1, type='n', 
          ylab='Energy sub metering',
          xlab='')
     lines(datetime, Sub_metering_1)
     lines(datetime, Sub_metering_2, col='red')
     lines(datetime, Sub_metering_3, col='blue')
     legend('topright', lty=1, col=c('black', 'red', 'blue'),
       legend=c('Sub_metering_1', 'Sub_metering_2', 'Sub_metering_3'))
     plot(datetime, Global_reactive_power, type='n')
     lines(datetime, Global_reactive_power)
})
dev.off()