After 4 days trying, I finally know how to use SAS to figure out which variable is associate with specific variable.

you can use following SAS program to figure out which variable between HAC1A~HAC1O
are associate with HFE7;
--------------------------
proc glm data=adult;
class HAC1A HAC1B HAC1D HAC1E HAC1F HAC1G HAC1H HAC1I HAC1J HAC1K HAC1L HAC1M HAC1N HAC1O;
model HAC1A HAC1B HAC1D HAC1E HAC1F HAC1G HAC1H HAC1I HAC1J HAC1K HAC1L HAC1M HAC1N HAC1O=HFE7;
run;
--------------------------


after you find something's P<0.05 (ex:HAC1E)

you can use following SAS program to see
what's the different between group 1 and 2 in HFE7
--------------------------
proc glm data=adult;
class HFE7;
model HAC1E=HFE7;
lsmeans HFE7;
run;
--------------------------