Ich habe eine Reihe von Beobachtungenobs.

enn ich das Histogramm der Beobachtung aufzeichne, sehe ich, dass sie aus einer Gammaverteilung stammen könnte

[counts,x] = hist(obs,[1:max(obs)]);

Ich möchte es mit Chi-Quadrat-Anpassungsgüte beweisen.

So schätze ich zuerst die Gamma-Parameter

paramEsts = fitdist(obs,'Gamma');  

und benutze chi2gof, um zu sehen, ob die Hypotesis wahr ist (h = 0).

[h,p] = chi2gof(obs,'CDF',paramEsts)

Mein Problem ist, dass ich ein p = NaN bekomme ... wie ist das möglich? Wo ist mein Fehler? Vielen Dan

Um es zu testen, laden Sie bitte die obs.mat-Datei herunterhttps: //drive.google.com/file/d/0B3vXKJ_zYaCJbHU2SHhac29MRms/view? usp = sharing