我正在尝试使用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)
,这应该会给出带有错误和p值的回归输出。
这是我对分位数回归的第一次尝试,所以也许我错过了之前的步骤或需要指定别的东西。有什么想法吗?我在RStudio版本0.98.501中运行它。