Survival Regression with Smoothed Error Distribution

Contains, as a main contribution, a function to fit a regression model with possibly right, left or interval censored observations and with the error distribution expressed as a mixture of G-splines. Core part of the computation is done in compiled C++ written using the Scythe Statistical Library Version 0.3.


24/09/2017 (version 2.0)

  • registration of C routines
  • correction of data.frame arithmetic inconsistency

27/07/2015 (version 1.6)

  • NAMESPACE updated by imports from standard packages graphics and stats

05/02/2015 (version 1.5)

  • CITATION file updated (desc changed to meta)


  • small piece of code added to smoothSurvReg function to treat non-zero fail of any type when doing the grid search for optimal lambda
  • confint.smoothSurvReg implemented


  • NAMESPACE added
  • MP.pseudoinv function ( La.eigen replaced by eigen.


  • erf and erfc in gauss.h and gauss.cpp changed to erfAK, erfcAK (there were some conflicts with /usr/include/bits/mathcalls.h reported by CRAN checks)


  • bug fixed in estimTdiff (estimTdiff did not work correctly when there was only 1 covariate in a model)


  • computation of confidence intervals added to estimTdiff
  • print.estimTdiff adjusted to give nicer output when there are no covariates in the model


  • lines fit.logn$failure <- 2; fit.logl$failure <- 2; fit.weib$failure <- 2 in smoothSurvReg.R commented due to changes in survival package in R 2.9.0 (argument 'failure' is no more included in survreg.control)


  • #include added to Scythe_Error.h to let abort() work also with g++ 4.3
  • estimTdiff.smoothSurvReg will no more work when logscale.cov1 and logscale.cov2 are specified


  • new macro AK_PRETTY_FUNCTION defined in Scythe_Error.h to avoid back compatibility problems with some compilers (thanks to B. Ripley)

  • all instances of PRETTY_FUNCTION replaced by AK_PRETTY_FUNCTION in files: Scythe_Matrix.h (8x) Scythe_Matrix_Iterator.h (4x) Scythe_IDE.cpp (20x) Scythe_LA.cpp (11x) Scythe_Math.cpp (7x) Scythe_Matrix.cpp (43x) AKMatrix.cpp (2x)

  • bug in Scythe_Matrix.cpp (in Submatrix operator) fixed

  • operations of the type "int / SCYTHE::Matrix" and "SCYTHE::Matrix / double" which were declared as illegal by Sun's compilers (thanks to B. Ripley) fixed

Reference manual

It appears you don't have a PDF plugin for this browser. You can click here to download the reference manual.


2.3 by Arnošt Komárek, a year ago

Browse source code at

Authors: Arnošt Komárek <[email protected]>

Documentation:   PDF Manual  

Task views: Survival Analysis

GPL (>= 2) license

Imports graphics, stats

Depends on survival

Imported by icensBKL.

Depended on by bayesSurv.

See at CRAN