Search

Paired organs like eyes, ears, and lungs in humans exhibit similarities, and data from these organs often display remarkable correlations. Accounting for these correlations could enhance classification models used in predicting disease phenotypes. To our knowledge, there is limited, if any, literature addressing this topic, and existing methods do not exploit such correlations. For example, the conventional approach treats each ear as an independent observation when predicting audiometric phenotypes and is agnostic about the correlation of data from the two ears of the same person. This approach may lead to information loss and reduce the model performance. In response to this gap, particularly in the context of audiometric phenotype prediction, this paper proposes new quadratic discriminant analysis (QDA) algorithms that appropriately deal with the dependence between ears. We propose two‐stage analysis strategies: (1) conducting data transformations to reduce data dimensionality before applying QDA; and (2) developing new QDA algorithms to partially utilize the dependence between phenotypes of two ears. We conducted simulation studies to compare different transformation methods and to assess the performance of different QDA algorithms. The empirical results suggested that the transformation may only be beneficial when the sample size is relatively small. Moreover, our proposed new QDA algorithms performed better than the conventional approach in both person‐level and ear‐level accuracy. As an illustration, we applied them to audiometric data from the Medical University of South Carolina Longitudinal Cohort Study of Age‐related Hearing Loss. In addition, we developed an R package, PairQDA, to implement the proposed algorithms. [ABSTRACT FROM AUTHOR]

While advances continue to be made in model-based clustering, challenges persist in modeling various data types such as panel data. Multivariate panel data present difficulties for clustering algorithms because they are often plagued by missing data and dropouts, presenting issues for estimation algorithms. This research presents a family of hidden Markov models that compensate for the issues that arise in panel data. A modified expectation–maximization algorithm capable of handling missing not at random data and dropout is presented and used to perform model estimation. [ABSTRACT FROM AUTHOR]

Copyright of Laboratory Animals is the property of Sage Publications, Ltd. and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

Clustered data are a special type of correlated data where units within a cluster are correlated while units between different clusters are independent. The number of units in a cluster can be associated with that cluster's outcome. This is called the informative cluster size (ICS), which is known to impact clustered data inference. However, when comparing the outcomes from multiple groups of units in clustered data, investigating ICS may not be enough. This is because the number of units belonging to a particular group in a cluster can be associated with the outcome from that group in that cluster, leading to an informative intra-cluster group size or IICGS. This phenomenon of IICGS can exist even in the absence of ICS. Ignoring the existence of IICGS can result in a biased inference for group-based outcome comparisons in clustered data. In this article, we mathematically formulate the concept of IICGS while distinguishing it from ICS and propose a nonparametric bootstrap-based statistical hypothesis-testing mechanism for testing any claim of IICGS in a clustered data setting. Through simulations and real data applications, we demonstrate that our proposed statistical testing method can accurately identify IICGS, with substantial power, in clustered data. [ABSTRACT FROM AUTHOR]

When integrating data from multiple sources, a common challenge is block-wise missing. Most existing methods address this issue only in cross-sectional studies. In this paper, we propose a method for variable selection when combining datasets from multiple sources in longitudinal studies. To account for block-wise missing in covariates, we impute the missing values multiple times based on combinations of samples from different missing pattern and predictors from different data sources. We then use these imputed data to construct estimating equations, and aggregate the information across subjects and sources with the generalized method of moments. We employ the smoothly clipped absolute deviation penalty in variable selection and use the extended Bayesian Information Criterion criteria for tuning parameter selection. We establish the asymptotic properties of the proposed estimator, and demonstrate the superior performance of the proposed method through numerical experiments. Furthermore, we apply the proposed method in the Alzheimer's Disease Neuroimaging Initiative study to identify sensitive early-stage biomarkers of Alzheimer's Disease, which is crucial for early disease detection and personalized treatment. [ABSTRACT FROM AUTHOR]

Rapid development in data science enables machine learning and artificial intelligence to be the most popular research tools across various disciplines. While numerous articles have shown decent predictive ability, little research has examined the impact of complex correlated data. We aim to develop a more accurate model under repeated measures or hierarchical data structures. Therefore, this study proposes a novel algorithm, the Generalized Estimating Equations Boosting (GEEB) machine, to integrate the gradient boosting technique into the benchmark statistical approach that deals with the correlated data, the generalized Estimating Equations (GEE). Unlike the previous gradient boosting utilizing all input features, we randomly select some input features when building the model to reduce predictive errors. The simulation study evaluates the predictive performance of the GEEB, GEE, eXtreme Gradient Boosting (XGBoost), and Support Vector Machine (SVM) across several hierarchical structures with different sample sizes. Results suggest that the new strategy GEEB outperforms the GEE and demonstrates superior predictive accuracy than the SVM and XGBoost in most situations. An application to a real-world dataset, the Forest Fire Data, also revealed that the GEEB reduced mean squared errors by 4.5% to 25% compared to GEE, XGBoost, and SVM. This research also provides a freely available R function that could implement the GEEB machine effortlessly for longitudinal or hierarchical data. [ABSTRACT FROM AUTHOR]

In this paper we derive estimating equations for modeling unbalanced correlated data sets in which the marginal distributions follow the one parameter unit-Lindley distributions with domain on the interval (0,1). A class of regressions models is proposed for modeling the location parameter and a reweighted iterative process is developed for the joint estimation of the regression coefficients and the correlation structure. Simulation studies are performed to assess the empirical properties of the derived estimators and diagnostic procedures, such as residual analysis and sensitivity studies based on conformal local influence are given. Finally, we analyze the proportion of people in households with inadequate water supply and sewage within federation units of Brazil by the procedures developed in the paper. [ABSTRACT FROM AUTHOR]

Understanding the relationships between exposure and disease incidence is an important problem in environmental epidemiology. Typically, a large number of these exposures are measured, and it is found either that a few exposures transmit risk or that each exposure transmits a small amount of risk, but, taken together, these may pose a substantial disease risk. Further, these exposure effects can be nonlinear. We develop a latent functional approach, which assumes that the individual effect of each exposure can be characterized as one of a series of unobserved functions, where the number of latent functions is less than or equal to the number of exposures. We propose Bayesian methodology to fit models with a large number of exposures and show that existing Bayesian group LASSO approaches are a special case of the proposed model. An efficient Markov chain Monte Carlo sampling algorithm is developed for carrying out Bayesian inference. The deviance information criterion is used to choose an appropriate number of nonlinear latent functions. We demonstrate the good properties of the approach using simulation studies. Further, we show that complex exposure relationships can be represented with only a few latent functional curves. The proposed methodology is illustrated with an analysis of the effect of cumulative pesticide exposure on cancer risk in a large cohort of farmers. [ABSTRACT FROM AUTHOR]

