Article,

Revisiting the pseudorandom number generator ran1 from the NUMERICAL RECIPES

.
Computational Statistics & Data Analysis, 27 (4): 487--495 (June 1998)
DOI: 10.1016/S0167-9473(98)80002-9

Abstract

In this contribution we show the danger which can arise from the use of the pseudorandom number generator ran1 (and partially ran2 as well) from the first edition of the monograph Numerical Recipes in C for some Monte Carlo simulations. This paper complements the paper of Baker (1997) published in SSN in CSDA recently. We show experimentally that the main problem is not in the implementation of these generators but in “their definition”. Two examples show what can happen when we try to establish critical values through simulations basing the necessary calculations on these generators. The third example, “casting doubts on the validity of the central limit theorem”, demonstrates another dangerous situation. Finally, a fourth example shows the high sensitivity of the random number generators studied to changes of the constants defining them. These examples again clearly show the imperative need to check pseudorandom generators before their use, as we do (or at least as we should do) with the theorems we apply. Moreover, it warns against the blind belief in the not-always-valid Czech proverb that … What is written is given .

Tags

Users

  • @earthfare

Comments and Reviews