> h = c(0,0,0,1,1,3,12,27,42,63,80,108,129,108,106,87,79,45,42,31,18,8,5,1,1,2,1)
> y = rep(0:26, h)
> n = length(y)
> n
[1] 1000
> mean(y)
[1] 13.05
> var(y)
[1] 11.80931
> var(y) / mean(y) # ポアソンならこれはほぼ1のはず
[1] 0.9049279
> var(y) / mean(y) * (n-1)
[1] 904.023
$y$ がポアソンなら，最後の量は近似的に自由度 $n - 1 = 999$ の $\chi^2$ 分布をするはずである。
> pchisq(904.023, 999)
[1] 0.01458752
であるから，ちょっと小さすぎる。よって，ポアソンよりも分散が小さい分布のようである。
ところで，ググってみたところ，E. L. Frome, Algorithm AS 171: Fisher's Exact Variance Test for the Poisson Distribution というのがあるようだ。Rで実装されていないだろうか？
（追記）シミュレーションによる計算でも同程度である：
> d = replicate(100000, {r=rpois(1000,13.05); var(r)/mean(r)*999})
> mean(d <= 904.023)
[1] 0.01416
Last modified: