Hello, This is my first time writing here. I have mortality data from a cohort study. I've been trying to calculate incidence rates for mortality with no success. I used proc genmod with lsmeans to estimate the IR using Poisson, but the results are very very different from those obtained by my colleague using STATA.This is the code I've been using:   proc genmod data=muj; class muerto diabetes(ref=first) educacion(ref=first) sexo(ref=last); model muerto=edad educacion  diabetes/dist=poisson                         link=log                         offset=logpyears; lsmeans diabetes/ilink exp cl diff OM; run;  The stata code used for this is the following:   poisson muerto i.diabetes i.cat_edad i.educacion, exposure(pyears) irr margins diabetes, predict(ir)   Is there some procedure or code that gives me adjusted incidence rates? Thanks to all in advance ... View more