function plotinterval(int,mu)
% plotinterval(int,mu) - plot confidence intervals
% plotinterval(INT,MU) where INT is an N x 2 matrix of confidence intervals
% and MU is a scalar, plots the N confidence intervals and the mean MU
[N,test]=size(int);
if test~=2, error('plotinterval: interval must be an N x 2 matrix of intervals'); end

for i=1:N
    plot([i i],[int(i,1) int(i,2)],'k');
    hold on
end;
plot([0 N+1], [mu mu],'r');
hold off;
set(gca,'Xlim',[0 N+1])