R/obtCount.R

In SPsimSeq: Semi-parametric simulation tool for bulk and single-cell RNA sequencing data

#' Calculates height and mid points of a distribution
#' 
#' @param Y a vector of gene expression data for a particular gene (in log CPM)
#' @param w a numeric value between 0 and 1 or NULL refering the number of classes to be created
# for the outcome data (if NULL the algorithm in graphics::hist() function will be used)
# 
#' @return a list object contating class breaks, mid points and counts
#' @importFrom graphics hist
obtCount <- function(Y, w){
  h = hist(Y, plot = FALSE, 
           nclass = if(is.null(w)) NULL else nclass = round(w*length(Y))) 
  list(breaks = h$breaks, mids = h$mids, counts = h$counts) 
}