Hi Guillaume,
Thank you so much for this! I just have another question: for example if I
have random factor A and B in both logistic model part and beta model part,
then after I fit the whole model and got variance component estimation of
random effect for factor A and B for both logistic model part and beta model
model part, will there be any way to combine variance together? I.e. I can
estimate a total variance from factor A, and a total variance from factor B
(i.e. only differ by factor, not model)? Something like variance
decomposition but I believe here is more complex as this is a mixture model.
Thank you again for all your help
Best regards,
Meng
On Sun, Jun 10, 2018 at 11:03 AM, Guillaume Chaumet
<guillaumechaumet at gmail.com> wrote:
To whom it may concern,
I am trying to fit a model for a data among which the response value is
within [0,1). I am thinking about fitting the zeros as a complete
separate
category from the non-zero data, i.e. a binomial (Bernoulli) model to
"==0
vs >0" and a Beta model to the >0 responses. Also, my data contains both
nested factors and crossed factors, which means I need to add nested
random
effects and crossed random effects to both logistic model part and beta
model model. However, I didn't find any R packages can do exactly what I
want (By far I found gamlss, glmmTMB, zoib but they either can only
assume
random zero or they can only fit repeated measures/clustered data but
not
nested and crossed design). Therefore, I am wondering if any one know if
there is any available package or function can do this.
Thank you very much for your help!
Best regards
Meng
[[alternative HTML version deleted]]