The Annals of Probability

Confidence intervals for nonhomogeneous branching processes and polymerase chain reactions

Didier Piau
Source: Ann. Probab. Volume 33, Number 2 (2005), 674-702.

Abstract

We extend in two directions our previous results about the sampling and the empirical measures of immortal branching Markov processes. Direct applications to molecular biology are rigorous estimates of the mutation rates of polymerase chain reactions from uniform samples of the population after the reaction. First, we consider nonhomogeneous processes, which are more adapted to real reactions. Second, recalling that the first moment estimator is analytically known only in the infinite population limit, we provide rigorous confidence intervals for this estimator that are valid for any finite population. Our bounds are explicit, nonasymptotic and valid for a wide class of nonhomogeneous branching Markov processes that we describe in detail. In the setting of polymerase chain reactions, our results imply that enlarging the size of the sample becomes useless for surprisingly small sizes. Establishing confidence intervals requires precise estimates of the second moment of random samples. The proof of these estimates is more involved than the proofs that allowed us, in a previous paper, to deal with the first moment. On the other hand, our method uses various, seemingly new, monotonicity properties of the harmonic moments of sums of exchangeable random variables.

First Page: Show Hide
Primary Subjects: 60J80, 92D20
Secondary Subjects: 60K37, 60J85
Full-text: Open access
Links and Identifiers

Permanent link to this document: http://projecteuclid.org/euclid.aop/1109868597
Digital Object Identifier: doi:10.1214/009117904000000775
Mathematical Reviews number (MathSciNet): MR2123207
Zentralblatt MATH identifier: 02164479

References

Athreya, K. B. (1994). Large deviation rates for branching processes. I. Single type case. Ann. Appl. Probab. 4 779--790.
Mathematical Reviews (MathSciNet): MR1284985
Athreya, K. B. and Vidyashankar, A. N. (1995). Large deviation rates for branching processes. II. The multitype case. Ann. Appl. Probab. 5 566--576.
Mathematical Reviews (MathSciNet): MR1336883
Higuchi, R., Dollinger, G., Walsh, P. S. and Griffith, R. (1992). Simultaneous amplification and detection of specific DNA sequences. Biotechnology 10 413--417.
Higuchi, R., Fockler, C., Dollinger, G. and Watson, R. (1993). Kinetic PCR analysis: Real-time monitoring of DNA amplification reactions. Biotechnology 11 1026--1030.
Jagers, P. and Klebaner, F. (2003). Random variations and concentration effects in PCR. J. Theoret. Biol. 224 299--304.
Mathematical Reviews (MathSciNet): MR2067238
Digital Object Identifier: doi:10.1016/S0022-5193(03)00166-8
Piau, D. (2002). Mutation-replication statistics for polymerase chain reactions. J. Comput. Biol. 9 831--847.
Piau, D. (2004a). Immortal branching Markov processes: Averaging properties and polymerase chain reaction applications. Ann. Probab. 32 337--367.
Mathematical Reviews (MathSciNet): MR2040785
Digital Object Identifier: doi:10.1214/aop/1078415838
Project Euclid: euclid.aop/1078415838
Zentralblatt MATH: 1047.60087
Piau, D. (2004b). Harmonic moments of branching processes. Preprint.
Saiki, R. K., Gelfand, D. H., Stoffel, S., Scharf, S. J., Higuchi, R., Horn, G. T., Mullis, K. B. and Erlich, H. A. (1988). Primer-directed enzymatic amplification of DNA with a thermostable DNA polymerase. Science 239 487--491.
Schnell, S. and Mendoza, C. (1997). Enzymological considerations for a theoretical description of the quantitative competitive polymerase chain reaction (QC-PCR). J. Theoret. Biol. 184 433--440.
Sun, F. (1995). The polymerase chain reaction and branching processes. J. Comput. Biol. 2 63--86.
Wang, D., Zhao, C., Cheng, R. and Sun, F. (2000). Estimation of the mutation rate during error-prone polymerase chain reaction. J. Comput. Biol. 7 143--158.
Weiss, G. and von Haeseler, A. (1995). Modeling the polymerase chain reaction. J. Comput. Biol. 2 49--61.
Weiss, G. and von Haeseler, A. (1997). A coalescent approach to the polymerase chain reaction. Nucleic Acids Research 25 3082--3087.

2012 © Institute of Mathematical Statistics

The Annals of Probability

The Annals of Probability