Home > Research > Publications & Outputs > Improved variational Bayes inference for transc...
View graph of relations

Improved variational Bayes inference for transcript expression estimation

Research output: Contribution to journalJournal article

<mark>Journal publication date</mark>04/2014
<mark>Journal</mark>Statistical Applications in Genetics and Molecular Biology
Issue number2
Number of pages14
Pages (from-to)203-216
Publication StatusPublished
Early online date10/01/14
<mark>Original language</mark>English


RNA-seq studies allow for the quantification of transcript expression by aligning millions of short reads to a reference genome. However, transcripts share much of their sequence, so that many reads map to more than one place and their origin remains uncertain. This problem can be dealt using mixtures of distributions and transcript expression reduces to estimating the weights of the mixture. In this paper, variational Bayesian (VB) techniques are used in order to approximate the posterior distribution of transcript expression. VB has previously been shown to be more computationally efficient for this problem than Markov chain Monte Carlo. VB methodology can precisely estimate the posterior means, but leads to variance underestimation. For this reason, a novel approach is introduced which integrates the latent allocation variables out of the VB approximation. It is shown that this modification leads to a better marginal likelihood bound and improved estimate of the posterior variance. A set of simulation studies and application to real RNA-seq datasets highlight the improved performance of the proposed method.