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