R 如何报告svyglm对象的拟合优度?
一个关于R调查包的问题。我知道这是非常基本的,但我什么也没发现。因此:R 如何报告svyglm对象的拟合优度?,r,survey,R,Survey,一个关于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 + mobilit
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
图书馆(调查)
数据(api)
dclus2 | t |)
(截距)819.279121.390038.3024.65e-13***
厄尔-0.5167 0.3240-1.595 0.139
膳食-3.12320.2781-11.2312.29e-07***
流动性-0.1689 0.4449-0.380 0.711
---
签名。代码:0'***'0.001'***'0.01'*'0.05'.'0.1''1
(高斯族的色散参数取3157.85)
Fisher评分迭代次数:2
有没有像R平方或任何其他统计数据那样的东西能让读者对模型与数据的拟合程度有一个好的印象
谢谢。2016年末,我发现了一个名为“poliscidata”的软件包,它有一个函数fit.svyglm,返回R平方和Adj。svyglm对象的R平方输出。我成功地使用了它。如果你不使用
摘要
你会得到更多的统计信息。。也许可以浏览一下你的svyglm
电话的属性:)谢谢@AnthonyDamico。我根据你的建议,通过取R2=1-object$deviance/object$null.deviance