See <a href="http://www.admb-project.org/svn/trunk/src/nh99/mod_hess.cpp">http://www.admb-project.org/svn/trunk/src/nh99/mod_hess.cpp</a><div><br></div><div><pre style="word-wrap:break-word;white-space:pre-wrap">void function_minimizer::hess_inv(void)</pre>
<br><div class="gmail_quote">On Tue, Feb 21, 2012 at 5:32 PM, James Thorson <span dir="ltr"><<a href="mailto:JimThor@uw.edu">JimThor@uw.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi all,<br><br>I've seen a question posted previously (e.g. <a href="http://groups.google.com/group/admb-users/browse_thread/thread/e8ae05a7919a4d19" target="_blank">http://groups.google.com/group/admb-users/browse_thread/thread/e8ae05a7919a4d19</a>), but I've never seen an answer: How does one extract the hessian from the .hes file in ADMB.<br>


<br>I'm attaching two simple examples.<br><br>First simple.tpl and simple.dat are taken from the ADMB examples webpage.  When run, the hessian can be extracted from the .hes in R using the following commands<br>
<br>x<-file('admodel.hes','rb')<br>  nopar<-readBin(x,"integer",1)<br>  H<-matrix(readBin(x,"numeric",nopar*nopar),nopar)<br><br>When inverted, the square-root of the diagonals of the inverted-Hessian are identical to the standard deviations in the STD file.  <br>


<br>However, the simple_bounded.tpl and simple_bounded.dat differ only in having bounds placed on the two parameters.  When run, the hessian in the HES file has no obvious relation to the standard deviations in the STD file.  I assume this is because of some transformation being done on bounded parameters internally to ADMB.  I cannot identify the transformation used by simply exploring the effect of changing the bounds.<br>


<br>It would be useful to extract the hessian from ADMB for two reasons:<br>1.  Using a generalized inverse for a non-positive-definite Hessian as a sampling kernal for Metropolis MCMC a la Gill and King (2004: <a href="http://artsci.wustl.edu/~jgill/papers/gill_king.pdf" target="_blank">http://artsci.wustl.edu/~jgill/papers/gill_king.pdf</a>) either in the REPORT_SECTION or externally<br>


2.  Using ADMB computations of the hessian when using ADMB within a larger model.<br><br>Does anyone know how to extract the "real" hessian from the HES file?  <br><br>Thanks in advance,<br>Jim<br><br><br><br>
<br>_______________________________________________<br>
Users mailing list<br>
<a href="mailto:Users@admb-project.org">Users@admb-project.org</a><br>
<a href="http://lists.admb-project.org/mailman/listinfo/users" target="_blank">http://lists.admb-project.org/mailman/listinfo/users</a><br>
<br></blockquote></div><br></div>