Rheumatology research often involves correlated and clustered data. A common error when analyzing these data occurs when instead we treat these data as independent observations. This can lead to incorrect statistical inference. The data used are a subset of the 2017 study from Raheel et al consisting of 633 patients with rheumatoid arthritis (RA) between 1988 and 2007. RA flare and the number of swollen joints served as our binary and continuous outcomes, respectively. Generalized linear models (GLM) were fitted for each, while adjusting for rheumatoid factor (RF) positivity and sex. Additionally, a generalized linear mixed model with a random intercept and a generalized estimating equation were used to model RA flare and the number of swollen joints, respectively, to take additional correlation into account. The GLM's β coefficients and their 95% confidence intervals (CIs) are then compared to their mixed-effects equivalents. The β coefficients compared between methodologies are very similar. However, their standard errors increase when correlation is accounted for. As a result, if the additional correlations are not considered, the standard error can be underestimated. This results in an overestimated effect size, narrower CIs, increased type I error, and a smaller P value, thus potentially producing misleading results. It is important to model the additional correlation that occurs in correlated data. [ABSTRACT FROM AUTHOR]

Entropy is widely used in ecological and environmental studies, where data often present complex interactions. Difficulties arise in linking entropy to available covariates or data dependence structures, thus, all existing entropy estimators assume independence. To overcome this limit, we take a Bayesian model-based approach which focuses on estimating the probabilities that compose the index, accounting for any data dependence and correlation. An estimate of entropy can be constructed from the model fitted values, returning an observation-specific measure of entropy rather than an overall index. This way, the latent heterogeneity of the system can be represented by a curve in time or a surface in space, according to the characteristics of the survey study at hand. An empirical study illustrates the flexibility and interpretability of our results over temporally and spatially correlated data. An application is presented about the biodiversity of spatially structured rainforest tree data. [ABSTRACT FROM AUTHOR]

The area under the true ROC curve (AUC) is routinely used to determine how strongly a given model discriminates between the levels of a binary outcome. Standard inference with the AUC requires that outcomes be independent of each other. To overcome this limitation, a method was developed for the estimation of the variance of the AUC in the setting of two-level hierarchical data using probit-transformed prediction scores generated from generalized estimating equation models, thereby allowing for the application of inferential methods. This manuscript presents an extension of this approach so that inference for the AUC may be performed in a three-level hierarchical data setting (e.g., eyes nested within persons and persons nested within families). A method that accounts for the effect of tied prediction scores on inference is also described. The performance of 95% confidence intervals around the AUC was assessed through the simulation of three-level clustered data in multiple settings, including ones with tied data and variable cluster sizes. Across all settings, the actual 95% confidence interval coverage varied from 0.943 to 0.958, and the ratio of the theoretical variance to the empirical variance of the AUC varied from 0.920 to 1.013. The results are better than those from existing methods. Two examples of applying the proposed methodology are presented. [ABSTRACT FROM AUTHOR]

Background: Chemotherapy-induced peripheral neuropathy can occur in the right and left hand. Studies on prevention treatments for chemotherapy-induced peripheral neuropathy have largely adopted either self-controlled designs or parallel designs to compare two preventive treatments. When three treatment options (two experimental treatments and a control treatment) are available, both designs can be extended. However, no clinical trials have adopted a self-controlled design to compare three prevention treatments for chemotherapy-induced peripheral neuropathy. The incomplete block crossover design for more than two treatments can be extended to compare three treatments in the self-controlled design. In simple extension, some of the participants receive two experimental treatments in both hands; however, it may be difficult to administer different experimental treatments in both hands for practical reasons, such as a concern for the different types of unexpected adverse events. This study proposes a design and analysis method appropriate for the situation where only one experimental treatment is provided to each participant. Methods: We assume clinical trials to compare each of the two experimental treatments (E 1 and E 2) with the control treatment (C) and between two experimental treatments only when both experimental treatments are superior to the control treatment. We propose a self-controlled design, which equally randomizes to four arms to adjust for the dominant hand effect: Arm 1: E 1 for right hand, C for left hand; Arm 2: C for right hand, E 1 for left hand; Arm 3: E 2 for right hand, C for left hand; and Arm 4: C for right hand, E 2 for left hand. We compare operating characteristics of the proposed design with the three-arm parallel design in which the same treatment is performed in both hands by participants. We also assess three proposed analysis methods for comparisons between experimental treatments in the self-controlled design under several conditions of correlations between right and left hands using simulation studies. Results: The simulation studies showed that the proposed design was more powerful than the three-arm parallel design when correlation was 0.3 or higher. For comparisons between experimental treatments, the methods based on the regression model, including the outcome of hands with C as a covariate, had the highest power under modest to high correlation among the analysis methods in the self-controlled design. Conclusion: The proposed design can improve the power for comparing between two experimental treatments and the control treatment. Our design is useful in situations where it is undesirable for participants to receive different experimental treatments in both hands for practical reasons. [ABSTRACT FROM AUTHOR]

In recent decades, the use of regression models with random effects has made great progress. Among these models' attractions is the flexibility to analyze correlated data. In various situations, the distribution of the response variable presents asymmetry or bimodality. In these cases, it is possible to use the normal regression with random effect at the intercept. In light of these contexts, i.e. the desire to analyze correlated data in the presence of bimodality or asymmetry, in this paper we propose a regression model with random effect at the intercept based onthe generalized inverse Gaussian distribution model with correlated data. The maximum likelihood is adopted to estimate the parameters and various simulations are performed for correlated data. A type of residuals for the new regression is proposed whose empirical distribution is close to normal. The versatility of the new regression is demonstrated by estimating the average price per hectare of bare land in 10 municipalities in the state of São Paulo (Brazil). In this context, various databases are constantly emerging, requiring flexible modeling. Thus, it is likely to be of interest to data analysts, and can make a good contribution to the statistical literature. [ABSTRACT FROM AUTHOR]

We develop and calibrate stochastic continuous models that capture crime dynamics in the city of Valencia, Spain. From the emergency phone, data corresponding to three crime events, aggressions, stealing and women alarms, are available from the year 2010 until 2020. As the resulting time series, with monthly counts, are highly noisy, we decompose them into trend and seasonality parts. The former is modeled by geometric Brownian motions, both uncorrelated and correlated, and the latter is accommodated by randomly perturbed sine-cosine waves. Albeit simple, the models exhibit high ability to simulate the real data and show promising for crimes-interaction identification and short-term predictive policing. [ABSTRACT FROM AUTHOR]

