R 包 quantreg : Extract p-values

Question

我有大约 250 个年度最大降雨量测量值 ma​​xima[,] 的数据系列，并且希望立即对所有系列应用分位数回归并获取 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 值。实际上是否可以自动将 p 值从 quantreg 提取到表中，而不是仅仅在每个模型的 summary() 中单独查看它们？

Answer 1

看看通过运行 summary-object 的 str() 生成的结构:

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