  In this note we describe experiments on an implementation of two methods
proposed in the literature for computing regions that correspond to a notion of
order statistics for multidimensional data. Our implementation, which works for
any dimension greater than one, is the only that we know of to be publicly
available. Experiments run using the software confirm that half-space peeling
generally gives better results than directly peeling convex hulls, but at a
computational cost.
