R中的简单变差函数，理解gstat:: variogram ()和object gstat

Question

我在R中有一个data.frame，它的变量表示位置，其观测值是这些位置中某个变量的度量。我想测量某些位置依赖于距离的衰减，所以变异函数对我的研究特别有用。

我正在尝试使用gstat库，但是我对某些参数有点困惑。据我所知，(经验性)变异函数只需要作为基本数据：

1. 变量的位置
2. 对这些变量的观察

还有其他参数，比如最大距离，方向，.

现在，gstat::variogram()函数首先需要输入类gstat的对象。检查函数gstat()的文档我发现它输出这个类的一个对象，但是这个函数需要一个formula参数，它被描述为：

将因变量定义为自变量的线性模型的公式；假设因变量的名称为z，对于普通的和简单的kriging，使用公式z~1；对于简单的kriging也定义β(见下文)；对于通用kriging，假设z是线性依赖于x和y的，则使用公式z~x+y。

有人能解释一下这个公式是干什么用的吗？

Answer 1

试一试

methods(variogram)

您将看到gstat有几个用于方差函数的方法，其中一个方法需要一个gstat对象作为第一个参数。

给定一个data.frame，最简单的方法是使用公式方法：

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

它指定在data中，z是感兴趣的观察变量，~1指定一个常量平均模型，~x+y指定坐标在data的列x和y中找到。