pb-refdist function

Calculate reference distribution using parametric bootstrap