目標分布 目標分布をベータ分布とします。 ベータ分布の描画のためにライブラリを使います。MCMC時には正規化定数は計算しません。 # 目標分布 def p(x): return (x ** (a - 1)) * ((1 - x) ** (b - 1)) # カーネルのみを計算 提案分布 ベータ分布の定義域がなので、 提案分布は正規分布の確率変数の実現値が0以上の箇所を使うことにします。 具体的には、平均0の正規分布で乱数を発生させたのち絶対値を取ります。イメージは以下です。 赤色の線は平均0の正規分布です。 緑色の線が今回採用する提案分布です。 ソース ####################…