需要从ggplot geom_histogram中提取数据 [英] Need to extract data from the ggplot geom_histogram

问题描述

您好我已经使用数据框中的ggplot2 geom_histogram函数绘制了一个图表，参见下面的示例并链接到ggplot直方图需要使用嵌套的ddply函数和facet wrap来标记每个geom_vline

现在我需要制作一个数据框，其中包含上图中使用的汇总数据。谁能帮忙？提前致谢。 G

  Sector2家庭年长度
 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 
  

解决方案

为了得到实际绘制的值，可以使用函数 ggplot_build（），其中参数就是您的绘图。 b
$ b

  p < -  ggplot（mtcars，aes（mpg））+ geom_histogram（）+ 
 facet_wrap（〜cyl）+ geom_vline = data.frame（x = c（20,30）），aes（xintercept = x））
 
 pg < -  ggplot_build（p）
  

这将使列表和一个子列表被命名为 data code>。这个子列表包含数据框，其中使用了plot中使用的值，例如，对于histrogramm它包含 y 值（与 count ）。如果您使用构面，则列 PANEL 显示使用哪些构面值。如果在你的图中有多个 geom _ ，那么数据将包含每个数据帧 - 在我的例子中，有一个数据帧用于histogramm，另一个用于vlines。

$ p $ 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]]）
x侦听PANEL组xnd 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





Hi I have made a plot using the ggplot2 geom_histogram function from a data frame see sample below and link to the ggplot histogram Need to label each geom_vline with the factors using a nested ddply function and facet wrap

Now I need to make a data frame that contains the summarized data used in the plot above. Can anyone help? Thanks in advance. G

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

解决方案

To get values actually plotted you can use function ggplot_build() where argument is your plot.

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)

This will make list and one of sublists is named data. This sublist contains dataframe with values used in plot, for example, for histrogramm it contains y values (the same as count). If you use facets then column PANEL shows in which facet values are used. If there are more than one geom_ in your plot then data will contains dataframes for each - in my example there is one dataframe for histogramm and another for 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

这篇关于需要从ggplot geom_histogram中提取数据的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持IT屋！