R中的简单变差函数,了解gstat :: variogram()和对象gstat

时间:2017-08-20 18:50:07

标签: r gstat covariogram

我在R中有data.frame,其变量代表位置,其观察值是这些位置中某个变量的度量。我想根据距离测量某些位置的依赖性衰减,因此变异函数对我的研究特别有用。

我正在尝试使用gstat库,但我对某些参数感到有点困惑。据我所知,(经验)变差函数应该只需要作为基本数据:

  1. 变量的位置
  2. 对这些变量的观察
  3. 然后是其他参数,例如最大距离,方向,......

    现在,gstat::variogram()函数需要第一个输入gstat类对象。检查函数gstat()的文档我看到它输出了这个类的对象,但是这个函数需要一个formula参数,其描述如下:

      

    公式,将因变量定义为自变量的线性模型;假设因变量具有名称z,对于普通和简单克里金法,使用公式z~1;简单克里金法也定义了beta(见下文);对于通用克里金法,假设z线性依赖于x和y,使用公式z~x + y

    有人可以解释一下这个公式的用途吗?

1 个答案:

答案 0 :(得分:0)

methods(variogram)

你会看到gstat有几种变异函数方法,一种需要gstat对象作为第一个参数。

给定data.frame,最简单的方法是使用公式方法:

variogram(z~1, ~x+y, data)

指定在data中,z是感兴趣的观察变量,~1指定常量均值模型,~x+y指定在列x中找到坐标{1}}和y的{​​{1}}。