? P&R Exercise 6.7

options crt, limwarn=0; 
freq q; 
smpl 47:1 95:3;
read(file='ex67.xls', format=excel) gcq gydq;
smpl 59:1 95:3;
olsq gcq c gydq gcq(-1);
print @dw @nob @ses(3);
set h=(1-(@dw/2))*sqrt(@nob/(1-@nob*(@ses(3)**2)));
print h;
end;stop;