Keeping the noise down: common random numbers for disease simulation modeling

Disease simulation models are used to conduct decision analyses of the comparative benefits and risks associated with preventive and treatment strategies. To address increasing model complexity and computational intensity, modelers use variance reduction techniques to reduce stochastic noise and improve computational efficiency. One technique, common random numbers, further allows modelers to conduct counterfactual-like analyses with direct computation of statistics at the individual level. This technique uses synchronized random numbers across model runs to induce correlation in model output thereby making differences easier to distinguish as well as simulating identical individuals across model runs. We provide a tutorial introduction and demonstrate the application of common random numbers in an individual-level simulation model of the epidemiology of breast cancer.

[1]  Averill M. Law,et al.  Simulation Modeling and Analysis , 1982 .

[2]  Pierre L'Ecuyer,et al.  An Object-Oriented Random-Number Package with Many Long Streams and Substreams , 2002, Oper. Res..

[3]  M C Weinstein,et al.  The cost effectiveness of combination antiretroviral therapy for HIV disease. , 2001, The New England journal of medicine.

[4]  Gordon M. Clark,et al.  Use of common random numbers in comparing alternatives , 1990, 1990 Winter Simulation Conference Proceedings.

[5]  Pierre L'Ecuyer,et al.  TestU01: A C library for empirical testing of random number generators , 2006, TOMS.

[6]  J. B. Peacock,et al.  Statistical Distributions: A Handbook for Students and Practitioners , 1975 .

[7]  J D Habbema,et al.  The MISCAN simulation program for the evaluation of screening for disease. , 1985, Computer methods and programs in biomedicine.

[8]  Bernard R. Rosner,et al.  Fundamentals of Biostatistics. , 1992 .

[9]  Natasha K. Stout,et al.  Chapter 7: The Wisconsin Breast Cancer Epidemiology Simulation Model , 2006 .

[10]  Bruce W. Schmeiser.  Some myths and common errors in simulation experiments , 2001, Proceeding of the 2001 Winter Simulation Conference (Cat. No.01CH37304).

[11]  Ruth Davies,et al.  Stream correlations in multiple recursive and congruential generators , 2007 .

[12]  Barry L. Nelson.  Robust multiple comparisons under common random numbers , 1993, TOMC.

[13]  Makoto Matsumoto,et al.  Common defects in initialization of pseudorandom number generators , 2007, TOMC.

[14]  M. Weinstein,et al.  A comprehensive natural history model of HPV infection and cervical cancer to estimate the clinical impact of a prophylactic HPV‐16/18 vaccine , 2003, International journal of cancer.

[15]  Takuji Nishimura,et al.  Mersenne twister: a 623-dimensionally equidistributed uniform pseudo-random number generator , 1998, TOMC.

[16]  Oguzhan Alagoz,et al.  Incorporating Biological Natural History in Simulation Models: Empirical Estimates of the Progression of End-Stage Liver Disease , 2005, Medical decision making : an international journal of the Society for Medical Decision Making.

[17]  Stephen E. Chick,et al.  A decision-theoretic approach to screening and selection with common random numbers , 1999, WSC '99.

[18]  A. Huitson,et al.  Statistical distributions: A handbook for students and practitioners , 1975 .

[19]  P. Hellekalek.  Good random number generators are (not so) easy to find , 1998 .

[20]  Ruth Davies,et al.  A Simulation to Evaluate Screening for Helicobacter Pylori Infection in the Prevention of Peptic Ulcers and Gastric Cancers , 2002, Health care management science.

[21]  Mark S Roberts,et al.  Increasing the Efficiency of Monte Carlo Cohort Simulations with Variance Reduction Techniques , 2006, Medical decision making : an international journal of the Society for Medical Decision Making.