Как объединить два базовых графика в R?

Вы можете использовать стандартные методы базы R, например. par(mfrow=c(2, 1)). Просто не забудьте установить graphics.reset=F для wt.image(), иначе при завершении будут сброшены исходные графические параметры, что, в частности, повторно применяет настройку mfrow, и следующий график будет помещен на новую страницу.

par(mfrow=c(2, 1))
wt.image(my.wt, main = "wavelet power spectrum",
         legend.params = list(lab = "wavelet power levels (quantiles)", 
                              lab.line = 3.5, 
                              label.digits = 2),
         periodlab = "period (days)",
         graphics.reset=F)
wt.image(my.wt, main = "wavelet power spectrum",
         legend.params = list(lab = "wavelet power levels (quantiles)", 
                              lab.line = 3.5, 
                              label.digits = 2),
         periodlab = "period (days)",
         graphics.reset=F)

Очевидно, вы не можете размещать "recordedplot" объекты на layout.

Не знаю, зачем вам вообще использовать grDevices::recordPlot, возможно, вы хотите имитировать какое-то поведение ggplot. Однако база plot или WaveletComp::wt.image, которая ее использует, всегда завершает построение графика, поэтому вы всегда будете видеть свои графики дважды.

Таким образом, вы можете придерживаться layout или использовать par(mfrow=), как в ответе @robert-hacken, и установить graphics.reset=FALSE согласно его совету.

Я рекомендую использовать устройство png, а не окно предварительного просмотра RStudio; это намного быстрее и создает жестко запрограммированный размер изображения. График хранится в рабочем каталоге (или используйте png('/<path>/foo.R')). Если вы не хотите его использовать, просто закомментируйте строки png и dev.off:

png('foo.png', 800, 600)  ## open png device

layout(matrix(c(1, 2), nrow=2))
par(mar=c(3, 3, 3, 3) + .6)  ## slightly reduce margins
wt.image(my.wt, main = "wavelet power spectrum",
         legend.params=list(lab = "wavelet power levels (quantiles)", 
                            lab.line=3.5, 
                            label.digits=2),
         periodlab = "period (days)",
         graphics.reset=FALSE)
wt.image(my.wt, main = "wavelet power spectrum",
         legend.params=list(lab = "wavelet power levels (quantiles)", 
                            lab.line=3.5, 
                            label.digits=2),
         periodlab = "period (days)",
         graphics.reset=FALSE)

dev.off()  ## close png device

Примечание. Не используйте layout вместе с опциями par(mfrow=, mfcol=), поскольку они несовместимы.

Обновлять

Чтобы объединить несколько анализов в графике, сначала сохраните анализы в list, а затем прокрутите список для построения графика. Пример:

my.data3 <- my.data2 <- my.data1 <- my.data

my.wt_lst <- Map(
  analyze.wavelet,
  my.data=list(my.data1, my.data2, my.data3),
  my.serie='x',
  loess.span=c(0, .25, .5),
  dt=1/24,
  dj=1/20,
  lowerPeriod=1/4,
  make.pval=TRUE,
  n.sim=10
)

Map перебирает элементы синхронно, т.е. если есть только один, как в dt=1/24, он всегда будет использовать 1/24, если их несколько, как в my.data=list(my.data1, my.data2, my.data3), или loess.span=c(0, .25, .5), он будет использовать элементы один за другим. Обратите внимание, что список или data.frame сам по себе представляет собой список с несколькими элементами, и вам нужно будет перечислить его снова, например. my.data=list(my.data1) вместо просто my.data=my.data1.

png('foo.png', 1200, 600)
layout(matrix(seq_along(my.wt_lst), nrow=length(my.wt_lst)))
par(mar=c(3, 3, 3, 3) + .6)  ## slightly reduce margins
for (i in seq_along(my.wt_lst)) {
  wt.image(my.wt_lst[[i]], main = "wavelet power spectrum",
           legend.params=list(lab = "wavelet power levels (quantiles)", 
                              lab.line=3.5, 
                              label.digits=2),
           periodlab = "period (days)",
           graphics.reset=FALSE)
}
dev.off()