Stable computation of maximum likelihood estimates in identity link Poisson regression

View Researcher's Other Codes

Disclaimer: The provided code links for this paper are external links. Science Nest has no responsibility for the accuracy, legality or content of these links. Also, by downloading this code(s), you agree to comply with the terms of use as set out by the author(s) of the code(s).

Authors Ian C. Marschner
Journal/Conference Name Journal of Computational and Graphical Statistics
Paper Category
Paper Abstract Identity link Poisson regression is useful when the mean of a count variable depends additively on a collection of predictor variables. It is particularly important in epidemiology, for modeling absolute differences in disease incidence rates as a function of covariates. A complication of such models is that standard computational methods for maximum likelihood estimation can be numerically unstable due to the nonnegativity constraints on the Poisson means. Here we present a straightforward and flexible method that provides stable maximization of the likelihood function over the constrained parameter space. This is achieved by conducting a sequence of maximizations within subsets of the parameter space, after which the global maximum is identified from among the subset maxima. The method adapts and extends EM algorithms that are useful in specialized applications involving Poisson deconvolution, but which do not apply in more general regression contexts. As well as allowing categorical and continuous cova...
Date of publication 2010
Code Programming Language R

Copyright Researcher 2022