*** Rerun fla.sas to get predictions for selected models;
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);
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;
run;
;
proc reg;
model y1=lpop1 whit1 lblac1 lhisp1 o651 hsed1 inco1 pbrow1 / cli;
run;
;
proc reg;
model y1=whit1 lblac1 lhisp1 hsed1 inco1 pbrow1 / cli;
run;
;
proc reg;
model y2=lpop2 whit2 lhisp2 hsed2 inco2 pbush2 pbrow2 / cli;
run;
;
proc reg;
model y2=lpop2 whit2 lhisp2 inco2 pbush2 / cli;
run;
;
proc reg;
model y3=lpop lhisp hsed inco pbush / cli;
run;
;
proc reg;
model y3=lhisp hsed inco pbush / cli;
run;
;