The 'checkJacobian' function can be used to check if the analytic gradient / Jacobian is aligned with the numerically approximated gradient / Jacobian, and the 'checkFisher' function can be used to check if the analytic Hessian is aligned with the numerically approximated Hessian.
checkJacobian(x, f = "default", jac = "default", transpose = FALSE,
plot = TRUE, perturbStart = FALSE, method = "Richardson")checkFisher(x, f = "default", fis = "default", transpose = FALSE,
plot = TRUE, perturbStart = FALSE)
Invisibly returns a list with two elements, analytic and numeric, containing the analytic and the numerically approximated gradient/Jacobian (checkJacobian) or Hessian/Fisher information (checkFisher) values being compared. A scatter plot of the two is drawn when plot = TRUE.
A 'psychonetrics' object
A custom fit function or the psychonetrics default fit function (default).
A custom Jacobian function or the psychonetrics default Jacobian function (default).
A custom Fisher information function or the psychonetrics default Fisher information function (default).
Should the numeric Jacobian be transposed?
Should a diagnostic plot be produced?
Should start values be perturbed (only used in development)
Numeric derivative method (default: Richardson)
Sacha Epskamp