计算 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)

输出显示置信区间外的点数。