fitDispersionFunction {DEXSeq} | R Documentation |
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
.
fitDispersionFunction(ecs)
ecs |
An ExonCountSet object. |
An ExonCountSet object with information of the fit included, as well as fData(ecs)$dispersion
filled.
data("pasillaExons", package="pasilla") pasillaExons <- estimateSizeFactors( pasillaExons ) pasillaExons <- estimateDispersions( pasillaExons ) pasillaExons <- fitDispersionFunction( pasillaExons )