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

【问题标题】：Whats the difference between z~1 and z~LON+LAT in plotting variogram in R?在 R 中绘制变异函数时，z~1 和 z~LON+LAT 有什么区别？
【发布时间】：2015-10-19 07:58:52
【问题描述】：

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 返回 v 和 v1 的不同 gamma 值。这两种操作的基本区别是什么？

help(variogram) 告诉我“定义响应向量和（可能的）回归量的公式，在没有回归量的情况下，使用例如 z~1;”。但是，我没有清楚地理解这句话！谁能详细告诉我什么时候应该使用z~1，或者什么时候应该使用z~LON+LAT？

【问题讨论】：

标签： r spatial sp gstat

【解决方案1】：

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

【讨论】：