{ "cells": [ { "cell_type": "markdown", "metadata": { "tags": [] }, "source": [ "# MCMC diagnostics via a case study: Artificial funnel of hell {#sec-mcmc-diagnostics}\n", "\n", "
\\n\"+\n", " \"BokehJS does not appear to have successfully loaded. If loading BokehJS from CDN, this \\n\"+\n", " \"may be due to a slow or bad network connection. Possible fixes:\\n\"+\n", " \"
\\n\"+\n", " \"\\n\"+\n",
" \"from bokeh.resources import INLINE\\n\"+\n",
" \"output_notebook(resources=INLINE)\\n\"+\n",
" \"\\n\"+\n",
" \"\\n\"+\n \"BokehJS does not appear to have successfully loaded. If loading BokehJS from CDN, this \\n\"+\n \"may be due to a slow or bad network connection. Possible fixes:\\n\"+\n \"
\\n\"+\n \"\\n\"+\n \"from bokeh.resources import INLINE\\n\"+\n \"output_notebook(resources=INLINE)\\n\"+\n \"\\n\"+\n \"<xarray.Dataset> Size: 16B\n",
"Dimensions: ()\n",
"Data variables:\n",
" theta float64 8B 1.054\n",
" v float64 8B 1.099\n",
"Attributes:\n",
" created_at: 2026-02-09T19:09:04.506296+00:00\n",
" arviz_version: 0.23.0\n",
" inference_library: cmdstanpy\n",
" inference_library_version: 1.3.0| \n", " | mean | \n", "sd | \n", "hdi_3% | \n", "hdi_97% | \n", "mcse_mean | \n", "mcse_sd | \n", "ess_bulk | \n", "ess_tail | \n", "r_hat | \n", "
|---|---|---|---|---|---|---|---|---|---|
| theta | \n", "0.060 | \n", "7.233 | \n", "-9.078 | \n", "8.971 | \n", "0.478 | \n", "1.507 | \n", "467.0 | \n", "298.0 | \n", "1.05 | \n", "
| v | \n", "0.667 | \n", "2.627 | \n", "-4.464 | \n", "5.145 | \n", "0.392 | \n", "0.313 | \n", "41.0 | \n", "18.0 | \n", "1.10 | \n", "
<xarray.Dataset> Size: 16B\n",
"Dimensions: ()\n",
"Data variables:\n",
" theta float64 8B 1.003\n",
" v float64 8B 1.005\n",
"Attributes:\n",
" created_at: 2026-02-09T19:17:39.526236+00:00\n",
" arviz_version: 0.23.0\n",
" inference_library: cmdstanpy\n",
" inference_library_version: 1.3.0| \n", " | mean | \n", "sd | \n", "hdi_3% | \n", "hdi_97% | \n", "mcse_mean | \n", "mcse_sd | \n", "ess_bulk | \n", "ess_tail | \n", "r_hat | \n", "
|---|---|---|---|---|---|---|---|---|---|
| theta | \n", "0.060 | \n", "7.233 | \n", "-9.078 | \n", "8.971 | \n", "0.478 | \n", "1.507 | \n", "467.0 | \n", "298.0 | \n", "1.05 | \n", "
| v | \n", "0.667 | \n", "2.627 | \n", "-4.464 | \n", "5.145 | \n", "0.392 | \n", "0.313 | \n", "41.0 | \n", "18.0 | \n", "1.10 | \n", "
<xarray.DataArray 'diverging' (chain: 4, draw: 1000)> Size: 4kB\n",
"array([[False, False, False, ..., False, False, False],\n",
" [False, False, False, ..., False, False, False],\n",
" [ True, False, True, ..., False, False, False],\n",
" [ True, False, False, ..., False, False, False]], shape=(4, 1000))\n",
"Coordinates:\n",
" * chain (chain) int64 32B 0 1 2 3\n",
" * draw (draw) int64 8kB 0 1 2 3 4 5 6 7 ... 993 994 995 996 997 998 999<xarray.DataArray 'tree_depth' (chain: 4, draw: 1000)> Size: 32kB\n",
"array([[2, 4, 4, ..., 2, 3, 4],\n",
" [3, 3, 2, ..., 2, 1, 2],\n",
" [1, 1, 4, ..., 4, 4, 4],\n",
" [1, 1, 1, ..., 4, 2, 2]], shape=(4, 1000))\n",
"Coordinates:\n",
" * chain (chain) int64 32B 0 1 2 3\n",
" * draw (draw) int64 8kB 0 1 2 3 4 5 6 7 ... 993 994 995 996 997 998 999<xarray.DataArray 'energy' (chain: 4, draw: 1000)> Size: 32kB\n",
"array([[ 2.38275 , 4.71404 , 4.36306 , ..., 1.57001 , 3.59906 ,\n",
" 3.32988 ],\n",
" [ 0.717477, 2.60399 , 3.26155 , ..., 0.396415, 2.11211 ,\n",
" 2.62692 ],\n",
" [-0.789693, -0.621038, 2.87434 , ..., 3.22471 , 1.53535 ,\n",
" 1.90595 ],\n",
" [ 0.604643, -0.27291 , 1.58246 , ..., 2.04745 , 0.59004 ,\n",
" 0.525632]], shape=(4, 1000))\n",
"Coordinates:\n",
" * chain (chain) int64 32B 0 1 2 3\n",
" * draw (draw) int64 8kB 0 1 2 3 4 5 6 7 ... 993 994 995 996 997 998 999