Add function to perform monte carlo bootstrapping

Based on excellent idea by @beckyfisher to deal with uncertainty due to censoring and/or uncertainty in ssd data points by using monte carlo sampling to draw from the data distributions and then perform parametric bootstrapping on the sampled data.