migpd function

Fit multiple independent generalized Pareto models