https://github.com/cran/bbmle
Tip revision: 526bb859fca17d9fc7d6640c84793e3d8254a248 authored by Ben Bolker on 22 May 2010, 05:22:54 UTC
version 0.9.5.1
version 0.9.5.1
Tip revision: 526bb85
BIC.R
require(bbmle)
x <- 0:10
y <- c(26, 17, 13, 12, 20, 5, 9, 8, 5, 4, 8)
fit <- mle2(y~dpois(lambda=ymax/(1+x/xhalf)), start=list(ymax=25,xhalf=3))
fit2 <- mle2(y~dpois(lambda=(x+1)*slope), start=list(slope=1))
BIC(fit,nobs=length(x))
BIC(fit,fit2,nobs=length(x))
