Rate/Mean Regression for Multiple-Sequence Recurrent Event Data with Missing Event Category

Abstract.  Censored recurrent event data frequently arise in biomedical studies. Often, the events are not homogenous, and may be categorized. We propose semiparametric regression methods for analysing multiple-category recurrent event data and consider the setting where event times are always known, but the information used to categorize events may be missing. Application of existing methods after censoring events of unknown category (i.e. 'complete-case' methods) produces consistent estimators only when event types are missing completely at random, an assumption which will frequently fail in practice. We propose methods, based on weighted estimating equations, which are applicable when event category missingness is missing at random. Parameter estimators are shown to be consistent and asymptotically normal. Finite sample properties are examined through simulations and the proposed methods are applied to an end-stage renal disease data set obtained from a national organ failure registry.     

Additive–multiplicative rates model for recurrent events

Recurrent events are frequently encountered in biomedical studies. Evaluating the covariates effects on the marginal recurrent event rate is of practical interest. There are mainly two types of rate models for the recurrent event data: the multiplicative rates model and the additive rates model. We consider a more flexible additive–multiplicative rates model for analysis of recurrent event data, wherein some covariate effects are additive while others are multiplicative. We formulate estimating equations for estimating the regression parameters. The estimators for these regression parameters are shown to be consistent and asymptotically normally distributed under appropriate regularity conditions. Moreover, the estimator of the baseline mean function is proposed and its large sample properties are investigated. We also conduct simulation studies to evaluate the finite sample behavior of the proposed estimators. A medical study of patients with cystic fibrosis suffered from recurrent pulmonary exacerbations is provided for illustration of the proposed method.     

An additive-multiplicative rates model for recurrent event data with informative terminal event

In this article, we propose an additive-multiplicative rates model for recurrent event data in the presence of a terminal event such as death. The association between recurrent and terminal events is nonparametric. For inference on the model parameters, estimating equation approaches are developed, and the asymptotic properties of the resulting estimators are established. The finite sample behavior of the proposed estimators is evaluated through simulation studies, and an application to a bladder cancer study is provided.     

A Marginal Additive Rates Model for Recurrent Event Data with a Terminal Event

Recurrent events data with a terminal event often arise in many longitudinal studies. Most of existing models assume multiplicative covariate effects and model the conditional recurrent event rate given survival. In this article, we propose a marginal additive rates model for recurrent events with a terminal event, and develop two procedures for estimating the model parameters. The asymptotic properties of the resulting estimators are established. In addition, some numerical procedures are presented for model checking. The finite-sample behavior of the proposed methods is examined through simulation studies, and an application to a bladder cancer study is also illustrated.     

Multiple imputation methods for recurrent event data with missing event category

Frequently in clinical and epidemiologic studies, the event of interest is recurrent (i.e., can occur more than once per subject). When the events are not of the same type, an analysis which accounts for the fact that events fall into different categories will often be more informative. Often, however, although event times may always be known, information through which events are categorized may potentially be missing. Complete‐case methods (whose application may require, for example, that events be censored when their category cannot be determined) are valid only when event categories are missing completely at random. This assumption is rather restrictive. The authors propose two multiple imputation methods for analyzing multiple‐category recurrent event data under the proportional means/rates model. The use of a proper or improper imputation technique distinguishes the two approaches. Both methods lead to consistent estimation of regression parameters even when the missingness of event categories depends on covariates. The authors derive the asymptotic properties of the estimators and examine their behaviour in finite samples through simulation. They illustrate their approach using data from an international study on dialysis.     

Additive Transformation Models for Recurrent Events

In this article, we propose a class of additive transformation models for recurrent event data, which includes the additive rates model as a special case. The new models offer great flexibility in formulating the effects of covariates on the mean function of recurrent events. Estimating equation approaches are developed for the model parameters, and asymptotic properties of the resulting estimators are established. In addition, a model checking procedure is presented to assess the adequacy of the model. The finite sample performance of the proposed estimators is examined through simulation studies, and an application to a bladder cancer study is presented.     

Semiparametric additive marginal regression models for multiple type recurrent events

Recurrent event data are often encountered in biomedical research, for example, recurrent infections or recurrent hospitalizations for patients after renal transplant. In many studies, there are more than one type of events of interest. Cai and Schaube (Lifetime Data Anal 10:121-138, 2004) advocated a proportional marginal rate model for multiple type recurrent event data. In this paper, we propose a general additive marginal rate regression model. Estimating equations approach is used to obtain the estimators of regression coefficients and baseline rate function. We prove the consistency and asymptotic normality of the proposed estimators. The finite sample properties of our estimators are demonstrated by simulations. The proposed methods are applied to the India renal transplant study to examine risk factors for bacterial, fungal and viral infections.     

Non parametric analysis of gap times for bivariate recurrent event data with cure fraction

