Rのfor文を回していると、途中でエラーが出て止まってしまいました。 その解決策としてエラーが出てもとりあえず最後まで回してくれるtry関数を紹介します。 以下のコードは10次元ベクトルの成分をlogで生成。 logを計算する際、はからのランダムサンプルで…

Rで最尤法を行おうと思ったのですが、ネットで調べてもいまいち分からなかったため自分で作ってみました。 自然な設定として線形回帰モデルを考え、最尤推定量を求めます。 モデル： ここで. 誤差項の分布として を仮定。 設定: N <- 1000 # サンプル数 mu <…