fitDispersionFunction {DEXSeq}R Documentation

Fit the mean-variance function.

Description

This function fits a parametric model of the mean-dispersion relationship to the per-gene estimates of mean \hat{μ} and dispersion \hat{α}. The parametric model is

α(μ) = \frac{α_1}{μ} + α_0,

where μ is the mean, α the dispersion and α_1 and α_0 are two parameters. After this, for each exon, the maximum between the per-gene estimate \hat{α} and the modelled value \hat{α}_1/\hat{μ} + \hat{α}_0 is stored in fData$dispersion.

Usage

fitDispersionFunction(ecs)

Arguments

ecs

An ExonCountSet object.

Value

An ExonCountSet object with information of the fit included, as well as fData(ecs)$dispersion filled.

Examples

	data("pasillaExons", package="pasilla")
	pasillaExons <- estimateSizeFactors( pasillaExons )
        pasillaExons <- estimateDispersions( pasillaExons )
        pasillaExons <- fitDispersionFunction( pasillaExons )

[Package DEXSeq version 1.0.2 Index]