Recurrent event data often arise in longitudinal studies. In many applications, subjects may experience two different types of events alternatively over time or a pair of subjects may experience recurrent events of the same type. Medical advances have made it possible for some patients to be cured such that the disease of interest does not recur. In this article, we consider non parametric analysis of bivariate recurrent event data with cure fraction. Using the inverse-probability weighted (IPW) approach, we propose non parametric estimators for the proportion of cured patients and for the joint distribution functions of bivariate recurrence times of the uncured ones. The asymptotic properties of the proposed estimators are established. Simulation study indicates that the proposed estimators perform well in finite samples.     

The analysis of multivariate recurrent events with partially missing event types

In many clinical studies, subjects are at risk of experiencing more than one type of potentially recurrent event. In some situations, however, the occurrence of an event is observed, but the specific type is not determined. We consider the analysis of this type of incomplete data when the objectives are to summarize features of conditional intensity functions and associated treatment effects, and to study the association between different types of event. Here we describe a likelihood approach based on joint models for the multi-type recurrent events where parameter estimation is obtained from a Monte-Carlo EM algorithm. Simulation studies show that the proposed method gives unbiased estimators for regression coefficients and variance–covariance parameters, and the coverage probabilities of confidence intervals for regression coefficients are close to the nominal level. When the distribution of the frailty variable is misspecified, the method still provides estimators of the regression coefficients with good properties. The proposed method is applied to a motivating data set from an asthma study in which exacerbations were to be sub-typed by cellular analysis of sputum samples as eosinophilic or non-eosinophilic.     

Nonparametric estimators for a survivor function of paired recurrent events

Recurrent event data arise in longitudinal studies where each study subject may experience multiple events during the follow-up. In many situations in survival studies, pairs of individuals can potentially experience recurrent events. The analysis of such data is not straightforward as it involves two kinds of dependences, namely, dependence between the individuals in the same pair and dependence among a sequence of pairs. In the present paper, we introduce a new stochastic model for the analysis of such recurrent event data. Nonparametric estimators for a bivariate survivor function are developed. Asymptotic properties of the estimators are discussed. Simulation studies are carried out to assess the finite sample properties of the estimator. We illustrate the procedure with real life data on eye disease.     

Simulated maximum likelihood estimation in joint models for multiple longitudinal markers and recurrent events of multiple types,in the presence of a terminal event

In medical studies we are often confronted with complex longitudinal data. During the follow-up period, which can be ended prematurely by a terminal event (e.g. death), a subject can experience recurrent events of multiple types. In addition, we collect repeated measurements from multiple markers. An adverse health status, represented by ‘bad’ marker values and an abnormal number of recurrent events, is often associated with the risk of experiencing the terminal event. In this situation, the missingness of the data is not at random and, to avoid bias, it is necessary to model all data simultaneously using a joint model. The correlations between the repeated observations of a marker or an event type within an individual are captured by normally distributed random effects. Because the joint likelihood contains an analytically intractable integral, Bayesian approaches or quadrature approximation techniques are necessary to evaluate the likelihood. However, when the number of recurrent event types and markers is large, the dimensionality of the integral is high and these methods are too computationally expensive. As an alternative, we propose a simulated maximum-likelihood approach based on quasi-Monte Carlo integration to evaluate the likelihood of joint models with multiple recurrent event types and markers.     

Joint model for bivariate zero-inflated recurrent event data with terminal events

Bivariate recurrent event data are observed when subjects are at risk of experiencing two different type of recurrent events. In this paper, our interest is to suggest statistical model when there is a substantial portion of subjects not experiencing recurrent events but having a terminal event. In a context of recurrent event data, zero events can be related with either the risk free group or a terminal event. For simultaneously reflecting both a zero inflation and a terminal event in a context of bivariate recurrent event data, a joint model is implemented with bivariate frailty effects. Simulation studies are performed to evaluate the suggested models. Infection data from AML (acute myeloid leukemia) patients are analyzed as an application.     

Marginal Regression of Multivariate Event Times Based on Linear Transformation Models

Multivariate event time data are common in medical studies and have received much attention recently. In such data, each study subject may potentially experience several types of events or recurrences of the same type of event, or event times may be clustered. Marginal distributions are specified for the multivariate event times in multiple events and clustered events data, and for the gap times in recurrent events data, using the semiparametric linear transformation models while leaving the dependence structures for related events unspecified. We propose several estimating equations for simultaneous estimation of the regression parameters and the transformation function. It is shown that the resulting regression estimators are asymptotically normal, with variance–covariance matrix that has a closed form and can be consistently estimated by the usual plug-in method. Simulation studies show that the proposed approach is appropriate for practical use. An application to the well-known bladder cancer tumor recurrences data is also given to illustrate the methodology.     

Nonparametric inference for the joint distribution of recurrent marked variables and recurrent survival time

