Events that may occur repeatedly for individual subjects are of interest in many medical studies. We review methods of analysis for repeated events, emphasizing that the approach taken in a given study should allow clinical questions to be addressed as directly as possible. Methods based on full models for event processes as well as on simpler ‘marginal’ assumptions are considered. The treatment of dependent terminating events related to the recurrent events is also discussed. We apply various methods of analysis to studies involving pulmonary exacerbations in persons with cystic fibrosis, and the occurrence of bone metastases and skeletal events in cancer patients, respectively. Most of the methodology considered can be implemented with existing software.
Get full access to this article
View all access options for this article.
References
1.
Hobson RW , Weiss DG , Fields WSet al.The Veterans Affairs Cooperative Study Group . Effect of carotid endarterectomy for asymptomatic carotid stenosis . New England Journal of Medicine1993; 328: 221-227 .
2.
Byar D , Kaihara R , Sylvester Ret al.Statistical analysis techniques and sample size determination for clinical trials of treatments for bladder cancer. Developments in bladder cancer. New York: Alan R. Liss , 1986, 49-64.
3.
Fuchs HJet al.The effect of aerosolized recombinant human DNase on respiratory exacerbations and pulmonary function in patients with cystic fibrosis . New England Journal of Medicine1994; 331: 637-642 .
4.
Sears MR , Taylor DR , Print CGet al.Regular inhaled beta-agonist treatment in bronchial asthma . Lancet1990; 336: 1391-1396 .
5.
Ng ETM , Cook RJ.Robust inference for bivariate point processes . The Canadian Journal of Statistics1999; 27: 509-524 .
6.
Tuli S , Drake J , Lawless J , Wigg M , Lamberti-Pasculli M.Risk factors for repeat cerebrospinal shunt failures in pediatric hydrocephalus . Journal of Neurosurgery2000; 92: 31-38 .
7.
Hortobagyi GN , Thierault RL , Porter Let al.Efficacy of pamidronate in reducing skeletal complications in patients with breast cancer and lytic bone metastases . New England Journal of Medicine1996; 335: 1785-1791 .
8.
Hortobagyi GN , Theriault RL , Liptonet al.Long-term prevention of skeletal complications of metastic breast cancer with pamidronate . Journal of Clinical Oncology1998; 16: 2031-2044 .
9.
Kalbfleisch JD , Prentice RL.The statistical analysis of failure time data. New York: Wiley , 1980.
10.
Andersen PK , Borgan O , Gill RD , Keiding N.Statistical models based on counting processes. New York: Springer-Verlag , 1993.
11.
Berman M , Turner TR.Approximating point process likelihoods . Applied Statistics1992; 41: 31-38 .
12.
Aalen OO , Husebye E.Statistical analysis of repeated events forming renewal processes . Statistics in Medicine1991; 10: 1227-1240 .
13.
Hougaard P.Analysis of multivariate survival data. New York: Springer-Verlag , 2000.
14.
Lawless JF.Regression methods for Poisson process data . Journal of the American Statistical Association1987; 82: 808-815 .
15.
Dean CB , Balshaw R.Efficiency lost by analysing counts rather than event times in Poisson and overdispersed Poisson regression models . Journal of the American Statistical Association1997; 92: 1387-1398 .
16.
Thall PF , Vail SC.Some covariance models for longitudinal count data with overdispersion . Biometrics1995; 46: 657-671 .
17.
Lawless JF , Nadeau JC.Some simple robust methods for the analysis of recurrent events . Technometrics1995; 37: 158-168 .
18.
Godambe VP , Kale BK.Estimating functions: an overview. In: Godambe VP , ed. Estimating functions. Oxford: Clarendon Press , 1991.
19.
Cook RJ , Lawless JF , Nadeau JC.Robust tests for treatment comparisons based on recurrent event responses . Biometrics1996; 52: 732-739 .
20.
Cook RJ , Lawless JF.Interim monitoring of longitudinal comparative studies with recurrent event responses . Biometrics1996; 52: 1311-1323 .
21.
Thall PF , Lachin JM.Analysis of recurrent events: nonparametric methods for random interval count data . Journal of the American Statistical Association1988; 83: 339-347 .
22.
Sun J , Kalbfleisch JD.The analysis of current status data on point processes . Journal of the American Statistical Association1993; 88: 1449-1454 .
23.
Sun J , Kalbfleisch JD.Estimation of the mean function of point processes based on panel count data . Statistics Sinica1995; 5: 279-290 .
24.
Thall PF.Mixed Poisson likelihood regression models for longitudinal interval count data . Biometrics1988; 44: 197-209 .
Lawless JF , Zhan M.Analysis of interval-grouped recurrent-event data using piecewise constant rate functions . Canadian Journal of Statistics1998; 4: 549-565 .
27.
Lawless JF , Fong D.State duration models in clinical and observational studies . Statistics in Medicine1999; 18: 2365-2376 .
28.
Cox DR.Regression models with lifetables (with discussion) . Journal of the Royal Statistical Society, Series B1972; 34: 187-220 .
29.
Gail M , Santner T , Brown C.An analysis of comparative carcinogenesis experiments based on multiple times to tumor . Biometrics1980; 36: 255-266 .
30.
Prentice RL , Williams BJ , Peterson AV.On the regression analysis of multivariate failure time data . Biometrika1981; 68: 373-379 .
31.
Dabrowska DM , Sun G , Horwitz MM.Cox regression in a Markov renewal model: an application to the analysis of bone marrow transplant data . Journal of the American Statistical Association1994; 89: 867-877 .
32.
Lawless JF , Wigg MB , Tuli S , Drake J , Lamberti-Pasculli M.Analysis of repeated shunt failures or durations, with applications to shunt failures for patients with paediatric hydrocephalus . Applied Statistics2001; 50: 449-465 .
33.
Wei LJ , Lin DY , Weissfeld L.Regression analysis of multivariate incomplete failure time data by modeling marginal distributions . Journal of the American Statistical Association1989; 84: 1065-1073 .
34.
Lee EW , Wei LJ , Amato DA.Cox-type regression analysis for large numbers of small groups of correlated failure time observations. In: Klein JP , Goel PK , eds. Survival analysis: state of the art. New York: Kluwer , 1992, pp 237-247.
35.
Visser M.Nonparametric estimation of the bivariate survival function with an application to vertically transmitted AIDS . Biometrika1996; 83: 507-518 .
36.
Wang W , Wells MT.Nonparametric estimation of successive duration times under dependent censoring . Biometrika1998; 85: 561-572 .
37.
Lin DY , Sun W , Ying Z.Nonparametric estimation of the gap time distribution for serial events with censored data . Biometrika1999; 86: 59-70 .
38.
Follman D , Goldberg M.Distinguishing heterogeneity from decreasing hazard rates . Technometrics1988; 30: 389-396 .
39.
Therneau TM , Grambsch PM.Modeling survival data: extending the cox model. New York: Springer , 2000.
40.
Lawless JF , Thiagarajah K.A point process model incorporating renewals and time trends . Technometrics1996; 38: 131-138 .
41.
Cox DR.The statistical analysis of dependencies in point processes. In: Lewis PAW , ed. Symposium on point processes. New York: Wiley , 1972.
42.
Cox DR , Isham V.Point processes. London: Chapman and Hall , 1980.
43.
Daley DJ , Vere-Jones O.An Introduction to the theory of point processes. New York: Springer , 1988.
44.
Lawless JF.The analysis of recurrent events for multiple subjects . Applied Statistics1995; 44: 487-498 .
45.
Cole EH , Cattran DC , Farewell VTet al.A comparison of rabbit anti-thymocyte serum and OKT3 as prophylaxis against renal allograft rejection . Transplantation1994; 57: 60-67 .
Cook RJ , Lawless JF.Marginal analysis of recurrent events and a terminating event . Statistics in Medicines1997; 16: 911-924 .
48.
Kaplan EL , Meier P.Non-parametric estimation from incomplete observations . Journal of the American Statistical Association1958; 53: 457-481 .
49.
Ghosh D , Lin DY.Nonparametric analysis of recurrent events and death . Biometrics2000; 56: 554-562 .
50.
Strawderman R.Estimating the mean of an increasing stochastic process at a censored stopping time . Journal of the American Statistical Association2000; 95: 1192-1208 .
51.
Lin DY , Feuer EJ , Etzioni R , Wax Y.Estimating medical costs from incomplete follow up data . Biometrics1997; 53: 419-434 .
52.
Zhao H , Tsiatis AA.A consistent estimator for the distribution of quality-adjusted survival time . Biometrika1997; 84: 339-348 .
53.
Robins JM , Rotnitzky A.Recovery of information and adjustment for dependent censoring using surrogate markers. In: Jewell N , Dietz K , Farewell V , eds. AIDS epidemiology-methodological issues. Boston: Birkaüuser , 1992, pp 297-331.
54.
Lawless JF.Statistical models and methods for lifetime data, 2nd edn.New York: Wiley , 2002.
55.
Gray RJ.A class of K-sample tests for comparing the cumulative incidence of a competing risk . Annals of Statistics1988; 16: 1141-1154 .
56.
Cook RJ , Major P.Methodology for treatment evaluation in patients with cancer metastatic to bone . Journal of the National Cancer Institute2001; 93: 534-538 .
57.
Cook RJ.The design and analysis of randomized trials with recurrent events . Statistics in Medicine1995; 14: 2081-2098 .
58.
Hughes MD.Power considerations for clinical trials using multivariate time-to-event data . Statistics in Medicine1997; 16: 865-882 .
