Recurrent
event time data and more general multiple event time data are commonly analyzed
using extensions of Cox regression, or proportional hazards regression, as used
with single event time data. These methods treat covariates, either
time-invariant or time-varying, as having multiplicative effects while general
dependence on time is left un-estimated. An adaptive approach is formulated for
analyzing multiple event time data. Conditional hazard rates are modeled in
terms of dependence on both time and covariates using fractional polynomials
restricted so that the conditional hazard rates are positive-valued and so that
excess time probability functions (generalizing survival functions for single
event times) are decreasing. Maximum likelihood is used to estimate parameters
adjusting for right censored event times. Likelihood cross-validation (LCV)
scores are used to compare models. Adaptive searches through alternate
conditional hazard rate models are controlled by LCV scores combined with
tolerance parameters. These searches identify effective models for the
underlying multiple event time data. Conditional hazard regression is
demonstrated using data on times between tumor recurrence for bladder cancer
patients. Analyses of theory-based models for these data using extensions of
Cox regression provide conflicting results on effects to treatment group and
the initial number of tumors. On the other hand, fractional polynomial analyses
of these theory-based models provide consistent results identifying significant
effects to treatment group and initial number of tumors using both model-based
and robust empirical tests. Adaptive analyses further identify distinct
moderation by group of the effect of tumor order and an additive effect to
group after controlling for nonlinear effects to initial number of tumors and
tumor order. Results of example analyses indicate that adaptive conditional
hazard rate modeling can generate useful insights into multiple event time
data.
References
[1]
Kalbfleisch, J.D. and Schaubel, D.E. (2023) Fifty Years of the Cox Model. Annual Review of Statistics and Its Application, 10, 1-23.
https://doi.org/10.1146/annurev-statistics-033021-014043
[2]
Andersen, P.K. and Gill, R.D. (1982) Cox’s Regression Model for Counting Processes: A Large Sample Study. The Annals of Statistics, 10, 1100-1120.
https://doi.org/10.1214/aos/1176345976
[3]
Lin, D.Y., Wei, L.J., Yang, I. and Ying, Z. (2000) Semiparametric Regression for the Mean and Rate Functions of Recurrent Events. Journal of the Royal Statistical Society Series B, 62, 711-730. https://doi.org/10.1111/1467-9868.00259
[4]
Prentice, R.L., Williams, B.J. and Peterson, A.V. (1981) On the Regression Analysis of Multivariate Failure Time Data. Biometrika, 68, 373-379.
https://doi.org/10.1093/biomet/68.2.373
[5]
Wei, L.J., Lin, D.Y. and Weissfeld, L. (1989) Regression Analysis of Multivariate Incomplete Failure Time Data by Modeling Marginal Distributions. Journal of the American Statistical Association, 84, 1065-1073.
https://doi.org/10.1080/01621459.1989.10478873
[6]
Royston, P. and Altman, D.G. (1994) Regression Using Fractional Polynomials of Continuous Covariates: Parsimonious Parametric Modeling. Applied Statistics, 43, 429-467. https://doi.org/10.2307/2986270
[7]
Knafl, G.J. (2018) Adaptive Fractional Polynomial Modeling. Open Journal of Statistics, 8, 159-186. https://doi.org/10.4236/ojs.2018.81011
[8]
Knafl, G.J. (2023) An Adaptive Approach for Hazard Regression Modeling. Open Journal of Statistics, 13, 300-315. https://doi.org/10.4236/ojs.2023.133016
[9]
Amorim, L. and Cai, J. (2014) Modelling Recurrent Events: A Tutorial for Analysis in Epidemiology. International Journal of Epidemiology, 44, 324-333.
https://doi.org/10.1093/ije/dyu222
[10]
Royston, P. and Sauerbrei, W. (2008) Multivariable Model-Building: A Practical Approach to Regression Analysis Based on Fractional Polynomials for Modelling Continuous Variables. John Wiley & Sons, Hoboken.
https://doi.org/10.1002/9780470770771
[11]
Knafl, G.J. and Ding, K. (2016) Adaptive Regression for Modeling Nonlinear Relationships. Springer International Publishing, Berlin.
[12]
Baron, R.M. and Kenny, D.A. (1986) The Moderator-Mediator Variable Distinction in Social Psychology Research: Conceptual, Strategic, and Statistical Considerations. Journal of Personality & Social Psychology, 51, 1173-1182.
https://doi.org/10.1037/0022-3514.51.6.1173
[13]
Knafl, G.J. (2018) A Reassessment of Birth Defects for Children of Participants of the Air Force Health Study. Open Journal of Epidemiology, 8, 187-200.
https://doi.org/10.4236/ojepi.2018.84015
[14]
Knafl, G.J. and Meghani, S.H. (2021) Modeling Individual Patient Count/Rate Data over Time with Applications to Cancer Pain Flares and Cancer Pain Medication Usage. Open Journal of Statistics, 11, 633-654.
https://doi.org/10.4236/ojs.2021.115038
[15]
Knafl, G.J. (2022) Adaptive Regression for Nonlinear Interrupted Time Series Analyses with Application to Birth Defects in Children of Vietnam War Veterans. Open Journal of Statistics, 12, 789-809. https://doi.org/10.4236/ojs.2022.126045
[16]
Knafl, G.J. and Meghani, S.H. (2021) Regression Modeling of Individual-Patient Correlated Discrete Outcomes with Applications to Cancer Pain Ratings. Open Journal of Statistics, 12, 456-485. https://doi.org/10.4236/ojs.2022.124029
[17]
Byar, D.P. (1980) The Veterans Administration Study of Chemoprophylaxis for Recurrent Stage I Bladder Tumors: Comparisons of Placebo, Pyridoxine, and Topical Thiotepa. In: Pavone-Macaluso, M., Smith, P.H. and Edsmyn, F., Eds., Bladder Tumors and Other Topics in Urological Oncology, Plenum, New York, 363-370.
https://doi.org/10.1007/978-1-4613-3030-1_74
[18]
SAS Institute Inc. (2004) SAS/STAT 9.1 User’s Guide: Volume 5. SAS Institute Inc., Cary.
[19]
Lee, E., Wei, L. and Amato, D. (1992) Cox-Type Regression Analyses for Large Numbers of Small Groups of Correlated Failure Time Observations. In: Klein, J.P. and Goel, P.K., Eds., Survival Analysis: State of the Art, Nato Science, Vol. 211, Springer, Dordrecht, 237-247. https://doi.org/10.1007/978-94-015-7983-4_14