Jacques, provided that X and Y are colocated (i.e., have exactly the same observation locations), you get the cross variogram right; the definition of this cross variogram is however: gamma(h)= E[(X(s)-X(s+h))*(Y(s)-Y(s+h))] also, where you select: cv <- v$gamma[1:14] you may be better off using the more general v$gamma[v$id == "X.Y"] Best regards, -- Edzer
Seemingly Similar Threads
- [R-sig-Geo] plot() and Jpeg() increase font size and resolution
- default value for cutoff in gstat variogram()
- gstat 0.9-12: cokriging cross validation and class name incompatibilities
- gstat 0.9-12: cokriging cross validation and class name incompatibilities
- Are package maintainers responsible for name-mangling class names?