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