Linear mixed effects models (LMMs) are a popular and powerful tool for analysing grouped or repeated observations for numeric outcomes. LMMs consist of a fixed and a random component, which are specified in the model through their respective design matrices. Verifying the correct specification of the two design matrices is important since mis-specifying them can affect the validity and efficiency of the analysis. We show how to use empirical stochastic processes constructed from appropriately ordered and standardized residuals from the model to test whether the design matrices of the fitted LMM are correctly specified. We define two different processes: one can be used to test whether both design matrices are correctly specified, and the other can be used only to test whether the fixed effects design matrix is correctly specified. The proposed empirical stochastic processes are smoothed versions of cumulative sum processes, which have a nice graphical representation in which model mis-specification can easily be observed. The amount of smoothing can be adjusted, which facilitates visual inspection and can potentially increase the power of the tests. We propose a computationally efficient procedure for estimating p-values in which refitting of the LMM is not necessary. Its validity is shown by using theoretical results and a large Monte Carlo simulation study. The proposed methodology could be used with LMMs with multilevel or crossed random effects. [ABSTRACT FROM AUTHOR]

Despite major advances in research and treatment, identifying important genotype risk factors for high blood pressure remains challenging. Traditional genome-wide association studies (GWAS) focus on one single nucleotide polymorphism (SNP) at a time. We aim to select among over half a million SNPs along with time-varying phenotype variables via simultaneous modeling and variable selection, focusing on the most dangerous blood pressure levels at high quantiles. Taking advantage of rich data from a large-scale public health study, we develop and apply a novel quantile penalized generalized estimating equations (GEE) approach, incorporating several key aspects including ultra-high dimensional genetic SNPs, the longitudinal nature of blood pressure measurements, time-varying covariates, and conditional high quantiles of blood pressure. Importantly, we identify interesting new SNPs for high blood pressure. Besides, we find blood pressure levels are likely heterogeneous, where the important risk factors identified differ among quantiles. This comprehensive picture of conditional quantiles of blood pressure can allow more insights and targeted treatments. We provide an efficient computational algorithm and prove consistency, asymptotic normality, and the oracle property for the quantile penalized GEE estimators with ultra-high dimensional predictors. Moreover, we establish model-selection consistency for high-dimensional BIC. Simulation studies show the promise of the proposed approach. for this article are available online. [ABSTRACT FROM AUTHOR]

We propose a model selection criterion for correlated survival data when the cluster size is informative to the outcome. This approach, called Resampling Cluster Survival Information Criterion (RCSIC), uses the Cox proportional hazards model that is weighted with the inverse of the cluster size. The RCSIC based on the within‐cluster resampling idea takes into account the possible variability of the within‐cluster subsampling and the possible informativeness of cluster sizes. The RCSIC allows for easy execution for the within‐cluster resampling idea without a large number of resamples of the data. In contrast with the traditional model selection method in survival analysis, the RCSIC has an additional penalization for the within‐cluster subsampling variability. Our simulations show the satisfactory results where the RCSIC provides a more robust power for variable selection in terms of clustered survival analysis, regardless of whether informative cluster size exists or not. Applying the RCSIC method to a periodontal disease studies, we identify the tooth loss in patients associated with the risk factors, Age, Filled Tooth, Molar, Crown, Decayed Tooth, and Smoking Status, respectively. [ABSTRACT FROM AUTHOR]

We introduce Gaussian orthogonal latent factor processes for modeling and predicting large correlated data. To handle the computational challenge, we first decompose the likelihood function of the Gaussian random field with a multi-dimensional input domain into a product of densities at the orthogonal components with lower-dimensional inputs. The continuous-time Kalman filter is implemented to compute the likelihood function efficiently without making approximations. We also show that the posterior distribution of the factor processes is independent, as a consequence of prior independence of factor processes and orthogonal factor loading matrix. For studies with large sample sizes, we propose a flexible way to model the mean, and we derive the marginal posterior distribution to solve identifiability issues in sampling these parameters. Both simulated and real data applications confirm the outstanding performance of this method. [ABSTRACT FROM AUTHOR]

In many situations, it is common to have more than one observation per experimental unit, thus generating the experiments with repeated measures. In the modeling of such experiments, it is necessary to consider and model the intra-unit dependency structure. In the literature, there are several proposals to model positive continuous data with repeated measures. In this paper, we propose one more with the generalization of the beta prime regression model. We consider the possibility of dependence between observations of the same unit. Residuals and diagnostic tools also are discussed. To evaluate the finite-sample performance of the estimators, using different correlation matrices and distributions, we conducted a Monte Carlo simulation study. The methodology proposed is illustrated with an analysis of a real data set. Finally, we create an R package for easy access to publicly available the methodology described in this paper. [ABSTRACT FROM AUTHOR]

Background: Studies that utilize RNA Sequencing (RNA-Seq) in conjunction with designs that introduce dependence between observations (e.g. longitudinal sampling) require specialized analysis tools to accommodate this additional complexity. This R package contains a set of utilities to fit linear mixed effects models to transformed RNA-Seq counts that properly account for this dependence when performing statistical analyses. Results: In a simulation study comparing lmerSeq and two existing methodologies that also work with transformed RNA-Seq counts, we found that lmerSeq was comprehensively better in terms of nominal error rate control and statistical power. Conclusions: Existing R packages for analyzing transformed RNA-Seq data with linear mixed models are limited in the variance structures they allow and/or the transformation methods they support. The lmerSeq package offers more flexibility in both of these areas and gave substantially better results in our simulations. [ABSTRACT FROM AUTHOR]

Many statistical models currently used in ecology and evolution account for covariances among random errors. Here, I address five points: (i) correlated random errors unite many types of statistical models, including spatial, phylogenetic and time‐series models; (ii) random errors are neither unpredictable nor mistakes; (iii) diagnostics for correlated random errors are not useful, but simulations are; (iv) model predictions can be made with random errors; and (v) can random errors be causal?These five points are illustrated by applying statistical models to analyse simulated spatial, phylogenetic and time‐series data. These three simulation studies are paired with three types of predictions that can be made using information from covariances among random errors: predictions for goodness‐of‐fit, interpolation, and forecasting.In the simulation studies, models incorporating covariances among random errors improve inference about the relationship between dependent and independent variables. They also imply the existence of unmeasured variables that generate the covariances among random errors. Understanding the covariances among random errors gives information about possible processes underlying the data.Random errors are caused by something. Therefore, to extract full information from data, covariances among random errors should not just be included in statistical models; they should also be studied in their own right. Data are hard won, and appropriate statistical analyses can make the most of them. [ABSTRACT FROM AUTHOR]

