/*   This is file Summary.GAU: Read indexes from disk and print summary table and figures.  
	It reads actual expenditure as well as Max and Min Laspeyres, EKS, CCD and Geary 
	         real incomes (created by Basics.gau) from disk.  
	It also reads GAIA real incomes for all 3 systems (created by GaiaQuai.GAU) from disk.
         Then it calculates all correlations etc. and prints summary table and Figures 4-6.  

(c) J.P. Neary 1999

*/
new;  /*************************************************************************/
/*          USER-SPECIFIED VARIABLES: ECONOMIC DATA AND PARAMETERS  */
n=11;                       /* n = # goods        */
nyear=1980;
mrel=60;             /* Index of country relative to which the results are printed */
/*************************************************************************/
/*    USER-SPECIFIED VARIABLES WHICH CONTROL OPERATION OF PROGRAM  */
library gearylib, pgraph;
printifS=1;         /* printifS=1 for printing of summary table */
GRAPHIF=1;    /* graphif=0: no graph; >1: print graphs */
/*************************************************************************/
/*  Retrieve m-by-n matrices of ZGAIA from disk: Location of files depends on demand system */
load path = c:\me\research\geary2\empirics\gaussout\HAIDS ;
   load ZGAIAALL; ZG_HALL = ZGAIAALL; 
load path = c:\me\research\geary2\empirics\gaussout\AIDS ;
   load ZGAIAALL; ZG_AALL = ZGAIAALL;  
load path = c:\me\research\geary2\empirics\gaussout\QUAIDS ;
   load ZGAIAALL; ZG_QALL = ZGAIAALL; 

load path = c:\me\research\geary2\empirics\gaussout\ALL ;  /* Empirical indexes read from here */
load  z  StarMax  StarMin  EKS  CCD  ZGeary; 
/*************************************************************************/
/* These variables are not needed; but for the present they should be defined, since
      otherwise there are conflicts with GearyLib.src (where they are assumed as globals).
      Good idea to fix this sometime.   */
testif=0;
alpha=zeros(n-1,1);    /* (n-1)-by-1 vector of ML alpha's:  n'th to be estimated residually */
beta=zeros(n-1,1);      /* (n-1)-by-1 vector of ML beta's:  n'th to be estimated residually */
lambda=zeros(n-1,1);  /* (n-1)-by-1 vector of ML lambdas: n'th to be estimated residually */
ggamma=zeros(n-1,n-1);   /* (n-1)-by-(n-1) matrix of ML gamma's:  n'th to be estimated residually */
/*************************************************************************/
/*                      PRINT FINAL SUMMARY INFORMATION               */
/*   Assemble the key matrix "results" to be printed and charted  */
/* De-comment whichever of the next lines is desired:
results = z ~StarMax ~StarMin ~EKS ~CCD ~ZGeary ;               /* Table 1 and Figure 1 */
results = EKS  ~ZGeary ~ZG_HALL ~ZG_AALL ~ZG_QALL ;      /* Basis for Figure 3 .xls file */
results = EKS  ~ZGeary ~ZG_HALL ;                                        /* Table A1: All GAIA-Haids */
results = EKS  ~ZGeary ~ZG_AALL ;                                         /* Table A3: All GAIA-Aids */
results = EKS  ~ZGeary ~ZG_QALL ;                                         /* Table A5: All GAIA-Quaids */
*/

results = EKS ~CCD ~ZGeary ~ZG_HALL ~ZG_AALL ~ZG_QALL ;     /* Figures 4-6 in paper */

/*************************************************************************/
results = results./results[mrel,.];   /*  Reexpress all results relative to country indexed "mrel"  */
call printout (results, EKS, ZGeary, mrel, printifS,graphif);
/*************************************************************************/
finish:
end;