I am trying to calculate Poisson-based 95% confidence intervals for rates. Is there a way to do this in SAS? The rates are based on small numbers of events, so the standard normal-based 95% CIs aren't appropriate.
Would you describe your problem with a little more detail? You speak of rates, so that leads me to believe that you have differing exposure time for different observations. But it would be good to know more precisely what data you have and what you want to estimate.
where log_PopSize is the (natural) log of the population size and birth_count should be self explanatory. The above code would provide the rate for a population of size 1. You probably don't want to specify the rate for a population of size 1. In order to specify the rate for a population of size K, compute log_PopSize as
log_PopSize = log(PopSize / K);
Of course, the variable log_PopSize needs to be constructed before executing the GENMOD procedure. You might also want to examine whether a negative binomial distribution specification provides a better fit than the Poisson.