r - 如何报告 svyglm 对象的拟合优度？

Question

关于 R 中调查包的问题。我知道这是非常基本的，但我什么也没找到。所以：

library(survey)
data(api)
dclus2 <- svydesign(id=~dnum, weights=~pw, data=apiclus1, fpc=~fpc)
summary(svyglm(api00 ~ ell + meals + mobility, design = dclus2))
Call:
svyglm(formula = api00 ~ ell + meals + mobility, design = dclus2)

Survey design:
svydesign(id = ~dnum, weights = ~pw, data = apiclus1, fpc = ~fpc)

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept) 819.2791    21.3900  38.302 4.65e-13 ***
ell          -0.5167     0.3240  -1.595    0.139    
meals        -3.1232     0.2781 -11.231 2.29e-07 ***
mobility     -0.1689     0.4449  -0.380    0.711    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for gaussian family taken to be 3157.85)

Number of Fisher Scoring iterations: 2

是否有任何类似于 R 平方或任何其他统计数据可以让读者了解模型对数据的拟合程度？

谢谢。

Answer 1

我发现 2016 年末推出了一个名为“poliscidata”的包，它有一个函数 fit.svyglm，它返回 R 平方和 Adj。svyglm 对象的 R Square 输出。我已经成功地使用了它。