我最近在分析NHANES数据时遇到了点疑问,想请大家帮忙解答一下,谢谢。
正常NHANES数据进行加权分析时,是用下面的步骤:
library(survey)
Design <- svydesign(id = ~psu,
strata = ~strata,
weights = ~persWeight,
nest = TRUE,
data = Analysis)
但是由于我清洗的数据包括1999年至2018年的数据,因此权重分析时有4年权重(WTMEC4YR)和2年权重(WTMEC2YR)
Design <- svydesign(id = ~SDMVPSU,
strata = ~SDMVSTRA,
weights = ~(1/5*WTMEC4YR+4/5*WTMEC2YR),
nest = TRUE,
data = Analysis)
这样计算会报错,请教一下大家该怎么处理,谢谢。