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", 
range=40e3,nugget=5))
rad_data_mod$exp = fit.variogram(rad_data_vgm,vgm(psill=60,model="Exp", 
range=40e3,nugget=5))
rad_data_mod$gau = fit.variogram(rad_data_vgm,vgm(psill=60,model="Gau", 
range=40e3,nugget=5))

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.

plot(rad_data_vgm,rad_data_mod$sph)
plot(rad_data_vgm,rad_data_mod$exp)
plot(rad_data_vgm,rad_data_mod$gau)

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
plot(x,y)
plot(x,y)
plot(x,y)

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]]


