Evaluating the Hessian

[tam10]

Hi, is it possible to retrieve the Hessian from the output in RMark? I'm running simulations and would like to simply evaluate whether the Hessian is positive definite.

many thanks!
Tom

[jlaake]

The beta var-cov matrix (inverse of hessian) is stored under results. Here is an example

data(dipper)
mod=mark(dipper)
mod$results$beta.vcv
# hessian
solve(mod$results$beta.vcv)

--jeff