ExData_Plotting1/Plot4.R at master · araujorlp/ExData_Plotting1 · GitHub

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
43
44
45
#Plot4

url <- "https://d396qusza40orc.cloudfront.net/exdata%2Fdata%2Fhousehold_power_consumption.zip"
destfile <- "/Users/raquelaraujo/Documents/Cursos/JHDS/4. Exploratory Data Analysis/Exploratory Data Analysis/data.txt"
download.file(url,destfile)

household_power_consumption <- read.csv("~/Documents/Cursos/JHDS/4. Exploratory Data Analysis/Exploratory Data Analysis/data.txt", sep=";")

#You may find it useful to convert the Date and Time variables to Date/Time classes in R using the
#strptime() and as.Date() functions.

household_power_consumption$Date <- as.Date(household_power_consumption$Date, format = "%d/%m/%Y")


#household_power_consumption$Time <- strptime(household_power_consumption$Time, format = "%H:%M:%S")

household_power_consumption$Date_Time <- as.POSIXct(paste(household_power_consumption$Date, household_power_consumption$Time), format="%Y-%m-%d %H:%M:%S")

household_power_consumption_ss <- subset(household_power_consumption, Date == "2007-02-01" | Date == "2007-02-02")

#Plot4


household_power_consumption_ss$Voltage <- as.numeric(household_power_consumption_ss$Voltage)
household_power_consumption_ss$Global_reactive_power <- as.numeric(household_power_consumption_ss$Global_reactive_power)


# set the plotting area into a 2*2 array

png(file = "Plot4.png")


par(mfrow=c(2,2))

with(household_power_consumption_ss, plot(Date_Time, Global_active_power, type = "l", ylab = "Global Active Power (kilowatts)", xlab = ""))
with(household_power_consumption_ss, plot(Date_Time, Voltage, type = "l", ylab = "Voltage", xlab = "datetime"))

with(household_power_consumption_ss, plot(Date_Time, Sub_metering_1, type = "l", ylab = "Energy sub metering", xlab = ""))
with(household_power_consumption_ss, lines(Date_Time, Sub_metering_2, col="red"))
with(household_power_consumption_ss, lines(Date_Time, 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), ncol=1)

with(household_power_consumption_ss, plot(Date_Time, Global_reactive_power, type = "l", ylab = "Global_reactive_power", xlab = "datetime"))

dev.off()