[英]Adding confidence intervals to a qq plot?
有没有办法在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 109 0.932 562 0.958 1.87 0.59 114 391 13.5 1.41 208 2.37 166 3.42]
我现在想绘制95%置信区间来检查哪些数据点在外面。 关于如何做到这一点的任何提示?
该库car
提供了功能qqPlot(...)
增加了一个逐点置信区间在默认情况下正常的QQ-情节:
library(car)
qqPlot(pca1$x)
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.