Returns a hilo central probability interval with probability coverage of
size. By default, the distribution's quantile() will be used to compute
the lower and upper bound for a centered interval
# S3 method for class 'distribution'
hilo(x, size = 95, ...)