分位数回归的r语言代码

来源:互联网 发布:免费区 无敌淘宝王 编辑:程序博客网 时间:2024/05/15 13:40
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)
0 0
原创粉丝点击