使用pivot_wider后空白转0失败，无法运行DGEList函数

一、从长格式直接转宽并自动补0（推荐方案） #

别先转宽再手动补值，一步到位更省心：

1. 先清理原始长数据，确保表达值是数值型，删掉无效的基因行（比如你宽表里的null行）
2. 用tidyr::pivot_wider转宽时直接指定空白填充为0

library(tidyr)
# 假设原始长数据框叫long_df
long_df$VALUE <- as.numeric(long_df$VALUE)
long_df <- long_df[long_df$Gene_Symbol != "null", ]

# 转宽并自动填充空白为0
wide_df <- pivot_wider(long_df, 
                       id_cols = Gene_Symbol, 
                       names_from = Sample_ID, 
                       values_from = VALUE, 
                       values_fill = 0)

# 把基因名设为行名，去掉多余列
rownames(wide_df) <- wide_df$Gene_Symbol
wide_df <- wide_df[, -match("Gene_Symbol", colnames(wide_df))]

二、修复已生成的错误宽表 #

如果已经有了带空白的宽表，别用apply处理——它会把数据转成字符型，导致全变0。正确操作：

# 把空白和NA统一替换成0
wide_df[wide_df == "" | is.na(wide_df)] <- 0

# 确保所有列都是数值型
wide_df <- as.data.frame(lapply(wide_df, as.numeric))

# 删掉无效的null行
wide_df <- wide_df[rownames(wide_df) != "null", ]

三、生成DGEList #

处理好的宽表直接用edgeR包生成：

library(edgeR)
dge <- DGEList(counts = wide_df)

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