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