MuMIN包装中的非标准斜坡

时间:2014-11-22 18:38:46

标签: r lme4 standardized mumin

我使用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。但是,由于这些斜率估算标准化,我想知道是否有任何方法可以得到非标准化的斜坡?

如果我的问题不明确,请告诉我。

由于

0 个答案:

没有答案