Nonparametric estimation of a Markov 'illness-death' process from interval-censored observations, with application to diabetes survival data

The nonparametric estimation of the cumulative transition intensityfunctions in a threestate time-nonhomogeneous Markov processwith irreversible transitions, an ‘illness-death’model, is considered when times of the intermediate transition,e.g. onset of a disease, are interval-censored. The times of‘death’ are assumed to be known exactly or to beright-censored. In addition the observed process may be left-truncated.Data of this type arise when the process is sampled periodically.For example, when the patients are monitored through periodicexaminations the observations on times of change in their diseasestatus will be interval-censored. Under the sampling schemeconsidered here the Nelson–Aalen estimator (Aalen, 1978)for a cumulative transition intensity is not applicable. Inthe proposed method the maximum likelihood estimators of someof the transition intensities are derived from the estimatorsof the corresponding subdistribution functions. The maximumlikelihood estimators are shown to have a self-consistency property.The self-consistency algorithm is developed for the computationof the estimators. This approach generalises the results fromTurnbull (1976) and Frydman (1992). The methods are illustratedwith diabetes survival data.     

Maximum likelihood estimation of ordered multinomial parameters

The pool adjacent violator algorithm Ayer et al. (1955, The Annals of Mathematical Statistics, 26, 641-647) has long been known to give the maximum likelihood estimator of a series of ordered binomial parameters, based on an independent observation from each distribution (see Barlow et al., 1972, Statistical Inference under Order Restrictions, Wiley, New York). This result has immediate application to estimation of a survival distribution based on current survival status at a set of monitoring times. This paper considers an extended problem of maximum likelihood estimation of a series of 'ordered' multinomial parameters p(i)= (p(1i),p(2i),.,p(mi)) for 1 相似文献   

Why we should take care of the competing risk bias in survival analysis: A phase II trial on the toxicity profile of radiotherapy for prostate cancer

AimThe aim of the present study is to evaluate and quantify the bias of competing risks in an Italian oncologic cohort comparing results from different statistical analysis methods.BackgroundCompeting risks are very common in randomized clinical trials and observational studies, in particular oncology and radiotherapy ones, and their inappropriate management causes results distortions widely present in clinical scientific articles.Materials and methodsThis is a single-institution phase II trial including 41 patients affected by prostate cancer and undergoing radiotherapy (IMRT-SIB) at the University Hospital of Udine.Different outcomes were considered: late toxicities, relapse, death.Death in the absence of relapse or late toxicity was considered as a competing event.ResultsThe Kaplan Meier method, compared to cumulative incidence function method, overestimated the probability of the event of interest (toxicity and biochemical relapse) and of the competing event (death without toxicity/relapse) by 9.36%. The log-rank test, compared to Gray's test, overestimated the probability of the event of interest by 5.26%.The Hazard Ratio's and cause specific hazard's Cox regression are not directly comparable to subdistribution hazard's Fine and Gray's modified Cox regression; nonetheless, the FG model, the best choice for prognostic studies with competing risks, found significant associations not emerging with Cox regression.ConclusionsThis study confirms that using inappropriate statistical methods produces a 10% overestimation in results, as described in the literature, and highlights the importance of taking into account the competing risks bias.     

Semiparametric models for cumulative incidence functions

In analyses of time-to-failure data with competing risks, cumulative incidence functions may be used to estimate the time-dependent cumulative probability of failure due to specific causes. These functions are commonly estimated using nonparametric methods, but in cases where events due to the cause of primary interest are infrequent relative to other modes of failure, nonparametric methods may result in rather imprecise estimates for the corresponding subdistribution. In such cases, it may be possible to model the cause-specific hazard of primary interest parametrically, while accounting for the other modes of failure using nonparametric estimators. The cumulative incidence estimators so obtained are simple to compute and are considerably more efficient than the usual nonparametric estimator, particularly with regard to interpolation of cumulative incidence at early or intermediate time points within the range of data used to fit the function. More surprisingly, they are often nearly as efficient as fully parametric estimators. We illustrate the utility of this approach in the analysis of patients treated for early stage breast cancer.     

The influence of competing-risks setting on the choice of hypothesis test for treatment effect

There is considerable debate regarding the choice of test for treatment difference in a randomized clinical trial in the presence of competing risks. This question arose in the study of standard and new antiepileptic drugs (SANAD) trial comparing new and standard antiepileptic drugs. This paper provides simulation results for the log-rank test comparing cause-specific hazard rates and Gray's test comparing cause-specific cumulative incidence curves. To inform the analysis of the SANAD trial, competing-risks settings were considered where both events are of interest, events may be negatively correlated, and the degree of correlation may differ in the 2 treatment groups. In settings where there are effects in opposite directions for the 2 event types, a likely situation for the SANAD trial, Gray's test has greater power to detect treatment differences than log-rank analysis. For the epilepsy application, conclusions were qualitatively similar for both log-rank and Gray's tests.