Uses function hessian() from numDeriv package to compute (approximate numerically) Hessian matrix
evaluated at arbitrary values of item easiness parameters.
Usage
tcl_hessian(data, eta, W, model = "RM")
Value
Hessian matrix evaluated at eta.
Arguments
data
data matrix.
eta
numeric vector of item easiness parameters.
W
design matrix.
model
RM, PCM, RSM, LLTM. Default is set to "RM".
References
Gilbert, P., Gilbert, M. P., & Varadhan, R. (2016). numDeriv: Accurate Numerical Derivatives. R package
version 2016.8-1.1. url: https://CRAN.R-project.org/package=numDeriv
if (FALSE) {
# Rasch model with beta_1 restricted to 0y <- eRm::raschdat1
res <- eRm::RM(X = y, sum0 = FALSE)
mat <- tcl_hessian(data = y, eta = res$etapar, model = "RM")
}