Dear all,
I have a problem with accessing class attributes. I was unable to solve this
yet, but someone may know how to solve it.
I'm trying to extract some information from the summary, and Akaike
weight has the desired value.
Object for a model fitted using the glmmML function from the
glmmML package:
result <- glmmML(cbind (y, n-y)~ x+a+b+c, family = binomial, data,
cluster）
library(MASS)
stepAIC(result)
Then calculated the delta AIC by hand (following is the best four).
model1 <- 0.0
model2 <- 1.8
model3 <- 4.2
model4 <- 6.2
Then followed equation as below:
*W <- exp(-0.5 * Delta) / sum(exp(-0.5 * Delta))*
However, result was always same value as [1] even each delta AIC is
different values.
I don't know why happened.
I've also tried Ben's AIC tab in the bbmle package under the Ben's
suggestion:
> summary <- stepAIC(result)
> AICtab(summary)
When I try to run the code from within a package, error came up
as UseMethod("logLik") no method to use "logLik".
I've tried adding > slot <-summary (result)
> slotName (slot)
but it didn't help, slotName is NULL.
Thanks for any kind of suggestions!
Odette
<r-help@r-project.org>
[[alternative HTML version deleted]]
______________________________________________
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.
