A probability distribution on Euclidean d-space can be described by its zonoid regions. These regions form a nested family of convex sets around the expectation, each being closed and bounded. The zonoid regions of an empirical distribution introduce an ordering of the data that has many applications in multivariate statistical analysis, e.g.cluster analysis, tests for multivariate location and scale, and risk analysis. An exact algorithm is developed to constructing the zonoid regions of a d-variate empirical distribution by their facets when d>=3. The vertices of the region and their adjacency are characterized, and a procedure is suggested by which all vertices and facets can be determined. The algorithm is available as an R-package.
