我使用mgcv中的gam函数在R中拟合模型,看起来有点像这样。
File
然后我绘制变异函数。我希望有一个与x作为空间坐标的相关性,以及组内的相关性,所以我做
m <- gamm(y ~ s(x))
然而,变异函数有一个名为“n.pairs”的列,它由非常大的数字组成。
这些数字代表什么?
例如,对于距离等于5的第i行,我目前看到的是20000 n.pairs。
然而,没有20,000对观察结果1)具有5的精确差异(相对于x)和2)来自同一组。
那么,如果不是那样,那个数字还代表什么?
不太明白这个: https://stat.ethz.ch/R-manual/R-devel/library/nlme/html/Variogram.lme.html
答案 0 :(得分:0)
对数是指距离等级中的对数,而不是距离为5
的对的数量。
在大多数空间数据集中,任何两对样本都不会以完全相同的距离分开。因此,距离被分组为组,并且是在Group
组之上求和的样本对的数量,它们在彼此的距离D_i
内。 D_i
是当前的距离等级,可能是(0 - 100)m。
nlme 并不估算Group
内的单独变异图,但它只计算Group
级别内的样本之间的距离,但如上所述,然后,它会将Group
。