fanno pannelli con stessi margini quando combinano grafici ggplot e basi

Question

ho generato una figura che combina ggplot e basi grafici:

t <- c(1:(24*14)) 
P <- 24 
A <- 10 
y <- A*sin(2*pi*t/P)+20 
#***************************************************************************** 
par(mfrow = c(2,1)) 
plot(y,type = "l",xlab = "Time (hours)",ylab = "Amplitude") 
aa <- par("mai") 
plot.new() 

require(gridBase) 
vps <- baseViewports() 
pushViewport(vps$figure) 
pushViewport(plotViewport(margins = aa)) ## I use 'aa' to set the margins 
#******************************************************************************* 
require(ggplot2) 
acz <- acf(y, plot = FALSE) 
acd <- data.frame(Lag = acz$lag, ACF = acz$acf) 
p <- ggplot(acd, aes(Lag, ACF)) + geom_area(fill = "grey") + 
    geom_hline(yintercept = c(0.05, -0.05), linetype = "dashed") + 
    theme_bw() 
grid.draw(ggplotGrob(p)) ## draw the figure 

uso il comando plotViewport e impostare le dimensioni del pannello secondo le dimensioni il primo pannello, ottenuto con il par ("mai"). La figura allegata mostra il risultato. Tuttavia, le dimensioni di entrambi i pannelli non corrispondono, vale a dire il secondo pannello sembra essere leggermente più largo del primo. Come posso superare questo senza dover impostare manualmente i margini con

pushViewport(plotViewport(c(4,1.2,0,1.2))) 

Answer 1

Questo dovrebbe dare alcuni suggerimenti:

library(grid) 
library(ggplot2) 
require(gridBase) 

par(mfrow = c(2,1)) 
plot(1:10) 
a <- par("mai") 
plot.new() 
vps <- baseViewports() 
pushViewport(vps$figure) 

p = qplot(1:10, 1:10) + theme_bw() 
g <- ggplotGrob(p) 

lw = unit(a[2], "inch") - sum(g$widths[1:3]) 

g$widths[[2]] <- as.list(lw + g$widths[[2]]) 
g$widths[[4]] <- as.list(unit(1, "npc") - unit(a[2] + a[4], "inch")) 
g$widths[[5]] <- unit(a[4], "inch") 
grid.draw(g) 

# draw a shaded vertical band to test the alignment 
grid.rect(unit(a[2], "inch"), unit(0, "inch"), 
      unit(1,"npc") - unit(a[2] + a[4], "inch"), 
      unit(2,"npc"), 
      gp=gpar(lty=2, fill="red", alpha=0.1), hjust=0, vjust=0) 

upViewport() 

ma, in realtà, perchè non fare tutto quanto in ggplot2?

Answer 2

L'idea principale è di premere 2 finestre di viste di base per ottenere le dimensioni del riquadro di stampa. La soluzione non è generale.

Per prima cosa tracciare la mia trama di base

t <- c(1:(24*14)) 
P <- 24 
A <- 10 
y <- A*sin(2*pi*t/P)+20 
#***************************************************************************** 
par(mfrow = c(2,1)) 
plot(t,y,type = "l",xlab = "Time (hours)",ylab = "Amplitude") 
plot.new() 

In secondo luogo ho le dimensioni del pannello trama. vpp verrà utilizzato solo per le dimensioni di grobs ggplot (simile all'idea di baptiste sopra)

require(gridBase) 
vps <- baseViewports() 
vpp <- pushViewport(vps$figure,vps$plot) ## here I add a new viewport 
vpp <- current.viewport() 
upViewport(2) 

ggplot2 trama SAVEC come una tabella grob:

require(ggplot2) 
p <- ggplot(acd, aes(Lag, ACF)) + geom_area(fill = "grey") + 
    geom_hline(yintercept = c(0.05, -0.05), linetype = "dashed") + 
    theme_bw() 
data <- ggplot_build(p) 
gtable <- ggplot_gtable(data) 

cambio le dimensioni dei grobs. (Qui il motivo per cui la soluzione non è generale)

gtable$heights[[2]] <- vpp$height 
gtable$heights[[4]] <- vpp$height 
gtable$widths[[4]] <- vpp$width 

ho tracciare

grid.draw(gtable)