Abstract
Current methods used to analyze real-time quantitative polymerase chain reaction (qPCR) data exhibit systematic deviations from the assumed model over the progression of the reaction. Slight variations in the amount of the initial target molecule or in early amplifications are likely responsible for these deviations. Commonly used 4- and 5-parameter sigmoidal models appear to be particularly susceptible to this issue, often displaying patterns of autocorrelation in the residuals. The presence of this phenomenon, even for technical replicates, suggests that these parametric models may be misspecified. Specifically, they do not account for the sequential dependent nature of the amplification process that underlies qPCR fluorescence measurements. We demonstrate that a Smooth Transition Autoregressive (STAR) model addresses this limitation by explicitly modeling the dependence between cycles and the gradual transition between amplification regimes. In summary, application of a STAR model to qPCR amplification data improves model fit and reduces autocorrelation in the residuals.
Benjamin Hsu
MS Student in Actuarial Science at Columbia University
Valeriia Sherina
Principal Statistician at GSK
