R包quantreg:提取p-值
R包quantreg:提取p-值
提问于 2011-10-27 09:45:53
我有一个大约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)以及其他类似的变体，但得到的是空值。实际上，是否可以将p值从quantreg自动提取到表中，而不只是在summary()中单独查看每个模型的p值呢？
回答 1
Stack Overflow用户
回答已采纳
发布于 2011-10-27 11:16:08
查看运行str() of summary-object生成的结构：
require(quantreg)
data(engel)
mod <- rq(foodexp ~ income, data = engel)
summ <- summary(mod, se = "boot")
summ
str(summ)
summ$coefficients[,4]页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/7914192
复制相关文章
相似问题
腾讯云开发者