Performing a genome-wide association study (GWAS) with a binary phenotype using family data is a challenging task. Using linear mixed effects models is typically unsuitable for binary traits, and numerical approximations of the likelihood function may not work well with rare genetic variants with small counts. Additionally, imbalance in the case-control ratios poses challenges as traditional statistical methods such as the Score test or Wald test perform poorly in this setting. In the last couple of years, several methods have been proposed to better approximate the likelihood function of a mixed effects logistic regression model that uses Saddle Point Approximation (SPA). SPA adjustment has recently been implemented in multiple software, including GENESIS, SAIGE, REGENIE and fastGWA-GLMM: four increasingly popular tools to perform GWAS of binary traits. We compare Score and SPA tests using real family data to evaluate computational efficiency and the agreement of the results. Additionally, we compare various ways to adjust for family relatedness, such as sparse and full genetic relationship matrices (GRM) and polygenic effect estimates. We use the New England Centenarian Study imputed genotype data and the Long Life Family Study whole-genome sequencing data and the binary phenotype of human extreme longevity to compare the agreement of the results and tools' computational performance. The evaluation suggests that REGENIE might not be a good choice when analyzing correlated data of a small size. fastGWA-GLMM is the most computationally efficient compared to the other three tools, but it appears to be overly conservative when applied to family-based data. GENESIS, SAIGE and fastGWA-GLMM produced similar, although not identical, results, with SPA adjustment performing better than Score tests. Our evaluation also demonstrates the importance of adjusting by full GRM in highly correlated datasets when using GENESIS or SAIGE. [ABSTRACT FROM AUTHOR]

High throughput technologies like microarrays, next generation sequencing and mass spectrometry enable the measurement of tens of thousands of biomarker candidates in pilot studies. Biological systems are often too complex to be based on simple single cause-effect associations and from the medical practice point of view, a single biomarker may not possess the desired sensitivity and/or specificity for disease classification and outcome prediction. Therefore, the efforts of researchers currently aims at combining biomarkers. The intention of biomarker pilot studies with small sample sizes is often to explore the possibility of finding good biomarker combinations and not to find and evaluate a final combination of biomarkers with high predictive value. The aim of the pilot study is to answer the question whether it is worthwhile to extend the study to a larger study and to obtain information about the required sample size. In this paper, we propose a method to judge the potential in a small biomarker pilot study without the need to explicitly identifying and confirming a specific subset of biomarkers. In addition, we provide a method for sample size estimation for an extended study when the results of the pilot study look promising. Abbreviations: ROC: receiver operating characteristic curve; AUC: Area Under the ROC curve; HAUCA: high AUC abundance; ER: Estrogen receptor; BCs: Biomarker candidates; w: with; wt: without [ABSTRACT FROM AUTHOR]

Summary: Differential privacy is a rigorous mathematical framework for evaluating and protecting data privacy. In most existing studies, there is a vulnerable assumption that records in a dataset are independent when differential privacy is applied. However, in real‐world datasets, records are likely to be correlated, which may lead to unexpected data leakage. In this survey, we investigate the issue of privacy loss due to data correlation under differential privacy models. Roughly, we classify existing literature into three lines: (1) using parameters to describe data correlation in differential privacy, (2) using models to describe data correlation in differential privacy, and (3) describing data correlation based on the framework of Pufferfish. First, a detailed example is given to illustrate the issue of privacy leakage on correlated data in real scenes. Then our main work is to analyze and compare these methods, and evaluate situations that these diverse studies are applied. Finally, we propose some future challenges on correlated differential privacy. [ABSTRACT FROM AUTHOR]

In ophthalmology, data from both eyes of a person are frequently included in the statistical evaluation. This violates the requirement of data independence for classical statistical tests (e.g. t-Test or analysis of variance (ANOVA)) because it is correlated data. Linear mixed models (LMM) were used as a possibility to include the data of both eyes in the statistical evaluation. The LMM is available for a variety of statistical software such as SPSS or R. The application was applied to a retrospective longitudinal analysis of an accelerated corneal cross-linking (ACXL (9*10)) treatment in progressive keratoconus (KC) with a follow-up period of 36 months. Forty eyes of 20 patients were included, whereas sequential bilateral CXL treatment was performed within 12 months. LMM and ANOVA for repeated measurements were used for statistical evaluation of topographical and tomographical data measured by Pentacam (Oculus, Wetzlar, Germany). Both eyes were classified into a worse and better eye concerning corneal topography. Visual acuity, keratometric values and minimal corneal thickness were statistically significant between them at baseline (p < 0.05). A significant correlation between worse and better eye was shown (p < 0.05). Therefore, analyzing the data at each follow-up visit using ANOVA partially led to an overestimation of the statistical effect that could be avoided by using LMM. After 36 months, ACXL has significantly improved BCVA and flattened the cornea. The evaluation of data of both eyes without considering their correlation using classical statistical tests leads to an overestimation of the statistical effect, which can be avoided by using the LMM. [ABSTRACT FROM AUTHOR]

Background: As the cost of RNA-sequencing decreases, complex study designs, including paired, longitudinal, and other correlated designs, become increasingly feasible. These studies often include multiple hypotheses and thus multiple degree of freedom tests, or tests that evaluate multiple hypotheses jointly, are often useful for filtering the gene list to a set of interesting features for further exploration while controlling the false discovery rate. Though there are several methods which have been proposed for analyzing correlated RNA-sequencing data, there has been little research evaluating and comparing the performance of multiple degree of freedom tests across methods.Methods: We evaluated 11 different methods for modelling correlated RNA-sequencing data by performing a simulation study to compare the false discovery rate, power, and model convergence rate across several hypothesis tests and sample size scenarios. We also applied each method to a real longitudinal RNA-sequencing dataset.Results: Linear mixed modelling using transformed data had the best false discovery rate control while maintaining relatively high power. However, this method had high model non-convergence, particularly at small sample sizes. No method had high power at the lowest sample size. We found a mix of conservative and anti-conservative behavior across the other methods, which was influenced by the sample size and the hypothesis being evaluated. The patterns observed in the simulation study were largely replicated in the analysis of a longitudinal study including data from intensive care unit patients experiencing cardiogenic or septic shock.Conclusions: Multiple degree of freedom testing is a valuable tool in longitudinal and other correlated RNA-sequencing experiments. Of the methods that we investigated, linear mixed modelling had the best overall combination of power and false discovery rate control. Other methods may also be appropriate in some scenarios. [ABSTRACT FROM AUTHOR]

