Работает ли команда "subset" в функции земснаряда MuMIn со случайными эффектами?
У меня есть следующая GLMM со случайным наклоном, на которой я провожу усреднение модели с использованием MuMIn "драги" (я не думаю, что мой набор данных актуален, поэтому я просто включаю здесь код):
sett_global2 <- glmmTMB(sett_dens~ depth_scale + sst_six_scale + exposure_six_month_daily_log_scale + (1|kelp)+ (0 + sst_six_scale|kelp) + (1|site) + (1|Year), data = sett_model_data, family=nbinom2)
Когда я изначально пытался углубить модель, она работала, но со следующими предупреждениями:
dredge(sett_global2, beta = "none", rank = "AIC", trace = 2)
Предупреждающие сообщения: 1: Случайные уклоны не являются фиксированными эффектами. Это искусственно раздувает условные дисперсии случайных эффектов. Решение: Уточните фиксированную конструкцию! 2: Случайные наклоны не являются фиксированными эффектами. Это искусственно раздувает условные дисперсии случайных эффектов. Решение: Уточните фиксированную конструкцию! 3: Случайные уклоны не являются фиксированными эффектами. Это искусственно увеличивает дисперсию условных случайных эффектов. Решение: Уточните фиксированную конструкцию!
это указывает на то, что даже когда фиксированный эффект "sst_six_scale" не был в модели, подходящей к полному набору моделей, модель все еще соответствовала случайному наклону "(0 + sst_six_scale|kelp)".
Поэтому я снова попытался углубить модель, используя аргумент "dc" в команде "subset", чтобы включить случайный уклон, только если соответствующий фиксированный эффект присутствовал в модели:
dredge(sett_global2, beta = "none", rank = "AIC", subset = dc("cond(sst_six_scale)" , "cond(0 + sst_six_scale | kelp)"), trace = 2)
Но я все еще получаю те же предупреждающие сообщения, что и раньше, указывающие мне, что команда подмножества не работает должным образом. Я убедился, что использовал "getAllTerms", чтобы получить соответствующие имена переменных для ввода в команду подмножества, и я не получаю ошибок при запуске приведенного выше кода, поэтому я действительно не могу придумать ни одной причины, по которой он не работает. Это заставляет меня задаться вопросом, работает ли команда "subset" вообще со случайными эффектами. Может ли кто-нибудь пролить свет на это для меня?