如何将geom_tile热图图例标签改为log10格式并解决NA报错

你的代码存在两个核心问题：

• 生成的Z值包含负数，直接对负数取log10()会产生NA，这是报错的根源之一；
• 在aes()里提前对Z做了log10转换，后续在scale_fill_gradient中使用log10相关的格式转换会逻辑混乱，导致NA生成。

下面提供两种可行的解决方案：

方法1：使用原始Z值，通过scale参数实现log10转换 #

这种方法更符合ggplot的语法逻辑，让scale自动处理数值转换，同时避免手动转换带来的NA问题：

library(tidyverse)

# 创建示例数据集
x <- LETTERS[1:20]
y <- paste0("var", seq(1,20))
data <- expand.grid(X=x, Y=y)
length_measurement <- 400
data$Z <- rnorm(400, 0.000025, 0.0005)

# 处理Z的负数和0：log10不接受非正数，将其设为NA
data$Z[data$Z <= 0] <- NA
# 添加随机NA
addNA <- sample(seq_len(length_measurement), length_measurement/3)
data$Z[addNA] <- NA

# 绘制热图
ggplot(data, aes(Y, X, fill = Z)) + 
  geom_tile() +
  scale_fill_gradient(
    name = "measurement of interest",
    low = "white", 
    high = "red",
    trans = "log10",  # 对fill值应用log10转换
    breaks = 10^c(-4.5, -4, -3.5, -3),  # 设置原始数值的breaks
    labels = scales::math_format(expr = 10^.x)  # 将breaks转换为10^x的数学格式
  ) +
  labs(x = "Y", y = "X")

方法2：提前对Z做log10转换，手动设置breaks和labels #

如果坚持在aes里使用log10(Z)，需要先确保Z都是正数，再对应设置log10后的breaks：

library(tidyverse)

# 创建示例数据集
x <- LETTERS[1:20]
y <- paste0("var", seq(1,20))
data <- expand.grid(X=x, Y=y)
length_measurement <- 400
data$Z <- rnorm(400, 0.000025, 0.0005)

# 处理Z的负数和0
data$Z[data$Z <= 0] <- NA
# 添加随机NA
addNA <- sample(seq_len(length_measurement), length_measurement/3)
data$Z[addNA] <- NA

# 提前计算log10(Z)
data$logZ <- log10(data$Z)

# 绘制热图
ggplot(data, aes(Y, X, fill = logZ)) + 
  geom_tile() +
  scale_fill_gradient(
    name = "log10(measurement of interest)",
    low = "white", 
    high = "red",
    breaks = c(-4.5, -4, -3.5, -3),  # 设置log10后的数值作为breaks
    labels = scales::math_format(expr = 10^.x)  # 将log值转换为10^x的格式
  ) +
  labs(x = "Y", y = "X")

关键要点 #

• 必须先处理Z的非正数：log10()无法处理0或负数，会返回NA，提前替换为NA可避免绘图异常；
• 两种方法的核心区别是数值转换的阶段：方法1由scale自动处理转换，方法2手动提前转换，对应breaks的设置逻辑不同；
• 使用scales::math_format(expr = 10^.x)可以让标签显示为数学表达式格式（如10⁻⁴.⁵），而非纯文本。

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