Você não explorou o suficiente o argumento "par". Veja se é isto aí! op <- par(mfrow=c(3,1), mai = c(2,1,0,0), oma = c(3,0,0,0), pin = c(5.5,1.6), xpd = T) tempo<-c(5,7,9,11,13) m1A<-rnorm(5,90,4) m2A<-rnorm(5,90,4) m3A<-rnorm(5,90,4) m4A<-rnorm(5,90,4) plot(m1A ~ tempo, axes = F, xlab = '', ylab ="%", main = 'Intercepto', type = 'b', xlim=c(5,13), ylim=c(75,100)) lines(m2A~tempo,col="red",type="b") lines(m3A~tempo,col="green",type="b") lines(m4A~tempo,col="black",type="b") axis(2) mtext('GRÁFICOS', line = 1) m1B<-rnorm(5,90,4) m2B<-rnorm(5,90,4) m3B<-rnorm(5,90,4) m4B<-rnorm(5,90,4) plot(m1B ~ tempo, axes = F, xlab = '', ylab ="%", main = 'Tratamento B', type = 'b', xlim=c(5,13), ylim=c(75,100)) lines(m2B~tempo,col="red",type="b") lines(m3B~tempo,col="green",type="b") lines(m4B~tempo,col="black",type="b") axis(2) legend(13.5, 100, letters[1:4], col = 1:4, lty = rep(1,4), ncol = 1) m1C<-rnorm(5,90,4) m2C<-rnorm(5,90,4) m3C<-rnorm(5,90,4) m4C<-rnorm(5,90,4) plot(m1C ~ tempo, axes = F, xlab = "medidas de tempo", ylab ="%", main = 'Tratamento C', type = 'b', xlim=c(5,13), ylim=c(75,100)) lines(m2C~tempo,col="red",type="b") lines(m3C~tempo,col="green",type="b") lines(m4C~tempo,col="black",type="b") axis(1, at = seq(5, 13, 2)) axis(2) mtext('medidas de tempo', line = -15, cex = .8) (s,f,p) Allaman