setwd("I:/研一课程/2.2回归分析/R/data")#设定当前的工作目录,重要!
d<-read.table("faithful.txt",header=T)
library("quantreg")
fit1=rq(waiting~eruptions,tau=0.5,data=d)
fit1
summary(fit1)
r1=resid(fit1) #残差序列
c1=coef(fit1) #模型系数
summary(fit1,se="nid") #系数的假设检验
fit2=rq(waiting~eruptions,tau=0.25,data=d)
fit3=rq(waiting~eruptions,tau=0.75,data=d)
anova(fit1,fit2,fit3)
d<-read.table("faithful.txt",header=T)
library("quantreg")
fit1=rq(waiting~eruptions,tau=0.5,data=d)
fit1
summary(fit1)
r1=resid(fit1) #残差序列
c1=coef(fit1) #模型系数
summary(fit1,se="nid") #系数的假设检验
fit2=rq(waiting~eruptions,tau=0.25,data=d)
fit3=rq(waiting~eruptions,tau=0.75,data=d)
anova(fit1,fit2,fit3)