[R] Residual deviance (cross-post from sci.stat.consult)
Serebrenik, A.
a.serebrenik at tue.nl
Wed Jan 7 07:43:55 CET 2009
Dear all,
I'm trying to fit a statistical model to series of measurements.
Unfortunately, my knowledge of statistics is rather limited, so I'm a
bit at loss of what is going on with the model.
First of all, I've prepared a histogram. Then, I've tried to fit a Poisson model to express the relation between the middle points of
classes (mids) and the corresponding frequencies (density). I've got the following Poisson models using R:
> summary(fmDP)
Deviance Residuals:
Min 1Q Median 3Q Max
-1.20831 -0.56363 -0.28010 0.08324 3.19099
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 1.2495 0.1528 8.18 2.84e-16 ***
hD$mids -3.3683 0.4420 -7.62 2.53e-14 ***
---
Signif. codes: 0 ~Q***~R 0.001 ~Q**~R 0.01 ~Q*~R 0.05 ~Q.~R 0.1 ~Q ~R 1
(Dispersion parameter for poisson family taken to be 1)
Null deviance: 129.59 on 99 degrees of freedom
Residual deviance: 55.67 on 98 degrees of freedom
AIC: Inf
Number of Fisher Scoring iterations: 5
> summary(fmDPAll)
Deviance Residuals:
Min 1Q Median 3Q Max
-1.3687 -0.5672 -0.3386 0.0513 4.9680
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 1.478 0.148 9.988 <2e-16 ***
hDAll$mids -4.327 0.501 -8.636 <2e-16 ***
---
Signif. codes: 0 ~Q***~R 0.001 ~Q**~R 0.01 ~Q*~R 0.05 ~Q.~R 0.1 ~Q ~R 1
(Dispersion parameter for poisson family taken to be 1)
Null deviance: 181.646 on 99 degrees of freedom
Residual deviance: 74.457 on 98 degrees of freedom
AIC: Inf
Number of Fisher Scoring iterations: 5
As 55.67 < 74.457 the first model seems to fit better than the second
one, but how good is it? Should I compare these residual deviances
with chi-square? Should I look for some other model with smaller
residual deviance?
Best regards,
Alexander
More information about the R-help
mailing list