Tengo un conjunto de observacionesobs.

Si trazo el histograma de la observación, veo que podrían provenir de una distribución gamma

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

Me gustaría probarlo usando chi cuadrado bondad de ajuste.

Así que primero calculo los parámetros gamma

paramEsts = fitdist(obs,'Gamma');  

y use chi2gof para ver si la hipotesis es verdadera (h = 0).

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

Mi problema es que obtengo un p = NaN ... ¿cómo es posible? ¿Dónde está mi error? Gracias

para probarlo, descargue el archivo obs.mathttps://drive.google.com/file/d/0B3vXKJ_zYaCJbHU2SHhac29MRms/view?usp=sharing