程序如何修改，或如何实现更好！

我觉得，楼主不会不明白我问的是什么。

简单说，“板凳”处给的程序的核心即，针对pdf文件指出的那个多项式，计算多项式中各次项（共有0、1、……、(n+1)n/2等次数项）的系数——个人以为此即本题的关键。

至于这些系数与可能的秩和及相应概率有什么关系，我想楼主应该是清楚的（抱歉，本人尚不清楚，但本人很想从楼主那里知道，虽然到现在也不知道）。楼主自己做出相应修改即可。

楼主虽然前面说了pdf文件“说得清楚”，但楼主前面似乎一直并没有按照pdf文件的表述及逻辑来说明这个问题（如果楼主已经知道了其他语言的程序代码，并且这些代码是有效的，我猜想Stata在program define这方面未必有优势）。

sungmoo 发表于 2012-12-29 19:12
我觉得，楼主不会不明白我问的是什么。

简单说，“板凳”处给的程序的核心即，针对pdf文件指出的那个多项 ...

谢谢sungmoo老师！
程序没问题！

sungmoo 发表于 2012-12-28 23:41
*上面的程序可以大幅减化（首数是样本量，末数是可能的秩和，求给定首数时末数对应的概率）
clear
cap pr ...

int(`e')!=`e' 这个条件可以去掉。