我想生成一個包含base和ggplot圖形組合的圖形。下面的代碼顯示了使用基礎繪製的R功能我的身影:在R圖窗口中合併base和ggplot圖形
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)")
產生
大多數面板看起來夠我在我的報表中。但是,顯示自相關的情節需要改進。這看起來使用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()
然而,看到ggplot不是將基本圖形,我們不能與佈局或參數(mfrow)結合ggplot。我怎樣才能將由基礎圖形生成的自相關圖替換爲由ggplot生成的自相關圖?我知道我可以使用grid.arrange,如果我的所有數字都是用ggplot製作的,但是如果只有其中一個圖形在ggplot中生成,我該如何做到這一點?
它可能是幾乎一樣簡單,而且看起來更一致,使用'polygon'與ACF的'輸出()'來構建類似於'ggplot'基本圖形情節一。 –
感謝您的回覆。這個問題的目的在於學習如何在數字窗口中結合ggplot和基礎圖形,我意識到可能有更有效的方式來生成所示的圖形,但爲了將來的目的,我想了解指定的方法。 – KatyB
查看'gridBase'軟件包... –