options crt;  ? 9e2
in klem;
freq a;
smpl 47,71;
lnqvql = log(qv/ql);
lnqkql = log(qk/ql);
olsq lnqvql c lnqkql;  ? (a)

lnqvqk = log(qv/qk);
lnpkpv = log(pk/pv);
lnplpv = log(pl/pv);
lnpkpl = log(pk/pl);
olsq lnqvqk c lnpkpv;
  set sigma1 = @coef(2);
  set sigma4 = @coef(2)/@rsq;  ? for (c)
  copy @rsq rsq1;
  set dft = @nob-2;
  proc test1;
    set t1 = (@coef(2)-1)/@ses(2);
    cdf(t,df=dft) t1;
  endproc;
  test1;
olsq lnqvql c lnplpv;
  set sigma2 = @coef(2);
  set sigma5 = @coef(2)/@rsq;
  test1;
lnpkplm = - lnpkpl;
olsq lnqkql c lnpkplm;  ? reverse sign of coef
  set sigma3 = @coef(2);
  set sigma6 = @coef(2)/@rsq;
  test1;
print sigma1-sigma3;  ? (b)

print sigma4-sigma6;
olsq lnpkpv c lnqvqk;
print rsq1 @rsq;
set sig4 = 1/@coef(2);
print sigma4 sig4;  ? (c)