R语言 使用gstat软件包的变差函数计算点对的规则是什么?

3yhwsihp  于 12个月前  发布在  其他
关注(0)|答案(1)|浏览(159)

在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是这个估计的点对数量。我想知道它的计算规则。
先谢谢你了!

vulvrdjw

vulvrdjw1#

np是在相应的距离间隔中具有距离的对的数量;使用widthcutoff来操作这些间隔。

相关问题