如何在R中使用Survey包运行线性回归

Question

我正在尝试分析国民健康采访调查的数据，以研究移民与薪资之间的关系。为此，我需要使用R中的调查包来适当考虑加权、分层等因素... 我已经使用调查函数编写了代码，但似乎无法运行回归，因为我收到以下错误消息:

"Error in onestrat(x[index, , drop = FALSE], clusters[index], nPSU[index][1], : Stratum (6001) has only one PSU at stage 1"

我不确定该如何修复，并希望有人能帮忙！这是我迄今为止运行的代码。

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))

如果有需要，以下是疾病控制与预防中心提供的更多关于如何使用此功能的文档：https://www.cdc.gov/nchs/data/nhis/2016var.pdf

谢谢！

- juliah0494

Answer 1

针对该错误，您可能希望使用此处列出的最保守选项

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

如需更多详细信息，请查看http://asdfree.com/national-health-interview-survey-nhis.html中的回归模型和关联测试部分。

- Anthony Damico