「統計学関連なんでもあり」の過去ログ---002

　306　スチューデント化された範囲　　青木繁伸　　1999/08/19 (木) 10:38
　　307　Re: スチューデント化された範囲　　青木繁伸　　1999/08/19 (木) 12:10

306.　スチューデント化された範囲　　青木繁伸　　1999/08/19 (木) 10:38

夏休み前に，堀さんのページ
http://www.ec.kagawa-u.ac.jp/~hori/delphistat/index.html
を参考にして，多重比較の Tukey の方法を JavaScript で書いてみたりしました。ルジャンドル多項式の定数の在処などについて，たいへんありがたかったです。

堀さんは，「qerf, qerfc については大浦拓哉（京都大学数理解析研究所）氏」のプログラムを pascal に移植したと書かれてありましたが，大浦さんはなぜこの二つを別のプログラムにしたのでしょうかというのがちょっと疑問でした。

JavaScript はかなり遅い（私の Macintosh が遅いだけかも知れないが）ので，元々のプログラムの計算精度がよすぎるのをちょっと改悪して，計算時間を半分くらいに短縮してみました。

　　　　　[このページのトップへ]

307.　Re: スチューデント化された範囲　　青木繁伸　　1999/08/19 (木) 12:10

Numerical recipes in C にあったものは以下のようになっていました。

function qerfc(x) { var t,z,ans; z=Math.abs(x); t=1.0/(1.0+0.5*z); ans=t*Math.exp(-z*z-1.26551223+t*(1.00002368 +t*(0.37409196+t*(0.09678418+t*(-0.18628806 +t*(0.27886807+t*(-1.13520398+t*(1.48851587 +t*(-0.82215223+t*0.17087277))))))))) return x >= 0.0 ? ans : 2.0-ans }
これを使って，

function qerf(x) { return 1-qerfc(x) }

　　　　　[このページのトップへ]