%
%
% (C) Michael Mann, 2016.
%
%
% FIGURE 1
%
figure(1)
subplot(2,3,1)
load Runs/Unbiased100Real100Site
hold
title('Unbiased')
set(gca,'FontSize',14);
hline1=plot(nflips,pvalmedianm,'b','Linewidth',2);
hline2=plot(nflips,pvallowerm,'b--','Linewidth',1.5);
hline3=plot(nflips,pvalupperm,'b--','Linewidth',1.5);
hline4=plot(nflips,pvalsig,'k--','Linewidth',1);
ylabel('{\itp} value');
legend([hline1,hline2,hline4],'median','upper/lower','p=0.05');
set(gca,'YLim',[0 0.8],'FontName','Courier'); grid on;
set(gca,'XLim',[0 64],'FontName','Courier'); grid on;
%
subplot(2,3,2)
load Runs/Biased100Real100Site0pt6
hold
title('Biased, \theta_0=0.6')
set(gca,'FontSize',14);
plot(nflips,pvalmedianm,'b','Linewidth',2);
plot(nflips,pvallowerm,'b--','Linewidth',1.5);
plot(nflips,pvalupperm,'b--','Linewidth',1.5);
plot(nflips,pvalsig,'k--','Linewidth',1);
set(gca,'YLim',[0 0.8],'FontName','Courier'); grid on;
set(gca,'XLim',[0 64],'FontName','Courier'); grid on;
%
subplot(2,3,3)
load Runs/Biased100Real100Site0pt75
hold
title('Biased, \theta_0=0.75')
set(gca,'FontSize',14);
plot(nflips,pvalmedianm,'b','Linewidth',2);
plot(nflips,pvallowerm,'b--','Linewidth',1.5);
plot(nflips,pvalupperm,'b--','Linewidth',1.5);
plot(nflips,pvalsig,'k--','Linewidth',1);
set(gca,'YLim',[0 0.8],'FontName','Courier'); grid on;
set(gca,'XLim',[0 64],'FontName','Courier'); grid on;
%
subplot(2,3,4)
load Runs/Unbiased100Real5Site
hold
title('Unbiased')
set(gca,'FontSize',14);
plot(nflips,pvalmedianm,'b','Linewidth',2);
plot(nflips,pvallowerm,'b--','Linewidth',1.5);
plot(nflips,pvalupperm,'b--','Linewidth',1.5);
plot(nflips,pvalsig,'k--','Linewidth',1);
xlabel('{\itN_{tot}}'); ylabel('{\itp} value');
set(gca,'YLim',[0 0.8],'FontName','Courier'); grid on;
set(gca,'XLim',[0 64],'FontName','Courier'); grid on;
%
subplot(2,3,5)
load Runs/Biased100Real5Site0pt6
hold
title('Biased, \theta_0=0.6')
set(gca,'FontSize',14);
plot(nflips,pvalmedianm,'b','Linewidth',2);
plot(nflips,pvallowerm,'b--','Linewidth',1.5);
plot(nflips,pvalupperm,'b--','Linewidth',1.5);
plot(nflips,pvalsig,'k--','Linewidth',1);
xlabel('{\itN_{tot}}');
set(gca,'YLim',[0 0.8],'FontName','Courier'); grid on;
set(gca,'XLim',[0 64],'FontName','Courier'); grid on;
%
subplot(2,3,6)
load Runs/Biased100Real5Site0pt75
hold
title('Biased, \theta_0=0.75')
set(gca,'FontSize',14);
plot(nflips,pvalmedianm,'b','Linewidth',2);
plot(nflips,pvallowerm,'b--','Linewidth',1.5);
plot(nflips,pvalupperm,'b--','Linewidth',1.5);
plot(nflips,pvalsig,'k--','Linewidth',1);
xlabel('{\itN_{tot}}');
set(gca,'YLim',[0 0.8],'FontName','Courier'); grid on;
set(gca,'XLim',[0 64],'FontName','Courier'); grid on;
%
% FIGURE 2 of article
%
figure(2)
subplot(2,3,1)
load Runs/Unbiased100Real100Site
hold
title('Unbiased')
set(gca,'FontSize',14);
hline1=plot(b,pfullmedm,'k','Linewidth',2); 
hline2=plot(b,pfulllowm,'k--','Linewidth',1); 
plot(b,pfullhighm,'k--','Linewidth',1); 
hline3=plot(b,pmidmedm,'g','Linewidth',2); 
hline4=plot(b,pmidlowm,'g--','Linewidth',1); 
plot(b,pmidhighm,'g--','Linewidth',1); 
hline5=plot(b,prior,'r','Linewidth',2);
ylabel('Probability \it{P}(\theta)');
set(gca,'YLim',[0 1.1],'FontName','Courier'); grid on;
legend([hline5,hline1,hline3],'prior','final','intermediate');
subplot(2,3,2)
load Runs/Biased100Real100Site0pt6
hold
title('Biased, \theta_0=0.6 ')
set(gca,'FontSize',14);
plot(b,pfullmedm,'k','Linewidth',2); 
plot(b,pfulllowm,'k--','Linewidth',1); 
plot(b,pfullhighm,'k--','Linewidth',1); 
plot(b,pmidmedm,'g','Linewidth',2); 
plot(b,pmidlowm,'g--','Linewidth',1); 
plot(b,pmidhighm,'g--','Linewidth',1); 
plot(b,prior,'r','Linewidth',2);
set(gca,'YLim',[0 1.1],'FontName','Courier'); grid on;
%
subplot(2,3,3)
load Runs/Biased100Real100Site0pt75
hold
title('Biased, \theta_0=0.75 ')
set(gca,'FontSize',14);
plot(b,pfullmedm,'k','Linewidth',2); 
plot(b,pfulllowm,'k--','Linewidth',1); 
plot(b,pfullhighm,'k--','Linewidth',1); 
plot(b,pmidmedm,'g','Linewidth',2); 
plot(b,pmidlowm,'g--','Linewidth',1); 
plot(b,pmidhighm,'g--','Linewidth',1); 
plot(b,prior,'r','Linewidth',2);
set(gca,'YLim',[0 1.1],'FontName','Courier'); grid on;
%
subplot(2,3,4)
load Runs/Unbiased100Real5Site
hold
set(gca,'FontSize',14);
hline1=plot(b,pfullmedm,'k','Linewidth',2); 
hline2=plot(b,pfulllowm,'k--','Linewidth',1); 
plot(b,pfullhighm,'k--','Linewidth',1); 
hline3=plot(b,pmidmedm,'g','Linewidth',2); 
hline4=plot(b,pmidlowm,'g--','Linewidth',1); 
plot(b,pmidhighm,'g--','Linewidth',1); 
hline5=plot(b,prior,'r','Linewidth',2);
xlabel('bias \theta'); ylabel('Probability \it{P}(\theta)');
set(gca,'YLim',[0 1.1],'FontName','Courier'); grid on;
%
subplot(2,3,5)
load Runs/Biased100Real5Site0pt6
hold
set(gca,'FontSize',14);
plot(b,pfullmedm,'k','Linewidth',2); 
plot(b,pfulllowm,'k--','Linewidth',1); 
plot(b,pfullhighm,'k--','Linewidth',1); 
plot(b,pmidmedm,'g','Linewidth',2); 
plot(b,pmidlowm,'g--','Linewidth',1); 
plot(b,pmidhighm,'g--','Linewidth',1); 
plot(b,prior,'r','Linewidth',2);
xlabel('bias \theta'); 
set(gca,'YLim',[0 1.1],'FontName','Courier'); grid on;
%
subplot(2,3,6)
load Runs/Biased100Real5Site0pt75
hold
set(gca,'FontSize',14);
plot(b,pfullmedm,'k','Linewidth',2); 
plot(b,pfulllowm,'k--','Linewidth',1); 
plot(b,pfullhighm,'k--','Linewidth',1); 
plot(b,pmidmedm,'g','Linewidth',2); 
plot(b,pmidlowm,'g--','Linewidth',1); 
plot(b,pmidhighm,'g--','Linewidth',1); 
plot(b,prior,'r','Linewidth',2);
xlabel('bias \theta'); 
set(gca,'YLim',[0 1.1],'FontName','Courier'); grid on;
%
% FIGURE 3 of article
%
figure(3)
subplot(2,3,1)
load Runs/Unbiased100Real100Site
hold
title('Unbiased')
set(gca,'FontSize',14);
hline1=plot(nflips,medianbiasnullm,'r','Linewidth',2);
hline2=plot(nflips,mean(medianbiasnullrej),'m','Linewidth',2);
hline3=plot(nflips,mean(medianbias),'b','Linewidth',2);
hline4=plot(nflips,mean(biaslower'),'b--','Linewidth',1.5);
plot(nflips,mean(biasupper'),'b--','Linewidth',1.5);
plot(nflips,mean(biasnullrejlower'),'m--','Linewidth',1.5);
hline5=plot(nflips,mean(biasnullrejupper'),'m--','Linewidth',1.5);
ylabel('bias, \theta_0');
set(gca,'YLim',[0.4 0.78],'FontName','Courier'); grid on;
set(gca,'XLim',[0 64],'FontName','Courier'); grid on;
legend([hline3,hline1,hline2],'median','null','adjusted null');
%
subplot(2,3,2)
load Runs/Biased100Real100Site0pt6
hold
title('Biased, \theta_0=0.6')
set(gca,'FontSize',14);
plot(nflips,medianbiasnullm,'r','Linewidth',2);
plot(nflips,mean(medianbiasnullrej),'m','Linewidth',2);
plot(nflips,mean(medianbias),'b','Linewidth',2);
plot(nflips,mean(biaslower'),'b--','Linewidth',1.5);
plot(nflips,mean(biasupper'),'b--','Linewidth',1.5);
plot(nflips,mean(biasnullrejlower'),'m--','Linewidth',1.5);
plot(nflips,mean(biasnullrejupper'),'m--','Linewidth',1.5);
set(gca,'YLim',[0.4 0.78],'FontName','Courier'); grid on;
set(gca,'XLim',[0 64],'FontName','Courier'); grid on;
%
subplot(2,3,3)
load Runs/Biased100Real100Site0pt75
hold
title('Biased, \theta_0=0.75')
set(gca,'FontSize',14);
plot(nflips,medianbiasnullm,'r','Linewidth',2);
plot(nflips,mean(medianbiasnullrej),'m','Linewidth',2);
plot(nflips,mean(medianbias),'b','Linewidth',2);
plot(nflips,mean(biaslower'),'b--','Linewidth',1.5);
plot(nflips,mean(biasupper'),'b--','Linewidth',1.5);
plot(nflips,mean(biasnullrejlower'),'m--','Linewidth',1.5);
plot(nflips,mean(biasnullrejupper'),'m--','Linewidth',1.5);
set(gca,'YLim',[0.4 0.78],'FontName','Courier'); grid on;
set(gca,'XLim',[0 64],'FontName','Courier'); grid on;
%
subplot(2,3,4)
load Runs/Unbiased100Real5Site
hold
set(gca,'FontSize',14);
plot(nflips,medianbiasnullm,'r','Linewidth',2);
plot(nflips,mean(medianbiasnullrej),'m','Linewidth',2);
plot(nflips,mean(medianbias),'b','Linewidth',2);
biasupperm=biasupper(:,1);
biaslowerm=biaslower(:,1);
biasnullrejupperm=biasnullrejupper(:,1);
biasnullrejlowerm=biasnullrejlower(:,1);
plot(nflips,mean(biaslower'),'b--','Linewidth',1.5);
plot(nflips,mean(biasupper'),'b--','Linewidth',1.5);;
plot(nflips,mean(biasnullrejlower'),'m--','Linewidth',1.5);
plot(nflips,mean(biasnullrejupper'),'m--','Linewidth',1.5);
xlabel('{\itN_{tot}}'); ylabel('bias, \theta_0');
set(gca,'YLim',[0.4 0.78],'FontName','Courier'); grid on;
set(gca,'XLim',[0 64],'FontName','Courier'); grid on;
%
subplot(2,3,5)
load Runs/Biased100Real5Site0pt6
hold
set(gca,'FontSize',14);
plot(nflips,medianbiasnullm,'r','Linewidth',2);
plot(nflips,mean(medianbiasnullrej),'m','Linewidth',2);
plot(nflips,mean(medianbias),'b','Linewidth',2);
biasupperm=biasupper(:,1);
biaslowerm=biaslower(:,1);
biasnullrejupperm=biasnullrejupper(:,1);
biasnullrejlowerm=biasnullrejlower(:,1);
plot(nflips,mean(biaslower'),'b--','Linewidth',1.5);
plot(nflips,mean(biasupper'),'b--','Linewidth',1.5);
plot(nflips,mean(biasnullrejlower'),'m--','Linewidth',1.5);
plot(nflips,mean(biasnullrejupper'),'m--','Linewidth',1.5);
xlabel('{\itN_{tot}}'); 
set(gca,'YLim',[0.4 0.78],'FontName','Courier'); grid on;
set(gca,'XLim',[0 64],'FontName','Courier'); grid on;
%
subplot(2,3,6)
load Runs/Biased100Real5Site0pt75
hold
set(gca,'FontSize',14);
plot(nflips,medianbiasnullm,'r','Linewidth',2);
plot(nflips,mean(medianbiasnullrej),'m','Linewidth',2);
plot(nflips,mean(medianbias),'b','Linewidth',2);
biasupperm=biasupper(:,1);
biaslowerm=biaslower(:,1);
biasnullrejupperm=biasnullrejupper(:,1);
biasnullrejlowerm=biasnullrejlower(:,1);
plot(nflips,mean(biaslower'),'b--','Linewidth',1.5);
plot(nflips,mean(biasupper'),'b--','Linewidth',1.5);
plot(nflips,mean(biasnullrejlower'),'m--','Linewidth',1.5);
plot(nflips,mean(biasnullrejupper'),'m--','Linewidth',1.5);
xlabel('{\itN_{tot}}'); 
set(gca,'YLim',[0.4 0.78],'FontName','Courier'); grid on;
set(gca,'XLim',[0 64],'FontName','Courier'); grid on;
%
% FIGURE 4 of article
%
figure(4)
subplot(2,3,1)
load Runs/Unbiased100Real100Site
hold
title('Unbiased')
set(gca,'FontSize',14);
hline1=plot(nflipsnew,rmsenullm,'r','Linewidth',2);
hline2=plot(nflipsnew,rmsenullrejm,'m','Linewidth',2);
hline3=plot(nflipsnew,rmsebayesm,'b','Linewidth',2);
hline4=plot(nflipsnew,rmsebayeslower,'b--','Linewidth',1.5);
plot(nflipsnew,rmsebayesupper,'b--','Linewidth',1.5);
plot(nflipsnew,rmsenullrejlower,'m--','Linewidth',1.5);
hline5=plot(nflipsnew,rmsenullrejupper,'m--','Linewidth',1.5);
ylabel('\epsilon');
set(gca,'YLim',[0.0 8.0],'FontName','Courier'); grid on;
set(gca,'XLim',[0 32],'FontName','Courier'); grid on;
legend([hline3,hline1,hline2],'Bayesian','null','adjusted null');
%
subplot(2,3,2)
load Runs/Biased100Real100Site0pt6
hold
title('Biased, \theta_0=0.6')
set(gca,'FontSize',14);
plot(nflipsnew,rmsenullm,'r','Linewidth',2);
plot(nflipsnew,rmsenullrejm,'m','Linewidth',2);
plot(nflipsnew,rmsebayesm,'b','Linewidth',2);
plot(nflipsnew,rmsebayeslower,'b--','Linewidth',1.5);
plot(nflipsnew,rmsebayesupper,'b--','Linewidth',1.5);
plot(nflipsnew,rmsenullrejlower,'m--','Linewidth',1.5);
plot(nflipsnew,rmsenullrejupper,'m--','Linewidth',1.5);
set(gca,'YLim',[0.0 8.0],'FontName','Courier'); grid on;
set(gca,'XLim',[0 32],'FontName','Courier'); grid on;
%
subplot(2,3,3)
load Runs/Biased100Real100Site0pt75
hold
title('Biased, \theta_0=0.75')
set(gca,'FontSize',14);
plot(nflipsnew,rmsenullm,'r','Linewidth',2);
plot(nflipsnew,rmsenullrejm,'m','Linewidth',2);
plot(nflipsnew,rmsebayesm,'b','Linewidth',2);
plot(nflipsnew,rmsebayeslower,'b--','Linewidth',1.5);
plot(nflipsnew,rmsebayesupper,'b--','Linewidth',1.5);
plot(nflipsnew,rmsenullrejlower,'m--','Linewidth',1.5);
plot(nflipsnew,rmsenullrejupper,'m--','Linewidth',1.5);
set(gca,'YLim',[0.0 8.0],'FontName','Courier'); grid on;
set(gca,'XLim',[0 32],'FontName','Courier'); grid on;
%
subplot(2,3,4)
load Runs/Unbiased100Real5Site
hold
set(gca,'FontSize',14);
plot(nflipsnew,rmsenullm,'r','Linewidth',2);
plot(nflipsnew,rmsenullrejm,'m','Linewidth',2);
plot(nflipsnew,rmsebayesm,'b','Linewidth',2);
plot(nflipsnew,rmsebayeslower,'b--','Linewidth',1.5);
plot(nflipsnew,rmsebayesupper,'b--','Linewidth',1.5);
plot(nflipsnew,rmsenullrejlower,'m--','Linewidth',1.5);
plot(nflipsnew,rmsenullrejupper,'m--','Linewidth',1.5);
xlabel('{\itN} '); ylabel('\epsilon');
set(gca,'YLim',[0.0 8.0],'FontName','Courier'); grid on;
set(gca,'XLim',[0 32],'FontName','Courier'); grid on;
%
subplot(2,3,5)
load Runs/Biased100Real5Site0pt6
hold
set(gca,'FontSize',14);
plot(nflipsnew,rmsenullm,'r','Linewidth',2);
plot(nflipsnew,rmsenullrejm,'m','Linewidth',2);
plot(nflipsnew,rmsebayesm,'b','Linewidth',2);
plot(nflipsnew,rmsebayeslower,'b--','Linewidth',1.5);
plot(nflipsnew,rmsebayesupper,'b--','Linewidth',1.5);
plot(nflipsnew,rmsenullrejlower,'m--','Linewidth',1.5);
plot(nflipsnew,rmsenullrejupper,'m--','Linewidth',1.5);
xlabel('{\itN} '); 
set(gca,'YLim',[0.0 8.0],'FontName','Courier'); grid on;
set(gca,'XLim',[0 32],'FontName','Courier'); grid on;
%
subplot(2,3,6)
load Runs/Biased100Real5Site0pt75
hold
set(gca,'FontSize',14);
plot(nflipsnew,rmsenullm,'r','Linewidth',2);
plot(nflipsnew,rmsenullrejm,'m','Linewidth',2);
plot(nflipsnew,rmsebayesm,'b','Linewidth',2);
plot(nflipsnew,rmsebayeslower,'b--','Linewidth',1.5);
plot(nflipsnew,rmsebayesupper,'b--','Linewidth',1.5);
plot(nflipsnew,rmsenullrejlower,'m--','Linewidth',1.5);
plot(nflipsnew,rmsenullrejupper,'m--','Linewidth',1.5);
xlabel('{\itN} '); 
set(gca,'YLim',[0.0 8.0],'FontName','Courier'); grid on;
set(gca,'XLim',[0 32],'FontName','Courier'); grid on;
%
%
% Supp Figure 1
%
%
figure(5)
subplot(2,3,1)
load Runs/Unbiased100Real100Site
hold
title('Unbiased')
set(gca,'FontSize',14);
hline1=plot(nflipsnew,newheadsmeanm,'k','Linewidth',2);
hline2=plot(nflipsnew,newheadslowerm,'k--','Linewidth',1.5);
plot(nflipsnew,newheadsupperm,'k--','Linewidth',1.5);
hline3=plot(nflipsnew,newheadsbayesmeanm,'b','Linewidth',2);
hline4=plot(nflipsnew,newheadsbayeslowerm,'b--');
plot(nflipsnew,newheadsbayesupperm,'b--');
hline5=plot(nflipsnew,newheadsnullmeanm,'r','Linewidth',2);
hline6=plot(nflipsnew,newheadsnullmeanrejm,'m','Linewidth',2);
newheadsnullrejlowerm=mean(newheadsnullrejlower');
newheadsnullrejupperm=mean(newheadsnullrejupper');
plot(nflipsnew,newheadsnullrejlowerm,'m--');
plot(nflipsnew,newheadsnullrejupperm,'m--');
ylabel('{\itN+}');
set(gca,'YLim',[0 28],'FontName','Courier'); grid on;
set(gca,'XLim',[0 32],'FontName','Courier'); grid on;
legend([hline1,hline3,hline5,hline6],'actual','Bayesian','null','adjusted null');
%legend([hline5,hline1,hline2,hline3,hline4],'prior','intermediate (median)','intermediate (upper/lower)', 'final (median)','final(upper/lower');
%
subplot(2,3,2)
load Runs/Biased100Real100Site0pt6
hold
title('Biased, \theta_0=0.6')
set(gca,'FontSize',14);
plot(nflipsnew,newheadsmeanm,'r','Linewidth',2);
plot(nflipsnew,newheadslowerm,'k--','Linewidth',1.5);
plot(nflipsnew,newheadsupperm,'k--','Linewidth',1.5);
plot(nflipsnew,newheadsbayesmeanm,'b','Linewidth',2);
plot(nflipsnew,newheadsbayeslowerm,'b--');
plot(nflipsnew,newheadsbayesupperm,'b--');
plot(nflipsnew,newheadsnullmeanm,'r','Linewidth',2);
plot(nflipsnew,newheadsnullmeanrejm,'m','Linewidth',2);
newheadsnullrejlowerm=mean(newheadsnullrejlower');
newheadsnullrejupperm=mean(newheadsnullrejupper');
plot(nflipsnew,newheadsnullrejlowerm,'m--');
plot(nflipsnew,newheadsnullrejupperm,'m--');
set(gca,'YLim',[0 28],'FontName','Courier'); grid on;
set(gca,'XLim',[0 32],'FontName','Courier'); grid on;
%
subplot(2,3,3)
load Runs/Biased100Real100Site0pt75
hold
title('Biased, \theta_0=0.75')
set(gca,'FontSize',14);
plot(nflipsnew,newheadsmeanm,'r','Linewidth',2);
plot(nflipsnew,newheadslowerm,'k--','Linewidth',1.5);
plot(nflipsnew,newheadsupperm,'k--','Linewidth',1.5);
plot(nflipsnew,newheadsbayesmeanm,'b','Linewidth',2);
plot(nflipsnew,newheadsbayeslowerm,'b--');
plot(nflipsnew,newheadsbayesupperm,'b--');
plot(nflipsnew,newheadsnullmeanm,'r','Linewidth',2);
plot(nflipsnew,newheadsnullmeanrejm,'m','Linewidth',2);
newheadsnullrejlowerm=mean(newheadsnullrejlower');
newheadsnullrejupperm=mean(newheadsnullrejupper');
plot(nflipsnew,newheadsnullrejlowerm,'m--');
plot(nflipsnew,newheadsnullrejupperm,'m--');
set(gca,'YLim',[0 28],'FontName','Courier'); grid on;
set(gca,'XLim',[0 32],'FontName','Courier'); grid on;
%
subplot(2,3,4)
load Runs/Unbiased100Real5Site
hold
set(gca,'FontSize',14);
plot(nflipsnew,newheadsmeanm,'r','Linewidth',2);
plot(nflipsnew,newheadslowerm,'k--','Linewidth',1.5);
plot(nflipsnew,newheadsupperm,'k--','Linewidth',1.5);
plot(nflipsnew,newheadsbayesmeanm,'b','Linewidth',2);
plot(nflipsnew,newheadsbayeslowerm,'b--');
plot(nflipsnew,newheadsbayesupperm,'b--');
plot(nflipsnew,newheadsnullmeanm,'r','Linewidth',2);
plot(nflipsnew,newheadsnullmeanrejm,'m','Linewidth',2);
newheadsnullrejlowerm=mean(newheadsnullrejlower');
newheadsnullrejupperm=mean(newheadsnullrejupper');
plot(nflipsnew,newheadsnullrejlowerm,'m--');
plot(nflipsnew,newheadsnullrejupperm,'m--');
xlabel('{\itN} '); ylabel('{\itN+}');
set(gca,'YLim',[0 28],'FontName','Courier'); grid on;
set(gca,'XLim',[0 32],'FontName','Courier'); grid on;
%
subplot(2,3,5)
load Runs/Biased100Real5Site0pt6
hold
set(gca,'FontSize',14);
plot(nflipsnew,newheadsmeanm,'r','Linewidth',2);
plot(nflipsnew,newheadslowerm,'k--','Linewidth',1.5);
plot(nflipsnew,newheadsupperm,'k--','Linewidth',1.5);
plot(nflipsnew,newheadsbayesmeanm,'b','Linewidth',2);
plot(nflipsnew,newheadsbayeslowerm,'b--');
plot(nflipsnew,newheadsbayesupperm,'b--');
plot(nflipsnew,newheadsnullmeanm,'r','Linewidth',2);
plot(nflipsnew,newheadsnullmeanrejm,'m','Linewidth',2);
newheadsnullrejlowerm=mean(newheadsnullrejlower');
newheadsnullrejupperm=mean(newheadsnullrejupper');
plot(nflipsnew,newheadsnullrejlowerm,'m--');
plot(nflipsnew,newheadsnullrejupperm,'m--');
xlabel('{\itN} '); 
set(gca,'YLim',[0 28],'FontName','Courier'); grid on;
set(gca,'XLim',[0 32],'FontName','Courier'); grid on;
%
subplot(2,3,6)
load Runs/Biased100Real5Site0pt75
hold
set(gca,'FontSize',14);
plot(nflipsnew,newheadsmeanm,'r','Linewidth',2);
plot(nflipsnew,newheadslowerm,'k--','Linewidth',1.5);
plot(nflipsnew,newheadsupperm,'k--','Linewidth',1.5);
plot(nflipsnew,newheadsbayesmeanm,'b','Linewidth',2);
plot(nflipsnew,newheadsbayeslowerm,'b--');
plot(nflipsnew,newheadsbayesupperm,'b--');
plot(nflipsnew,newheadsnullmeanm,'r','Linewidth',2);
plot(nflipsnew,newheadsnullmeanrejm,'m','Linewidth',2);
newheadsnullrejlowerm=mean(newheadsnullrejlower');
newheadsnullrejupperm=mean(newheadsnullrejupper');
plot(nflipsnew,newheadsnullrejlowerm,'m--');
plot(nflipsnew,newheadsnullrejupperm,'m--');
xlabel('{\itN} '); 
set(gca,'YLim',[0 28],'FontName','Courier'); grid on;
set(gca,'XLim',[0 32],'FontName','Courier'); grid on;
%
