пространственная автокорреляция для данных временных рядов

11

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

Я уверен, что в данных имеется пространственная автокорреляция, которая обязательно повлияет на регрессионный анализ агрегированных данных. У меня вопрос - как запустить тест SAC для данных временных рядов? Нужно ли мне смотреть на SAC остатков от моей регрессии для каждого года (глобальный Моран I)? Или я могу провести один тест за все годы?

После того, как я проверил, что да, есть SAC, легко ли было это решить? Мой фон статистики минимален, и все, что я прочитал о пространственно-временном моделировании, звучит очень сложно. Я знаю, что в R есть автоковариатная функция, взвешенная по расстоянию - это вообще просто в использовании?

Я действительно совершенно не понимаю, как оценить / добавить SAC для этой проблемы, и буду очень признателен за любые предложения, ссылки или ссылки. Заранее спасибо!

Rozza
источник
1
Вы хотите смоделировать лежащие в основе пространственные процессы, или вы хотите скорректировать оценку матрицы дисперсии-ковариации для учета пространственной автокорреляции?
generic_user
1
Или вы хотите сделать и то, и другое ...
gregmacfarlane
Спасибо за ваши ответы! Я думаю, что мои данные действительно имеют пространственную автозависимость - это биологическое явление, и очень вероятно, что количество видов в соседних единицах повлияет на окружающие единицы. Мои юниты довольно большие, поэтому я намерен просто использовать опцию «соединение по краям и углам», чтобы установить интервал расстояния для проверки SAC. Функции R для моделей SAC выглядят выполнимо (все еще над моей головой!). Еще раз спасибо.
Роза
Добро пожаловать на сайт Rozza, это нужно оставить как комментарий к одному из ответов. Хотя не забудьте опротестовать существующие ответы, если они предоставили полезные ответы, и пометьте один из них как ответ на ваш вопрос, если вы считаете, что он был удовлетворительным.
Энди W

Ответы:

5

Согласно этой статье , OLS является последовательным при наличии пространственной автокорреляции, но стандартные ошибки неверны и должны быть скорректированы. Соломон Сян предоставляет для этого код stata и matlab . К сожалению, я не знаком с любым R-кодом для этого.

Есть, конечно, другие подходы к такого рода проблемам в пространственной статистике, которые явно моделируют пространственные процессы. Этот просто раздувает стандартные ошибки.

К сожалению, теоретические эконометристы получают удовольствие от запутывания. Связанный документ действительно трудно читать. В основном это говорит о том, что нужно выполнить любую регрессию, которую вы хотите, а затем исправить стандартные ошибки, т. Е. Использовать код из Сян. Пространство не входит в него, пока вы не попытаетесь оценить дисперсию вашей оценки. Интуитивно понятно, что если все различия близки друг к другу, вы менее уверены, что ваша оценка - это не просто след от некоторого ненаблюдаемого пространственного шока.

Обратите внимание, что вам нужно указать пропускную способность ядра, на которой, по вашему мнению, может работать пространственный процесс.

Этот ответ является в основном копией / вставкой повторяющего подобного ответа, который я сделал здесь

generic_user
источник
5

Yзнак равноИксβ+U,Uзнак равноρWU+ε

Yзнак равноρWY+Иксβ+εYзнак равноρWY+Иксβ+WИксλ+ε

Пакет spdep для R содержит множество функций, которые вычисляют матрицы пространственных весов, оценивают пространственные регрессии и делают другие вещи. У меня большой опыт работы с lagsarlmфункциями, но в документации к пакету видно, что есть sacsarlmфункция, которая, кажется, больше того, что вы ищете.

WTT-1T+1

WИкс

gregmacfarlane
источник
1
Хороший пост. Я думаю, что различие, которое ОП должен учитывать при выборе между двумя подходами, заключается в том, повлияет ли «результат» в одном полигоне на результат его соседей. Если это так, используйте подход gmacfarlane. Если нет, то я предлагаю проще.
generic_user