-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathpollutantmean.r
58 lines (42 loc) · 1.03 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
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
formfilename <- function(id){
flname <- NA
counter <- 0
for (i in id)
{
if (i < 10){
fname <- paste("00",i,sep="")
}
else{
if (i <100){
fname <- paste("0",i,sep="")
}
else{
fname <- toString(i)
}
}
counter <- counter+1
flname[counter] <- paste(fname,".csv",sep="")
}
return(flname)
}
pollutantmean <- function(directory, pollutant, id){
if(missing(id))
id<-1:332
flnames <- formfilename(id)
flnames <- flnames[!is.na(flnames)]
finaldata <- NA
polluData <- NA
count <- 0
file.names <- list.files(paste(directory,"/", sep=""), full.names=FALSE)
for(i in 1:length(flnames)){
for(file in file.names){
if(flnames[i]==file) {
file.data <- read.csv(paste(paste(directory,"/", sep=""),file,sep=""))
polluData <- file.data[pollutant]
polluData <- polluData[!is.na(polluData)]
finaldata <- c(finaldata[!is.na(finaldata)], polluData)
}
}
}
return(mean(finaldata))
}