为qq图添加置信区间？

有没有办法将置信区间添加到qqplot？为qq图添加置信区间？

我有基因表达值的数据集，我一直在使用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 109 0.932 562 0.958 1.87 0.59 114 391 13.5 1.41 208 2.37 166 3.42]

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

来源

2013-10-11 user2846211

所以，你想从理论正态分布减去你的样本分布？听起来你想要做的就是使用'nls'来使你的数据适合正常的dist函数，并从'nls'的输出中获取置信度数据。 –

如果您提供[最小，可重现的数据集]，您更可能收到有用答案（http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example/5963610＃5963610）以及您尝试过的代码。谢谢！ – Henrik

我用一些数据编辑了我的初始文章。是否可以从qqnorm的输出中获取置信度数据？ – user2846211

库car提供了默认情况下增加了一个逐点置信区间，以正常的QQ积函数qqPlot(...)：

library(car) 
qqPlot(pca1$x)

来源

2013-10-11 13:23:48 sieste

为qq图添加置信区间？

回答

相关问题