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]]
_______________________________________________
R-sig-Geo mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
