[R] Generating a model fitness when score using svyglm?

Thomas Lumley tlumley at u.washington.edu
Wed Apr 28 21:11:45 CEST 2010

On Wed, 28 Apr 2010, Brad Fulton wrote:

> Does anyone know how to calculated a BIC score (or an equivalent model
> fitness score) when using svyglm for logistic regressions?

No.  That is, the model is not fitted by maximum likelihood, so BIC doesn't approximate posterior probabilities.

Now, the deviance returned by svyglm() is scaled to the sample size, so if the survey design isn't informative it should be more or less in the same ballpark as a deviance from an independent sample, and the usual BIC calculation might give somewhat helpful results.


Thomas Lumley			Assoc. Professor, Biostatistics
tlumley at u.washington.edu	University of Washington, Seattle

More information about the R-help mailing list