R中均匀QQ图的同时零带答案

【问题标题】：simultaneous null band for uniform QQ plot in RR中均匀QQ图的同时零带
【发布时间】：2013-10-15 22:17:13
【问题描述】：

这个问题是关于基于 R 中的逐点分位数估计绘制同时空带的便捷方法。假设我获得了一个 p 值向量，并且我想通过使用均匀 QQ 图来查看它的分布是否接近均匀分布（......如果不是，哪些点“脱颖而出”，即落在乐队之外）。这是一个示例（其中u_obs 表示不一定来自均匀分布的 p 值——仅用于说明目的）：

library(gap)
u_obs = qqrunif(1000)
qqunif(u_obs,pch=21,bg="blue",bty="n")

根据第 5 页讨论的非参数方法，本文的图 1（左下图）：http://stat.wharton.upenn.edu/~buja/PAPERS/paper-sim.pdf，如果需要 90% 的同时界限，我们可以进行 10,000 次模拟并获得下限和上限 5% 分位数.

我的问题是：

是否有能够有效执行此任务的 R 包/代码 sn-ps？我认为这属于“功能数据分析”领域，但我不太确定。
我们可以从均匀 QQ 图的同时空界得出什么推论？检查那些超出范围的点？如果 p 值的真实分布是均匀的（即来自空值），我期望有多少点合理地在带外？

【问题讨论】：

我不确定它是否是同时的，但你应该看看car::qqPlot()和mgcv::qq.gam()

标签： r simultaneous quantile

【解决方案1】：

看看car package 中的?qqPlot 函数。这将为您绘制一个置信带（您只需告诉它您对哪个分布感兴趣）。

library(car)
set.seed(9)
u_obs <- runif(100)
windows()
  qqPlot(u_obs, distribution="unif")

【讨论】：

感谢您的评论！是的，它非常接近我想要的。当我查看qqPlot 的帮助手册，特别是qqPlot.lm 时，我注意到simulate 参数如果为真，则通过参数引导计算置信包络，但它仅适用于lm 对象。该方法归功于 Atkinson (1985)。我希望使用参数引导程序获得qqPlot.default 工作以获得（同时）置信带。有什么建议？非常感谢！
理论线不经过[0,0]是最烦人的。
我想你可以尝试从头开始编写代码，@alittleboy，但我不知道有什么函数可以做到这一点。我试过mod <- lm(u_obs~1)，但qqPlot.lm 将接受的唯一分布是“norm”和“t”。 OTOH，我相信，绘制的信心信封确实来自理论，但似乎将这条线作为参考。根据 DWin 的评论，使用line="robust" 选项可能会更好。
@gung: 是的，line=robust 可以让参考线通过(0,0)。这是一个更基本的问题：我们如何解释生成的 QQ 图？预计 5% 会超出范围？顺便说一句，这支乐队是“同时”乐队吗？作者确实出于某种原因叠加了乐队，但我不确定如何使用它。谢谢！