1. 学ぶ
    2. /
  2. コース
    3. /
  3. RJAGS によるベイズモデリング
Connected

演習

マルコフ連鎖のトレースプロット

トレースプロットは、マルコフ連鎖の時系列的な挙動を可視化するグラフです。具体的には、\(m\) 連鎖のトレースプロットは、x 軸に反復番号、y 軸に連鎖の観測値をプロットします。

\(m\) 連鎖のトレースプロットを、2 つの方法で作成しましょう。1 つ目は、組み込みの plot() 関数を mcmc.list オブジェクトの sleep_sim に適用する方法です。2 つ目は、グラフをより細かく制御するために（後の章での分析でも役立ちます）、data.frame オブジェクトの sleep_chains に ggplot() を適用する方法です。sleep_sim と sleep_chains はどちらもワークスペースに用意されています。

sleep_sim <- coda.samples(model = sleep_jags, variable.names = c("m", "s"), n.iter = 10000)
sleep_chains <- data.frame(sleep_sim[[1]], iter = 1:10000)

指示

100 XP

  • sleep_sim に対して density = FALSE を指定して plot() を適用し、\(m\) 連鎖と \(s\) 連鎖のトレースプロットを作成します。注意：記録された 10,000 回の反復は、サンプルが破棄される「バーンイン」期間の後から始まります。そのため、反復カウントは 1 からスタートしません。

  • sleep_chains に対して geom_line() レイヤーを加えた ggplot() を適用し、\(m\) 連鎖のトレースプロットを再作成します。

  • ズームイン：\(m\) 連鎖の最初の 100 反復に絞ったトレースプロットを ggplot() で作成します。