使用R survey包的svyby函数无法获取分组比例的置信区间

问题出在默认调用svyby时，仅提取了svyciprop返回的比例系数，没有捕获置信区间。可以通过自定义函数，让svyby同时输出比例、置信区间下限和上限。

修改后的代码示例 #

library(haven)
library(survey)

# 读取数据
Data_STEPs2016 <- read_dta("Downloads/New Data/Data_STEPs2016.dta")

# 构建抽样设计（注意语法优化：直接用变量名而非数据框$变量名）
steps1395 <- svydesign(id = ~code, weights = ~weight_Q, data = Data_STEPs2016)

# 自定义函数提取比例和置信区间，svyby会自动按分组整合结果
svyby(~t2, ~i07, design = steps1395, 
      FUN = function(x, design) {
        # 计算带置信区间的比例，可指定method（如logit/likelihood/asin）
        prop_result <- svyciprop(x, design = design, method = "logit")
        # 提取比例值、置信区间上下限并返回数据框
        data.frame(
          比例 = coef(prop_result),
          置信区间下限 = confint(prop_result)[1],
          置信区间上限 = confint(prop_result)[2]
        )
      })

关键说明 #

• 自定义函数中调用svyciprop时，需明确指定method参数（推荐logit，适合比例类数据的置信区间计算，避免边界值失真）。
• 通过coef()提取比例值，confint()提取置信区间的上下限，将这些结果整合成数据框返回，svyby就会按分组输出所有列。

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