摘要 quantreg backsolve 中的错误

Question

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

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

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

Answer 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无法处理的对角矩阵。