谁用SAS做过负二项分布的拟合，求助！急急急

shiyiming

举例：钉螺分布（摘自杨树勤的中国医学百科全书。医学统计学）
每框螺数X   相应出现的框数即频数f
0           3836
1           1170

已经建立好数据集如下：
data a;
input x f@@;
lines;
0 3836 1 1170 2 444 3 174 4 71 5 31 6 18 7 7 8 3 9 2 10 2 11 1 12 0 13 0 14 1
;
run;

通过该数据集如何拟合负二项分布，看是否符合该分布，非常感谢！！

shiyiming

手工算，不行吗？

利用卡方拟合优度检验。try it!  <!-- s:D --><img src="{SMILIES_PATH}/icon_biggrin.gif" alt=":D" title="Very Happy" /><!-- s:D -->

shiyiming

行，但是数据行数比较多时就太累啦 <!-- s:P --><img src="{SMILIES_PATH}/icon_razz.gif" alt=":P" title="Razz" /><!-- s:P -->
书上的数据量仅10几行，计算量比较少，但是现在有上百行，太恐怖啦
我折腾了半天，每搞定： 问题是计算概率是当x=0时，概率是比较简单的，但是当x^=0时需要用递推公式，每一次都要用到前一个概率，而这所有的概率应是在一列而不是一行上，有点麻烦 <!-- s:oops: --><img src="{SMILIES_PATH}/icon_redface.gif" alt=":oops:" title="Embarassed" /><!-- s:oops: -->
有兴趣可以试试帮忙解决，T U