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
______________________________________________ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
