Thank you very much!
I'm using  par(new = TRUE)  because in my real case, the 2 plots have 
different ylim (different y-scale).
I got what I wanted by using the same xlim in the barplot and in the 
plot, as suggested by Peter.

My real case:

par(mar= c(7, 4, 5, 5) + 0.1)
area   <- c(136, 3426, 5594, 29268, 19080, 31461, 72629, 225443)/ 10000
nfires <- c(48, 134, 80, 131, 27, 24, 23, 14)

t <- barplot(area, xlim=c(0.5,9.5), ylim=c(0, 25), ylab="Area burnt (x 
1000 ha)", xlab="Fire size class")
box()
par(new= T, yaxs = "i", xpd=TRUE)
plot(t, nfires, type="o", xlim=c(0.5,9.5), pch=19, ylim=c(0, 140), 
axes=F, xlab="", ylab="")
axis(4); mtext("Number of fires", side=4, line=3, col=1)
legend(2, 155, "Number", lty= 1, pch= 19, bty="n", cex=0.8)
legend(2, 163, "    Area", fill= 2, bty="n", cex=0.8)


Cheers

Juli


Martin Maechler wrote:

>>>>>>"juli" == juli g pausas <juli@ceam.es>
>>>>>>    on Thu, 18 Dec 2003 10:54:08 +0100 writes:
>>>>>>            
>>>>>>
>
>    juli> Dear colleges,
>    juli> I'm trying to combine a barplot and a plot in a single figure as follows:
>
>    juli> data <- 1:6
>    juli> t <- barplot(data, axes=F)
>    juli> par(new= T)
>    juli> plot(t, data, type="b")
>
>    juli> However, as you can see in the example, the dots of
>    juli> the second plot do not fall in the midpoint of the
>    juli> bars in the first. Any trick for setting the 2 plots
>    juli> at the same scale?
>
>yes, use
>
>   bd <- barplot(data)
>   points(bd, data, type = "b")
>
>instead.  
>A general recommendation: 
>   Try to *not* use par(new = TRUE) if you can.
>
>Martin
>
>  
>

-- 
Juli G. Pausas
Centro de Estudios Ambientales del Mediterraneo (CEAM)
C/ Charles R. Darwin 14, Parc Tecnologic,
46980 Paterna, Valencia, SPAIN
Tel: (+ 34) 96 131 8227; Fax: (+ 34) 96 131 8190
mailto:juli@ceam.es
http://www.gva.es/ceam

GCTE Fire Network - http://www.gva.es/ceam/FireNetwork


	[[alternative HTML version deleted]]


