DiscoverMathematical, Statistical and Computational Aspects of the New Science of MetagenomicsApplication of Bayesian model averaging and population Monte Carlo to inference from metagenomic mixture
Application of Bayesian model averaging and population Monte Carlo to inference from metagenomic mixture

Application of Bayesian model averaging and population Monte Carlo to inference from metagenomic mixture

Update: 2014-03-31
Share

Description

Co-author: Vincent Plagnol (University College London Genetics Institute)

For many practical applications, for example to uncover the pathogen that caused an infection after the acute phase, very deep short read sequencing can be effective provided that we can reliably assign short sequencing reads to species. This problem of assignment of reads to species is complicated by the fact that, in the absence of very large contigs, most short reads reads match to multiple species. This is essentially a mixture model, where the complete knowledge of all species present in the mixture provides information about the assignment of each read individually. However, metagenomic data analysis rarely formulates the problem in these terms because the very large number of potential species typically makes the inference intractable. Here, we propose a Bayesian model averaging strategy designed to explore the high dimensional space of species present in a metagenomic mixture. We use approximate Bayesian computation and a Monte Carlo strategy to implement the search o f the most appropriate mixture models. Owing to the computationally intensive aspects of the work, we used a population Monte Carlo Markov Chain to leverage the use of parallel computing. We find that the methodolgy is effective to provide a full Bayesian inference for samples with > 10M reads, hence providing interpretable Bayes Factors and posterior probabilities for practical problems that regularly arise in a clinical context.
Comments 
In Channel
Why mixture modelling?

Why mixture modelling?

2014-04-1001:28:00

loading
00:00
00:00
1.0x

0.5x

0.8x

1.0x

1.25x

1.5x

2.0x

3.0x

Sleep Timer

Off

End of Episode

5 Minutes

10 Minutes

15 Minutes

30 Minutes

45 Minutes

60 Minutes

120 Minutes

Application of Bayesian model averaging and population Monte Carlo to inference from metagenomic mixture

Application of Bayesian model averaging and population Monte Carlo to inference from metagenomic mixture

Vincenzo Abete

We and our partners use cookies to personalize your experience, to show you ads based on your interests, and for measurement and analytics purposes. By using our website and our services, you agree to our use of cookies as described in our Cookie Policy.