[R-sig-ME] varIdent in glmmTMB?
Ben Bolker
bbo|ker @end|ng |rom gm@||@com
Tue Jan 14 15:18:36 CET 2025
The glmmTMB analogue of varIdent(form = ~ Temperature) would be
dispformula = ~ Temperature
or (if you want the variance model to be parameterized as a
log-dispersion for each temperature rather than a log-dispersion for the
first temperature level plus differences in log-dispersion)
dispformula = ~ 0 + Temperature
cheers
Ben Bolker
On 2025-01-14 7:14 a.m., Davide Nespoli wrote:
> Dear development team,
>
> this is Davide Nespoli, researcher at the University of Florence (Italy).
> After building a zero-inflated negative binomial GLMM of the form:
>
> m<-glmmTMB(Time~Temperature*Sex+Position+(1|Individual),family=nbinom1(),
> ziformula=~1, data=Perm)
>
> I have a good model, but with a bit of heteroschedasticity due to the
> different levels of Temperature (which is a 3-level factor).
> I am wondering if it is possible to integrate a varIdent covariance
> structure in the model.
> I have already tried
> "+(1|Temperature), dispformula=~0,"
> but it doesn't work (neither for ZINB nor for ordered beta on proportional
> data).
> Is there a way to integrate it in such a model in glmmTMB?
>
> Thank you very much.
> Best regards,
>
>
>
> --
> Davide Nespoli, PhD
> Department of Biology
> University of Florence
> via Madonna del Piano 6, I-50019
> Sesto Fiorentino (FI), Italy
> email: davide.nespoli using unifi.it
> Tel. +390554574722
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-mixed-models using r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
More information about the R-sig-mixed-models
mailing list