Bayesian Estimation for Quantification by Real-time Polymerase Chain Reaction Under a Branching Process Model of the DNA Molecules Amplification Process

Authors: Lalam, Nadia¹; Jacob, Christine²

Source: Mathematical Population Studies, Volume 14, Number 2, April 2007 , pp. 111-129(19)

Publisher: Routledge, part of the Taylor & Francis Group

Buy & download fulltext article:

Price: $50.43 plus tax (Refund Policy)

Abstract:

The aim of Quantitative Polymerase Chain Reaction is to determine the initial amount X0 of specific nucleic acids from an observed trajectory of the amplification process, the amplification being achieved through successive replication cycles. This process depends on the efficiency {pn}n of replication of the DNA molecules, pn being the probability that a molecule will duplicate at replication cycle n. Assuming pn = p for all n, Bayesian estimators of the unknown parameter θ = (p, X0) are constructed by Markov Chain Monte Carlo methods under a Bienaymé-Galton-Watson branching model of the amplification process. The Bayesian approach takes into account some prior information on the parameter. Relying on simulated data, the proposed Bayesian estimators and their credibility sets are shown to be quite accurate.

Keywords: 2000 mathematics subject classification: 60J85, 62; bayesian inference; branching processes; Markov Chain Monte Carlo; population dynamics; quantitative Polymerase Chain Reaction

Document Type: Research article

DOI: http://dx.doi.org/10.1080/08898480701298418

Affiliations: 1: EURANDOM, The Netherlands 2: INRA, Unit of Applied Mathematics and Informatics, France

Publication date: 2007-04-01