具有两个Y轴的分组框图 [英] Grouped boxplot with two Y axis

问题描述

我想使用GGPLOT2为同一类别制作两个具有不同信息的框图，并使用两个Y轴。 如下图所示：

我有两个区域(AM和AR)的类别信息(速度和旋转角度)，我想将这两个运动模式组合在一个根据区域分组的图形中。

到目前为止，我的代码如下所示：

areaveleang <- ggplot(dados2, aes(x = area)) +
        geom_boxplot(aes(y = speed, color="Bruto" ), fill = "darkblue", alpha=0.7) +
  geom_boxplot(aes(y = angulo_modulo), color="#a50026", position = "dodge") +
        scale_y_continuous(name = "speed (km/h)", sec.axis = sec_axis(~ . *6,name=expression("turn angle (graus)"))) +
                           #breaks = seq(0, 175, 25),
                           #limits=c(0, 175)) +
        scale_x_discrete(name = "Area") +
        #ggtitle("Boxplot of mean ozone by month") +
        theme_bw() +
        theme(plot.title = element_text(size = 14, family = "Tahoma", face = "bold"),
              text = element_text(size = 12, family = "Tahoma"),
              axis.title = element_text(face="bold"),
              axis.text.x=element_text(size = 11),
              legend.position = "bottom") +
        scale_fill_brewer(palette = "Accent")

我的部分数据：

structure(list(lon = c(-38.8745, -38.8676, -38.9445, -38.9811, 
-38.9322, -38.9386), lat = c(-14.09437, -14.25664, -14.34518, 
-14.51857, -14.95994, -15.41684), bmode = c(1.396, 1.375, 1.346, 
1.248, 1.074, 1.289), bmode.5 = c(1, 1, 1, 1, 1, 1), ID = structure(c(53L, 
53L, 53L, 53L, 53L, 53L), .Label = c("10946.05", "20162.03", 
"21792.03", "21800.03", "21810.03", "24640.03", "24641.05", "24642.03", 
"27258.05", "27259.03", "27261.03", "27261.05", "27261.07", "33000.05", 
"37231.07", "37234.05", "37288.07", "42521.07", "50682.07", "50687.07", 
"60004.07", "81122.09", "81123.09", "81124.09", "81125.09", "81126.09", 
"84497.1", "87632.12", "87759.08", "87761.08", "87762.08", "87763.08", 
"87764.08", "87765.08", "87766.08", "87768.08", "87768.11", "87769.11", 
"87770.08", "87773.08", "87773.1", "87773.11", "87774.09", "87774.11", 
"87775.08", "87776.08", "87776.11", "87777.1", "87778.08", "87783.09", 
"88724.1", "88727.09", "111868.11", "111871.12", "112702.12", 
"112712.12", "112714.12", "120947.12", "121189.12", "121192.12", 
"121196.12", "1211931.12"), class = "factor"), sex = structure(c(3L, 
3L, 3L, 3L, 3L, 3L), .Label = c("F", "Fc", "M"), class = "factor"), 
    speed = c(0.921691675, 2.990097297, 2.137609076, 3.259383146, 
    8.173674567, 8.413534172), area = structure(c(3L, 3L, 3L, 
    3L, 3L, 3L), .Label = c("AA", "AM", "AR"), class = "factor"), 
    angulo = c(41.027, -43.41, 29.056, 18.241, -7.125, -4.702
    ), angle_rel = c(0.7160619, -0.7576522, 0.5071251, 0.3183724, 
    -0.1243479, -0.08206201), bat.depth = c(-45L, -397L, -32L, 
    -29L, -21L, -3L), dist = c(5.53015005, 17.94058378, 12.82565446, 
    19.55629887, 49.0420474, 50.48120503), angulo_modulo = c(41.027, 
    43.41, 29.056, 18.241, 7.125, 4.702), curva = structure(c(1L, 
    1L, 1L, 1L, 1L, 1L), .Label = c("curva", "reta"), class = "factor")), row.names = 2:7, class = "data.frame")

PS：对不起，我的英语不好，希望你能理解我的问题

推荐答案

不确定我是否正确，但可能这就是您要查找的内容：

1. 若要获得正确的次轴比例，必须对数据应用反变换，即乘以6时，必须将数据中的angulo_modulo除以6。

2. 要避开这两个框图，我建议使用tidyr::pivot_longer，这样您的两个变量就变成了一个变量的两个类别，该变量默认命名为name，而值则为value。这样就可以将value映射到y和name映射到color、alpha和fill上。默认情况下，ggplot2将为每个区域提供两个并排的框图。

library(ggplot2)
library(tidyr)
library(dplyr)

dados2 %>% 
  mutate(angulo_modulo = angulo_modulo / 6) %>% 
  pivot_longer(cols = c(speed, angulo_modulo)) %>% 
  mutate(name = factor(name, levels = c("speed", "angulo_modulo"))) %>% 
  ggplot(aes(x = area)) +
  geom_boxplot(aes(y = value, color = name, alpha = name, fill = name)) +
  scale_y_continuous(name = "speed (km/h)", sec.axis = sec_axis(~ . * 6, name = expression("turn angle (graus)"))) +
  scale_x_discrete(name = "Area") +
  scale_color_manual(values = c(angulo_modulo = "darkblue", speed = "#a50026")) +
  scale_alpha_manual(values = c(0.7, 1)) +
  theme_bw() +
  theme(plot.title = element_text(size = 14, family = "sans", face = "bold"),
        text = element_text(size = 12, family = "sans"),
        axis.title = element_text(face="bold"),
        axis.text.x=element_text(size = 11),
        legend.position = "bottom") +
  scale_fill_brewer(palette = "Accent")

这篇关于具有两个Y轴的分组框图的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持IT屋！