気になったのでパラパラ調べてみたんだが、ここら辺のトピックもまだまだやることあるんだなあ。 実証へ繋がる可能性も多い面白い分野でもある。 暇があったらRで実装されていない単位根検定パッケージを書いてみたい。難しそうなのは、極限分布を作るとコラ…

と言うか、数値計算ライブラリ一般に言えることなんだけど、中身がよく分からなくて怖いことがある。 例えば、固有値計算とかするときに複素固有値とか出てきしまったときの処理とかどうやってるんだろうなあとか。アルゴリズム通りに実装できてるのかとか、…

Rのoptimを使った汎用最適化は非常に遅い。 本当は、最尤推定法の関数を使いたいんだけど、出てくる値が怪しい。 とか、色々書こうと思ったけど面倒くさいからやめる。 取りあえず、特殊な組み込み関数で実装されている統計モデルの推定に使われている最適化…