在这个问题之后,我使用该MuMIN
包根据信息标准进行模型平均。
options(na.action = "na.fail")
创建一个包含所有变量和双向交互的全局模型:
global.model<-lmer(yld.res ~ rain + brk +
onset + wid + (1|state),data=data1,REML="FALSE")
标准化全局模型,因为变量的尺度不同
stdz.model <- standardize(global.model,standardize.y = TRUE)
创建模型的所有可能组合
model.set <- dredge(stdz.model)
根据 deltaAICc<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
。但是,由于这些斜率估计是标准化的,我想知道是否有任何方法可以获得非标准化的斜率?
如果我的问题不清楚,请告诉我。
谢谢