Table of Contents
Fetching ...

Semiparametric estimation of GLMs with interval-censored covariates via an augmented Turnbull estimator

Andrea Toloba, Klaus Langohr, Guadalupe Gómez Melis

TL;DR

This paper addresses regression with interval-censored covariates in generalized linear models by introducing GELc, a likelihood-based estimator that augments Turnbull's nonparametric estimator for the censored covariate distribution. GELc jointly estimates the GLM parameters and a nonparametric covariate distribution using augmented Turnbull intervals and self-consistent updating, with the authors proving consistency and asymptotic normality under mild regularity. The method is validated via extensive simulations showing favorable finite-sample performance and good coverage, and demonstrated on two real datasets (carotenoids and ACTG359) to illustrate practical applicability; an R package ICenCov provides implementation. Overall, GELc offers robust inference for GLMs with interval-censored covariates by avoiding parametric misspecification of the covariate distribution while providing standard errors and feasible computation.

Abstract

Interval-censored covariates are frequently encountered in biomedical studies, particularly in time-to-event data or when measurements are subject to detection or quantification limits. Yet, the estimation of regression models with interval-censored covariates remains methodologically underdeveloped. In this article, we address the estimation of generalized linear models when one covariate is subject to interval censoring. We propose a likelihood-based approach, GELc, that builds upon an augmented version of Turnbull's nonparametric estimator for interval-censored data. We prove that the GELc estimator is consistent and asymptotically normal under mild regularity conditions, with available standard errors. Simulation studies demonstrate favorable finite-sample performance of the estimator and satisfactory coverage of the confidence intervals. Finally, we illustrate the method using two real-world applications: the AIDS Clinical Trials Group Study 359 and an observational nutrition study on circulating carotenoids. The proposed methodology is available as an R package at github.com/atoloba/ICenCov.

Semiparametric estimation of GLMs with interval-censored covariates via an augmented Turnbull estimator

TL;DR

This paper addresses regression with interval-censored covariates in generalized linear models by introducing GELc, a likelihood-based estimator that augments Turnbull's nonparametric estimator for the censored covariate distribution. GELc jointly estimates the GLM parameters and a nonparametric covariate distribution using augmented Turnbull intervals and self-consistent updating, with the authors proving consistency and asymptotic normality under mild regularity. The method is validated via extensive simulations showing favorable finite-sample performance and good coverage, and demonstrated on two real datasets (carotenoids and ACTG359) to illustrate practical applicability; an R package ICenCov provides implementation. Overall, GELc offers robust inference for GLMs with interval-censored covariates by avoiding parametric misspecification of the covariate distribution while providing standard errors and feasible computation.

Abstract

Interval-censored covariates are frequently encountered in biomedical studies, particularly in time-to-event data or when measurements are subject to detection or quantification limits. Yet, the estimation of regression models with interval-censored covariates remains methodologically underdeveloped. In this article, we address the estimation of generalized linear models when one covariate is subject to interval censoring. We propose a likelihood-based approach, GELc, that builds upon an augmented version of Turnbull's nonparametric estimator for interval-censored data. We prove that the GELc estimator is consistent and asymptotically normal under mild regularity conditions, with available standard errors. Simulation studies demonstrate favorable finite-sample performance of the estimator and satisfactory coverage of the confidence intervals. Finally, we illustrate the method using two real-world applications: the AIDS Clinical Trials Group Study 359 and an observational nutrition study on circulating carotenoids. The proposed methodology is available as an R package at github.com/atoloba/ICenCov.
Paper Structure (26 sections, 7 theorems, 105 equations, 9 figures, 20 tables, 1 algorithm)

This paper contains 26 sections, 7 theorems, 105 equations, 9 figures, 20 tables, 1 algorithm.

Key Result

Lemma A.1

Consider $W_1, W_2$ two possible distribution functions for $Z$, and let $a,b$ be two arbitrary constants. Then, $L_i(aW_1+bW_2) = a L_i(W_1) + b L_i(W_2)$.

Figures (9)

  • Figure 1: Simulated observed intervals for $n=300$ under the four censoring scenarios: the true $Z_i$ values ($\mu=0$) and the three scenarios of censoring process with increasing mean widths. Observations are ordered by the true $Z_i$ values to facilitate comparison of how censoring interval widths vary across scenarios.
  • Figure 2: Results for the gamma regression with true parameters $\alpha=10$, $\gamma = -0.05$ and $\phi = 1$, showing relative bias (with Monte Carlo error bars) and empirical standard error. From left to right: $\alpha$, $\gamma$, and $\phi$ parameters.
  • Figure B.3: Relative bias of $\gamma$ for increasing mean width of observed intervals (x-axis) and sample size, in the Gamma regression model for varying $\gamma_0$ and $\phi_0$.
  • Figure B.4: Relative bias of $\gamma$ for increasing mean width of observed intervals (x-axis) and sample size, in the Logistic regression model for varying $\gamma_0$.
  • Figure B.5: Relative bias of $\alpha$ for increasing mean width of observed intervals (x-axis) and sample size, in the Gamma regression model for $\alpha_0=10$, and varying $\gamma_0$ and $\phi_0$.
  • ...and 4 more figures

Theorems & Definitions (15)

  • proof
  • Lemma A.1
  • proof
  • Corollary A.2
  • proof
  • Lemma A.3
  • proof
  • Lemma A.4
  • proof
  • Lemma A.5
  • ...and 5 more