% Copyright (C) 2001 Michel Juillard
%
function disp_moments(y,var_list)
  global lgy_ endo_nbr ykmin_ options_
  
  warning off
  nvar = size(var_list,1);
  if nvar == 0
    nvar = endo_nbr;
    ivar = [1:nvar]';
  else
    ivar=zeros(nvar,1);
    for i=1:nvar
      i_tmp = strmatch(var_list(i,:),lgy_,'exact');
      if isempty(i_tmp)
      	error (['One of the variable specified does not exist']) ;
      else
	ivar(i) = i_tmp;
      end
    end
  end
  y = y(ivar,options_.drop+ykmin_+1:end)';
  
  m = mean(y);
  y = y - repmat(m,size(y,1),1);
  s2 = mean(y.*y);
  s = sqrt(s2);
  
  labels = deblank(lgy_(ivar,:));
  
  if options_.nomoments == 0
    z = [ m' s' s2' (mean(y.^3)./s2.^1.5)' (mean(y.^4)./(s2.*s2)-3)' ];
  
    title='MOMENTS OF SIMULATED VARIABLES';
    if options_.hp_filter == 1
      title = [title ' (HP filter, lambda = ' ...
	       int2str(options_.hp_filter) ')'];
    end
    headers=strvcat('VARIABLE','MEAN','STD. DEV.','VARIANCE','SKEWNESS', ...
		    'KURTOSIS');
    table(title,headers,labels,z,size(labels,2)+2,16,6);
  end
  
  if options_.nocorr == 0
    corr = (y'*y/size(y,1))./(s'*s);
    title = 'CORRELATION OF SIMULATED VARIABLES';
    if options_.hp_filter == 1
      title = [title ' (HP filter, lambda = ' ...
	       int2str(options_.hp_filter) ')'];
    end
    headers = strvcat('VARIABLE',lgy_(ivar,:));
    table(title,headers,labels,corr,size(labels,2)+2,8,4);
  end
  
  ar = options_.ar;
  options_ = set_default_option(options_,'ar',5);
  ar = options_.ar;
  if ar > 0
    autocorr = [];
    for i=1:ar
      autocorr = [ autocorr mean(y(1:end-i,:).*y(i+1:end,:))'./s2' ];
    end
    title = 'AUTOCORRELATION OF SIMULATED VARIABLES';
    if options_.hp_filter == 1
      title = [title ' (HP filter, lambda = ' ...
	       int2str(options_.hp_filter) ')'];
    end
    headers = strvcat('VARIABLE',int2str([1:ar]'));
    table(title,headers,labels,autocorr,size(labels,2)+2,8,4);
  end
  
  warning on
% 10/03/02 MJ corrected order std. dev var in printed report.
% 01/02/03 MJ added correlation and autocorrelation
% 01/19/03 MJ corrected variable name truncation
% 02/18/03 MJ added subtitle for HP filter
% 03/02/03 MJ added ykmin_ to the number of entries of y
% 04/28/03 MJ modified handling of options_
% 06/23/03 MJ added warning off
