当我在R中运行分位数回归时,使用quantreg
包,然后运行summary(quantregObject)
,我收到以下错误消息:
base :: backsolve中的错误(r,x,k = k,upper.tri = upper.tri,transpose = transpose,:奇异矩阵在' backsolve'。对角线中的第一个零[1]
有任何建议我如何解决这个问题?
答案 0 :(得分:1)
简而言之,尝试:
summary(quantregObject, se = "iid")
这对您的残差有很强的假设。或者,如果您需要准确性,请使用引导带来获取标准错误:
summary(quantregObject, se = "boot")
如果对由quantreg:rq返回的对象调用摘要
summary(quantregObject)
这将称为summary.rq
。
来自?summary.rq
。
您可以看到有4个选项可以计算标准误差(se)。
根据样本大小(N <1000),默认值为se = "rank"
或se = "nid"
。
"nid
”听起来有些复杂,然后可能会产生backsolve
无法处理的对角矩阵。