[R] Negative Binomial Model

From: Mike Ryckman <ryckman_at_arizona.edu>
Date: Wed, 14 May 2008 02:01:52 -0700


I am trying to run a negative binomial regression model in R and can't get the
standard errors to match the output I get from the Stata nbreg command. I've tried a few different options but haven't had much luck. The closest I've found

gamlss(formula, family = NBI, sigma.formula = ~ 1,data=dataframe)

...But this is still a little off most of the time and pretty far off at other
times (compared with the Stata output). The glm.nb from the MASS package produces the correct coefficients, but different (usually very different) standard errors.

Could anybody explain this and point me in the right direction? I'd really appreciate it.

Thank you,

Mike Ryckman
Department of Political Science
University of Arizona

R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Wed 14 May 2008 - 10:56:11 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Wed 14 May 2008 - 12:30:37 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive