I have TB case surveillance data within a county from 2006-2012. I want to run a poisson regression and calculate incidence rates to determine which years had the highest incident rates and to identify any years with large outbreaks. I just want to make sure I am doing this right in sas, I was told to use the estimate statement and have never used it before. This is my code: Proc genmod data = TB;
class year(param=ref ref='2006');
model CASES = YEAR/ dist=poisson OFFSET=LOGN ;
estimate '2006' year 1 0 0 0 0 0 0 0 0 0 0 0/ exp;
estimate '2007' year 0 1 0 0 0 0 0 0 0 0 0 0/ exp;
estimate '2008' year 0 0 1 0 0 0 0 0 0 0 0 0/ exp;
estimate '2009' year 0 0 0 1 0 0 0 0 0 0 0 0 / exp;
estimate '2010' year 0 0 0 0 1 0 0 0 0 0 0 0 / exp;
estimate '2011' year 0 0 0 0 0 1 0 0 0 0 0 0 / exp;
estimate '2012' year 0 0 0 0 0 0 1 0 0 0 0 0/ exp;
estimate '2013' year 0 0 0 0 0 0 0 1 0 0 0 0/ exp;
estimate '2014' year 0 0 0 0 0 0 0 0 1 0 0 0 / exp;
estimate '2015' year 0 0 0 0 0 0 0 0 0 1 0 0/ exp;
estimate '2016' year 0 0 0 0 0 0 0 0 0 0 1 0 / exp;
estimate '2017' year 0 0 0 0 0 0 0 0 0 0 0 1 / exp;
run; LOGN is a variable I created earlier that is equal to LOG(population). I have a variable for population data of each year. Please let me know if I am using the estimate statement correctly? The output doesn't seem to be matching the case count numbers--for instance 2012 had the highest case count numbers but it is not the highest incidence rate... This is my output from the contrast statements Output from contrast statment Any help is appreciated! Thank you!
... View more
