I used the MuMIN
package to do a model averaging based on information criterion, following this question.
options(na.action = "na.fail")
Create a global model with all variables and two way interaction:
global.model<-lmer(yld.res ~ rain + brk +
onset + wid + (1|state),data=data1,REML="FALSE")
Standardise the glboal model since the variables are on different scale
stdz.model <- standardize(global.model,standardize.y = TRUE)
Create all possible combinations of the model
model.set <- dredge(stdz.model)
Get the best model based on deltaAICc<2 criteria
top.models <- get.models(model.set, subset= delta<2)
Average the models to calculate the effect size (standardised slopes of the input variables)
s<-model.avg(top.models)
summary(s);confint(s)
The effect size of the variables are as follows:
Variable slope estimate
brk -0.28
rain 0.13
wid 0.10
onset 0.09
As you can see, I had standardize my model in step 3 so I can compare these slope estimates i.e. I can say slope estimate of brk
is greater (in negative direction) than rain
. However, since these slope estimates where standardised, I wanted to know if there is any way I can get the unstandardised slopes?
Please let me know if my question is not clear.
Thanks