[R-sig-Geo] gstat anisotropy : multipanel plotting of experiment and model

[Ashton Shortridge]

Just a thought, but how about using the 'regular' plotting system, and not the 
lattice plotting system?

Then something like:
par(mfrow=c(2,3))
for (i in 6): {
	PLOT EMPIRICAL VARIOGRAM i
	OVERLAY MODEL i
}

Might do the trick? I suggest this because, being efficient (lazy!), I have 
often found it easier to avoiding lattice than to spend time figuring out how 
to use it properly.

Ashton

On Thursday, October 17, 2013 11:47:38 AM efisio solazzo wrote:
> 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]]
> 
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo

-----
Ashton Shortridge
Associate Professor			ashton at msu.edu
Dept of Geography			http://www.msu.edu/~ashton
235 Geography Building		ph (517) 432-3561
Michigan State University		fx (517) 432-1671