ExData_Plotting1/plot3.R at master · evgaster/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
rawDataFileName <- "../household_power_consumption.txt"

# Read the data from the file.
# For now keep the date and time as a  character strings.
d <- read.table(rawDataFileName,
                header = TRUE,
                sep = ";",
                na.strings = "?",
                colClasses = c("character",
                               "character",
                               "numeric",
                               "numeric",
                               "numeric",
                               "numeric",
                               "numeric",
                               "numeric",
                               "numeric"
                ),
                comment.char = ""
)

# Keep only the needed dates.
# The dates appear to be in the D/M/YYYY format,
# without leading zeros.
d <- subset(d, Date %in% c("1/2/2007", "2/2/2007"))

# Convert the Date and Time character strings to a timestamp.
d$Timestamp <- as.POSIXct(paste(d$Date, d$Time, sep = " "),
                          format = "%d/%m/%Y %H:%M:%S"
)
# Drop the Date and Time character strings.
d <- d[, ! names(d) %in% c("Date", "Time")]

# Draw the Sub_metering_[1-3] over time.
png("plot3.png")
plot(d$Timestamp, d$Sub_metering_1,
     type = "l",
     col = "black",
     xlab = "",
     ylab = "Energy sub metering"
)
lines(d$Timestamp, d$Sub_metering_2,
      type = "l",
      col = "red"
)
lines(d$Timestamp, d$Sub_metering_3,
      type = "l",
      col = "blue"
)
legend("topright",
       legend = c("Sub_metering_1", "Sub_metering_2", "Sub_metering_3"),
       col = c("black", "red", "blue"),
       lty = 1
)
dev.off()