統計とか,研究とか,日常についての備忘録的なもの。この広告は、90日以上更新していないブログに表示しています。2017-12-21JAGSで対数尤度を得る方法(Stanも)RWAICなどを計算するには,対数尤度のMCMCサンプルが必要。個々のMCMCサンプルからがんばって計算する方法もあるが,より簡単にはMCMCサンプラーのなかで対数尤度を得ると楽。quantityが容易されているが,JAGSでは少し工夫が必要。参考https://sourceforge.net/p/mcmc-jags/discussion/610037/thread/a3bef062/https://sourceforge.net/p/mcmc-jags/discussion/610037/thread/87902622/その他参考ito-hi.blog.so-net.ne.jphamada.hatenablog.jpJAGSのモデルの中でのようなものでも対数尤度が得られる一変量の分布に関しては,密度関数などがRのように用意されている。logdensity.**系の関数が用意されているが,これはマニュアルに書かれていないので必要に応じて調べる。相関係数については,以下の頁のものを参考にした。logdensity.mnormで多変量正規分布の対数密度を得られる。ベクトルの配列について注意。y[i]で得られるのはi番目の長さが2のベクトル。y[i][1]でi番目の配列にあるベクトルの1番目の要素を取り出すということになる。値は大体一致している。https://mc-stan.org/docs/2_18/stan-users-guide/multi...2019-06-16とある論文を読んで考えたこと下のリンクの論文を読んで,ベイズ統計について,少し気になる...2018-07-02ベイズ統計学について最近思うこと最近,ベイズ統計学が流行っていることに関して,個人的に思っ...2017-12-23JAGSで予測分布の作成その2IRTモデルで予測分布や対数尤度の計算などをおこなってみる。putについて月別アーカイブはてなブログをはじめよう!kazzstatさんは、はてなブログを使っています。ブログを報告するスターをつけました引用をストックしましたストック一覧を見る閉じる引用するにはまずログインしてくださいログイン閉じる引用をストックできませんでした。再度お試しください閉じる限定公開記事のため引用できません。