This article demonstrates how data from a biology paper, which analyzes the relationship between mass and metabolic rate for two species of marine bryozoan, can be used to teach a variety of regression topics to both introductory and advanced students. A thorough analysis requires intelligent data wrangling, variable transformations, and accounting for correlation between observations. The bryozoan data can be used as a valuable class example throughout the semester, or as a dataset for extended homework assignments and class projects. for this article are available online. [ABSTRACT FROM AUTHOR]

Longitudinal assessments are crucial in evaluating the disease state and trajectory in patients with neurodegenerative diseases. Neuropsychological outcomes measured over time often have a non-linear trajectory with autocorrelated residuals and a skewed distribution. We propose the adjusted local linear trend model, an extended state-space model in lieu of the commonly used linear mixed-effects model in modeling longitudinal neuropsychological outcomes. Our contributed model has the capability to utilize information from the stochasticity of the data while accounting for subject-specific trajectories with the inclusion of covariates and unequally spaced time intervals. The first step of model fitting involves a likelihood maximization step to estimate the unknown variances in the model before parsing these values into the Kalman filter and Kalman smoother recursive algorithms. Results from simulation studies showed that the adjusted local linear trend model is able to attain lower bias, lower standard errors, and high power, particularly in short longitudinal studies with equally spaced time intervals, as compared to the linear mixed-effects model. The adjusted local linear trend model also outperforms the linear mixed-effects model when data is missing completely at random, missing at random, and, in certain cases, even in data with missing not at random. [ABSTRACT FROM AUTHOR]

Collinearity can be a problem in regression models. When examining the effects of an exposure at different time points, constrained distributed lag models can alleviate some of the problems caused by collinearity. Still, some consequences of collinearity may remain and they are often unexplored. We aimed to illustrate the effects of collinearity in the context of distributed lag models, and to provide a tool to assess whether the results of a study could be influenced by collinearity. We used simulations under different scenarios of hypothesized effects of an exposure to visualize the resulting curves of lagged effects. We analysed three real datasets: a cohort study looking for windows of vulnerability to air pollution, a time series study examining the linear association of air pollution with hospital admissions, and a time series study examining the non-linear association between temperature and mortality. We showed that collinearity could be the explanation for some unexpected results, e.g. for statistically significant associations in the opposite direction from that expected, or for wrongly suggesting that some periods are more important than others. We implemented the collin R package to explore the potential consequences of collinearity in the context of distributed lag models. Our visual tool can be a useful way to assess if the results of an analysis may be influenced by collinearity. [ABSTRACT FROM AUTHOR]

In studies following selective sampling protocols for secondary outcomes, conventional analyses regarding their appearance could provide misguided information. In the large type 1 diabetes prevention and prediction (DIPP) cohort study monitoring type 1 diabetes-associated autoantibodies, we propose to model their appearance via a multivariate frailty model, which incorporates a correlation component that is important for unbiased estimation of the baseline hazards under the selective sampling mechanism. As further advantages, the frailty model allows for systematic evaluation of the association and the differences in regression parameters among the autoantibodies. We demonstrate the properties of the model by a simulation study and the analysis of the autoantibodies and their association with background factors in the DIPP study, in which we found that high genetic risk is associated with the appearance of all the autoantibodies, whereas the association with sex and urban municipality was evident for IA-2A and IAA autoantibodies. [ABSTRACT FROM AUTHOR]

Clinical visit data are clustered within people, which complicates prediction modeling. Cluster size is often informative because people receiving more care are less healthy and at higher risk of poor outcomes. We used data from seven health systems on 1,518,968 outpatient mental health visits from January 1, 2012 to June 30, 2015 to predict suicide attempt within 90 days. We evaluated true performance of prediction models using a prospective validation set of 4,286,495 visits from October 1, 2015 to September 30, 2017. We examined dividing clustered data on the person or visit level for model training and cross‐validation and considered a within cluster resampling approach for model estimation. We evaluated optimism by comparing estimated performance from a left‐out testing dataset to performance in the prospective dataset. We used two prediction methods, logistic regression with least absolute shrinkage and selection operator (LASSO) and random forest. The random forest model using a visit‐level split for model training and testing was optimistic; it overestimated discrimination (area under the curve, AUC = 0.95 in testing versus 0.84 in prospective validation) and classification accuracy (sensitivity = 0.48 in testing versus 0.19 in prospective validation, 95th percentile cut‐off). Logistic regression and random forest models using a person‐level split performed well, accurately estimating prospective discrimination and classification: estimated AUCs ranged from 0.85 to 0.87 in testing versus 0.85 in prospective validation, and sensitivity ranged from 0.15 to 0.20 in testing versus 0.17 to 0.19 in prospective validation. Within cluster resampling did not improve performance. We recommend dividing clustered data on the person level, rather than visit level, to ensure strong performance in prospective use and accurate estimation of future performance at the time of model development. [ABSTRACT FROM AUTHOR]

In a cluster randomized trial (CRT), groups of people are randomly assigned to different interventions. Existing parametric and semiparametric methods for CRTs rely on distributional assumptions or a large number of clusters to maintain nominal confidence interval (CI) coverage. Randomization-based inference is an alternative approach that is distribution-free and does not require a large number of clusters to be valid. Although it is well-known that a CI can be obtained by inverting a randomization test, this requires testing a non-zero null hypothesis, which is challenging with non-continuous and survival outcomes. In this article, we propose a general method for randomization-based CIs using individual-level data from a CRT. This approach accommodates various outcome types, can account for design features such as matching or stratification, and employs a computationally efficient algorithm. We evaluate this method's performance through simulations and apply it to the Botswana Combination Prevention Project, a large HIV prevention trial with an interval-censored time-to-event outcome. [ABSTRACT FROM AUTHOR]

Before performing the inversion process, the original measured data set is often transformed (corrected, smoothed, Fourier-transformed, interpolated etc.). These preliminary transformations may make the original (statistically independent) noisy measurement data correlated. The noise correlation on transformed data must be taken into account in the parameter fitting procedure (inversion) by proper derivation of likelihood function. The covariance matrix of transformed data system is no longer diagonal, so the likelihood based metrics, which determines the fitting process is also changed as well as the results of inversion. In the practice, these changes are often neglected using the "customary" estimation procedure (simple least square method) resulting wrong uncertainty estimation and sometimes biased results. In this article the consequence of neglected correlation is studied and discussed by decomposing the inversion functional to "customary" and additional part which represents the effect of correlation. The ratio of two components demonstrates the importance and justification of the inversion method modification. [ABSTRACT FROM AUTHOR]

