
clear all
close all
%--------------------------------
  lw = 2.5; 
  set(0, 'DefaultAxesFontSize', 16);
  fs = 16;
  msize = 10;
%----------------------------------------
% Needs m-file histo.m
disp('Output From BUGS: alpha, beta, tau, and sigma')
load 'C:\Springer\BUGS\BUGSdat\alpha.txt'
load 'C:\Springer\BUGS\BUGSdat\beta.txt'
load 'C:\Springer\BUGS\BUGSdat\tau.txt'
load 'C:\Springer\BUGS\BUGSdat\sigma.txt'
alph = alpha(:,2); bet = beta(:,2); ta = tau(:,2); sig = sigma(:,2);
figure(1)
histo(alph(alph < 6 & alph > 0),100,0,1)
axis tight
%print -depsc 'C:\Brani\Courses\Bayes\Handouts\Working14\Figs\alpha_bugs.eps'

figure(2)
histo(bet(bet < 3 & bet > -1),100,0,1)
axis tight
%print -depsc 'C:\Brani\Courses\Bayes\Handouts\Working14\Figs\beta_bugs.eps'
figure(3)
histo(ta(ta < 10),100,0,1)
axis tight
%print -depsc 'C:\Brani\Courses\Bayes\Handouts\Working14\Figs\tau_bugs.eps'
figure(4)
histo(sig(sig<5),100,0,1)
axis tight
%print -depsc 'C:\Brani\Courses\Bayes\Handouts\Working14\Figs\sigma_bugs.eps'
