/*
** Lesson 3.2: Residual Analysis
*/
use gpe2; @ using GPE package (version 2) @
output file = gpe\output3.2 reset;
load data[17,7] = gpe\longley.txt;
data = data[2:17,.];
    
PGNP = data[.,2];
GNP = data[.,3]/1000;               
EM = data[.,7]/1000;
RGNP = 100*GNP./PGNP;

call reset; 	@ initialize global variables @
_rstat = 1;	@ report residual statistics @
_rlist = 1;	@ list residuals @
_rplot = 1;	@ plot residuals @
_names = {"EM","RGNP"};
call estimate(EM,RGNP);
end;