A simulation procedure based on copulas to generate clustered multi-state survival data

Authors:
Federico Rotolo;Catherine Legrand;Ingrid Van Keilegom
Affiliations:
Dipartimento di Scienze Statistiche, Universití di Padova, Via Cesare Battisti 241, 35121 Padova, Italy;Institut de Statistique, Biostatistique et Sciences Actuarielles, Université catholique de Louvain, Voie du Roman Pays, 20, 1348 Louvain-la-Neuve, Belgium;Institut de Statistique, Biostatistique et Sciences Actuarielles, Université catholique de Louvain, Voie du Roman Pays, 20, 1348 Louvain-la-Neuve, Belgium
Venue:
Computer Methods and Programs in Biomedicine
Year:
2013

Citing 5
Cited 0

The shared frailty model and the power for heterogeneity tests in multicenter trials

Computational Statistics & Data Analysis
An Introduction to Copulas (Springer Series in Statistics)

An Introduction to Copulas (Springer Series in Statistics)
Comparison of different estimation procedures for proportional hazards model with random effects

Computational Statistics & Data Analysis
The mstate package for estimation and prediction in non- and semi-parametric multi-state and competing risks models

Computer Methods and Programs in Biomedicine
Statistical models for recurrent events and death: Application to cancer events

Mathematical and Computer Modelling: An International Journal

Quantified Score

Hi-index	0.00

Visualization

Abstract

Generating survival data with a clustered and multi-state structure is useful to study finite sample properties of multi-state models, competing risks models and frailty models. We propose a simulation procedure based on a copula model for each competing events block, allowing to introduce dependence between times of different transitions and between those of grouped subjects. The effect of simulated frailties and covariates can be added in a proportional hazards way. In order to mimic information from real data, we also propose a method for the tuning of parameters via numerical minimization of a criterion function based on the ratios of target and observed values of median times and of probabilities of competing events. An example is provided on simulation of data mimicking those from a multicenter study on head and neck cancer, where the interest is in studying both time to local relapses and to distant metastases before death. The results demonstrated that data simulated according to our proposed method have characteristics very close to the target values.