--reml to loop over multiple phenotypes?