A explicação do script a seguir está no artigo R na Prática – Parte 18.

P18_1.R

def.par=par(no.readonly=TRUE)
#setwd("C:\\geoStats\\figuras\\chap3")
#jpeg("Fig3_1.jpeg",width=5,height=5,
#units="in",res=300)
layout(matrix(1:4,2,2))
layout.show(4)
#dev.off()
par(def.par)

P18_2.R

def.par=par(no.readonly=TRUE)
#setwd("C:\\geoStats\\figuras\\chap3")
#jpeg("Fig3_2.jpeg",width=5,height=5,
#units="in",res=300)
x=matrix(c(1:3,3),2,2)
layout(x); layout.show(3)
#dev.off()
par(def.par)

P18_3.R

#setwd("C:\\geoStats\\figuras\\chap3")
#jpeg("Fig3_3.jpeg",width=9.055,height=4.33,
#units="in",res=300)
plot(NA,NA,xlim=c(0,25),ylim=c(0,12),type="n",
frame=FALSE,xaxt="n",yaxt="n",xlab="",ylab="")
stx=c(1,7,7,1,1); sty=c(5,5,7,7,5)
sgx=c(13,19,19,13,13); sgy=sty
slx=c(9,15,15,9,9); sly=c(1,1,3,3,1)
ggx=c(17,23,23,17,17); ggy=c(sly)
grx=c(7,13,13,7,7); gry=c(9,9,11,11,9)
lx1=c(10,10,4,4,16,16); ly1=c(9,8,8,7,8,7)
lx2=c(16,16,12,12,20,20); ly2=c(5,4,4,3,4,3)
corRGB=c(253,185,83,221,254,98,207,223,2,220,
192,67,223,234,103)
corRGB=corRGB/255
segmentos=function(sx,sy,corRGB,cor){
for (i in 1:4){
  segments(sx[i],sy[i],sx[i+1],sy[i+1])}
r=corRGB[1+(cor-1)*3];g=corRGB[2+(cor-1)*3];
b=corRGB[3+(cor-1)*3]
polygon(sx,sy,col=rgb(r,g,b))
}
setas=function(sx,sy){
segments(sx[1],sy[1],sx[2],sy[2])
arrows(sx[3],sy[3],sx[4],sy[4],length=0.10,angle=20)
arrows(sx[5],sy[5],sx[6],sy[6],length=0.10,angle=20)
segments(sx[3],sy[3],sx[5],sy[5])
}
segmentos(stx,sty,corRGB,2)
segmentos(sgx,sgy,corRGB,3)
segmentos(slx,sly,corRGB,4)
segmentos(ggx,ggy,corRGB,5)
segmentos(grx,gry,corRGB,1)
text(((stx[1]+stx[2])/2),((sty[2]+sty[3])/2),
"tradicional",cex=1.25,adj=0.5)
text(((sgx[1]+sgx[2])/2),((sgy[2]+sgy[3])/2),
"grid",cex=1.25,adj=0.5)
text(((slx[1]+slx[2])/2),((sly[2]+sly[3])/2),
"lattice",cex=1.25,adj=0.5)
text(((ggx[1]+ggx[2])/2),((ggy[2]+ggy[3])/2),
"ggplot2",cex=1.25,adj=0.5)
text(((grx[1]+grx[2])/2),((gry[2]+gry[3])/2),
"Gr�fico no R",cex=1.25,adj=0.5)
setas(lx1,ly1); setas(lx2,ly2)
#dev.off()