我有一个来自格子库的图(levelplot),以及来自oce库的棒图,我想将它们叠加在一起。我正在堆积这些图,因为我通常使用基R,我认为这是因为这些都是专门的图。是否可以在同一页面上堆叠下面的两个图,以便时间戳(x轴)相互对齐?
library(oce)
library(lattice)
library(grid)
第一个图(来自OCE的帮助文档的示例)
# Oceanographic example
data(met)
t <- met[["time"]]
u <- met[["u"]]
v <- met[["v"]]
p <- met[["pressure"]]
dev.new(width=15, height=4)
plotSticks(t, 99, u, v, yscale=25)
第二个情节:
AllT<-NULL
for(i in seq(1,length(t))){
myTime<-t[i]
myTime2<-rep(myTime,10)
AllT<-append(AllT,myTime2)
}
myY<-NULL
for(j in seq(1,length(t))){
mySeq<-seq(1,10)
myY<-append(myY,mySeq)
}
Temp<-seq(1,7200)
MyDF<-data.frame(Temp,myY,AllT)
#Plot code
dev.new(width=15, height=6)
p1<-levelplot(Temp ~ AllT * myY,
data = MyDF,ylim=c(10,1),
xlab = "Time", ylab = "y]",
aspect=0.4,
)
p1
trellis.focus("legend", side="right", clipp.off=TRUE, highlight=FALSE)
grid.text('[Temp]', 1.9, .5, hjust=0.5, vjust=1,rot=270)
trellis.unfocus()
理想情况下,最终的绘图将在顶部绘制棒图(如果x轴排成一行,则不需要在贴图上打勾标签),下面的水平图最小
这是我现在能想到的一样好。必须有一种方法来实际修改grobs
,就像我所知道的那样,但是当我尝试时它不喜欢gTrees
。所以这会让你接近。
这些链接可能有一些用途。
https://www.andrewheiss.com/blog/2016/12/08/save-base-graphics-as-pseudo-objects-in-r/
force a regular plot object into a Grob for use in grid.arrange
#libraries
library(gridGraphics)
library(grid)
library(gridExtra)
library(oce)
library(gggplot2)
library(ggplotify)
#plots
data(met)
t <- met[["time"]]
u <- met[["u"]]
v <- met[["v"]]
p <- met[["pressure"]]
grab_grob <- function(){
grid.echo()
grid.grab()
}
plotSticks(t, 99, u, v, yscale=25)
p1 <- grab_grob() #grab the last plot as a grob
AllT <- NULL
for(i in seq(1, length(t))){
myTime <- t[i]
myTime2 <- rep(myTime, 10)
AllT <- append(AllT, myTime2)
}
myY <- NULL
for(j in seq(1, length(t))){
mySeq <- seq(1, 10)
myY <- append(myY, mySeq)
}
Temp <- seq(1, 7200)
MyDF <- data.frame(Temp, myY, AllT)
#Plot code
dev.new(width = 15, height = 6)
p2 <- levelplot(Temp ~ AllT * myY,
data = MyDF, ylim=c(10, 1),
xlab = "Time", ylab = "y]",
aspect = 0.4,
)
p2
trellis.focus("legend", side="right", clipp.off = TRUE, highlight = FALSE)
grid.text('[Temp]', 1.9, .5, hjust = 0.5, vjust = 1, rot = 270)
trellis.unfocus()
p2a <- as.grob(p2)
grid.arrange(p1, p2a, ncol = 1)