Quantum backsolve总结错误

Posted 2021-04-12

当我在R中运行分位数回归时，使用quantreg包，然后运行summary(quantregObject)，我收到以下错误消息：

base :: backsolve中的错误（r，x，k = k，upper.tri = upper.tri，transpose = transpose，：'backsolve'中的奇异矩阵。对角线中的第一个零[1]

任何建议我怎么能解决这个问题？

答案

简而言之，尝试：

summary(quantregObject, se = "iid")

这会对您的残差做出强有力的假设。或者，如果您需要准确性，请使用引导带来获取标准错误：

summary(quantregObject, se = "boot")

如果在quantreg：rq返回的对象上调用summary

summary(quantregObject)

这将调用summary.rq。

来自?summary.rq。

您可以看到有4个选项来计算标准错误（se）。根据样本大小（N <1000），默认值为se = "rank"或se = "nid"。 "nid“做了一些听起来很复杂的事情，然后可能产生backsolve无法处理的对角线矩阵。