freq n;
read(file='h:hetero.xls');
print y x;
graph x y;
olsq y c x;
ste=@res/@s;
graph x ste;
yy=y/x;
xx=1/x;
olsq yy c xx;
ste=@res/@s;
graph xx ste;
end;