R：警告：當我嘗試使用我的功能時得到NA

Question

我是編程的新手，我試圖在R中編寫一個函數來計算指定的監視器列表中的污染物（硝酸鹽或硫酸鹽）的平均值每個監視器在文件夾「specdata」中都有自己的.csv文件）。我已經構建了以下功能：

pollutantmean <- function(directory="specdata", pollutant="sulfate", id=1:332) 
{ 
    files_f<-list.files(directory,full.names=TRUE) 
    d <- data.frame() 

    for(i in 1:332){ 
     d <- rbind(d,read.csv(files_f[i])) 
    } 
    if(pollutant=="sulfate"){ 
     mean(d$Sulfate[which(d$ID==id)], na.rm=TRUE) 
    } 
    else{ 
     mean(d$Nitrate[which(d$ID==id)], na.rm=TRUE) 
    } 
} 

然後我試圖與測試功能： pollutantmean(directory="specdata",pollutant="sulfate", id=1:10)

然後我得到以下錯誤：

[1] NA Warning messages: 
1: In d$ID == id : 
longer object length is not a multiple of shorter object length 
2: In mean.default(d$Sulfate[which(d$ID == id)], na.rm = TRUE) : 
argument is not numeric or logical: returning NA 

這是什麼意思？我已經多次查看了代碼，但無法確定問題所在。

謝謝。

Answer 1

在這裏，我認爲我已經在評論中實施了建議，縮短了代碼，甚至推廣了這個功能，以防您想調查其他污染物（只要確保將它們拼寫成與csv相同，包括大寫）：

pollutantmean <- function(directory="specdata", 
           pollutant="Sulfate", 
           id=1:332){ 
     files_f <- list.files(directory,full.names=TRUE) 
     d <- do.call(rbind, lapply(files_f, read.csv, stringsAsFactors=FALSE)) 
     mean(d[[pollutant]][which(d$ID %in% id)], na.rm=TRUE) 
    } 

希望這樣的作品，好運與污染物監測