New technologies have produced increasingly complex and massive datasets, such as next generation sequencing and microarray data in biology, dynamic treatment regimes in clinical trials and long-term wide-scale studies in the social sciences. Each study exhibits its unique data structure within individuals, clusters and possibly across time and space. In order to draw valid conclusion from such large dimensional data, we must account for intracluster correlations, varying cluster sizes, and outliers in response and/or covariate domains to achieve valid and efficient inferences. A weighted rank-based method is proposed for selecting variables and estimating parameters simultaneously. The main contribution of the proposed method is four fold: (1) variable selection using adaptive lasso is extended to robust rank regression so that protection against outliers in both response and predictor variables is obtained; (2) within-subject correlations are incorporated so that efficiency of parameter estimation is improved; (3) the computation is convenient via the existing function in statistical software R. (4) the proposed method is proved to have desirable asymptotic properties for fixed number of covariates (p). Simulation studies are carried out to evaluate the proposed method for a number of scenarios including the cases when p equals to the number of subjects. The simulation results indicate that the proposed method is efficient and robust. A hormone dataset is analyzed for illustration. By adding additional redundant variables as covariates, the penalty approach and weighting schemes are proven to be effective. [ABSTRACT FROM AUTHOR]

An appropriate technique for monitoring a stochastic system is to utilize the control charts. In statistical process, control a drift in characteristics of output may be due to one or several assignable causes. In the establishment of X ¯ charts in statistical process control, an assumption is made that there is no correlation within the samples. However, in practice there are many cases where the correlation does exist within the samples. It would be more appropriate to assume that each sample is a realization of a multivariate normal random vector. Although some research works have been done on the economic design of control charts with single assignable cause with correlated data, the economic statistical design of X ¯ control chart under Weibull shock model with multiple assignable causes and correlated samples have not been presented yet. Based on the optimization of the average cost per unit of time and taking into account the different combination values of Weibull distribution parameters, optimal design values of sample size, sampling interval, and control limit coefficient were derived and calculated. Then the cost models under non-uniform and uniform sampling scheme were compared. The results revealed that the model under multiple assignable causes with correlated samples with non-uniform sampling has a lower cost than that with uniform sampling. [ABSTRACT FROM AUTHOR]

Continuous publication of statistics collected from various location-based applications may compromise users' privacy as the statistics could be procured from users' private data. Differential Privacy (DP) is a new privacy notion that offers a strong privacy guarantee to all users who participate in the statistics. However, the existing DP mechanism for continuous publication of location statistics provides a privacy guarantee with the assumption that the data-points of users stream at consecutive timestamps are independent. In reality, users' data-points may be temporally correlated, resulting in more privacy leakage due to an inadequate supply of privacy budget to the timestamps where the data-points are correlated. In this paper, we present a reformulated differential privacy definition to quantify the impact of temporal correlation on privacy leakage. Then, we introduce a privacy budget allocation method for allocating an adequate amount of privacy budget to each successive timestamps under the protection of differential privacy. Our solution adopts w-event privacy for continuously releasing statistics over infinite streams. The main idea is to check the dissimilarity between statistics at each timestamp and decide whether to publish current statistics or last release statistics. Finally, we evaluate the data utility of our proposed method by presenting experimental results for real and synthetic data sets. [ABSTRACT FROM AUTHOR]

We propose a multivariate regression model to deal with multiple continuous bounded data. The proposed model is based on second-moment assumptions, only. We adopted the quasi-score and Pearson estimating functions for estimation of the regression and dispersion parameters, respectively. Thus, the proposed approach does not require a multivariate probability distribution for the variable response vector. The multivariate quasi-beta regression model can easily handle multiple continuous bounded outcomes taking into account the correlation between the response variables. Furthermore, the model allows us to analyze continuous bounded data on the interval [0, 1], including zeros and/or ones. Simulation studies were conducted to investigate the behavior of the NORmal To Anything (NORTA) algorithm and to check the properties of the estimating function estimators to deal with multiple correlated response variables generated from marginal beta distributions. The model was motivated by a data set concerning the body fat percentage, which was measured at five regions of the body and represent the response variables. We analyze each response variable separately and compare it with the fit of the multivariate proposed model. The multivariate quasi-beta regression model provides better fit than its univariate counterparts, as well as allows us to measure the correlation between response variables. Finally, we adapted diagnostic tools to the proposed model. In the supplementary material, we provide the data set and R code. [ABSTRACT FROM AUTHOR]

Count data have emerged in many applied research areas. In recent years, there has been a considerable interest in models for count data. In modelling such data, it is common to face a large frequency of zeroes. The data are regarded as zero-inflated when the frequency of observed zeroes is larger than what is expected from a theoretical distribution such as Poisson distribution, as a standard model for analysing count data. Data analysis, using the simple Poisson model, may lead to over-dispersion. Several classes of different mixture models were proposed for handling zero-inflated data. But they do not apply to cases when inflated counts happen at some other points, in addition to zero. In these cases, a doubly-inflated Poisson model has been suggested which only be used for cross-sectional data and cannot consider correlations between observations. However, correlated count data have a large application, especially in the health and medical fields. The present study aims to introduce a Doubly-Inflated Poisson models with random effect for correlated doubly-inflated data. Then, the best performance of the proposed method is shown via different simulation scenarios. Finally, the proposed model is applied to a dental study. [ABSTRACT FROM AUTHOR]

Copyright of Afrika Statistika is the property of Statistics & Probability African Society and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

We present an online expectation maximization (EM) based algorithm for the problem of model parameter estimation and hypothesis testing based on observations from a network of heterogeneous sensors. The sensors’ measurements are assumed to be correlated and copula theory is used to model this correlation. Moreover, it is assumed that the statistical model for the sensors’ data is not completely known. We first develop the batch-mode EM to estimate the parameters of the underlying distributions of the sensors’ measurements and to detect the hypothesis. Next, we develop an online EM-based method that processes the data on a sample-by-sample basis. Numerical results are presented for both simulation data and two different real-world datasets. These results demonstrate the efficacy of the proposed method showing significant improvements for both parameter estimation and hypothesis testing compared to the method that ignores the correlation in sensors’ measurements. Hypothesis testing results from the real-world data are also compared with several unsupervised and supervised learning methods. It is shown that the proposed methods outperform the unsupervised and even some supervised methods. [ABSTRACT FROM AUTHOR]

