options crt;  ? 3e1
in poly;
freq a;
smpl 60,72;
lncp = log(cump);
lnuc = log(ucostp);
lny = log(prodp);
?graph ucostp cump;
graph(preview) cump ucostp;  ? graphics version
?graph lnuc lncp;  ? (a)
graph(preview) lncp lnuc;  ? (a) - graphics version
olsq lnuc c lncp;
frml eqd d = 2**lncp;
analyz eqd;       ? (b) - d and its standard error
?
in tio2;
freq a;
smpl 55,70;
lncp = log(cumt);
ucostd = ucostt/defl;
lnuc = log(ucostd);
lny = log(prodt);
?graph ucostd cumt;
graph(preview) cumt ucostd;  ? graphics version
?graph lnuc lncp;  ? (a)
graph(preview) lncp lnuc;  ? (a) - graphics version
olsq lnuc c lncp;
frml eqd d = 2**lncp;  ? We don't actually have to repeat this FRML from above
analyz eqd;       ? (b)