Assignment 3

Assigment 3

library(readxl)
library(data.table)
library(leaflet)
directory = "C:/Users/Erik/Documents"

mataro = read_excel(paste(directory,"mataro.xlsx",sep = "/"))
mataro <- data.frame("Date"=mataro$Date,"PM10"=mataro$PM10) 
mataro <- subset(mataro,!is.na(PM10))
mataro <- subset(mataro,PM10!="No data")

svin = read_excel(paste(directory,"santvi.xls",sep = "/")) 
svin  <- data.frame("Date"=svin$Date,"PM10"=svin$PM10) 
svin  <- subset(svin,!is.na(PM10)) 
svin  <- subset(svin,PM10!="No data") 
svin$Date  <- as.POSIXct(svin$Date,format="%d/%m/%Y %H:%M:%S",tz="CET")
svin$Day   <- as.factor(weekdays(svin$Date))
svin$Time  <- format(svin$Date,"%T")

eixample = read_excel(paste(directory,"eixemple.xls",sep = "/")) 
eixample <- data.frame("Date"=eixample$Date,"PM10"=eixample$PM10) 
eixample <- subset(eixample,!is.na(PM10)) 
eixample <- subset(eixample,PM10!="No data") 
eixample$Date <- as.POSIXct(eixample$Date,format="%d/%m/%Y %H:%M:%S",tz="CET")
eixample$Day  <- as.factor(weekdays(eixample$Date))
eixample$Time <- format(eixample$Date,"%T")
eixample$PM10 <- as.integer(eixample$PM10)

work <- subset(eixample, eixample$Day=="lunes" | eixample$Day=="martes"| eixample$Day=="miércoles"| eixample$Day=="jueves"| eixample$Day=="viernes")
work <- data.frame("Date"=work$Date,"PM10"=work$PM10) 
dayoff <- subset(eixample, eixample$Day=="sábado" | eixample$Day=="domingo")
dayoff <- data.frame("Date"=dayoff$Date,"PM10"=dayoff$PM10)

Are the PM10 concentrations from working days close to the weekends?

cat("\neixample:","\n")

## 
## eixample:

mean(work[["PM10"]])

## [1] 53.94444

mean(dayoff[["PM10"]])

## [1] 47.69325

cat("\nmataro:","\n")

## 
## mataro:

mean(mataro$PM10)

## [1] 19.48444