glmmTMB с бета-семейством

Question

Я пытаюсь смоделировать ответ фотосинтеза c (fvfm), который принимает значения от 0 до 0,7, я использую бета-распределение, поскольку у меня нулевая инфляция, и я рассматриваю модель дисперсии, потому что мои значения отходят от среднее значение по дням, но при запуске этой модели:

beta7<-glmmTMB(fvfm~Temperature*dia+(1|ID), zi=~1, disp=~dia, data=ds, family=beta_family(link="logit"))

, и это выглядит так:

Error in eval(family$initialize) : y values must be 0 < y < 1

Я добавляю ко всем значениям моей переменной отклика 0,001, поэтому я больше не иметь нулевую инфляцию и модель работает. Это правильно?? или есть способ запустить без добавления 0,001 ?? Спасибо!! С уважением, Фрэн

Answer 1

Я думаю, что бета с нулевым уровнем инфляции может сработать, если ваш результат включает 0 - последняя версия glmmTMB на CRAN (1.0.0) может соответствовать бета-регрессии с нулевым уровнем инфляции.

Иначе, вы также можете подумайте о «сжатии» вашей переменной ответа, используя normalize() -функцию и установив include_bounds в FALSE.