Uses the parametric bootstrap to calculate the p-value for a REML ratio test to compare two models.
Useful links