将颜色作为变量传递给aes_string

时间:2013-12-18 05:33:47

标签: r ggplot2

我想在一列上绘制数据矩阵的几列,所以不是单独绘制每一列,而是使用循环。问题是如果要使用循环,要绘制的列名称和颜色必须是可变的。我尝试了以下方法:

allDs <- sort(unique(plotdata$D))
p <- ggplot(plotdata, aes(SpaceWidth))
for (thisD in allDs) {
    tlColName <- paste("M2D", thisD, "Tl", sep="")
    colorName <- paste("D", thisD, sep="")
    p <- p + geom_line(data = plotdata[!is.na(plotdata[[tlColName]]),], aes_string(y = tlColName, color = colorName))
}
p <- p + scale_colour_manual("Legend", values = c("D2" = "blue", "D3" = "red", "D4" = "green", "D6" = "violet", "D7" = "yellow" ))
p <- p + scale_x_log10(breaks = composite$SpaceWidth)
p <- p + facet_wrap(~ Drawn, ncol = 3)
p <- p + labs(title = "Fu plot", y = "MTN")
p

但是当我运行这个时,我收到以下错误:

Error in eval(expr, envir, enclos) : object 'D2' not found

如何将颜色(或任何其他值)作为变量传递给aes_string?提前谢谢。

测试数据可用here

1 个答案:

答案 0 :(得分:3)

错误的原因非常明显:原始数据集中不存在D2。请注意,您可以将颜色直接映射到变量D,因此您的colorName构造是多余的。看看这个:

allDs <- sort(unique(plotdata$D))
plotdata$D <- as.factor(plotdata$D)
p <- ggplot(plotdata, aes(SpaceWidth, color=D))
for (thisD in allDs) {
   tlColName <- paste("M2D", thisD, "Tl", sep="")
   p <- p + geom_line(data = plotdata[!is.na(plotdata[[tlColName]]),], 
            aes_string(y = tlColName))
}
p <- p + scale_colour_manual("Legend", 
         values = c("blue", "red", "green", "violet", "yellow"))
p <- p + scale_x_log10(breaks = plotdata$SpaceWidth)
p <- p + facet_wrap(~ D, ncol = 3)
p <- p + labs(title = "Fu plot", y = "MTN")
p

请注意,要正确映射颜色,您需要先将其转换为因子。 enter image description here

<强> UPD: 好吧,让我告诉你如何摆脱for循环,这通常不是一个好习惯。

library(reshape2)
melt.plotdata <- melt(plotdata, id.vars=c("SpaceWidth", "D"))
melt.plotdata <- melt.plotdata[order(melt.plotdata$SpaceWidth), ]
melt.plotdata <- na.omit(melt.plotdata)
q <- ggplot(melt.plotdata, aes(SpaceWidth, value, colour=variable)) + geom_path()
q + scale_colour_manual("Legend", 
                        values = c("blue", "red", "green", "violet", "yellow")) +
  scale_x_log10(breaks = melt.plotdata$SpaceWidth) +
  facet_wrap(~ D, ncol = 3) + 
  labs(title = "Fu plot", y = "MTN")

情节将与我上面发布的情节相同。