Skip to contents

Estimate Marginal distribution and Gaussian Copula for Gene Expression Matrix

Usage

Est_ModelPara(
  expr,
  anno,
  sim_method = c("ind", "copula"),
  region = NULL,
  ncores = 1,
  subsetN = NULL
)

Arguments

expr

Expression levels of input data

anno

Cell type annotation of input data

sim_method

c('ind', 'copula'). Use 'ind' if later plans to simulate independent genes (fast). Use 'Coplua" if later plans to simulated gene-gene correlations.

region

Region annotation of the input data. NULL if input has no region information.

ncores

No of cores for parallel computing.

subsetN

If the number of cells is large, one can use it to estimate expression distributions from subset of cells (e.g. subsetN=2500) per cell type.