2017-08-20 111 views
1

R中有一個data.frame,它的變量代表位置,其觀測值是這些位置中某個變量的度量值。我想根據距離來測量某些位置的相關衰減,所以變異函數對我的研究特別有用。R中的簡單變異函數,瞭解gstat :: variogram()和對象gstat

我想使用gstat庫,但我對某些參數有些困惑。據我理解(經驗)變差函數應該只需要基本的數據:

  1. 變量的位置
  2. 觀察這些變量

,然後在其它參數,如maximun距離,方向,...

現在,gstat::variogram()功能需要作爲第一個輸入類的對象gstat。檢查的功能gstat()文檔我看到它輸出這個類的對象,但該功能需要一個formula參數,它被描述爲:

式定義從屬變量作爲獨立變量的線性模型;假設因變量具有名稱z,對於普通和簡單的克里格使用公式z〜1;對於簡單的克里格也定義了beta(見下文);對於通用克里格法,假設z線性依賴於x和y,使用公式z〜x + y

有人能解釋我這個公式的用途嗎?

回答

0

嘗試

methods(variogram) 

,你會看到gstat有變差的幾種方法,一個需要gstat對象作爲第一個參數。

給定一個data.frame,最簡單的是使用式方法:

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

其指定在dataz是感興趣的觀察變量,~1指定一個恆定的平均模型,~x+y指定在data的列xy中找到座標。