numb
by numb
grid
from a Potts distribution with four colours and a four neighbour
structure, relying on niter
iterations of a standard Gibbs sampler.pottsgibbs(niter, numb, beta)
pottshm
ex=pottsgibbs(100,15,.4)
image(ex)
Run the code above in your browser using DataLab