Longitudinal outcomes of high-risk human papillomavirus (HPV) infections as competing-risks events following cervical HPV test at baseline visit in the *NIS-LAMS** cohort

Background: The complex natural history of human papillomavirus (HPV) infections following a single HPV test can be modeled as competing-risks events (i.e., no-, transient- or persistent infection) in a longitudinal setting. The covariates associated with these compet ng events have not been previously assessed using competing-risks regression models. Objectives: To gain further insights in the outcomes of cervical HPV infections, we used univariate- and multivariate competing-risks regression models to assess the covariaies associated with these competing events. Study Design and Methods: Covariates associated with three competing outcomes (no-, transient- or persistent HR-HPV infection) were analysed in a sub-cohort of 1,865 women prospectively followed-up in the NIS (n = 3,187) and LAMS Study (n = 12,114). Results: In multivariate competing-risks models (with two other outcomes as competing events), permanently HR-HPV negative outcome was significantly predicted only by the clearance of ASCUS+Pap during FU, while three independent covariates predicted transient HR-HPV infections: i) number of recent (< 12 months) sexual partners (risk increased), ii) previous Pap screening history (protective), and history of previous CIN (increased risk). The two most powerful predictors of persistent HR-HPV infections were persistent ASCUS+Pap (risk increased), and previous Pap screening history (protective). In pair-wise comparisons, number of recent sexual partners and previous CIN history increase the probability of transient HR-HPV infection against the HR-HPV negative competing event, while previous Pap screening history is protective. Persistent ASCUS+Pap during FU and no previous Pap screening history are significantly associated with the persistent HR-HPV outcome (compared both with i) always negative, and ii) transient events), whereas multiparity is protective. Conclusions: Different covariates are associated with the three main outcomes of cervical HPV infections. The most significant covariates of each competing events are probably distinct enough to enable constructing of a risk-profile for each main outcome.

Covariates of high-risk human papillomavirus (HPV) infections are distinct for incident CIN1, CIN2 and CIN3 as disclosed by competing-risks regression models

Background: In addition to the oncogenic human papillomavirus (HPV), several cofactors are needed in cervical carcinogenesis, but whether the HPV covariates associated with incident i) CIN1 are different from those of incident ii) CIN2 and iii) CIN3 needs further assessment. Objectives: To gain further insights into the true biological differences between CIN1, CIN2 and CIN3, we assessed HPV covariates associated with incident CIN1, CIN2, and CIN3. Study Design and Methods: HPV covariates associated with progression to CIN1, CIN2 and CIN3 were analysed in the combined cohort of the NIS (n = 3,187) and LAMS study (n = 12,114), using competing-risks regression models (in panel data) for baseline HR-HPV-positive women (n = 1,105), who represent a sub-cohort of all 1,865 women prospectively followed-up in these two studies. Results: Altogether, 90 (4.8%), 39 (2.1%) and 14 (1.4%) cases progressed to CIN1, CIN2, and CIN3, respectively. Among these baseline HR-HPV-positive women, the risk profiles of incident GIN I, CIN2 and CIN3 were unique in that completely different HPV covariates were associated with progression to CIN1, CIN2 and CIN3, irrespective which categories (non-progression, CIN1, CIN2, CIN3 or all) were used as competing-risks events in univariate and multivariate models. Conclusions: These data confirm our previous analysis based on multinomial regression models implicating that distinct covariates of HR-HPV are associated with progression to CIN1, CIN2 and CIN3. This emphasises true biological differences between the three grades of GIN, which revisits the concept of combining CIN2 with CIN3 or with CIN1 in histological classification or used as a common end-point, e.g., in HPV vaccine trials.

An EM-based Semi-Parametric Mixture Model Approach to the Regression Analysis of Competing-Risks Data

We consider a mixture model approach to the regression analysis of competing-risks data. Attention is focused on inference concerning the effects of factors on both the probability of occurrence and the hazard rate conditional on each of the failure types. These two quantities are specified in the mixture model using the logistic model and the proportional hazards model, respectively. We propose a semi-parametric mixture method to estimate the logistic and regression coefficients jointly, whereby the component-baseline hazard functions are completely unspecified. Estimation is based on maximum likelihood on the basis of the full likelihood, implemented via an expectation-conditional maximization (ECM) algorithm. Simulation studies are performed to compare the performance of the proposed semi-parametric method with a fully parametric mixture approach. The results show that when the component-baseline hazard is monotonic increasing, the semi-parametric and fully parametric mixture approaches are comparable for mildly and moderately censored samples. When the component-baseline hazard is not monotonic increasing, the semi-parametric method consistently provides less biased estimates than a fully parametric approach and is comparable in efficiency in the estimation of the parameters for all levels of censoring. The methods are illustrated using a real data set of prostate cancer patients treated with different dosages of the drug diethylstilbestrol. Copyright (C) 2003 John Wiley Sons, Ltd.

