Подсчет количества точек вне доверительного интервала (CI) на графике qqp (график квантиля-сравнения) - PullRequest
0 голосов
/ 09 марта 2019

Я хотел бы найти любую функцию для вычисленного количества точек вне доверительного интервала (CI 95%) на графике qqp (график квантиля-сравнения).

В моем примере:

Пакеты

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 plot

Теперь я хотел бы использовать любую функцию для создания вектора с количеством точек вне доверительного интервала (CI) на графике qqp (Quantile-Comparison Plot). Это возможно? Спасибо

1 Ответ

0 голосов
/ 10 марта 2019

qqp () не считает количество точек за пределами доверительного интервала, но вычисляет информацию, необходимую для получения этого количества. Вы можете просто изменить код (car:::qqPlot.default), если измените:

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

вывод показывает количество точек за пределами доверительной границы.

...