Monte Carlo Methods Flashcards
what are monte carlo methods
methods for generating random variables –
i.e. samples of numbers which
behave as if they are drawn
from some particular pdf (e.g.
uniform, Gaussian, Poisson etc).
algorithms only generate pseudo-random numbers which are
very long (deterministic) sequences of numbers which are
approximately random (i.e. no discernible pattern).
the better the random number generator, the
better it approximates U[0,1]
phase portraits
scatterplots of the ith value against against the (i+1)th value
we can compute the auto-correlation function where j is known as
the lag
if the sequence is uniformly random, we expect
p(j)=1 for j=0
p(j)=0 otherwise
generating random numbers from other pdfs can be done by
transforming random numbers drawn from simpler pdfs
probability integral transform steps
- sample a random variable y~U[0,1]
- compute x such that y=p(x)
- then x~p(x) ie x is drawn from the pdf corresponding to the cdf
MCMC provides a
simple metropolis algorithm for generating random samples of points from L(a,b)
MCMC steps
- sample random initial point
- create a new pdf Q called the proposal density on p1
- sample tentative new point P’ from Q
- compute R
- if R>1 we accept next point
- If R<1 we reject
markov chain
acceptance probability depends only on the previous point
If the variance of the proposal density is too small (σ=0.01) then
not all of the distribution is sampled
If the variance of the proposal density is too large (σ=50) then
values jump around significantly
if the variance of the proposal density is correct (σ=1) then
the entire distribution is sampled