On modifications to the long-term survival mixture model in the presence of competing risks

A mixture model for long-term survivors has been adopted in various fields such as biostatistics and criminology where some individuals may never experience the type of failure under study. It is directly applicable in situations where the only information available from follow-up on individuals who will never experience this type of failure is in the form of censored observations. In this paper, we consider a modification to the model so that it still applies in the case where during the follow-up period it becomes known that an individual will never experience failure from the cause of interest. Unless a model allows for this additional information, a consistent survival analysis will not be obtained. A partial maximum likelihood (ML) approach is proposed that preserves the simplicity of the long-term survival mixture model and provides consistent estimators of the quantities of interest. Some simulation experiments are performed to assess the efficiency of the partial ML approach relative to the full ML approach for survival in the presence of competing risks.

Health investment complementarities under competing risks

Applying the competing--risks model to multi--cause mortality,this paper provides a theoretical and empirical investigation of the positive complementarities that occur between disease--specific policy interventions. We argue that since an individual cannot die twice, competing risks imply that individuals will not waste resources on causes that are not the most immediate, but will make health investments so as to equalize cause--specific mortality. However, equal mortality risk from a variety of diseases does not imply that disease--specific public health interventions are a waste. Rather, a cause--specific intervention produces spillovers to other disease risks, so that the overall reduction in mortality will generally be larger than the direct effect measured on the targeted disease. The assumption that mortality from non--targeted diseases remains the same after a cause--specific intervention under--estimates the true effect of such programs, since the background mortality is also altered as a result of intervention. Analyzing data from one of the most important public health programs ever introduced, the Expanded Program on Immunization (EPI) of the United Nations, we find evidence for the existence of such complementarities, involving causes that are not biomedically, but behaviorally, linked.

Competing risks to breast cancer mortality in Catalonia

Background: Breast cancer mortality has experienced important changes over the last century. Breast cancer occurs in the presence of other competing risks which can influence breast cancer incidence and mortality trends. The aim of the present work is: 1) to assess the impact of breast cancer deaths among mortality from all causes in Catalonia (Spain), by age and birth cohort and 2) to estimate the risk of death from other causes than breast cancer, one of the inputs needed to model breast cancer mortality reduction due to screening or therapeutic interventions. Methods: The multi-decrement life table methodology was used. First, all-cause mortality probabilities were obtained by age and cohort. Then mortality probability for breast cancer was subtracted from the all-cause mortality probabilities to obtain cohort life tables for causes other than breast cancer. These life tables, on one hand, provide an estimate of the risk of dying from competing risks, and on the other hand, permit to assess the impact of breast cancer deaths on all-cause mortality using the ratio of the probability of death for causes other than breast cancer by the all-cause probability of death. Results: There was an increasing impact of breast cancer on mortality in the first part of the 20th century, with a peak for cohorts born in 1945–54 in the 40–49 age groups (for which approximately 24% of mortality was due to breast cancer). Even though for cohorts born after 1955 there was only information for women under 50, it is also important to note that the impact of breast cancer on all-cause mortality decreased for those cohorts. Conclusion: We have quantified the effect of removing breast cancer mortality in different age groups and birth cohorts. Our results are consistent with US findings. We also have obtained an estimate of the risk of dying from competing-causes mortality, which will be used in the assessment of the effect of mammography screening on breast cancer mortality in Catalonia.

The Conway-Maxwell-Poisson-generalized gamma regression model with long-term survivors

In this paper, we proposed a flexible cure rate survival model by assuming the number of competing causes of the event of interest following the Conway-Maxwell distribution and the time for the event to follow the generalized gamma distribution. This distribution can be used to model survival data when the hazard rate function is increasing, decreasing, bathtub and unimodal-shaped including some distributions commonly used in lifetime analysis as particular cases. Some appropriate matrices are derived in order to evaluate local influence on the estimates of the parameters by considering different perturbations, and some global influence measurements are also investigated. Finally, data set from the medical area is analysed.

The Geometric Birnbaum-Saunders regression model with cure rate

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

The geometric Birnbaum-Saunders regression model with cure rate

In this paper, we propose a cure rate survival model by assuming the number of competing causes of the event of interest follows the Geometric distribution and the time to event follow a Birnbaum Saunders distribution. We consider a frequentist analysis for parameter estimation of a Geometric Birnbaum Saunders model with cure rate. Finally, to analyze a data set from the medical area. (C) 2011 Elsevier B.V. All rights reserved.

Semiparametric Methods for Semi-competing Risks Problem with Censoring and Truncation

