R/genSampDs.R

In SysBioChalmers/DSAVE-R: DSAVE: Tools to Investigate Variation and Purity of Subpopulations in Single-Cell RNA-Seq Data

#' genSampDs
#'
#' Help function for DSAVEGetGeneVariation. Generates a SNO dataset for the gene-wise metric,
#' where the counts per gene is constant instead of counts per cell.
#'
#' genSampDs
#'
#' @param data numeric matrix, the input dataset (cell population)
#' @param countsPerGene counts per gene
#' @param edges edges
#' @importFrom graphics hist
#' @author Juan Inda, <inda@@chalmers.se>
#' @return list(logCV, variances)
#'
genSampDs <- function(data, countsPerGene, edges){
  new_data <- data
  new_data[,] <- 0
  numGenes <- dim(data)[1]
  for(i in 1:numGenes){
    if(countsPerGene[i] > 0){
      r <- runif(countsPerGene[i])
      a <- hist(r,edges, plot = F)$counts
      new_data[i,] <- a
    }
  }
  return(new_data)
}