R help - Feb 2009 - Extracting Coefficients and Such from mle2 Output

The mle2 function (bbmle library) gives an example something like the
following in its help page. How do I access the coefficients, standard
errors, etc in the summary of "a"?
 
> x <- 0:10
> y <- c(26, 17, 13, 12, 20, 5, 9, 8, 5, 4, 8)
> LL <- function(ymax=15, xhalf=6)
+   -sum(stats::dpois(y, lambda=ymax/(1+x/xhalf),
log=TRUE))
> a <- mle2(LL, fixed=list(xhalf=6))
> summary(a)
Maximum likelihood estimation

Call:
mle2(minuslogl = LL, fixed = list(xhalf = 6))

Coefficients:
     Estimate Std. Error z value     Pr(z)    
ymax  19.2881     1.7115  11.269 < 2.2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05
'.' 0.1 ' ' 1

-2 log L: 60.39244 


Tom


-- 
View this message in context:
http://www.nabble.com/Extracting-Coefficients-and-Such-from-mle2-Output-tp21770860p21770860.html
Sent from the R help mailing list archive at Nabble.com.

On Jan 31, 2009, at 9:13 PM, Tom La Bone wrote:
>
> The mle2 function (bbmle library) gives an example something like the
> following in its help page. How do I access the coefficients, standard
> errors, etc in the summary of "a"?
?coef
?vcov

eeep. Further comment on "etc" not possible at this time. Mind reading
equipment malfunction.

-- 
David Winsemius

>
>
>> x <- 0:10
>> y <- c(26, 17, 13, 12, 20, 5, 9, 8, 5, 4, 8)
>> LL <- function(ymax=15, xhalf=6)
> +   -sum(stats::dpois(y, lambda=ymax/(1+x/xhalf), log=TRUE))
>> a <- mle2(LL, fixed=list(xhalf=6))
>> summary(a)
> Maximum likelihood estimation
>
> Call:
> mle2(minuslogl = LL, fixed = list(xhalf = 6))
>
> Coefficients:
>     Estimate Std. Error z value     Pr(z)
> ymax  19.2881     1.7115  11.269 < 2.2e-16 ***
> ---
> Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05
'.' 0.1 ' ' 1
>
> -2 log L: 60.39244
>
>
> Tom
>
>
> -- 
> View this message in context:
http://www.nabble.com/Extracting-Coefficients-and-Such-from-mle2-Output-tp21770860p21770860.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

Tom La Bone <booboo <at> gforcecable.com> writes:
> 
> 
> The mle2 function (bbmle library) gives an example something like the
> following in its help page. How do I access the coefficients, standard
> errors, etc in the summary of "a"?
> 
example(mle2)
coef(summary(fit1))

  Ben Bolker