0

このドキュメントに従って、一連の頻度データを R のポアソン分布に適合させようとしましたが、goodfit多くの警告がスローされます。

>library(vcd)
>freq
[1]  71   9   1   7  27 158   1   6   1  28   1   7  33  10   2  14   3   5  17  73   8   248
[23]   6 129  11  15
>gf<-goodfit(freq,type="poisson",method="MinChisq")
Warning messages:
1: In optimize(chi2, range(count)) :
  NA/Inf replaced by maximum positive value
2: In optimize(chi2, range(count)) :
  NA/Inf replaced by maximum positive value
3: In optimize(chi2, range(count)) :
  NA/Inf replaced by maximum positive value
4: In optimize(chi2, range(count)) :
  NA/Inf replaced by maximum positive value
5: In optimize(chi2, range(count)) :
  NA/Inf replaced by maximum positive value
6: In optimize(chi2, range(count)) :
  NA/Inf replaced by maximum positive value
7: In optimize(chi2, range(count)) :
  NA/Inf replaced by maximum positive value
8: In optimize(chi2, range(count)) :
  NA/Inf replaced by maximum positive value
9: In optimize(chi2, range(count)) :
  NA/Inf replaced by maximum positive value

freqは整数ベクトルです。ここで問題を説明するために、より大きな整数ベクトルから取得されます。より大きい整数ベクトルは、ポアソンに似た分布を持ちます。

の何が問題なのですかfreq、それともパッケージのバグvcdですか?

4

0 に答える 0