scNoiseSim: simulation to generate scRNA-seq data with varying level of...

In quon-titative-biology/bfa: A dimensionality reduction tool using gene detection pattern to mitigate noisy expression profile of scRNA-seq

Description

simulation to generate scRNA-seq data with varying level of gene detection noise versus gene count noise

Usage

scNoiseSim(zinb, celltype, disper, var_dropout = 1, var_count = 1, delta)

Arguments

zinb	a ZINB-WaVE object representing ZINB-WaVE fit to real data to get realistic simulation parameters
celltype	a factor to specify the ground-truth cell types in the original dataset that the parameter of zinb object is fit to. Since we filter out some simulated cells due to low amount of genes detected in that cell, we subset the ground truth cell types correspondingly
disper	numeric value, parameter to control the size factor r in NB(μ, r). r is varied in the set 0.5,1,5 in our simulation(as outlined in our paper)
var_dropout	numeric value, parameter to control the noise level added to a common embedding space for to generate gene detection matrix. This parameter is formulated as σ_π and in the paper is selected from the set 0.1, 0.5, 1, 2, 3
var_count	numeric value, parameter to control the noise level added to a common embedding space to generate gene count matrix. This parameter is formulated as σ_μ and and in the paper is selected from the set 0.1, 0.5, 1, 2, 3
delta	intercept to control the overall gene detection rate. and in the paper is selected from the set -2, -0.5, 1,2.5,4

Value

GeneExpr,a count matrix with rows number of genes and columns number of cells

celltype,a vector specify the corresponding celltype after QC measures.

Examples

## raw counts matrix with rows are genes and columns are cells
data("zinb_toy",package = "scBFA", envir = environment())
## a vector specify the ground truth of cell types provided by conquer database
data("celltype_toy",package = "scBFA",envir = environment())

scData = scNoiseSim(zinb = zinb_toy,
         celltype = celltype_toy,
         disper = 1,
         var_dropout =1,
         var_count = 1,
         delta = 1)

quon-titative-biology/bfa documentation built on March 18, 2020, 11:53 a.m.