使用svydesign时遇Stratum(0)仅含一个PSU错误求助

错误原因 #

报错的核心问题是：数据中编号为0的分层（stratum）仅包含1个抽样单位（PSU）。在复杂抽样的方差计算逻辑里，每个分层至少需要2个PSU才能有效估计抽样误差，单个PSU无法提供方差计算的基础，因此触发错误。

可行解决方案 #

• 合并单个PSU的分层
  先统计各分层的PSU数量，再将只有1个PSU的分层与业务逻辑相似的分层合并（比如同国家、同人口特征组的分层）：

  # 统计每个stratum的PSU数量
  psu_count <- steps_3jan_matrix %>% 
    group_by(stratum) %>% 
    summarise(n_psu = n_distinct(psu))
  
  # 筛选出仅含1个PSU的分层
  single_psu_strata <- psu_count$stratum[psu_count$n_psu == 1]
  

  后续可根据数据特征（如国家、年龄组）完成分层合并，再重新构建调查设计。

• 移除单个PSU的分层
  若单PSU分层的样本占比极低，对整体分析结果影响可忽略，直接剔除这些分层：

  # 过滤掉仅含1个PSU的分层
  steps_filtered <- steps_3jan_matrix %>%
    group_by(stratum) %>%
    filter(n_distinct(psu) > 1) %>%
    ungroup()
  
  # 重新创建调查设计并计算均值
  steps_design <- svydesign(id = ~psu, 
                            strata = ~stratum,
                            nest = TRUE,
                            weights = ~wstep1,
                            data = steps_filtered)
  
  svymean(~age, steps_design, na.rm = TRUE)
  

• 调整svydesign的参数保留分层
  若必须保留这些单PSU分层，可在创建调查设计时设置lonely.psu参数调整方差计算逻辑：

  steps_3jan_design <- svydesign(id = ~psu, 
                                 strata = ~stratum,
                                 nest = TRUE,
                                 weights = ~wstep1,
                                 data = steps_3jan_matrix,
                                 lonely.psu = "adjust")
  
  svymean(~age, steps_3jan_design, na.rm = TRUE)
  

  lonely.psu="adjust"会通过合并方差或使用全局方差估计的方式，绕过单个PSU无法计算方差的问题，具体调整逻辑由survey包自动处理。

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