Hello,

I was looking for a way to evaluate the goodness-of-fit of a logistic regression model. After googling, I found that I could use "resid(fit, 'gof')" method implemented in the rms package. However, since I am not used to the "le Cessie-van Houwelingen normal test statistic," I do not know which statistic from the returned from the "resid(fit, 'gof')" call that I could use to evaluate the goodness of fit.

When I ran the "resid(fit, 'gof')", I got the following results:
##############################################
Sum of squared errors   Expected value|H0             SD
      6844.684594       6805.672315         2.790969
            Z              P
       13.978043         0.000000
##############################################

I tried to read the le Cessie and van Houwelingen's original paper, but I found that it required prerequisite knowledge I don't current have.
Could someone explain how to interpret the results from "resid(fit, 'gof') call?

Any help would be much appreciated.

Young-Jin Lee
______________________________________________
R-help@(protected)
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.