経度と緯度が関連付けられたデータがあります。ポイント間の大圏距離に基づいて、このデータのバリオグラムを取得するにはどうすればよいですか?
この単純な例には、赤道上のすべてのデータが含まれています。
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 度離れていると考えてしまいます。