The hazard function is a key component in the inferential process in survival analysis and relevant for describing the pattern of failures. However, it is rarely shown in research papers due to the difficulties in nonparametric estimation. We developed the bshazard package to facilitate the computation of a nonparametric estimate of the hazard function, with data-driven smoothing. The method accounts for left truncation, right censoring and possible covariates. B-splines are used to estimate the shape of the hazard within the generalized linear mixed models framework. Smoothness is controlled by imposing an autoregressive structure on the baseline hazard coefficients. This perspective allows an ’automatic’ smoothing by avoiding the need to choose the smoothing parameter, which is estimated from the data as a dispersion parameter. A simulation study demonstrates the capability of our software and an application to estimate the hazard of Non-Hodgkin’s lymphoma in Swedish population data shows its potential.
muhaz, flexsurv, bshazard, Epi, survival, splines
Survival, ClinicalTrials, Distributions, Econometrics, SocialSciences
Text and figures are licensed under Creative Commons Attribution CC BY 4.0. The figures that have been reused from other sources don't fall under this license and can be recognized by a note in their caption: "Figure from ...".
For attribution, please cite this work as
Rebora, et al., "The R Journal: bshazard: A Flexible Tool for Nonparametric Smoothing of the Hazard Function", The R Journal, 2015
BibTeX citation
@article{RJ-2014-028,
author = {Rebora, Paola and Salim, Agus and Reilly, Marie},
title = {The R Journal: bshazard: A Flexible Tool for Nonparametric Smoothing of the Hazard Function},
journal = {The R Journal},
year = {2015},
note = {https://doi.org/10.32614/RJ-2014-028},
doi = {10.32614/RJ-2014-028},
volume = {6},
issue = {2},
issn = {2073-4859},
pages = {114-122}
}