Questions tagged [monte-carlo]

Using (pseudo-)random numbers and the Law of Large Numbers to simulate the random behavior of a real system.

Monte Carlo methods are used to mimic the random behavior of real random phenomena by drawing sequences of random or pseudo-random numbers, usually via a computer code. These numbers are incorporated as inputs to the model describing the phenomena of interest, as proxy to a true random sample. They require a (pseudo-)random uniform generator; the most common implementations provide a number from 0 to 1 that mimics the random draw of a $U(0,1)$ variate, let us call this (pseudo-)random variate $U$.

For instance, to simulate a fair coin toss, draw a random number $U$ and associate values in $[0,0.5)$ with heads, and values in $[0.5,1)$, with tails. More complex models require a more involved use of the random uniform generators, as well as smart ways of transforming the basic $U[0,1]$ variates into random variables with the targeted distribution, e.g., using acceptance-rejection methods or Metropolis algorithm. In general, a random generator will be a function $\Psi$ of a random number of uniform variates, $X=\Psi(U_1,U_2,\ldots)$.

In Monte Carlo studies where time and budget are an issue, an efficient use of variance reduction techniques is often warranted. These include antithetic samples, importance sampling, control variates, etc.

Besides the direct use of Monte Carlo methods to model random phenomena, they are also used in mathematical physics and statistics to evaluate multidimensional integrals. Such integrals can be represented as expectations of an integrand against a wide range of (importance) distributions. A point is drawn randomly from such a distribution (often with support the support of the integrand); the value of the integrand is added to the running sum; the process is repeated until the desired accuracy is achieved. The approximation is converging by virtue of the Law of Large Numbers. The advantage of the method is that the error decreases at the rate $O( N^{-1/2})$, independent of the dimension of the space. Low discrepancy, or quasi-Monte Carlo, sequences are arguably better suited for this purpose.

References:

Lemieux, C. (2009) Monte Carlo and Quasi-Monte Carlo Sampling. Springer. Amazon link.

Liu, J. S. (2001) Monte Carlo Strategies in Scientific Computing. Springer. Amazon link.

Robert, C.P. and Casella, G. (2004) Monte Carlo Statistical Methods. Springer. Amazon link

See also: tags sampling, simulation.

1313 questions
39
votes
7 answers

Are all simulation methods some form of Monte Carlo?

Is there a simulation method that is not Monte Carlo? All simulation methods involve substituting random numbers into the function to find a range of values for the function. So are all simulation methods in essence Monte Carlo methods?
Victor
  • 6,565
29
votes
5 answers

Why use Monte Carlo method instead of a simple grid?

when integrating a function or in complex simulations, I have seen the Monte Carlo method is widely used. I'm asking myself why one doesn't generate a grid of points to integrate a function instead of drawing random points. Wouldn't that bring more…
7
votes
3 answers

Estimating Pi using Monte Carlo simulations but without using fractional numbers

Euler's number can be estimated by using an array of size $n$, randomly permuting it and checking whether it is a derangement or not. The simulation repeats this process several times and keeps track of how many times this permutation was not a…
Vk1
  • 73
6
votes
1 answer

Monte Carlo estimation of convex hull overlap probability

This is a statistical version of my Math.SE post. Given natural numbers $b$ and $r$, uniformly randomly choose $b+r$ points within a unit square. Call the $b$ points the blue points and the $r$ points red points. Let $p(b,r)$ be the probability…
PKG
  • 461
5
votes
1 answer

Monte-Carlo estimation of the mean chord length in a polygon

This is a more "statistical" version of my Math.SE problem . Question Let $\mathbf{P}$ be a convex polygon of $m$ sides. Pick two of its edges at random, and further pick a random point on each of these edges, and compute the length of the chord…
PKG
  • 461
5
votes
1 answer

Run a Monte Carlo simulation of a startup while conceiving it? Has it been done?

As an entrepreneur you try to find a model of your business. If you simulate this via Monte Carlo you could see sensitive areas in e.g. the marketing strategy and how it affects e.g. sales. I like to think of it as a prediction device. My intuition…
Roland Kofler
  • 681
  • 1
  • 6
  • 16
4
votes
6 answers

What free tool can I use to do simple Monte Carlo simulations on OS X?

What free tool can I use to do simple Monte Carlo simulations on OS X?
sorin
  • 191
4
votes
1 answer

Confusion related to calculation of pi

I was referring to this lecture http://videolectures.net/mlss09uk_murray_mcmc/. However, I didn't get how the pi value was calculated. Here is a screenshot As far as I know the integral gives the area under the curve. So I didn't get how the double…
user34790
  • 6,757
  • 10
  • 46
  • 69
4
votes
1 answer

Does the number of Required Monte Carlo simulations increase for including more input variables

This question describes a method to calculate the number of Monte Carlo simulation runs required. Another method checks the convergence of the mean of a particular output variable. Both of these methods focus on the output variables without regard…
Patrick
  • 53
4
votes
1 answer

Is there an intuitive way to understand a sequential Monte Carlo Markov Chain?

I was wondering if there was an intuitive way to understand what a sequential Monte Carlo Markov Chain is. It also goes by the name of particle filtering. Is there an intuitive way to think about it/teach it to someone? Thanks!
user123276
  • 2,077
4
votes
1 answer

What is a good Monte Carlo algorithm for generating an archipelago in a particular map?

Sort of like the archipelago map in Age of Empires II (But this is actually a real question that I do need for my exoplanet research)
3
votes
2 answers

Proof of Marsaglia polar method

I studied Polar method and I can use it very well to simulate to Standard Normal Variable. But I can't figure it out that how it works! So is there any proof/theorem to learn reasoning behind Polar method? Thanks! Regards,
2
votes
1 answer

A Monte Carlo simulation of the "winners' curse" in a poorly powered experimental study

I am trying to carry out a Monte Carlo simulation that results in N pairings of N p-values from Student's t-test and N values of Cohen's d, such that the resulting vectors of p-values and d-values can be joined in a data.frame I can subset to…
2
votes
1 answer

Calculate certainty of Monte Carlo simulation

(Hi, sorry, this is probably a very entry level question for this site. Let me know if something is not OK.) Let's say that we use the Monte Carlo method to estimate the area of an object, in the exact same way you'd use it to estimate the value of…
2
votes
2 answers

Metropolis Algorithm for approximating integrals

I've been trying to approximate $\int_a^b f(x)\hspace{1pt} dx$ with the weight function $W(x)=A(e^{x^2}-1)$, $A$ is a constant of normalization (when $x\in[0,1]$ then $A=2.16145...$). I don't have any problems when the integration limits are $x \in…
Leinad
  • 21
1
2 3