R语言使用if else循环为dataframe新增评分列赋值的问题排查

你写的代码存在3个明确错误：

• 循环逻辑完全错误：for(oP.Res in Merged)是按列遍历数据框，你的数据集有22列，这段代码只会循环22次，根本不是逐行判断的逻辑
• 没有赋值动作：所有条件分支里只写了"0""1"这类字符串，既没有指定存储位置，也没有给目标列赋值，运行后不会对数据集产生任何修改
• 判断逻辑不符合R的语法规则：if()函数要求传入单个逻辑值，你直接在判断条件里调用Merged$MEAS_ANAL_METH_CODE是取整列所有值做判断，会触发长度不匹配警告，运行结果完全不可控

R处理这类按条件分箱赋值的场景不需要写显式for循环，用dplyr包的mutate+case_when写法可读性最高、运行速度最快，2000行数据可以瞬间完成计算，逻辑完全匹配你设定的评分规则：

# 首次使用先运行安装命令：install.packages("dplyr")
library(dplyr)

Merged <- Merged %>%
  mutate(
    # 计算OP_METH列评分
    OP_METH = case_when(
      MEAS_ANAL_METH_CODE == 24 & oP.Res < 0.01 ~ 0,
      MEAS_ANAL_METH_CODE == 24 & oP.Res >= 0.01 & oP.Res < 0.0375 ~ 1,
      MEAS_ANAL_METH_CODE == 24 & oP.Res >= 0.0375 & oP.Res < 0.0725 ~ 2,
      MEAS_ANAL_METH_CODE == 24 & oP.Res >= 0.0725 & oP.Res < 0.1 ~ 3,
      MEAS_ANAL_METH_CODE == 24 & oP.Res >= 0.1 ~ 4,
      MEAS_ANAL_METH_CODE == 21 & oP.Res < 0.01 ~ 4,
      MEAS_ANAL_METH_CODE == 21 & oP.Res >= 0.01 & oP.Res < 0.0375 ~ 3,
      MEAS_ANAL_METH_CODE == 21 & oP.Res >= 0.0375 & oP.Res < 0.0725 ~ 2,
      MEAS_ANAL_METH_CODE == 21 & oP.Res >= 0.0725 & oP.Res < 0.1 ~ 1,
      # 方法码非21/24、值缺失等其余情况统一赋值0
      TRUE ~ 0
    ),
    # 计算TP_METH列评分，仅需把oP.Res替换为TP.Res，规则完全一致
    TP_METH = case_when(
      MEAS_ANAL_METH_CODE == 24 & TP.Res < 0.01 ~ 0,
      MEAS_ANAL_METH_CODE == 24 & TP.Res >= 0.01 & TP.Res < 0.0375 ~ 1,
      MEAS_ANAL_METH_CODE == 24 & TP.Res >= 0.0375 & TP.Res < 0.0725 ~ 2,
      MEAS_ANAL_METH_CODE == 24 & TP.Res >= 0.0725 & TP.Res < 0.1 ~ 3,
      MEAS_ANAL_METH_CODE == 24 & TP.Res >= 0.1 ~ 4,
      MEAS_ANAL_METH_CODE == 21 & TP.Res < 0.01 ~ 4,
      MEAS_ANAL_METH_CODE == 21 & TP.Res >= 0.01 & TP.Res < 0.0375 ~ 3,
      MEAS_ANAL_METH_CODE == 21 & TP.Res >= 0.0375 & TP.Res < 0.0725 ~ 2,
      MEAS_ANAL_METH_CODE == 21 & TP.Res >= 0.0725 & TP.Res < 0.1 ~ 1,
      TRUE ~ 0
    )
  )

注意：如果后续需要对评分做数值计算，代码里直接写0-4的数字即可，不要加引号；如果需要存为文本分类，再给数字加双引号转为字符串类型。

如果一定要用原生for循环实现，需要按行索引遍历，逐行取值、逐行给目标列对应位置赋值，参考代码如下：

# 先初始化两列，填充默认值0
Merged$OP_METH <- 0
Merged$TP_METH <- 0

# 按行号遍历，从第1行到最后1行
for(i in seq_len(nrow(Merged))){
  # 提取当前行的所需字段值
  current_meth <- Merged$MEAS_ANAL_METH_CODE[i]
  current_op <- Merged$oP.Res[i]
  current_tp <- Merged$TP.Res[i]
  
  # 给OP_METH赋值
  if(current_meth == 24){
    if(current_op < 0.01) Merged$OP_METH[i] <- 0
    else if(current_op < 0.0375) Merged$OP_METH[i] <- 1
    else if(current_op < 0.0725) Merged$OP_METH[i] <- 2
    else if(current_op < 0.1) Merged$OP_METH[i] <- 3
    else Merged$OP_METH[i] <- 4
  }else if(current_meth == 21){
    if(current_op < 0.01) Merged$OP_METH[i] <- 4
    else if(current_op < 0.0375) Merged$OP_METH[i] <- 3
    else if(current_op < 0.0725) Merged$OP_METH[i] <- 2
    else if(current_op < 0.1) Merged$OP_METH[i] <- 1
    else Merged$OP_METH[i] <- 0
  }
  
  # 给TP_METH赋值
  if(current_meth == 24){
    if(current_tp < 0.01) Merged$TP_METH[i] <- 0
    else if(current_tp < 0.0375) Merged$TP_METH[i] <- 1
    else if(current_tp < 0.0725) Merged$TP_METH[i] <- 2
    else if(current_tp < 0.1) Merged$TP_METH[i] <- 3
    else Merged$TP_METH[i] <- 4
  }else if(current_meth == 21){
    if(current_tp < 0.01) Merged$TP_METH[i] <- 4
    else if(current_tp < 0.0375) Merged$TP_METH[i] <- 3
    else if(current_tp < 0.0725) Merged$TP_METH[i] <- 2
    else if(current_tp < 0.1) Merged$TP_METH[i] <- 1
    else Merged$TP_METH[i] <- 0
  }
}

内容的提问来源于stack exchange，提问作者Dani.Elizabeth