Quantiles Dirichlet distribution

Quantiles of multivariate distributions are not unique, so you’ll get the same error for other multivariate distributions:

julia> quantile(MvNormal(2, 1.0), 0.5)
ERROR: MethodError: no method matching iterate(::ZeroMeanIsoNormal{Tuple{Base.OneTo{Int64}}})