[R矢量不打印預期輸出

Question

corr <- function(directory, threshold = 0){ 

    #get all the cases that have non-NA values 
    complete_cases <- complete(directory) 
    #get all the observations over the threshold amount 
    filter_cases <- complete_cases[complete_cases[["nobs"]] > threshold, ] 

    #The returned data frame contains two columns "ID" and "nobs" 

    #get all file names in a vector 
    all_files <- list.files(directory, full.names=TRUE) 

    correlation <- vector("numeric") 

    for(i in as.numeric(filter_cases[["ID"]])){ 
    #get all the files that are in the filter_cases 
    output <- read.csv(all_files[i]) 
    #remove all NA values from the data 
    output <- output[complete.cases(output), ] 
    #get each of the correlations and store them 
    correlation[i] <- cor(output[["nitrate"]], output[["sulfate"]]) 
    } 

    correlation 
} 

我預計將從這一點說就是這樣的：

corr("directory", 200) 

[1] -1.023 0.0456 0.8231 etc 

我得到的是：

NA NA -1.023 NA NA 
NA NA NA 0.0456 NA 
0.8231 NA NA NA NA etc 

我覺得這是一件簡單的我作爲print（cor（輸出[[「硝酸鹽」]，輸出[[「硫酸鹽」]]））缺少基本上得到我所期望的。輸出必須是一個向量，因爲我打算在其他函數中使用該函數。

Answer 1

在我看來，你的問題可能是由於你的循環索引。這導致相關矢量的一些條目被跳過並因此被設置爲NAs。如果沒有訪問你的數據的信息，很難確定，但看起來上面的行的目的是讓你只能循環訪問某些文件。如果是這種情況，由於您爲了兩個目的而使用for循環，因此使關聯索引使用明確的計數器可能是有意義的，如下所示。

cor_index = 0 
for(i in as.numeric(filter_cases[["ID"]])){ 
    #get all the files that are in the filter_cases 
    output <- read.csv(all_files[i]) 
    #remove all NA values from the data 
    output <- output[complete.cases(output), ] 
    #get each of the correlations and store them 
    cor_index = cor_index + 1 
    correlation[cor_index] <- cor(output[["nitrate"]], output[["sulfate"]]) 
}