This version introduces a complete overhaul of the package, with substantial improvements to both the core computational engine and the user-facing interface.
The main change in this version is the exploitation of the banded structure of the penalization matrix, which allows for tremendous savings in both computation time and memory footprint.
-
The package now uses Rcpp and calls directly Lapack functions for banded matrices when available.
-
Some additional functions have also been coded using Rcpp. The
backsolvefunction for banded matrix was also recoded although it already exists in Lapack because the Lapack version overwrites the Cholesky factor so a copy has to be made beforehand. -
Parameter estimation now rely on backward-forward solves from the Cholesky factor instead of forming the variance-covariance matrix explicitly.
-
By default, only the diagonal of the variance-covariance matrix is now computed (which is sufficient for credibility intervals). The full matrix can still be retrieved using the newly introduced
vcov()method if necessary.
-
The two main functions from earlier versions have been merged into a single unified function:
WH(). -
Performance iteration and rank reduction are no longer available, as they no longer provide meaningful benefits with the new implementation. Indeed, the full-rank version is now efficient enough to handle several thousand observation points without optimization tricks.
- Some minor changes have been made to the example datasets to better align the characteristics of the mortality and long-term care datasets
- Greatly increased tolerance from 1e-6 to 1e-2 to remove a lone failing test which only occurs when using MKL BLAS
- Increased tolerance to remove a lone failing test which only occurs when using MKL BLAS
- Made significant changes to the extrapolation methods provided by be
predict.WH_1dandpredict.WH_2dfunctions. It turns out that the formula used for the extrapolation ignored the innovation error caused by the prior on the extrapolated region, resulting in smaller credibility intervals than they should have been. This now has been fixed.
- Removed tests of the form
expect_equal(f(x), f(x))after confirmation they do not work with MKL BLAS
-
Further improved test robustness
-
Added tests of the form
expect_equal(f(x), f(x))for testing purposes
- Further improved tests robustness
-
Simplified computation of the matrix
tUWUby using thecrossprodfunction, which should reduce memory usage (by half) and computation time (slightly) -
Fixed an issue with the
WH_2dplot and the what = "edf" argument -
Improved tests robustness
- Replaced backquotes by normal quotes in the description field of the DESCRIPTION file
- Updated the package description to link it to the paper whose methods are used in the package
- Fixed an issue in the author field of the DESCRIPTION file
- Official candidate for CRAN release
- Added a
NEWS.mdfile to track changes to the package