Come rendere i grafici ggplot2 più belli?

Question

Mi hanno generato la trama di seguito utilizzando il codice R che lo segue:

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

Ora mi piacerebbe fare la trama un po 'più carina:

1. rendere lo spettacolo di lunghezza asse x ogni unità 5000 (anziché ogni 20000)
2. Aggiungere i valori x sopra i tre picchi (circa 3000,5000 e 35000).

Come posso farlo?

aggiornamento in risposta a James:

Answer 1

ne dite:

(prima creare un esempio riproducibile)

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)))) 

(roba carina per trovare le posizioni di picco e le altezze)

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 

(disegnare la trama)

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) 

voi probabilmente vorrai modificare leggermente l'altezza verticale delle etichette

Answer 2

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

si potrebbe anche mettere la dichiarazione xlim qui, utilizzando limits, ad esempio,

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

2: Per le etichette è possibile utilizzare annotate() o geom_text(). Vedi il post this per gli esempi. Dovresti comunque calcolare i valori per questo.