proc phreg data=d;
class drug diabetes;
model time_primary_ep*primary_ep(0) = bmi drug diabetes drug*bmi / risklimits ties=EFRON;
hazardratio drug / at (bmi =20) diff=ALL;
hazardratio drug / at (bmi =25) diff=ALL;
run;