Для доступа к коэффициентам используйте функцию coef
, см. ?bbmle::mle2
для следующего примера:
x <- 0:10
y <- c(26, 17, 13, 12, 20, 5, 9, 8, 5, 4, 8)
d <- data.frame(x,y)
LL <- function(ymax=15, xhalf=6)
-sum(stats::dpois(y, lambda=ymax/(1+x/xhalf), log=TRUE))
fit1F <- mle2(LL, fixed=list(xhalf=6))
Получите первый коэффициент:
coef(fit1F)[1]
Или используйте @
:
fit1F@fullcoef[1]