-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathpollutantmean.R
30 lines (30 loc) · 1.07 KB
/
pollutantmean.R
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
pollutantmean <- function(directory, pollutant, id = 1:332) {
setwd(directory)
for (i in id) {
if (i < 10) {
if (exists("pollution")) {
temp <- read.csv(paste("00",i,".csv",sep=""))
pollution <- rbind(pollution, temp)
} else {
pollution <- read.csv(paste("00",i,".csv",sep=""))
}
} else if (i >= 10 & i < 100) {
if (exists("pollution")) {
temp <- read.csv(paste("0",i,".csv",sep=""))
pollution <- rbind(pollution, temp)
} else {
pollution <- read.csv(paste("0",i,".csv",sep=""))
}
} else if (i >= 100) {
if (exists("pollution")) {
temp <- read.csv(paste(i,".csv",sep=""))
pollution <- rbind(pollution, temp)
} else {
pollution <- read.csv(paste(i,".csv",sep=""))
}
}
}
colNum = if(pollutant == "sulfate"){ 2 } else if (pollutant == "nitrate") { 3 } else NA
polMean <- mean(pollution[,colNum], na.rm = TRUE)
polMean
}