2011-10-27 17 views
3

私は約250の年間最大降水量のデータ系列、maxima [、]を持っており、すべてのシリーズに一括して分数回帰を適用し、各回帰モデルの意義を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値を取得しました。実際には、モデルごとにsummary()で個別に表示するのではなく、quantregからテーブルに自動的にp値を抽出することは可能ですか?

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

答えて

9

は()要約対象のSTRを見て
関連する問題