Simulate sparse counting data with a mixture of Poisson parameters

Simulate sparse counting data with a mixture of Poisson parameters

Usage

mmutil_simulate_poisson_mixture(
  r_mu_list,
  Ncell,
  output,
  dir_alpha = 1,
  gam_alpha = 2,
  gam_beta = 2,
  rseed = 42L,
  MAX_COL_WORD = 100L,
  COL_WORD_SEP = "@"
)

Arguments

r_mu_list

a list of gene x individual matrices

Ncell

the total number of cells (may not make it if too sparse)

output

a file header string for output files

dir_alpha

a parameter for Dirichlet(alpha * 1, ..., 1)

gam_alpha

a parameter for Gamma(alpha, beta)

gam_beta

a parameter for Gamma(alpha, beta)

rseed

random seed