R语言基于行值与列名替换列值及列子集越界报错解决

• 运算符优先级错误：R中冒号:的运算优先级高于加号+，代码中polydata$s_mean[x]:187+10的逻辑是先生成从s_mean[x]到187的整数序列，再给序列每个值加10，而非先计算187+10作为序列终点。比如某行s_mean值为294时，这段代码会生成304、303……197的列索引，而数据集总共只有197列，自然会提示找不到304及以上的列。
• 列偏移计算错误：原代码备注doy列从第11列开始，但从样例数据看，前16列分别是ID、doy、FireID、Year、sE、NAME、L1NAME、ID_2、area、s_count、s_mean、s_median、s_stdev、s_min、doydiff、ID_E5，数字命名的doy列实际从第17列（列名32）开始，偏移量10的计算完全错误。
• 实现思路偏差：原代码逐行切片选列的逻辑，只能筛选保留符合条件的列，无法实现“指定单元格填NA”的需求，和预期最终效果不匹配。

直接用列名匹配逻辑替换逐行切片写法，无需手动计算列位置，自动识别目标列，从根源避免索引越界问题：

library(readr)
library(dplyr)

# 读入并预处理数据
polydata <- read_csv("path/E15.csv")
polydata$s_mean <- round(polydata$s_mean)

# 自动识别所有数字命名的doy列（列名为32到212的数值列）
doy_cols <- polydata |> 
  select(where(~ is.numeric(.) && as.numeric(cur_column()) %in% 32:212)) |> 
  names()
doy_col_values <- as.numeric(doy_cols)

# 遍历行，将列名小于对应行s_mean的单元格赋值为NA
for (row_idx in 1:nrow(polydata)) {
  row_smean <- polydata$s_mean[row_idx]
  # 找出当前行需要设为NA的列
  target_na_cols <- doy_cols[doy_col_values < row_smean]
  polydata[row_idx, target_na_cols] <- NA
}

# 提取最终需要的输出列
final_output <- polydata |> 
  select(ID, s_mean, all_of(doy_cols))

代码说明：

• 无需手动计算列偏移，自动识别所有doy列，避免列索引越界问题
• 赋值逻辑完全匹配需求：仅将列名数值小于该行s_mean的单元格设为NA，其余值保留
• 最终输出结构为ID、s_mean + 所有32~212命名的doy列，和示例结构完全一致，运行效率远高于逐行切片拼接列表的写法。

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