76

base グラフィックと 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()

ここに画像の説明を入力

ただし、ggplot はベース グラフィックではないため、ggplot をレイアウトまたは par(mfrow) と組み合わせることはできません。ベースグラフィックスから生成された自己相関プロットをggplotによって生成されたものに置き換えるにはどうすればよいですか? すべての図が ggplot で作成されている場合、 grid.arrange を使用できることはわかっていますが、ggplot で生成されたプロットが 1 つだけの場合、どうすればよいですか?

4

4 に答える 4

55

gridBase パッケージを使用すると、2 行追加するだけで実行できます。グリッドで面白いプロットをしたい場合は、ビューポートを理解して習得するだけでよいと思います。これは実際には grid パッケージの基本オブジェクトです。

vps <- baseViewports()
pushViewport(vps$figure) ##   I am in the space of the autocorrelation plot

baseViewports() 関数は、3 つのグリッド ビューポートのリストを返します。ここでは図のビューポートを使用します現在のプロットの図の領域に対応するビューポート。

最終的な解決策は次のとおりです。

ここに画像の説明を入力

library(gridBase)
library(grid)

par(mfrow=c(2, 2))
plot(y,type = "l",xlab = "Time (hours)",ylab = "Amplitude",main = "Time series")
plot(wt.t1, plot.cb=FALSE, plot.phase=FALSE,main = "Continuous wavelet transform",
     ylab = "Period (hours)",xlab = "Time (hours)")
spectrum(y,method = "ar",main = "Spectral density function", 
         xlab = "Frequency (cycles per hour)",ylab = "Spectrum")
## the last one is the current plot
plot.new()              ## suggested by @Josh
vps <- baseViewports()
pushViewport(vps$figure) ##   I am in the space of the autocorrelation plot
vp1 <-plotViewport(c(1.8,1,0,1)) ## create new vp with margins, you play with this values 
require(ggplot2)
acz <- acf(y, plot=F)
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()+labs(title= "Autocorrelation\n")+
  ## some setting in the title to get something near to the other plots
  theme(plot.title = element_text(size = rel(1.4),face ='bold'))
print(p,vp = vp1)        ## suggested by @bpatiste
于 2013-01-02T16:29:27.390 に答える
8

私は gridGraphics パッケージのファンです。何らかの理由で、gridBase に問題がありました。

library(ggplot2)
library(gridGraphics)
data.frame(x = 2:10, y = 12:20) -> dat
plot(dat$x, dat$y)
grid.echo()
grid.grab() -> mapgrob
ggplot(data = dat) + geom_point(aes(x = x, y = y)) 
pushViewport(viewport(x = .8, y = .4, height = .2, width = .2))    
grid.draw(mapgrob)

ここに画像の説明を入力

于 2015-09-28T21:21:38.857 に答える