Gaussian process modelling of transcription factor networks using Markov Chain Monte-Carlo
published: April 17, 2008, recorded: March 2008, views: 4460
Report a problem or upload filesIf you have found a problem with this lecture or would like to send us extra material, articles, exercises, etc., please use our ticket system to describe your request and upload the data.
Enter your e-mail into the 'Cc' field, and we will keep you updated with your request's status.
Ordinary differential equations (ODEs) can provide an useful framework for modelling the dynamics of biological networks. In this study, we focus on a small biological sub-system where a set of target genes are regulated by one transcription factor protein. The concentration of the protein and the gene specific kinetic parameters such as basal rates, decay rates and sensitivities are typically unknown. The objective of modelling is to estimate these quantities by making use of a set of observed gene expression levels. We consider a Bayesian framework for modelling the system of ODEs that is based on Gaussian processes. The Gaussian process is used as the prior for the transcription factor protein and allows us to infer the concentration of the protein in a time continuous manner. We present a Markov chain Monte Carlo algorithm for a full Bayesian statistical inference. The essential property of our MCMC algorithm is that we efficiently infer the protein concentration by applying a novel sampling algorithm for Gaussian process models. We apply our technique to linear and non-linear models.
Link this pageWould you like to put a link to this lecture on your homepage?
Go ahead! Copy the HTML snippet !