A Concave Pairwise Fusion Approach to Clustering of Multi-Response Regression and Its Robust Extensions