[R] Using filled.contour and contour functions together

Xavier Bodin xavier.bodin at univ-savoie.fr
Fri Feb 11 10:56:57 CET 2011


Dear R help contributors,

I'd like to plot ground temperature with time on X-axis and depth on Y-axis
on this datasets ( http://r.789695.n4.nabble.com/file/n3301033/NEdaily.csv
NEdaily.csv ), and to do so I use the following commands:

	library(RSEIS) 

	xNE	<- seq(1, as.numeric(as.Date(max(NEdaily[[1]])) -
as.Date(min(NEdaily[[1]]))), 1)
	yNE <- rev(c(-0.3, -0.5, -0.7, -0.9, -1.1, -1.4, -1.7, -2, -2.5, -3, -4,
-5, -7, -9, -10))
	zNE <-
mirror.matrix(as.matrix(NEdaily[1:(nrow(NEdaily)-1),2:length(NEdaily)]))

	filled.contour(xNE,yNE,zNE
	, col = myPal(20)
	, zlim = c(-20,20)
	, ylab = "Depth [m]", 
	, xlab = paste("Days since ", as.Date(min(NEdaily[[1]]), format
="%d.%m.%Y"))
	)
	contour(xNE,yNE,zNE, lty = 3, add = T)
	contour(xNE,yNE,zNE, nlevels = 1, level = 0, add = T, lwd = 1.5)

I get this graph ( http://r.789695.n4.nabble.com/file/n3301033/NEdaily.png
NEdaily.png ) and don't understand why filled.contour and contour plots are
no set on the same dimensions and why they don't exactly overlay. Does
anyone have an idea and a solution ?!

Thanks in advance,

Xavier
-- 
View this message in context: http://r.789695.n4.nabble.com/Using-filled-contour-and-contour-functions-together-tp3301033p3301033.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list