- Fixed C code compatibility issues with
_R_USE_STRICT_R_HEADERS_
(upcoming default in R 4.5.0) (#13):- Defined
STRICT_R_HEADERS
. - Replaced
Calloc
andFree
withR_Calloc
andR_Free
. - Included
float.h
forFLT_EPSILON
.
- Defined
- Fix "lost braces" check notes on r-devel.
- Export S3 method per guidelines in roxygen2 7.3.0.
- Use GitHub Actions to build the pkgdown site.
- Fix code linting issues.
- Use all samples instead of a much smaller set of samples in some code examples
as a temporary fix for reverse dependency check errors which result in the
null model. A proper fix will involve updating the logic to use
ncvsurv()
in a pathwise manner instead of with a single value of lambda.
- Removed the dependency on
survAUC
. Ported the essential C code for computing time-dependent AUC and fixed the build issues in r-devel.
This version is a major refactor of the package, with several technical adjustments to improve the functional interface, code structure, and execution performance. As a result, a few critical API-breaking changes have been made. Please update your previous code that calls hdnom accordingly. For the detailed changes, please check the updated items below.
- Renamed exported functions. Most of the exported function have been renamed to be more meaningful and succinct. For example,
hdcox.*()
are renamed asfit_*()
,hdnom.nomogram()
is renamed asas_nomogram()
,hdnom.validate()
is renamed asvalidate()
, and so on. - Removed the dependency on
rms
, by reusing a minimal set of code fromrms
for nomogram construction and plotting. This results in clearer code structure, better maintainability, and faster package installation/loading speed. Also removed other non-essential package dependencies. - The first argument for
print
functions are now returned invisbily, to make it easier to use them in a pipe.
- The components in the model fitting function returns are now unified across model types. For example, the model object can all be accessed by
fit$model
, and the selected "optimal" hyperparameters can be accessed byfit$lambda
. The model type is now stored explicitly asfit$type
.
as_nomogram
(previouslyhdnom.nomogram()
) now accepts the fitted model objects directly instead of the$model
component. It now will recognize the model type automatically, thus the previous argumentsmodel.type
has been deprecated. so that it is easier to chain the function calls together usingmagrittr
.- In
as_nomogram
, the previousddist
argument is not needed anymore and has been removed. There is also no more need to set adatadist
object as a into the global options variable (which was required in therms
user flow). - The new nomogram implementation prints and plots the nomogram for the penalized regression models directly. This supersedes the old implementation, which fits an OLS model to regress the linear predictors on the same set of predictors selected by the penalized Cox regression model, aiming to approximate the penalized model. The numerical or visual difference is minimal, though.
- Add a new ggplot2 theme
theme_hdnom()
and applies it to most of the validation, calibration, and comparison plots for a consistent, cleaner look across plots within the package.
- Exported the survival curve prediction functions (
glmnet.survcurve()
,ncvreg.survcurve()
,penalized.survcurve()
) and Breslow baseline hazard estimator functions (glmnet.basesurv()
,ncvreg.basesurv()
,penalized.basesurv()
). - New URL for the documentation website: https://nanx.me/hdnom/.
- Use system font stack instead of Google Fonts in vignettes to avoid pandoc SSL issue.
- Reduced example running time to less than 10s for
hdnom.calibrate()
.
- Better code indentation style.
- Reduced example running time for fused lasso.
- Updated gallery images in
README.md
. - HTTPS enabled for the website.
- Fixed issues in parameter tuning and cross-validation procedures for
fused lasso models (nanxstats/hdnom@afc49c9).
The user-visible change is that two parameters
lambda1
andlambda2
instead of a single "lambda" are now required to fit, validate, and calibrate fused lasso models.
- The argument
lambda
inhdnom.nomogram
is no longer needed and has been deprecated. - Allow users to specify
eps
andmax.iter
for MCP and SCAD penalty related models. Setting the default values to be1e-4
and10000
, which is consistent with ncvreg 3.8-0.
- Fixed vanishing axis problem in Kaplan-Meier plot
hdnom.kmplot()
under ggplot2 2.2.0, which is caused by a previous workaround for a bug introduced in ggplot2 2.1.0. - Fixed potential convergence issues for examples under ncvreg >= 3.7-0 new
convergence criterion, by increasing
max.iter
forncvsurv
to a substantially higher value (5e+4). - Fixed single lambda support issues in
ncvsurv
under ncvreg >= 3.7-0.
- Added Windows continuous integration using AppVeyor.
- New website design for project website: consistent with the web application.
- More concrete examples for several functions.
- Introduce argument
ylim
forplot.hdnom.validate()
,plot.hdnom.external.validate()
, andplot.hdnom.compare.validate()
(#4).
- Removed one redundant color from the lancet color palette.
- Added 4 new color palettes (JCO, Lancet, NPG, AAAS) for all plots. See the vignette for details.
- Fixed vanishing axis problem in Kaplan-Meier plot due to ggplot2 2.1.0 update.
- New CSS style for the HTML vignette.
- New function
hdnom.compare.validate()
for model comparison by validation - New function
hdnom.compare.calibrate()
for model comparison by calibration - New function
hdnom.external.validate()
for external validation - New function
hdnom.external.calibrate()
for external calibration - New
predict
andprint
methods forhdcox.model
objects - New function
hdnom.kmplot()
: Kaplan-Meier analysis for risk groups using internal/external calibration results - New function
hdnom.logrank()
: Log-rank test for risk groups using internal/external calibration results
- The web application is substantially improved to reflect the new features
- Record random seeds in the generated reports to improve reproducibility
- Allow users to download the model objects for further exploration in
R
- Improvements on random seed handling for internal validation and calibration
- The vignette is extended to reflect the new features
- Fixed an error in internal calibration which mistakenly used testing data when should use training data in computation
- New documentation website.
- New Shiny web application.
- Added exception handling for null models in all
hdcox.*()
functions. Make examples compatible with ncvreg 3.5-0, which refined CV implementation for survival models and improved computation speed.
-
Support five more high-dimensional penalized Cox model types:
- Fused lasso
- MCP
- Mnet
- SCAD
- Snet
- Reduced example running time for
hdnom.validate()
,hdnom.calibrate()
,hdcox.aenet()
, andhdcox.enet()
by reducing resampling times.
- Added argument
parallel
tohdcox.aenet()
andhdcox.enet()
to enable or disable the use of parallel parameter tuning.
- Initial version
- Nomograms for glmnet models
- Validation for glmnet models
- Calibration for glmnet models