即使在每分钟扫描速率的预测值(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。