Time between recurrent medical events may be correlated with the cost incurred at each event. As a result, it may be of interest to describe the relationship between recurrent events and recurrent medical costs by estimating a joint distribution. In this paper, we propose a nonparametric estimator for the joint distribution of recurrent events and recurrent medical costs in right-censored data. We also derive the asymptotic variance of our estimator, a test for equality of recurrent marker distributions, and present simulation studies to demonstrate the performance of our point and variance estimators. Our estimator is shown to perform well for a wide range of levels of correlation, demonstrating that our estimators can be employed in a variety of situations when the correlation structure may be unknown in advance. We apply our methods to hospitalization events and their corresponding costs in the second Multicenter Automatic Defibrillator Implantation Trial (MADIT-II), which was a randomized clinical trial studying the effect of implantable cardioverter-defibrillators in preventing ventricular arrhythmia.     

An additive marginal regression model for clustered recurrent event in the presence of a terminal event

In this article, an additive rate model is proposed for clustered recurrent event with a terminal event. The subjects are clustered by some property. For the clustered subjects, the recurrent event is precluded by the death. An estimating equation is developed for the model parameter and the baseline rate function. The asymptotic properties of the resulting estimators are established. In addition, a goodness-of-fit test is presented to assess the adequacy of the model. The finite-sample behavior of the proposed estimators is evaluated through simulation studies, and an application to a bladder cancer data is illustrated.     

Joint modeling approach for semicompeting risks data with missing nonterminal event status

Semicompeting risks data, where a subject may experience sequential non-terminal and terminal events, and the terminal event may censor the non-terminal event but not vice versa, are widely available in many biomedical studies. We consider the situation when a proportion of subjects’ non-terminal events is missing, such that the observed data become a mixture of “true” semicompeting risks data and partially observed terminal event only data. An illness–death multistate model with proportional hazards assumptions is proposed to study the relationship between non-terminal and terminal events, and provide covariate-specific global and local association measures. Maximum likelihood estimation based on semiparametric regression analysis is used for statistical inference, and asymptotic properties of proposed estimators are studied using empirical process and martingale arguments. We illustrate the proposed method with simulation studies and data analysis of a follicular cell lymphoma study.     

Nonparametric inference for panel count data with competing risks

In survival and reliability studies, panel count data arise when we investigate a recurrent event process and each study subject is observed only at discrete time points. If recurrent events of several types are possible, we obtain panel count data with competing risks. Such data arise frequently from transversal studies on recurrent events in demography, epidemiology and reliability experiments where the individuals cannot be observed continuously. In the present paper, we propose an isotonic regression estimator for the cause specific mean function of the underlying recurrent event process of a competing risks panel count data. Further, a nonparametric test is proposed to compare the cause specific mean functions of the panel count competing risks data. Asymptotic properties of the proposed estimator and test statistic are studied. A simulation study is conducted to assess the finite sample behaviour of the proposed estimator and test statistic. Finally, the procedures developed are applied to a real data arising from skin cancer chemo prevention trial.     

A semiparametric additive rates model for recurrent event data

Recurrent event data often arise in biomedical studies, with examples including hospitalizations, infections, and treatment failures. In observational studies, it is often of interest to estimate the effects of covariates on the marginal recurrent event rate. The majority of existing rate regression methods assume multiplicative covariate effects. We propose a semiparametric model for the marginal recurrent event rate, wherein the covariates are assumed to add to the unspecified baseline rate. Covariate effects are summarized by rate differences, meaning that the absolute effect on the rate function can be determined from the regression coefficient alone. We describe modifications of the proposed method to accommodate a terminating event (e.g., death). Proposed estimators of the regression parameters and baseline rate are shown to be consistent and asymptotically Gaussian. Simulation studies demonstrate that the asymptotic approximations are accurate in finite samples. The proposed methods are applied to a state-wide kidney transplant data set.     

Nonparametric analysis of recurrent gap time data

Abstract

Recurrent event data are frequently encountered in longitudinal studies. In many applications, the times between successive recurrent events (gap times) are often of interest and lead to problems that have received much attention recently. In this article, using the approach of inverse probability-of-censoring weights (IPCW), we propose nonparametric estimators for the estimation of the bivariate distribution and survival functions for gap times of recurrent event data. We also consider the estimation of Kendall’s tau for two gap times by expressing it as an integral functional of the bivariate survival function. The asymptotic properties of the proposed estimators are established. Simulation studies are conducted to investigate their finite sample performance.     

Estimation of the joint survival function for three successive duration times under double truncation and right censoring

In incident cohort studies, it is common to include subjects who have experienced a certain event within a calendar time window. For all the included individuals, the time of the previous events is retrospectively confirmed and the occurrence of subsequent events is observed during the follow-up periods. During the follow-up periods, subjects may undergo three successive events. Since the second/third duration process becomes observable only if the first/second event has occurred, the data is subject to double truncation and right censoring. We consider two cases: the case when the first event time is subject to double truncation and the case when the second event time is subject to double truncation. Using the inverse-probability-weighted approach, we propose nonparametric and semiparametric estimators for the estimation of the joint survival function of three successive duration times. We establish the asymptotic properties of the proposed estimators and conduct a simulation study to investigate the finite sample properties of the proposed estimators.