<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=us-ascii">
<META NAME="Generator" CONTENT="MS Exchange Server version 6.5.7650.5">
<TITLE>mceval output inconsistent with mcmc display</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/rtf format -->
<P><FONT SIZE=2 FACE="Arial">mcmc stdout displays the current value of the (negative log) likelihood (twice for some reason) then "mcmc sim x acceptance rate y1 y2" like this:</FONT></P>
<P><FONT SIZE=2 FACE="Arial">757.732553269085 757.732553269085</FONT>
<BR><FONT SIZE=2 FACE="Arial"> mcmc sim 401 acceptance rate 0.0698254364089776 0.14</FONT>
</P>
<P><FONT SIZE=2 FACE="Arial">I've checked the code to make sure it is actually the likelihood being displayed and it seems to be but tell me if I'm wrong. when I use mceval the parameter values are inconsistent with this displayed likelihood (their likelihood is much lower). How could this be happening? by the end of a million iteration run the displayed likelihood was 603 which seems within reason. the corresponding parameter values obtained from mceval were around -700! and needless to say a really terrible fit. </FONT></P>
<P><FONT SIZE=2 FACE="Arial">any idea what could be going wrong? Am I interpreting the displayed output correctly? The current likelihood should be the current *accepted* likelihood, so should translate to valid parameters in the saved to psv and read during mceval.? Code and input files are rather big - I won't spam people with it. I can send directly if anyone volunteers. This is what I'm doing for the mceval section if that sheds any light (I doubt it). </FONT></P>
<P><FONT SIZE=2 FACE="Arial"> if (mceval_phase()) {</FONT>
<BR><FONT SIZE=2 FACE="Arial"> cout << a_sr << " " << b_sr << " " << c_sr << " " << iota << " " << vartheta << " " << vartheta_ce << " " << nexpp << " ";</FONT></P>
<P><FONT SIZE=2 FACE="Arial"> for (int i=1;i<=nyears;i++) {</FONT>
<BR><FONT SIZE=2 FACE="Arial"> cout << xi(i) << " ";</FONT>
<BR><FONT SIZE=2 FACE="Arial"> } </FONT>
<BR><FONT SIZE=2 FACE="Arial"> cout << log_sigma_xi << " ";</FONT>
<BR><FONT SIZE=2 FACE="Arial"> for (int i=1;i<=(nyears_spatial*ncells);i++) {</FONT>
<BR><FONT SIZE=2 FACE="Arial"> cout << psi(i) << " ";</FONT>
<BR><FONT SIZE=2 FACE="Arial"> }</FONT>
<BR><FONT SIZE=2 FACE="Arial"> cout << log_gamma << " " << rho_psi << " ";</FONT>
<BR><FONT SIZE=2 FACE="Arial"> cout << log_sigma_psi << " ";</FONT>
<BR><FONT SIZE=2 FACE="Arial"> cout << sigma_C << " " << sigma_E << " ";</FONT>
<BR><FONT SIZE=2 FACE="Arial"> cout << early_kappa << " ";</FONT>
<BR><FONT SIZE=2 FACE="Arial"> //cout << rho_psi_cl << " ";</FONT>
<BR><FONT SIZE=2 FACE="Arial"> //cout << log_sigma_psi_cl << " ";</FONT>
<BR><FONT SIZE=2 FACE="Arial"> cout << E1980 << " " << E1985 << " " << nexpp2 << endl;</FONT>
<BR><FONT SIZE=2 FACE="Arial"> //cout << mu_psi_cl << " ";</FONT>
<BR><FONT SIZE=2 FACE="Arial"> }</FONT>
</P>
<BR>
<P><FONT SIZE=2 FACE="Arial">This has occurred both using admb 10.0 on a windoze box and admb 9.1 (I think) on a bigger, faster, 64-bit windoze box. I've also discovered more linux-windoze inconsistencies but that's another story for another day..</FONT></P>
<P><FONT SIZE=2 FACE="Arial">-alex</FONT>
</P>
<BR>
<BR>
<br><br><table bgcolor=white style="color:black"><tr><td><br>********************************DISCLAIMER****************************<br>
The information contained in the above e-mail message or messages <br>
(which includes any attachments) is confidential and may be legally <br>
privileged. It is intended only for the use of the person or entity <br>
to which it is addressed. If you are not the addressee any form of <br>
disclosure, copying, modification, distribution or any action taken <br>
or omitted in reliance on the information is unauthorised. Opinions <br>
contained in the message(s) do not necessarily reflect the opinions <br>
of the Queensland Government and its authorities. If you received <br>
this communication in error, please notify the sender immediately <br>
and delete it from your computer system network.</td></tr></table></BODY>
</HTML>