[英]R Package quantreg: Extract p-values
我有一个包含大约 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 值。 实际上是否可以将 quantreg 中的 p 值自动提取到表格中,而不是仅在每个模型的summary()
单独查看它们?
看看结构产生通过运行str()
的的summary
-object:
require(quantreg)
data(engel)
mod <- rq(foodexp ~ income, data = engel)
summ <- summary(mod, se = "boot")
summ
str(summ)
summ$coefficients[,4]
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.