R语言获取104×182数据集各列最大值行索引的报错解决

核心错误分析 #

• 循环内错误使用了Z0[, p]：p是1到182的整列范围，这会导致每次循环都取所有列，而非当前要处理的第i列，应改为Z0[, i]。
• 数据类型问题：报错提示'list' object cannot be coerced to type 'double'，说明Z0可能是tibble或list类型，直接用[,i]取列会返回子数据框而非向量，需改用[[i]]提取向量，或者将Z0转为矩阵。

修正后的循环代码 #

library(tibble)

p <- 1:182
DF <- tibble(i = numeric(), max = numeric())
for (i in p) {
  # 用[[i]]提取单列向量，避免list类型报错
  max_idx <- which.max(Z0[[i]])
  DF <- DF %>% add_row(i = i, max = max_idx)
}

DF %>% filter(max != 1) %>% distinct(max, .keep_all = TRUE)

更高效的非循环写法（推荐） #

没必要用循环，用apply函数可直接按列处理，代码更简洁且效率更高：

# 按列提取每列最大值的行索引
max_indices <- apply(Z0, 2, which.max)

# 转换为目标tibble格式
DF <- tibble(i = 1:182, max = max_indices)

# 执行后续过滤与去重
DF %>% filter(max != 1) %>% distinct(max, .keep_all = TRUE)

补充提示 #

• 如果Z0包含非数值列，先过滤处理：Z0_numeric <- Z0 %>% select(where(is.numeric))，再执行上述代码。
• apply的第二个参数2代表按列操作，若需按行处理则改为1。

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