options crt;  ? 3e2
in poly;
freq a;
smpl 60,72;
lncp = log(cump);
lnuc = log(ucostp);
lny = log(prodp);
olsq lnuc c lncp lny;
 unmake @coef b0 b1 b2;
 rename @rsq r2m;
 rename @ssr ssrm;
 rename @fst fstm;
 set ssrym = @sdev*@sdev*(@nob-1);
 set t1 = @t(2);
 set t2 = @t(3);
frml eqr r = 1/(1+lny);
analyz eqr;  ? (a)
olsq lnuc c lncp;
 unmake @coef a0 a1;
 rename @rsq r2s;
olsq lny c lncp;
 unmake @coef d0 d1;
set a1id = b1 + d1*b2;
print a1,a1id;  ? (b)
print r2s,r2m;  ? (c)
set nob3 = @nob-3;
set fstssr = (ssrym-ssrm)/2/(ssrm/nob3);
set fstr2 = r2m/2/((1-r2m)/nob3);
print fstm fstssr fstr2 t1 t2;   ? (d)
cdf(f,df1=2,df2=nob3) fstssr;
?
in tio2;
freq a;
smpl 55,70;
lncp = log(cumt);
ucostd = ucostt/defl;
lnuc = log(ucostd);
lny = log(prodt);
olsq lnuc c lncp lny;
 unmake @coef b0 b1 b2;
 rename @rsq r2m;
 rename @ssr ssrm;
 rename @fst fstm;
 set ssrym = @sdev*@sdev*(@nob-1);
 set t1 = @t(2);
 set t2 = @t(3);
frml eqr r = 1/(1+lny);
analyz eqr;  ? (a)
olsq lnuc c lncp;
 unmake @coef a0 a1;
 rename @rsq r2s;
olsq lny c lncp;
 unmake @coef d0 d1;
set a1id = b1 + d1*b2;
print a1,a1id;  ? (b)
print r2s,r2m;  ? (c)
set nob3 = @nob-3;
set fstssr = (ssrym-ssrm)/2/(ssrm/nob3);
set fstr2 = r2m/2/((1-r2m)/nob3);
print fstm fstssr fstr2 t1 t2;   ? (d)
cdf(f,df1=2,df2=nob3) fstssr;