Как нарисовать QQ сюжет в R?

Question

Задача:

Чтобы нарисовать график квантиль-квантиль (Q-Q), чтобы проверить, является ли гамма-распределение хорошей моделью для моих данных, не полагаясь на qqplot. Чтобы оценить правильность соответствия на этом графике Q-Q, нарисуйте графики Q-Q для трех наборов из 150 наблюдений, полученных из вашего подобранного гамма-распределения.

Итак, у меня есть 150 наблюдений, которые я смоделировал как гамма-распределение и нашел оценки для k и лямбда:

pp = c(1:150)/151
qq= qgamma(pp,k,rate=lambda)
layout(matrix(c(1,1,2,3), 2, 2, byrow = TRUE))
plot(qq,sort(x),main="Q-Q plot with gamma distribution",xlab='Quantiles of gamma distribution',ylab="Ordered data")

Однако, когда я строю qq индивидуально 3 раза, получается точно такой же график?

Answer 1

Я предполагаю, что у вас есть вектор x с вашими наблюдениями, распределенными по гамма-распределению; Кроме того, я предполагаю, что вам известны параметры shape и rate базового гамма-распределения и что вы хотите сравнить наблюдаемые квантили с теоретическими квантилями из гамма-распределения с такими же параметрами shape и rate.

Вот пример, где я генерирую данные x из гамма-распределения с параметром shape = 2 и rate = 1.

# Set seed for reproducibility
set.seed(2017);

# Generate some Gamma distributed data
x <- rgamma(100, shape = 2, rate = 1);

# Sort x values
x <- sort(x);

# Theoretical distribution
x0 <- qgamma(ppoints(length(x)), shape = 2, rate = 1);

plot(x = x0, y = x, xlab = "Theoretical quantiles", ylab = "Observed quantiles");
abline(a = 0, b = 1, col = "red");