2

我有一个图，直方图上有几个观察结果。观察结果是按顺序收集的，我需要观察它们收集的顺序。使用scale_colour_brewer. 问题是，连续 brewer 调色板的最大长度为 9。我有多达 20 个观察值的示例，但我不确定如何使用插值颜色。下面是一些代码，它用少于 10 个点演示了我想要的输出。

# Setting this to be > 9 will cause a warning and not produce the desired result.
observations = 9
subset <-1:observations
res = data.frame(x_data = rnorm(5000),TestID=1:5000)
ggplot(res,aes(x=x_data)) + 
  stat_bin(aes(y=..density..))+
  stat_density(colour="blue", fill=NA)+
  geom_point(data = res[res$TestID %in% subset,], 
             aes(x = x_data, 
                 y = 0, 
                 colour = as.factor(res$TestID[res$TestID %in% subset])
             ),
             size = 5) +scale_colour_brewer("Fancy title", type="seq", palette='Reds')

我知道随着观察次数的增加，这个图将变得难以阅读。但是，我相信使用多达 20 种颜色，应该可以在我的应用程序中解释结果。

4

1 回答 1

5

扩展我的评论，您需要使用colorRampPalette：

library(RColorBrewer)
blues_fun <- colorRampPalette(brewer.pal(9,"Blues"))
> blues_fun(20)
 [1] "#F7FBFF" "#ECF4FB" "#E1EDF8" "#D7E6F4" "#CDE0F1" "#C1D9ED" "#B0D2E7" "#A0CAE1" "#8BBFDC" "#75B3D8" "#62A8D2" "#519CCB"
[13] "#4090C5" "#3282BD" "#2474B6" "#1966AD" "#0E59A2" "#084B94" "#083D7F" "#08306B"

然后通过以下方式构建规模scale_colour_manual：

ggplot(res,aes(x=x_data)) + 
  stat_bin(aes(y=..density..))+
  stat_density(colour="blue", fill=NA)+
  geom_point(data = res[res$TestID %in% subset,], 
             aes(x = x_data, 
                 y = 0, 
                 colour = as.factor(res$TestID[res$TestID %in% subset])
             ),
             size = 5) + 
  scale_colour_manual("Fancy title",values = blues_fun(9))

您只需将生成的颜色交给values参数即可。

于 2013-06-04T16:58:19.303 回答