PRIST (version 0.925)

ebayes: Estimate hyperparameters for hierarchical variance model for continuous component

Description

ebayesControl is a named list with (optional) components 'method' (one of 'MOM' or 'MLE') and 'model' (one of 'H0' or 'H1') method MOM uses a method-of-moments estimator, while MLE using the marginal likelihood. H0 model estimates the precisions using the intercept alone in each gene, while H1 fits the full model specified by formula

Usage

ebayes(sca, ebayesControl, Formula, truncate = Inf)

Arguments

sca

SingleCellAssay

ebayesControl

list with (optional) components 'method', 'model'. See details.

Formula

a formula (using variables in cData(sca) used when model='H1'.

truncate

Genes with sample precisions exceeding this value are discarded when estimating the hyper parameters

Value

numeric of length two, giving the hyperparameters in terms of a variance (v) and prior observations (df), inside a structure, with component hess, giving the Fisher Information of the hyperparameters.