我有一个大约250个年度最大降雨量测量的数据系列,最大值[,]并且想要一次性对所有系列应用分位数回归并获得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值.实际上是否可以自动从quantreg中将p值提取到表中,而不是仅在每个模型的summary()中单独查看它们？