1+ # setwd("C:/Users/hariz_000/Dropbox/R/exdata-002/")
2+ data <- read.table(" household_power_consumption.txt" ,header = T , sep = " ;" ,as.is = T ,na.strings = " ?" )
3+ data $ Date <- as.Date(data $ Date ," %d/%m/%Y" )
4+ data <- data [data $ Date == ' 2007-02-01' | data $ Date == ' 2007-02-02' ,]
5+ for (i in 3 : 9 ) data [,i ] <- as.numeric(data [,i ])
6+
7+ png(" plot4.png" )
8+ par(mfrow = c(2 ,2 ))
9+ with(data ,plot(strptime(paste(Date ,Time )," %Y-%m-%d %H:%M:%S" ),
10+ Global_active_power ,type = ' l' ,xlab = " " ,
11+ ylab = ' Global Active Power (kilowatts)' ))
12+ with(data ,plot(strptime(paste(Date ,Time )," %Y-%m-%d %H:%M:%S" ),Voltage ,
13+ type = ' l' ,xlab = ' datetime' ,ylab = ' Voltage' ))
14+ with(data ,plot(strptime(paste(Date ,Time )," %Y-%m-%d %H:%M:%S" ),Sub_metering_1 ,type = ' l' ,xlab = " " ,ylab = ' Energy sub metering' ))
15+ with(data ,lines(strptime(paste(Date ,Time )," %Y-%m-%d %H:%M:%S" ),Sub_metering_2 ,col = ' red' ))
16+ with(data ,lines(strptime(paste(Date ,Time )," %Y-%m-%d %H:%M:%S" ),Sub_metering_3 ,col = ' blue' ))
17+ legend(" topright" ,legend = c(' Sub_metering_1' ,' Sub_metering_2' ,' Sub_metering_3' ),col = c(' black' ,' red' ,' blue' ),lwd = 1 )
18+ with(data ,plot(strptime(paste(Date ,Time )," %Y-%m-%d %H:%M:%S" ),
19+ Global_reactive_power ,type = ' l' ,xlab = ' datetime' ))
20+ dev.off()
0 commit comments