search

bbolker / bbmle

maximum likelihood estimation package
GNU General Public License v3.0
25 stars 13 forks source link

serious qAIC glitches #22

Open bbolker opened 5 years ago

bbolker commented 5 years ago

for example, consider

library(bbmle)

set.seed(101)
z = rpois(100,lambda=5)

m1 = mle2(z~dpois(lambda=L),start=list(L=4),data=data.frame(z))
qAICc(m1,nobs=100,dispersion=1.2)
qAICc(m1,m1,nobs=100,dispersion=1.2) ## !!