forked from rdpeng/ExData_Plotting1
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathplot2.R
42 lines (40 loc) · 1.29 KB
/
plot2.R
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
download.file("https://d396qusza40orc.cloudfront.net/exdata%2Fdata%2Fhousehold_power_consumption.zip",
"household_power_consumption.zip",
"curl")
unzip("household_power_consumption.zip")
df <- read.table(
"household_power_consumption.txt",
sep=";",
na.strings="?",
col.names=c("Date",
"Time",
"Global_active_power",
"Global_reactive_power",
"Voltage",
"Global_intensity",
"Sub_metering_1",
"Sub_metering_2",
"Sub_metering_3"),
colClasses=c("character",
"character",
rep("numeric",7)),
comment.char="",
nrows=3000,
# skip forward to the start of the data of interest
skip=tail(grep("31/1/2007",
readLines("household_power_consumption.txt")),1)
)
# convert date column from character class to date class
df[, "Date"] <- as.Date(df[, "Date"],"%d/%m/%Y")
# remove rows outside of relevant date range
df <- subset(df,Date < "2007-02-03")
# create a new column dt that combines the date and time elements
dt <- strptime(paste(df$Date,df$Time),"%Y-%m-%d %T")
df <- cbind(df,dt)
png("plot2.png")
plot(df$dt,
df$Global_active_power,
type="l",
ylab="Global Active Power (kilowatts)",
xlab="")
dev.off()