DESeq2调用DESeqDataSetFromMatrix时报计数矩阵为逻辑型错误如何解决

Error in DESeqDataSet(se, design = design, ignoreRank) : counts matrix should be numeric, currently it has mode: logical

核心诱因 #

报错的本质是你传入DESeqDataSetFromMatrix的countData参数df经过代码处理后已经不是数值型的计数矩阵，变成了逻辑型，由以下3个代码逻辑错误导致：

逐项修改建议 #

• filter条件语法错误 #

你在筛选coldata时给condition加了双引号，R会将其识别为纯字符串而非数据框列，字符串"condition"和"Tumor"/"Normal"/"At_risk"永远不相等，筛选后coldata为空表，后续用空行名选df列会返回0列/全NA的逻辑型结果。
修改后代码：

coldata <- coldata %>%
  dplyr::filter(condition == "Tumor"| condition == "Normal" | condition == "At_risk")
# 筛选后先验证有有效样本
stopifnot(nrow(coldata) > 0)

• 计数表未设置基因名为行名 #

read.csv读入计数表后没有把第一列的基因ID设为行名，df的第一列是字符型基因名，会导致后续数据类型异常；原有代码geneID <- rownames(df)只能拿到1、2、3的数字索引，无法获取真实基因ID。
修改后读入代码：

# 读入时指定第一列为行名，确保剩下的列都是数值型计数
df <- read.csv(df, row.names = 1)

• 样本信息表未设置样本名为行名+列筛选前未验证匹配关系 #

原有代码读入coldata时未指定样本ID为行名，rownames(coldata)默认生成1、2、3的数字索引，和df的样本列名完全不匹配，dplyr::select会返回全NA列，全NA的数值列会被R自动识别为逻辑型。
修改后增加验证步骤：

# 读入coldata时指定第一列样本ID为行名
coldata <- read.csv(coldata_path, row.names = 1)
# 确认coldata行名全部存在于df列名中
stopifnot(all(rownames(coldata) %in% colnames(df)))
# 筛选列
df <- df[, rownames(coldata)]
# 强制确认df为数值型整数矩阵
df <- as.matrix(df)
storage.mode(df) <- "integer"

修正后完整可运行代码 #

# 读入计数表
df_path <- "C:\\Users\\reach\\OneDrive\\Desktop\\REU Summer 2021\\Autotaxin Cancer\\mRNA_counts_and_annotation\\mRNA_counts.csv"
df <- read.csv(df_path, row.names = 1)

# 读入样本信息表，指定第一列为样本ID行名
coldata_path <- "C:\\Users\\reach\\OneDrive\\Desktop\\REU Summer 2021\\Autotaxin Cancer\\mRNA_counts_and_annotation\\coldata.txt"
coldata <- read.csv(coldata_path, row.names = 1)

# 筛选样本
coldata <- coldata %>%
  dplyr::filter(condition == "Tumor"| condition == "Normal" | condition == "At_risk")
stopifnot(nrow(coldata) > 0)
coldata$condition <- as.factor(coldata$condition)

# 标准化列名，替换横杠为下划线
colnames(df) <- stringr::str_replace_all(colnames(df), "-", "_")

# 匹配计数表和样本信息表
stopifnot(all(rownames(coldata) %in% colnames(df)))
df <- df[, rownames(coldata)]
stopifnot(all(rownames(coldata) == colnames(df)))

# 强制转换为数值型整数矩阵，避免类型异常
df <- as.matrix(df)
storage.mode(df) <- "integer"

# 构建dds对象
ddsall <- DESeq2::DESeqDataSetFromMatrix(countData = df,
                                         colData = coldata,
                                         design = ~ condition)

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