У меня есть множество наблюденийobs.

Если я построю гистограмму наблюдения, я вижу, что они могут быть получены из гамма-распределения

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

Я хотел бы доказать это, используя добротность подбора.

Поэтому я сначала оцениваю параметры гаммы

paramEsts = fitdist(obs,'Gamma');  

и используйте chi2gof, чтобы увидеть, является ли гипотеза истинной (h = 0).

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

Моя проблема в том, что я получаю р = NaN ... как это возможно? Где моя ошибка? Спасибо

чтобы проверить это, пожалуйста, скачайте файл obs.mathttps://drive.google.com/file/d/0B3vXKJ_zYaCJbHU2SHhac29MRms/view?usp=sharing