Scalable random number generation for truncated log-concave distributions

04/04/2022
by   Lambardi di San Miniato, et al.
0

Inverse transform sampling is an exceptionally general method to generate non-uniform-distributed random numbers, but can be rather unstable when simulating extremely truncated distributions. Many famous probability models share a property called log-concavity, which is not affected by truncation, so they can all be simulated via rejection sampling using Devroye's approach. This sampler is based on rejection and thus more stable than inverse transform, and uses a very simple envelope whose acceptance rate is guaranteed to be at least 20%. The aim of this paper is threefold: firstly, to warn against the risk of wrongly simulating from truncated distributions; secondly, to motivate a more extensive use of rejection sampling to mitigate the issues; lastly, to motivate Devroye's automatic method as a practical standard in the case of log-concave distributions. We illustrate the proposal by means of simulations based on some Tweedie distributions, for their relevance in regression analysis.

READ FULL TEXT

page 9

page 15

page 16

page 17

page 21

research
12/17/2021

Moments and random number generation for the truncated elliptical family of distributions

This paper proposes an algorithm to generate random numbers from any mem...
research
08/19/2019

Probability Estimation with Truncated Inverse Binomial Sampling

In this paper, we develop a general theory of truncated inverse binomial...
research
11/01/2018

Entropy versus variance for symmetric log-concave random variables and related problems

We show that the uniform distribution minimises entropy among all symmet...
research
06/12/2019

Flexible Modeling of Diversity with Strongly Log-Concave Distributions

Strongly log-concave (SLC) distributions are a rich class of discrete pr...
research
05/17/2021

On log-concave approximations of high-dimensional posterior measures and stability properties in non-linear inverse problems

The problem of efficiently generating random samples from high-dimension...
research
10/21/2020

Optimal dual quantizers of 1D log-concave distributions: uniqueness and Lloyd like algorithm

We establish for dual quantization the counterpart of Kieffer's uniquene...
research
12/24/2021

Concave-Convex PDMP-based sampling

Recently non-reversible samplers based on simulating piecewise determini...

Please sign up or login with your details

Forgot password? Click here to reset