%% compute the expected and the lower and upper bounds for parameters

post_lam1=sort(posterior_sim(:,1));
post_lam2=sort(posterior_sim(:,2));
post_al=sort(posterior_sim(:,3));

display('true -- lower - expected - upper for lambda1')
disp([p(1) post_lam1(int32((MM*.05))) mean(post_lam1) post_lam1(int32((MM*.95)))])
display('true -- lower - expected - upper for lambda2')
disp([p(2) post_lam2(int32((MM*.05))) mean(post_lam2) post_lam2(int32((MM*.95)))])
display('true -- lower - expected - upper for alpha')
disp([p(3) post_al(int32((MM*.05))) mean(post_al) post_al(int32((MM*.95)))])