[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