A slightly different approach would be to use the n/d format in GENMOD. You could enter the "rate", numerator/denominator for each year, so that the model statement would be: model n/d=year/d=p etc. The output will give you the calculated poisson rates along with standard errors for each year. Then you could compare the yearly rates using LSMEANS.
... View more