[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