GRM for REML run on subset