[R] [Q] Goodness-of-fit test of a logistic regression model using rms package

GMail (KU) youngjin.lee.ku at gmail.com
Wed Sep 1 17:57:56 CEST 2010


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


More information about the R-help mailing list