我使用后面的R代码生成了以下图:
ggplot(lengths, aes(length, fill = library)) + geom_density(alpha = 0.2) + coord_cartesian(xlim = c(0, 60000))
现在我想让情节变得更漂亮:
我该怎么做?
更新 回应詹姆斯:
答案 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))
。
你也可以使用xlim
将limits
声明放在这里,例如,
+ scale_x_continuous(breaks=rep(5000,12),limits=c(0,60000))
2:对于标签,您可以使用annotate()
或geom_text()
。有关示例,请参阅this帖子。你必须自己计算这个值。