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)
扫一扫
1954
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
