A Rejection Technique for Sampling from Log-Concave Multivariate Distributions

Leydold, Josef (1998) A Rejection Technique for Sampling from Log-Concave Multivariate Distributions. Preprint Series / Department of Applied Statistics and Data Processing, 21. Department of Statistics and Mathematics, Abt. f. Angewandte Statistik u. Datenverarbeitung, WU Vienna University of Economics and Business, Vienna.


Download (400kB)


Different universal methods (also called automatic or black-box methods) have been suggested to sample from univariate log-concave distributions. The description of a suitable universal generator for multivariate distributions in arbitrary dimensions has not been published up to now. The new algorithm is based on the method of transformed density rejection. To construct a hat function for the rejection algorithm the multivariate density is tranformed by a proper transformation T into a concave function (in the case of log-concave density T(x) = log(x).) Then it is possible to construct a dominating function by taking the minimum of several tangent hyperplanes which are transformed back by $T^(-1)$ into the original scale. The domains of different pieces of the hat function are polyhedra in the multivariate case. Although this method can be shown to work, it is too slow and complicated in higher dimensions. In this paper we split the $R^n$ into simple cones. The hat function is constructed piecewise on each of the cones by tangent hyperplanes. The resulting function is not continuous any more and the rejection constant is bounded from below but the setup and the generation remains quite fast in higher dimensions, e.g. n=8. The paper describes the details how this main idea can be used to construct algorithm TDRMV that generates random tuples from multivariate log-concave distribution with a computable density. Although the developed algorithm is not a real black box method it is adjustable for a large class of log-concave densities. (author's abstract)

Item Type: Paper
Additional Information: published in: ACM Transactions on Modeling and Computer Simulation 8(3), 1998, pp. 254-280
Keywords: Rejection method / multivariate log-concave distributions / universal method
Classification Codes: CCS_G.3
Divisions: Departments > Finance, Accounting and Statistics > Statistics and Mathematics
Depositing User: Repository Administrator
Date Deposited: 10 Jul 2006 10:13
Last Modified: 22 Oct 2019 00:41
Related URLs:
URI: https://epub.wu.ac.at/id/eprint/946


View Item View Item


Downloads per month over past year

View more statistics