The Loess Plateau is the largest loess accumulation area around the world, onto which loess landslides have been occurring frequently each year, thus bringing significant threats to communities there. To mitigate or manage the risk brought by loess landslides, many methods have been developed to gain insights into mechanisms that trigger loess landslides, or to identify regions that are susceptible to landslides through landslide susceptibility mapping. However, none of these methods can be used to quantitatively evaluate possible impact of potentially unstable slopes, which offers important information for risk management, especially for regions with high susceptibility to landslides. This study aims to fill this gap by constructing a loess landslide database first from field investigation. Then, a multivariate model for loess landslide data, including its height, width, area, and length, is developed considering correlation among these parameters. Subsequently, the multivariate model is used to predict statistically and quantitatively impact of a potentially unstable loess slope, in terms of slide width, length, and area, given height of the potentially unstable loess slope. The proposed method is applied to loess landslides occurred in Baoji City for illustration. Results show that the proposed method works reasonably well. In addition, some key equations are provided using results from the multivariate model. With these equations, geotechnical engineers or decision-makers can evaluate possible impact of a potentially unstable loess slope with minimal effort. [ABSTRACT FROM AUTHOR]

Purpose: To describe and demonstrate methods for analyzing longitudinal correlated eye data with a continuous outcome measure. Methods: We described fixed effects, mixed effects and generalized estimating equations (GEE) models, applied them to data from the Complications of Age-Related Macular Degeneration Prevention Trial (CAPT) and the Age-Related Eye Disease Study (AREDS). In CAPT (N = 1052), we assessed the effect of eye-specific laser treatment on change in visual acuity (VA). In the AREDS study, we evaluated effects of systemic supplement treatment among 1463 participants with AMD category 3. Results: In CAPT, the inter-eye correlations (0.33 to 0.53) and longitudinal correlations (0.31 to 0.88) varied. There was a small treatment effect on VA change (approximately one letter) at 24 months for all three models (p =.009 to 0.02). Model fit was better with the mixed effects model than the fixed effects model (p <.001). In AREDS, there was no significant treatment effect in all models (p >.55). Current smokers had a significantly greater VA decline than non-current smokers in the fixed effects model (p =.04) and the mixed effects model with random intercept (p =.0003), but marginally significant in the mixed effects model with random intercept and slope (p =.08), and GEE models (p =.054 to 0.07). The model fit was better with the fixed effects model than the mixed effects model (p <.0001). Conclusion: Longitudinal models using the eye as the unit of analysis can be implemented using available statistical software to account for both inter-eye and longitudinal correlations. Goodness-of-fit statistics may guide the selection of the most appropriate model. [ABSTRACT FROM AUTHOR]

Generalised estimating equations with the sandwich standard-error estimator provide a promising method of analysis for stepped wedge cluster randomised trials. However, they have inflated type-one error when used with a small number of clusters, which is common for stepped wedge cluster randomised trials. We present a large simulation study of binary outcomes comparing bias-corrected standard errors from Fay and Graubard; Mancl and DeRouen; Kauermann and Carroll; Morel, Bokossa, and Neerchal; and Mackinnon and White with an independent and exchangeable working correlation matrix. We constructed 95% confidence intervals using a t -distribution with degrees of freedom including clusters minus parameters (DFC-P), cluster periods minus parameters, and estimators from Fay and Graubard (DFFG), and Pan and Wall. Fay and Graubard and an approximation to Kauermann and Carroll (with simpler matrix inversion) were unbiased in a wide range of scenarios with an independent working correlation matrix and more than 12 clusters. They gave confidence intervals with close to 95% coverage with DFFG with 12 or more clusters, and DFC-P with 18 or more clusters. Both standard errors were conservative with fewer clusters. With an exchangeable working correlation matrix, approximated Kauermann and Carroll and Fay and Graubard had a small degree of under-coverage. [ABSTRACT FROM AUTHOR]

Although data analysis and mining technologies can efficiently provide intelligent and personalized services to us, data owners may not always be willing to share their true data because of privacy concerns. Recently, differential privacy (DP) technology has achieved a good trade-off between data utility and privacy guarantee by publishing noisy outputs. Nonetheless, DP still has a risk of privacy leakage when handling correlated data directly. Current schemes attempt to extend DP to publish correlated data, but are faced with the challenge of violating DP or low-level data utility. In this paper, we try to explore the essential cause of this inapplicability. Specifically, we suppose that this inapplicability is caused by the different correlations between noise and original data. To verify our supposition, we propose the notion of Correlation-Distinguishability Attack (CDA) to separate IID (Independent and Identically Distributed) noise from correlated data. Furthermore, taking time series as an example, we design an optimum filter to realize CDA in practical applications. Experimental results support our supposition and show that, the privacy degree of current approaches has a degradation under CDA. [ABSTRACT FROM AUTHOR]

RNA sequencing data have been abundantly generated in biomedical research for biomarker discovery and other studies. Such data at the exon level are usually heavily tailed and correlated. Conventional statistical tests based on the mean or median difference for differential expression likely suffer from low power when the between-group difference occurs mostly in the upper or lower tail of the distribution of gene expression. We propose a tail-based test to make comparisons between groups in terms of a specific distribution area rather than a single location. The proposed test, which is derived from quantile regression, adjusts for covariates and accounts for within-sample dependence among the exons through a specified correlation structure. Through Monte Carlo simulation studies, we show that the proposed test is generally more powerful and robust in detecting differential expression than commonly used tests based on the mean or a single quantile. An application to TCGA lung adenocarcinoma data demonstrates the promise of the proposed method in terms of biomarker discovery. [ABSTRACT FROM AUTHOR]

