Autoregressive modeling and diagnostics for qPCR amplification.
Ontology highlight
ABSTRACT: Motivation
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.Results
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.Availability and implementation
R scripts to reproduce all the analyses and results described in this manuscript can be found at: https://github.com/bhsu4/GAPDH.SO.Supplementary information
Supplementary data are available at Bioinformatics online.
SUBMITTER: Hsu B
PROVIDER: S-EPMC8016497 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Measurement error is omnipresent in psychological data. However, the vast majority of applications of autoregressive time series analyses in psychology do not take measurement error into account. Disregarding measurement error when it is present in the data results in a bias of the autoregressive parameters. We discuss two models that take measurement error into account: An autoregressive model with a white noise term (AR+WN), and an autoregressive moving average (ARMA) model. In a simulation study we compare the parameter recovery performance of these models, and compare this performance for both a Bayesian and frequentist approach. We find that overall, the AR+WN model performs better. Furthermore, we find that for realistic (i.e., small) sample sizes, psychological research would benefit from a Bayesian approach in fitting these models. Finally, we illustrate the effect of disregarding measurement error in an AR(1) model by means of an empirical application on mood data in women. We find that, depending on the person, approximately 30-50% of the total variance was due to measurement error, and that disregarding this measurement error results in a substantial underestimation of the autoregressive parameters.
| S-EPMC4516825 | biostudies-literature
Project description:Every year, Agroscope examines nearly 300,000 tubers for the presence of viruses, as regulated for the certification of seed potatoes intended for Swiss growers. Since 2016, this examination has been performed via RT-qPCR on dormant tubers directly after harvest. This method offers fast results and eliminates the need for the use of Rindite, which is a toxic and polluting gaseous compound previously used in Switzerland to break the dormancy of seed tubers. The implementation of this molecular analytical method for the routine diagnosis of regulated viruses makes it possible to conduct additional analyses via Illumina sequencing to assess the conformity of the primers and probes used with the sequences of the different viral isolates. This form of quality control in routine diagnosis is a source of information that can answer more fundamental scientific questions related to the epidemiology of viral strains related to certification. The datasets produced in this framework can also be used to explore the diversity of rare or unknown virus species in potato crops. Supplementary Information The online version contains supplementary material available at 10.1007/s11540-021-09491-3.
| S-EPMC8572825 | biostudies-literature
Project description:Predicting health outcomes from longitudinal health histories is of central importance to healthcare. Observational healthcare databases such as patient diary databases provide a rich resource for patient-level predictive modeling. In this paper, we propose a Bayesian hierarchical vector autoregressive (VAR) model to predict medical and psychological conditions using multivariate time series data. Compared to the existing patient-specific predictive VAR models, our model demonstrated higher accuracy in predicting future observations in terms of both point and interval estimates due to the pooling effect of the hierarchical model specification. In addition, by adopting an elastic-net prior, our model offers greater interpretability about the associations between variables of interest on both the population level and the patient level, as well as between-patient heterogeneity. We apply the model to two examples: 1) predicting substance use craving, negative affect and tobacco use among college students, and 2) predicting functional somatic symptoms and psychological discomforts.
| S-EPMC6294362 | biostudies-literature
Project description:Generative models emerge as promising candidates for novel sequence-data driven approaches to protein design, and for the extraction of structural and functional information about proteins deeply hidden in rapidly growing sequence databases. Here we propose simple autoregressive models as highly accurate but computationally efficient generative sequence models. We show that they perform similarly to existing approaches based on Boltzmann machines or deep generative models, but at a substantially lower computational cost (by a factor between 102 and 103). Furthermore, the simple structure of our models has distinctive mathematical advantages, which translate into an improved applicability in sequence generation and evaluation. Within these models, we can easily estimate both the probability of a given sequence, and, using the model's entropy, the size of the functional sequence space related to a specific protein family. In the example of response regulators, we find a huge number of ca. 1068 possible sequences, which nevertheless constitute only the astronomically small fraction 10-80 of all amino-acid sequences of the same length. These findings illustrate the potential and the difficulty in exploring sequence space via generative sequence models.
| S-EPMC8490405 | biostudies-literature
Project description:Naturally derived bioactive peptides with antihypertensive activities serve as promising alternatives to pharmaceutical drugs. There are few relevant reports on the mapping relationship between the EC50 value of antihypertensive peptide activity (AHTPA-EC50) and its corresponding amino acid sequence (AAS) at present. In this paper, we have constructed two group series based on sorting natural logarithm of AHTPA-EC50 or sorting its corresponding AAS encoding number. One group possesses two series, and we find that there must be a random number series in any group series. The random number series manifests fractal characteristics, and the constructed series of sorting natural logarithm of AHTPA-EC50 shows good autocorrelation characteristics. Therefore, two non-linear autoregressive models with exogenous input (NARXs) were established to describe the two series. A prediction method is further designed for AHTPA-EC50 prediction based on the proposed model. Two dynamic neural networks for NARXs (NARXNNs) are designed to verify the two series characteristics. Dipeptides and tripeptides are used to verify the proposed prediction method. The results show that the mean square error (MSE) of prediction is about 0.5589 for AHTPA-EC50 prediction when the classification of AAS is correct. The proposed method provides a solution for AHTPA-EC50 prediction.
| S-EPMC8787326 | biostudies-literature
Project description:Influenza poses a serious health threat and creates an economic burden for people around the world. The accurate diagnosis of influenza is critical to the timely clinical treatment of patients and the control of outbreaks to protect public health. Commercially available rapid influenza diagnostic tests (RIDTs) that are operated by visual readout are widely used in clinics to screen influenza infections, but RIDTs suffer from imperfect analytical sensitivity, especially when the virus concentration in the sample is low. Fortunately, the sensitivity can be simply improved through an add-on signal amplification step, i.e., thermal contrast amplification (TCA). To demonstrate the advantage of TCA for influenza diagnosis, we conducted a prospective cohort study on 345 clinical specimens collected for influenza A and B testing during the 2017-2018 influenza season. All samples were tested using the Quidel QuickVue Influenza A + B test, followed by a TCA readout, and then confirmatory polymerase chain reaction testing. Through the TCA detecting sub-visual weak positives, TCA reading improved the overall influenza sensitivity by 53% for influenza A and 33% for influenza B over the visual RIDTs readings. Even though the specificity was compromised slightly by the TCA protocol (relative decrease of 0.09% for influenza A and 0.01% for influenza B), the overall performance was still better than that achieved by visual readout based on comparison of their plots in receiver operating characteristic space and F1 scores (relative increase of 14.5% for influenza A and 12.5% for influenza B). Performing a TCA readout on wet RIDTs also improved the overall TCA performance (relative increase in F1 score of 48%). Overall, the TCA method is a simple and promising way to improve the diagnostic performance of commercial RIDTs for infectious diseases, especially in the case of specimens with low target analytes.
| S-EPMC7999055 | biostudies-literature
Project description:The reverse transcription quantitative polymerase chain reaction (RT-qPCR) is an established tool for the diagnosis of RNA pathogens. Its potential for automation has caused it to be used as a presence/absence diagnostic tool even when RNA quantification is not required. This technology has been pushed to the forefront of public awareness by the COVID-19 pandemic, as its global application has enabled rapid and analytically sensitive mass testing, with the first assays targeting three viral genes published within days of the publication of the SARS-CoV-2 genomic sequence. One of those, targeting the RNA-dependent RNA polymerase gene, has been heavily criticised for supposed scientific flaws at the molecular and methodological level, and this criticism has been extrapolated to doubts about the validity of RT-qPCR for COVID-19 testing in general. We have analysed this assay in detail, and our findings reveal some limitations but also highlight the robustness of the RT-qPCR methodology for SARS-CoV-2 detection. Nevertheless, whilst our data show that some errors can be tolerated, it is always prudent to confirm that the primer and probe sequences complement their intended target, since, when errors do occur, they may result in a reduction in the analytical sensitivity. However, in this case, it is unlikely that a mismatch will result in poor specificity or a significant number of false-positive SARS-CoV-2 diagnoses, especially as this is routinely checked by diagnostic laboratories as part of their quality assurance.
| S-EPMC8395416 | biostudies-literature
Project description:Physiological time series are affected by many factors, making them highly nonlinear and nonstationary. As a consequence, heart rate time series are often considered difficult to predict and handle. However, heart rate behavior can indicate underlying cardiovascular and respiratory diseases as well as mood disorders. Given the importance of accurate modeling and reliable predictions of heart rate fluctuations for the prevention and control of certain diseases, it is paramount to identify models with the best performance in such tasks. The objectives of this study were to compare the results of three different forecasting models (Autoregressive Model, Long Short-Term Memory Network, and Convolutional Long Short-Term Memory Network) trained and tested on heart rate beats per minute data obtained from twelve heterogeneous participants and to identify the architecture with the best performance in terms of modeling and forecasting heart rate behavior. Heart rate beats per minute data were collected using a wearable device over a period of 10 days from twelve different participants who were heterogeneous in age, sex, medical history, and lifestyle behaviors. The goodness of the results produced by the models was measured using both the mean absolute error and the root mean square error as error metrics. Despite the three models showing similar performance, the Autoregressive Model gave the best results in all settings examined. For example, considering one of the participants, the Autoregressive Model gave a mean absolute error of 2.069 (compared to 2.173 of the Long Short-Term Memory Network and 2.138 of the Convolutional Long Short-Term Memory Network), achieving an improvement of 5.027% and 3.335%, respectively. Similar results can be observed for the other participants. The findings of the study suggest that regardless of an individual's age, sex, and lifestyle behaviors, their heart rate largely depends on the pattern observed in the previous few minutes, suggesting that heart rate can be reasonably regarded as an autoregressive process. The findings also suggest that minute-by-minute heart rate prediction can be accurately performed using a linear model, at least in individuals without pathologies that cause heartbeat irregularities. The findings also suggest many possible applications for the Autoregressive Model, in principle in any context where minute-by-minute heart rate prediction is required (arrhythmia detection and analysis of the response to training, among others).
| S-EPMC8747593 | biostudies-literature
Project description:To model data from multi-item scales, many researchers default to a confirmatory factor analysis (CFA) approach that restricts cross-loadings and residual correlations to zero. This often leads to problems of measurement-model misfit while also ignoring theoretically relevant alternatives. Existing research mostly offers solutions by relaxing assumptions about cross-loadings and allowing residual correlations. However, such approaches are critiqued as being weak on theory and/or indicative of problematic measurement scales. We offer a theoretically-grounded alternative to modeling survey data called an autoregressive confirmatory factor analysis (AR-CFA), which is motivated by recognizing that responding to survey items is a sequential process that may create temporal dependencies among scale items. We compare an AR-CFA to other common approaches using a sample of 8,569 people measured along five common personality factors, showing how the AR-CFA can improve model fit and offer evidence of increased construct validity. We then introduce methods for testing AR-CFA hypotheses, including cross-level moderation effects using latent interactions among stable factors and time-varying residuals. We recommend considering the AR-CFA as a useful complement to other existing approaches and treat AR-CFA limitations.
| S-EPMC6763968 | biostudies-literature
Project description:Gait asymmetry is often observed in populations with varying degrees of neuromuscular control. While changes in vertical ground reaction force (vGRF) peak magnitude are associated with altered limb loading that can be observed during asymmetric gait, the challenge is identifying techniques with the sensitivity to detect these altered movement patterns. Autoregressive (AR) modeling has successfully delineated between healthy and pathological gait during running; but has been little explored in walking. Thus, AR modeling was implemented to assess differences in vGRF pattern dynamics during symmetric and asymmetric walking. We hypothesized that the AR model coefficients would better detect differences amongst the symmetric and asymmetric walking conditions than the vGRF peak magnitude mean. Seventeen healthy individuals performed a protocol that involved walking on a split-belt instrumented treadmill at different symmetric (0.75m/s, 1.0 m/s, 1.5 m/s) and asymmetric (Side 1: 0.75m/s-Side 2:1.0 m/s; Side 1:1.0m/s-Side 2:1.5 m/s) gait conditions. Vertical ground reaction force peaks extracted during the weight-acceptance and propulsive phase of each step were used to construct a vGRF peak time series. Then, a second order AR model was fit to the vGRF peak waveform data to determine the AR model coefficients. The resulting AR coefficients were plotted on a stationarity triangle and their distance from the triangle centroid was computed. Significant differences in vGRF patterns were detected amongst the symmetric and asymmetric conditions using the AR modeling coefficients (p = 0.01); however, no differences were found when comparing vGRF peak magnitude means. These findings suggest that AR modeling has the sensitivity to identify differences in gait asymmetry that could aid in monitoring rehabilitation progression.
| S-EPMC7714243 | biostudies-literature