remMap is developed for fitting multivariate response regression models under the high-dimension-low-sample-size setting