Computing multiple-output regression quantile regions

A procedure relying on linear programming techniques is developed to compute (regression) quantile regions that have been defined recently. In the location case, this procedure allows for computing halfspace depth regions even beyond dimension two. The corresponding algorithm is described in detail, and illustrations are provided both for simulated and real data. The efficiency of a Matlab implementation of the algorithm1 is also investigated through extensive simulations.