如何使用 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
的 回归模型和关联测试 部分
我正在尝试分析来自全国健康访谈调查的数据,以了解出生与薪水之间的关联。为此,我需要使用 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
的 回归模型和关联测试 部分