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 call11.stat.1; s11=call11;
load callss.stat.1; sl=callss;
load callss.stat.2; st=callss;
load callfiml.stat.1; sfl=callfiml;
load callfiml.stat.2; sfm=callfiml;
load callfiml.stat.3; sft=callfiml;
s11 = s11(:,[1:21,31:54]);
for i=1:45

  disp(['Statistic',int2str(i)])
  disp(['True = ',num2str(true(i))])
  x=sort(s11(:,i));
  lb11 =x(25);
  ub11 =x(975);
  x=sort(sl(:,i));
  lbsl =x(25);
  ubsl =x(975);
  x=sort(st(:,i));
  lbst =x(25);
  ubst =x(975);
  x=sort(sfl(:,i));
  lbfl =x(25);
  ubfl =x(975);
  x=sort(sfm(:,i));
  lbfm =x(25);
  ubfm =x(975);
  x=sort(sft(:,i));
  lbft =x(25);
  ubft =x(975);
  disp(sprintf('SS-R-T:  %7.3f %7.3f %7.3f ',lbft,mean(sft(:,i)),ubft))
  disp(sprintf('SS-R-M:  %7.3f %7.3f %7.3f ',lbfm,mean(sfm(:,i)),ubfm))
  disp(sprintf('SS-R-L:  %7.3f %7.3f %7.3f ',lbfl,mean(sfl(:,i)),ubfl))
  disp(sprintf('SS-UR-T: %7.3f %7.3f %7.3f ',lbst,mean(st(:,i)),ubst))
  disp(sprintf('SS-UR-L: %7.3f %7.3f %7.3f ',lbsl,mean(sl(:,i)),ubsl))
  disp(sprintf('VARMA:   %7.3f %7.3f %7.3f ',lb11,mean(s11(:,i)),ub11))


end;



