R包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()中单独查看它们？

- Mari

Answer 1

看一下运行str()函数后，summary对象所生成的结构：

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

- EDi