如何在R中基于已有相关矩阵文件绘制相关图（解决“矩阵不在[-1,1]范围”报错）

看起来你遇到的问题根源是输入的矩阵不是对称方阵，corrplot要求输入的是标准的相关矩阵（行数=列数、对称结构，且所有值在[-1,1]区间），而你的原始数据是一个6行3列的非完整矩阵，这才是触发报错的核心原因（虽然报错提示是值范围问题，但本质是矩阵结构不符合要求）。

下面是一步步的解决方案：

步骤1：正确读取数据 #

你的.txt文件里，第一列是基因名（行名），第一行是前3个基因的列名，剩下的数值是对应基因与列名基因的相关系数。首先要正确读取数据，指定第一列为行名：

# 如果文件是空格分隔，用默认sep即可；如果是制表符分隔，加上sep="\t"
m <- read.table("file.txt", header = TRUE, row.names = 1, check.names = FALSE)

读取后你会得到一个6行3列的数据框，行名是6个基因，列名是CXCL9、IL2RG、TAP1。

步骤2：构建完整的对称相关矩阵 #

我们需要把这个不完整的矩阵补全为6x6的对称方阵：

# 获取所有基因名称
all_genes <- rownames(m)

# 创建空的对称矩阵，行名和列名都用所有基因
full_corr_matrix <- matrix(NA, nrow = length(all_genes), ncol = length(all_genes),
                           dimnames = list(all_genes, all_genes))

# 填充对角线（每个基因与自身的相关系数为1）
diag(full_corr_matrix) <- 1

# 填充已知的相关系数
full_corr_matrix[rownames(m), colnames(m)] <- as.matrix(m)

# 对称填充矩阵（相关矩阵是对称的，A与B的相关系数等于B与A的）
full_corr_matrix[lower.tri(full_corr_matrix)] <- t(full_corr_matrix)[lower.tri(full_corr_matrix)]

现在你得到的full_corr_matrix就是一个标准的对称相关矩阵，所有值都在0.5~1之间，完全符合corrplot的要求。

步骤3：用corrplot绘制相关图 #

现在就可以正常使用corrplot绘制你想要的相关图了：

library(corrplot)
corrplot(full_corr_matrix, type = "lower", method = "number")

这个命令会绘制下三角区域的相关图，同时显示具体的相关系数数值。

备选方案：用热图展示 #

如果你想用热图来展示相关关系，推荐使用pheatmap包，可视化效果更灵活：

library(pheatmap)
# 关闭聚类（因为我们的矩阵是已知的相关关系，不需要聚类）
pheatmap(full_corr_matrix, display_numbers = TRUE, 
         cluster_rows = FALSE, cluster_cols = FALSE)

这个热图会用颜色深浅表示相关系数的大小，同时在每个单元格显示具体数值。

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