我已使用其后的R代码生成了以下图:
ggplot(lengths, aes(length, fill = library)) + geom_density(alpha = 0.2) + coord_cartesian(xlim = c(0, 60000))
现在,我想使情节更漂亮:
5000个单位(而不是每20000个单位)
峰值(大约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/