r - 使用 quantreg 获取模型的显着性水平（p 值）

Question

set.seed(5)
library(quantreg)
d <- data.frame(x=seq(-5, 5, len=51))

d$y <- 50 - 0.3*d$x^2 + rnorm(nrow(d))

Taus <- c(0.1,0.3,0.5,0.7,0.9)

QR<-rq(y ~ 1 + x + I(x^2)+I(x^3), tau=Taus, data=d)

我需要找到这个模型的估计系数的 p 值