Xi (i = 1~n)を正規母集団N(μ、σ2)からのi.i.d.サンプルとする。Rを用いたGibbsサンプラーによるμとσ2のベイズ推定例を示す。 (参考文献)中妻照雄著「入門ベイズ統計学」、朝倉書店、p.145。 1.事前分布 2.尤度関数 3.事後分布 4.未知パラメータの全条件付事後分布 5.Gibbsサンプラーのアルゴリズム 6.Rコードの一例 library(MCMCpack) # rinvgamma()を使用するため X <- c(0.0459, 0.0436, 0.0207, 0.0867, 0.1678, 0.1748) #データ n <- NROW(X) Xbar <- mean(X) Sum2 <- sum((X-Xbar)^2) mu0 <- 0.0 # ハイパーパラメータ sig20 <- 0.01 # ハイパーパラメータ (sig20はτ02を表す。) nu0