updated

Author: abdullahau

03 - Geocentric Models.ipynb

@@ -229516,7 +229516,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 139,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -229548,8 +229548,8 @@
     "    axs[1].plot(data.year, B_mat * wp)\n",
     "    axs[1].set(xlabel='Year', ylabel='Basis * Weight')    \n",
     "    \n",
-    "    az.plot_hdi(d2.year, trace_quap.posterior['mu'], smooth=False, color='k', ax=axs[2])\n",
-    "    axs[2].plot(d2.year, d2.doy, 'o', color='pink', alpha=0.6)\n",
+    "    az.plot_hdi(data.year, trace_quap.posterior['mu'], smooth=False, color='k', ax=axs[2])\n",
+    "    axs[2].plot(data.year, data.doy, 'o', color='pink', alpha=0.6)\n",
     "    axs[2].set(xlabel='Year', ylabel='Day in Year')"
    ]
   },

03b - Geocentric Models.ipynb

@@ -80038,7 +80038,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 41,
+   "execution_count": null,
    "metadata": {},
    "outputs": [
     {
@@ -82048,16 +82048,16 @@
     "fig, ax = plt.subplots()\n",
     "\n",
     "hdi_1 = az.plot_hdi(weight_seq, pred.predictions[\"height\"], \n",
-    "                     hdi_prob=0.89, smooth=False)  \n",
+    "                     hdi_prob=0.89, smooth=False, color='C1')  \n",
     "\n",
     "hdi_2 = az.plot_hdi(weight_seq, pred.predictions[\"mu\"], \n",
-    "                     hdi_prob=0.89, smooth=False) \n",
+    "                     hdi_prob=0.89, smooth=False, color='green') \n",
     "\n",
     "mean_line, = ax.plot(weight_seq, pred.predictions[\"mu\"].mean(axis=1)[0], 'k', label='Mean')\n",
     "obs_data, = ax.plot(d2.weight, d2.height, 'o', fillstyle='none', label='Observed Data', color='k')\n",
     "\n",
-    "hdi_patch1 = plt.Line2D([0], [1], color='C0', lw=3, label=\"Posterior Predictive Sample\")\n",
-    "hdi_patch2 = plt.Line2D([0], [1], color='C1', lw=3, label=\"Uncertainty in Mean\")\n",
+    "hdi_patch1 = plt.Line2D([0], [1], color='C1', lw=3, label=\"Posterior Predictive Sample\")\n",
+    "hdi_patch2 = plt.Line2D([0], [1], color='green', lw=3, label=\"Uncertainty in Mean\")\n",
     "\n",
     "ax.set(xlabel=\"Weight (kg)\",\n",
     "       ylabel=\"Height (cm)\",\n",