rdpeng / exdata_plotting1 Goto Github PK

hpc <- read.csv("household_power_consumption.txt", sep=";",
na.strings="?",stringsAsFactors=FALSE)
hpc2<-subset(hpc,hpc$Date==c("1/2/2007","2/2/2007"))
hpc2$Timestamp = strptime(paste(hpc2$Date, hpc2$Time),
format = "%d/%m/%Y %H:%M:%S")
Sys.setlocale("LC_TIME", "C")
with(hpc2,plot(Timestamp,Global_active_power, type ="l", ylab = "Global active power (kilowatts)"
,xlab=" ",ylim=c(0,6)))
dev.copy(png,file="plot2.png",width=480,height=480)
dev.off()

https://github.com/SurabhiShah/Project_1.git

So nice!

hpc <- read.csv("household_power_consumption.txt", sep=";",
na.strings="?",stringsAsFactors=FALSE)
hpc2<-subset(hpc,hpc$Date==c("1/2/2007","2/2/2007"))
hpc2$Timestamp = strptime(paste(hpc2$Date, hpc2$Time),
format = "%d/%m/%Y %H:%M:%S")
Sys.setlocale("LC_TIME", "C")
par(mfrow=c(2,2))
with(hpc2,plot(Timestamp,Global_active_power,type="l",ylim=c(0,6),xlab=" ",ylab="Global Active Power"))
with(hpc2,plot(Timestamp,Voltage,type="l",ylim=c(234,246),xlab="datetime",ylab="Voltage"))
with(hpc2,plot(Timestamp,Sub_metering_1,type="l",col="black", xlab=" ",ylab="Energy sub metering"))
points(hpc2$Timestamp,hpc2$Sub_metering_2,type="l",col="red")
points(hpc2$Timestamp,hpc2$Sub_metering_3,type="l",col="blue")
legend("topright",cex=0.3,lty=1:1,bty="n",col=c("black","red","blue"), legend=c("Sub_metering_1","Sub_metering_2","Sub_metering_3"))
with(hpc2,plot(Timestamp,Global_reactive_power,type="l",ylim=c(0.0,0.5),xlab="datetime",ylab="Global_reactive_power"))
dev.copy(png,file="plot4.png")
dev.off()

plot(Sub_metering_1, xlab=" ", ylab="Energy Sub Metering",
type="l", col="black")
lines(Sub_metering_2, col="red")
lines(Sub_metering_3, col="blue")
legend("topright", legend = c("Sub_metering_1","Sub_metering_2","Sub_metering_3"), col = c("black", "red", "blue"), lty = c(1,1,1))
dev.copy(png, file = "plot3.png", width=480,height=480) ## Copy plot to a PNG file
dev.off() ## close the PNG device

For the 4th plot:

The y-label of the 4th subplot (Global Reactive Power vs datetime) has not been formatted yet. The prettier y-label should be like, 'Global Reactive Power'.

Is this a bug or intention by the author?

plot2 <- plot(Global_active_power, type ="l", ylab = "Global active power (kilowatts)", xlab = "") ## Create plot on screen device

dev.copy(png, file = "plot2.png") ## Copy plot to a PNG file
dev.off() ## close the PNG device

plot1.R

data
+#plot 1
+PowerConsumption <- read.table("household_power_consumption.txt",sep=";",header=TRUE,stringsAsFactors=FALSE)
+head(PowerConsumption)
+dim(PowerConsumption)
+class(PowerConsumption)
+
+#Plot 1 - using data from the dates 02-01-2007 and 02-02-2007
+PowerConsumption$Date=strptime(PowerConsumption$Date, format="%d/%m/%Y")
+PowerConsumption$Global_active_power=as.numeric(PowerConsumption$Global_active_power)
+PowerConsumption$Global_intensity=as.numeric(PowerConsumption$Global_intensity)
+PowerConsumption$Global_reactive_power=as.numeric(PowerConsumption$Global_reactive_power)
+PowerConsumption$Sub_metering_1=as.numeric(PowerConsumption$Sub_metering_1)
+PowerConsumption$Sub_metering_2=as.numeric(PowerConsumption$Sub_metering_2)
+PowerConsumption$Sub_metering_3=as.numeric(PowerConsumption$Sub_metering_3)
+PowerConsumption$Voltage=as.numeric(PowerConsumption$Voltage)
+
+PowerConsumption=subset(PowerConsumption,PowerConsumption$Date>as.POSIXlt("2007-01-31") & PowerConsumption$Date<as.POSIXlt("2007-02-03"))
+
+png(file="plot1.png",width=480,height=480)
+hist(PowerConsumption$Global_active_power,col="red",main="Global Active Power",xlab="Global Active Power (kilowatts)")
+dev.off()

plot1 <- hist(Global_active_power, col ="red", main = paste("Global active power"), xlab = "Global active power (kilowatts)")
plot1 ## Create plot on screen device
dev.copy(png, file = "plot1.png") ## Copy plot to a PNG file
dev.off() ## close the PNG device

hpc <- read.csv("household_power_consumption.txt", sep=";",
na.strings="?",stringsAsFactors=FALSE)
hpc2<-subset(hpc,hpc$Date==c("1/2/2007","2/2/2007"))
hpc2$Timestamp = strptime(paste(hpc2$Date, hpc2$Time),
format = "%d/%m/%Y %H:%M:%S")
Sys.setlocale("LC_TIME", "C")
with(hpc2,plot(Timestamp,Sub_metering_1,type="l",col="black",
xlab=" ",ylab="Energy sub metering"))
points(hpc2$Timestamp,hpc2$Sub_metering_2,type="l",col="red")
points(hpc2$Timestamp,hpc2$Sub_metering_3,type="l",col="blue")
legend("topright",cex=0.6,lty=1:1,col=c("black","red","blue"),
legend=c("Sub_metering_1","Sub_metering_2","Sub_metering_3"))
dev.copy(png,file="plot3.png")
dev.off()

hpc <- read.csv("household_power_consumption.txt", sep=";", na.strings="?",
stringsAsFactors=FALSE)
hpc2 <- hpc[min(which(hpc$Date == "1/2/2007")):max(which(hpc$Date == "2/2/2007")),]
hist(hpc2$Global_active_power,main="Global Active Power",xlab="Global Active Power(kilowatts)",col="red",breaks=12,ylim=c(0,1200))
dev.copy(png,width=480,height=480,file="plot1.png")
dev.off()