我已使用其后的R代码生成了以下图:

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

现在,我想使情节更漂亮:
  • 使x轴每次显示长度
    5000个单位(而不是每20000个单位)
  • 在三个值的顶部添加x值
    峰值(大约3000,5000和35000)。

  • 我怎样才能做到这一点?

    更新
    回应詹姆斯:

    最佳答案

    怎么样:

    (首先创建一个可复制的示例)

    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)
    

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

    关于r - 如何使ggplot2绘图更漂亮?,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/3606879/

    10-12 19:11