论坛 › 算法及编程语言 › [求助]卡方分布函数？？

noahhu 发表于 2006-4-4 23:00

[求助]卡方分布函数？？

以下是我参照数理统计实用算法编写的卡方分布函数，这个程序在df为偶数时正确，但df为奇数时计算结果不正确，请各位高手帮忙，看看什么原因。另外，论坛中<font size="-1"><font color="#cc0033">Handbook Of Mathematical Functions找不到，如果哪位能下载请帮我下载一下， 我的邮箱是：dreams_wing@163.com 谢谢了！</font></font><br>double CDFChi(double X2, int df){<br>    int n=df, IAI;<br>    double F, u, pp;<br>    if (X2==0) return 0.0;<br>    double PIS=sqrt(PI());<br>    double x=X2/2;<br>    double CHS=sqrt(x);<br><br>    if (n%2==0) {<br>      u=x*exp(-x);<br>      F=1-exp(-x);<br>      IAI=2;<br>    }<br>    else {<br>      u=sqrt(x)*exp(-x)/PIS;<br>      pp=CDFNor(CHS);<br>      F=2*(pp-0.5);<br>      IAI=1;<br>    };<br><br>    if (IAI==n) return F;<br>    double n2=n-2;<br>    for (double i=IAI; i&lt;=n2; i+=2) {<br>      F=F-2*u/i;<br>      u=X2*u/i;<br>    };<br>    return F;<br>};

johhan 发表于 2006-4-5 15:03

C语言编的,俺们看不懂!!!!!!!!!!!

查看完整版本: [求助]卡方分布函数？？