Chciałbym wygenerować figurkę, która ma kombinację grafiki bazowej i graficznej. Poniższy kod pokazuje moją figurę przy użyciu funkcji drukowania bazowego R:

t <- c(1:(24*14)) 
P <- 24 
A <- 10 
y <- A*sin(2*pi*t/P)+20

par(mfrow=c(2,2))
plot(y,type = "l",xlab = "Time (hours)",ylab = "Amplitude",main = "Time series")
acf(y,main = "Autocorrelation",xlab = "Lag (hours)", ylab = "ACF")
spectrum(y,method = "ar",main = "Spectral density function", 
         xlab = "Frequency (cycles per hour)",ylab = "Spectrum")
require(biwavelet)
t1 <- cbind(t, y)
wt.t1=wt(t1)
plot(wt.t1, plot.cb=FALSE, plot.phase=FALSE,main = "Continuous wavelet transform",
     ylab = "Period (hours)",xlab = "Time (hours)")

Który generuje

Większość z tych paneli wydaje mi się wystarczająca do uwzględnienia w moim raporcie. Jednak wykres pokazujący autokorelację wymaga poprawy. Wygląda to znacznie lepiej za pomocą ggplot:

require(ggplot2)
acz <- acf(y, plot=F)
acd <- data.frame(lag=acz$lag, acf=acz$acf)
ggplot(acd, aes(lag, acf)) + geom_area(fill="grey") +
  geom_hline(yintercept=c(0.05, -0.05), linetype="dashed") +
  theme_bw()

Jednak widząc, że ggplot nie jest podstawową grafiką, nie możemy łączyć ggplot z układem lub par (mfrow). W jaki sposób mogę zastąpić wykres autokorelacji wygenerowany z grafiki bazowej wykresem wygenerowanym przez ggplot? Wiem, że mogę użyć grid.arrange, jeśli wszystkie moje dane zostały utworzone za pomocą ggplot, ale jak mam to zrobić, jeśli tylko jeden z wykresów jest generowany w ggplot?