在 R 中绘制变差函数时，z~1 和 z~LON+LAT 有什么区别？

Question

library(sp)
data(meuse)
coordinates(meuse)<-~x+y
v<-variogram(log(zinc)~1,meuse)
v
         np       dist     gamma dir.hor dir.ver   id
    1   57   79.29244 0.1234479       0       0 var1
    2  299  163.97367 0.2162185       0       0 var1
    3  419  267.36483 0.3027859       0       0 var1
    4  457  372.73542 0.4121448       0       0 var1
    .
    .

v1<-variogram(log(zinc)~x+y,meuse)
v1
            np       dist     gamma dir.hor dir.ver   id
        1   57   79.29244 0.1060834       0       0 var1
        2  299  163.97367 0.1829983       0       0 var1
        3  419  267.36483 0.2264256       0       0 var1
        4  457  372.73542 0.2847192       0       0 var1
        .
        .

从上面的代码和输出我可以发现 log(zinc)~1 和 log(zinc)~x+y returns v 和 v1 的不同伽马值。这两个操作的基本区别是什么？

help(variogram)告诉我"formula defining the response vector and (possible) regressors, in case of absence of regressors, use e.g. z~1;"。但是，我没看清楚这句话！谁能详细告诉我什么时候应该使用 z~1 或者什么时候应该使用 z~LON+LAT?

Answer 1

variogram 的第一个参数定义了平均结构的模型。在常数均值的情况下，模型仅包含一个截距，因此 log(zinc)~1。如果均值在 x 和 y 中建模为线性回归模型，则使用 log(zinc)~x+y。在这种情况下，计算普通最小二乘残差作为半变异函数值的基础，而不是测量的 log(zinc) 值。

在 R 中绘制变差函数时，z~1 和 z~LON+LAT 有什么区别？

Whats the difference between z~1 and z~LON+LAT in plotting variogram in R?

r

spatial

gstat

sp