bootstrap: Bootstrap.

In mnem: Mixture Nested Effects Models

Description

Run bootstrap simulations on the components (phi) of an object of class mnem.

Usage

bootstrap(x, size = 1000, p = 1, logtype = 2, complete = FALSE, ...)

Arguments

x	mnem object
size	size of the booststrap simulations
p	percentage of samples (e.g. for 100 E-genes p=0.5 means sampling 50)
logtype	logarithm type of the data (e.g. 2 for log2 data or exp(1) for natural)
complete	if TRUE, complete data log likelihood is considered (for very large data sets, e.g. 1000 cells and 1000 E-genes)
...	additional parameters for hte nem function

Value

returns bootstrap support for each edge in each component (phi); list of adjacency matrices

Author(s)

Martin Pirkl

Examples

sim <- simData(Sgenes = 3, Egenes = 2, Nems = 2, mw = c(0.4,0.6))
data <- (sim$data - 0.5)/0.5
data <- data + rnorm(length(data), 0, 1)
result <- mnem(data, k = 2, starts = 1)
boot <- bootstrap(result, size = 2)

mnem documentation built on Nov. 18, 2020, 2 a.m.