A explicação do script a seguir está no artigo R na Prática – Parte 21.
def.par=par(no.readonly=TRUE)
#setwd("C:\\geoStats\\figuras\\chap3")
#jpeg("Fig3_7.jpeg",width=5,height=5,
#units="in",res=300)
cor=c("red","blue","darkgreen","maroon")
tit=c("FUNÇÃO SENO","FUNÇÃO COSENO",
"FUNÇÃO LOG","FUNÇÃO LOG10")
par(mfrow=c(2,2),col=cor[1],cex=1)
x=seq(-pi,pi,0.1)
plot(x,sin(x),type="h",main=tit[1],
xlab="x em radianos",ylab="Função Seno (x)")
par(col=cor[2])
plot(x,cos(x),main=tit[2],
xlab="x em radianos",ylab="Função Coseno (x)")
x=seq(1:51)
par(col=cor[3])
plot(x,log(x),main=tit[3],
xlab="X",ylab="log(x)")
par(col=cor[4])
plot(x,log10(x),main=tit[4],log="x",
xlab="Escala logaritmica",ylab="log10(x)")
par(def.par)
#dev.off()
#setwd("C:\\geoStats\\figuras\\chap3")
#jpeg("Fig3_5.jpeg",width=5,height=5,
#units="in",res=300)
plot.new()
for (i in 1:5){
for (j in 1:5){
kb=j+(i-1)*5
xp=(j-1)*0.20
yp=(i-1)*0.20
points(xp,yp,pch=kb,cex=1.5)
text((xp+0.025*xp),yp,pos=4,cex=1,
labels=as.character(kb))}
}
#dev.off()
#fonte: https://kktg.net/sgr/otw-portfolio
#/fig-15-7/
#setwd("C:\\geoStats\\figuras\\chap3")
#jpeg("Fig3_6.jpeg",width=7.5,height=3,
#units="in",res=300)
def.par=par(no.readonly=TRUE)
par(mai = c(0, 0, 0, 0))
plot.new()
segments(x0=0.1,y0=seq(.9,.1,-.16),
x1=.6,y1=seq(.9,.1,-.16),
lty = c(1:6),lwd = 2)
text(
x=.65,y=seq(.9,.1,-.16),pos = 4,
cex = 1.25,
label = c("lty = 1 ou 'sólido'",
"lty = 2 ou 'tracejado'",
"lty = 3 ou 'pontilhado'",
"lty = 4 ou 'ponto-traço'",
"lty = 5 ou 'traço longo'",
"lty = 6 ou 'dois traços'"))
#dev.off()
#par(def.par)