- 162 名前：nobodyさん mailto:sage [2009/11/29(日) 13:23:41 ID:???]
- >>150
JavaScriptでの実装を見つけました。
正規分布の確率計算
econom01.cc.sophia.ac.jp/stat/normprob.htm
確率密度関数をテイラー展開して積分しているようです。
下記サイトにその説明があります。
誤差関数の近似式（その１）
www.geocities.jp/ikuro_kotaro/koramu/gosa.htm
同じようにPerlで実装すると、>>143よりさらに簡単になりました。
収束も早いようです。
sub normdist {
my %param = @_;
use Math::Trig qw( pi );
my $precision = 10**-$param{precision};
my $z = ( $param{x} - $param{mean} ) / $param{sigma};
my $z2 = $z * $z;
my $p = my $s = $z * exp( -0.5 * $z2 ) / sqrt( 2*pi() );
for ( my $n = 3; abs($s) > $precision; $n += 2 ) {
$s *= $z2 / $n;
$p += $s;
}
return 0.5 + $p;
}
|
 |
