データへのポアソン分布を適合させる

Question

ヒストグラムをプロットし、ヒストグラムにポアソン分布をフィットさせたいと思います。これを行うために、xとyのヒストグラム座標ベクトルをラムダを推定するpoissfit()関数に渡しました。例えば

、ここで私がやったことです：

expecteddata = cat(2,x,y) 
[l,lci] = poissfit(expecteddata) 

私の出力はそうのようになります。

l = 

    44.3766 0.0130 


lci = 

    42.8887 0.0003 
    45.8645 0.0724 

私は考えをプロットするために興味があるラムダを想定しています0.013（私のヒストグラムは周波数ヒストグラムなので私のラムダはとても小さいと思います）。私が使用してポアソンPDFをプロットします

x = 0:50 
y = poisspdf(x,0.013); 
plot(x,y,'r') 

そして、私は結果のオーバーレイプロットを得る：

はしかし、私はこのフィット分布は少し奇妙に見えると思います。それは非常に "ポアソン"のようには見えません。誰かが私が何かを間違ってやっているのか知っていますか

Answer 1

"私はヒストグラムをプロットしており、ヒストグラムにポアソン分布をフィットさせたいと思います。"

私が理解することは、ポアソン分布を測定データの既存のヒストグラムに合わせる必要があることです。私はあなたがfitdist()機能を使用できると信じています。

たとえば、データがxの場合。

[n,bin]=hist(x,100); 
m=n/trapz(bin,n); 
bar(bin,m,'w'); 
hold on 
pd=fitdist(x,'poisson'); 
y=pdf(pd,bin); 
plot(bin,y,'k'); 
hold off; 

は、ポアソン分布曲線がフィットするヒストグラムを表示します。