R Package quantreg：提取p值

Question

我有一個大約250年的最大降雨量測量數據系列，maxima [，]，並且想要一次對所有系列應用分位數迴歸，並獲得每個迴歸模型的顯着性R.

library(quantreg) 


    qmag <- array(NA, c(250,4)) 
    taus <- c(0.05, 0.1, 0.95, 0.975) 

    for(igau in 1:250){ 
    qure <- rq(maxima[,igau+1]~maxima[,1], tau=taus) 
    qmag[igau,] <- coef(qure)[2,] 

    } 

我已經試過

summary(qure, se="boot")$p.value 
    ci(qure) 

和其他類似的變化，但得到NULL值。實際上是否可以自動將quantreg中的p值提取到表中，而不是僅爲每個模型在summary（）中單獨查看它們？

Answer 1

看看在STR（）的概要對象的：

require(quantreg) 
data(engel) 
mod <- rq(foodexp ~ income, data = engel) 
summ <- summary(mod, se = "boot") 
summ 
str(summ) 
summ$coefficients[,4]