plot3.R

# setwd("~/Documents/github_repositories/exploratory_data_analysis")
library(lubridate)
data <- read.table('household_power_consumption.txt', header = TRUE, na.strings = '?', sep = ';')

## clean data for downstream analysis
names(data) <- tolower(names(data))
data$all_time <- paste(data$date, data$time, sep = ' ')
data$time_2 <- dmy_hms(data$all_time)
data$date <- dmy(data$date)
new.data <- data[data$date == ymd(20070201) | data$date == ymd(20070202),]

## look for NA values in new.data
sapply(new.data, FUN = function(x) {sum(is.na(x))}) ## there are 0 NA values in new.data

## create plot
png('plot3.png', height = 480, width = 480)
plot(new.data$time_2, new.data$sub_metering_1, type = 'l', xlab = '', ylab = 'Energy sub metering')
lines(new.data$time_2, new.data$sub_metering_2, col = 'red')
lines(new.data$time_2, new.data$sub_metering_3, col = 'blue')
legend('topright', c('sub_metering_1','sub_metering_2','sub_metering_3'), lty = 1, col = c('black', 'red', 'blue'))

dev.off()