calc.post {SAVER}

Calculates SAVER posterior

Description

Given prediction and prior variance, calculates the Gamma posterior distribution parameters for a single gene.

Usage

calc.post(y, mu, sf, scale.sf)

Arguments

y

A vector of observed gene counts.

mu

A vector of prior means.

sf

Vector of normalized size factors.

scale.sf

Mean of the original size factors.

Details

Let \alpha be the shape parameter and \beta be the rate parameter of the prior Gamma distribution. Then, the posterior Gamma distribution will be

Gamma(y + \alpha, sf + \beta),

where y is the observed gene count and sf is the size factor.

Value

A list with the following components

estimate

Recovered (normalized) expression

se

Standard error of expression estimate