R 包 quantreg:提取 p 值
我有一个大约 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 值。实际上是否可以自动将 p 值从 quantreg 提取到表中,而不是仅仅在每个模型的 summary()
中单独查看它们?
I have a data series of around 250 annual maximum rainfall measurements, maxima[,] and want to apply quantile regression to all series at once and obtain the significance of each regression model in 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,]
}
I've tried
summary(qure, se="boot")$p.value
ci(qure)
and other similar variations but get NULL values. Is it actually possible to automatically extract the p-values from quantreg to a table, rather than just viewing them individually in summary()
for each model?
如果你对这篇内容有疑问,欢迎到本站社区发帖提问 参与讨论,获取更多帮助,或者扫码二维码加入 Web 技术交流群。
绑定邮箱获取回复消息
由于您还没有绑定你的真实邮箱,如果其他用户或者作者回复了您的评论,将不能在第一时间通知您!
发布评论
评论(1)
看一下通过运行
summary
对象的str()
生成的结构:have a look at the structure produced by running
str()
of thesummary
-object: