true=[ ...
    0.5820
    0.5020
    0.3510
    0.2660
   -1.5210
   -1.0630
    0.8810
   -1.8770
   -3.5240
    0.758
   27.034
  113.423
    0.4520
    0.0340
    0.0890
    0.3580
    0.6910
    0.1910
    0.1900
    0.2750
    0.7200
             1.920
             2.400
             6.922
             0.698
             0.690
             0.690
             1.000
             0.891
             0.906
             0.891
             1.000
             0.922
             0.906
             0.922
             1.000
             0.325
             0.021
             0.105
             0.471
             0.668
             0.295
             0.232
             0.332
             0.623];



load ejunk.dat; s11=ejunk;
load call11.stat.1; s11=call11;
load callss.stat.2; ss=callss;
load callfiml.stat.1; sfa=callfiml;
load callfiml.stat.2; sfb=callfiml;
load callfiml.stat.3; sfc=callfiml;

s11 = s11(:,[1:21,31:54]);
X11 = zeros(45,100);
Y11 = zeros(45,100);
Xss = zeros(45,100);
Yss = zeros(45,100);
Xfa = zeros(45,100);
Xfb = zeros(45,100);
Xfc = zeros(45,100);
Yfa = zeros(45,100);
Yfb = zeros(45,100);
Yfc = zeros(45,100);
for i=1:45
[x11,y11,c11,n11]=plotpdf(s11(:,i));
[xss,yss,css,nss]=plotpdf(ss(:,i));
[xfa,yfa]=plotpdf(sfa(:,i));
[xfb,yfb]=plotpdf(sfb(:,i));
[xfc,yfc]=plotpdf(sfc(:,i));
%plot(x11,y11,xss,yss,xfa,yfa,xfb,yfb,xfc,yfc, ...
%     [true(i),true(i)],[0,max(yfc)],'k:') 
%title(['Statistic ',int2str(i)])
X11(i,:)=x11;
Y11(i,:)=y11;
Xss(i,:)=xss;
Yss(i,:)=yss;
Xfa(i,:)=xfa;
Yfa(i,:)=yfa;
Xfb(i,:)=xfb;
Yfb(i,:)=yfb;
Xfc(i,:)=xfc;
Yfc(i,:)=yfc;
end;
