optimize_variance: Optimizes variance

In mohuangx/SAVER: Single-Cell RNA-Seq Gene Expression Recovery

Optimizes variance

Description

Finds the prior parameter that maximizes the marginal likelihood given the prediction.

Usage

calc.a(y, mu, sf)

calc.b(y, mu, sf)

calc.k(y, mu, sf)

Arguments

y	A vector of observed gene counts.
mu	A vector of predictions from expr.predict.
sf	Vector of normalized size factors.

Details

calc.a returns a prior alpha parameter assuming constant coefficient of variation. calc.b returns a prior beta parameter assuming constant Fano factor. calc.k returns a prior variance parameter assuming constant variance.

Value

A vector with the optimized parameter and the negative log-likelihood.

mohuangx/SAVER documentation built on June 16, 2022, 2:42 p.m.