2011-10-27 3 views
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] 
:

답변

9

는() 요약 객체의 STR를 보라