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.
______________________________________________
R-help@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.
