Ggplot2: добавление вторичной трансформированной оси x поверх графика

Я работаю с масштабированными данными x, и вам нужно добавить немасштабированную ось x в начало графика для упрощения интерпретации. Я нашел подход для добавления вторичной оси y в Как я могу поместить преобразованный масштаб в правую сторону ggplot2?. Однако я не могу заставить его работать правильно для оси x. Я уверен, что не понимаю какую-то часть кода, но я не могу понять, что это такое. Я пробовал искать в файлах справки ggplot2, а также книгу Wickham ggplot2: Элегантная графика для анализа данных, но если кто-нибудь может указать мне на какую-то соответствующую документацию, я бы очень признателен!

Я работаю с данными о температуре, но я буду использовать данные о озере по ссылке выше, поскольку для этого был написан код. Вот код orignal из этой ссылки:

Чтобы проверить метод добавления оси x вместо оси y, я переключил оси x и y и изменил axis-l на axis-b, чтобы дать:

Это создает новую ось x, но есть пара проблем: 1) находится в середине участка 2) он не создает новую x-метку для "Elevation (ft)"

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

Кроме того, как упоминалось выше, я работаю с масштабированными данными температуры, поэтому в идеале тики не будут выравниваться по верхней и нижней осям, как в этом примере. Есть ли способ сделать это в ggplot2? Произвольный пример из Интернета:

Ответы

Ответ 1

Корень вашей проблемы состоит в том, что вы изменяете столбцы, а не строки.

Настройка с масштабированными метками на оси X второго графика:

## 'base' plot
p1 <- ggplot(data=LakeLevels) + geom_line(aes(x=Elevation,y=Day)) + 
    scale_x_continuous(name="Elevation (m)",limits=c(75,125))

## plot with "transformed" axis
p2<-ggplot(data=LakeLevels)+geom_line(aes(x=Elevation, y=Day))+
    scale_x_continuous(name="Elevation (ft)", limits=c(75,125),
                   breaks=c(90,101,120),
                   labels=round(c(90,101,120)*3.24084) ## labels convert to feet
                   )

## extract gtable
g1 <- ggplot_gtable(ggplot_build(p1))
g2 <- ggplot_gtable(ggplot_build(p2))

## overlap the panel of the 2nd plot on that of the 1st plot
pp <- c(subset(g1$layout, name=="panel", se=t:r))

g <- gtable_add_grob(g1, g2$grobs[[which(g2$layout$name=="panel")]], pp$t, pp$l, pp$b, 
               pp$l)

ИЗМЕНИТЬ, чтобы линии сетки выровнялись с тиками нижней оси, замените указанную выше строку следующим образом: g <- gtable_add_grob(g1, g1$grobs[[which(g1$layout$name=="panel")]], pp$t, pp$l, pp$b, pp$l)

## steal axis from second plot and modify
ia <- which(g2$layout$name == "axis-b")
ga <- g2$grobs[[ia]]
ax <- ga$children[[2]]

Теперь вам нужно убедиться, что вы изменяете правильное измерение. Поскольку новая ось горизонтальна (строка, а не столбец), whatever_grob$heights - это вектор для изменения, чтобы изменить количество вертикального пространства в данной строке. Если вы хотите добавить новое пространство, обязательно добавьте строку, а не столбец (т.е. Используйте gtable_add_rows()).

Если вы сами изменяете сами гниль (в этом случае мы изменяем вертикальное обоснование тиков), обязательно измените y (вертикальное положение), а не x (горизонтальное положение).

## switch position of ticks and labels
ax$heights <- rev(ax$heights)
ax$grobs <- rev(ax$grobs)
ax$grobs[[2]]$y <- ax$grobs[[2]]$y - unit(1, "npc") + unit(0.15, "cm")

## modify existing row to be tall enough for axis
g$heights[[2]] <- g$heights[g2$layout[ia,]$t]

## add new axis
g <- gtable_add_grob(g, ax, 2, 4, 2, 4)

## add new row for upper axis label
g <- gtable_add_rows(g, g2$heights[1], 1)
g <- gtable_add_grob(g, g2$grob[[6]], 2, 4, 2, 4)

# draw it
grid.draw(g)

Замечу, что gtable_show_layout() - очень и очень удобная функция для выяснения того, что происходит.

Ответ 2

Вы можете использовать sec_axis или dup_axis.

library(ggplot2)

ggplot(mtcars, aes(x=disp, y=mpg, color=as.factor(cyl))) +
  geom_point() +
  labs(x="displacement (cubic inches)", col="# of cylinders") +
  scale_x_continuous(sec.axis=sec_axis(trans=~ . * 0.0163871, name="displacement (L)"))

ggplot2 версия 3.1.1