[R-sig-ME] glmmADMB fails to fit poisson data
dave fournier
davef at otter-rsch.com
Sat Jun 11 16:38:19 CEST 2016
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,
More information about the R-sig-mixed-models
mailing list