删除`qqPlot`中的置信度区间

Question

考虑以下方法来测试来自假设分布的一组观察数据。为了简单起见，

x<-rnorm(200)
qqPlot(x,  distribution="norm")

在这里我使用了来自car包的qqPlot。我得到一个Q-Q图，如下所示：

​

有没有人能告诉我怎样才能去掉这个手环？作为一个后续问题，有没有人能给我一些参考，这个乐队是干什么的？

Answer 1

我建议使用ggplot2，因为ggplot对象更容易定制。不过，数据x必须是data.frame的列。

set.seed(1)
df <- data.frame(x = rnorm(200))

现在您可以使用stat_qq()添加Q-Q图，并使用stat_qq_line()添加这条线。

library(ggplot2)

ggplot(df, aes(sample = x)) + 
  stat_qq() + 
  stat_qq_line()

​

​

如果您对置信区间的计算感兴趣，我推荐ggpubr包中的ggqqplot()或qqplotr包中的stat_qq_band()的源代码。

Answer 2

我们可以使用envelope=FALSE。

x <- rnorm(200)
car::qqPlot(x, distribution="norm", envelope=FALSE)

​

​

“置信度包络基于来自比较分布的独立随机样本的顺序统计量的SEs (参见Fox，2016)”(?car::qqPlot)