如何在R中绘制QQ图？

Question

任务：

在不依赖 qqplot 的情况下，绘制分位数-分位数 (Q-Q) 图来检查 gamma 分布是否适合我的数据。要判断此 QQ 图中的拟合优度，请为您的拟合 Gamma 分布生成的三组 150 个观测值绘制 QQ 图。

所以我有 150 个观察值，我将它们建模为 Gamma 分布，并找到了 k 和 lambda 的估计值：

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 与你观察到的伽马分布观察；我进一步假设您知道基础 Gamma 分布的 shape 和 rate 参数，并且您希望将观察到的分位数与具有相同 shape 和 rate 参数的 Gamma 分布的理论分位数进行比较。强>

这是一个示例，我使用 shape = 2 和 rate = 1 参数从 Gamma 分布生成数据 x。

# 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");