4

私は約 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)

および他の同様のバリエーションですが、NULL 値を取得します。summary()モデルごとに個別に表示するだけでなく、quantreg からテーブルに p 値を自動的に抽出することは実際に可能ですか?

4

1 に答える 1

10

str()オブジェクトの実行によって生成された構造を見てsummaryください。

require(quantreg)
data(engel)
mod <- rq(foodexp ~ income, data = engel)
summ <- summary(mod, se = "boot")
summ
str(summ)
summ$coefficients[,4]
于 2011-10-27T11:16:08.877 に答える