Exploratory_Data_Analysis/plot4.R at master · demydd/Exploratory_Data_Analysis · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
#set working directory with the source data file
setwd("D:/Dimon/Coursera/Exploratory Data Analysis/project1")
#read the source data file
x<-read.csv2("household_power_consumption.txt",sep=";")
png(filename = "plot4.png", width = 480, height = 480)


#extract data for two days
data$myTime = strptime(paste(data$Date, data$Time), "%d/%m/%Y %H:%M:%S")

Sys.setlocale("LC_TIME","en_US.UTF-8")     # Depends on the system, works for me on my OS

par(mfrow = c(2, 2))

# Plot 1
plot(data$myTime, data$Global_active_power, type = "n", ylab = "Global Active Power", xlab="")
lines(data$myTime, data$Global_active_power)

# Plot 2
plot(data$myTime, data$Voltage, type = "n", ylab = "Voltage", xlab="datetime")
lines(data$myTime, data$Voltage)

# Plot 3
plot(c(data$myTime,data$myTime,data$myTime), c(data$Sub_metering_1, data$Sub_metering_2, data$Sub_metering_3), type = "n", ylab = "Energy sub metering", xlab="")
lines(data$myTime, data$Sub_metering_1, col="black")
lines(data$myTime, data$Sub_metering_2, col="red")
lines(data$myTime, data$Sub_metering_3, col="blue")
legend("topright", c("Sub_metering_1", "Sub_metering_2", "Sub_metering_3"), lty=c(1,1), col=c("black", "red", "blue"))

# Plot 4
plot(data$myTime, data$Global_reactive_power, type = "n", ylab = "Global_reactive_power", xlab="datetime")
lines(data$myTime, data$Global_reactive_power)


dev.off()