我使用MuMIN
包根据this question之后的信息标准进行模型平均。
options(na.action = "na.fail")
创建一个包含所有变量和双向交互的全局模型:
global.model<-lmer(yld.res ~ rain + brk +
onset + wid + (1|state),data=data1,REML="FALSE")
标准化glboal模型,因为变量的规模不同
stdz.model <- standardize(global.model,standardize.y = TRUE)
创建模型的所有可能组合
model.set <- dredge(stdz.model)
根据deltaAICc&lt; 2条件
获取最佳模型top.models <- get.models(model.set, subset= delta<2)
平均模型以计算效果大小(输入变量的标准化斜率)
s<-model.avg(top.models)
summary(s);confint(s)
变量的效果大小如下:
Variable slope estimate
brk -0.28
rain 0.13
wid 0.10
onset 0.09
正如您所看到的,我已经在步骤3中对我的模型进行了标准化,因此我可以比较这些斜率估计值,即我可以说brk
的斜率估计值(负方向)大于rain
。但是,由于这些斜率估算标准化,我想知道是否有任何方法可以得到非标准化的斜坡?
如果我的问题不明确,请告诉我。
由于