我一直在努力编写一个函数来计算相关性,这些函数的条件是完整个案的数量大于一定数量(否则返回0表示相关性)。
问题是如何添加阈值。请参阅下面的公式,该公式仅计算每个文件的相关性并打印它们。
我补充说 如果(完成> treshhold)但它似乎没有做任何事情。
提前谢谢! 苏珊
corr<-function(directory,treshhold=0){
files_list <- list.files(directory, full.names=TRUE)
mydata <- data.frame()
treshhold<-vector(,1)
cr<-vector(,332)
for (i in 1:332){
mydata <- read.csv(files_list[i])
complete<-sum(complete.cases(mydata))
if (complete > treshhold) {
x<- cor(mydata$sulfate,mydata$nitrate,use="complete.obs")
cr<-c(x)
print(cr)
} else{
x<-0
}
}
}