turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

Find a Community

- Home
- /
- Analytics
- /
- Stat Procs
- /
- Confidence intervals for rates using Poisson distr...

Topic Options

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

10-29-2010 01:57 PM

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.

Thanks.

Thanks.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

10-29-2010 03:58 PM

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.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

10-29-2010 06:14 PM

We are using birth and population data to calculate birth rates. The confidence intervals will be used to compare rates by time intervals and locations.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

10-29-2010 08:00 PM

Assuming that the data are distributed as Poisson conditional on the population size, then you can obtain confidence intervals for the Poisson rate using the GENMOD procedure as follows:

proc genmod data=mydata;

model birth_count = / dist=poisson offset=log_PopSize;

estimate "log(rate)" intercept 1;

run;

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.

proc genmod data=mydata;

model birth_count = / dist=poisson offset=log_PopSize;

estimate "log(rate)" intercept 1;

run;

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.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

11-01-2010 03:20 PM

See this usage note that discusses the modeling or rates and computing rate estimates and confidence limits:

http://support.sas.com/kb/24188

http://support.sas.com/kb/24188