5

我是ggplot的新手,所以请耐心等待。我正在绘制 35 个小区域的增长预测,即使使用了很棒的directlabels图书馆,这对于一个地块来说也是一个不健康的数字。但是,我需要所有系列进行初步筛选。

挑战在于使其具有可读性。我找到了@Ben Bolker 使用大量不同颜色的修复程序,但在改变线型时遇到了麻烦。35 系列不需要是唯一的,但我想使用 12 种不同的类型来使单个系列更易于阅读。

我的计划是创建一个包含 12 种可能类型的 35 个元素的随机列表,并将其作为 linetype 参数传递,但我无法让它工作,并出现错误:

Error: Aesthetics must either be length one, or the same length as the dataProblems:lty

线型列表中有 35 个值。当然,我希望类型、颜色和所有内容都反映在图例中。

融化的数据看起来像这样;35 个系列中每个系列的 9 年观察结果:

> simulation_long_index[16:24,]    
      year    geography    value
16    2018    sfr_2    101.1871
17    2019    sfr_2    101.1678
18    2020    sfr_2    101.2044
19    2012    sfr_3    100.0000
20    2013    sfr_3    100.1038
21    2014    sfr_3    100.2561
22    2015    sfr_3    100.0631
23    2016    sfr_3    100.8071
24    2017    sfr_3    101.2405    

到目前为止,这是我的代码:

lty <- data.frame(lty=letters[1:12][sample(1:12, 35,replace=T)])

g3<-ggplot(data=simulation_long_index,
   aes(
     x=as.factor(year), 
     y=value, 
     colour=geography,
     group=geography,
     linetype=lty$lty))+
       geom_line(size=.65) + 
       scale_colour_manual(values=manyColors(35)) +
     geom_point(size=2.5) +
     opts(title="growth")+
     xlab("Year") + 
     ylab(paste("Indexed Value (Rel. to 2012")) + 
     opts(axis.text.x=theme_text(angle=90, hjust=0))

print(g3)

添加

    scale_linetype_manual("",values=lty$lty) +

在 scale_color_manual 而不是 linetype 参数生成图表之后,但线条都是相同的。那么,如何让线条因大系列计数而变化?

示例图

4

1 回答 1

10

using 的技巧scale_..._manual通常是发送一个命名向量作为value参数。这个setNames功能很好

一、一些虚拟数据

## some dummy data 
simulations<- expand.grid(year = 2012:2020, geography = paste0('a',1:35))
library(plyr)
library(RColorBrewer)
simulation_long_index <- ddply(simulations, .(geography), mutate, 
  value = (year-2012) * runif(1,-2, 2) + rnorm(9, mean = 0, sd = runif(1, 1, 3)))
## create a manyColors function 
manyColors <- colorRampPalette(brewer.pal(name = 'Set3',n=11))

接下来我们创建一个向量,它是 1:12 的随机样本(有替换),并将名称设置为与geography 变量相同

lty <- setNames(sample(1:12,35,T), levels(simulation_long_index$geography))

这就是它的样子

lty
## a1  a2  a3  a4  a5  a6  a7  a8  a9 a10 a11 a12 a13 a14 a15 a16 
## 7   5   8  11   2  10   3   2   5   4   6   6  11   8   2   2 
## a17 a18 a19 a20 a21 a22 a23 a24 a25 a26 a27 a28 a29 a30 a31 a32 
## 12   7   6   8  11   5   1   1   8  12   8   1  12   2   3   5 
## a33 a34 a35 
#7   1   3 

现在您可以line_type = geography结合使用scale_linetype_manual(values = lty)

ggplot(data=simulation_long_index,
        aes(
          x=as.factor(year), 
          y=value, 
          colour=geography,
          group=geography,
          linetype = geography))+
            geom_line(size=.65) + 
            scale_colour_manual(values=manyColors(35)) +
            geom_point(size=2.5) +
            opts(title="growth")+
            xlab("Year") + 
            ylab(paste("Indexed Value (Rel. to 2012")) + 
            opts(axis.text.x=theme_text(angle=90, hjust=0)) +
            scale_linetype_manual(values = lty)

这给了你

在此处输入图像描述

顺便说一句,您真的想将年份绘制为因子变量吗?

于 2012-07-30T02:34:59.887 回答