如何从 R 中的 facet_grid 图中获取真实比例？

Question

我正在尝试添加此中显示的字幕。

因此，当我使用 facet_grid 时，我需要绘图的真实比例（x 和 y 轴）。我知道我可以使用 layer_data，因为它保存了情节中的所有内容......但是，它并不准确，因为当我尝试使用 min 和 max 建立限制时根据该输出，情节发生了变化。

这里有一个例子：

library(ggplot2)
library(dplyr)

val1 <- c(2.1490626,2.2035281,1.5927854,3.1399245,2.3967338,3.7915825,4.6691277,3.0727319,2.9230937,2.6239759,3.7664386,4.0160378,1.2500835,4.7648343,0.0000000,5.6740227,2.7510256,3.0709322,2.7998003,4.0809085,2.5178086,5.9713330,2.7779843,3.6724801,4.2648527,3.6841084,2.5597235,3.8477471,2.6587736,2.2742209,4.5862788,6.1989269,4.1167091,3.1769325,4.2404515,5.3627032,4.1576810,4.3387921,1.4024381,0.0000000,4.3999099,3.4381837,4.8269218,2.6308474,5.3481382,4.9549753,4.5389650,1.3002293,2.8648220,2.4015338,2.0962332,2.6774765,3.0581759,2.5786137,5.0539080,3.8545796,4.3429043,4.2233248,2.0434363,4.5980727)
val2 <- c(3.7691229,3.6478055,0.5435826,1.9665861,3.0802654,1.2248374,1.7311236,2.2492826,2.2365337,1.5726119,2.0147144,2.3550348,1.9527204,3.3689502,1.7847986,3.5901329,1.6833872,3.4240479,1.8372175,0.0000000,2.5701453,3.6551315,4.0327091,3.8781182)

df1 <- data.frame(value = val1)   
df2 <- data.frame(value = val2)   

data <- bind_rows(lst(df1, df2), .id = 'id')
data$Sex <- rep(c("Male", "Female"), times=84/2)

p <- data %>% 
  ggplot(aes(value)) +
  geom_density(lwd = 1.2, colour="red", show.legend = FALSE) +
  geom_histogram(aes(y=..density.., fill = id), bins=10, col="black", alpha=0.2) +
  facet_grid(id ~ Sex ) +
  xlab("type_data") + 
  ylab("Density") +
  ggtitle("title") +
  guides(fill=guide_legend(title="legend_title")) +
  theme(strip.text.y = element_blank())
  
p

plot_info <- layer_data(p)
> min(plot_info$density)
[1] 7.166349e-09
> max(plot_info$density)
[1] 0.5738021

正如您在图中看到的那样，y 轴从 0 开始，如果结束则减少 0.7 左右。但是，最大密度是0.57.

如果我尝试使用来自 layer_data 的信息：

p + coord_cartesian(clip="off", ylim=c(min(plot_info$density), max(plot_info$density)), 
                  xlim = c(min(plot_info$x), max(plot_info$x)))

剧情完全改变

有谁知道如何获得 ggplot2 和 facet_grid 正在使用的秤？我需要密度信息 (y_axis) 和来自 x_axis.

的信息

Answer 1

是的，要直接获取比例，请使用 layer_scales(p)，它会为您提供轴的范围，而不仅仅是数据的范围，这是您从 layer_data(p)[= 中获得的16=]

p + coord_cartesian(clip = "off", 
                    ylim = layer_scales(p)$y$range$range, 
                    xlim = layer_scales(p)$x$range$range)

或者，将这个问题与你的最后一个问题结合起来，在绘图面板之外添加文本标签，你的结果可能类似于：

p + coord_cartesian(clip = "off", 
                    ylim = layer_scales(p)$y$range$range, 
                    xlim = layer_scales(p)$x$range$range) +
  geom_text(data = data.frame(value = c(0, 6), id = c("df2", "df2"),
                              Sex = c('Female', 'Male')),
            aes(y = -0.15, label = c('Female', 'Male')))

Answer 2

这有帮助吗？

?layer_data
summary(layer_data(p, i = 2))

i是你要的图层return

可以最小 xmin 和最大 xmax 等等

如何从 R 中的 facet_grid 图中获取真实比例？

How can I get the real scale from a facet_grid plot in R?

r

scale

ggplot2

density-plot

facet-grid