Given a set of predictive quantiles from a distribution, estimate the
distribution and create d, p, q,
and r functions to evaluate its density function,
distribution function, and quantile function, and generate random
samples. On the interior of the provided quantiles, an interpolation
method such as a monotonic cubic spline is used; the tails are
approximated by a location-scale family.