[R-sig-ME] mixed lognormal hurdle model with multiple grouping factors

[Guillaume Adeux]

Hi everyone,

I am looking for a package which can handle "hurdle.lognormal" distribution
family and multiple grouping factors.

GLMMadaptive seemed as the way to go but unfortunately, to the best of my
knowledge, it does not handle multiple grouping factors (random effects).

You may ask why? I am analyzing plant diversity and one of the treatments
led to plots which were dominated by one species. Hence, certain diversity
indices are estimated as zero in these plots, and produces a mass at zero.
All other values are positive and continuous.

Anyone have an idea of a package/function which can handle this? Or any
alternative approach?

In lmer syntax, the model is the following:

mod=lmer(diversity~block+syst+(1|plot)+(1|year)+(1|plot:year),data=density,REML=F)

Thank you for your time and help.

Sincerely,

Guillaume ADEUX

	[[alternative HTML version deleted]]