将置信区间添加到 qq 图？答案

【问题标题】：Adding confidence intervals to a qq plot?将置信区间添加到 qq 图？
【发布时间】：2013-10-19 11:02:22
【问题描述】：

有没有办法在 qqplot 中添加置信区间？

我有一个基因表达值数据集，我使用 PCA 对其进行了可视化：

pca1 = prcomp(data, scale. = TRUE)

我现在通过检查数据的分布与正态分布来寻找异常值：

qqnorm(pca1$x,pch = 20, col = c(rep("red", 73), rep("blue", 33)))

qqline(pca1$x)

这是我的数据：

数据= [2.48 104 4.25 219 0.682 0.302 1.09 0.586 90.7 344 13.8 1.17 305 2.8 79.7 3.18 3.18 109 0.932 562 0.958 1.87 0.59 114 391 13.5 1.41 208 2.37 166 3.42]

我现在想绘制 95% 的置信区间来检查哪些数据点位于外部。关于如何做到这一点的任何提示？

【问题讨论】：

所以你想从理论正态分布中减去你的样本分布？听起来您想要做的是使用 nls 将您的数据拟合到正常的 dist 函数并从 nls 的输出中获取置信度数据。
如果您提供minimal, reproducible data set 以及您尝试过的代码，您更有可能收到有用的答案。谢谢！
我已经用一些数据编辑了我的初始帖子。 qqnorm的输出是否可以抓取置信度数据？
请阅读我发布的链接并相应地格式化您的示例数据。谢谢。
尝试将数据剪切并粘贴到 R 命令行中。不工作。使用dput 并分配给data <-

标签： r confidence-interval

【解决方案1】：

库car 提供了qqPlot(...) 函数，它默认为普通qq-plot 添加一个逐点置信包络：

library(car)
qqPlot(pca1$x)

【讨论】：