Я пытаюсь смоделировать данные подсчета птиц в лесных обломках разного размера.Поскольку графики, на которых проводились опросы, также различаются по размеру между фрагментами, я хотел бы добавить размер графика опроса в качестве смещенного члена, чтобы преобразовать счет в плотности.
Как я понял из предыдущих вопросов на этом сайте, это обычно делается для моделей Пуассона, поскольку они имеют ссылку на журнал.Модель GAM (пакет mgcv), в которой я работаю с семейством ziP, имеет link = "identity".Насколько я понимаю, в таких случаях срок смещения будет вычитаться из ответа, а не в результате в желаемой скорости ответа / смещения.
Однако, когда я запускаю модель со смещенным членом и строю график, результаты, похоже, дают желаемый результат (я сравнил график для модели Пуассона с моделью ziP).
Это модель, которую я использовал, в которой Гильдия отражает различные гильдии кормления, logArea - это журнал размера фрагмента, а Study - мой случайный эффект (данные получены из нескольких исследований).
gam1 <- gam(Count ~ Guild + s(logArea, by=Guild) + s(Study,bs="re"), offset=lnTotalPlotsize, family=ziP(),data=Data_ommited2)
Может кто-нибудь объяснить, как GAM обрабатывает смещенные термины в этом случае (модель ziP с идентификационной связью)?Это действительно приводит к желаемому ответу / скорости смещения, или это делает что-то еще?
Спасибо за вашу помощь!
С уважением, Роберт