ggplot 在同一面板中组合线和区域(平面网格)

Posted

技术标签:

【中文标题】ggplot 在同一面板中组合线和区域(平面网格)【英文标题】:ggplot combine line and area in same panel (facet-grid) 【发布时间】:2016-09-17 03:37:35 【问题描述】:

我正在尝试创建一个构面网格。它将有三行四列。我希望两行是折线图,另一行是堆叠的 geom_area。除非我绘制这个,它总是出现四行——两行,两个区域。我想把这两个区域放在一起。我该如何解决这个问题?

library(ggplot2)
rawdata_Fig1 <- read.table(text="
Date    Number  mMP Day Rep Parameter   FeType  Value
13/03/2016  1   0   2   1   pH  NA  3.11
13/03/2016  2   0   2   2   pH  NA  3.13
13/03/2016  3   0   2   3   pH  NA  3.43
13/03/2016  4   0   2   4   pH  NA  3.49
21/04/2016  0   10  41  1   pH  NA  2.87
21/04/2016  0   10  41  2   pH  NA  2.87
21/04/2016  0   10  41  3   pH  NA  2.87
21/04/2016  0   10  41  4   pH  NA  2.85
13/03/2016  1   0   2   1   Fe2 Fe2 0.965538462
13/03/2016  57  1   2   1   Fe2 Fe2 0.217846154
13/03/2016  113 5   2   1   Fe2 Fe2 0.300153846
13/03/2016  169 10  2   1   Fe2 Fe2 0.226307692
21/04/2016  0   0   41  1   Fe2 Fe2 9.277692308
21/04/2016  0   1   41  1   Fe2 Fe2 13.92
21/04/2016  0   5   41  1   Fe2 Fe2 4.181538462
21/04/2016  0   10  41  1   Fe2 Fe2 1.066153846
13/03/2016  1   0   2   1   'PO4 uM'    NA  9.39
13/03/2016  2   0   2   2   'PO4 uM'    NA  0.39
13/03/2016  3   0   2   3   'PO4 uM'    NA  -0.06
13/03/2016  4   0   2   4   'PO4 uM'    NA  -0.06
13/03/2016  57  1   2   1   'PO4 uM'    NA  0.09
13/03/2016  58  1   2   2   'PO4 uM'    NA  -0.06
13/03/2016  59  1   2   3   'PO4 uM'    NA  0.54
13/03/2016  60  1   2   4   'PO4 uM'    NA  0.69
13/03/2016  113 5   2   1   'PO4 uM'    NA  113.19
13/03/2016  114 5   2   2   'PO4 uM'    NA  120.54
13/03/2016  115 5   2   3   'PO4 uM'    NA  144.24
13/03/2016  116 5   2   4   'PO4 uM'    NA  101.34
13/03/2016  169 10  2   1   'PO4 uM'    NA  235.17
13/03/2016  170 10  2   2   'PO4 uM'    NA  263.52
13/03/2016  171 10  2   3   'PO4 uM'    NA  252.72
13/03/2016  172 10  2   4   'PO4 uM'    NA  261.72
21/04/2016  0   10  41  1   'PO4 uM'    NA  17.04
21/04/2016  0   10  41  2   'PO4 uM'    NA  26.79
21/04/2016  0   10  41  3   'PO4 uM'    NA  29.49
21/04/2016  0   10  41  4   'PO4 uM'    NA  23.64
13/03/2016  1   0   2   1   TotFe   TotFe   1.316027906
13/03/2016  57  1   2   1   TotFe   TotFe   0.233872429
13/03/2016  113 5   2   1   TotFe   TotFe   0.273511824
13/03/2016  169 10  2   1   TotFe   TotFe   0.158557579
21/04/2016  0   0   41  1   TotFe   TotFe   12.51812087
21/04/2016  0   1   41  1   TotFe   TotFe   18.02799674
21/04/2016  0   5   41  1   TotFe   TotFe   6.167889825
21/04/2016  0   10  41  1   TotFe   TotFe   1.625215185
",header=TRUE)

rawdata_Fig1_Fe <- read.table(text="
Date    Number  mMP Day Rep Parameter Value
13/03/2016  1   0   2   1   Fe2     0.965538462
13/03/2016  57  1   2   1   Fe2     0.217846154
13/03/2016  113 5   2   1   Fe2     0.300153846
13/03/2016  169 10  2   1   Fe2     0.226307692
21/04/2016  0   0   41  1   Fe2     9.277692308
21/04/2016  0   1   41  1   Fe2     13.92
21/04/2016  0   5   41  1   Fe2     4.181538462
21/04/2016  0   10  41  1   Fe2     1.066153846
13/03/2016  1   0   2   1   TotFe       1.316027906
13/03/2016  57  1   2   1   TotFe       0.233872429
13/03/2016  113 5   2   1   TotFe       0.273511824
13/03/2016  169 10  2   1   TotFe       0.158557579
21/04/2016  0   0   41  1   TotFe       12.51812087
21/04/2016  0   1   41  1   TotFe       18.02799674
21/04/2016  0   5   41  1   TotFe       6.167889825
21/04/2016  0   10  41  1   TotFe       1.625215185
",header=TRUE)

Fig1 <-ggplot(rawdata_Fig1,aes(x=Day,y=Value, group=Rep))
Fig1Pan=Fig1 + theme_bw()

Fig1panel=Fig1Pan + facet_grid(Parameter ~ mMP, scales="free") +
  geom_line(data=subset(rawdata_Fig1, Parameter == "pH")) +
  geom_line(data=subset(rawdata_Fig1, Parameter == "PO4 uM")) +
  geom_area(data=rawdata_Fig1_Fe, aes(y=pH, x=Day, fill=Parameter, stat='identity'))

Fig1panel

【问题讨论】:

你有8列9列标签,是多了一个还是一个两个字? 有一个额外的,现在已修复。谢谢 所以你想把两个区域的地块合并成一个地块,基本上是叠起来的?不是简单地重新排列行以使两个区域行彼此重叠? 正确,堆叠面积图。当不是 facet_grid 的一部分时,我可以自己执行此操作。但在结合 facet 和 stacking 时不会。我试过 position="stack",但可能不在正确的位置?我是 ggplot 的新手。 【参考方案1】:

您可以组合两个构面行并创建新的fillgroup 变量以匹配构面层

rawdata_Fig1$Parameter[rawdata_Fig1$Parameter %in% 'Fe2'] <- 'TotFe'
# rawdata_Fig1 <- droplevels(rawdata_Fig1)
rawdata_Fig1_Fe$Parameter2 <- rawdata_Fig1_Fe$Parameter
rawdata_Fig1_Fe$Parameter[rawdata_Fig1_Fe$Parameter %in% 'Fe2'] <- 'TotFe'

Fig1 <-ggplot(rawdata_Fig1,aes(x=Day,y=Value, group=Rep))
Fig1Pan=Fig1 + theme_bw()

Fig1Pan + facet_grid(Parameter ~ mMP, scales="free") +
  geom_line(data=subset(rawdata_Fig1, Parameter == "pH")) +
  geom_line(data=subset(rawdata_Fig1, Parameter == "'PO4 uM'")) +

  ## I only changed these lines
  geom_area(data=rawdata_Fig1_Fe, aes(y=Value, x=Day, group = Parameter2, fill = Parameter2)) +
  labs(fill = 'Parameter')

【讨论】:

谢谢!!我花了一些时间弄清楚(了解你做了什么,代码很好),但它有效!

以上是关于ggplot 在同一面板中组合线和区域(平面网格)的主要内容,如果未能解决你的问题,请参考以下文章

R语言ggplot2可视化:ggplot2可视化基本散点图(scatter plot)通过在theme_bw中指定参数base_size来改变轴标签的大小并控制网格线和轴标签的大小

ggplot:生成具有多个系列的平面网格图

网格面板中的 EXT.NET 组合框显示值而不是文本

R语言ggplot2可视化:默认情况下ggplot2在x轴和y轴的刻度线和轴之间保留了一些空间设置ggplot2可视化去除可视化结果与坐标轴之间的空间可视化结果与坐标轴紧紧贴合,没有空白区域

ggplot:使用一个共同的 y 轴由多个变量(而不是变量中的多个类别)分隔的多面板/平面散点图

如何网格化面板图,使其在 R 中的 ggplot 中具有分类变量和不同的 x 变量