Stata
Products Purchase Support Company
Search
   >> Home >> Products >> Capabilities >> Count variables >> Poisson regression Bookmark and Share

Poisson regression

Stata’s poisson command fits maximum-likelihood models of the number of occurrences (counts) of an event. In a Poisson regression model, the incidence rate for the jth observation is assumed to be given by

   r_j = exp(b_0 + b_1*x_(1,j) + ... + b_k*x_(k,j)

If E_j is the exposure, the expected number of events C_j will be

   C_j = E_j * r_j
       = exp[ ln(E_j) + b_0 + b_1*x_(1,j) + ... + b_k*x_(k,j) ]

This is the model fitted by poisson. E_j may be specified or, if not specified, is assumed to be 1.

. poisson deaths smokes i.agecat, exposure(pyears) irr

Iteration 0:   log likelihood = -33.823284  
Iteration 1:   log likelihood = -33.600471  
Iteration 2:   log likelihood = -33.600153  
Iteration 3:   log likelihood = -33.600153  

Poisson regression                                Number of obs   =         10
                                                  LR chi2(5)      =     922.93
                                                  Prob > chi2     =     0.0000
Log likelihood = -33.600153                       Pseudo R2       =     0.9321

------------------------------------------------------------------------------
      deaths |        IRR   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
      smokes |   1.425519   .1530638     3.30   0.001     1.154984    1.759421
             |
      agecat |
          2  |   4.410584   .8605197     7.61   0.000     3.009011    6.464997
          3  |    13.8392   2.542638    14.30   0.000     9.654328    19.83809
          4  |   28.51678   5.269878    18.13   0.000     19.85177    40.96395
          5  |   40.45121   7.775511    19.25   0.000     27.75326    58.95885
      pyears | (exposure)
------------------------------------------------------------------------------

The syntax of all estimation commands is the same: the name of the dependent variable is followed by the names of the independent variables, which are followed by a comma and any options. In this case, we controlled for the exposure (person-years recorded in the variable pyears) and asked that results be displayed as incidence-rate ratios rather than as coefficients.

Stata can perform stepwise Poisson regression, too. Also, Stata provides Cox regression, exponential, Weibull, and other parametric survival models, as well as logistic regression and stepwise variants of these procedures.

See New in Stata 11 for more about what was added in Stata Release 11.

Stata 11
Overview: Why use Stata?
Stata/MP
64-bit Stata
Capabilities
Overview
Data management
Graphics
Basic statistics
Linear models
Binary and discrete outcomes
Panel data
Survey methods
Time series
Survival analysis
Epidemiology tools
Mixed models
GLM
ANOVA / MANOVA
Multiple imputation
Exact statistics
Nonparametric methods
Multivariate methods
Cluster analysis
Resampling
Model testing
Maximum likelihood
Other statistical methods
Programming
Matrix programming—Mata
Internet capabilities
Accessibility
Sample session
User-written commands
New in Stata 11
Supported platforms
Which Stata package?
Technical support
User comments
Products
Stata 11
Order Stata
Upgrade
Training
Bookstore
Stata Journal
Stata Press
Stata News
STB
Stat/Transfer
Gift Shop

Site overview
Products
Resources & support
Company
Site index

© Copyright 1996–2009 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   Site index