%Use results to construct means, standard deviations, and RMSE 
%for Tables 1-4
%

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.2; ss=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]);
disp('RESULTS FOR TABLES 1-4')
disp('----------------------')
disp(' ')
disp('           Mean   St.Dev.  RMSE')
disp('          ----------------------')

for i=1:45

  disp(['Statistic ',int2str(i)])
  disp(sprintf('True  =  %7.3f ',true(i)))
  disp(sprintf('SS-R-T:  %7.3f %7.3f %7.3f', ...
      mean(sft(:,i)),std(sft(:,i)),rmse(sft(:,i),true(i))))
  disp(sprintf('SS-R-M:  %7.3f %7.3f %7.3f', ...
      mean(sfm(:,i)),std(sfm(:,i)),rmse(sfm(:,i),true(i))))
  disp(sprintf('SS-R-L:  %7.3f %7.3f %7.3f', ...
      mean(sfl(:,i)),std(sfl(:,i)),rmse(sfl(:,i),true(i))))
  disp(sprintf('SS-UR:   %7.3f %7.3f %7.3f', ...
      mean(ss(:,i)),std(ss(:,i)),rmse(ss(:,i),true(i))))
  disp(sprintf('VARMA:   %7.3f %7.3f %7.3f', ...
      mean(s11(:,i)),std(s11(:,i)),rmse(s11(:,i),true(i))))
  disp(' ')

end;



