# NOT RUN {
library(bbmle)
x <- 0:10
y <- c(26, 17, 13, 12, 20, 5, 9, 8, 5, 4, 8)
d <- data.frame(x, y)
fit <- mle2(y ~ dpois(lambda = ymean),
start = list(ymean = mean(y)), data = d
)
tidy(fit)
# }
Run the code above in your browser using DataLab