R中统一QQ图的同时零带

Question

这个问题是关于在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％分位数是理想的。

我的问题是：

1. 是否有能够高效执行此任务的R包/代码片段？我认为这属于“功能数据分析”领域，但我不太确定。

2. 我们可以从统一的QQ图的同时空界限中做出什么推论？检查那些落在乐队之外的点？如果p值的真实分布是均匀的（即来自零），我期望合理地在乐队之外多少分？

Answer 1

查看?qqPlot中的car package功能。这将为您绘制一个置信区间（您只需告诉它您感兴趣的分布）。

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