<title>Introduction</title> <title>Digital PCR and Digital Array</title> <p>Digital PCR conventionally utilizes sequential limiting dilutions of target DNA, followed by amplification using the polymerase chain reaction (PCR) ##UREF##0##[1]##, ##REF##18301740##[2]##. As a result, it is possible to quantitate single DNA target molecules. We utilize the digital array, which is a novel nanofluidic biochip ##REF##18301740##[2]##, ##UREF##1##[3]## where digital PCR reactions can be performed (##FIG##0##Figure 1##) by partitioning DNA molecules, instead of diluting them. This chip utilizes integrated channels and valves that partition mixtures of sample and reagents into 765 nanolitre volume reaction chambers. DNA molecules in each mixture are randomly partitioned into the 765 chambers of each panel (the total volume of the PCR mix in each panel: 6 nl×765 = 4.59 µl). The chip is then thermocycled and imaged on Fluidigm's BioMark real-time PCR system and the positive chambers that originally contained 1 or more molecules can be counted by the digital array analysis software (##FIG##1##Figure 2##).</p> <title>Copy number variation</title> <p>Copy number variations (CNVs) are the gains or losses of genomic regions which range from 500 bases on upwards in size. Whole genome studies have revealed the presence of large numbers of CNV regions in human and a broad range of genetic diversity among the general population ##REF##15286789##[4]##, ##REF##15273396##[5]##, ##REF##17122850##[6]##. CNVs have been the focus of many recent studies because of their roles in human genetic disorders ##REF##17160897##[7]##, ##REF##17668371##[8]##, ##REF##17597781##[9]##.</p> <p>Current whole-genome scanning technologies use array-based platforms (array-CGH and high-density SNP microarrays) to study CNVs. They are high throughput but lack resolution and sensitivity. Real-time PCR is a sequence-specific technique which is easy to perform, but is limited in its discriminating power beyond a 2-fold difference ##REF##17597776##[11]##, ##REF##17664418##[12]##.</p> <p>CNV determination on the digital array is based upon its ability to partition DNA sequences. Given the number of molecules per panel and the dilution factor, the concentration of the target sequence in a DNA sample can be accurately calculated. In a multiplex PCR reaction with 2 or more assays, multiple genes can be quantitated simultaneously and independently, effectively eliminating any pipetting errors if separate reactions have to be set up for different genes. When a single copy reference gene (RNase P in this study, ##REF##2308839##[10]##) is used in the reaction, the ratio of the target gene to the reference gene would reflect the copy number per haploid genome of the target gene.</p> <title>Primary contribution of this paper</title> <p>In this paper we will show that the digital array provides a robust and easy-to-use platform to study CNVs. We have derived a mathematical framework to calculate the true concentration of molecules from the observed positive reactions in a panel. We also show how to perform statistical analysis to find the 95% confidence intervals of the true concentrations and the ratio of two concentrations in a CNV experiment using the digital array with multiplex PCR.</p> <p>The copy number variation problem can be stated as follows. <italic>Given two counts h</italic>
<sub>1</sub>
<italic>and h</italic>
<sub>2</sub>
<italic>of positive chambers for two genes in a digital array panel, how can one estimate a ratio of true concentrations r = λ</italic>
<sub>1</sub>
<italic>/λ</italic>
<sub>2</sub>
<italic>of the two genes and a confidence interval</italic> [<italic>r<sub>Low</sub></italic>, <italic>r<sub>High</sub></italic>] <italic>on the estimation?</italic>
</p> <p>Our approach is built on well-known tools and techniques from statistics. It decomposes the problem into two parts.</p> <p>Given a count <italic>h</italic> of positive chambers, how can one estimate the true concentration <italic>λ</italic> of target molecules in the DNA sample and a confidence interval [<italic>λ<sub>Low</sub></italic>, <italic>λ<sub>High</sub></italic>] on this estimation?</p> <p>Given estimated true concentrations <italic>λ</italic>
<sub>1</sub> and <italic>λ</italic>
<sub>2</sub> of the reference gene and the target gene, respectively, in the DNA sample and their respective confidence intervals, how can one estimate the ratio <italic>r</italic> = <italic>λ</italic>
<sub>1</sub>/<italic>λ</italic>
<sub>2</sub> and a confidence interval [<italic>r<sub>Low</sub></italic>, <italic>r<sub>High</sub></italic>] on this estimation?</p> <p>It turns out that the first question can be answered by applying sampling and estimation theories from statistics and probability, and the second question can be answered by a numerical algorithm based on generalization of a mathematical theorem.</p> <p>For related work on answering the first question, using Bayesian approach, see unpublished preprint by Warren et al. titled “The Digital Array Response Curve” dated March 2007 at <ext-link ext-link-type="uri" xlink:href="http://thebigone.stanford.edu/papers.htm">http://thebigone.stanford.edu/papers.htm</ext-link>. Warren et al. assumed a uniform probability distribution of number of molecules, with maximum number assumed to be 4000, and using Bayesian and combinatorial methods, presented a solution. The confidence interval obtained using Bayesian probabilistic framework, is often referred to as <italic>credible interval or Bayesian confidence interval</italic> which requires one to incorporate problem-specific contextual information from the prior distribution.</p> <p>This paper differs from this prior work by Warren et al. in two different ways. First, we consider the parameter <italic>λ</italic> to be a fixed constant, unlike having a probability distribution as in Bayesian approach. Second, in addition to providing the answer to the first question, we are interested in estimating the confidence interval of the ratio of two concentrations which is new work. For difference between credible interval and confidence interval, see ##UREF##2##[13]##. Both approaches give good results depending upon the question one is trying to answer.</p> <p>We will prove mathematical correctness of our results in this paper and present simulation results to help the reader build useful insight. Finally, we present actual CNV experiments on the digital array with known ratios and show the results using the techniques developed in this paper.</p>
<title>Materials and Methods</title> <p>A 10-µl reaction mix is normally prepared for each panel. It contains 1× TaqMan Universal master mix (Applied Biosystems, Foster City, CA), 1× RNase P-VIC TaqMan assay, 1× TaqMan assay for the target gene (900 nM primers and 200 nM FAM-labeled probe), 1× sample loading reagent (Fluidigm, South San Francisco, CA) and DNA with about 1,100–1,300 copies of the RNase P gene. 4.59 µl of the 10-µl reaction mix was uniformly partitioned into the 765 reaction chambers of each panel and the digital array was thermocycled on the BioMark system. Thermocycling conditions included a 95°C, 10 minute hot start followed by 40 cycles of two-step PCR: 15 seconds at 95°C for denaturing and 1 minute at 60°C for annealing and extension. Molecules of the two genes were independently amplified. FAM and VIC signals of all chambers were recorded at the end of each PCR cycle. After the reaction was completed, Digital PCR Analysis software (Fluidigm, South San Francisco, CA) was used to process the data and count the numbers of both FAM-positive chambers (target gene) and VIC-positive chambers (RNase P) in each panel.</p> <p>A spike-in experiment was performed using a synthetic construct to explore the digital array's feasibility as a robust platform for the CNV study. A 65-base oligonucleotide was ordered from Integrated DNA Technologies (Coralville, IA) that is identical to a fragment of the human RPP30 gene. The sequences of the primers and FAM-BHQ probe used to amplify this construct are from Emery et al ##REF##15030703##[22]##. The primers and probe were ordered from Biosearch Technologies (Novato, CA).</p> <p>Both RPP30 synthetic construct and human genomic DNA NA10860 (Coriell Cell Repositories Camden, NJ) were quantitated using the RPP30 assay on a digital array. Different amounts of RPP30 synthetic construct was then added into the genomic DNA so that mixtures with ratios of RPP30 to RNase P of 1∶1 (no spike-in), 1∶1.5, 1∶2, 1∶2.5, 1∶3, and 1∶3.5 were made simulating DNA samples containing 2 to 7 copies of the RPP30 gene per diploid cell.</p> <p>These DNA mixtures were analyzed on the digital arrays as described above. Five panels were used for each mixture and 400–500 RNase P molecules were present in each panel. The ratios of RPP30/RNase P of all samples were calculated using the techniques developed in this paper. For each ratio, we did pooled analysis by adding the numbers of positive chambers in the first <italic>P</italic> = 1,2,3,4,5 panels. The results are summarized in the previous section and in ##FIG##7##Figure 8##.</p>
<title>Results</title> <p>DNA quantitation in the digital array is based on the partitioning of a PCR reaction into an array of several hundreds or even few thousands of chambers or wells. One panel of the digital array consists of 765 chambers and one can use up to 12 panels at a time. If the concentration of the target molecules is low in the DNA sample, most of the chambers capture either one or no molecules and the number of positive chambers at the end point of the PCR yields close approximation to the true concentration of the target. However, if the number of molecules is large, then there is greater probability of several molecules being in the same chamber, and therefore the number of positive chambers would be significantly lower compared with the number of molecules in the chambers.</p> <p>We are interested in estimating the true concentration of the molecules in the DNA sample from which we extracted 6 nl×765 = 4.59 µl of sample for each panel.</p> <p>Consider the universe of infinite number of the digital array chambers filled with an infinite amount of the DNA sample where the true concentration of the target molecules is <italic>λ</italic> per chamber (per 6 nl). The true concentration is an unknown population parameter of this infinite DNA sample. If a chamber gets no molecule then it constitutes failure in the sense of Bernoulli experiment. If it gets one or more molecules, that is, if it gets a “hit” and is therefore positive, then it constitutes success. Let the probability of success be <italic>p</italic>. Note that <italic>p</italic> is an unknown population parameter. We will use the standard hat notation to denote sample estimators of population parameters. For example, <italic>pˆ</italic> and  will denote the estimators of <italic>p</italic> and <italic>λ</italic>, respectively.</p> <title>Relationship between <italic>p</italic> and <italic>λ</italic>
</title> <p>One can model <italic>K</italic>, the number of molecules in each chamber as a Poisson process, and this gives the relationship between <italic>p</italic> and <italic>λ</italic> as followsAlternatively, consider <italic>M</italic> molecules randomly distributed in <italic>C</italic> chambers. The probability of any given molecule being in any given chamber is . So the probability <italic>p</italic> of a given chamber having at least one molecule isSince <italic>M</italic> = <italic>λC</italic>, we haveAs number of chambers becomes arbitrarily large, the above approaches <italic>e</italic>
<sup>−<italic>λ</italic></sup>. Therefore,which establishes the relationship between <italic>λ</italic> and <italic>p</italic>.</p> <title>Confidence Intervals for estimation of <italic>p</italic> and <italic>λ</italic>
</title> <p>A chamber getting a hit or no hit is a binomial process, same as toss of a coin, with success probability <italic>p</italic>. Let the number of positive chambers in the panel be <italic>H</italic>. Consider  as an estimator of <italic>p</italic>. It is well known that <italic>pˆ</italic> is an unbiased estimator of <italic>p</italic> and has expectation <italic>p</italic> and standard deviation  and its sampling distribution <italic>f</italic>(<italic>pˆ</italic>) is approximately normal for large <italic>C</italic>. See ##FIG##2##Figure 3## for illustration of the above ideas.</p> <p>See ##UREF##2##[13]##, ##UREF##3##[14]##, ##UREF##4##[15]## for extensive literature on obtaining confidence interval for the estimation of binomial probability. It is referred to as binomial sign test when the test statistic can be approximated with the chi-square distribution, specifically through the use of the chi-square goodness-of-fit. An alternative and equivalent approximation is obtained by using the normal distribution and then the test is referred to as the population proportion test, see ##UREF##4##[15]##.</p> <p>If <italic>C</italic> is large enough, then the confidence limits are approximately given byFor 95% confidence interval, <italic>z<sub>c</sub></italic> = 1.96. For the digital array, <italic>C</italic> is an integral multiple of 765 and is comfortably large enough for the above approximation.</p> <p>Define the estimator of <italic>λ</italic> asSince the probabilities in any given differential area of a probability density function are preserved under change of variables, the 95% confidence interval [, ] is directly given as follows See ##FIG##3##Figure 4## for illustration.</p> <title>Expectation of estimation of <italic>λ</italic>
</title> <p>Let a random variable <italic>X</italic> have probability density function <italic>f<sub>X</sub></italic>(<italic>x</italic>). If <italic>h</italic>(<italic>x</italic>) is either increasing or decreasing in <italic>x</italic>, then <italic>U</italic> = <italic>h</italic>(<italic>X</italic>) has density function given bywhich follows from the fact that probabilities in any given differential area have to be invariant under change of variables, see ##UREF##5##[16]##. Furthermore,which can be expanded using Taylor series expansion of <italic>h</italic>(<italic>x</italic>) around the mean <italic>ξ</italic> = <italic>E</italic>(<italic>X</italic>), as followsSince in our case, we have the followingtherefore, in above, we have <italic>x</italic> = <italic>pˆ</italic>,  and <italic>h</italic>(<italic>x</italic>) = −ln(1−<italic>x</italic>). Since  is a monotonically increasing function of <italic>pˆ</italic>, one can get the sampling distribution of  from the sampling distribution of <italic>pˆ</italic> asNote that due to nonlinear relationship between  and <italic>pˆ</italic>, one can not make assumptions about <italic>g</italic>. In general, <italic>g</italic> is not normal and .</p> <p>Now we derive an approximation for  from the Taylor series expansion shown above. Higher order central moments of Gaussian function <italic>f</italic>(<italic>pˆ</italic>) with mean <italic>p</italic> areFor proof see ##UREF##6##[17]##. Since <italic>f</italic>(<italic>pˆ</italic>) has very small <italic>σ</italic> due to very large number of chambers, the higher order terms for all <italic>n</italic>&gt;0 in the Taylor expansion are small,and therefore the only contributing term is when <italic>n</italic> = 0, which impliesAs <italic>C</italic>→∞,  approaches the true concentration <italic>λ</italic> of molecules.</p> <title>Simulation Results on Estimation of <italic>p</italic>
</title> <p>It is informative and useful to run a simulation experiment on the computer to see how the real-world matches with the theory developed above. For this purpose, one can use a random number generator and a computer program to simulate the universe of the digital array chambers.</p> <p>If a panel has <italic>C</italic> chambers, consider a universe of <italic>C</italic>×<italic>K</italic> many chambers where <italic>K</italic> is a large number chosen for simulation. Choose some value of <italic>λ</italic> as the true concentration of molecules in one chamber. Therefore, in total, there will be <italic>λ</italic>×<italic>C</italic>×<italic>K</italic> molecules. Assign each of these molecules randomly to one of the chambers.</p> <p>Extract <italic>K</italic> panels out of this universe and for each of the panels, compute  as an estimator of <italic>p</italic> and plot its histogram over all the <italic>K</italic> panels. The mean should be <italic>p</italic> = 1−<italic>e</italic>
<sup>−<italic>λ</italic></sup> and standard deviation should be . For each of these panels, estimate <italic>λ</italic> and compute the 95% confidence interval. In 95% of the <italic>K</italic> panels, the true value of <italic>λ</italic> should lie within the confidence interval.</p> <p>For our simulation experiments we chose <italic>M</italic> = 400, that is, . We chose <italic>K</italic> = 70000. In ##FIG##4##Figure 5## we show the histogram of <italic>H</italic> which is really same as distribution of <italic>P</italic> scaled by a factor of 765. In ##TAB##0##Table 1## we show how the predicted values match with the actual simulation values. In the same way, the sampling distribution of number of molecules matched with what is predicted by theory. Though the results of the simulation follow from elementary probability, we conducted these simulations in order to build more advanced simulations for ratios of concentrations later. They also illustrate the meaning of the confidence interval.</p> <title>Determination of Ratio of Concentrations</title> <p>In previous section, we established a method for estimating the true concentration of the target molecules in the DNA sample from the count of positive chambers as well as the 95% confidence interval for this estimation. We also showed how the sampling distribution  is related to the sampling distribution <italic>f</italic>(<italic>pˆ</italic>).</p> <p>In CNV, the goal is to determine ratio of true concentrations of two genes, one being reference gene and the other being test gene, and associated confidence interval, which we now accomplish in next subsections.</p> <title>Fieller's Theorem and its geometric interpretation</title> <p>Let the sampling distributions of the test gene and the reference gene be  and , respectively. If these distributions were normal, then one can make use of Fieller's Theorem ##UREF##7##[18]##, ##UREF##8##[19]##.</p> <p>However, as mentioned in previous section, one can not make this assumption in general. It is useful to go through the geometric interpretation of Fieller's theorem so that one can solve the problem for arbitrary sampling distributions. See ##FIG##5##Figure 6## for geometric interpretation of Fieller's Theorem ##UREF##9##[20]##, ##UREF##10##[21]##.</p> <p>Assume  and  are normal. For  and , the ratio  can be shown as the slope of the line in the two-dimensional plane which passes through the origin and the 2-D point (, ). Luxburg et al. show in ##UREF##9##[20]##, ##UREF##10##[21]## how a confidence ellipse in the two-dimensional plane can be constructed. Consider the two lines which pass through the origin and are tangents to this ellipse. The intersection of these lines with the vertical line at  gives the desired confidence interval.</p> <title>Ratio of concentrations</title> <p>In this paper we have presented data in a controlled experimental system, where a synthetic DNA construct was spiked into human cell line DNA at different concentrations. In this case, the synthetic construct, which was to the RPP30 gene, was used as the target, and the RNase P gene which was endogenous to the human cell line, was used as the reference gene. The two genes were identified using two separate PCR reactions, using separate PCR primers and probes. Since there is no reason to assume that the amplification and detection of the target and reference genes are linked,  and  are independent variables.</p> <p>It is easy to see from the proof of Fieller's theorem and its geometric interpretation that one can compute sampling distribution <italic>q</italic> of the ratio estimator  as follows:This can be interpreted as cutting out thin wedges in the joint distribution of  and  and accumulating the probabilities inside the wedge to compute the function <italic>q</italic> in the corresponding thin interval of the ratio. This is the basis of our numerical algorithm which implements integration in order to compute <italic>q</italic>(<italic>rˆ</italic>):</p> <p>Build histograms of sampling distributions  and . The tails of the histograms where probabilities become very small are approximated by zero.</p> <p>Build a histogram of sampling distribution <italic>q</italic>(<italic>rˆ</italic>) of  by considering each bin [<italic>r</italic>
<sub>1</sub>, <italic>r</italic>
<sub>2</sub>] and by adding all the joint probabilities of different values of concentrations which give a ratio <italic>rˆ</italic>∈[<italic>r</italic>
<sub>1</sub>, <italic>r</italic>
<sub>2</sub>].</p> <p>Compute the mean and the 95% confidence interval from the ratio histogram.</p> <p>See ##FIG##6##Figure 7## for illustration of the above algorithm.</p> <p>One can still use direct formulas, as an approximation, to compute confidence interval as follows.</p> <p>The means of  and  are <italic>λ</italic>
<sub>1</sub> and <italic>λ</italic>
<sub>2</sub> respectively. Let the standard deviations be <italic>σ<sub>x</sub></italic> and <italic>σ<sub>y</sub></italic> respectively. For given estimations  and , assuming that distributions are normal, it follows from the analysis in ##UREF##9##[20]##, ##UREF##10##[21]## that the boundary of the confidence ellipse for a given confidence level <italic>z<sub>c</sub></italic> would be defined byIt is easy to generalize this to the case, under a reasonably close approximation, when we have asymmetric distributions which are assumed to be normal in each of the four quadrants of the coordinate system centered at (, ). Then the confidence region is made of union of four quadrant-wise elliptic regions.</p> <p>Let the asymmetric confidence intervals for specified <italic>z<sub>c</sub></italic> and the two concentrations be [, ] and [, ].</p> <p>If <italic>W<sub>R</sub></italic> = <italic>W<sub>L</sub></italic> = <italic>z<sub>c</sub>σ<sub>x</sub></italic> and <italic>H<sub>T</sub></italic> = <italic>H<sub>B</sub></italic> = <italic>z<sub>c</sub>σ<sub>y</sub></italic>, it is symmetric case ##UREF##9##[20]##, ##UREF##10##[21]##. Using simple algebraic manipulations, it can be shown, as in symmetric case, that the slopes of lines that will be tangents to this union of four quadrant-wise ellipses will beThe above equations can be used as an approximation though numerical algorithm will give more accurate results as the algorithm does not make any assumptions and works with arbitrary sampling distributions.</p> <p>One detail has to be mentioned. Special care has to be taken if the confidence region gets too close to  axis when  is small. If it touches  axis, then <italic>r<sub>High</sub></italic> = ∞. If either  or  is too small, one can build respective histogram with smaller bin size to get more accurate results.</p> <p>See ##TAB##1##Table 2## for summary of equations derived in order to solve the copy number variation problem. Though the numerical approach based on histograms is recommended as it does not make assumptions, these direct formulas can be used as close approximation.</p> <title>Simulation Results on Estimation of Ratio</title> <p>We conducted simulation studies, using a random number generator and a computer program as in previous section, by choosing a ratio of 2 of concentrations of two genes, which are independent of each other, and building a distribution of estimated ratios over 50 thousand panels. In 94.9% of the panels, the true chosen ratio did lie in the computed confidence intervals thereby showing the correctness of our mathematical analysis.</p> <title>Actual Results on Estimation of Ratio</title> <p>The copy number variation results for known ratios of 1, 1.5, 2, 2.5, 3 and 3.5 are shown in ##FIG##7##Figure 8##. <xref ref-type="sec" rid="s3">Materials and methods</xref> for this experiment are discussed in next section. As the number of panels <italic>P</italic> increases, then the number of chambers <italic>C</italic> = 765 <italic>P</italic>increases and therefore the estimation of the ratio becomes more accurate as well as the confidence interval shrinks. When only 1 panel is used, there is significant overlap between 95% confidence intervals of certain ratios e.g. between ratio 2 and 2.5. There is no overlap when 3 or more panels are used. In all cases the known ratio lies within the computed 95% confidence interval. Note that using mathematical analysis one can find optimal numbers of positive chambers for each ratio which give smallest confidence intervals and which will therefore improve the results.</p> <p>In summary, Fluidigm's digital array is capable of accurately quantitating DNA samples and is a valuable platform for studying copy number variation. It is a robust technology that is sequence-specific, easy-to-use, and extremely flexible. We have presented mathematical and algorithmic solutions to analyze CNV on a digital array. The solution is an elegant application of statistical sampling and estimation theories to such an important real-world problem. We have shown how one can compute the true concentration of a target sequence in a DNA sample and the associated confidence interval on this estimation, and how one can compute the ratio of true concentrations of multiple sequences and the associated confidence interval on the estimation of this ratio.</p>
<p>Conceived and designed the experiments: JQ RR. Performed the experiments: JQ. Analyzed the data: SD. Wrote the paper: SD JQ RR. Worked out the mathematical analysis and algorithms and conducted simulation experiments: SD.</p> <p>Copy Number Variations (CNVs) of regions of the human genome have been associated with multiple diseases. We present an algorithm which is mathematically sound and computationally efficient to accurately analyze CNV in a DNA sample utilizing a nanofluidic device, known as the digital array. This numerical algorithm is utilized to compute copy number variation and the associated statistical confidence interval and is based on results from probability theory and statistics. We also provide formulas which can be used as close approximations.</p>
<p>The authors would like to thank the members of the Research and Development division of Fluidigm Corporation, in particular the software and molecular biology teams for facilitating this work. The first author would like to thank Gang Sun for bringing the problem to his attention, and to Gang Sun and Hou-Pu Chou for discussions leading to clarification of the problem statement and for initial feedback on the solution. The authors would also like thank an anonymous referee for constructive feedback on the paper which led to its significant improvement.</p>
<fig id="pone-0002876-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002876.g001</object-id><label>Figure 1</label><caption><title>A digital array has 12 panels of 765 reaction chambers each.</title><p>PCR mixes are loaded into each panel and single DNA molecules are randomly partitioned into the chambers. The digital array can be thermocycled, imaged on a BioMark instrument, and the data analyzed using the Digital PCR Analysis software.</p></caption></fig> <fig id="pone-0002876-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002876.g002</object-id><label>Figure 2</label><caption><title>Human genomic DNA NA10860 (left 5 panels) and the RPP30 synthetic construct (right 5 panels) were quantitated using the RPP30 (FAM) assay on this digital array.</title><p>The two bottom panels are NTC (no template control). Digital PCR Analysis software can count the number of positive chambers in each panel. When two assays with two fluorescent dyes are used in a multiplex digital PCR reaction, two genes can be independently quantitated. This is the basis of the CNV study using the digital array.</p></caption></fig> <fig id="pone-0002876-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002876.g003</object-id><label>Figure 3</label><caption><title>Consider an infinite universe of chambers.</title><p>A digital array panel is a finite sampling of this universe. The goal is to determine <italic>λ</italic>, the mean number of the target molecules per chamber in the DNA sample. The number of positive chambers, which have hits of one or more molecules, shown as filled green squares in the panel with <italic>C</italic>( = 765) chambers is <italic>H</italic>.</p></caption></fig> <fig id="pone-0002876-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002876.g004</object-id><label>Figure 4</label><caption><title>From the sampling distribution of estimation of p, one can obtain the sampling distribution of estimation of λ.</title></caption></fig> <fig id="pone-0002876-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002876.g005</object-id><label>Figure 5</label><caption><title>Histogram of number of positive chambers <italic>H</italic> = <italic>P</italic>×<italic>C</italic> obtained by choosing <italic>M</italic> = 400 as the mean number of molecules per panel over 70 thousand panels and running a simulation using a random number generator.</title><p>The green curve is the sampling distribution predicted by the theory.</p></caption></fig> <fig id="pone-0002876-g006" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002876.g006</object-id><label>Figure 6</label><caption><title>Geometric interpretation of Fieller's Theorem to compute confidence interval of ratio of two normally distributed random variables  and  in which confidence ellipse of the joint sampling distribution is projected on a vertical line.</title></caption></fig> <fig id="pone-0002876-g007" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002876.g007</object-id><label>Figure 7</label><caption><title>Illustration of a numerical projection algorithm to compute the sampling distribution of ratio of two random variables with arbitrary probability distributions by slicing the 2-D space into thin wedges and accumulating the joint probabilities in the wedges.</title><p>Most of the contribution would come from the confidence ellipse region.</p></caption></fig> <fig id="pone-0002876-g008" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002876.g008</object-id><label>Figure 8</label><caption><title>Results of actual CNV experiments on the digital array with varying number of copies of the target gene.</title><p>In total, 6 different known ratios were estimated by running the experiments for varying number of panels. The graphs for different numbers of copies are slightly staggered to allow visual comparison of overlap of the 95% confidence intervals.</p></caption></fig>
<table-wrap id="pone-0002876-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002876.t001</object-id><label>Table 1</label><caption><title>Comparison of the metrics of histogram, shown in ##FIG##4##Figure 5##, of number of positive chambers obtained in simulation with those predicted by the theory.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Theoretical Predictions</td><td align="left" rowspan="1" colspan="1">Simulation Results</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Mean</td><td align="left" rowspan="1" colspan="1">311.5</td><td align="left" rowspan="1" colspan="1">311.48</td></tr><tr><td align="left" rowspan="1" colspan="1">Standard Deviation</td><td align="left" rowspan="1" colspan="1">13.59</td><td align="left" rowspan="1" colspan="1">13.58</td></tr><tr><td align="left" rowspan="1" colspan="1">Percent of times <italic>M</italic> = 400 lies in the computed 95% confidence interval</td><td align="left" rowspan="1" colspan="1">95%</td><td align="left" rowspan="1" colspan="1">94.44%</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0002876-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002876.t002</object-id><label>Table 2</label><caption><title>Given number of chambers <italic>C</italic> and counts <italic>H</italic>
<sub>1</sub> and <italic>H</italic>
<sub>2</sub> of the positive chambers in a digital array for the target gene and the reference gene, respectively, list of formulas needed to analyze copy number variation.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/></colgroup><tbody><tr><td align="left" rowspan="1" colspan="1">
, 
</td></tr><tr><td align="left" rowspan="1" colspan="1">
, 
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>pˆ</italic>
<sub>1,<italic>Low</italic></sub> = <italic>pˆ</italic>
<sub>1</sub>−1.96<italic>S</italic>
<sub>1</sub>, <italic>pˆ</italic>
<sub>1,<italic>High</italic></sub> = <italic>pˆ</italic>
<sub>1</sub>+1.96<italic>S</italic>
<sub>1</sub>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>pˆ</italic>
<sub>2,<italic>Low</italic></sub> = <italic>pˆ</italic>
<sub>2</sub>−1.96<italic>S</italic>
<sub>2</sub>, <italic>pˆ</italic>
<sub>2,<italic>High</italic></sub> = <italic>pˆ</italic>
<sub>2</sub>+1.96<italic>S</italic>
<sub>2</sub>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
, , 
</td></tr><tr><td align="left" rowspan="1" colspan="1">
, , 
</td></tr><tr><td align="left" rowspan="1" colspan="1">
, 
</td></tr><tr><td align="left" rowspan="1" colspan="1">
, 
</td></tr><tr><td align="left" rowspan="1" colspan="1">

