[R] Extracting Coefficients and Such from mle2 Output

BooBoo booboo at gforcecable.com
Sun Feb 1 13:50:04 CET 2009


I had asked this question once before about a function in the NADA 
package, and you provided this neat response:

[Begin quote]

An approach that may yield somewhat more self-documenting code would be 
to examine either the fit object or the summary object with str and then 
to access results by extracting named elements. Since I don't have the 
package in question, let me use the lm object on its help page as an 
example:

  ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
  trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
  group <- gl(2,10,20, labels=c("Ctl","Trt"))
  weight <- c(ctl, trt)
  lm.D9 <- lm(weight ~ group)

  str(lm.D9)
  str(summary(lm.D9))
  summary(lm.D9)$coefficients
  summary(lm.D9)$coefficients["groupTrt", "Pr(>|t|)"]

# to get the p-value
 > summary(lm.D9)$coefficients["groupTrt","Pr(>|t|)"]
[1] 0.2490232

[End quote]

However, this does not seem to work with mle2, and the methods you 
suggested this time don't allow me access to the "etc" (the standard 
error and its probability for example). Is there a way to do this 
similar to what you suggested for lm, where I can at anything reported 
by mle2?

Tom


David Winsemius wrote:
> 
> 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.
>




More information about the R-help mailing list