Dear Paul, Thanks for helping. Is there a way to call grid() first? The problem seems to be that everything drawn before grid() is overplotted.
Cheers, Marius require(grid) require(gridBase) pdf(file="Rplot.pdf", width=8, height=8, onefile=FALSE) ## set up the grid layout plot.new() # start (empty) new page with 'graphics' gl <- grid.layout(5, 5, widths=unit(c(1.8, 8, 0.8, 8, 0.8), "cm"), heights=unit(c(0.8, 8, 0.8, 8, 1.5), "cm")) pushViewport(viewport(layout=gl)) ## plot data par. <- par(no.readonly=TRUE) # save plot settings for(i in 1:2) { # rows i. <- if(i > 1) i+2 else i+1 # jumping over gaps for(j in 1:2) { # columns j. <- if(j > 1) j+2 else j+1 # jumping over gaps pushViewport(viewport(layout.pos.row=i., layout.pos.col=j.)) grid.rect(gp=gpar(col=NA, fill="gray90")) # background par(plt=gridPLT()) ## plot par(new=TRUE) # always do this before each new 'graphics' plot plot(1:5, if(i==1 && j==2) c(1, 1, 2, 3, 1) else 1:5, log="y", xlab="", ylab="", frame.plot=FALSE, xaxt="n", yaxt="n") grid(col="white", lty="solid", lwd=1.6, equilogs=FALSE) # background grid upViewport() } } par(par.) dev.off() Paul Murrell <p...@stat.auckland.ac.nz> writes: > Hi > > On 24/09/12 18:06, Marius Hofert wrote: >> Dear Paul, >> >> Thank you for helping. This works great. >> >> I then tried to put in a grid (via grid()). Why does that fail? > > Because grid() is used to add lines to an existing plot; just put the grid() > call AFTER the plot() call and it should work ok. > > Paul > >> Cheers, >> >> Marius >> >> >> require(grid) >> require(gridBase) >> >> pdf(file="Rplot.pdf", width=8, height=8, onefile=FALSE) >> >> ## set up the grid layout >> plot.new() # start (empty) new page with 'graphics' >> gl <- grid.layout(5, 5, widths=unit(c(1.8, 8, 0.8, 8, 0.8), "cm"), >> heights=unit(c(0.8, 8, 0.8, 8, 1.5), "cm")) >> pushViewport(viewport(layout=gl)) >> >> ## plot data >> par. <- par(no.readonly=TRUE) # save plot settings >> for(i in 1:2) { # rows >> i. <- if(i > 1) i+2 else i+1 # jumping over gaps >> for(j in 1:2) { # columns >> j. <- if(j > 1) j+2 else j+1 # jumping over gaps >> pushViewport(viewport(layout.pos.row=i., layout.pos.col=j.)) >> grid.rect(gp=gpar(fill="gray90")) # background >> par(plt=gridPLT()) >> ## plot >> par(new=TRUE) # always do this before each new 'graphics' plot >> grid(col=1) >> plot(1:10, 1:10, log="y", xlab="", ylab="", >> xaxt=if(i==2) "s" else "n", yaxt=if(j==1) "s" else "n") >> upViewport() >> } >> } >> par(par.) >> dev.off() >> >> >> >> Paul Murrell <p...@stat.auckland.ac.nz> writes: >> >>> Hi >>> >>> On 24/09/12 09:36, Marius Hofert wrote: >>>> Hi, >>>> >>>> Why does the upper left panel (in the plot below) not have a gray >>>> background? >>> >>> That is a symptom of the conflict that occurs when the 'graphics' package >>> and >>> the 'grid' package both try to initialise a new page. >>> A good rule of thumb is to start a new page with 'graphics' first and THEN >>> add >>> grid' stuff ('grid' is better at sharing), so a minor adjustment to your >>> code >>> would be (#PAUL marks the changes) ... >>> >>> >>> pdf(file="Rplot.pdf", width=8, height=8, onefile=FALSE) >>> #PAUL >>> # Start (empty) new page with 'graphics' >>> plot.new() >>> ## set up the grid layout >>> gl <- grid.layout(5, 5, widths=unit(c(1.8, 8, 0.8, 8, 0.8), "cm"), >>> heights=unit(c(0.8, 8, 0.8, 8, 1.5), "cm")) >>> pushViewport(viewport(layout=gl)) >>> ## plot data >>> par. <- par(no.readonly=TRUE) # save plot settings >>> for(i in 1:2) { # rows >>> i. <- if(i > 1) i+2 else i+1 # jumping over gaps >>> for(j in 1:2) { # columns >>> j. <- if(j > 1) j+2 else j+1 # jumping over gaps >>> pushViewport(viewport(layout.pos.row=i., layout.pos.col=j.)) >>> grid.rect(gp=gpar(fill="gray90")) # background >>> par(plt=gridPLT()) >>> ## plot >>> #PAUL >>> # ALWAYS do this before each new 'graphics' plot >>> par(new=TRUE) >>> plot(1:10, 1:10, log="y", xlab="", ylab="", >>> xaxt=if(i==2) "s" else "n", yaxt=if(j==1) "s" else "n") >>> upViewport() >>> } >>> } >>> par(par.) >>> dev.off() >>> >>> >>> Hope that helps. >>> >>> Paul >>> >>>> Cheers, >>>> >>>> Marius >>>> >>>> >>>> require(grid) >>>> require(gridBase) >>>> >>>> pdf(file="Rplot.pdf", width=8, height=8, onefile=FALSE) >>>> >>>> ## set up the grid layout >>>> gl <- grid.layout(5, 5, widths=unit(c(1.8, 8, 0.8, 8, 0.8), "cm"), >>>> heights=unit(c(0.8, 8, 0.8, 8, 1.5), "cm")) >>>> if(FALSE) grid.show.layout(gl) >>>> pushViewport(viewport(layout=gl)) >>>> >>>> ## plot data >>>> par. <- par(no.readonly=TRUE) # save plot settings >>>> for(i in 1:2) { # rows >>>> i. <- if(i > 1) i+2 else i+1 # jumping over gaps >>>> for(j in 1:2) { # columns >>>> j. <- if(j > 1) j+2 else j+1 # jumping over gaps >>>> pushViewport(viewport(layout.pos.row=i., layout.pos.col=j.)) >>>> grid.rect(gp=gpar(fill="gray90")) # background >>>> par(plt=gridPLT()) >>>> ## plot >>>> plot(1:10, 1:10, log="y", xlab="", ylab="", >>>> xaxt=if(i==2) "s" else "n", yaxt=if(j==1) "s" else "n") >>>> par(new=TRUE) # to be run after first plot >>>> upViewport() >>>> } >>>> } >>>> par(par.) >>>> dev.off() >>>> >>>> ______________________________________________ >>>> R-help@r-project.org mailing list >>>> https://stat.ethz.ch/mailman/listinfo/r-help >>>> PLEASE do read the posting guide >>>> http://www.R-project.org/posting-guide.html >>>> and provide commented, minimal, self-contained, reproducible code. >>>> >> ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.