#r #lme4 #mixed-models #poisson
Вопрос:
Я использую predictInterval() функцию из merTools пакета. Моя модель подходит для использования спецификации семейства Пуассона, как показано ниже:
glmer(y ~ (1|key) x, data = dat, family = poisson())
Когда я использую predictInterval() для вычисления интервала прогнозирования, связанного с моей моделью, я получаю следующее предупреждающее сообщение:
Warning message:
Prediction for NLMMs or GLMMs that are not mixed binomial regressions is not tested. Sigma set at 1.
Я понимаю это так, что predictInterval() у этого нет реализации для моделей, соответствующих распределению Пуассона. Поэтому я не доверяю полученному интервалу.
Верна ли моя интерпретация? Я искал подобные проблемы, но ничего не нашел.
Любая помощь будет очень признательна.
Комментарии:
1. Пожалуйста, предоставьте достаточно кода, чтобы другие могли лучше понять или воспроизвести проблему.
2.
metToolsработайте в основном в рамках линейных моделей со смешанным эффектом . Мир становится намного сложнее в обобщенных моделях, так как нет замкнутых решений, и для оценки параметров используются различные приближения, поэтому предупреждение (как вы правильно интерпретируете) заключается в том, что им нельзя доверять. Если вы ищете доверительные интервалы и интервалы прогнозирования, различные методы начальной загрузки или оценки с использованием байесовских альтернатив обеспечивают гораздо более надежные оценки.3. Спасибо за ваш очень вдумчивый вклад.