bootstrap_rho function

Perform a parametric bootstrapping correction on an estimated rho vector