[R-br] formatação de gráfico

[Maurício Lordêlo]

Saudações a todos.
No CMR abaixo construo seis gráficos e estou tendo dificuldade para
acrescentar uma legenda única e de posicionar um "título" no centro de cada
dupla de gráficos.
Fiz um desenho manualmente para verem como gostaria que ficasse o resultado
final:
https://www.dropbox.com/s/o7snb3i206bz4dr/graficos.jpg
Consegui ter acesso ao livro "R Graphics - Paul Murrel" e encontrei esta
página que talvez tenha a informação que preciso:
https://www.dropbox.com/s/n96r3h1aw9az41q/pag91_92Rgraphics.jpg
Porém, fiz várias tentativas (nas linhas com "comentário" estão as últimas
mudanças) juntamente com as informações contidas no help da função "par" e
não obtive sucesso.
Grato se puderem me ajudar.
Maurício
par(mfrow=c(3,2)
#     ,plt = c(-10,-10,0,0)
# ,pty="s"
# ,mai = c(-5,1,-7,1)
# ,mar = c(-5, 4, -4, 2)
# ,oma = c(1,0,1,0)
# ,pin = c(-3,-6)
# ,xpd = T
# ,din=c(0.5,0.5)
# ,fin=c(0.5,0.5)
)
#par(mfrow=c(3,2))
#layout(matrix(c(1, 2, 3, 4, 5, 6), byrow=TRUE, ncol=2))
tempo<-c(4,5,6,7,8,9,10,11)
A1A10<-rnorm(8,90,4)
A3A10<-rnorm(8,90,4)
plot(A1A10 ~ tempo,axes = F,xlab = '',ylab ="medida",main =
'Intercepto',cex.main=1,type = 'b',xlim=c(4,11),ylim=c(80,100))
mtext('10 repetições',line = 2.9,cex = .8)
lines(A3A10~tempo, type="o",pch="+",lty=2)
axis(2)
A1A20<-rnorm(8,90,4)
A3A20<-rnorm(8,90,4)
plot(A1A20 ~ tempo,axes = F,xlab = '',ylab ="medida",main =
'Intercepto',cex.main=1,type = 'b',xlim=c(4,11),ylim=c(80,100))
mtext('20 repetições',line = 2.8,cex = .8)
lines(A3A20~tempo,type="o",pch="+",lty=2)
axis(2)
A1B10<-rnorm(8,90,4)
A3B10<-rnorm(8,90,4)
plot(A1B10 ~ tempo,axes = F,xlab = '',ylab ="medida",main = 'Tratamento
B',type = 'b',xlim=c(4,11),ylim=c(80,100))
lines(A3B10~tempo,type="o",pch="+",lty=2)
axis(2)
A1B20<-rnorm(8,90,4)
A3B20<-rnorm(8,90,4)
plot(A1B20 ~ tempo,axes = F,xlab = '',ylab ="medida",main = 'Tratamento
B',type = 'b',xlim=c(4,11),ylim=c(80,100))
lines(A3B20~tempo,type="o",pch="+",lty=2)
axis(2)
A1C10<-rnorm(8,90,4)
A3C10<-rnorm(8,90,4)
plot(A1C10 ~ tempo,axes = F,xlab = "tempo",ylab ="medida",main =
'Tratamento C',type = 'b',xlim=c(4,11),ylim=c(80,100))
lines(A3C10~tempo,type="o",pch="+",lty=2)
axis(1, at = seq(4, 11, 1))
axis(2)
#legend(9.5,60,c("Opção 1",""),lty = c(1,2),ncol = 1,horiz=TRUE)
A1C20<-rnorm(8,90,4)
A3C20<-rnorm(8,90,4)
plot(A1C20 ~ tempo,axes = F,xlab = "tempo",ylab ="medida",main =
'Tratamento C',type = 'b',xlim=c(4,11),ylim=c(80,100))
lines(A3C20~tempo,type="o",pch="+",lty=2)
axis(1, at = seq(4, 11, 1))
axis(2)
#legend(0.5,60,c("","Opção 3"),lty = c(1,2),ncol = 1,horiz=TRUE)
#legend(0,60,c("Opção 1","Opção 3"),lty = c(1,2),ncol = 1,horiz=TRUE)
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140203/4b129453/attachment.html>