「ち」は「ち」でも値の話

誰か検算してくれ、と書いたが、よく考えればMathematicaに検算させればいいんじゃないか。

俺はどうも数式で展開するのが苦手で、実際に数字を一つずつ見ていかないと納得できないのだけど、おかげで計算プロセスを電卓で追いかけることが出来るのだ。

面倒な手続きではあるが、一つずつ計算してみると、自分のプログラムと近いようでずれているところがあったりする。
やっぱり精度の問題みたいだ。

Mathematicaで計算すると、イコール・ゼロになるべき数値が、-0.01235と出た。うむむむぅー、これはどうやら正解っぽいですぞ。となると、自分のソフトの精度問題だなぁ。

問題は二変数正規分布のCDFの精度のようです。まぁこれの精度の上げ方は、原典に書いてあったからなんとかなるかも。参考になりそうな本もわかったし(Mathematicaのヘルプで出てきた)。

こういうかすかな光に導かれて、やる気を出しちゃうんだよな。

余談ですが、フリーの数値計算ソフトとして有名なMaximaOctaveには、多変量正規分布関数のCDFを求めるような関数は用意されていないみたいだ。いや、探し出せてないだけかもしれないけど、ヘルプ等々が十分でないといういみではやっぱりこれらのソフトはまだまだ、で。Mathematicaってすげぇなぁ。お金があったら、新しいバージョンを買いたいです。

日記
連休、やってます

令和になりましたな。 世間は10連休。わたしは初日の土曜日に研究会を入れたし、最後の月曜日は授業日な …

日記
バルサン焚きました

去年の今頃もやったなぁと思いつつ。 家の中にムカデが出るんですよ。家族はみんな虫嫌いで、中でもムカデ …

日記
さようなら平成30年度

大学という年度単位で生きている私どもにとっては,今日がある意味一年の最後の日です。明日は新元号も,科 …