Optimizing Code to Fit a Generalized Pareto Distribution to Data

Btw, Distributions.jl currently doesn’t fit a GPD.
It’d be great to add this: Distributions.jl/generalizedpareto.jl at master · JuliaStats/Distributions.jl · GitHub