在qqp（分位数比较图）图中计算置信区间（CI）之外的点数

Question

我想在qqp（分位数比较图）图中找到置信区间（CI 95％）以外的计算点数的任何函数。

在我的示例中：

包裹

require(MASS)
require(car)

模拟60个泊松值

Resp<-rpois(60,1)

拟合二项式负分布

nbinom <- fitdistr(Resp, "Negative Binomial")

使用qqp绘制

qqp(Resp, "nbinom", size = nbinom$estimate[[1]], mu = nbinom$estimate[[2]])

现在，我想使用任何函数来创建一个向量，该向量的点数超出qqp（分位数比较图）图中的置信区间（CI）。这个有可能？谢谢

Answer 1

qqp（）不会计算置信度包络之外的点数，但会计算获得此计数所需的信息。如果您进行更改，则只需修改代码（car:::qqPlot.default）

：

outerpoints <- sum(ord.x > upper |  ord.x < lower)
print(outerpoints)
if (length(outerpoints) > -1) 
    outerpoints
else invisible(NULL)

输出显示置信度包络之外的点数。