Joint modelling of pre-randomisation event counts and multiple post-randomisation survival times with cure rates: application to data for early epilepsy and single seizures. (English) Zbl 07265814

Summary: In this paper, we consider the analysis of recurrent event data that examines the differences between two treatments. The outcomes that are considered in the analysis are the pre-randomisation event count and post-randomisation times to first and second events with associated cure fractions. We develop methods that allow pre-randomisation counts and two post-randomisation survival times to be jointly modelled under a Poisson process framework, assuming that outcomes are predicted by (unobserved) event rates. We apply these methods to data that examine the difference between immediate and deferred treatment policies in patients presenting with single seizures or early epilepsy. We find evidence to suggest that post-randomisation seizure rates change at randomisation and following a first seizure after randomisation. We also find that there are cure rates associated with the post-randomisation times to first and second seizures. The increase in power over standard survival techniques, offered by the joint models that we propose, resulted in more precise estimates of the treatment effect and the ability to detect interactions with covariate effects.


62-XX Statistics
Full Text: DOI Link


[1] Akacha, M. and Benda, N. 2010. The impact of dropouts on the analysis of dose-finding studies with recurrent event data. Stat. Med., 29: 1635-1646.
[2] Berg, A. T. and Shinnar, S. 1991. The risk of seizure recurrence following a first unprovoked seizure: A quantitative review. Neurology, 41: 965-972. (doi:10.1212/WNL.41.7.965)
[3] Chandra, B. 1992. First seizure in adults: To treat or not to treat. Clin. Neurol. Neurosurg., 94(Suppl.): S61-S63. (doi:10.1016/0303-8467(92)90024-W)
[4] Cockerell, O. C., Johnson, A. L., Sander, J. W.A.S., Hart, Y. M. and Shorvon, S. D. 1995. Remission of epilepsy: Results from the national general practice study of epilepsy. Lancet, 346: 140-144. (doi:10.1016/S0140-6736(95)91208-8)
[5] Cook, R. J. and Lawless, J. F. 2007. The Statistical Analysis of Recurrent Events, New York, NY: Statistics for Biology and Health, Springer. · Zbl 1159.62061
[6] Cook, R. J. and Wei, W. 2003. Conditional analysis of mixed Poisson processes with baseline counts: Implications for trial design and analysis. Biostatistics, 4: 479-494. (doi:10.1093/biostatistics/4.3.479) · Zbl 1154.62387
[7] Cowling, B. J. 2003. Survival models for censored point processes, Coventry: Department of Statistics, University of Warwick. PhD thesis
[8] Cowling, B. J., Hutton, J. L. and Shaw, J. E.H. 2006. Joint modelling of event counts and survival times. J. R. Statist. Soc. C, 55: 31-39. (doi:10.1111/j.1467-9876.2005.00529.x) · Zbl 05188725
[9] Diggle, P. J., Sousa, I. and Chetwynd, A. G. 2008. Joint modelling of repeated measurements and time-to-event outcomes: The fourth Armitage lecture. Stat. Med., 27: 2981-2998. (doi:10.1002/sim.3131)
[10] Duchateau, L., Janssen, P., Kezic, I. and Fortpied, C. 2003. Evolution of recurrent asthma event rate over time in frailty models. J. R. Statist. Soc. C, 52: 355-363. (doi:10.1111/1467-9876.00409) · Zbl 1111.62336
[11] Filliben, J. J. 1975. The probability plot correlation coefficient test for normality. Technometrics, 11: 111-117. (doi:10.1080/00401706.1975.10489279) · Zbl 0295.62049
[12] Greenwood, M. and Yule, G. U. 1920. An inquiry into the nature of frequency distributions representative of multiple happenings with particular reference to the occurrence of multiple attacks of disease or of repeated accidents. J. R. Statist. Soc., 83: 255-279. (doi:10.2307/2341080)
[13] Hougaard, P., Lee, M. and Whitmore, G. A. 1997. Analysis of overdispersed count data by mixtures of Poisson variables and Poisson processes. Biometrics, 53: 1225-1238. (doi:10.2307/2533492) · Zbl 0911.62101
[14] Hutton, J. L. and Monaghan, P. F. 2002. Choice of parametric accelerated life and proportional hazards models for survival data: Asymptotic results. Lifetime Data Anal., 8: 375-393. (doi:10.1023/A:1020570922072) · Zbl 1116.62423
[15] ILAE Commission on Antiepileptic Drugs. 1998. Considerations on designing clinical trials to evaluate the place of new antiepileptic drugs in the treatment of newly diagnosed and chronic patients with epilepsy. Epilepsia, 39: 799-803. (doi:10.1111/j.1528-1157.1998.tb01167.x)
[16] Kim, L., Johnson, A., Marson, A. and Chadwick, D. 2006. Prediction of risk of seizure recurrence after a single seizure and early epilepsy: Further results from the MESS trial. Lancet Neurol., 5: 317-322. (doi:10.1016/S1474-4422(06)70383-0)
[17] Lin, D. Y., Sun, W. and Ying, Z. 1999. Nonparametric estimation of the gap time distributions for serial events with censored data. Biometrika, 86: 59-70. (doi:10.1093/biomet/86.1.59) · Zbl 0917.62030
[18] Little, R. J.A. and Rubin, D. B. 2002. Statistical Analysis with Missing Data, Hoboken, NJ: Wiley Series in Probability and Statistics, John Wiley and Sons.
[19] Maller, R. and Zhou, X. 1996. Survival Analysis with Long-Term Survivors, Chichester, , West Sussex: Wiley Series in Probability and Statistics, John Wiley and Sons. · Zbl 1151.62350
[20] Marson, A., Jacoby, A., Johnson, A., Kim, L., Gamble, C. and Chadwick, D. 2005. Immediate versus deferred antiepileptic drug treatment for early epilepsy and single seizures: A randomised control trial. Lancet, 365: 2007-2013. (doi:10.1016/S0140-6736(05)66694-9)
[21] McCullagh, P. and Nelder, J. A. 1989. Generalised Linear Models, Boca Raton, FL: Monographs on Statistics and Applied Probability, Chapman and Hall/CRC. · Zbl 0588.62104
[22] Nayak, T. K. 1987. Multivariate Lomax distribution: Properties and usefulness in reliability theory. J. Appl. Probab., 24: 170-177. (doi:10.2307/3214068) · Zbl 0615.60082
[23] Y. Peng, gfcure. Available at http://post.queensu.ca/ pengp/software.html (downloaded 14 July 2009).
[24] Rogers, J. K. 2011. Statistical models for censored point processes with cure fractions, Coventry: Department of Statistics, University of Warwick. PhD thesis
[25] J.K. Rogers, J.L. Hutton, and K. Hemming, Joint modelling of event counts and survival times, CRiSM Working Paper 44, University of Warwick, 2009. Available at http://www2.warwick.ac.uk/fac/sci/statistics/crism/research/2009/paper09-44.
[26] Rogers, J. K., Hutton, J., Marson, A. G. and Chadwick, D. W. 2012. Assessing the risk of subsequent tonic-clonic seizures in patients with a history of simple or complex partial seizures. J. Neurol. Neurosurg. Psychiatry, 83: 803-809. (doi:10.1136/jnnp-2011-300917)
[27] Verity, C. M., Hosking, G. and Easter, D. J. 1995. A multicentre comparative trial of sodium valproate and carbamazepine in paediatric epilepsy. Dev. Med. Child Neurol., 37: 97-108. (doi:10.1111/j.1469-8749.1995.tb11978.x)
[28] Visser, M. 1996. Nonparametric estimation of the bivariate survival function with an application to vertically transmitted AIDS. Biometrika, 83: 507-518. (doi:10.1093/biomet/83.3.507) · Zbl 0866.62076
[29] Warrell, D. A., Cox, T. M., Firth, J. D. and Benz, E. J. 2003. Oxford Textbook of Medicine, Edited by: Warrell, D. A., Cox, T. M., Firth, J. D. and Benz, E. J. Oxford: Oxford University Press.
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.