R SUR回归使用系统拟合导致误差: LU计算单数:极端项的比率

Question

我使用R中的systemfit包运行SUR回归，数据集包含80家银行的回报，这些银行在虚拟变量上回归，在特定日期等于1，否则等于0。在运行这个程序时，我总是获得相同的错误：

.solve.dgC.lu中的错误(as(a，"dgCMatrix")，b= b，tol = tol)：LU计算单数:极端条目的比率(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不会产生任何误差，也不会给出相同的估计，但不会计算协方差矩阵，因此，标准误差将是不同的。