У меня есть data.frame
в R, чьи переменные представляют местоположения, а чьи наблюдения являются мерами определенной переменной в этих местоположениях. Я хочу измерить ослабление зависимости для определенных мест в зависимости от расстояния, поэтому вариограмма особенно полезна для моих исследований.
Я пытаюсь использовать библиотеку gstat
, но меня немного смущают некоторые параметры. Насколько я понимаю (эмпирическая) вариограмма должна быть нужна только в качестве исходных данных:
- Расположение переменных
- Наблюдения за этими переменными
А затем другие параметры, такие как максимальное расстояние, направления, ...
Теперь функция gstat::variogram()
требует в качестве первого входа объект класса gstat. Проверяя документацию функции gstat()
, я вижу, что она выводит объект этого класса, но эта функция требует аргумента formula
, который описывается как:
формула, определяющая зависимую переменную как линейную модель независимых переменных; предположим, что зависимая переменная имеет имя z, для обычного и простого кригинга используется формула z~1; для простого кригинга также определите бета (см. ниже); для универсального кригинга предположим, что z линейно зависит от x и y, используйте формулу z~x+y
Может ли кто-нибудь объяснить мне, для чего эта формула?