Dear Gregory R. Warnes, I'm a phd student in statistics, at the University of Milano Bicocca. I'm interested to the methods of estimate the parameters of the Dirichlet distribution. Do you have implemented an algorithm in R? If so, can you give me the script? or, in general, some helps? Best Regards Gianna --------------------------------------------------------------- Gianna Monti Dottorato di Ricerca in Statistica Università degli Studi di Milano Bicocca Via Bicocca degli Arcimboldi, 8 20126 Milano e-mail: gianna.monti@unimib.it [[alternative HTML version deleted]]
Gianna Monti <gianna.monti <at> unimib.it> writes:> > Dear Gregory R. Warnes, > I'm a phd student in statistics, at the University of Milano Bicocca. > I'm interested to the methods of estimate the parameters of the Dirichletdistribution.> Do you have implemented an algorithm in R? > If so, can you give me the script? or, in general, some helps?Scrolling through the results of RSiteSearch("dirichlet") suggests some useful tools in the VGAM package. The gtools package and MCMC packages also have ddirichlet() functions that you could use to construct a (negative log) likelihood function and optimize with optim/nlmin/etc. Ben Bolker
Reasonably Related Threads
- Rpart - best split selection for class method and Gini splitting index
- Problems with APC Smart-UPS 1500 USB and, newhidups
- Problems with APC Smart-UPS 1500 USB and newhidups driver
- User defined split function in Rpart
- ordered Twoing criterion in classification trees