R help - Apr 2010 - Generating a model fitness when score using svyglm?

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

Thanks

Brad
-- 
View this message in context:
http://r.789695.n4.nabble.com/Generating-a-model-fitness-when-score-using-svyglm-tp2069280p2069280.html
Sent from the R help mailing list archive at Nabble.com.

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

So are you saying that one way to estimate goodness of fit would be to run
each of models using glm() and compare their BIC scores?

Is there a recommended way to demonstrate improvements in goodness of fit
when using svyglm?

Thanks

Brad
-- 
View this message in context:
http://r.789695.n4.nabble.com/Generating-a-model-fitness-when-score-using-svyglm-tp2069280p2073835.html
Sent from the R help mailing list archive at Nabble.com.