Dear all,
I used the glm.nb with the default values from the MASS package to run a
negative binomial regression. Here is a simple example:
set.seed(123)
y <- c( rep(0, 30), rpois(70, lambda=2) )
smoke <- factor( sample( c("NO", "YES"), 100, replace=T ) )
height <- c( rnorm(30, mean=100, sd=20), rnorm(70, mean=150, sd=20) )
fit <- glm.nb( y ~ smoke + height )
coef(summary(fit))
Estimate Std. Error z value Pr(>|z|)
(Intercept) -2.34907191 0.537610710 -4.3694664 1.245505e-05
smokeYES -0.03479730 0.197627539 -0.1760751 8.602349e-01
height 0.01942373 0.003527538 5.5063142 3.664243e-08
The question now is how do I report the results, say, for height? Do I
simply take the anti logs. i.e. 1.019613 = exp(0.019423) ?
I have seen one paper where they report using anti log base 10 instead
of natural base but they use STATA though.
Please kindly advise. Thank you.
Regards, Adai
______________________________________________
R-help@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.