7

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

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

有什么建议我该如何解决这个问题?

4

1 回答 1

8

简而言之,尝试:

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

于 2018-12-20T10:31:12.387 回答