Abstract:Event of the same type occurring several times for one individual (recurrent events) are present in various domains (industrial systems reliability, episodes of unemployment, political conflicts, chronic diseases episodes). Analysis of such kind of data should account for the whole recurrence process dynamics rather than only focusing on the number of observed events. Statistical models for recurrent events analysis are developed in the counting process probabilistic framework. One of the often-used models is the Andersen-Gill model, a generalization of the well-known Cox model for durations, which assumes that the baseline intensity of the recurrence process is time-dependent and is adjusted for covariates. For an individual i with covariates Xi, the intensity is as follows: $\lambda$_{ik}(t;$\theta$) = $\lambda$_0(t) exp (X_i $\beta$). The baseline intensity can be specified parametrically, in a form of Weibull: $\lambda$_0 (t) = $\gamma$_{1} $\gamma$_{2} t^{$\gamma$_2-1}, with $\gamma$1 scale parameter et $\gamma$2 shape parameter. However, the observed covariates are often insufficient to explain the observed heterogeneity in data. This is often the case of clinical trials data containing information on patients. In this article a mixture model for recurrent events analysis is proposed. This model allows to account for unobserved heterogeneity and to cluster individuals according to their recurrence process. The intensity of the process is parametrically specified within each class and depend on observed covariates. Thus, the intensity becomes specific to class k: $\lambda$_{ik} (t; $\theta$_k) = $\gamma$_{1k} $\gamma$_{2k} t^{$\gamma$_{2k}-1} exp (X_i $\beta$_k). The model parameters are estimated by the Maximum Likelihood method, using the EM algorithm. The BIC criterion is employed to choose the optimal number of classes. Model feasibility is verified by Monte Carlo simulations. An application to real data concerning hospital readmissions of elderly patients is proposed. The proposed model feasibility is empirically verified (the optimization algorithm converges, providing non-biased estimates). The real data application allows to identify two clinically relevant classes of patients.

Some applications of phase-type distributions in recurrent events

Identifying Markov chain models from time-to-event data: an algebraic approach

Heavy-tailed phase-type distributions: A unified approach

Classification non supervis{é}e des processus d'{é}v{é}nements r{é}currents

Bayesian inference on the number of recurrent events: A joint model of recurrence and survival

Quantile Regression Modeling of Recurrent Event Risk

Recurrent Events Analysis in the Presence of Terminal Event and Zero-recurrence Subjects

Accelerated failure time models for recurrent event data analysis and joint modeling

Copula-Frailty Models for Recurrent Event Data Based on Monte Carlo EM Algorithm

Heterogeneous individual risk modelling of recurrent events

Proportional rates model for recurrent event data with intermittent gaps and a terminal event

Bayesian Nonparametric Modelling of Joint Gap Time Distributions for Recurrent Event Data

Sample size estimation for recurrent event data using multifrailty and multilevel survival models

Conditional modeling of recurrent event data with terminal event

Advanced models for predicting event occurrence in event-driven clinical trials accounting for patient dropout, cure and ongoing recruitment

Advanced method and application on Markov state transfer chain

Bayesian semi-parametric inference for clustered recurrent events with zero-inflation and a terminal event/4163305

Regression analysis of multivariate recurrent event data allowing time-varying dependence with application to stroke registry data

Joint modelling of recurrent and terminal events with discretely-distributed non-parametric frailty: application on re-hospitalizations and death in heart failure patients

Incorporating delayed entry into the joint frailty model for recurrent events and a terminal event

A competing risks approach for nonparametric estimation of transition probabilities in a non-Markov illness-death model