1

这个问题是关于基于 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，我们可以进行 10,000 次模拟和如果需要 90% 的同时界限，请获取下 5% 和上 5% 的分位数。

我的问题是：

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

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

4

1 回答 1

2

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

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

在此处输入图像描述

于 2013-10-15T22:38:08.847 回答