3
나는 최대 250 개의 연간 최대 강우 측정 데이터 시리즈를 가지고 있습니다. maxima [,] 모든 시리즈에 동시에 quantile regression을 적용하고 각 회귀 모델의 중요성을 R.R 패키지 quantreg : 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 값을 가져했습니다. quantify에서 테이블로 p 값을 추출하는 것이 실제로 가능합니까? 각 모델에 대해 summary()로 개별적으로 보는 것보다?
require(quantreg)
data(engel)
mod <- rq(foodexp ~ income, data = engel)
summ <- summary(mod, se = "boot")
summ
str(summ)
summ$coefficients[,4]
: