如何在R中绘制QQ情节？

Question

任务：

绘制分位数 - 分位数（Q-Q）图以检查伽马分布是否是我的数据的良好模型，而不依赖于qqplot。要在此Q-Q图中判断拟合优度，请绘制Q-Q图，以便根据拟合的Gamma分布生成三组150个观测值。

所以我有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分布式观测值;我进一步假设您知道基础Gamma分布的shape和rate参数，并且您希望将观察到的分位数与Gamma分布中的理论分位数进行比较，并使用相同的shape和{ {1}}参数。

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

rate = 1