ggplot2：修复对数色标的图例

Question

我正在尝试在ggplot2绘制拟合模型效果，以替代effects包返回的图，并且在将连续分组因子映射到scale_color_gradient 。 这些问题源于分组因子具有偏斜分布的事实。 使用默认的颜色映射，大多数颜色是无法区分的，但是对数转换会混淆图例。 我看过一些相关的SO答案（ 是否有内置的方法可以在ggplot2中执行对数色标？ ， 用于对数数据的R色标？以及ggplot2中的对数色标会挤压某些图例编号 ）不太合适。

这是我的数据：

myEffs <- structure(list(PrimeShiftIndex = c(-4, -0.2, 4, -4, -0.2, 4,
-4, -0.2, 4, -4, -0.2, 4), PrimeVowelDur = c(0.03, 0.03, 0.03,
0.06, 0.06, 0.06, 0.09, 0.09, 0.09, 1.59, 1.59, 1.59), fit = c(-0.184306629528313,
-0.164313919815862, -0.142216714344205, -0.200749305969527, -0.178039844592615,
-0.152939913597082, -0.210367655099129, -0.186068995874736, -0.159212583047775,
-0.278488972243709, -0.242934925102426, -0.203638346683111),
    se = c(0.0437103286485701, 0.0342751848548937, 0.0446524040373885,
    0.0417352317881704, 0.0340007746839495, 0.042093900962637,
    0.0441609220226782, 0.0341565687974652, 0.0442166991273061,
    0.0995662189943997, 0.041203801253227, 0.0993299532144987
    ), lower = c(-0.269979086288845, -0.231493448847753, -0.229735643449126,
    -0.282550563276701, -0.24468152835471, -0.235444164230796,
    -0.296923277064622, -0.253016036857112, -0.245877528409642,
    -0.473639245768052, -0.323694575976204, -0.398325538129639
    ), upper = c(-0.0986341727677806, -0.0971343907839703, -0.0546977852392849,
    -0.118948048662354, -0.11139816083052, -0.0704356629633676,
    -0.123812033133635, -0.119121954892359, -0.0725476376859078,
    -0.0833386987193667, -0.162175274228647, -0.00895115523658357
    )), class = "data.frame", row.names = c(NA, -12L), .Names = c("PrimeShiftIndex",
"PrimeVowelDur", "fit", "se", "lower", "upper"))

在此，分组因子PrimeVowelDur歪斜右，其值在0.03 ， 0.06 ， 0.09 ，和1.59 。 这是获得可区分的颜色和可读图例的一些失败尝试（带有令人讨厌的宽线以突出显示颜色对比或缺少颜色对比）。

p <- ggplot(myEffs, aes(x=PrimeShiftIndex, y=fit, group=PrimeVowelDur, color=PrimeVowelDur)) +
  geom_line(size=6)
##Legend fine but line colors indistinguishable
p
##Missuse's suggestion yields the same issue
p + scale_color_gradientn(colors = colorRampPalette(colors = c("#132B43", "#56B1F7"))(nrow(myEffs)), 
                          values = scales::rescale(log(sort(myEffs$PrimeVowelDur))))


##Colors distinguishable but legend messed up
p + scale_color_gradient(trans="log")
##Using trans="log" with pre-defined breaks as per Gregor doesn't make legend much better
brks <- seq(0, 1.6, length.out=5)
p + scale_color_gradient(trans="log", breaks=brks, labels=brks)
##Nor does S Rivero's suggestion
p + scale_color_gradient(trans="log", breaks=brks, labels=brks, guide="legend")

我的直觉是ggplot应该能够自然地处理p + scale_color_gradient(trans="log")解决方案，而不会弄乱图例。 无论如何，我有一个建议的解决方案，但是我想看看是否有任何我想念的东西。

Answer 1

如果将基于连续分组因子的值的间断传递给breaks并为scale_color_gradient labels参数scale_color_gradient labels ，则图例将按预期显示。 无论您使用默认的颜色栏指南还是图例指南，此方法均有效：

brks2 <- sort(unique(myEffs$PrimeVowelDur))
p + scale_color_gradient(trans="log", breaks=brks2, labels=brks2)
p + scale_color_gradient(trans="log", breaks=brks2, labels=brks2, guide="legend")