Computing for Data Analysis Programming assignment 2 (part 3)

corr.R

corr <- function(directory, threshold = 0) {
  source("getmonitor.R")
  source("complete.R")
  all.cases <- complete("specdata", id = 1:332)
  needed.cases <- subset(all.cases, all.cases$nobs > threshold)
  result <- numeric(0)
  id <- needed.cases$id
  for (i in id) {
    data <- getmonitor(i, "specdata")
    crl <- cor(data$sulfate, data$nitrate, use = "pairwise.complete.obs")
    result <- c(result, crl)
  }
  return (result)
}