R SUR регрессия с использованием systemfit, приводящая к ошибке: LU в вычислительном отношении единственное число: отношение экстремальных входов

Question

Я запускаю регрессию SUR с использованием пакета systemfit в R. Набор данных содержит возвраты для 80 банков, которые регрессируют по фиктивным переменным, равным 1 в определенные даты и 0 в противном случае. При запуске этого я всегда получаю ту же ошибку:

Ошибка в .solve.dgC.lu (as (a, "dgCMatrix"), b = b, tol = tol): LU в вычислительном отношении единственное число: отношение экстремальных элементов в | diag (U) | знак равно 3.703e-20

Я загрузил набор данных в pastebin, поэтому я надеюсь, что вы легко сможете повторить ошибку, запустив этот код:

library("systemfit")
library("plm")
den <- read.table("https://pastebin.com/raw.php?i=WF3vn1G8", sep=";", header=TRUE)
denpanel<-pdata.frame(den, c("id", "t"))
densur<-systemfit(returns ~ Price + Pre + Event + Post, method = "SUR",data = denpanel)

Регресс SUR работает до 78 банков. Когда я добавляю 79-й банк, он больше не работает.

Я использую R версии 3.5.1 (64 бит).

Буду очень признателен за вашу помощь! Это мой первый пост, поэтому, пожалуйста, дайте мне знать, если я что-то забыл.

Answer 1

Это происходит потому, что у нас G = 80 банков, T = 83 наблюдения на банк и K = 5 параметров на банк, что дает 83-5 = 78 степеней свободы на уравнение, что является проблемой, поскольку SUR оценивает остатки GxGковариационная матрица.Между тем, OLS не будет выдавать никаких ошибок и будет давать те же оценки, но не будет вычислять ковариационную матрицу, и, следовательно, стандартные ошибки будут другими.