quantreg分位数回归：summary()与stargazer()计算标准误的差异疑问

我使用quantreg::rq()估计分位数回归时，调用summary()返回如下错误：

Error in base::backsolve(r, x, k = k, upper.tri = upper.tri, transpose = transpose, : singular matrix in 'backsolve'. First zero in diagonal

查阅资料后，我找到解决方案：将标准误计算方法设为iid，即执行summary(df, se="iid")，此方法有效。但当我使用stargazer()并设置stargazer(df, rq.se="iid")时，却再次触发相同错误：

Error in base::backsolve(r, x, k = k, upper.tri = upper.tri, transpose = transpose, :
singular matrix in 'backsolve'. First zero in diagonal

因此我想请教：summary()与stargazer()在计算标准误时是否存在差异？

以下是复现代码：

library(quantreg)
library(stargazer)
df <- rq(logincome ~ ., data=rfmale1995,tau=c(0.2))
summary(df,se="iid")
stargazer(df,rq.se="iid")

编辑：我查看了stargazer的内部代码，发现rq.se参数本质是直接调用summary(model, se="")提取结果，理论上两者的标准误计算逻辑应该完全一致，这让我更加困惑。

内容的提问来源于stack exchange，提问作者Fang