Dirichlet regression in R/c for translation to Python task