Identify genes exhibiting evidence for splice variants (alternative exon usage/transcript isoforms) from exon-level count data using negative binomial generalized linear models.

1 2 |

`y` |
either a matrix of exon-level counts or a |

`geneID` |
vector of length equal to the number of rows of |

`dispersion` |
scalar (in future a vector will also be allowed) supplying the negative binomial dispersion parameter to be used in the negative binomial generalized linear model. |

`group` |
factor supplying the experimental group/condition to which each sample (column of |

`estimate.genewise.disp` |
logical, should genewise dispersions (as opposed to a common dispersion value) be computed if the |

`trace` |
logical, whether or not verbose comments should be printed as function is run. Default is |

This function can be used to identify genes showing evidence of splice variation (i.e. alternative splicing, alternative exon usage, transcript isoforms). A negative binomial generalized linear model is used to assess evidence, for each gene, given the counts for the exons for each gene, by fitting a model with an interaction between exon and experimental group and comparing this model (using a likelihood ratio test) to a null model which does not contain the interaction. Genes that show significant evidence for an interaction between exon and experimental group by definition show evidence for splice variation, as this indicates that the observed differences between the exon counts between the different experimental groups cannot be explained by consistent differential expression of the gene across all exons. The function `topTags`

can be used to display the results of `spliceVariants`

with genes ranked by evidence for splice variation.

`spliceVariants`

returns a `DGEExact`

object, which contains a table of results for the test of differential splicing between experimental groups (alternative exon usage), a data frame containing the gene identifiers for which results were obtained and the dispersion estimate(s) used in the statistical models and testing.

Davis McCarthy, Gordon Smyth

`estimateExonGenewiseDisp`

for more information about estimating genewise dispersion values from exon-level counts. `DGEList`

for more information about the `DGEList`

class. `topTags`

for more information on displaying ranked results from `spliceVariants`

. `estimateCommonDisp`

and related functions for estimating the dispersion parameter for the negative binomial model.

1 2 3 4 5 6 |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.