[R] R and Poisson
Paul Mcgeoghan
Mcgeoghan at Cardiff.ac.uk
Mon Aug 18 20:46:41 CEST 2003
Hi, I wonder if anyone can answer the following or point me in the direction of
how to obtain answers to the questions. Below is Output from R and further down
are the questions raised and explanation of the study.
Output from R:
glm(formula = CB95TO00 ~ URB + INC, family = poisson)
Deviance Residuals:
Min 1Q Median 3Q Max
-1.2272 -1.1290 0.2709 0.4272 2.1376
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -0.30621 0.13499 -2.268 0.0233 *
URB2 0.02253 0.16826 0.134 0.8935
URB3 -0.00936 0.15263 -0.061 0.9511
INC2 -0.14430 0.12342 -1.169 0.2423
INC3 -0.55092 0.31351 -1.757 0.0789 .
---
Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1
(Dispersion parameter for poisson family taken to be 1)
Null deviance: 403.97 on 420 degrees of freedom
Residual deviance: 399.61 on 416 degrees of freedom
AIC: 883.8
Number of Fisher Scoring iterations: 4
Explanation and Questions raised.
The dependent variable is:
Number of children born in last 5 years: (values range from 0 to 3).
Distribution of dependent variable (named CB95TO00)
0 203
1 157
2 59
3 2
Predictors are:
Level of Urbanisation (3 categories 1: Rural; 2:Semi-Urban; 3: Urban)
Income Level (3 categories: 1: Low; 2:Medium; 3: High)
The questions are (1) how does one interpret the coefficients in the output:
Our interpretation is Urb2 compared to Urb1 gives an estimate of .02253;
Urb3 compared to Urb1 gives a parameter estimate of -.00936 etc. Neither of
these shows significance. How does one interpret this exactly with regards to
the dependent variable which is Number of children?
2) How does one interpret the intercept which shows significance?
3) What does the Null Deviance tell us and the Residual Deviance?
4) What does the AIC tell us?
5) Is it possible to obtain goodness of fit statistics such as Pearson
ChiSquare and Log-Likelihood similar to what SAS statistical software gives?
6) Is it possible to find out if Urbanisation and Income are significant
overall in R?
Thanks in advance for any assistance,
Regards,
Paul
==================
Paul McGeoghan,
Application support specialist (Statistics and Databases),
Information Services,
Cardiff University.
Tel. 02920 (875035).
More information about the R-help
mailing list