我想生成一个图形,有一个基本图形和ggplot图形的组合.以下代码显示了我使用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)")
这就产生了
这些小组中的大多数看起来足以让我将其纳入我的报告中.然而,显示自相关性的曲线图需要改进.通过使用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()
然而,作为GGTRAP不是基础图形,我们不能结合GGPrice与布局或PAR(MFROW).如何用ggplot生成的自相关图替换从基本图形生成的自相关图?我知道我可以使用网格.如果我所有的图形都是用ggplot绘制的,那我该怎么做呢?如果ggplot中只生成了一个图形呢?