Построение нескольких переменных X для одной непрерывной переменной Y в одном графике

Question

Я пытаюсь нарисовать линейный график с двумя x переменными в x-axis с одной непрерывной y переменной в y-axis.Количество x1 и x2 различно.df выглядит следующим образом -

df <- structure(list(val = c(3817,2428,6160,6729,7151,7451,6272,7146,7063,6344,5465,6169,7315,6888,7167,6759,4903,6461,7010,7018,6920,3644,6541,31862,31186,28090,28488,29349,28284,25815,23529,20097,19945,22118), type = c("1wt", "1wt", "3wt", "3wt", "3wt", "5wt", "5wt", "7wt", "7wt", "7wt","10wt","10wt","10wt","15wt","15wt","20wt","20wt","25wt","25wt","25wt","30wt","30wt","30wt","20m","20m","15m","15m","15m","10m","10m","5m", "5m", "5m", "5m"), group = c("A", "A", "A", "A", "A", "A", "A", "A", "A", "A", "A", "A", "A", "A", "A", "A", "A", "A", "A", "A", "A", "A", "A", "B", "B", "B", "B", "B", "B", "B", "B", "B", "B", "B")), row.names = c(NA, 34L), class = "data.frame")

, где x переменные -

x1 <- factor(df$type, levels = c("1wt", "3wt", "5wt", "7wt", "10wt", "15wt", "20wt", "25wt", "30wt")) и

x2 <- factor(df$type, levels = c("20m", "15m","10m","5m"))

Я хочу иметь отдельные строки для x1 и x2 с различными цветами и легендами в соответствии с df$group на оси x с df$val наось Y.не могли бы вы помочь мне сделать это?Заранее спасибо.

Answer 1

РЕДАКТИРОВАТЬ: добавлено ниже

Вот подход, который предполагает, что цель состоит в том, чтобы сопоставить диапазон возможных значений типа из группы A с диапазоном возможных значений из группы B.

Маркировка может быть добавлена ​​вручную, но я не думаю, что есть какой-либо простой способ использовать две категориальные оси x вместе на одном графике.

df2 <- df %>%
  mutate(x = case_when(type == "1wt" ~ 0,
                       type == "3wt" ~ 1,
                       type == "5wt" ~ 2,
                       type == "7wt" ~ 3,
                       type == "10wt" ~ 4,
                       type == "15wt" ~ 5,
                       type == "20wt" ~ 6,
                       type == "25wt" ~ 7,
                       type == "30wt" ~ 8,

                       type == "20m"  ~ 0/3 * 8,
                       type == "15m"  ~ 1/3 * 8,
                       type == "10m"  ~ 2/3 * 8,
                       type == "5m"   ~ 3/3 * 8))


ggplot(df2, aes(x, val, color = group, group = group)) + 
  geom_point() +
  geom_smooth(method = lm)

2-й подход

Похоже, что ОП хотел бы использовать значения type численно в некотором роде.Если они неразрывно связаны друг с другом описанным способом, я подозреваю, что вводить их в заблуждение так, как если бы они были введены в заблуждение.(См. здесь для обсуждения того, почему это проблема.)

Тем не менее, вот как вы можете это сделать.Во-первых, вот подход, который просто использует числовую часть type как есть.Обратите внимание, что «m», связанный с группой B, находится внизу, а «wt» вверху, связанный с группой A, как в примере, добавленном в комментарии OP ниже.Я добавил цвета к осям, чтобы прояснить это.Это немного нелогично визуально, так как точки, связанные с верхней осью, находятся снизу, и наоборот.

df2 <- df %>%
  # First, let's take the number used in "type" without adjustment
  mutate(x_unadj = parse_number(type))

ggplot(df2, aes(x_unadj, val, color = group, group = group)) + 
  geom_point() +
  geom_smooth(method = lm) + # Feel free to use other smoothing method, but
                             # not obvious to me what would be improvement.
  scale_x_continuous("m", sec.axis = sec_axis(~., name = "wt")) +
  theme(axis.text.x.bottom  = element_text(color = "#00BFC4"),
        axis.title.x.bottom = element_text(color = "#00BFC4"),
        axis.text.x.top     = element_text(color = "#F8766D"),
        axis.title.x.top    = element_text(color = "#F8766D"))

Если этонеудовлетворительно, мы могли бы изменить порядок обеих осей, используя

scale_x_reverse("m", sec.axis = sec_axis(~., name = "wt")) +

Используя ggplot 3.1.0 (с октября 2018 г.), я не смогзаставить вторичную ось x сместиться в направлении, противоположном основной оси. Этот пример из 2017 года, похоже, больше не работает с этой версией.По состоянию на декабрь 2018 года предлагается предложенное исправление , которое рассматривается для решения этой проблемы.