如何创建以泊松作为理论分布的 QQ 图

Question

我需要创建一个 QQ 图来检查我观察到的数据是否符合泊松分布。

这是我的数据框：

df = read.table(text = 'Var1 Freq
 1975   10
 1976   12
 1977    9
 1978   14
 1979   14
 1980   11
 1981    8
 1982    7
 1983   10
 1984    8
 1985   12
 1986    9
 1987   10
 1988    9
 1989   10
 1990    9
 1991   11
 1992   12
 1993    9
 1994   10', header = TRUE)

这df$Freq专栏是我感兴趣的专栏，因为观察结果代表了每年的事件数量。

我知道我必须使用该qqplot函数以及qpois创建理论分位数的函数，但是如何呢？

Answer 1

此外，该fitdistrplus包可以用更少的代码来完成此操作。比较经验和理论密度以及 CDF。

library('fitdistrplus')
plot(fitdist(df$Freq,"pois"))

您可以获取 lambda 等并检查其他发行版。不如该ggplot方法灵活，但适合快速检查。