предсказание в glmm с glmulti - PullRequest
Новая
       34

предсказание в glmm с glmulti

0 голосов
/ 29 октября 2018

Мне нужно проверить, играет ли роль размер территории [га], даже в случае прогнозируемых значений скорости сканирования в минуту (spm, непрерывная переменная). Кроме того, мне нужно запустить блеск. Мое решение состояло в том, чтобы создать новый фрейм данных с дополнительными переменными, которые должны быть проверены с предсказанным 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")

Мой вопрос здесь: нужно ли преобразовывать spm в прогнозируемое значение перед использованием glmer или после него? Я не уверен, выбрал ли я правильную последовательность, чтобы сделать glmer и glmulti с предсказанным spm.

...