按共同列名筛选矩阵列表子集时返回空列表的问题

你的代码返回空列表的核心问题是矩阵子集筛选的语法错误：

• 矩阵用x[col]会默认按一维向量的方式索引（列优先），而非筛选列。如果col是列名，这种写法会试图匹配行名，自然找不到对应元素，最终返回空结果。

修复后的代码： #

mat.list <- list(mRNA=as.matrix(mrna.deg), Methylation=as.matrix(diff.meth), Protein=as.matrix(protein.dep))
# 获取所有矩阵的共享列名
common_cols <- Reduce(intersect, lapply(mat.list, colnames))
# 正确筛选列：用[, common_cols]指定列维度
mat.list <- lapply(mat.list, function(x) x[, common_cols])

额外检查点（若修复后仍为空）： #

• 先验证common_cols是否为空：运行print(common_cols)，如果输出是character(0)，说明三个矩阵本来就没有共享列名，需检查原始数据列名是否存在拼写差异、大小写不一致（比如GeneID和geneid）或完全无重叠的情况。
• 确认原始对象mrna.deg、diff.meth、protein.dep的列名格式一致，比如是否带有多余空格。

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