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

[Thomas Lumley]

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

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