如何使用 R 中的 Survey Package 进行 运行 线性回归

How to use Survey Package in R to run linear regression

我正在尝试分析来自全国健康访谈调查的数据,以了解出生与薪水之间的关联。为此,我需要使用 R 中的 survey 包来正确考虑权重、分层等......我已经用 survey 函数编写了代码,但似乎无法回归到 运行 因为我收到错误消息:

“onestrat 出错(x[index, , drop = FALSE], clusters[index], nPSU[index][1], : Stratum (6001) 在第 1 阶段只有一个 PSU"

我不确定要解决什么问题,希望有人能提供帮助!这是我 运行 到目前为止的代码

nhissvy <- svydesign(id=~PSU, strata=~STRATA,
                 nest = TRUE,
                 weights=~PERWEIGHT,
                 data=NHIS1)

 summary(svyglm(adj_SALARY ~ AGE+IMMIGRANT+FEMALE+MARRIED+NEWEDUC+YEAR_2015+YEAR_2016+YEAR_2017+YEAR_2018, design = nhissvy))

以下是来自 CDC 的更多文档,说明如何使用此功能(如果有帮助的话):https://www.cdc.gov/nchs/data/nhis/2016var.pdf

谢谢!

对于该错误,您可能需要使用 the most conservative option listed here

options( survey.lonely.psu = "adjust" )

有关更多详细信息,请查看 http://asdfree.com/national-health-interview-survey-nhis.html

回归模型和关联测试 部分