I am writing a post about fitting Dirichlet regression with
brms. I am wondering if it would be possible to estimate category-specific \phi. I tried this way
bf(Y ~ X, phi ~ 1, family = Dirichlet)
get_prior only proposes me one intercept and not one for each category.
Is it possible to do this in some way?
I was wondering if it would be generally a good idea to do so, but I already used this with a hand-written stan code, and it did not seem to cause an identifiability problem.
Thank you very much!
Please also provide the following information in addition to your question:
- Operating System:
- brms Version: 2.10.0