R-optimal designs for multi-response regression models with multi-factors