options ls=77 ps=200 nonumber label;
data fl;
infile 'fldat1.sasdat';
input co lat lon npop whit blac hisp o65 hsed coll inco bush gore 
brow nade harr hage buch mcre phil moor;
ydot=127.2651;
tvot=bush+gore+brow+nade+harr+hage+mcre+phil+moor;
pbush=100*bush/tvot;
pgore=100*gore/tvot;
pbrow=100*brow/tvot;
pnade=100*nade/tvot;
pbuch=100*buch/tvot;
inco=inco/1000;
lblac=log(blac);
lhisp=log(hisp);
lpop=log(npop);
y1=buch;
y2=2*sqrt(ydot*buch);
y3=ydot*log(buch/tvot);
y4=3*exp(2*log(ydot)/3)*exp(log(buch)/3);
tvot3=exp(log(tvot)/3);
lpop1=tvot*lpop;
whit1=tvot*whit;
lblac1=tvot*lblac;
lhisp1=tvot*lhisp;
o651=tvot*o65;
hsed1=tvot*hsed;
coll1=tvot*coll;
inco1=tvot*inco;
pbush1=tvot*pbush;
pbrow1=tvot*pbrow;
pnade1=tvot*pnade;
lpop2=sqrt(tvot)*lpop;
whit2=sqrt(tvot)*whit;
lblac2=sqrt(tvot)*lblac;
lhisp2=sqrt(tvot)*lhisp;
o652=sqrt(tvot)*o65;
hsed2=sqrt(tvot)*hsed;
coll2=sqrt(tvot)*coll;
inco2=sqrt(tvot)*inco;
pbush2=sqrt(tvot)*pbush;
pbrow2=sqrt(tvot)*pbrow;
pnade2=sqrt(tvot)*pnade;
lpop3=tvot3*lpop;
whit3=tvot3*whit;
lblac3=tvot3*lblac;
lhisp3=tvot3*lhisp;
o653=tvot3*o65;
hsed3=tvot3*hsed;
coll3=tvot3*coll;
inco3=tvot3*inco;
pbush3=tvot3*pbush;
pbrow3=tvot3*pbrow;
pnade3=tvot3*pnade;
run;
;
proc reg;
model y1=lpop1 whit1 lblac1 lhisp1 o651 hsed1 coll1 inco1 pbush1 pbrow1 pnade1 / 
selection=cp;
run;
;
proc reg;
model y1=lpop1 whit1 lblac1 lhisp1 o651 hsed1 coll1 inco1 pbush1 pbrow1 pnade1 / 
selection=backward slstay=0.05;
run;
;
proc reg;
model y2=lpop2 whit2 lblac2 lhisp2 o652 hsed2 coll2 inco2 pbush2 pbrow2 pnade2 / 
selection=cp;
run;
;
proc reg;
model y2=lpop2 whit2 lblac2 lhisp2 o652 hsed2 coll2 inco2 pbush2 pbrow2 pnade2 / 
selection=backward slstay=0.05;
run;
;
proc reg;
model y3=lpop whit lblac lhisp o65 hsed coll inco pbush pbrow pnade / 
selection=cp;
run;
;
proc reg;
model y3=lpop whit lblac lhisp o65 hsed coll inco pbush pbrow pnade / 
selection=backward slstay=0.05;
run;
;
proc reg;
model y4=lpop3 whit3 lblac3 lhisp3 o653 hsed3 coll3 inco3 pbush3 pbrow3 pnade3 / 
selection=cp;
run;
;
proc reg;
model y4=lpop3 whit3 lblac3 lhisp3 o653 hsed3 coll3 inco3 pbush3 pbrow3 pnade3 / 
selection=backward slstay=0.05;
run;
;