pb__modcomp function

Model comparison using parametric bootstrap methods.