分位数回归不产生输出

时间:2014-05-07 15:59:30

标签: r regression quantreg

我正在尝试使用rq()包中的quantreg函数运行一些分位数回归。我已经设置了我的数据并创建了模型等(按照其他网站上的说明),这一切看起来都不错。见下文; REWTPC是连续变量,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中运行它。

0 个答案:

没有答案