P値（統計学）を求めるプログラム

質問

統計学でP値を求めるプログラムで、通常ならば、台形公式やシンプソン公式などを用いて、微小ずつ積分していくと考え作ってみたのですが、先生から毎回積分するのではなく近似値を出しておいて利用するようにと言われ、積分をしないプログラムを組めと言われています。
おそらく、積分はしたという前提で計算するという事でしょうか？

フリーの統計ソフトR　のコードを参考にと見ていたんですが理解し難く、どなたかにご教授願いたいと考えました。
よろしくお願い致します。
今回は、カイ二乗の値からP値を求めようと考えています。
その際の、カイ二乗分布に従う式は

f(x) = x^(n-2/2) * e^(-x/2) / 2^(n/2)gamma(n/2) は分かりました。

ガンマ関数に関しては、
double gamma(double x, int *ier) /* ガンマ関数の計算 */
{
double err, g, s, t, v, w, y;
long k;

*ier = 0;

if (x > 5.0) {
v = 1.0 / x;
s = ((((((-0.000592166437354 * v + 0.0000697281375837) * v +
0.00078403922172) * v - 0.000229472093621) * v -
0.00268132716049) * v + 0.00347222222222) * v +
0.0833333333333) * v + 1.0;
g = 2.506628274631001 * exp(-x) * pow(x,x-0.5) * s;
}

else {

err = 1.0e-20;
w = x;
t = 1.0;

if (x < 1.5) {

if (x < err) {
k = (long)x;
y = (double)k - x;
if (fabs(y) < err || fabs(1.0-y) < err)
*ier = -1;
}

if (*ier == 0) {
while (w < 1.5) {
t /= w;
w += 1.0;
}
}
}

else {
if (w > 2.5) {
while (w > 2.5) {
w -= 1.0;
t *= w;
}
}
}

w -= 2.0;
g = (((((((0.0021385778 * w - 0.0034961289) * w +
0.0122995771) * w - 0.00012513767) * w + 0.0740648982) * w +
0.0815652323) * w + 0.411849671) * w + 0.422784604) * w +
0.999999926;
g *= t;
}

return g;
}

回答No.1

SE-1 · ベストアンサー

私は特殊関数の計算や統計処理をするプログラムを組む際、以下のサイトに大変お世話になりました。

ガンマ関数：function.c
カイ二乗：dist.c

参考URL：

http://www5.airnet.ne.jp/tomy/cpro/science.htm