01.08.2013 Views

Information Theory, Inference, and Learning ... - MAELabs UCSD

Information Theory, Inference, and Learning ... - MAELabs UCSD

Information Theory, Inference, and Learning ... - MAELabs UCSD

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Copyright Cambridge University Press 2003. On-screen viewing permitted. Printing not permitted. http://www.cambridge.org/0521642981<br />

You can buy this book for 30 pounds or $50. See http://www.inference.phy.cam.ac.uk/mackay/itila/ for links.<br />

368 29 — Monte Carlo Methods<br />

(a)<br />

12<br />

10<br />

8<br />

6<br />

4<br />

2<br />

0<br />

40<br />

35<br />

30<br />

25<br />

20<br />

15<br />

10<br />

5<br />

0<br />

90<br />

80<br />

70<br />

60<br />

50<br />

40<br />

30<br />

20<br />

10<br />

0<br />

(b) Metropolis<br />

100 iterations<br />

0 5 10 15 20<br />

400 iterations<br />

0 5 10 15 20<br />

1200 iterations<br />

0 5 10 15 20<br />

12<br />

10<br />

8<br />

6<br />

4<br />

2<br />

0<br />

40<br />

35<br />

30<br />

25<br />

20<br />

15<br />

10<br />

5<br />

0<br />

90<br />

80<br />

70<br />

60<br />

50<br />

40<br />

30<br />

20<br />

10<br />

0<br />

(c) Independent sampling<br />

100 iterations<br />

0 5 10 15 20<br />

400 iterations<br />

0 5 10 15 20<br />

1200 iterations<br />

0 5 10 15 20<br />

Figure 29.12. Metropolis method<br />

for a toy problem. (a) The state<br />

sequence for t = 1, . . . , 600.<br />

Horizontal direction = states from<br />

0 to 20; vertical direction = time<br />

from 1 to 600; the cross bars mark<br />

time intervals of duration 50. (b)<br />

Histogram of occupancy of the<br />

states after 100, 400, <strong>and</strong> 1200<br />

iterations. (c) For comparison,<br />

histograms resulting when<br />

successive points are drawn<br />

independently from the target<br />

distribution.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!