Quasi-Monte Carlo Methods: What, Why, and How?

Abstract

Many questions in quantitative finance, uncertainty quantification, and other disciplines are answered by computing the population mean, μ := E(Y), where instances of Y:=f(X) may be generated by numerical simulation and X has a simple probability distribution. The population mean can be approximated by the sample mean, μn := n-1 Σi=0n-1 f(xi) for a well chosen sequence of nodes, \x0, x1, …\ and a sufficiently large sample size, n. Computing μ is equivalent to computing a d-dimensional integral, ∫ f(x) (x) \, d x, where is the probability density for X. Quasi-Monte Carlo methods replace independent and identically distributed sequences of random vector nodes, \xi \i = 0∞, by low discrepancy sequences. This accelerates the convergence of μn to μ as n ∞. This tutorial describes low discrepancy sequences and their quality measures. We demonstrate the performance gains possible with quasi-Monte Carlo methods. Moreover, we describe how to formulate problems to realize the greatest performance gains using quasi-Monte Carlo. We also briefly describe the use of quasi-Monte Carlo methods for problems beyond computing the mean, μ.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…