r - 如何估计地球数据的变异函数？

我有相关经度和纬度的数据。如何根据点之间的大圆距离获得此数据的变异函数？

这个简单的例子包含赤道上的所有数据:

require(geoR)

long <- seq(-179, 180)
x <- sin(pi * long / 180) + rnorm(length(long))
V <- variog(data=x, coords=cbind(long, 0))
# variog: computing omnidirectional variogram
plot(V)

第一个和最后一个点实际上仅相隔 1 度，但我天真的尝试导致 variog 认为它们相隔 359 度。

最佳答案

您应该使用 Semi-variogram from nmle 。它允许您指定一个距离矩阵，您可以自己轻松计算。