[R-sig-ME] glmmADMB fails to fit poisson data

[Ben Bolker]

  Thanks Dave. Would you suggest replacing the previous Poisson code
entirely?  (That seems sensible; unfortunately I don't think we have the
test data for the previous [May 2013] example lying around ... )

   Ben

   case 0:   // Poisson
	if (poiss_prob_bound==0) {
	    tmpl =  log_density_poisson(y(_i,1),lambda);
	} else {
          if (cph<5)
	    tmpl = log(e3+exp(log_density_poisson(y(_i,1),lambda)));
          else
	    tmpl = log(e4+exp(log_density_poisson(y(_i,1),lambda)));
	}
      break;

On 16-06-11 10:38 AM, dave fournier wrote:
>   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;