Thank you, Koen. That was very helpful. You are right, my reference was Female group, not Male. I made an error in my post. I never thought about trying a new variable for sex*age. It was interesting to see that I had different result for this variable to the one where I use param=glm and compare all sex*age interaction groups to Female*12-17. I ended up using proc logistic and oddsratio for age_grps AT (gender='1') to figure out if that result matches my proc genmod output for gender*age_grps. It provided ORs for age_grps xx vs 12-17 at gender=1 and they matched my estimates from proc genmod. That means to me that my estimations for interaction in proc genmod are for Males in each age group vs Males in 12-17.
... View more
