The acceptance-rejection method for low-discrepancy sequences

Nguyet Nguyen and Giray \"Okten

Papers from arXiv.org

Abstract: Generation of pseudorandom numbers from different probability distributions has been studied extensively in the Monte Carlo simulation literature. Two standard generation techniques are the acceptance-rejection and inverse transformation methods. An alternative approach to Monte Carlo simulation is the quasi-Monte Carlo method, which uses low-discrepancy sequences, instead of pseudorandom numbers, in simulation. Low-discrepancy sequences from different distributions can be obtained by the inverse transformation method, just like for pseudorandom numbers. In this paper, we will present an acceptance-rejection algorithm for low-discrepancy sequences. We will prove a convergence result, and present error bounds. We will then use this acceptance-rejection algorithm to develop quasi-Monte Carlo versions of some well known algorithms to generate beta and gamma distributions, and investigate the efficiency of these algorithms numerically. We will also consider the simulation of the variance gamma model, a model used in computational finance, where the generation of these probability distributions are needed. Our results show that the acceptance-rejection technique can result in significant improvements in computing time over the inverse transformation method in the context of low-discrepancy sequences.

Date: 2014-03
New Economics Papers: this item is included in nep-cmp
References: View references in EconPapers View complete reference list from CitEc
Citations:

Downloads: (external link)
http://arxiv.org/pdf/1403.5599 Latest version (application/pdf)

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:arx:papers:1403.5599

Access Statistics for this paper

More papers in Papers from arXiv.org
Bibliographic data for series maintained by arXiv administrators ().