histograms - stat_bin

Necesidad de extraer datos del ggplot geom_histogram

Question

Hola, he hecho una gráfica usando la función ggplot2 geom_histogram de un marco de datos. Vea la muestra a continuación y vincule al histograma ggplot. Necesite etiquetar cada geom_vline con los factores usando una función ddply anidada y un ajuste de facetas.

Ahora necesito hacer un marco de datos que contenga los datos resumidos utilizados en el diagrama anterior. ¿Alguien puede ayudar? Gracias por adelantado. GRAMO

Sector2 Family Year Length
BUN Acroporidae 2010 332.1300496
BUN Poritidae 2011 141.1467966
BUN Acroporidae 2012 127.479
BUN Acroporidae 2013 142.5940556
MUR Faviidae 2010 304.0405
MUR Faviidae 2011 423.152
MUR Pocilloporidae 2012 576.0295
MUR Poritidae 2013 123.8936667
NTH Faviidae 2010 60.494
NTH Faviidae 2011 27.427
NTH Pocilloporidae 2012 270.475
NTH Poritidae 2013 363.4635

Answer 1

Para obtener los valores realmente trazados, puede usar la función ggplot_build() donde argumento es su trazado.

p <- ggplot(mtcars,aes(mpg))+geom_histogram()+
facet_wrap(~cyl)+geom_vline(data=data.frame(x=c(20,30)),aes(xintercept=x))
pg <- ggplot_build(p)

Esto hará que la lista y una de las sublistas se denomine data . Esta sublista contiene el marco de datos con los valores utilizados en la trama, por ejemplo, para histrogramm contiene valores y (igual que el count ). Si usa facetas, la columna PANEL muestra en qué valores de faceta se utilizan. Si hay más de un geom_ en su trazado, los datos contendrán marcos de datos para cada uno; en mi ejemplo, hay un marco de datos para histogramm y otro para vlines.

head(pg$data[[1]])
y count x ndensity ncount density PANEL group ymin ymax
1 0 0 9.791667 0 0 0 1 1 0 0
2 0 0 10.575000 0 0 0 1 1 0 0
3 0 0 11.358333 0 0 0 1 1 0 0
4 0 0 12.141667 0 0 0 1 1 0 0
5 0 0 12.925000 0 0 0 1 1 0 0
6 0 0 13.708333 0 0 0 1 1 0 0
xmin xmax
1 9.40000 10.18333
2 10.18333 10.96667
3 10.96667 11.75000
4 11.75000 12.53333
5 12.53333 13.31667
6 13.31667 14.10000
head(pg$data[[2]])
xintercept PANEL group xend x
1 20 1 1 20 20
2 30 1 1 30 30
3 20 2 2 20 20
4 30 2 2 30 30
5 20 3 3 20 20
6 30 3 3 30 30

Answer 2

Si solo necesita datos, parece que layer_data está diseñado precisamente para esto:

layer_data(p, 1)

le dará los datos de la primera capa ( ggplot_build(p)$data[[1]] como ggplot_build(p)$data[[1]] )