Блог

У меня есть ряд данных из примерно 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)
 

и другие подобные варианты, но получите нулевые значения. Действительно ли возможно автоматически извлекать p-значения из quantreg в таблицу, а не просто просматривать их по отдельности summary() для каждой модели?

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

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