将Stata的INCA3调查设计命令转为R的svydesign时遇权重缺失错误

背景 #

我正在处理INCA3营养数据集，分析前需完成调查设计声明，用户文档第69页提供了对应的Stata抽样设计命令：

svyset zae [pweight=pond_indiv_adu_pop3], ///
    strata(strate) fpc(fpc1) vce(linearized) singleunit(scaled) ///
    || NOMEN, fpc(fpc2) || NOIND, fpc(fpc3)

R代码转换尝试 #

由于不使用Stata且缺乏此类调查数据的处理经验，我尝试用R的survey包中的svydesign命令转换上述Stata代码，实现如下：

svydesign(id=~zae+NOMEN+NOIND,
          strata=~strate,
          weights=~pond_indiv_adu_pop3,
          fpc=~fpc1+fpc2+fpc3,
          data=conso_alim)

报错信息 #

运行上述R代码时出现错误：

Error in (function (object, ...)  : missing values in `weights'

请问我的转换实现是否存在问题？

内容的提问来源于stack exchange，提问作者DSMLSC