我一直在努力編寫一個函數,該函數需要一個數據文件的目錄和一個完整的案例的閾值,並計算硫酸鹽和硝酸鹽之間的相關性,以監測完全觀察到的病例數量對所有變量)大於閾值。該函數應返回滿足閾值要求的監視器的相關向量。如果沒有監視器滿足閾值要求,則該函數應該返回一個長度爲0的數字向量。產生了多個錯誤,因此我在這裏沒有提及它們。R函數產生錯誤
的代碼數據文件的位置: https://d396qusza40orc.cloudfront.net/rprog%2Fdata%2Fspecdata.zip
代碼
corr<-function(directory, threshold=0){
files.list=list.files(directory, full.names=TRUE, pattern=".csv")
comp.sum<-numeric()
num<-numeric()
for(i in 1:332){
data<-read.csv(files.list[i])
data.cor<-na.omit(data[,2:3])
comp.sum<-sum(data.cor)
if
{
comp.sum>threshold
cor.var<-cor(data.cor, use="all.obs")
}
else
{
num
}
}
cor.var
}
至少缺少''if''的條件。我認爲它意外地在括號內。 – Phann