如何使ggplot2图更漂亮?

时间:2010-08-31 07:26:26

标签: r plot ggplot2

我使用后面的R代码生成了以下图: alt text

ggplot(lengths, aes(length, fill = library)) + geom_density(alpha = 0.2) + coord_cartesian(xlim = c(0, 60000)) 

现在我想让情节变得更漂亮:

  1. 使x轴显示每个长度 5000单位(而不是每20000)
  2. 在三者之上添加x值     峰值(约3000,5000和35000)。
  3. 我该怎么做?

    更新 回应詹姆斯: alt text

2 个答案:

答案 0 :(得分:11)

怎么样:

(首先创建一个可重现的例子)

set.seed(1001)
lengths <- data.frame(length=c(rgamma(1000,shape=10,scale=500),
                    10000+rgamma(1000,shape=5,scale=700),
                    rnorm(500,mean=30000,sd=2000)),
                  library=factor(rep(2:1,c(2000,500))))

(寻找高峰位置和高度的可爱东西)

peakfun <- function(x) {
  d <- density(x$length)
  peaks <- which(diff(sign(diff(d$y)))==-2)
  data.frame(x=d$x[peaks],y=d$y[peaks])
}

peakdat <- ddply(lengths,.(library),peakfun)
peakdat <- peakdat[-1,] ## drop spurious peak

(绘制情节)

library(ggplot2)
ggplot(lengths, aes(length, fill = library)) +
  geom_density(alpha = 0.2) +
  scale_x_continuous(limits = c(0,60000),
                     breaks = seq(0,60000,by=5000))+
  geom_text(data=peakdat,aes(x=x,y=y,label=round(x)),vjust=1)

您可能想稍微调整标签的垂直高度

答案 1 :(得分:6)

1:+ scale_x_continuous(breaks=rep(5000,12))

你也可以使用xlimlimits声明放在这里,例如,

+ scale_x_continuous(breaks=rep(5000,12),limits=c(0,60000))

2:对于标签,您可以使用annotate()geom_text()。有关示例,请参阅this帖子。你必须自己计算这个值。