如何使用tbl_svysummary汇总5级Likert量表有序数据的中位数？

步骤1：将Likert变量转换为有序因子 #

首先必须确保5个Likert问题变量是有序因子，这样R才能识别选项的顺序逻辑，正确计算中位数。假设你的变量名为q1到q5，代码如下：

# 定义Likert选项的固定顺序
likert_levels <- c("very bad", "bad", "moderate", "good", "very good")

# 将目标变量批量转换为有序因子
your_survey_data <- your_survey_data %>%
  mutate(across(q1:q5, ~factor(., levels = likert_levels, ordered = TRUE)))

步骤2：在tbl_svysummary中指定中位数统计量 #

默认情况下，tbl_svysummary对分类变量仅展示频数和占比，你需要通过statistic参数自定义输出中位数。注意先创建正确的调查设计对象：

# 创建复杂抽样调查设计对象（替换weight_var为你的权重变量）
survey_design <- svydesign(id = ~1, weights = ~weight_var, data = your_survey_data)

# 生成包含中位数的汇总表
tbl_svysummary(
  survey_design,
  include = q1:q5, # 指定要汇总的5个Likert问题
  statistic = all_categorical() ~ "{median}", # 对有序分类变量输出中位数
  label = q1 ~ "问题1", q2 ~ "问题2", q3 ~ "问题3", q4 ~ "问题4", q5 ~ "问题5" # 可选：给变量添加友好标签
)

补充技巧 #

• 若需同时展示中位数与频数/占比，可修改statistic参数为：all_categorical() ~ "{median}\n(n = {n} [{p}%])"
• 务必确认调查设计对象的权重、集群等抽样参数设置正确，否则中位数计算会偏离实际抽样特征
• 有序因子的中位数计算基于其内部整数编码（"very bad"=1，"bad"=2……），最终会映射回对应的选项标签展示

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