</td></tr><tr><td align="left" rowspan="1" colspan="1">

</td></tr><tr><td align="left" rowspan="1" colspan="1">

</td></tr></tbody></table></alternatives></table-wrap>
<inline-formula></inline-formula> <disp-formula></disp-formula> <inline-formula></inline-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <disp-formula></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <disp-formula></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <disp-formula></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <disp-formula></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula>
<table-wrap-foot><fn id="nt101"><p>See the details in the paper for assumptions made so that these equations are close approximations to actual values.</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>The authors are employees of Fludigim Corporation and for this research project they used a technology platform which is being commercialized by Fluidigm.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>The authors are employees of Fludigim Corporation and this work was undertaken at Fluidigm.</p></fn></fn-group>
<graphic xlink:href="pone.0002876.g001"/> <graphic xlink:href="pone.0002876.g002"/> <inline-graphic xlink:href="pone.0002876.e001.jpg" mimetype="image"/> <graphic xlink:href="pone.0002876.e002.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pone.0002876.e003.jpg" mimetype="image"/> <graphic xlink:href="pone.0002876.e004.jpg" mimetype="image" position="float"/> <graphic xlink:href="pone.0002876.e005.jpg" mimetype="image" position="float"/> <graphic xlink:href="pone.0002876.e006.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pone.0002876.e007.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e008.jpg" mimetype="image"/> <graphic xlink:href="pone.0002876.g003"/> <graphic xlink:href="pone.0002876.e009.jpg" mimetype="image" position="float"/> <graphic xlink:href="pone.0002876.e010.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pone.0002876.e011.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e012.jpg" mimetype="image"/> <graphic xlink:href="pone.0002876.e013.jpg" mimetype="image" position="float"/> <graphic xlink:href="pone.0002876.g004"/> <graphic xlink:href="pone.0002876.e014.jpg" mimetype="image" position="float"/> <graphic xlink:href="pone.0002876.e015.jpg" mimetype="image" position="float"/> <graphic xlink:href="pone.0002876.e016.jpg" mimetype="image" position="float"/> <graphic xlink:href="pone.0002876.e017.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pone.0002876.e018.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e019.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e020.jpg" mimetype="image"/> <graphic xlink:href="pone.0002876.e021.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pone.0002876.e022.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e023.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e024.jpg" mimetype="image"/> <graphic xlink:href="pone.0002876.e025.jpg" mimetype="image" position="float"/> <graphic xlink:href="pone.0002876.e026.jpg" mimetype="image" position="float"/> <graphic xlink:href="pone.0002876.e027.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pone.0002876.e028.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e029.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e030.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e031.jpg" mimetype="image"/> <graphic xlink:href="pone.0002876.g005"/> <graphic id="pone-0002876-t001-1" xlink:href="pone.0002876.t001"/> <inline-graphic xlink:href="pone.0002876.e032.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e033.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e034.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e035.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e036.jpg" mimetype="image"/> <graphic xlink:href="pone.0002876.g006"/> <inline-graphic xlink:href="pone.0002876.e037.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e038.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e039.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e040.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e041.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e042.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e043.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e044.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e045.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e046.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e047.jpg" mimetype="image"/> <graphic xlink:href="pone.0002876.e048.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pone.0002876.e049.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e050.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e051.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e052.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e053.jpg" mimetype="image"/> <graphic xlink:href="pone.0002876.g007"/> <inline-graphic xlink:href="pone.0002876.e054.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e055.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e056.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e057.jpg" mimetype="image"/> <graphic xlink:href="pone.0002876.e058.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pone.0002876.e059.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e060.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e061.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e062.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e063.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e064.jpg" mimetype="image"/> <graphic xlink:href="pone.0002876.e065.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pone.0002876.e066.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e067.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e068.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e069.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e070.jpg" mimetype="image"/> <graphic id="pone-0002876-t002-2" xlink:href="pone.0002876.t002"/> <inline-graphic xlink:href="pone.0002876.e071.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e072.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e073.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e074.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e075.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e076.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e077.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e078.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e079.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e080.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e081.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e082.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e083.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e084.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e085.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e086.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002876.e087.jpg" mimetype="image"/> <graphic xlink:href="pone.0002876.g008"/>
[{"label": ["1"], "element-citation": ["\n"], "surname": ["Vogelstein", "Kinzler"], "given-names": ["B", "KW"], "year": ["1999"], "article-title": ["Digital PCR. Proc. Natl. Acad."], "source": ["Sci U S A"], "volume": ["96"], "fpage": ["9236"], "lpage": ["9241"]}, {"label": ["3"], "element-citation": ["\n"], "surname": ["Sindelka", "Jonak", "Hands", "Bustin", "Kubista"], "given-names": ["R", "J", "R", "SA", "M"], "year": ["2007"], "article-title": ["Intracellular expression profiles measured by real-time PCR tomography in the Xenopus laevis oocyte."], "source": ["Nucleic Acids Research"], "fpage": ["1"], "lpage": ["6"], "comment": ["doi:10.1093/nar/gkm1024"]}, {"label": ["13"], "element-citation": ["\n"], "surname": ["Kapadia", "Chan", "Moy\u00e9"], "given-names": ["AS", "W", "LA"], "year": ["2005"], "source": ["Mathematical Statistics with Applications"], "publisher-name": ["CRC Press"], "fpage": ["617"]}, {"label": ["14"], "element-citation": ["\n"], "surname": ["Motulsky"], "given-names": ["H"], "year": ["1995"], "source": ["Intuitive Biostatistics"], "publisher-loc": ["New York"], "publisher-name": ["Oxford University Press"], "fpage": ["408"]}, {"label": ["15"], "element-citation": ["\n"], "surname": ["Sheskin"], "given-names": ["DJ"], "year": ["2000"], "source": ["Handbook of Parametric and Nonparametric Statistical Procedures"], "publisher-loc": ["New York"], "publisher-name": ["Chapman & Hall/CRC Press"], "fpage": ["982"]}, {"label": ["16"], "element-citation": ["\n"], "surname": ["Scheaffer"], "given-names": ["RL"], "year": ["1990"], "source": ["Introduction to Probability and its Applications."], "publisher-name": ["PWS-KENT Publishing Company"], "fpage": ["363"]}, {"label": ["17"], "element-citation": ["\n"], "surname": ["Papoulis"], "given-names": ["A"], "year": ["1965"], "source": ["Probability, Random Variables, and Stochastic Processes"], "publisher-loc": ["New York"], "publisher-name": ["McGraw-Hill"], "fpage": ["583"]}, {"label": ["18"], "element-citation": ["\n"], "surname": ["Fieller"], "given-names": ["E"], "year": ["1932"], "article-title": ["The distribution of the index in a normal bivariate population."], "source": ["Biometrika"], "volume": ["24(3/4)"], "fpage": ["428"], "lpage": ["440"]}, {"label": ["19"], "element-citation": ["\n"], "surname": ["Fieller"], "given-names": ["E"], "year": ["1954"], "article-title": ["Some problems in interval estimation."], "source": ["J Roy Statist Soc Ser: B"], "volume": ["16"], "fpage": ["175"], "lpage": ["185"]}, {"label": ["20"], "element-citation": ["\n"], "surname": ["von Luxburg", "Franz"], "given-names": ["U", "VH"], "year": ["2002"], "article-title": ["Confidence Sets for Ratios: A Purely Geometric Approach To Fieller's Theorem."], "comment": ["Max Plank Institute for Biological Cybernetics, Technical Report No. TR-133"]}, {"label": ["21"], "element-citation": ["\n"], "surname": ["von Luxburg", "Franz"], "given-names": ["U", "VH"], "year": ["2007"], "article-title": ["A Geometric Approach to Confidence Sets for Ratios: Fieller's Theorem, Generalizations, and Boostrap."], "comment": ["ArXiv e-prints, arXiv:0711.0198v1. See: "], "ext-link": ["http://arxiv.org/abs/0711.0198"]}]
PLoS One. 2008 Aug 6; 3(8):e2876