时间戳

首页 > 解决方案 > 在 qqp（分位数比较图）图中计算置信区间 (CI) 之外的点数

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

问题描述

我想在 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图

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

标签： rperformanceplotdistributiondata-fitting

解决方案

qqp() 不计算置信范围之外的点数，但它计算获得此计数所需的信息。如果您更改，您可以简单地修改代码 ( car:::qqPlot.default)：

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

输出显示置信范围之外的点数。

推荐阅读