OPTIONS PS=40 lS=64;

data a;
infile 'nym';
input nn ST CO EXPEN WEALTH POP PINTERG DENS INCOME ID GROWR ;
lex = log(expen);
lw = log(wealth);
lpop = log(pop);
ldens = log(dens);
lincome = log(income);
pint2 = pinterg**2;
pint3 = pinterg**3;
lpop2 = lpop**2;
ldens2 = ldens**2;
lpop3 = lpop**3;
ldens3 = ldens**3;

if growr < 0 then lgrowr = - log(-growr);
if growr > 0 then lgrowr =  log(growr);
if _N_ ne 887;
if _N_ ne 475;
run;


run;

proc univariate plot;
var  ST CO EXPEN WEALTH POP PINTERG DENS INCOME ID GROWR
    lex lw lpop ldens lincome;

proc plot;
plot  EXPEN*(growr wealth pop dens income PINTERG GROWR);
plot lex*(lgrowr lw lpop ldens lincome);
run;

proc reg;
model lEX = lW lPOP lpop2 lpop3 PINTERG pint2 
       pint3 lDENS ldens2 ldens3 lINCOME lGROWR/ P R;
run;

proc reg data = a;
model lEX = lW lPOP lpop2 lpop3 PINTERG pint2 pint3 lDENS 
       ldens2 ldens3 lINCOME lGROWR/ METHOD=stepwise;

proc reg data = a;
model lEX = lW lPOP lpop2 lpop3 PINTERG pint2 pint3 lDENS 
       ldens2 ldens3 lINCOME lGROWR/ METHOD=backward;

proc reg data = a;
model lEX = lW lPOP lpop2 lpop3 PINTERG pint2 pint3 lDENS 
       ldens2 ldens3 lINCOME lGROWR/ METHOD=adjrsq;

run;