Hi All,

I'm trying to show two cross variograms and their fitted models side-by-side.  Models and variograms are produced with the gstat function fit.lmc.  I would like to do this without also showing the single variograms.  I'm not very familiar with lattice or gstat, and I'm stuck.  What I would like to do from the following data set would be to produce a plot of only the zn.cu graph.  If someone could demonstrate how to do this, I think I would be able to take it from there.


data(meuse)
coordinates(meuse)=~x+y
data(meuse.grid)
gridded(meuse.grid)=~x+y
meuse.g <- gstat(id = "zn", formula = log(zinc) ~ + 1, data = meuse, nmax = 10)
meuse.g <- gstat(meuse.g, "cu", log(copper) ~ 1, meuse, nmax = 10)
meuse.g <- gstat(meuse.g, model = vgm(1, "Sph", 900, + 1), fill.all = T)
meuse.g
meuse.g.var<- variogram(meuse.g, cutoff = 1000)
plot(meuse.g.var)
names(meuse.g)
meuse.g.fit<-fit.lmc(meuse.g.var, meuse.g)
plot(meuse.g.var, model=meuse.g.fit)


I get close to what I want with this extra line:
plot(meuse.g.var, subset=id=="zn.cu", meuse.g.fit)

but xyplot is still providing space for the other plots.  Not sure how to create just the single zn.cu plot.

Thanks,

Doug Collins
Small Farms Extension Specialist
Washington State University
Center for Sustaining Agriculture and Natural Resources
2606 W Pioneer
Puyallup, WA  98371


	[[alternative HTML version deleted]]

