エントリーの編集
![loading...](https://b.st-hatena.com/bdefb8944296a0957e54cebcfefc25c4dcff9f5f/images/v4/public/common/loading@2x.gif)
エントリーの編集は全ユーザーに共通の機能です。
必ずガイドラインを一読の上ご利用ください。
StanでWiener分布からの乱数を得る方法(+累積分布関数ほか)
記事へのコメント0件
- 注目コメント
- 新着コメント
このエントリーにコメントしてみましょう。
注目コメント算出アルゴリズムの一部にLINEヤフー株式会社の「建設的コメント順位付けモデルAPI」を使用しています
![アプリのスクリーンショット](https://b.st-hatena.com/bdefb8944296a0957e54cebcfefc25c4dcff9f5f/images/v4/public/entry/app-screenshot.png)
- バナー広告なし
- ミュート機能あり
- ダークモード搭載
関連記事
StanでWiener分布からの乱数を得る方法(+累積分布関数ほか)
Stanには,Drift Diffusion Modelを実装するのにも使えるWiener分布の確率密度関数が実装されています。... Stanには,Drift Diffusion Modelを実装するのにも使えるWiener分布の確率密度関数が実装されています。サンプリングステートメントとしてのy ~ wiener(alpha, tau, beta, delta)と対数確率密度を計算するwiener_lpdf(y | alpha, tau, beta, delta)を使えばデータをWiener分布に当てはめることは容易に行えますが,Wiener分布からの乱数を発生させる関数wiener_rng()は実装されていません。Stanで得たMCMCサンプルを使ってR側でRWiener::rwiener()関数を使えば乱数を得ることはできますが,事後予測分布を作りたいときにはStan側で実装できたほうが便利だと思います。そんな訳で,Stan用のwiener_rng()を自作してみました。2通りのやり方を思い付いたので2つ紹介します