[R-sig-ME] zero one inflated beta mixed model
bbo|ker @end|ng |rom gm@||@com
Wed Feb 12 15:56:44 CET 2020
At present glmmTMB doesn't do zero-one-inflated betas, only
zero-inflated betas. As far as I know your options are (1) use brms,
(2) squish your 1 values to something slightly less than 1, or (3) do
the hurdle model manually (i.e. fit two separate models, one for the
probability that the response== 1, and another (conditional) model for
the zero-inflated beta distribution applied only to the responses <1).
Others on the list may have other suggestions ... (e.g. does INLA
does zero-one-inflated betas?)
On Wed, Feb 12, 2020 at 8:42 AM Morgane Brachet
<morgane.brachet using hotmail.com> wrote:
> I am writing to you following one of the posts on GitHub (https://github.com/glmmTMB/glmmTMB/issues/355). I am trying to fit proportion data with lots of 0s and a few 1s into a hurdle model using glmmTMB. Is this possible? Would you have any example code please?
> Thank you!
> zero inflation for beta distribution model · Issue #355 · glmmTMB/glmmTMB · GitHub<https://github.com/glmmTMB/glmmTMB/issues/355>
> Hi Ben, Thanks for your reply. To fit hurdle model using glmmTMB, do you have any example code? will it still have the same issue? Actually i thought using the ziformula option is the way to fit hurdle model, but since it shows such errors of inappropriate values, i guess i may misse some options for hurdle model.
> [[alternative HTML version deleted]]
> R-sig-mixed-models using r-project.org mailing list
More information about the R-sig-mixed-models