I have been trying to do logistic regression with interactions. Since this would have required a lot of dummy coding in proc logistic, I used proc genmod.
proc genmod;
class ahsex(ref=first) ah66a qtype &smok_ &alco_;
model jncht= wstrs wstrs*ahsex ahsex ah66a ahage &smok_ &alco_ WLTHINDF QTYPE/ error=bin link=logit type3;
run;
Another way to encode the interaction term and the main effects would be using "wstrs|ahsex". This equals "wstrs wstrs*ahsex ahsex".
The ref option could be replaced with (ref="1") . Specifying more than one REF= variable option in the CLASS statement could be a problem.
In case you are wondering, &smok_ refer to the macro variable I created earlier in the code as follows:
%let &sm0k= smok1 smok2 smok3 smokm;
Showing posts with label genmod logistic interaction dummy macro variable class. Show all posts
Showing posts with label genmod logistic interaction dummy macro variable class. Show all posts
Saturday, March 31, 2007
Subscribe to:
Posts (Atom)