R Package Quantreg: экстракт р-значения

Question

У меня есть ряд данных около 250 ежегодных измерений максимального количества осадков, максимум [,] , и я хочу применить квантильную регрессию ко всем рядам одновременно и получить значение каждой модели регрессии в 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

взгляните на структуру, созданную при запуске str() объекта summary:

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