Samples from SAVER

Source: R/sample_saver.R

Samples from the posterior distribution output by SAVER.

sample.saver(x, rep = 1, efficiency.known = FALSE, seed = NULL)

Arguments

x

A saver object.
rep

Number of sampled datasets. Default is 1.
efficiency.known

Whether the efficiency is known. Default is FALSE.
seed

seed used in set.seed.

Value

A matrix of expression values sampled from SAVER posterior. If rep > 1, a list of matrices is returned

Details

The SAVER method outputs a posterior distribution, which we can sample from for downstream analysis. The posterior distribution accounts for uncertainty in the SAVER estimation procedure. If the efficiency is known, negative binomial sampling is performed; otherwise, gamma sampling is performed.

Examples

data("linnarsson_saver")

samp1 <- sample.saver(linnarsson_saver, seed = 50)