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.
David Winsemius
2009-Feb-01 04:11 UTC
[R] Extracting Coefficients and Such from mle2 Output
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