f <- function(directory,threshold){
files_list <- list.files(directory, full.names=T)
u<-numeric()
for (i in 1:332){
dat <- read.csv(files_list[i])
m <- na.omit(dat)
c <- nrow(m)
if (c>threshold){
x<-m["nitrate"]
y<-m["sulfate"]
u[i]<-cor(x,y)
}
i<-i+1
}
u
}
这是求cor的一个函数,输入f("specdata",150),显示结果为
[1] NA -0.018957541 -0.140512544 -0.043897372 -0.068159562 -0.123506666
[7] -0.075888144 -0.159673652 -0.086841940 NA 0.763128837 NA
[13] NA NA NA NA -0.157828603 N A
[19] -0.156998919 NA -0.044898818 NA 0.117249264 0.259057178
想把里面的NA去掉,当c<threshold的时候,不出现NA,直接进行下一个计算,不知道怎么改