bshazard-package: Nonparametric Smoothing of the Hazard Function

Description Details Author(s) References Examples

Description

The function estimates the hazard function non parametrically from a survival object (possibly adjusted for covariates). The smoothed estimate is based on B-splines from the perspective of generalized linear mixed models. Left truncated and right censoring data are allowed.

Details

The DESCRIPTION file: This package was not yet installed at build time.

Index: This package was not yet installed at build time.

Author(s)

Paola Rebora, Agus Salim, Marie Reilly Maintainer: Paola Rebora <paola.rebora@unimib.it>

References

Rebora P, Salim A, Reilly M (2014) bshazard: A Flexible Tool for Nonparametric Smoothing of the Hazard Function.The R Journal Vol. 6/2:114-122.

Lee Y, Nelder JA, Pawitan Y (2006). Generalized Linear Models with Random Effects: Unified Analysis via H-likelihood, volume 106. Chapman & Hall/CRC.

Pawitan Y (2001). In All Likelihood: Statistical Modelling and Inference Using Likelihood. Oxford University Press

Examples

1
2
3
data(lung,package="survival")
  fit<-bshazard(Surv(time, status==2) ~ 1,data=lung)
  plot(fit)

Example output

Loading required package: splines
Loading required package: survival
Loading required package: Epi

Attaching package: 'Epi'

The following object is masked from 'package:base':

    merge.data.frame

Iterations: relative error in phi-hat = 1e-04 
phi= 1.227307   sv2= 0.02626932   df= 9.117773   lambda= 46.72016 
phi= 1.243423   sv2= 0.006915531   df= 6.543105   lambda= 179.8016 
phi= 1.262371   sv2= 0.002194549   df= 4.957587   lambda= 575.23 
phi= 1.276969   sv2= 0.0008269603   df= 3.94415   lambda= 1544.172 
phi= 1.285215   sv2= 0.00044136   df= 3.283988   lambda= 2911.944 
phi= 1.289474   sv2= 0.0003237787   df= 2.941378   lambda= 3982.579 
phi= 1.291689   sv2= 0.0002773528   df= 2.793532   lambda= 4657.204 
phi= 1.292859   sv2= 0.0002555827   df= 2.724854   lambda= 5058.475 
phi= 1.293494   sv2= 0.000244403   df= 2.689998   lambda= 5292.465 
phi= 1.293847   sv2= 0.0002383763   df= 2.671347   lambda= 5427.749 
phi= 1.294045   sv2= 0.0002350399   df= 2.661065   lambda= 5505.64 
phi= 1.294158   sv2= 0.0002331652   df= 2.655302   lambda= 5550.389 
phi= 1.294222   sv2= 0.000232103   df= 2.652041   lambda= 5576.066 
phi= 1.294258   sv2= 0.0002314983   df= 2.650186   lambda= 5590.79 
phi= 1.294279   sv2= 0.000231153   df= 2.649128   lambda= 5599.231 
phi= 1.294291   sv2= 0.0002309556   df= 2.648523   lambda= 5604.068 
phi= 1.294298   sv2= 0.0002308427   df= 2.648177   lambda= 5606.84 
phi= 1.294302   sv2= 0.000230778   df= 2.647979   lambda= 5608.428 
phi= 1.294304   sv2= 0.000230741   df= 2.647865   lambda= 5609.338 
phi= 1.294305   sv2= 0.0002307198   df= 2.6478   lambda= 5609.859 

bshazard documentation built on May 2, 2019, 5:58 a.m.