data a;
infile 'ipums_2009.dat';
input hhwt 1-4  perwt 5-8  age 9-11 sex 12  race 13
incearn 14-19  poverty 20-22;
if incearn lt 1 then delete;
data b;
   set a;
l_earn = log(incearn);
if race eq 1 then white=1; else white=0;
a2 = age*age;
w_age = white*age;
w_age2 = white*a2;
proc means;
proc reg;
   model l_earn= white age a2 w_age w_age2;
run;