Dear,
I would like to create a multipanel plot of anisotropic variogram and 
relative model.
So far I have managed to create either 1. individual plot showing 
experiment and model, or 2. multipanel directional plots where all the 
panels show the correct experiment but the same model.

Say "vm" is the results of variogram function (six horizontal 
directions) and "auto_Sph" is the result of fitting.

The function: plot(vm, group.id=TRUE, auto.key=TRUE) returns multipanel 
plots of experimental data, and that is fine.

Now I'd like to overlay to each panel the respective fitted model. I 
tried with

for (i in 1:6){  plot(vm,model=auto_Sph[[i]]$var_model, ...,  multipanel=TRUE, auto.key=TRUE)},

but in this case the last model curve (auto_Sph[[6]]$var_model) is 
overlaid to all panels.
I then tried with:

vms <- split(vm,vm$dir.hor)

for (i in 1:6){  plot(vms[[i]],model=auto_Sph[[i]]$var_model, ...,  multipanel=TRUE, auto.key=TRUE)},

but still got the same results.

I can create each plot separately and then combine them, but I'd like to 
know if there a more efficient way.

Thanks.

	[[alternative HTML version deleted]]

