Answer by Jarle Tufto for How to fit a finite mixture of Dirichlet distributions
Reparameterizing in terms of theta1=log(p1/p3), theta2=log(p2/p3) and logs of all 9 alpha parameters, and then maximising the log likelihood using optim() with method="BFGS" seems to work if using...
View ArticleHow to fit a finite mixture of Dirichlet distributions
I have a compositional sample and I would like to fit a finite mixture of Dirichlet distributions. To be more precise, consider the following example:library(gtools)set.seed(1)PROB = c(0.25, 0.15,...
View Article
