The R Journal: article published in 2014, volume 6:2

bshazard: A Flexible Tool for Nonparametric Smoothing of the Hazard Function
Paola Rebora, Agus Salim and Marie Reilly , The R Journal (2014) 6:2, pages 114-122.

Abstract 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.

Received: 2014-05-25; online 2015-01-09
CRAN packages: muhaz, flexsurv, bshazard, Epi, survival, splines
CRAN Task Views implied by cited CRAN packages: Survival, ClinicalTrials, Distributions, Econometrics, SocialSciences

This article is licensed under a Creative Commons Attribution 3.0 Unported license .

@article{RJ-2014-028,
  author = {Paola Rebora and Agus Salim and Marie Reilly},
  title = {{bshazard: A Flexible Tool for Nonparametric Smoothing of the
          Hazard Function}},
  year = {2014},
  journal = {{The R Journal}},
  doi = {10.32614/RJ-2014-028},
  url = {https://doi.org/10.32614/RJ-2014-028},
  pages = {114--122},
  volume = {6},
  number = {2}
}