使用glmulti在glmm中进行预测

时间:2018-10-29 10:48:58

标签: r prediction lme4 glmulti

即使在每分钟扫描速率的预测值(spm,连续变量)的情况下,我也需要测试区域[ha]的大小是否起作用。此外,我需要浏览一下。 我的解决方案是使用其他变量构建新的数据框架,并应使用预测的spm对其进行测试。 在这里,我的原木将原木转换为领土尺寸,例如觅食(连续变量):

terri_spm<-lmer(logterrisize~spm + (1|studyarea/teriid),
        data = Data_table_for_analysis_Character_studyarea,
        control=lmerControl(optimizer="Nelder_Mead",
                             optCtrl=list(maxfun=1e4)))

forage_spm<-lmer(foraging~spm + (1|studyarea/teriid),
         data = Data_table_for_analysis_Character_studyarea, REML = TRUE,
         control = lmerControl(optimizer="bobyqa",optCtrl=list(maxfun=1000)))

和预测公式

t.size<-predict(terri_spm, newdata = Data_table_for_analysis_Character_studyarea, newparams = NULL,
    re.form = NULL, terms = NULL,
    type = c("response"), allow.new.levels = FALSE,
    na.action = na.pass)

forage<-predict(forage_spm, newdata = Data_table_for_analysis_Character_studyarea, newparams = NULL,
                   re.form=NULL,
                   type = c("response"), allow.new.levels = FALSE,
                   na.action = na.pass)

新的预测值数据框:

all_preds<-data.frame(t.size,forage)

然后我需要运行glmulti:

behave<-glmulti(logterrisize~forage+t.size,data=all_preds,crit="aicc")

我的问题是:在使用glmer之前或之后,我是否必须将spm转换为预测值?我不确定,是否选择了正确的序列来使用预测的spm进行glmer和glmulti。

0 个答案:

没有答案