% run_me



[p1,p2,y,s2]=setup('ahead');

[Mp1,Mp2]=gibbs(y,s2,10000);

plot(p2-p1,mean(Mp2-Mp1),'o');

xlabel('observed difference')

ylabel('diff. in posterior means')

print -deps 'c:\matlab\bayes\hierarch\plot.eps'

