R生成贫困状态指标时报「替换项数非长度倍数」错误如何解决

• 第一处语法错误：你在定义poverty向量长度时，误将赋值符号<-写为小于号<，导致poverty没有被正确设置为和数据集行数一致的长度，后续赋值时两边长度不匹配，直接触发你遇到的「number of items to replace is not a multiple of length」报错。
• 第二处逻辑错误：给poverty的子集赋值时，右侧使用了完整的data$inc向量，左侧是满足筛选条件的部分观测，两者长度不一致，也会引发长度不匹配问题。
• 额外问题：原代码逐行赋值的写法冗余度极高，后续调整参数时维护成本很高，可通过映射表批量匹配大幅简化。

# 1. 家庭收入编码映射为实际收入值
inc_map <- c(
  "1" = NA, "2" = 500, "3" = 2000, "4" = 4000, "5" = 6250,
  "6" = 8750, "7" = 12500, "8" = 17500, "9" = 22500, "10" = 30000,
  "11" = 42500, "12" = 62500, "13" = 87500, "14" = 150000, "15" = 200000,
  "98" = NA, "99" = NA
)
data$inc <- as.numeric(inc_map[data$BYFAMINC])
table(data$inc)

# 2. 家庭规模清洗
data$fmlysz <- as.numeric(data$BYFAMSIZ)
data$fmlysz[data$fmlysz %in% c(98,99)] <- NA
table(data$fmlysz)

# 3. 不同家庭规模对应贫困线映射，直接计算收入/贫困线比值
pov_threshold <- c(
  "1" = NA, "2" = 7704, "3" = 9435, "4" = 12092, "5" = 14304,
  "6" = 16146, "7" = 18232, "8" = 20253, "9" = 24129, "10" = 24129,
  "98" = NA, "99" = NA
)
data$poverty <- data$inc / pov_threshold[as.character(data$fmlysz)]
data$poverty <- as.numeric(data$poverty)

# 如需生成二分类贫困标识（收入低于贫困线为1，否则为0）可补充以下代码
# data$pov_flag <- ifelse(data$poverty < 1, 1, 0)

table(data$poverty)
summary(data$poverty)

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