MathGroup Archive 2003

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Warning Long Post: Help Speeding up algorithm


Mark Coleman <mark at markscoleman.com> wrote...

> drawing random numbers and random vectors from nonstandard distributions. 
...
> Acceptance-rejection algorithms repeatedly draw a random number from a 
> standard distribution, ...
> Given that one  (obviously) cannot tell in advance how many passes are 
> required to generate a satisfactory random draw, it is not clear to me
> how to code  this more efficiently in Mathematica, ideally without
> the explicit loop. 

See section 2.6 of:
  Rose/Smith, Mathematical Statistics with Mathematica -> pages 72-79. 

At pp.77-79, they specifically discuss the rejection method in Mathematica, 
and how to resolve the 'loop' problem you refer to (in 1 line). Recommended.

Have a good one

Simon


  • Prev by Date: Re: Integration
  • Next by Date: NDSolve question
  • Previous by thread: Warning Long Post: Help Speeding up algorithm
  • Next by thread: Trigonometric math functions