ExData_Plotting1/plot4.R at master · ESKA17/ExData_Plotting1 · 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
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
library(tibble)

## Checking if code was run before
if (file.exists("extrData.txt")) {
        extrData <- read.table("extrData.txt", header = T, sep = ";")
} else {
        ## Downloading and unzipping file if it does not exist
        if (!file.exists("household_power_consumption.zip") &
            !file.exists("household_power_consumption.txt")) {
                myURL <- "https://d396qusza40orc.cloudfront.net/exdata%2Fdata%2Fhousehold_power_consumption.zip"
                download.file(myURL, destfile = "household_power_consumption.zip")
                unzip("household_power_consumption.zip")
        }

        ## Reading data
        data <- read.table("household_power_consumption.txt", sep = ";", header = T,
                           na.strings = "?")
        ## Subseting required dates and removing stored data
        extrData <- subset(data, data$Date == "1/2/2007" | data$Date == "2/2/2007")
        rm(list = "data")

        ## Writing new table
        write.table(extrData, file = "extrData.txt", row.names = F, sep = ";")
}

## Adding  date with time column
Date_Time <- paste(extrData$Date, extrData$Time, sep = " ") %>%
        strptime(format = "%d/%m/%Y %H:%M:%S")
extrData <- add_column(extrData, Date_Time = Date_Time)
rm(list = "Date_Time")

## writing to png file "plot1.png"
png(file="plot1.png", width=480, height=480)
hist(extrData$Global_active_power, col = "red", main = "Global Active Power",
     xlab = "Global Active Power (kilowatts)")
dev.off()


## writing to png file "plot4.png"
png(file="plot4.png", width=480, height=480)
par(mfcol = c(2, 2))
plot(extrData$Date_Time, extrData$Global_active_power, xlab = "",
     ylab = "Global Active Power (kilowatts)", type = "l")
plot(extrData$Date_Time, extrData$Sub_metering_1, xlab = "",
     ylab = "Global Active Power (kilowatts)", type = "l")
lines(extrData$Date_Time, extrData$Sub_metering_2, xlab = "",
      ylab = "Global Active Power (kilowatts)", type = "l", col = "red")
lines(extrData$Date_Time, extrData$Sub_metering_3, xlab = "",
      ylab = "Global Active Power (kilowatts)", type = "l", col = "blue")
legend("topright", c("Sub_metering_1","Sub_metering_2", "Sub_metering_3"),
       col = c("black", "red", "blue"), lty = 1, lwd  = 2, box.lty = 1, cex=0.7)
plot(extrData$Date_Time, extrData$Voltage, xlab = "datetime",
     ylab = "Voltage", type = "l")
plot(extrData$Date_Time, extrData$Global_reactive_power, xlab = "datetime",
     ylab = "Global_reactive_power", type = "l")
dev.off()