Пространственные регрессии с панельными данными в R
У меня есть набор данных панели с несколькими сотнями регионов, ~10 лет и пространственными данными для регионов. Я создал весовую матрицу с spdepпакет (через стандартный способ, а затем, nb2listw). Таким образом, у меня есть матрица с весами для каждого региона (по отношению к другим регионам), но каждый регион представлен только один раз.
Я хотел бы запустить некоторые из пространственных регрессий из spdepпакет (lagsarlm, errorsarlm), но я получаю ошибку:
Error in subset.listw(listw, subset, zero.policy = zero.policy) :
Not yet able to subset general weights lists
а также
Error in lagsarlm(y ~ x1 + x2: Input data and weights have different dimensions
Я предполагаю, что это потому, что матрица весов имеет только одну строку на регион (и тогда можно рассчитать только один год). Есть ли у вас какие-либо предложения, как решить проблему? Мои идеи вращаются вокруг следующего:
- Расширить матрицу пространственного веса ИЛИ
- Скажите spdep, что регионы будут повторяться в том же порядке (но как?)
Ждем ваших предложений.