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
MCMC (Markov Chain Monte Carlo) provides a
simple metropolis algorithm for generating random samples of points from L(a,b)
MCMC steps
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