C4W1 Course Project 1 files

Author: PythonWillRule

getPlotData.R

@@ -0,0 +1,25 @@
+# By Richard Wheeler
+# Script for C4W1 project
+# getPlotData.R
+
+getPlotData <- function() {
+    inputfile <- file("household_power_consumption.txt", open="rt")
+    outputfile <- file("output.txt", open = "wt")
+    colnames <- readLines(inputfile,n=1, warn = FALSE)
+    colnames <- unlist(strsplit(colnames,split = ";"))
+    colnames <- c(paste0(colnames[1], colnames[2],collaspe=""), colnames[3:9])
+    while(length(data <- readLines(inputfile,n=1000, warn = FALSE)) > 0) {
+        for (i in 1:length(data)) {
+            if(any(sapply(as.Date(strsplit(data[i],";")[[1]][1],"%d/%m/%Y"),grepl,"2007-02-01 2007-02-02"))) {
+                dline <- paste(format(as.Date(strsplit(data[i],";")[[1]][1],"%d/%m/%Y"),"%Y-%m-%d"),paste(strsplit(data[i],";")[[1]][2:9], collapse = ';'),sep = " ")
+                writeLines(dline,outputfile)
+            }# Write only data at the sepecified dates
+        }
+        if(as.Date(strsplit(data[i],";")[[1]][1],"%d/%m/%Y") >= "2007-02-03") break
+    }#while LOOP
+    close(inputfile)
+    close(outputfile)
+    
+    data <- read.table("output.txt",sep=';',header = FALSE, col.names = colnames, blank.lines.skip = TRUE, stringsAsFactors = FALSE)
+    return(data)
+}#getPlotData()

plot1.R

@@ -0,0 +1,15 @@
+# By Richard Wheeler
+# Script for C4W1 project
+# plot1.R
+
+plot1 <- function() {
+    if(!exists("getPlotData", mode="function")) source("getPlotData.R")
+    pdata <- getPlotData()
+    hist(pdata$Global_active_power, breaks = 36, xaxt='n', col = "red",
+         main = "Global Active Power",
+         xlab = "Global Active Power (kilowatts)",
+         ylim = c(0,1200))
+    axis(side = 1, at=seq(0,6,2), labels=seq(0,6,2))
+    dev.copy(png,file="plot1.png")
+    dev.off()
+}#plot1()

plot1.png

@@ -0,0 +1,19 @@
+# By Richard Wheeler
+# Script for C4W1 project
+# plot2.R
+
+if(!exists("getPlotData", mode="function")) source("getPlotData.R")
+
+plot2 <- function(dispOnly = FALSE, pdata = NULL) {
+    if(dispOnly == FALSE)  pdata <- getPlotData()
+
+    plot.ts(pdata$Global_active_power, axes=FALSE, ann=FALSE)
+    box()
+    title(ylab="Global Active Power (Kilowatts)")
+    axis(2, at=c(0,2,4,6))
+    axis(1,at=c(0,1440,2880),labels=c("Thu","Fri","Sat"))
+    if(dispOnly == FALSE) {
+        dev.copy(png,file="plot2.png")
+        dev.off()
+    }
+}#plot2()

plot2.R

@@ -0,0 +1,29 @@
+# By Richard Wheeler
+# Script for C4W1 project
+# plot3.R
+
+if(!exists("getPlotData", mode="function")) source("getPlotData.R")
+
+plot3 <- function(dispOnly = FALSE, pdata = NULL) {
+    boxtype <- "o"
+    textwidth <- 600
+    if(dispOnly == FALSE) {  
+        pdata <- getPlotData()
+    } else {
+        boxtype <- "n"
+        textwidth <- 1050
+    }
+    
+    plot(pdata$Sub_metering_1, type="n", xaxt="n", ylab="Energy sub metering")
+    lines(pdata$Sub_metering_1,type="l")
+    lines(pdata$Sub_metering_2,type="l",col="red")
+    lines(pdata$Sub_metering_3,type="l",col="blue")
+    axis(1,at=c(0,1440,2880),labels=c("Thu","Fri","Sat"))
+    legend("topright",c("Sub_metering_1","Sub_metering_2","Sub_metering_3"),col=c("black","red","blue"),lty=1,
+           text.width = textwidth, cex=.75, y.intersp = .75, bty = boxtype)
+
+    if(dispOnly == FALSE) {
+        dev.copy(png,file="plot3.png")
+        dev.off()
+    }
+}#plot3()

plot2.png

@@ -0,0 +1,27 @@
+# By Richard Wheeler
+# Script for C4W1 project
+# plot4.R
+
+if(!exists("getPlotData", mode ="function")) source("getPlotData.R")
+if(!exists("plot2", mode = "function")) source("plot2.R")
+if(!exists("plot3", mode = "function")) source("plot3.R")
+
+plot4 <- function() {
+    par(mfrow=c(2,2))
+    pdata <- getPlotData()
+    plot2(TRUE, pdata)
+
+    plot(pdata$Voltage, type="n", xaxt="n", ylab="Voltage", xlab = "datetime")
+    lines(pdata$Voltage,type="l")
+    axis(1,at=c(0,1440,2880),labels=c("Thu","Fri","Sat"))
+    #text(1, labels = "datetime", cex = 1)
+    
+    plot3(TRUE, pdata)
+    
+    plot(pdata$Global_reactive_power, type="n", xaxt="n", ylab="Global Reactive Power", xlab = "datetime")
+    lines(pdata$Global_reactive_power,type="l")
+    axis(1,at=c(0,1440,2880),labels=c("Thu","Fri","Sat"))
+
+    dev.copy(png,file="plot4.png")
+    dev.off()
+}#plot4()