fastrml function

Rapidly converging algorithm for restricted maximum-likelihood (RML) estimation in linear mixed models