How to overlay contourplot of a dataset A and a levelplot of a dataset B?

Dear all,

Just to inform that I have solved my problem in a very elegant way, thanks
to the layer approach given by package latticeExtra.

After loading my grids like before:

x<-1:10
y<-1:10
grid1<-expand.grid(x=x,y=y)
grid2<-expand.grid(x=x,y=y)
z1<-grid1$x^2 + grid1$y^2
z2<-2*grid2$x^2 - grid2$y^2
grid1$z<-z1
grid2$z<-z2

I just need to do:

levelplot(z ~ x*y,grid2) + as.layer(contourplot(z ~ x*y,grid1));

The result is a beautiful plot with a contourplot over a levelplot. There
are some examples of overlay of trellis plots in the manual of
latticeExtra: http://latticeextra.r-forge.r-project.org/

Thank you for your attention,

Best,

Charles


On Sat, Aug 2, 2014 at 4:01 PM, Charles Novaes de Santana <
charles.santana at gmail.com> wrote:

Dear all,

Does anyone know a way to overlay a contourplot and a levelplot of
different datasets, both datasets with the same dimension?

Let's say I have 2 10x10 grids, like those below:

library(lattice)

x<-1:10
y<-1:10
grid1<-expand.grid(x=x,y=y)
grid2<-expand.grid(x=x,y=y)
z1<-grid1$x^2 + grid1$y^2
z2<-2*grid2$x^2 - grid2$y^2
grid1$z<-z1
grid2$z<-z2

I would like to plot z1 and z2 in the same plot: z1 as a contourline and
z2 as a levelplot. I tried to do this in two ways, without success:

plot.new();
contourplot(z ~ x*y,grid1)
par(new=T)
levelplot(z ~ x*y,grid2)

and

levelplot(z ~ x*y,grid1,region=FALSE,contour=TRUE)
par(new=T)
levelplot(z ~ x*y,grid2,region=TRUE,contour=FALSE)

Any clue?

 Thank you very much for your time and any help!

Charles

--
Um ax?! :)

--
Charles Novaes de Santana, PhD
http://www.imedea.uib-csic.es/~charles

--
Um ax?! :)

--
Charles Novaes de Santana, PhD
http://www.imedea.uib-csic.es/~charles

        [[alternative HTML version deleted]]