Я создал два вектора и «очистил» данные. Я собираюсь просмотреть каждый элемент в каждом векторе в виде пар, чтобы рассчитать корреляцию каждой пары и сохранить в третьем векторе.
Я могу видеть каждый элемент и печатать - однако функция 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()