使用rqpd-package在面板上运行分位数回归时,我可以获得上限/下限吗?

时间:2017-03-05 20:38:46

标签: r quantile quantreg

使用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")

但是这不适用于我目前正在查看的paneldata。你们有什么指针吗?

0 个答案:

没有答案