Estimates Hessian of a scalar-valued function, and returns it
in a sparse Matrix format. The sparsity pattern must be known in advance. The
algorithm is especially efficient for hierarchical models with a large number of
heterogeneous units. See Braun, M. (2017)
Added implementation of the complex step method.
The 'direct' argument to the sparseHessianFD initializer was removed (defunct).
An even more major rewrite of the package. All ACM code was removed, and replaced with original R/C++ implementations.
The sparseHessianFD class is now implemented as an R reference
class, and not as an Rcpp module. The
function is deprecated. Instead, use
initialize an object. Initialization once again takes place in a single step.
The 'direct' computation method has been removed. All computation uses the 'indirect' triangular substitution method. The 'direct' argument in the initializer for the sparseHessianFD class is now deprecated, and remains solely for compatibility with older versions of the package.
There is a new vignette with a lot more detail about what the package does and how it works.
New matrix helper functions
Coord.to.Pattern.Matrix function is deprecated.
sparseMatrix functions in the Matrix
With the removal of ACM-copyrighted code, this package is now licensed under the MPL 2.0.
Essentially a complete re-write of the package.
New vignette, using a binary choice model as an example. Functions for this model are in binary.R. Access sample simluated data with data(binary).
Documentation written using roxygen2
Added unit tests using testthat
Core class has been renamed sparseHessianFD. Construction and initialization are now two separate steps.
New function sparseHessianFD.new is a wrapper around the construction and initialization steps.
Deprecated functions in Version 0.2.0
Coord.to.Pattern.Matrixinstead, with option
Sym.CSC.to.Matrix. Use the
spMatrixfunction in the Matrix package instead.