论文信息 - Algorithm 802: an automatic generator for bivariate log-concave distributions

Algorithm 802: an automatic generator for bivariate log-concave distributions

Different automatic (also called universal or black-box) methods have been suggested to sample from univariate log-concave distributions. Our new automatic algorithm for bivariate log-concave distributions is based on the method of transformed density rejection. In order to construct a hat function for a rejection algorithm the bivariate density is transformed by the logarithm into a concave function. Then it is possible to construct a dominating function by taking the minimum of several tangent planes, which are by exponentiation transformed back into the original scale. The choice of the points of contact is automated using adaptive rejection sampling. This means that points that are rejected by the rejection algorithm can be used as additional points of contact. The article describes the details how this main idea can be used to construct Algorithm ALC2D that can generate random pairs from all bivariate log-concave distributions with known domain, computable density, and computable partial derivatives.

Wolfgang Hörmann | W. Hörmann

[1] Wolfgang Hörmann,et al. A rejection technique for sampling from T-concave distributions , 1995, TOMS.

[2] Wolfgang Hörmann,et al. A sweep-plane algorithm for generating random tuples in simple polytopes , 1998, Math. Comput..

[3] Wolfgang Hörmann,et al. Universal Generators for Correlation Induction , 1994 .

[4] J. Dagpunar. Principles of Random Variate Generation , 1988 .

[5] Josef Leydold,et al. A rejection technique for sampling from log-concave multivariate distributions , 1998, TOMC.

[6] Mark E. Johnson. Multivariate Statistical Simulation: Johnson/Multivariate , 1987 .

[7] L. Devroye. Non-Uniform Random Variate Generation , 1986 .

[8] Luc Devroye,et al. Random variate generation for multivariate unimodal densities , 1997, TOMC.

[9] Jennie Malboeuf. Algorithm , 1994, Neurology.

[10] M. E. Johnson,et al. Multivariate Statistical Simulation , 1988 .

[11] Ronald W. Golland,et al. Principles of Random Variate Generation , 1990 .

[12] J. G. Mauldon,et al. A GENERALIZATION OF THE BETA-DISTRIBUTION' , 1959 .

[13] W. Gilks,et al. Adaptive Rejection Sampling for Gibbs Sampling , 1992 .

[14] Adrian F. M. Smith,et al. Efficient generation of random variates via the ratio-of-uniforms method , 1991 .

[15] M. Evans,et al. Random Variable Generation Using Concavity Properties of Transformed Densities , 1998 .