R/weight_by_empirical_counts.R

In anthony-aylward/chenimbalance: Allelic imbalance mapping procedure from Chen et al. A uniform survey of allele-specific binding and expression over 1000-Genomes-Project individuals (2016)

#===============================================================================
# weight_by_empirical_counts.R
#===============================================================================

#' @title Weight by empirical counts
#'
#' @param total integer vector containing coverage per variant
#' @return A one-column matrix of weights for each coverage level
#' @export
weight_by_empirical_counts <- function(total) {
  t <- as.data.frame(table(total), stringsAsFactors = FALSE)
  w <- matrix(0, max(total), 1)
  for (j in 1:nrow(t)) {
    w[as.integer(t[j, 1]), 1] <- t[j, 2]
  }
  w
}