/*
** Lesson 3.6: Residual Diagnostics
*/
use gpe2;
output file = gpe\output3.6;
load data[40,6] = gpe\cjx.txt;

year = data[2:40,1];
X = ln(data[2:40,2]);
L = ln(data[2:40,3]);
K = ln(data[2:40,5]);
names = {"X", "L", "K"};

call reset;
_rstat = 1; 
_rlist = 2;	@ check influential obs. @
_bjtest = 1;	@ normality test @

call estimate(X,L~K);
end;