R: Помощь, необходимая для понимания RMtrend (RandomFields)
При работе с пакетом R RandomFields Вы указываете свою модель, описывающую пространственное распределение ваших данных, добавляя различные компоненты.
Например, точки в ваших данных ближе друг к другу, чем точки, поэтому они следуют экспоненциальной вариограмме. Но все еще есть некоторый шум, в результате чего ваша вариограмма имеет эффект самородка. Вы бы написали:
RMexp(var = partial sill, scale = range) + RMnugget(var = nugget)
добавляя самородок отдельно.
Представьте, что ваши значения данных дополнительно увеличиваются как линейная функция от ваших x-координат. Я уверен, что вы все равно добавите + RMtrend() к приведенной выше модели, указав тренд.
И это сейчас, где я застреваю. Я действительно не понимаю help(RMtrend)-страница - особенно примеры оставляют меня в затруднении (что это R.p-функцию делать? Это страница справки для меня совершенно бессмысленно, не зная C-язык).
Тем не мение, RMtrend принимает только один аргумент: mean, который не может быть формулой, только числовой вектор или другой RMmodel, Мне кажется, как будто не было RMmodel определить линейные или полиномиальные тренды, чтобы я мог передать их RMtrend, Я тоже пытался обмануть RMtrend чтобы дать ему формулу через
RMtrend(RFformula(response Variable ~ x))
который дает ошибку позже, говоря, что RFformula был неправильно использован. И кажется, что невозможно дать RMtrend() вектор с длиной числа узлов сетки сетки моделирования, определяющий "среднее" (= смоделированный тренд) для каждого из них.
Например, моя сетка должна иметь 62250 очков. И если я симулирую без тренда, я действительно получаю матрицу с 249 строками и 250 столбцами, но если я передаю вектор с 62250 числовыми записями RMtrend() RFsimulate() все равно выдаст мне ошибку (или приведет к фатальной ошибке R):
Error in rfInit(model = list("Simulate", setseed = eval(parse(text = "quote(set.seed(seed=seed))")), :
'RFsimulate' : Number of data (62250) not a multiple of the number of locations (62567 x 0)
Так как это можно дать RMtrend линейный или полиномиальный тренд?