使用 rqpd-package 在面板上运行分位数回归时，我可以获得系数的上限/下限吗？
以下是设置（很好，我是 R 的初学者）：
s = as.factor(rep(1:30,rep(63,30)))
s = as.matrix(s)
fit = rqpd(formula = y ~ x | s, panel = panel(lambda = 1, taus=c(0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 0.95), tauw=c(0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 0.95)), data = pdata, control = list(tmpmax = 500))
capture.output(summary(fit), file = "QR.txt")
代码运行良好，输出为我提供了通常的系数、t 值、p 值等。但我希望找出我的 SE 的上限和下限，以获得 95% 的显着性。
使用 quantreg 包中的 rq 我可以通过捕获以下输出来提取上限/下限：
fit = rq(y~z, tau = taeus)
plot.fit = summary(fit)
capture.output(plot.fit, file = "BOUNDS.txt")
但是，这不适用于我目前正在查看的面板数据。各位有什么指点吗？