Phylogenetic Dirichlet regression

Hi Tommy,

Did you try your models yet? I can’t really comment on Dirichlet Regression in particular, but you should be able to add the phylogenetic correlation (or covariance) matrix into basically any brms model as a group-level effect using the same notation (the data2 = list(...) part) as a gaussian or poisson regression (both are in that vignette I think). I have used the phylo matrix in ordinal models, binomial models, etc and everything works as it should! HERE is a manuscript that uses the brms phylo models in categorical regressions.

The same thing goes for the intra-specific variation part. It’s just an additional group-level effect.

One thing to consider is that if you are interested in the Pagel’s lambda value, that may not be an option in non-gaussian regression (or regressions that don’t include a sigma variable). There is some discussion of that here in this post.

Good luck!

[edit: spelling]