https://github.com/cran/bbmle
Tip revision: 07cd713b4f52d0cde5dbb27602e2432d222df3a8 authored by Ben Bolker on 29 July 2012, 00:00:00 UTC
version 1.0.5.2
version 1.0.5.2
Tip revision: 07cd713
optimize.R
## try to reconstruct error reported by Hofert Jan Marius
## (simpler version)
Lfun <- function(x) {
(x-5)^2
}
library(bbmle)
lb <- 6
## first try with L-BFGS-B and bounds
m1 <- mle2(Lfun,start=list(x=7),lower=6,method="L-BFGS-B")
coef(m1)
p1 <- profile(m1)
plot(p1)
(c1 <- confint(m1,quietly=TRUE))
## all OK
m2 <- mle2(Lfun,start=list(x=7),optimizer="optimize",
lower=lb,upper=10)
coef(m2)
p2 <- profile(m2)
(c2 <- confint(m2))
(c2 <- confint(m2))
plot(p2,show.points=TRUE)