R/get_ts.R

In somspace: Spatial Analysis with Self-Organizing Maps

globalVariables(c(".", "region"))

get_ts <- function(x, nregions){
  ids <- data.table(id = x$regions[[1]], region = x$regions[[9 + nregions]])
  out <- merge(x$input_dt, ids, by = "id")  
  out <- unique(out[, .(variable = mean(variable)), .(time, region)])
  return(out)
}

cor_mat <- function(x) {
  mat <- as.matrix(dcast(x, time~region, value.var = "variable")[, -1])
  out <- cor(mat)
  return(out)
}

cor_regs <- function(x, n, ...){
  reg_ts <- get_ts(x, n)
  mat <- cor_mat(reg_ts)
  return(mat)
}