Generating gamma variates by a modified rejection technique

  • Authors:
  • J. H. Ahrens;U. Dieter

  • Affiliations:
  • Univ. of Kiel, West Germany;Technical Univ., Graz, Austria

  • Venue:
  • Communications of the ACM
  • Year:
  • 1982

Quantified Score

Hi-index 48.23

Visualization

Abstract

A suitable square root transformation of a gamma random variable with mean a ≥ 1 yields a probability density close to the standard normal density. A modification of the rejection technique then begins by sampling from the normal distribution, being able to accept and transform the initial normal observation quickly at least 85 percent of the time (95 percent if a ≥ 4). When used with efficient subroutines for sampling from the normal and exponential distributions, the resulting accurate method is significantly faster than competing algorithms.