Background: There is high rate of under-five mortality in West Africa with little effort made to study determinants that significantly increase or decrease its risk across the West African sub-region. This is important since it will help in the design of effective intervention programs for each country or the entire region. The overall objective of this research evaluates the determinants of under-five mortality prior to the end of the 2015 Millennium Development Goals, to guide West African countries implement strategies that will aid them achieve the Sustainable Development Goal 3 by 2030. Method: This study used the Demographic and Health Survey (DHS) data from twelve (12) out of the eighteen West African countries; Ghana, Benin, Cote d' Ivoire, Guinea, Liberia, Mali, Niger, Nigeria, Sierra Leone, Burkina Faso, Gambia and Togo. Data were extracted from the children and women of reproductive age files as provided in the DHS report. The response or outcome variable of interest is under-five mortality rate. A Bayesian exponential, Weibull and Gompertz regression models via a gamma shared frailty model were used for the analysis. The deviance information criteria and Bayes factors were used to discriminate between models. These analyses were carried out using Stata version 15 software. Results: The study recorded 101 (95% CI: 98.6–103.5) deaths per 1000 live births occurring among the twelve countries. Burkina Faso (124.4), Cote D'lvoire (110.1), Guinea (116.4), Nigeria (120.6) and Niger (118.3) recorded the highest child under-5 mortality rate. Gambia (48.1), Ghana (60.1) and Benin (70.4) recorded the least unde-5 mortality rate per 1000 livebirths. Multiple birth children were about two times more likely to die compared to singleton birth, in all except Gambia, Nigeria and Sierra Leone. We observed significantly higher hazard rates for male compared to female children in the combined data analysis (HR: 1.14, 95% CI: [1.10–1.18]). The country specific analysis in Benin, Cote D'lvoire, Guinea, Liberia, Mali and Nigeria showed higher under-5 mortality hazard rates among male children compared to female children whilst Niger was the only country to report significantly lower hazard rate of males compared to females. Conclusion: There is still quite a substantial amount of work to be done in order to meet the Sustainable Development Goal 3 in 2030 in West Africa. There exist variant differences among some of the countries with respect to mortality rates and determinants which require different interventions and policy decisions. [ABSTRACT FROM AUTHOR]

Background: There is high rate of under-five mortality in West Africa with little effort made to study determinants that significantly increase or decrease its risk across the West African sub-region. This is important since it will help in the design of effective intervention programs for each country or the entire region. The overall objective of this research evaluates the determinants of under-five mortality prior to the end of the 2015 Millennium Development Goals, to guide West African countries implement strategies that will aid them achieve the Sustainable Development Goal 3 by 2030.Method: This study used the Demographic and Health Survey (DHS) data from twelve (12) out of the eighteen West African countries; Ghana, Benin, Cote d' Ivoire, Guinea, Liberia, Mali, Niger, Nigeria, Sierra Leone, Burkina Faso, Gambia and Togo. Data were extracted from the children and women of reproductive age files as provided in the DHS report. The response or outcome variable of interest is under-five mortality rate. A Bayesian exponential, Weibull and Gompertz regression models via a gamma shared frailty model were used for the analysis. The deviance information criteria and Bayes factors were used to discriminate between models. These analyses were carried out using Stata version 15 software.Results: The study recorded 101 (95% CI: 98.6-103.5) deaths per 1000 live births occurring among the twelve countries. Burkina Faso (124.4), Cote D'lvoire (110.1), Guinea (116.4), Nigeria (120.6) and Niger (118.3) recorded the highest child under-5 mortality rate. Gambia (48.1), Ghana (60.1) and Benin (70.4) recorded the least unde-5 mortality rate per 1000 livebirths. Multiple birth children were about two times more likely to die compared to singleton birth, in all except Gambia, Nigeria and Sierra Leone. We observed significantly higher hazard rates for male compared to female children in the combined data analysis (HR: 1.14, 95% CI: [1.10-1.18]). The country specific analysis in Benin, Cote D'lvoire, Guinea, Liberia, Mali and Nigeria showed higher under-5 mortality hazard rates among male children compared to female children whilst Niger was the only country to report significantly lower hazard rate of males compared to females.Conclusion: There is still quite a substantial amount of work to be done in order to meet the Sustainable Development Goal 3 in 2030 in West Africa. There exist variant differences among some of the countries with respect to mortality rates and determinants which require different interventions and policy decisions. [ABSTRACT FROM AUTHOR]

We propose a joint analysis of recurrent and nonrecurrent event data subject to general types of interval censoring. The proposed analysis allows for general semiparametric models, including the Box–Cox transformation and inverse Box–Cox transformation models for the recurrent and nonrecurrent events, respectively. A frailty variable is used to account for the potential dependence between the recurrent and nonrecurrent event processes, while leaving the distribution of the frailty unspecified. We apply the pseudolikelihood for interval‐censored recurrent event data, usually termed as panel count data, and the sufficient likelihood for interval‐censored nonrecurrent event data by conditioning on the sufficient statistic for the frailty and using the working assumption of independence over examination times. Large sample theory and a computation procedure for the proposed analysis are established. We illustrate the proposed methodology by a joint analysis of the numbers of occurrences of basal cell carcinoma over time and time to the first recurrence of squamous cell carcinoma based on a skin cancer dataset, as well as a joint analysis of the numbers of adverse events and time to premature withdrawal from study medication based on a scleroderma lung disease dataset. [ABSTRACT FROM AUTHOR]

Estimating equations for analyzing correlated Birnbaum–Saunders (BS) data are derived in this paper. A regression model is proposed for modeling the median of the life time until the failure, and a reweighted iterative process is developed for the joint estimation of the regression coefficients and the shape and correlation parameters. Diagnostic procedures, such as residual analysis and sensitivity studies based on case deletion and local influence, are given. Simulation studies are performed to assess the empirical distributions of the derived estimators and of a Pearson-type residual for correlated data. Finally, a longitudinal data set is analyzed by the procedures developed in the paper and extensions for the double case in which the median and the shape parameter are jointly modeled are discussed. [ABSTRACT FROM AUTHOR]

Background: A recently proposed method for estimating qPCR amplification efficiency E analyzes fluorescence intensity ratios from pairs of points deemed to lie in the exponential growth region on the amplification curves for all reactions in a dilution series. This method suffers from a serious problem: The resulting ratios are highly correlated, as they involve multiple use of the raw data, for example, yielding ~ 250 E estimates from ~ 25 intensity readings. The resulting statistics for such estimates are falsely optimistic in their assessment of the estimation precision. Results: Monte Carlo simulations confirm that the correlated pairs method yields precision estimates that are better than actual by a factor of two or more. This result is further supported by estimating E by both pairwise and Cq calibration methods for the 16 replicate datasets from the critiqued work, and then comparing the ensemble statistics for these methods. Conclusion: Contrary to assertions in the proposing work, the pairwise method does not yield E estimates a factor of 2 more precise than estimates from Cq calibration fitting (the standard curve method). On the other hand, the statistically correct direct fit of the data to the model behind the pairwise method can yield E estimates of comparable precision. Ways in which the approach might be improved are discussed briefly. [ABSTRACT FROM AUTHOR]