How to fit a Dirichlet Distribution using MLE in Distributions.jl?

You could just write a call to Optim.jl to optimize the logpdf.