StatModeling https://blog.hatena.ne.jp/StatModeling/ StatModeling Memorandum https://statmodeling.hatenablog.com/ Stan R 前回の理論的なまとめを踏まえてStanでWAICとWBICを計算してみます。 今回は例題として混合正規分布から発生させたデータ100個を用いて、2種類のモデルで推定を行い、それぞれに対してWAICとWBICを求めてみます。まずはデータ生成部分のRコードは以下です。 N <- 100 a_true <- 0.4 mean1 <- 0 mean2 <- 3 sd1 <- 1 sd2 <- 1 set.seed(1) Y <- c(rnorm((1-a_true)*N, mean1, sd1), rnorm(a_true*N, mean2, sd2)) 次にモデルその1の説明です。ここでは2つの正規… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Fstatmodeling.hatenablog.com%2Fentry%2Fcalc-waic-wbic" title="WAICとWBICを事後分布から計算する - StatModeling Memorandum" class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> https://blog.st-hatena.com/images/theme/og-image-1500.png Hatena Blog https://hatena.blog 2015-09-12 00:00:00 rich https://statmodeling.hatenablog.com/entry/calc-waic-wbic 1.0 100%