Make ML-refits optional

Minor change to enable RLR-based tests instead of LR-tests in PBmodcomp for testing variance components. Adds a `REML` argument (default: `FALSE`, so the standard behavior does not change) to PBmodcomp and the functions it calls.
Sorry about the huge diff, seems like there were encoding issues.