www.phidot.org

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

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