corr <- function(directory, threshold = 0){
source("complete.R")
ord <- complete(directory, id = 1:332)
nn <- list.files(directory, full.names = TRUE)
ma <- max(ord$"nobs")
if(threshold > ma){
numeric()
} else{
orde <- subset(ord, nobs >= threshold)
tar <- as.vector(as.matrix(orde$id))
tab <- nn[tar]
res_c <- numeric(length(tar))
for (i in 1 : length(tar)){
b <- read.csv(tab[i])
bb <- na.omit(b)
bbb <- bb[c("sulfate", "nitrate")]
res_c[i] <- cor(bbb[,1], bbb[,2])
}
res <- na.omit(res_c)
}
}


雷达卡
京公网安备 11010802022788号







