仕様上はできそうなので試してみたらできたと言うだけの話なのですが、Rで呼び出したFortranのプロシージャからRの関数を呼ぶ方法を記しておきます。技術的にどうこうと言うのは無いのですが、流れはややこしくて以下のようになります。 .CallでRのコールバ…

Rもだいたい線形代数の範囲で片付けている限りは、そこそこ速いことがわかります。Fortranとの速度差が出ないというか、R vs R/FortranだとRの方が微妙に速く*1。遅くなるのになぜFortranでRの対数尤度関数とその一階微分を用意したかと言うと、Fortranでハ…