The patchwork rejection technique for sampling from unimodal distributions

  • Authors:
  • Ernst Stadlober;Heinz Zechner

  • Affiliations:
  • Technical Univ. Graz, Graz, Austria;Technical Univ. Graz, Graz, Austria

  • Venue:
  • ACM Transactions on Modeling and Computer Simulation (TOMACS)
  • Year:
  • 1999

Quantified Score

Hi-index 0.00

Visualization

Abstract

We report on both theoretical developments and comutational experience with the patchwork rejection technique in Zechner and Stadlober [1993] and Zechner [1997]. The basic approach is due to Minh [1988], who suggested a special sampling method for the gamma distribution. The method's general objective is to rearrange the area below the density of histogram f (x) in the body of the distribution by certain point reflections such that variates may be generated efficiently within a large center interval. This is carried out via uniform hat functions, combined with minorizing rectangles for immediate acceptance of one transformed uniform deviate. The remaining tails of f(x) are covered by exponential functions. Experiments show that patchwork rejection algorithms are in general faster than their competitors at the cost of higher set-up times.