0

在这个问题之后,我使用该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。但是,由于这些斜率估计是标准化的,我想知道是否有任何方法可以获得非标准化的斜率?

如果我的问题不清楚,请告诉我。

谢谢

4

0 回答 0