Studies of chronic life-threatening diseases often involve both mortality and morbidity. In observational studies, the data may also be subject to administrative left truncation and right censoring. Since mortality and morbidity may be correlated and mortality may censor morbidity, the Lynden-Bell estimator for left truncated and right censored data may be biased for estimating the marginal survival function of the non-terminal event. We propose a semiparametric estimator for this survival function based on a joint model for the two time-to-event variables, which utilizes the gamma frailty specification in the region of the observable data. Firstly, we develop a novel estimator for the gamma frailty parameter under left truncation. Using this estimator, we then derive a closed form estimator for the marginal distribution of the non-terminal event. The large sample properties of the estimators are established via asymptotic theory. The methodology performs well with moderate sample sizes, both in simulations and in an analysis of data from a diabetes registry.

A log-extended Weibull regression model

A bathtub-shaped failure rate function is very useful in survival analysis and reliability studies. The well-known lifetime distributions do not have this property. For the first time, we propose a location-scale regression model based on the logarithm of an extended Weibull distribution which has the ability to deal with bathtub-shaped failure rate functions. We use the method of maximum likelihood to estimate the model parameters and some inferential procedures are presented. We reanalyze a real data set under the new model and the log-modified Weibull regression model. We perform a model check based on martingale-type residuals and generated envelopes and the statistics AIC and BIC to select appropriate models. (C) 2009 Elsevier B.V. All rights reserved.

A Log-Linear Regression Model for the Beta-Weibull Distribution

We introduce the log-beta Weibull regression model based on the beta Weibull distribution (Famoye et al., 2005; Lee et al., 2007). We derive expansions for the moment generating function which do not depend on complicated functions. The new regression model represents a parametric family of models that includes as sub-models several widely known regression models that can be applied to censored survival data. We employ a frequentist analysis, a jackknife estimator, and a parametric bootstrap for the parameters of the proposed model. We derive the appropriate matrices for assessing local influences on the parameter estimates under different perturbation schemes and present some ways to assess global influences. Further, for different parameter settings, sample sizes, and censoring percentages, several simulations are performed. In addition, the empirical distribution of some modified residuals are displayed and compared with the standard normal distribution. These studies suggest that the residual analysis usually performed in normal linear regression models can be extended to a modified deviance residual in the proposed regression model applied to censored data. We define martingale and deviance residuals to evaluate the model assumptions. The extended regression model is very useful for the analysis of real data and could give more realistic fits than other special regression models.

Likelihood-based estimation of the regression model with scrambled responses

A significant problem in the collection of responses to potentially sensitive questions, such as relating to illegal, immoral or embarrassing activities, is non-sampling error due to refusal to respond or false responses. Eichhorn & Hayre (1983) suggested the use of scrambled responses to reduce this form of bias. This paper considers a linear regression model in which the dependent variable is unobserved but for which the sum or product with a scrambling random variable of known distribution, is known. The performance of two likelihood-based estimators is investigated, namely of a Bayesian estimator achieved through a Markov chain Monte Carlo (MCMC) sampling scheme, and a classical maximum-likelihood estimator. These two estimators and an estimator suggested by Singh, Joarder & King (1996) are compared. Monte Carlo results show that the Bayesian estimator outperforms the classical estimators in almost all cases, and the relative performance of the Bayesian estimator improves as the responses become more scrambled.

Bayesian Model Averaging in the Instrumental Variable Regression Model

This paper considers the instrumental variable regression model when there is uncertainty about the set of instruments, exogeneity restrictions, the validity of identifying restrictions and the set of exogenous regressors. This uncertainty can result in a huge number of models. To avoid statistical problems associated with standard model selection procedures, we develop a reversible jump Markov chain Monte Carlo algorithm that allows us to do Bayesian model averaging. The algorithm is very exible and can be easily adapted to analyze any of the di¤erent priors that have been proposed in the Bayesian instrumental variables literature. We show how to calculate the probability of any relevant restriction (e.g. the posterior probability that over-identifying restrictions hold) and discuss diagnostic checking using the posterior distribution of discrepancy vectors. We illustrate our methods in a returns-to-schooling application.

An alternative method for dating unknown tephras based on a segmented regression model

In CoDaWork’05, we presented an application of discriminant function analysis (DFA) to 4 differentcompositional datasets and modelled the first canonical variable using a segmented regression modelsolely based on an observation about the scatter plots. In this paper, multiple linear regressions areapplied to different datasets to confirm the validity of our proposed model. In addition to dating theunknown tephras by calibration as discussed previously, another method of mapping the unknown tephrasinto samples of the reference set or missing samples in between consecutive reference samples isproposed. The application of these methodologies is demonstrated with both simulated and real datasets.This new proposed methodology provides an alternative, more acceptable approach for geologists as theirfocus is on mapping the unknown tephra with relevant eruptive events rather than estimating the age ofunknown tephra.Kew words: Tephrochronology; Segmented regression