如何解决R包EnhancedVolcano报"log2FoldChange is not numeric"错误

1. 解决log2FoldChange is not numeric!报错 #

从你提供的str(results)输出可以看到，baseMean、log2FoldChange、lfcSE、stat四列的类型都是字符型（chr），这是报错的核心原因。
这类问题通常是因为你用read.table加fill=T读入csv时，文件中存在字段数不匹配的行、数值列混入了非数字内容（比如空格、特殊符号、空值），导致整列被自动识别为字符型。你单独转as.numeric仍然报错的原因是转换过程中会生成NA值，没有提前清理。
修复后的完整代码如下：

library(EnhancedVolcano)
# 优先用read.csv读取csv文件，避免read.table的参数适配问题
results <- read.csv("Volcano_DEGlist_Age12_vs_Age0.csv", header = TRUE, stringsAsFactors = FALSE)
# 备份原始基因名，避免行名加后缀后丢失原始名称
results$raw_gene <- results[,1]
# 处理重复基因名生成唯一行名
row.names(results) <- make.unique(results[,1])
results[,1] <- NULL
# 批量转换所有数值列的类型，非数字内容会转为NA
num_cols <- c("baseMean", "log2FoldChange", "lfcSE", "stat", "pvalue", "padj")
results[num_cols] <- lapply(results[num_cols, drop=F], function(x) as.numeric(as.character(x)))
# 清理含NA的无效行
results <- na.omit(results)
# 绘图
pdf("Volcano_plot.pdf", width=10, height=8)
EnhancedVolcano(results,
                lab = rownames(results),
                x = "log2FoldChange",
                y = "padj")
dev.off()

如果不想图上的基因标签显示.1这类后缀，只需要将lab参数替换为lab = results$raw_gene即可，行名仅作为后台索引使用，不影响显示效果。

2. 非唯一基因名导入问题 #

R的data.frame底层规则要求行名必须唯一，EnhancedVolcano也要求输入数据的行名不能重复，因此无法直接导入非唯一基因名作为行名。除了make.unique之外，你也可以根据自己的需求自定义重复基因名的后缀规则，但必须保证行名唯一。
另外建议：如果是直接从DESeq2输出的差异结果，不需要导出csv再导入R，直接用results()函数生成的结果对象即可，既不会出现数据类型转换问题，也默认是唯一行名。

内容的提问来源于stack exchange，提问作者Rodriguez J Mathew