R/sim_H.R

In ayurovsky/rnaGinesis: Simulate gene expression data

#' Take in direchelet parameters
#' Return a mixture matrix
#'
#' @export
#' @import gtools

sim_H = function(S = 10,
                 seed = 1234, 
                 d.params = c("Tumor"   = .3,
                              "Stromal" = .5,
                              "Immune"  = .1,
                              "Normal"  = .1))
{
   if (sum(d.params)!=1)
  {
    stop("proportions do not add up to 1")
  }
  set.seed(seed)
  H = t(rdirichlet(S , d.params))
  rownames(H) = names(d.params)
  colnames(H) = paste0("sample_", 1:S)
  return(H)
}