[ Изменить апрель 2016: решение в этом потоке больше не отображает добавленную ось правильно - новый поток по этой проблеме был открыт в ggplot2 2.1.0 сломал мой код? Вторичная трансформированная ось теперь отображается неправильно]
Я работаю с масштабированными данными x, и вам нужно добавить немасштабированную ось x в начало графика для упрощения интерпретации. Я нашел подход для добавления вторичной оси y в Как я могу поместить преобразованный масштаб в правую сторону ggplot2?. Однако я не могу заставить его работать правильно для оси x. Я уверен, что не понимаю какую-то часть кода, но я не могу понять, что это такое. Я пробовал искать в файлах справки ggplot2, а также книгу Wickham ggplot2: Элегантная графика для анализа данных, но если кто-нибудь может указать мне на какую-то соответствующую документацию, я бы очень признателен!
Я работаю с данными о температуре, но я буду использовать данные о озере по ссылке выше, поскольку для этого был написан код. Вот код orignal из этой ссылки:
library(ggplot2)
library(gtable)
library(grid)
LakeLevels<-data.frame(Day=c(1:365),Elevation=sin(seq(0,2*pi,2*pi/364))*10+100)
p1 <- ggplot(data=LakeLevels) + geom_line(aes(x=Day,y=Elevation)) +
scale_y_continuous(name="Elevation (m)",limits=c(75,125))
p2<-ggplot(data=LakeLevels)+geom_line(aes(x=Day, y=Elevation))+
scale_y_continuous(name="Elevation (ft)", limits=c(75,125),
breaks=c(80,90,100,110,120),
labels=c("262", "295", "328", "361", "394"))
#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)
ia <- which(g2$layout$name == "axis-l")
ga <- g2$grobs[[ia]]
ax <- ga$children[[2]]
ax$widths <- rev(ax$widths)
ax$grobs <- rev(ax$grobs)
ax$grobs[[1]]$x <- ax$grobs[[1]]$x - unit(1, "npc") + unit(0.15, "cm")
g <- gtable_add_cols(g, g2$widths[g2$layout[ia, ]$l], length(g$widths) - 1)
g <- gtable_add_grob(g, ax, pp$t, length(g$widths) - 1, pp$b)
# draw it
grid.draw(g)
Чтобы проверить метод добавления оси x вместо оси y, я переключил оси x и y и изменил axis-l на axis-b, чтобы дать:
library(ggplot2)
library(gtable)
library(grid)
LakeLevels<-data.frame(Day=c(1:365),Elevation=sin(seq(0,2*pi,2*pi/364))*10+100)
p1 <- ggplot(data=LakeLevels) + geom_line(aes(x=Elevation,y=Day)) +
scale_x_continuous(name="Elevation (m)",limits=c(75,125))
p2<-ggplot(data=LakeLevels)+geom_line(aes(x=Elevation, y=Day))+
scale_x_continuous(name="Elevation (ft)", limits=c(75,125),
breaks=c(80,90,100,110,120),
labels=c("262", "295", "328", "361", "394"))
#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)
ia <- which(g2$layout$name == "axis-b")
ga <- g2$grobs[[ia]]
ax <- ga$children[[2]]
ax$widths <- rev(ax$widths)
ax$grobs <- rev(ax$grobs)
ax$grobs[[1]]$x <- ax$grobs[[1]]$x - unit(1, "npc") + unit(0.15, "cm")
g <- gtable_add_cols(g, g2$widths[g2$layout[ia, ]$l], length(g$widths) - 1)
g <- gtable_add_grob(g, ax, pp$t, length(g$widths) - 1, pp$b)
# draw it
grid.draw(g)
Это создает новую ось x, но есть пара проблем: 1) находится в середине участка 2) он не создает новую x-метку для "Elevation (ft)"
Мне нужна ось, чтобы она отображалась в верхней части графика, и мне нужна соответствующая метка оси. Может ли кто-нибудь сказать мне, что я делаю неправильно?
Кроме того, как упоминалось выше, я работаю с масштабированными данными температуры, поэтому в идеале тики не будут выравниваться по верхней и нижней осям, как в этом примере. Есть ли способ сделать это в ggplot2? Произвольный пример из Интернета:
