R语言ifelse判断小于极小Bonferroni校正p值时输出全yes异常排查

该问题的核心触发原因是df$score列被锁定为字符类型，后续的数值大小比较被R隐式转换为字符串字典序比较规则，而非原本预期的数值大小比较。

• 代码中在循环计算p值前，对score列做了空字符串初始化（原代码中预分配列空间的步骤），直接将该列的存储类型定义为字符型。后续循环中即使往该列赋值计算得到的数值，R也会自动将数值转换为字符串存入，不会自动还原为数值类型。
• R的比较运算存在强制类型转换规则：当比较操作的两个值中任意一个为字符类型时，系统会将两边的值全部转换为字符串，按照ASCII码的字典顺序逐位比较大小，而非按照数值逻辑判断。
• 结合运行场景，具体的错误逻辑可以完全对应复现：
  • Bonferroni阈值6.52912e-06转换为字符串为"6.52912e-06"，首个字符为ASCII码54对应的'6'
  • 实际存储的score值对应的字符串分别为"1"（首字符为ASCII49的'1'）、"0.9644..."（首字符为ASCII48的'0'），两类字符串的首字符ASCII值均小于'6'
  • 按照字典序规则，所有score值都会被判定为“小于”Bonferroni阈值，因此所有样本被统一标记为离群点（outler = yes）
• 提到的“替换为0.05的未校正p值后判断流程可正常运行”属于偶然出现的假象：0.05转换为字符串为"0.05"，首字符为ASCII48的'0'，此时字典序的比较结果刚好和数值大小的判断逻辑重合，造成了流程正常的错觉，实际上比较逻辑的类型错误始终存在。

• 修正score列的存储类型：直接将整列强制转换为数值型，从根源上消除类型不匹配问题，对应代码为：df$score <- as.numeric(df$score)
• 类型修正完成后，重新执行离群点判断逻辑，此时比较双方均为数值类型，会按照正常的数值大小规则运算。
• 流程优化建议：预分配data frame列空间时，避免使用空字符串初始化数值计算列的写法；如果需要提前预留列位置，可直接初始化为数值型空值，例如df$score <- numeric(nrow(df))，从设计上规避类型锁定问题。

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