4

如何根据数据集中的列覆盖图例指南的aes大小值?ggplot2

请参阅此示例(编辑 2:添加了试用 C,并更改了行大小以使用对数刻度):

library(data.table)
set.seed(26798)

dt<-rbind(data.table(Trial="A",Value=rweibull(1000,1.0,0.5)),
      data.table(Trial="B",Value=rweibull(100,1.2,0.75)),
      data.table(Trial="C",Value=rweibull(10,1.3,0.8)))

# Add a count and something like a cumulative distribution:
dt2<-dt[order(Trial,Value),list(Value,N=.N),by=Trial][,list(Value,N,y=1-cumsum(N)/sum(N)),by=Trial]
dt2
##      Trial        Value    N     y
##   1:     A 0.0003628745 1000 0.999
##   2:     A 0.0013002615 1000 0.998
##   3:     A 0.0017002173 1000 0.997
##   4:     A 0.0022597343 1000 0.996
##   5:     A 0.0026608082 1000 0.995
##  ---                              
##1096:     B 1.6821827814  100 0.040
##1097:     B 2.2431595707  100 0.030
##1098:     B 2.5122479833  100 0.020
##1099:     B 2.5519954416  100 0.010
##1100:     B 2.6848412995  100 0.000

ggplot(dt2) +
  geom_line(aes(x=Value,y=y,group=Trial,color=Trial,size=N)) +
  scale_size(range=c(0.1, 2), trans="log") +
  guides(size=F, color=guide_legend(override.aes=list(size=2)))

三个试验的情节

我希望指南图例中每个 Trial 值的线条粗细与图中的线条相匹配(即“A”应该是粗的,“B”应该是细的)。 编辑1: @Arun 和@ChelseaE 为手动调整每个厚度给出了很好的建议,但是我的实际数据集有很多因子水平并且不断变化,所以我需要它是“动态的”。

@DidzisElferts 对类似问题的回答(控制 ggplot2 图例外观而不影响绘图)显示了如何将大小设置为静态值。上面示例的最后一行中的size=2部分让我可以更改图例的线条大小,但我希望它与图中线条的大小相匹配。改用似乎合乎逻辑,但它给出size=N了错误“找不到对象'N'”。什么是正确的语法?

期望的输出:

带有所需图例的三个试验的情节

4

2 回答 2

6

您应该为 A 和 B 相应地设置尺寸。您只设置了 1 个尺寸。尝试这个:

p <- ggplot(dt2) +
geom_line(aes(x=Value,y=y,group=Trial,color=Trial,size=N)) +
scale_size(range=c(0.1, 2)) +
guides(size=FALSE, color=guide_legend(override.aes=list(size=c(2, .1))))

按照OP的评论:

好的,在这种情况下,您将不得不做更多的工作(可能有更简单的方法;我现在想不出它们,如果有的话)。

scales <- c(0.1, 2) # the range you want: min, max
vals <- summary(lm(scales ~ c(min(dt2$N), max(dt2$N))))$coefficients[,1]
sizes <- vals[2] * unique(dt2$N) + vals[1]

ggplot(dt2) +
geom_line(aes(x=Value,y=y,group=Trial,color=Trial,size=N)) +
scale_size(range=scales) +
guides(size=FALSE, color=guide_legend(override.aes=list(size=sizes)))

这应该有效。试试看,如果你有问题,请告诉我。

于 2013-08-21T17:59:02.670 回答
1

尝试在 guide_legend 中添加尺寸范围:

ggplot(dt2) + geom_line(aes(x = Value, y = y, group = Trial, color = Trial, size = N)) + scale_size(range = c(0.1, 2)) + guides(size = F, color = guide_legend(override.aes = list(size = range(0.1,2))))

编辑:(可能有效,不确定)

您也可以尝试为 N ( N <- dt2$N) 创建一个向量,然后使用size = N

希望这可以帮助。

于 2013-08-21T19:06:29.050 回答