non_null_sim: Simulate count matrix

In jhsiao999/ashbun: Methods and evaluation tools for single cell RNA-seq count data

Usage

non_null_sim(counts, condition, pi0, beta_args = list(betapi, betamu, betasd))

Arguments

counts	G*N null count matrix, rows are genes and columns are samples
args	List of arguments, including beta_args: parameters for normal mixture distributions. betapi Probability vector for the k components in the normal mixture. betamu Mean vector for k mixture components. betasd Standard deviation vector for k mixture components. pi0: null proportion. If pi0=="random" then pi0 will be randomly selected from U(0,1)

Examples

library(singleCellRNASeqMouseZeiselBrain)
eset <- get(data(MouseZeiselBrain))
counts <- exprs(eset)

counts_null <- makeSimCount2groups(counts,
                                   Ngenes = 1000,
                                   Nsamples = 20,
                                   sample_method = "all_genes")
counts_sim <- non_null_sim(counts_null$counts,
                           counts_null$condition,
                           pi0 = .5,
                           beta_args = args.big_normal())

jhsiao999/ashbun documentation built on May 8, 2019, 11:17 p.m.