Functional Cox model for interval-censored data
Pith reviewed 2026-05-19 15:56 UTC · model grok-4.3
The pith
A functional Cox model links continuous covariate trajectories to interval-censored event times with efficient estimators.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We formulate the effects of both scalar and functional covariates on the interval-censored event time through a functional Cox model. Penalized maximum likelihood estimation with an EM algorithm yields estimators that are consistent and asymptotically normal with limiting covariance matrices attaining the semiparametric efficiency bound, which can be estimated via the profile likelihood method. This supports a global test for the functional covariate effect.
What carries the argument
Penalized maximum likelihood estimation in the functional Cox model, computed via EM algorithm, for handling interval-censored data with functional covariates.
If this is right
- Regression parameters and functionals of the coefficient function can be consistently estimated.
- Confidence intervals and tests can be constructed using the asymptotic normality and efficient covariance estimates.
- A global test for the overall effect of the functional covariate is available.
- The approach applies to real datasets like those from Alzheimer's Disease Neuroimaging Initiative.
Where Pith is reading between the lines
- This framework may be adaptable to other semiparametric survival models with functional predictors.
- Potential extensions could include time-varying functional effects or integration with machine learning for high-dimensional functional data.
- Validation on additional datasets with different censoring patterns would strengthen generalizability.
Load-bearing premise
The functional Cox model correctly captures the relationship between the functional covariate trajectory and the hazard of the interval-censored event time, with the censoring mechanism independent of the event time conditional on the covariates.
What would settle it
Empirical evidence that the estimators do not achieve the semiparametric efficiency bound in large simulated samples from the model, or that the global test has incorrect size under the null hypothesis of no functional effect.
Figures
read the original abstract
Interval-censored data arise frequently in scientific studies, where the event of interest is known only to occur within a specific time interval. In such studies, functional covariates taking the form of continuous curves or spatial profiles are increasingly encountered, and it is of substantial scientific relevance to investigate how the trajectory of a functional covariate affects the event time. We formulate the effects of both scalar and functional covariates on the interval-censored event time through a functional Cox model. We consider penalized maximum likelihood estimation for this model and devise an EM algorithm to stably compute the parameter estimators. The resulting estimators for the regression parameters and linear functionals of the coefficient function are shown to be consistent and asymptotically normal, with limiting covariance matrices that attain the semiparametric efficiency bound and can be readily estimated through the profile likelihood method. Building upon these results, we construct a global test for the overall effect of the functional covariate. Finally, we assess the performance of the proposed methods through extensive simulation studies and present an application to data from the Alzheimer's Disease Neuroimaging Initiative.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes a functional Cox proportional hazards model for interval-censored event times that incorporates both scalar covariates and a functional covariate trajectory. Penalized maximum likelihood estimation is developed and implemented via an EM algorithm. The authors establish consistency and asymptotic normality of the estimators for the regression coefficients and for linear functionals of the coefficient function, with the limiting covariance attaining the semiparametric efficiency bound and estimable via the profile likelihood. A global test for the overall effect of the functional covariate is constructed. Performance is assessed in simulation studies and illustrated on Alzheimer's Disease Neuroimaging Initiative data.
Significance. If the asymptotic claims hold under the stated regularity conditions, the work supplies a practical and theoretically grounded method for a common data structure in longitudinal medical studies. The explicit attainment of the efficiency bound together with a profile-likelihood variance estimator and a global test constitute clear methodological contributions. The real-data example further demonstrates applicability.
major comments (2)
- [§3.2, Theorem 2] §3.2, Theorem 2: The proof of asymptotic normality and efficiency attainment for the estimator of the functional coefficient relies on the penalty parameter λ_n converging to zero at a specific rate; the manuscript should state the admissible range for λ_n explicitly in the theorem statement rather than only in the proof sketch, as this rate is load-bearing for the claimed √n-consistency.
- [Table 2] Table 2, n=200 row: the reported coverage probabilities for the pointwise confidence bands of the functional coefficient fall to 0.82–0.87 under moderate censoring; this under-coverage undermines the practical utility of the efficiency-bound variance estimator and requires either a larger-sample study or a finite-sample correction.
minor comments (3)
- [§2.3] The EM algorithm description in §2.3 does not specify the convergence criterion or the maximum number of iterations used in the reported simulations; adding these details would improve reproducibility.
- [Figure 3] Figure 3 caption: the label 'proposed' is ambiguous because two penalized variants are compared; please distinguish them as 'proposed (L2 penalty)' and 'proposed (L1 penalty)'.
- [§4] The global test in §4 is presented only for the functional covariate; a brief remark on how the test extends when both scalar and functional covariates are present would clarify its scope.
Simulated Author's Rebuttal
We are grateful to the referee for the thorough review and the recommendation for minor revision. The comments are helpful in improving the clarity and presentation of our work. Below we provide point-by-point responses to the major comments.
read point-by-point responses
-
Referee: [§3.2, Theorem 2] §3.2, Theorem 2: The proof of asymptotic normality and efficiency attainment for the estimator of the functional coefficient relies on the penalty parameter λ_n converging to zero at a specific rate; the manuscript should state the admissible range for λ_n explicitly in the theorem statement rather than only in the proof sketch, as this rate is load-bearing for the claimed √n-consistency.
Authors: We agree with the referee that the rate condition on the penalty parameter λ_n is essential for establishing the √n-consistency and semiparametric efficiency in Theorem 2. In the revised manuscript, we will explicitly state the admissible range for λ_n directly in the theorem statement, moving it from the proof sketch to the main theorem for better clarity. revision: yes
-
Referee: [Table 2] Table 2, n=200 row: the reported coverage probabilities for the pointwise confidence bands of the functional coefficient fall to 0.82–0.87 under moderate censoring; this under-coverage undermines the practical utility of the efficiency-bound variance estimator and requires either a larger-sample study or a finite-sample correction.
Authors: The referee correctly observes the finite-sample under-coverage in the coverage probabilities for n=200 under moderate censoring in Table 2. While the asymptotic results guarantee good performance for large n, we acknowledge that at moderate sample sizes the coverage can be below nominal levels. In the revision, we will include additional simulation results for larger sample sizes to demonstrate the convergence to nominal coverage, and add a brief discussion on the finite-sample behavior of the profile likelihood variance estimator. revision: yes
Circularity Check
No significant circularity
full rationale
The derivation proceeds from the functional Cox model specification under interval censoring to penalized MLE and EM algorithm, then invokes standard semiparametric regularity conditions to establish consistency, asymptotic normality, and attainment of the efficiency bound via profile likelihood. These steps rely on external mathematical results for Cox-type models extended to functional covariates rather than any self-definition, fitted-input renaming, or self-citation chain that reduces the central claims to the paper's own inputs by construction. The abstract and described results contain no equations or arguments that equate a prediction or estimator directly to a fitted quantity or prior self-cited ansatz.
Axiom & Free-Parameter Ledger
free parameters (1)
- penalty parameter
axioms (2)
- domain assumption The hazard rate follows a functional Cox form that additively incorporates scalar covariates and a linear functional of the coefficient function applied to the functional covariate.
- domain assumption Censoring is non-informative conditional on the observed covariates.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We formulate the effects of both scalar and functional covariates on the interval-censored event time through a functional Cox model. We consider penalized maximum likelihood estimation for this model and devise an EM algorithm...
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Theorem 1. The penalized maximum likelihood estimator for β is bβ(s) = sum dj ξj(s) + sum ci ∫ Zi(w) K1(w,s) dw
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
D. R. Cox , Journal =
-
[2]
D. Y. Lin and L. J. Wei and Z. Ying , journal =. Checking the
-
[3]
S. A. Murphy and A. W. van der Vaart , journal =. On Profile Likelihood , urldate =
-
[4]
Farrington, C. P. , title =. Biometrics , volume =. doi:https://doi.org/10.1111/j.0006-341X.2000.00473.x , year =
-
[5]
Biometrical Journal , volume =
Wang, Lianming and Sun, Liuquan and Sun, Jianguo , title =. Biometrical Journal , volume =. doi:https://doi.org/10.1002/bimj.200510269 , year =
-
[6]
Zeng, Donglin and Mao, Lu and Lin, D. Y. , title = ". Biometrika , volume =. 2016 , month =
work page 2016
-
[7]
Barlow, William E. and Prentice, Ross L. , title = ". Biometrika , volume =. 1988 , month =
work page 1988
-
[8]
Therneau, Terry M. and Grambsch, Patricia M. and Fleming, Thomas R. , title = ". Biometrika , volume =. 1990 , month =
work page 1990
-
[9]
Lagakos, S. W. , title = ". Biometrika , volume =. 1981 , month =
work page 1981
-
[10]
The ARIC (Atherosclerosis Risk In Communities) Study: JACC Focus Seminar 3/8 , journal =. 2021 , issn =. doi:https://doi.org/10.1016/j.jacc.2021.04.035 , url =
-
[11]
Xu, Yangjianchen and Zeng, Donglin and Lin, D Y , title = ". Biometrika , volume =. 2022 , month =
work page 2022
-
[12]
Zeng, Donglin and Gao, Fei and Lin, D. Y. , title = ". Biometrika , volume =. 2017 , month =
work page 2017
-
[13]
Xu, Yangjianchen and Zeng, Donglin and Lin, Dan-Yu , title = ". Biometrics , volume =. 2024 , month =
work page 2024
-
[14]
Zeng, Donglin and Lin, D Y , title = ". Biometrika , volume =. 2020 , month =
work page 2020
-
[15]
Zhang, Ying and Liu, Wei and Zhan, Yihui , title = ". Biometrika , volume =. 2001 , month =
work page 2001
- [16]
-
[17]
Cox, D. R. , title = ". Biometrika , volume =. 1975 , month =
work page 1975
-
[18]
Breslow, N. E. , Journal =
-
[19]
P. K. Andersen and R. D. Gill , title =. The Annals of Statistics , number =
-
[20]
The Annals of Statistics , number =
Jian Huang , title =. The Annals of Statistics , number =
-
[21]
Interval Censored Survival Data: A Review of Recent Progress
Huang, Jian and Wellner, Jon A. Interval Censored Survival Data: A Review of Recent Progress. Proceedings of the First Seattle Symposium in Biostatistics. 1997
work page 1997
-
[22]
Weak Convergence and Empirical Processes: with Applications to Statistics , author=. 1996 , publisher=
work page 1996
-
[23]
New England Journal of Medicine , volume =
Dan-Yu Lin and Yu Gu and Bradford Wheeler and Hayley Young and Shannon Holloway and Shadia-Khan Sunny and Zack Moore and Donglin Zeng , title =. New England Journal of Medicine , volume =
-
[24]
Lagakos, S. W. , title = ". Biometrika , volume =. 1988 , month =
work page 1988
-
[25]
New England Journal of Medicine , volume =
Dan-Yu Lin and Yu Gu and Yangjianchen Xu and Donglin Zeng and Bradford Wheeler and Hayley Young and Shadia Khan Sunny and Zack Moore , title =. New England Journal of Medicine , volume =
-
[26]
Zeng, D. and Lin, D. Y. , title =. Journal of the Royal Statistical Society Series B: Statistical Methodology , volume =. 2007 , month =
work page 2007
-
[27]
S. W. Lagakos and D. A. Schoenfeld , journal =. Properties of Proportional-Hazards Score Tests under Misspecified Regression Models , urldate =
-
[28]
DiRienzo, A. G. and Lagakos, S. W. , title =. Journal of the Royal Statistical Society Series B: Statistical Methodology , volume =. 2002 , month =
work page 2002
-
[29]
Gail, M. H. and Wieand, S. and Piantadosi, S. , title =. Biometrika , volume =. 1984 , month =
work page 1984
-
[30]
Struthers, C. A. and Kalbfleisch, J. D. , title =. Biometrika , volume =. 1986 , month =
work page 1986
-
[31]
Garnet L. Anderson and Thomas R. Fleming , journal =. Model Misspecification in Proportional Hazards Regression , urldate =
-
[32]
The Annals of Statistics , number =
Ker-Chau Li and Naihua Duan , title =. The Annals of Statistics , number =
-
[33]
Kosorok and Bee Leng Lee and Jason P
Michael R. Kosorok and Bee Leng Lee and Jason P. Fine , title =. The Annals of Statistics , number =. 2004 , doi =
work page 2004
-
[34]
Holme,. Effect of Flexible Sigmoidoscopy Screening on Colorectal Cancer Incidence and Mortality: A Randomized Clinical Trial , journal =. 2014 , month =
work page 2014
-
[35]
Donald M. Stablein and I. A. Koutrouvelis , journal =. A Two-Sample Test Sensitive to Crossing Hazards in Uncensored and Singly Censored Data , urldate =
-
[36]
Autologous Hematopoietic Stem Cell Transplantation vs Intravenous Pulse Cyclophosphamide in Diffuse Cutaneous Systemic Sclerosis: A Randomized Clinical Trial , journal =. 2014 , month =
work page 2014
- [37]
-
[38]
Timothy M. Morgan and S. W. Lagakos and D. A. Schoenfeld , journal =. Omitting Covariates from the Proportional Hazards Model , urldate =
-
[39]
Inconsistency of the misspecified proportional hazards model , journal =. 1986 , issn =. doi:10.1016/0167-7152(86)90091-X , url =
-
[40]
Anderson, Garnet L. and Fleming, Thomas R. , title =. Biometrika , volume =. 1995 , month =
work page 1995
-
[41]
Gerds, T. A. and Schumacher, M. , title =. Biometrika , volume =. 2001 , month =
work page 2001
-
[42]
D. Y. Lin and L. J. Wei , title =. Journal of the American Statistical Association , volume =. 1989 , publisher =
work page 1989
-
[43]
The Annals of Statistics , number =
Peter Sasieni , title =. The Annals of Statistics , number =. 1993 , doi =
work page 1993
-
[44]
Bhatt and Michael Szarek and P
Deepak L. Bhatt and Michael Szarek and P. Gabriel Steg and Christopher P. Cannon and Lawrence A. Leiter and Darren K. McGuire and Julia B. Lewis and Matthew C. Riddle and Adriaan A. Voors and Marco Metra and Lars H. Lund and Michel Komajda and Jeffrey M. Testani and Christopher S. Wilcox and Piotr Ponikowski and Renato D. Lopes and Subodh Verma and Pablo ...
work page 2021
-
[45]
Rosenberg, Eli S. and Dufort, Elizabeth M. and Udo, Tomoko and Wilberschied, Larissa A. and Kumar, Jessica and Tesoriero, James and Weinberg, Patti and Kirkwood, James and Muse, Alison and DeHovitz, Jack and Blog, Debra S. and Hutton, Brad and Holtgrave, David R. and Zucker, Howard A. , title =. JAMA , volume =. 2020 , month =
work page 2020
-
[46]
Jean-Claude Tardif and Simon Kouz and David D. Waters and Olivier F. Bertrand and Rafael Diaz and Aldo P. Maggioni and Fausto J. Pinto and Reda Ibrahim and Habib Gamra and Ghassan S. Kiwan and Colin Berry and José López-Sendón and Petr Ostadal and Wolfgang Koenig and Denis Angoulvant and Jean C. Grégoire and Marc-André Lavoie and Marie-Pierre Dubé and Dav...
work page 2019
-
[47]
Yangjianchen Xu and Donglin Zeng and D. Y. Lin , title =. Journal of the American Statistical Association , volume =. 2025 , publisher =. doi:10.1080/01621459.2025.2520460 , URL =
-
[48]
Bickel, Peter J. and Klaassen, Chris A. J. and Ritov, Ya'acov and Wellner, Jon A. , TITLE =. 1993 , PAGES =
work page 1993
-
[49]
Bycroft, Clare and Freeman, Colin and Petkova, Desislava and Band, Gavin and Elliott, Lloyd T. and Sharp, Kevin and Motyer, Allan and Vukcevic, Damjan and Delaneau, Olivier and O'Connell, Jared and Cortes, Adrian and Welsh, Samantha and Young, Alan and Effingham, Mark and McVean, Gil and Leslie, Stephen and Allen, Naomi and Donnelly, Peter and Marchini, J...
-
[50]
Dey, Rounak and Zhou, Wei and Kiiskinen, Tuomo and Havulinna, Aki and Elliott, Amanda and Karjalainen, Juha and Kurki, Mitja and Qin, Ashley and Lee, Seunggeun and Palotie, Aarno and Neale, Benjamin and Daly, Mark and Lin, Xihong and FinnGen , date =. Efficient and accurate frailty model approach for genome-wide survival association analysis in large-scal...
-
[51]
Spline Models for Observational Data , author=. 1990 , Publisher =
work page 1990
-
[52]
Yuan, Ming and Cai, T Tony , journal=. A reproducing kernel. 2010 , publisher=
work page 2010
-
[53]
Optimal penalized function-on-function regression under a reproducing kernel
Sun, Xiaoxiao and Du, Pang and Wang, Xiao and Ma, Ping , journal=. Optimal penalized function-on-function regression under a reproducing kernel. 2018 , publisher=
work page 2018
-
[54]
Optimal estimation for the functional
Qu, Simeng and Wang, Jane-Ling and Wang, Xiao , journal=. Optimal estimation for the functional
-
[55]
Statistical Methods in Medical Research , volume=
A functional proportional hazard cure rate model for interval-censored data , author=. Statistical Methods in Medical Research , volume=. 2022 , publisher=
work page 2022
-
[56]
Semiparametric inference for the functional
Hao, Meiling and Liu, Kin-yat and Xu, Wei and Zhao, Xingqiu , journal=. Semiparametric inference for the functional. 2021 , publisher=
work page 2021
-
[57]
The Annals of Statistics , volume=
Penalized log-likelihood estimation for partly linear transformation models with current status data , author=. The Annals of Statistics , volume=
-
[58]
Minimax nonparametric classification
Yang, Yuhong , journal=. Minimax nonparametric classification. 1999 , publisher=
work page 1999
-
[59]
Estimating time-varying effects for overdispersed recurrent events data with treatment switching , author=. Biometrika , volume=. 2013 , publisher=
work page 2013
- [60]
-
[61]
Gaussian Processes for Machine Learning , author=. 2006 , publisher=
work page 2006
-
[62]
Approximation Theory and Applications , editor=
Design problems for optimal surface interpolation , author=. Approximation Theory and Applications , editor=. 1981 , publisher=
work page 1981
-
[63]
R. C. Petersen and P. S. Aisen and L. A. Beckett and M. C. Donohue and A. C. Gamst and D. J. Harvey and C. R. Jack and W. J. Jagust and L. M. Shaw and A. W. Toga and J. Q. Trojanowski and M. W. Weiner , title =. Neurology , volume =. 2010 , doi =
work page 2010
- [64]
-
[65]
Piecewise-polynomial approximations of functions of the classes
Birman, M. Piecewise-polynomial approximations of functions of the classes. Mathematics of the USSR-Sbornik , volume=
-
[66]
Applied Functional Analysis: Applications to Mathematical Physics , author=. 1995 , publisher=
work page 1995
-
[67]
Bracketing metric entropy rates and empirical central limit theorems for function classes of
Nickl, Richard and P. Bracketing metric entropy rates and empirical central limit theorems for function classes of. Journal of Theoretical Probability , volume=. 2007 , publisher=
work page 2007
-
[68]
Focal decline of cortical thickness in
Lerch, Jason P and Pruessner, Jens C and Zijdenbos, Alex and Hampel, Harald and Teipel, Stefan J and Evans, Alan C , journal=. Focal decline of cortical thickness in. 2005 , publisher=
work page 2005
-
[69]
Neural network on interval-censored data with application to the prediction of
Sun, Tao and Ding, Ying , journal=. Neural network on interval-censored data with application to the prediction of. 2023 , publisher=
work page 2023
-
[70]
Dai, Ning and Kang, Hakmook and Jones, Galin L and Fiecas, Mark B and Alzheimer's Disease Neuroimaging Initiative , journal=. A. 2021 , publisher=
work page 2021
-
[71]
Event time analysis of longitudinal neuroimage data , author=. NeuroImage , volume=. 2014 , publisher=
work page 2014
-
[72]
Efficient estimation for functional accelerated failure time models , author=. Biometrika , volume=. 2025 , publisher=
work page 2025
-
[73]
Journal of the American Statistical Association , number=
Semiparametric Joint Modeling for Survival Analysis with Longitudinal Covariates , author=. Journal of the American Statistical Association , number=. 2026 , publisher=
work page 2026
-
[74]
Jian Huang , journal =. Maximum Likelihood Estimation for Proportional Odds Regression Model with Current Status Data , urldate =
-
[75]
Tianxi Cai and Rebecca A. Betensky , journal =. Hazard Regression for Interval-Censored Data with Penalized Spline , urldate =
-
[76]
Wang, Lianming and McMahan, Christopher S. and Hudgens, Michael G. and Qureshi, Zaina P. , title =. Biometrics , volume =. 2016 , month =
work page 2016
-
[77]
A Baseline-free Procedure for Transformation Models Under Interval Censorship , volume =
Gu, Ming Gao and Sun, Liuquan and Zuo, Guoxin , date =. A Baseline-free Procedure for Transformation Models Under Interval Censorship , volume =. Lifetime Data Analysis , number =. doi:10.1007/s10985-005-5235-x , isbn =
-
[78]
Zhigang Zhang and Liuquan Sun and Xingqiu Zhao and Jianguo Sun , journal =. Regression Analysis of Interval-Censored Failure Time Data with Linear Transformation Models , urldate =
-
[79]
Empirical likelihood for linear transformation models with interval-censored failure time data , journal =. 2013 , issn =. doi:https://doi.org/10.1016/j.jmva.2013.01.003 , author =
-
[80]
Journal of the American Statistical Association , volume =
Kun Chen and Kehui Chen and Hans-Georg Müller and Jane-Ling Wang , title =. Journal of the American Statistical Association , volume =. 2011 , publisher =
work page 2011
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.