Message-ID: <575C225B.9040703@otter-rsch.com>
Date: 2016-06-11T14:38:19Z
From: dave fournier
Subject: glmmADMB fails to fit poisson data
In-Reply-To: <575B5B57.2030101@gmail.com>
On the other hand maybe the code never will write itself. The following
change to the Poisson
model in glmmadmb using the log(1_+y_i) transformation to get good
initial parameter estimates
converges nicely on this data set.
case 0: // Poisson
if (cph<2)
tmpl=-square(log(1.0+y(_i,1))-log(1.0+lambda));
else
tmpl= log_density_poisson(y(_i,1),lambda);
break;
case 1: // Binomial: y(_i,1)=#successes, y(_i,2)=#failures,