在R中,使用gstat包的变差函数计算点对的规则是什么?计算半变差函数时,随机生成30个点进行计算,每次点对的总和都不一样。代码如下:
library(gstat)
n <- 40
X <- runif(n, 0, 100)
Y <- runif(n, 0, 100)
ref_exponential <- exp(runif(n, 0, 1))
data_exponential <- data.frame(X, Y, ref = ref_exponential)
coordinates(data_exponential) <- ~X + Y
variogram_model_exponential <- variogram(ref ~ 1, data_exponential)
print(sum(variogram_model_exponential$np))
字符串
我试着通过阅读txt文件中给出的数据进行计算,当有30个点时,共有435个点对,符合随机抽样(30*29/2=435)
的规则,当计算40个点时,共有528个点对(40*39/2=780)
,不符合随机抽样的规则。
同时,我参考了gstat variogram的文档,其中解释了np是这个估计的点对数量。我想知道它的计算规则。
先谢谢你了!
1条答案
按热度按时间vulvrdjw1#
np
是在相应的距离间隔中具有距离的对的数量;使用width
和cutoff
来操作这些间隔。