R Package quantreg：提取p值

Question

我有一个大约250个年度最大降雨量测量值的数据系列，最大值[，] ，并希望同时对所有序列应用分位数回归，并获得R中每个回归模型的显着性。 / p>

    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]