如何在ggplot2中为离散变量z分组设置连续变量w的独立配色方案？

要实现为z的每个分组（A/B/C）分配不同的viridis配色方案来展示连续变量w，这里提供两种实用方法：

方法1：用ggnewscale添加多组独立颜色刻度 #

这是最直观的方案，借助ggnewscale包可以在同一张图里添加多套独立的颜色比例尺：

首先安装并加载依赖包：

install.packages("ggnewscale")
library(ggplot2)
library(ggnewscale)

然后按分组分层绘图，每组用不同的viridis选项，用new_scale_colour()分隔不同的颜色刻度：

ggplot() +
  # 分组A：用magma配色
  geom_point(data = subset(d, z == "A"), aes(x, y, colour = w)) +
  scale_colour_viridis_c(option = "magma", name = "w (分组A)") +
  # 开启新的颜色刻度组
  new_scale_colour() +
  # 分组B：用plasma配色
  geom_point(data = subset(d, z == "B"), aes(x, y, colour = w)) +
  scale_colour_viridis_c(option = "plasma", name = "w (分组B)") +
  # 再次开启新刻度
  new_scale_colour() +
  # 分组C：用默认viridis配色
  geom_point(data = subset(d, z == "C"), aes(x, y, colour = w)) +
  scale_colour_viridis_c(option = "viridis", name = "w (分组C)") +
  labs(x = "x", y = "y") +
  theme_minimal()

方法2：手动映射颜色（无需额外包） #

如果不想安装新包，可以手动将每个分组的w值映射到对应viridis调色板的颜色：

library(ggplot2)
library(viridis)

# 给每个分组指定对应的viridis选项
group_opts <- c(A = "magma", B = "plasma", C = "viridis")

# 为每条数据生成对应颜色
d$point_col <- mapply(function(w_val, z_group) {
  # 把当前分组的w值归一化到0-1区间
  norm_w <- (w_val - min(d$w[d$z == z_group])) / (max(d$w[d$z == z_group]) - min(d$w[d$z == z_group]))
  # 从对应调色板中取色
  viridis_pal(option = group_opts[z_group])(100)[cut(norm_w, 100, labels = FALSE)]
}, d$w, d$z)

# 绘图
ggplot(d, aes(x, y)) +
  geom_point(colour = d$point_col) +
  # 手动构建图例（需要自行调整细节）
  guides(colour = guide_legend(
    override.aes = list(colour = c(
      viridis_pal(option = "magma")(5),
      viridis_pal(option = "plasma")(5),
      viridis_pal(option = "viridis")(5)
    )),
    labels = c(rep("分组A",5), rep("分组B",5), rep("分组C",5))
  )) +
  labs(x = "x", y = "y") +
  theme_minimal()

注意：方法2的图例需要手动调整，不如方法1灵活，优先推荐ggnewscale方案，它能自动生成清晰的分组颜色刻度，可读性更强。

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