testGeneForDEU {DEXSeq} | R Documentation |
This function first fits a GLM for the null model, then a GLM for the full model for each exon counting bin. Then, p values are derived with a chi-squared test from the deviance differences between the models.
testGeneForDEU( ecs, gene, formula0=NULL, formula1=NULL )
ecs |
An ExonCountSet object. |
gene |
The ID of the gene to be tested for differential exon usage. |
formula0 |
Formula for the null model. If NULL, the default
" |
formula1 |
Formula for the full model. If NULL, the default
" |
The terms in the formulas must be columns of design(ecs)
. In addition,
in formula1
, the variable exonID
is set to the
ID of the currently tested exon counting bin, looping through all the counting bins.
The GLMs are of the negative binomial family, using the dispersions from the dispersion
column in fData(ecs)
.
A data frame with columns "deviance", "df" (degrees of freedom) and pvalues from the test.
testForDEU
data("pasillaExons", package="pasilla") pasillaExons <- estimateSizeFactors( pasillaExons ) pasillaExons <- estimateDispersions( pasillaExons ) pasillaExons <- fitDispersionFunction( pasillaExons ) testGeneForDEU(pasillaExons, "FBgn0085442")