2012-09-14 45 views
19

He estado tratando de utilizar ggplot2 para producir una trama similar a este gráfica R:Como sombrear una región bajo una curva utilizando 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") 

Esto es por lo que yo he conseguido 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)) 

Me gustaría sombrear la región de la cola más allá de x = 0,12. ¿Cómo haría esto usando ggplot o qplot?

En general, ¿cómo se sombrea un subconjunto debajo de la curva, ya sea una cola, o entre dos líneas arbitrarias que dividen la región en áreas distintas?

Gracias por cualquier consejo.

Respuesta

14

crear un polígono con el área que desea sombrear

#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 respuesta ayuda de otra forma también. No estaba pensando en el modo de ggplot2 y creando un subconjunto de datos explícitos. Intentaba hacer esto desde un punto de vista de objeto puramente gráfico. –

+0

Creo que quiso decir (tenga en cuenta la ubicación del último paréntesis de cierre) ... shade <- rbind (c (0.12,0), subconjunto (MyDF, x> 0.12), c (MyDF [nrow (MyDF), " X "], 0)) –

+0

@ceilingcat Tienes razón, voy a editar la respuesta –

Cuestiones relacionadas