2012-09-14 14 views
19

Ho cercato di usare ggplot2 per produrre una trama simile a questo R grafica:Come ombreggiare una regione sotto una curva con ggplot2

xv<-seq(0,4,0.01) 
yv<-dnorm(xv,2,0.5) 
plot(xv,yv,type="l") 
polygon(c(xv[xv<=1.5],1.5),c(yv[xv<=1.5],yv[xv==0]),col="grey") 

Questo è quanto ho ottenuto con ggplot2:

x<-seq(0.0,0.1699,0.0001) 
ytop<-dnorm(0.12,0.08,0.02) 
MyDF<-data.frame(x=x,y=dnorm(x,0.08,0.02)) 
p<-qplot(x=MyDF$x,y=MyDF$y,geom="line") 
p+geom_segment(aes(x=0.12,y=0,xend=0.12,yend=ytop)) 

Vorrei ombreggiare la regione della coda oltre x = 0,12. Come lo farei usando ggplot o qplot?

In generale, in che modo si ombreggia qualsiasi sottoinsieme della curva, se una coda o tra due linee arbitrarie che dividono la regione in aree distinte?

Grazie per qualsiasi consiglio.

risposta

14

creare un poligono con l'area che si desidera ombra

#First subst the data and add the coordinates to make it shade to y = 0 
shade <- rbind(c(0.12,0), subset(MyDF, x > 0.12), c(MyDF[nrow(MyDF), "X"], 0)) 

#Then use this new data.frame with geom_polygon 
p + geom_segment(aes(x=0.12,y=0,xend=0.12,yend=ytop)) + 
    geom_polygon(data = shade, aes(x, y)) 

enter image description here

+0

La risposta aiuta in un altro modo troppo. Non stavo pensando nella modalità di ggplot2 e creando un sottoinsieme di dati esplicito. Stavo cercando di farlo funzionare da un punto di vista puramente grafico. –

+0

Penso che intendessi (annota la posizione dell'ultima parentesi chiusa) ... shade <- rbind (c (0,12,0), sottoinsieme (MyDF, x> 0,12), c (MyDF [nrow (MyDF), " X "], 0)) –

+0

@ceilingcat Hai ragione, modificherò la risposta –

Problemi correlati