1

我有一组观察结果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

4

0 回答 0