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]