Simulate2SCE: Simulate the data for two-group comparison; e.g., treatment...

In suke18/POWSC: Simulation, power evaluation, and sample size recommendation for single cell RNA-seq

Description

Simulate the data for two-group comparison; e.g., treatment v.s. control It simulates the DE changes in two forms corresponding two types of DE genes

Usage

Simulate2SCE(n = 100, perDE = 0.05, estParas1, estParas2)

Arguments

n	the number of total cells for two groups
perDE	percentage of DE genes
estParas1	the set of parameters corresponding to cell type I
estParas2	the set of parameters corresponding to cell type II

Value

a list of metrics recording the changes in the generated data: such as the DE gene indices for Form I and II DE genes, and simulated expression data in singlecellexperiment format.

Examples

data("es_mef_sce")
sce = es_mef_sce[, colData(es_mef_sce)$cellTypes == "fibro"]
set.seed(123)
rix = sample(1:nrow(sce), 500)
sce = sce[rix, ]
estParas = Est2Phase(sce)
simData = Simulate2SCE(n=100, estParas1 = estParas, estParas2 = estParas)

suke18/POWSC documentation built on April 2, 2021, 4:34 a.m.