我正在尝试使用包rq()
中的函数运行一些分位数回归quantreg
。我已经设置了我的数据并创建了模型等(按照另一个站点上的说明),一切看起来都很好。见下文; 是连续变量,RE
定义如下。WTPC
Armored
# Convert Treatment to 1/0 Binary factor where Armored=1
Armored<-as.numeric(Treatment=="A")
# Define variables
Y<-WTPC
X<-cbind(RE,Armored)
# Quantile Regression @ 25th
QR25<-rq(Y~X,tau=0.25)
summary(QR25)
但是,当我调用模型时summary()
,rq()
我得到了这个
> summary(QR25)
Call: rq(formula = Y ~ X, tau = 0.25)
tau: [1] 0.25
Coefficients:
coefficients lower bd upper bd
(Intercept) 11.09456 9.34170 18.22058
XRE -1.81530 -3.46350 -0.78062
XArmored -3.68480 -14.29227 -0.74389
这是很好的信息,但我想要标准误差、p 值等,就像正常的线性回归一样。我检查了文档和其他示例,看起来这是执行此操作的方法。QR25
建立summary(QR25)
一个模型(
这是我第一次尝试分位数回归,所以我可能错过了前面的步骤或需要指定其他内容。有什么想法吗?我在 RStudio 版本 0.98.501 中运行它。