R绘制物种丰度热图报错：‘x’需为数值矩阵问题求助

这个错误看起来特别矛盾——明明提示要求数值矩阵，但又说你的输入就是数值矩阵，其实问题出在R对"数值矩阵"的严格定义和你输入对象的实际类型不匹配上，大概率是你的log_mean_wide是数据框（data frame）而非真正的矩阵类型，或者数据框里藏着非数值元素。下面一步步帮你排查解决：

1. 先确认输入对象的真实类型 #

base R的heatmap()函数只认matrix类型，不认数据框，哪怕数据框里全是数值也不行。先跑这两行代码看看：

class(log_mean_wide)
str(log_mean_wide)

如果输出里有data.frame，那这就是问题根源了。

2. 把数据框转换为数值矩阵 #

根据你的数据结构（样地×物种），分两种情况处理：

• 如果你的数据框已经把样地设置为行名（没有单独的样地列）：

  log_mean_matrix <- as.matrix(log_mean_wide)
  

• 如果数据框第一列是样地名（作为普通列存在）：

  # 先把样地列设为行名
  rownames(log_mean_wide) <- log_mean_wide[, 1]
  # 再把剩余的物种丰度列转成矩阵
  log_mean_matrix <- as.matrix(log_mean_wide[, -1])
  

转换后再验证一下：

is.numeric(log_mean_matrix)

返回TRUE就说明没问题了。

3. 排查是否存在隐藏的非数值元素 #

有时候数据框里的列看起来是数值，但实际是因子类型（比如导入数据时自动识别错了），可以用这行检查：

any(sapply(log_mean_wide, function(x) !is.numeric(x)))

如果返回TRUE，先把这些列转成数值：

log_mean_wide[] <- lapply(log_mean_wide, function(x) {
  if (is.factor(x)) as.numeric(as.character(x)) else x
})

再重复步骤2转成矩阵。

4. 重新运行热图代码 #

用转换后的矩阵替换原来的数据框：

lrc <- rainbow(nrow(log_mean_matrix), start = 0, end = .3)
lcc <- rainbow(ncol(log_mean_matrix), start = 0, end = .3)
logmap <- heatmap(log_mean_matrix, col = cm.colors(256), scale = "column")

额外小建议 #

如果经常画热图，推荐用pheatmap包，它支持直接输入数据框，还自带更多美化和功能：

install.packages("pheatmap")
library(pheatmap)
pheatmap(log_mean_wide, color = cm.colors(256), scale = "column")

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