我有一组观察结果obs
。
如果我绘制观察的直方图,我发现它们可能来自伽马分布
[counts,x] = hist(obs,[1:max(obs)]);
我想用卡方拟合优度来证明它。
所以我首先估计伽玛参数
paramEsts = fitdist(obs,'Gamma');
并使用 chi2gof 查看假设是否为真(h = 0)。
[h,p] = chi2gof(obs,'CDF',paramEsts)
我的问题是我得到 ap = NaN ...这怎么可能?我的错误在哪里?谢谢
要对其进行测试,请下载 obs.mat 文件 https://drive.google.com/file/d/0B3vXKJ_zYaCJbHU2SHhac29MRms/view?usp=sharing