<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
<link href="chrome://translator/skin/floatingPanel.css"
type="text/css" rel="stylesheet">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Hi all<br>
<br>
I'm not sure if this is a problem with glmmabmb or the coda package,
or (most likely) my code. I ran this model:<br>
<br>
fit1.mcmc <- glmmadmb(Y ~ time*X1 + time*X2 + ( 1| Subject),
data=final, family="nbinom2", zeroInflation=TRUE,<br>
mcmc=TRUE,mcmc.opts=mcmcControl(mcmc=10000))<br>
<br>
and summary(zi1.mcmc) gave me the estimates I was expecting, that
is:<br>
<br>
Estimate Std. Error z value Pr(>|z|)
<br>
(Intercept) .......<br>
time ......<br>
X1 ......<br>
X2 ......<br>
time:X1 ......<br>
time:X2 ......<br>
<br>
However, when running head(HPDinterval(as.mcmc(fit1.mcmc$mcmc))) I
got<br>
<br>
lower upper<br>
pz ........... ...........<br>
(Intercept) ........... ...........<br>
time ........... ...........<br>
X1 ........... ...........<br>
X2 ........... ...........<br>
time:X1 ........... ...........<br>
<br>
The interval for time:X2 is absent.<br>
<br>
Am I missing something (excuse the pun) ?<br>
<br>
Thanks<br>
WRL<br>
<br>
<br>
<br>
<br>
<br>
<div style="bottom: auto; left: 89px; right: auto; top: 324px;
display: none;" class="translator-theme-default"
id="translator-floating-panel">
<div title="Click to translate"
id="translator-floating-panel-button"></div>
</div>
</body>
</html>