R Under development (unstable) (2023-12-07 r85661 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2023 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(bbmle)
Loading required package: stats4
> 
> set.seed(101)
> z = rpois(100,lambda=5)
> 
> m1 = mle2(z~dpois(lambda=L),start=list(L=4),data=data.frame(z))
> q1 <- qAICc(m1,nobs=100,dispersion=1.2)
> qAICc(m1,m1,nobs=100,dispersion=1.2) ## !!
      AICc df
1 374.2214  2
2 374.2214  2
> i1 <- ICtab(m1,type="qAICc",dispersion=1.2,nobs=100, base=TRUE)
> 
> m2 = glm(z~1,family=poisson)
> q2 <- qAICc(m2,nobs=100,dispersion=1.2)
> 
> ## test that dAIC ignores
> m3 <- glm(z~1,family=quasipoisson)
> aa <- AICtab(m1,m2,m3,weights=TRUE)
> stopifnot(any(!is.na(aa$dAIC)),
+           any(!is.na(aa$weight)))
> 
> set.seed(101)
> x <- rnorm(100)
> dd <- data.frame(y=rnorm(100,2+3*x,sd=1),x)
> m4A <- lm(y~x,dd)
> m4B <- mle2(y~dnorm(mean=a+b*x,sd=exp(logsd)),
+            data=dd,
+            start=list(a=1,b=1,logsd=0))
> ## cosmetic differences only
> stopifnot(all.equal(AIC(m4A,m4B)[,"AIC"],
+                     AIC(m4B,m4A)[,"AIC"]))
> 
> 
> proc.time()
   user  system elapsed 
   1.21    0.07    1.29