Простая вариограмма в R, понимание gstat::variogram() и объекта gstat

Question

У меня есть data.frame в R, чьи переменные представляют местоположения, а чьи наблюдения являются показателями определенной переменной в этих местоположениях. Я хочу измерить затухание зависимости для определенных мест в зависимости от расстояния, поэтому вариограмма особенно полезна для моих занятий.

Я пытаюсь использовать gstat библиотека, но я немного смущен определенными параметрами. Насколько я понимаю (эмпирическая) вариограмма должна быть нужна только в качестве исходных данных:

1. Расположение переменных
2. Наблюдения за этими переменными

А затем другие параметры, такие как максимальное расстояние, направления, ...

Сейчас, gstat::variogram() Функция требует в качестве первого ввода объекта класса gstat. Проверка документации функции gstat() Я вижу, что он выводит объект этого класса, но эта функция требует formula аргумент, который описывается как:

формула, определяющая зависимую переменную как линейную модель независимых переменных; предположим, что зависимая переменная имеет имя z, для обычного и простого кригинга используется формула z~1; для простого кригинга также определите бета-версию (см. ниже); для универсального кригинга предположим, что 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 указать, что координаты находятся в столбцах x а также y из data,