bootVcov1: Bootstrap a zlmfit

In RGLab/MAST: Model-based Analysis of Single Cell Transcriptomics

Bootstrap a zlmfit

Description

Sample cells with replacement to find bootstrapped distribution of coefficients

Usage

pbootVcov1(cl, zlmfit, R = 99)

bootVcov1(zlmfit, R = 99, boot_index = NULL)

Arguments

cl	a cluster object created by makeCluster
zlmfit	class ZlmFit
R	number of bootstrap replicates
boot_index	list of indices to resample. Only one of R or boot_index can be offered.

Value

array of bootstrapped coefficients

array of bootstrapped coefficients

Functions

• pbootVcov1(): parallel version of bootstrapping

Examples

data(vbetaFA)
zlmVbeta <- zlm(~ Stim.Condition, subset(vbetaFA, ncells==1)[1:5,])
#Only run 3 boot straps, which you wouldn't ever want to do in practice...
bootVcov1(zlmVbeta, R=3)

RGLab/MAST documentation built on Sept. 30, 2023, 1:08 p.m.