make_normalmix: Generate beta (effects) from normal mixture prior

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

Usage

make_normalmix(Ngenes, pi0, beta_args = list(betapi, betamu, betasd))

Arguments

Ngenes	Total number of genes (null + non-null in the data).
pi0	Proportion of null genes.Default = "random", which selects a random number between 0 and 1 from a uniform distribution.
pi	K-component probability vector. K = number of components in the mixture prior.
mu	Length-K mean vector of the mixture components.
sd	Length-K standard error vector of the mixture components.

Value

null Binary vector of length Ngenes. 1 = null, 0 = no null.

scalar value pi0 proportion of null

Examples

betas <- make_normalmix(100, pi = .5,
                        beta_args = args.big_normal())

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