我在R中有data.frame
,其变量代表位置,其观察值是这些位置中某个变量的度量。我想根据距离测量某些位置的依赖性衰减,因此变异函数对我的研究特别有用。
我正在尝试使用gstat
库,但我对某些参数感到有点困惑。据我所知,(经验)变差函数应该只需要作为基本数据:
然后是其他参数,例如最大距离,方向,......
现在,gstat::variogram()
函数需要第一个输入gstat类的对象。检查函数gstat()
的文档我看到它输出了这个类的对象,但是这个函数需要一个formula
参数,其描述如下:
公式,将因变量定义为自变量的线性模型;假设因变量具有名称z,对于普通和简单克里金法,使用公式z~1;简单克里金法也定义了beta(见下文);对于通用克里金法,假设z线性依赖于x和y,使用公式z~x + y
有人可以解释一下这个公式的用途吗?
答案 0 :(得分:0)
试
methods(variogram)
你会看到gstat有几种变异函数方法,一种需要gstat对象作为第一个参数。
给定data.frame,最简单的方法是使用公式方法:
variogram(z~1, ~x+y, data)
指定在data
中,z
是感兴趣的观察变量,~1
指定常量均值模型,~x+y
指定在列x
中找到坐标{1}}和y
的{{1}}。