Added single_pairs.R convenience function and updated reg_fs.R with covariate support (still not done)