2 
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 的不同伽马值。这两种操作的基本区别是什么？

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

4

1 回答 1

1

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

于 2015-08-13T08:54:15.913 回答