Нужно ли проходить через векторные элементы в R для вычисления корреляции? - PullRequest
0 голосов
/ 29 августа 2018

Я создал два вектора и «очистил» данные. Я собираюсь просмотреть каждый элемент в каждом векторе в виде пар, чтобы рассчитать корреляцию каждой пары и сохранить в третьем векторе.

Я могу видеть каждый элемент и печатать - однако функция cor () выдает результат NA.

Код ниже, любой совет приветствуется.

corr <- function(directory, threshold = 0) {
  files_list <- list.files(directory, full.names = TRUE)

  dat <- data.frame()
  cleandat <- data.frame()
  correlation <- c()

  for (count in 1:length(files_list)) {
    dat <- rbind(dat, read.csv(files_list[count]))
    cleandat <- dat[complete.cases(dat),]
  }

  if (nrow(cleandat) <= threshold) {
    print("Nope.")
  } else {
    sulfate_data <- cleandat$sulfate
    nitrate_data <- cleandat$nitrate

    for (element in 1:length(sulfate_data)) {
      print("Sulfate: ")
      print(sulfate_data[element])
      print("Nitrate: ")
      print(nitrate_data[element])
      print(cor(sulfate_data[element], y = nitrate_data[element]))
      correlation <- cor(sulfate_data[element], y = nitrate_data[element])

    }
  }
  print(correlation)
} #end corr()
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...