》》》》》》》》》》》》》》》》》》》》》》》》》》》》》》》》》》
pollutantmean <- function(directory, pollutant, id = 1:332) {
total = 0
observations = 0
for (i in id)
{
if (i <10) {
data <- read.csv(paste("0","0", as.character(i), ".csv", sep=""), ## for example, if 'id' =7, we get 007.csv
header = T,
na.strings=c("NA","NaN", " "))
}
else if (i>=10 & i<100) {
data <- read.csv(paste("0", as.character(i), ".csv", sep=""), ## for example, if 'id' = 17, we get 017.csv
header = T,
na.strings=c("NA","NaN", " ")
)
}
else {
data <- read.csv(paste(as.character(i), ".csv", sep=""), ## Normal
header = T,
na.strings=c("NA","NaN", " ")
)
}
data = na.omit(data)
observations = observations + nrow(data)
if (pollutant == "sulfate") {total = total + sum(data$sulfate)}
else {total = total + sum(data$nitrate)}
}
return(total/observations)