R Ошибка в chol.default(A) krigeST из пакета gstat

Question

Я работаю с почасовым набором данных о температуре воздуха, зарегистрированных на ~200 станциях на относительно небольшой территории. Я выбрал пространственно-временную вариограмму (например, суммарную метрику), чтобы соответствовать моим данным, и сейчас пытаюсь делать прогнозы по тем же станциям, чтобы заполнить пробелы в NA (отсутствующее значение). При использовании функции krigeST() для ежедневных агрегированных данных все кажется гладким, но когда я использую его с исходным часовым разрешением, я всегда получаю следующую ошибку:

Error in chol.default(A)
the leading minor of order 68 is not positive definite

Я прогуглил это и обнаружил, что это связано с тем, что матрица не является полностью положительно определенной. Тем не менее, я не уверен, почему это происходит, и мне было интересно, если кто-нибудь из вас знает способ исправить это (обходной путь, чтобы избежать этого).

Answer 1

Есть несколько возможностей, которые приводят к сингулярной ковариационной матрице. Два общих:

1. дубликаты наблюдений (идентичное местоположение и отметка времени),
2. модель вариограммы, которая недостаточно различает наблюдения, что приводит к почти идеально коррелированным наблюдениям.