[R-sig-Geo] Variogram does not plot correctly when using layout()
Edzer J. Pebesma
e.pebesma at geo.uu.nl
Tue Oct 24 23:34:03 CEST 2006
Paul, please note that at least 5 R packages calculate variograms, and
not everyone can see from your email address which one you are using
plot is a generic method; to find which instance is called try
so the plot returns is actually a trellis object; the default action on
returned objects is to print it, using the print.trellis function. Read
to find out how to arrange multiple trellis objects on a single page.
Trellis objects don't obey layout().
More likely you don't want separate plots but rather a single trellis
object; you may want to work along these lines:
g = gstat(formula = log(zinc)~1, data = meuse)
g = gstat(g, formula = log(zinc)~1, data = meuse) # define it twice
v = variogram(log(zinc)~1,meuse)
vgm1 = fit.variogram(v,vgm(1, "Sph", 900,1))
vgm2 = fit.variogram(v,vgm(1, "Exp", 300,1))
Paul Hiemstra wrote:
> Dear All,
> I want to make a plot that compares several different models fitted to the
> same experimental variogram.
> # Try different models
> rad_data_mod$sph = fit.variogram(rad_data_vgm,vgm(psill=60,model="Sph",
> rad_data_mod$exp = fit.variogram(rad_data_vgm,vgm(psill=60,model="Exp",
> rad_data_mod$gau = fit.variogram(rad_data_vgm,vgm(psill=60,model="Gau",
> Then I use layout to divide the plotting area into 4 different areas:
> layout(matrix(c(1,2,3,4),2,2,byrow=TRUE), TRUE)
> Now I want to fill the plot area by calling the plot() command three times.
> I expected that the sub-areas would be filled with the three plots. But the
> variogram plot fills the entire plot area and not just the small sub-area.
> A test function:
> x = seq(1:1000)
> y = x^2
> performed as expected producing one plot with four subplots.
> Does anybody know how to make such a trellis plot for variogram plots?
> kind regards,
> Paul Hiemstra
> [[alternative HTML version deleted]]
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
More information about the R-sig-Geo