Пример данных:
##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 или смежные R ^ 2!