在 R 中计算相关性时的 NA 答案

Question

我应该在 R 中编写一个函数，以便在每个文件的完整案例数高于指定阈值时计算“硝酸盐”和“硫酸盐”两种污染物之间的相关性，有 332 个单独的 CSV 文件，这里是到目前为止我得到了什么，但是当我运行它而不是矢量时，我会得到 NA 答案。我将非常感谢任何帮助，谢谢

corr <- function(directory, threshold=0) {
    data <- setwd("C:\\Users.....specdata")
    files <-list.files(data, full.names=FALSE)
    nobs <- c()
    threshold <- c() 
    correlation_list <- c()
    for (i in 1:332){
            dat <- read.csv(files[i], sep=",") 
            complete_cases <- dat[complete.cases(dat),]
            nobs <- c(nobs,nrow(complete_cases))
            good <- complete_cases [(nobs> threshold),]
            #if (nobs> threshold) {
            correlation_list <- c(correlation_list, cor(good["nitrate"], good["sulfate"]))      
            #}
            return(correlation_list)
    }

Answer 1

相关性是协方差除以方差，因此如果您的一个维度的方差为零，您将得到除以零。