Hello, 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! Morgane [https://avatars1.githubusercontent.com/u/13640228?s=400&v=4]<https://github.com/glmmTMB/glmmTMB/issues/355> 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. github.com
zero one inflated beta mixed model
2 messages · Morgane Brachet, Ben Bolker
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 at hotmail.com> wrote:
Hello, 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! Morgane [https://avatars1.githubusercontent.com/u/13640228?s=400&v=4]<https://github.com/glmmTMB/glmmTMB/issues/355> 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. github.com [[alternative HTML version deleted]]
_______________________________________________ R-sig-mixed-models at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models