R: svyglm(调查包)分类表

Question

我试着用'survey'包做逻辑回归分析。

我正在努力与我的模型的分类表。

design1 <- svydesign(ids = ~num, weights = ~weights, data = mydata, fpc=~nm, strata = ~strata1)
glm_m<- svyglm(var1 ~ sex + age1 + educat, design = design1, family = quasibinomial())

接下来我应该如何继续做分类表?其次，在此之后，用调查数据计算Hosmer-Lemeshow拟合优度(GOF)检验?

Answer 1

您需要一个分类来获得分类表。假设您希望阈值为0.5。然后svytable得到一个人口表

design1<-update(design1, fittedp=fitted(design1))
design1<-update(design1, fittedclass= fittedp>0.5)
svytable(~var1+fittedclass, design1)

Hosmer-Lemeshow测试至少有两种不同的扩展来调查数据，但我不知道有R实现。