All,
I am working on a multiple-regression meta-analysis and have too many
alternative models to fit by hand. I am using the "metafor" package in
R, which generates AIC scores among other metrics. I'm using a simple
formula to define these models. For example,
rma(Effect_size,variance, mods=~Myco_type + N.type +total,
method="ML")->mod where Effect_size is the mean response for each
experiment, variance is the SE around the mean, mods are the variables,
and the method used to fit the AIC score is maximum likelihood
estimation. I thought I might be able to use the function glmulti(mod,
level=1) to rank all alternative models using AICc scores. If I fit a
glm in in this way (e.g. glm(Effect_size~Myco_type+N.type+total)->mod),
the glmulti function works beautifully. Is it possible to extract AIC
scores from the rma models as well? Is there another package that
might assimilate AIC scores from many alternative models more easily?
Thank you for any insight,
Meghan