A Dynamic Risk Model for Multitype Recurrent Events

Recurrent events can occur more than once in the same individual; such events may be of different types, known as multitype recurrent events. They are very common in longitudinal studies. Often there is a terminating event, after which no further events can occur. The risk of any event, including terminating events such as death or cure, is typically affected by prior events. We propose a flexible joint multitype recurrent-events model that explicitly provides estimates of the change in risk for each event due to subject characteristics, including number and type of prior events and the absolute risk for every event type (terminating and nonterminating), and predicts event-free survival probability over a desired time period. The model is fully parametric, and therefore a standard likelihood function and robust standard errors can be constructed. We illustrate the model with applications to the Antihypertensive and Lipid-Lowering Treatment to Prevent Heart Attack Trial (1994-2002) and provide discussion of the results and model features.

Multistate modeling and structure selection for multitype recurrent events and terminal event data

In cardiovascular disease studies, a large number of risk factors are measured but it often remains unknown whether all of them are relevant variables and whether the impact of these variables is changing with time or remains constant. In addition, more than one kind of cardiovascular disease events can be observed in the same patient and events of different types are possibly correlated. It is expected that different kinds of events are associated with different covariates and the forms of covariate effects also vary between event types. To tackle these problems, we proposed a multistate modeling framework for the joint analysis of multitype recurrent events and terminal event. Model structure selection is performed to identify covariates with time-varying coefficients, time-independent coefficients, and null effects. This helps in understanding the disease process as it can detect relevant covariates and identify the temporal dynamics of the covariate effects. It also provides a more parsimonious model to achieve better risk prediction. The performance of the proposed model and selection method is evaluated in numerical studies and illustrated on a real dataset from the Atherosclerosis Risk in Communities study.