powered by
Simple Metropolis algorithm demonstration using a bivariate exponential target from Chapter 1 (pages 27-30).
expo.metrop(m,x,y,L1,L2,L,B)
number of iterations
starting point for the x vector
starting point for the y vector
event intensity for the x dimension
event intensity for the y dimension
shared event intensity
upper bound
Jeff Gill
expo.metrop(m=5000, x=0.5, y=0.5, L1=0.5, L2=0.1, L=0.01, B=8)
Run the code above in your browser using DataLab