Hello,
I am currently working on a project where I want to see the OR of having a cardiac event for: underweight, overweight and obese bmi categories at each of the age values I listed in my code (ages 50 & 60). For some reason my code below only outputs the OR for age and bmicat individually and the OR for age at 50 and 60. Is there something with my code? How can I get the output I want for the odds ratio(listed below)?
Code used:
proc logistic data=final;
class bmicat(ref="normal")/param=ref;
effect agespline=spline(age/naturalcubic); *Wanted to allow for curvature by using spline;
model cardiac(event="yes")= agespline bmicat/ rsq ;
oddsratio age/at (age= 50 60);
oddsratio bmicat/at (age= 50 60);
run;
Output I am trying to get for the Odds Ratios:
underweight vs. normal at age=50
overweight vs. normal at age= 50
obese vs. normal at age=50
underweight vs. normal at age=60
overweight vs. normal at age= 60
obese vs. normal at age=60