forked from rdpeng/ExData_Plotting1
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathplot4.R
17 lines (17 loc) · 1.01 KB
/
plot4.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
temp <- tempfile()
download.file('http://d396qusza40orc.cloudfront.net/exdata%2Fdata%2Fhousehold_power_consumption.zip',temp)
all_data<-read.csv(unzip(temp),sep = ";",header = T,quote = "",na.strings='?')
unlink(temp)
my_data<-subset(all_data,Date=='1/2/2007' | Date=='2/2/2007')
my_data$dt<-strptime(paste(my_data$Date,my_data$Time),"%d/%m/%Y %T")
png('plot4.png',bg = "transparent")
par(mfrow = c(2,2) )
plot(my_data$dt,my_data$Global_active_power,ylab='Global Active Power',xlab='',t='l')
plot(my_data$dt,my_data$Voltage,ylab='Voltage',xlab='datetime',t='l')
plot(my_data$dt,my_data$Sub_metering_1,ylab='Energy Submetering',xlab='',t='n')
lines(my_data$dt,my_data$Sub_metering_1,col='black')
lines(my_data$dt,my_data$Sub_metering_2,col='red')
lines(my_data$dt,my_data$Sub_metering_3,col='blue')
legend('topright',col=c("black","red","blue"),legend = c("Sub_metering_1","Sub_metering_2","Sub_metering_3"),lwd=1,bty ='n' )
plot(my_data$dt,my_data$Global_reactive_power,ylab='Global_reactive_power',xlab='datetime',t='l')
dev.off()