options crt;  ? 4e1
in waugh;
smpl 1,200;
olsq price c green nostalks disperse;  ? (a)
msd(cova) price green nostalks disperse;  ? (b,c)
mat cov200 = @cova*(199/200);  ? Apparently, Waugh divided by 200 (vs. 199)
print cov200;  ? (c)
?
set unitch = 2.782*(@coef(2)+@coef(3)+@coef(4));
set unitchw = 2.782*(.13826-1.53394-.27553);
print unitch unitchw;  ? (d)
?
?  Make coefficients from Waugh's covariance matrix
?  (much closer to his published answers, so his errors
?   were primarily in computing the data covariance matrix)
?  It's not a data error, because at least one number from
?  every row is correct to 4 digits or more.
?  (adding up 200 products by hand for each element
?    must have been tedious...)
?
read(nrow=3,ncol=1) wxy;    ? first row of his matrix, except for (1,1)
3430.89  -100.92   -82.35;
read(type=sym,nrow=3) wxx;  ? last 3 rows of his matrix, transposed
24317.19
  -17.01   61.33
 -154.54   25.51    83.07;
mat wbeta = wxx"wxy;
print wbeta;  ? (e)