Abstract:PURPOSE: A substantial body of epidemiologic literature addresses risks associated with occupational radiation exposure but comparing results between studies is often difficult as different statistical models are commonly used. It is unclear whether different methods produce similar results for estimates of radiation risk when applied to the same data. The goal of this study was to compare the radiation risk estimates for leukemia other than chronic lymphocytic leukemia (non-CLL) and ischemic heart disease (IHD) produced by both Cox and Poisson regression models for time-dependent dose-response analyses of occupational exposure.MATERIALS AND METHODS: For brevity, this methods paper presents the results from one cohort, the Nuclear Power Plant workers (NPP), though the evaluation considered five cohorts of varying size and exposure as part of the Million Worker Study. Cox Proportional Hazards models, with age as the underlying timescale for hazard, were conducted using three computer software programs: SAS, R, and Epicure. Doses lagged 2 years for non-CLL and 10 years for ischemic heart disease were treated as time-dependent exposures at the annual level and were examined both in categories and as a continuous term. Hazard ratios (HR) and 95% confidence intervals (CI) were reported for each model in SAS and R, while the Peanuts program of Epicure was utilized to produce Excess Relative Risk (ERR) estimates and 95% CI. All models were adjusted for gender and year of birth. Four piece-wise exponential Poisson models (log-linear regression for rate) were developed with varying cutpoints for age strata from very fine to broad categories using both R and the Amfit program in Epicure for ERR estimates.RESULTS: Comparable estimates of risk (both RR and ERR) were observed from Cox and Poisson models, regardless of software utilized, as long as appropriately narrow categories of age were utilized to control the confounding of age in Poisson models. The ERR estimates produced in Epicure tended to agree very closely with the HR or RR estimates, and the statistical software program used had no impact to risk estimates for the same model.CONCLUSIONS: As computational power is no longer the burden today as it has been in the past, the results of this evaluation support the use of the Cox proportional hazards or the ungrouped Poisson approach to analyzing time-dependent dose-response relationships to ensure that maximum control over the confounding of age is achieved in studies of mortality for cohorts occupationally exposed to radiation.

Bootstrap confidence interval estimation in generalized nonlinear models

Bootstrap method-based estimation on the confidence interval for additive interaction in cohort studies

Modified Poisson Regression Model for Data of Prospective Studies with Common Outcomes

Overview of model validation for survival regression model with competing risks using melanoma study data

A comparative study with bootstrap resampling technique to uncover behavior of unconditional hazards and survival functions for gamma and inverse Gaussian frailty models

A Bootstrap Approach for Evaluating Uncertainty in the Number of Groups Identified by Latent Class Growth Models

Estimation of the accelerated failure time frailty model under generalized gamma frailty

Bootstrap inference in functional linear regression models with scalar response under heteroscedasticity

Bootstrapping Inference of Average Treatment Effect in Completely Randomized Experiments with High-Dimensional Covariates

Benchmark Dose Analysis via Nonparametric Regression Modeling

Nonparametric estimation of the incubation time distribution

Evaluation of statistical modeling approaches for epidemiologic studies of low-dose radiation health effects

Shared dosimetry error in epidemiological dose-response analyses

Bootstrapping generalized linear models to accommodate overdispersed count data

Computationally Efficient Estimation for the Generalized Odds Rate Mixture Cure Model with Interval-Censored Data.

On models for the estimation of the excess mortality hazard in case of insufficiently stratified life tables

Smooth Semi‐nonparametric Analysis for Mixture Cure Models and Its Application to Breast Cancer

Learning Gaussian Graphical Models from Correlated Data

Radiotherapy toxicity prediction using knowledge-constrained generalized linear model

Model diagnostics for censored regression via randomized survival probabilities

Novel bootstrap-based ellipsoidal convex model for non-probabilistic reliability-based design optimization with insufficient input data