plot - Combine base and ggplot graphics in R figure window -

-

i generate figure has combination of base , ggplot graphics. following code shows figure using base plotting functions of r:

t <- c(1:(24*14))  p <- 24  <- 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)")

which generates enter image description here

most of these panels sufficient me include in report. however, plot showing autocorrelation needs improved. looks better using 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()

enter image description here

however, seeing ggplot not base graphic, cannot combine ggplot layout or par(mfrow). how replace autocorrelation plot generated base graphics 1 generated ggplot? know can use grid.arrange if of figures made ggplot how do if 1 of plots generated in ggplot?

using gridbase package, can adding 2 lines. think if want funny plot grid need understand , master viewports. basic object of grid package. 

vps <- baseviewports() pushviewport(vps$figure) ##   in space of autocorrelation plot

the baseviewports() function returns list of 3 grid viewports. use here figure viewport viewport corresponding figure region of current plot.

here how looks final solution:

enter image description here

library(gridbase) 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") ## last 1 current plot plot.new()              ## suggested @josh vps <- baseviewports() pushviewport(vps$figure) ##   in space of autocorrelation plot vp1 <-plotviewport(c(1.8,1,0,1)) ## create new vp margins, play 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")+   ## setting in title near other plots   theme(plot.title = element_text(size = rel(1.4),face ='bold')) print(p,vp = vp1)        ## suggested @bpatiste

Comments

Post a Comment

Popular posts from this blog

ruby - Trying to change last to "x"s to 23 -

-
i have been playing around irb , have string looks this: irb(main):072:0> puts ["[\"4354 5432 5432 xxxxx\", \"6547 6547 8543 xxxxx\", \"2344 6543 6674 xxxxx\", \"2346 6236 7543 xxxxx\", \"1273 5585 5587 xxxxx\"]"] => nil irb(main):073:0> what want gsub last 2 "x" s of each 17 digit combination, 23 . what i've tried far (i've been using rubular): a.gsub(/\d[x]/,"23") <= grabs "x" s i'm close a.gsub(/\w[x{2}]/,"23") <= grabs "x" s grabs two digits each combination. is there easier way i'm not understanding? a.gsub(/\d{4} \d{4} \d{4} xxx\kxx/, '23') \d{4} = 4 digits what \k drop matched far. can use don't replace entire number, still able validate have in desired format.
Read more

jquery - Clone last and append item to closest class -

-
this bare-bones code on page: <div class="list-container"> <div class="single-list-item"> <button>remove</button> </div> <div class="single-list-item"> <button>remove</button> </div> <div class="single-list-item"> <button>remove</button> </div> </div> <p><button class="add-list-item">add</button></p> <div class="list-container"> <div class="single-list-item"> <button>remove</button> </div> <div class="single-list-item"> <button>remove</button> </div> <div class="single-list-item"> <button>remove</button> </div> </div> <p><button class="add-list-item">add</button></p> <di
Read more

c - Unrecognised emulation mode: elf_i386 on MinGW32 -

-
i'm trying make kernel, , cannot link c output assembly. ld . i'm getting error: unrecognized emulation mode: elf_i386 i'm using windows 10 professional mingw32 , msys. code using: link.ld /* * link.ld */ output_format(elf32-i386) entry(start) sections { . = 0x100000; .text : { *(.text) } .data : { *(.data) } .bss : { *(.bss) } } kernel.c /* * kernel.c */ void kmain(void) { const char *str = "my first kernel"; char *vidptr = (char*)0xb8000; //video mem begins here. unsigned int = 0; unsigned int j = 0; /* loops clears screen * there 25 lines each of 80 columns; each element takes 2 bytes */ while(j < 80 * 25 * 2) { /* blank character */ vidptr[j] = ' '; /* attribute-byte - light grey on black screen */ vidptr[j+1] = 0x07; j = j + 2; } j = 0; /* loop writes string video memory */ while(str[j] != '\0') { /
Read more