R/match.patterns.R

In robyn-ball/PMCA: Permutation-based Maximum Covariance Analysis

#' Title
#'
#' @param score an array of distance scores
#' @param w a numeric vector denoting the optimal window widths for each component
#' @param rowterms a character vector denoting the elements of Y (rownames(Y))
#'
#' @return a list that maps rowterms of Y onto rowterms of X for each component
#' @export
#'
#' @examples see README
match.patterns <- function(score,w,rowterms) {
  n <- length(w)
  if (length(dim(score))==2) {
    g <- vector("list",n)
    g[[1]] <- rowterms[which(score[,1]<=w[1])]
    for (j in 2:n) {
      g[[j]] <- intersect(g[[j-1]],rowterms[which(score[,j]<=w[j])])
    }
  } else {
    p <- dim(score)[3]
    g <- vector("list",p)
    for (i in 1:p) {
      g[[i]] <- vector("list",n)
      g[[i]][[1]] <- rowterms[which(score[,1,i]<=w[1])]
      for (j in 2:n) {
        g[[i]][[j]] <- intersect(g[[i]][[j-1]],rowterms[which(score[,j,i]<=w[j])])
      }
    }
  }
  g
}