使用systemfit的R SUR回归导致错误：LU计算奇异：极端条目的比率

Question

我使用R中的systemfit包运行SUR回归。数据集包含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不会引起任何错误并且会给出相同的估计值，但是不会计算协方差矩阵，因此标准误差会不同。