样本数据:
##Data preperation
set.seed(5)
d <- data.frame(x=seq(-5, 5, len=51))
d$y <- 50 - 0.3*d$x^2 + rnorm(nrow(d))
功能:
Qr=function(x,y,taus) {
QUA5<-quantreg::rq(y ~ 1 + x + I(x^2)+I(x^3)+I(x^4)+I(x^5), tau=0.5, data=d)
return(summary(QUA5,se="boot"))
}
如果可能的话,我不仅要参数符号,还要模型符号和R ^ 2或adj R ^ 2!