{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Applying a discrete-time, non-stationary nucleotide model\n",
    "\n",
    "We fit a discrete-time Markov nucleotide model. This corresponds to a Barry and Hartigan 1987 model."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<style>\n",
       "tr:last-child {border-bottom: 1px solid #000;} tr > th {text-align: center !important;} tr > td {text-align: left !important;}\n",
       "</style>\n",
       "<caption style=\"color: rgb(250, 250, 250); background: rgba(30, 140, 200, 1); align=top;\"><span style=\"font-weight: bold;\">BH</span><span></span></caption>\n",
       "<thead style=\"background: rgba(161, 195, 209, 0.75); font-weight: bold; text-align: center;\">\n",
       "<th>key</th>\n",
       "<th>lnL</th>\n",
       "<th>nfp</th>\n",
       "<th>DLC</th>\n",
       "<th>unique_Q</th>\n",
       "</thead>\n",
       "<tbody>\n",
       "<tr>\n",
       "<td></td>\n",
       "<td style=\"font-family: monospace !important;\">-6941.4684</td>\n",
       "<td style=\"font-family: monospace !important;\">135</td>\n",
       "<td>True</td>\n",
       "<td></td>\n",
       "</tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/plain": [
       "BH\n",
       "============================================\n",
       "key           lnL    nfp     DLC    unique_Q\n",
       "--------------------------------------------\n",
       "       -6941.4684    135    True            \n",
       "--------------------------------------------\n",
       "\n",
       "1 rows x 5 columns"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from cogent3.app import io, evo\n",
    "\n",
    "loader = io.load_aligned(format=\"fasta\", moltype=\"dna\")\n",
    "aln = loader(\"../data/primate_brca1.fasta\")\n",
    "model = evo.model(\"BH\", tree=\"../data/primate_brca1.tree\")\n",
    "result = model(aln)\n",
    "result"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "**NOTE:** DLC stands for diagonal largest in column and the value is a check on the identifiability of the model. `unique_Q` is not applicable to a discrete-time model and so remains as `None`.\n",
    "\n",
    "Looking at the likelihood function, you will"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<h4>BH</h4>\n",
       "<p>log-likelihood = -6941.4684</p>\n",
       "<p>number of free parameters = 135</p>\n",
       "<table>\n",
       "<style>\n",
       "tr:last-child {border-bottom: 1px solid #000;} tr > th {text-align: center !important;} tr > td {text-align: left !important;}\n",
       "</style>\n",
       "<caption style=\"color: rgb(250, 250, 250); background: rgba(30, 140, 200, 1); align=top;\"><span style=\"font-weight: bold;\">Edge motif motif2 params</span><span></span></caption>\n",
       "<thead style=\"background: rgba(161, 195, 209, 0.75); font-weight: bold; text-align: center;\">\n",
       "<th>edge</th>\n",
       "<th>motif</th>\n",
       "<th>motif2</th>\n",
       "<th>psubs</th>\n",
       "</thead>\n",
       "<tbody>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">Galago</td>\n",
       "<td>T</td>\n",
       "<td>T</td>\n",
       "<td style=\"font-family: monospace !important;\">0.8751</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">Galago</td>\n",
       "<td>T</td>\n",
       "<td>C</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0649</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">Galago</td>\n",
       "<td>T</td>\n",
       "<td>A</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0409</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">Galago</td>\n",
       "<td>T</td>\n",
       "<td>G</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0192</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">Galago</td>\n",
       "<td>C</td>\n",
       "<td>T</td>\n",
       "<td style=\"font-family: monospace !important;\">0.1126</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">...</td>\n",
       "<td>...</td>\n",
       "<td>...</td>\n",
       "<td>...</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">edge.3</td>\n",
       "<td>A</td>\n",
       "<td>G</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0055</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">edge.3</td>\n",
       "<td>G</td>\n",
       "<td>T</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0000</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">edge.3</td>\n",
       "<td>G</td>\n",
       "<td>C</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0011</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">edge.3</td>\n",
       "<td>G</td>\n",
       "<td>A</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0039</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">edge.3</td>\n",
       "<td>G</td>\n",
       "<td>G</td>\n",
       "<td style=\"font-family: monospace !important;\">0.9950</td>\n",
       "</tr>\n",
       "</tbody>\n",
       "</table>\n",
       "\n",
       "<table>\n",
       "<style>\n",
       "tr:last-child {border-bottom: 1px solid #000;} tr > th {text-align: center !important;} tr > td {text-align: left !important;}\n",
       "</style>\n",
       "<caption style=\"color: rgb(250, 250, 250); background: rgba(30, 140, 200, 1); align=top;\"><span style=\"font-weight: bold;\">Motif params</span><span></span></caption>\n",
       "<thead style=\"background: rgba(161, 195, 209, 0.75); font-weight: bold; text-align: center;\">\n",
       "<th>A</th>\n",
       "<th>C</th>\n",
       "<th>G</th>\n",
       "<th>T</th>\n",
       "</thead>\n",
       "<tbody>\n",
       "<tr>\n",
       "<td style=\"font-family: monospace !important;\">0.3774</td>\n",
       "<td style=\"font-family: monospace !important;\">0.1763</td>\n",
       "<td style=\"font-family: monospace !important;\">0.2058</td>\n",
       "<td style=\"font-family: monospace !important;\">0.2404</td>\n",
       "</tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/plain": [
       "<cogent3.evolve.parameter_controller.AlignmentLikelihoodFunction at 0x1111da5f8>"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "result.lf"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Get a tree with branch lengths as paralinear\n",
    "\n",
    "This is the only possible length metric for a discrete-time process."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "        <script type=\"text/javascript\">\n",
       "        window.PlotlyConfig = {MathJaxConfig: 'local'};\n",
       "        if (window.MathJax) {MathJax.Hub.Config({SVG: {font: \"STIX-Web\"}});}\n",
       "        if (typeof require !== 'undefined') {\n",
       "        require.undef(\"plotly\");\n",
       "        requirejs.config({\n",
       "            paths: {\n",
       "                'plotly': ['https://cdn.plot.ly/plotly-latest.min']\n",
       "            }\n",
       "        });\n",
       "        require(['plotly'], function(Plotly) {\n",
       "            window._Plotly = Plotly;\n",
       "        });\n",
       "        }\n",
       "        </script>\n",
       "        "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "application/vnd.plotly.v1+json": {
       "config": {
        "plotlyServerURL": "https://plot.ly"
       },
       "data": [
        {
         "line": {
          "color": "black",
          "shape": "spline",
          "smoothing": 1.3,
          "width": 1.25
         },
         "mode": "lines",
         "showlegend": false,
         "type": "scatter",
         "x": [
          null,
          null,
          0,
          0,
          null,
          0,
          0.03139195169925202,
          null,
          0.03139195169925202,
          0.03139195169925202,
          null,
          0.03139195169925202,
          0.08159484122097149,
          null,
          0.08159484122097149,
          0.08159484122097149,
          null,
          0.08159484122097149,
          0.09617900440577909,
          null,
          0.09617900440577909,
          0.09617900440577909,
          null,
          0.09617900440577909,
          0.0961878725442622,
          null,
          0.0961878725442622,
          0.0961878725442622,
          null,
          0.0961878725442622,
          0.10959571362620046,
          null,
          0.0961878725442622,
          0.0961878725442622,
          null,
          0.0961878725442622,
          0.12185690868725763,
          null,
          0.09617900440577909,
          0.09617900440577909,
          null,
          0.09617900440577909,
          0.1057597283305971,
          null,
          0.08159484122097149,
          0.08159484122097149,
          null,
          0.08159484122097149,
          0.11501597595635671,
          null,
          0.03139195169925202,
          0.03139195169925202,
          null,
          0.03139195169925202,
          0.12464469563427283,
          null,
          0,
          0,
          null,
          0,
          0.1930734042426634,
          null,
          0,
          0,
          null,
          0,
          0.7165441919466353,
          null
         ],
         "y": [
          null,
          null,
          -2.6,
          -0.08125000000000004,
          null,
          -0.08125000000000004,
          -0.08125000000000004,
          null,
          -0.08125000000000004,
          1.1375,
          null,
          1.1375,
          1.1375,
          null,
          1.1375,
          2.275,
          null,
          2.275,
          2.275,
          null,
          2.275,
          3.25,
          null,
          3.25,
          3.25,
          null,
          3.25,
          3.9000000000000004,
          null,
          3.9000000000000004,
          3.9000000000000004,
          null,
          3.25,
          2.6,
          null,
          2.6,
          2.6,
          null,
          2.275,
          1.3,
          null,
          1.3,
          1.3,
          null,
          1.1375,
          0,
          null,
          0,
          0,
          null,
          -0.08125000000000004,
          -1.3,
          null,
          -1.3,
          -1.3,
          null,
          -2.6,
          -2.6,
          null,
          -2.6,
          -2.6,
          null,
          -2.6,
          -3.9000000000000004,
          null,
          -3.9000000000000004,
          -3.9000000000000004,
          null
         ]
        },
        {
         "hoverinfo": "text",
         "marker": {
          "color": "black",
          "size": 3,
          "symbol": "circle"
         },
         "mode": "markers",
         "showlegend": false,
         "text": [
          "root",
          "edge.3",
          "edge.2",
          "edge.1",
          "edge.0",
          "Chimpanzee",
          "Human",
          "Gorilla",
          "Orangutan",
          "Rhesus",
          "HowlerMon",
          "Galago"
         ],
         "type": "scatter",
         "x": [
          0,
          0.03139195169925202,
          0.08159484122097149,
          0.09617900440577909,
          0.0961878725442622,
          0.10959571362620046,
          0.12185690868725763,
          0.1057597283305971,
          0.11501597595635671,
          0.12464469563427283,
          0.1930734042426634,
          0.7165441919466353
         ],
         "y": [
          -2.6,
          -0.08125000000000004,
          1.1375,
          2.275,
          3.25,
          3.9000000000000004,
          2.6,
          1.3,
          0,
          -1.3,
          -2.6,
          -3.9000000000000004
         ]
        }
       ],
       "layout": {
        "annotations": [
         {
          "font": {
           "family": "Inconsolata, monospace",
           "size": 12
          },
          "showarrow": false,
          "text": "Chimpanzee",
          "textangle": 0,
          "x": 0.12750931842486635,
          "xanchor": "left",
          "xref": "x",
          "y": 3.9000000000000004,
          "yref": "y"
         },
         {
          "font": {
           "family": "Inconsolata, monospace",
           "size": 12
          },
          "showarrow": false,
          "text": "Human",
          "textangle": 0,
          "x": 0.1397705134859235,
          "xanchor": "left",
          "xref": "x",
          "y": 2.6,
          "yref": "y"
         },
         {
          "font": {
           "family": "Inconsolata, monospace",
           "size": 12
          },
          "showarrow": false,
          "text": "Gorilla",
          "textangle": 0,
          "x": 0.12367333312926299,
          "xanchor": "left",
          "xref": "x",
          "y": 1.3,
          "yref": "y"
         },
         {
          "font": {
           "family": "Inconsolata, monospace",
           "size": 12
          },
          "showarrow": false,
          "text": "Orangutan",
          "textangle": 0,
          "x": 0.1329295807550226,
          "xanchor": "left",
          "xref": "x",
          "y": 0,
          "yref": "y"
         },
         {
          "font": {
           "family": "Inconsolata, monospace",
           "size": 12
          },
          "showarrow": false,
          "text": "Rhesus",
          "textangle": 0,
          "x": 0.1425583004329387,
          "xanchor": "left",
          "xref": "x",
          "y": -1.3,
          "yref": "y"
         },
         {
          "font": {
           "family": "Inconsolata, monospace",
           "size": 12
          },
          "showarrow": false,
          "text": "HowlerMon",
          "textangle": 0,
          "x": 0.21098700904132928,
          "xanchor": "left",
          "xref": "x",
          "y": -2.6,
          "yref": "y"
         },
         {
          "font": {
           "family": "Inconsolata, monospace",
           "size": 12
          },
          "showarrow": false,
          "text": "Galago",
          "textangle": 0,
          "x": 0.7344577967453012,
          "xanchor": "left",
          "xref": "x",
          "y": -3.9000000000000004,
          "yref": "y"
         },
         {
          "ax": 0,
          "ay": 0,
          "showarrow": false,
          "text": "0.07",
          "x": 0.8240258207386306,
          "xref": "x",
          "y": 3.9000000000000004,
          "yref": "y"
         }
        ],
        "autosize": true,
        "font": {
         "family": "Balto",
         "size": 14
        },
        "height": 500,
        "hovermode": "closest",
        "margin": {
         "b": 50,
         "l": 50,
         "pad": 4,
         "r": 50,
         "t": 50
        },
        "shapes": [
         {
          "line": {
           "color": "black",
           "width": 1.25
          },
          "type": "line",
          "x0": 0.7165441919466353,
          "x1": 0.7881986111412989,
          "y0": 3.9000000000000004,
          "y1": 3.9000000000000004
         }
        ],
        "showlegend": false,
        "template": {
         "data": {
          "bar": [
           {
            "error_x": {
             "color": "#2a3f5f"
            },
            "error_y": {
             "color": "#2a3f5f"
            },
            "marker": {
             "line": {
              "color": "#E5ECF6",
              "width": 0.5
             }
            },
            "type": "bar"
           }
          ],
          "barpolar": [
           {
            "marker": {
             "line": {
              "color": "#E5ECF6",
              "width": 0.5
             }
            },
            "type": "barpolar"
           }
          ],
          "carpet": [
           {
            "aaxis": {
             "endlinecolor": "#2a3f5f",
             "gridcolor": "white",
             "linecolor": "white",
             "minorgridcolor": "white",
             "startlinecolor": "#2a3f5f"
            },
            "baxis": {
             "endlinecolor": "#2a3f5f",
             "gridcolor": "white",
             "linecolor": "white",
             "minorgridcolor": "white",
             "startlinecolor": "#2a3f5f"
            },
            "type": "carpet"
           }
          ],
          "choropleth": [
           {
            "colorbar": {
             "outlinewidth": 0,
             "ticks": ""
            },
            "type": "choropleth"
           }
          ],
          "contour": [
           {
            "colorbar": {
             "outlinewidth": 0,
             "ticks": ""
            },
            "colorscale": [
             [
              0,
              "#0d0887"
             ],
             [
              0.1111111111111111,
              "#46039f"
             ],
             [
              0.2222222222222222,
              "#7201a8"
             ],
             [
              0.3333333333333333,
              "#9c179e"
             ],
             [
              0.4444444444444444,
              "#bd3786"
             ],
             [
              0.5555555555555556,
              "#d8576b"
             ],
             [
              0.6666666666666666,
              "#ed7953"
             ],
             [
              0.7777777777777778,
              "#fb9f3a"
             ],
             [
              0.8888888888888888,
              "#fdca26"
             ],
             [
              1,
              "#f0f921"
             ]
            ],
            "type": "contour"
           }
          ],
          "contourcarpet": [
           {
            "colorbar": {
             "outlinewidth": 0,
             "ticks": ""
            },
            "type": "contourcarpet"
           }
          ],
          "heatmap": [
           {
            "colorbar": {
             "outlinewidth": 0,
             "ticks": ""
            },
            "colorscale": [
             [
              0,
              "#0d0887"
             ],
             [
              0.1111111111111111,
              "#46039f"
             ],
             [
              0.2222222222222222,
              "#7201a8"
             ],
             [
              0.3333333333333333,
              "#9c179e"
             ],
             [
              0.4444444444444444,
              "#bd3786"
             ],
             [
              0.5555555555555556,
              "#d8576b"
             ],
             [
              0.6666666666666666,
              "#ed7953"
             ],
             [
              0.7777777777777778,
              "#fb9f3a"
             ],
             [
              0.8888888888888888,
              "#fdca26"
             ],
             [
              1,
              "#f0f921"
             ]
            ],
            "type": "heatmap"
           }
          ],
          "heatmapgl": [
           {
            "colorbar": {
             "outlinewidth": 0,
             "ticks": ""
            },
            "colorscale": [
             [
              0,
              "#0d0887"
             ],
             [
              0.1111111111111111,
              "#46039f"
             ],
             [
              0.2222222222222222,
              "#7201a8"
             ],
             [
              0.3333333333333333,
              "#9c179e"
             ],
             [
              0.4444444444444444,
              "#bd3786"
             ],
             [
              0.5555555555555556,
              "#d8576b"
             ],
             [
              0.6666666666666666,
              "#ed7953"
             ],
             [
              0.7777777777777778,
              "#fb9f3a"
             ],
             [
              0.8888888888888888,
              "#fdca26"
             ],
             [
              1,
              "#f0f921"
             ]
            ],
            "type": "heatmapgl"
           }
          ],
          "histogram": [
           {
            "marker": {
             "colorbar": {
              "outlinewidth": 0,
              "ticks": ""
             }
            },
            "type": "histogram"
           }
          ],
          "histogram2d": [
           {
            "colorbar": {
             "outlinewidth": 0,
             "ticks": ""
            },
            "colorscale": [
             [
              0,
              "#0d0887"
             ],
             [
              0.1111111111111111,
              "#46039f"
             ],
             [
              0.2222222222222222,
              "#7201a8"
             ],
             [
              0.3333333333333333,
              "#9c179e"
             ],
             [
              0.4444444444444444,
              "#bd3786"
             ],
             [
              0.5555555555555556,
              "#d8576b"
             ],
             [
              0.6666666666666666,
              "#ed7953"
             ],
             [
              0.7777777777777778,
              "#fb9f3a"
             ],
             [
              0.8888888888888888,
              "#fdca26"
             ],
             [
              1,
              "#f0f921"
             ]
            ],
            "type": "histogram2d"
           }
          ],
          "histogram2dcontour": [
           {
            "colorbar": {
             "outlinewidth": 0,
             "ticks": ""
            },
            "colorscale": [
             [
              0,
              "#0d0887"
             ],
             [
              0.1111111111111111,
              "#46039f"
             ],
             [
              0.2222222222222222,
              "#7201a8"
             ],
             [
              0.3333333333333333,
              "#9c179e"
             ],
             [
              0.4444444444444444,
              "#bd3786"
             ],
             [
              0.5555555555555556,
              "#d8576b"
             ],
             [
              0.6666666666666666,
              "#ed7953"
             ],
             [
              0.7777777777777778,
              "#fb9f3a"
             ],
             [
              0.8888888888888888,
              "#fdca26"
             ],
             [
              1,
              "#f0f921"
             ]
            ],
            "type": "histogram2dcontour"
           }
          ],
          "mesh3d": [
           {
            "colorbar": {
             "outlinewidth": 0,
             "ticks": ""
            },
            "type": "mesh3d"
           }
          ],
          "parcoords": [
           {
            "line": {
             "colorbar": {
              "outlinewidth": 0,
              "ticks": ""
             }
            },
            "type": "parcoords"
           }
          ],
          "scatter": [
           {
            "marker": {
             "colorbar": {
              "outlinewidth": 0,
              "ticks": ""
             }
            },
            "type": "scatter"
           }
          ],
          "scatter3d": [
           {
            "line": {
             "colorbar": {
              "outlinewidth": 0,
              "ticks": ""
             }
            },
            "marker": {
             "colorbar": {
              "outlinewidth": 0,
              "ticks": ""
             }
            },
            "type": "scatter3d"
           }
          ],
          "scattercarpet": [
           {
            "marker": {
             "colorbar": {
              "outlinewidth": 0,
              "ticks": ""
             }
            },
            "type": "scattercarpet"
           }
          ],
          "scattergeo": [
           {
            "marker": {
             "colorbar": {
              "outlinewidth": 0,
              "ticks": ""
             }
            },
            "type": "scattergeo"
           }
          ],
          "scattergl": [
           {
            "marker": {
             "colorbar": {
              "outlinewidth": 0,
              "ticks": ""
             }
            },
            "type": "scattergl"
           }
          ],
          "scattermapbox": [
           {
            "marker": {
             "colorbar": {
              "outlinewidth": 0,
              "ticks": ""
             }
            },
            "type": "scattermapbox"
           }
          ],
          "scatterpolar": [
           {
            "marker": {
             "colorbar": {
              "outlinewidth": 0,
              "ticks": ""
             }
            },
            "type": "scatterpolar"
           }
          ],
          "scatterpolargl": [
           {
            "marker": {
             "colorbar": {
              "outlinewidth": 0,
              "ticks": ""
             }
            },
            "type": "scatterpolargl"
           }
          ],
          "scatterternary": [
           {
            "marker": {
             "colorbar": {
              "outlinewidth": 0,
              "ticks": ""
             }
            },
            "type": "scatterternary"
           }
          ],
          "surface": [
           {
            "colorbar": {
             "outlinewidth": 0,
             "ticks": ""
            },
            "colorscale": [
             [
              0,
              "#0d0887"
             ],
             [
              0.1111111111111111,
              "#46039f"
             ],
             [
              0.2222222222222222,
              "#7201a8"
             ],
             [
              0.3333333333333333,
              "#9c179e"
             ],
             [
              0.4444444444444444,
              "#bd3786"
             ],
             [
              0.5555555555555556,
              "#d8576b"
             ],
             [
              0.6666666666666666,
              "#ed7953"
             ],
             [
              0.7777777777777778,
              "#fb9f3a"
             ],
             [
              0.8888888888888888,
              "#fdca26"
             ],
             [
              1,
              "#f0f921"
             ]
            ],
            "type": "surface"
           }
          ],
          "table": [
           {
            "cells": {
             "fill": {
              "color": "#EBF0F8"
             },
             "line": {
              "color": "white"
             }
            },
            "header": {
             "fill": {
              "color": "#C8D4E3"
             },
             "line": {
              "color": "white"
             }
            },
            "type": "table"
           }
          ]
         },
         "layout": {
          "annotationdefaults": {
           "arrowcolor": "#2a3f5f",
           "arrowhead": 0,
           "arrowwidth": 1
          },
          "colorscale": {
           "diverging": [
            [
             0,
             "#8e0152"
            ],
            [
             0.1,
             "#c51b7d"
            ],
            [
             0.2,
             "#de77ae"
            ],
            [
             0.3,
             "#f1b6da"
            ],
            [
             0.4,
             "#fde0ef"
            ],
            [
             0.5,
             "#f7f7f7"
            ],
            [
             0.6,
             "#e6f5d0"
            ],
            [
             0.7,
             "#b8e186"
            ],
            [
             0.8,
             "#7fbc41"
            ],
            [
             0.9,
             "#4d9221"
            ],
            [
             1,
             "#276419"
            ]
           ],
           "sequential": [
            [
             0,
             "#0d0887"
            ],
            [
             0.1111111111111111,
             "#46039f"
            ],
            [
             0.2222222222222222,
             "#7201a8"
            ],
            [
             0.3333333333333333,
             "#9c179e"
            ],
            [
             0.4444444444444444,
             "#bd3786"
            ],
            [
             0.5555555555555556,
             "#d8576b"
            ],
            [
             0.6666666666666666,
             "#ed7953"
            ],
            [
             0.7777777777777778,
             "#fb9f3a"
            ],
            [
             0.8888888888888888,
             "#fdca26"
            ],
            [
             1,
             "#f0f921"
            ]
           ],
           "sequentialminus": [
            [
             0,
             "#0d0887"
            ],
            [
             0.1111111111111111,
             "#46039f"
            ],
            [
             0.2222222222222222,
             "#7201a8"
            ],
            [
             0.3333333333333333,
             "#9c179e"
            ],
            [
             0.4444444444444444,
             "#bd3786"
            ],
            [
             0.5555555555555556,
             "#d8576b"
            ],
            [
             0.6666666666666666,
             "#ed7953"
            ],
            [
             0.7777777777777778,
             "#fb9f3a"
            ],
            [
             0.8888888888888888,
             "#fdca26"
            ],
            [
             1,
             "#f0f921"
            ]
           ]
          },
          "colorway": [
           "#636efa",
           "#EF553B",
           "#00cc96",
           "#ab63fa",
           "#FFA15A",
           "#19d3f3",
           "#FF6692",
           "#B6E880",
           "#FF97FF",
           "#FECB52"
          ],
          "font": {
           "color": "#2a3f5f"
          },
          "geo": {
           "bgcolor": "white",
           "lakecolor": "white",
           "landcolor": "#E5ECF6",
           "showlakes": true,
           "showland": true,
           "subunitcolor": "white"
          },
          "hoverlabel": {
           "align": "left"
          },
          "hovermode": "closest",
          "mapbox": {
           "style": "light"
          },
          "paper_bgcolor": "white",
          "plot_bgcolor": "#E5ECF6",
          "polar": {
           "angularaxis": {
            "gridcolor": "white",
            "linecolor": "white",
            "ticks": ""
           },
           "bgcolor": "#E5ECF6",
           "radialaxis": {
            "gridcolor": "white",
            "linecolor": "white",
            "ticks": ""
           }
          },
          "scene": {
           "xaxis": {
            "backgroundcolor": "#E5ECF6",
            "gridcolor": "white",
            "gridwidth": 2,
            "linecolor": "white",
            "showbackground": true,
            "ticks": "",
            "zerolinecolor": "white"
           },
           "yaxis": {
            "backgroundcolor": "#E5ECF6",
            "gridcolor": "white",
            "gridwidth": 2,
            "linecolor": "white",
            "showbackground": true,
            "ticks": "",
            "zerolinecolor": "white"
           },
           "zaxis": {
            "backgroundcolor": "#E5ECF6",
            "gridcolor": "white",
            "gridwidth": 2,
            "linecolor": "white",
            "showbackground": true,
            "ticks": "",
            "zerolinecolor": "white"
           }
          },
          "shapedefaults": {
           "line": {
            "color": "#2a3f5f"
           }
          },
          "ternary": {
           "aaxis": {
            "gridcolor": "white",
            "linecolor": "white",
            "ticks": ""
           },
           "baxis": {
            "gridcolor": "white",
            "linecolor": "white",
            "ticks": ""
           },
           "bgcolor": "#E5ECF6",
           "caxis": {
            "gridcolor": "white",
            "linecolor": "white",
            "ticks": ""
           }
          },
          "title": {
           "x": 0.05
          },
          "xaxis": {
           "automargin": true,
           "gridcolor": "white",
           "linecolor": "white",
           "ticks": "",
           "zerolinecolor": "white",
           "zerolinewidth": 2
          },
          "yaxis": {
           "automargin": true,
           "gridcolor": "white",
           "linecolor": "white",
           "ticks": "",
           "zerolinecolor": "white",
           "zerolinewidth": 2
          }
         }
        },
        "width": 500,
        "xaxis": {
         "autorange": true,
         "range": [
          -0.052283017591692425,
          0.8569868535681758
         ],
         "type": "linear",
         "visible": false
        },
        "yaxis": {
         "autorange": true,
         "range": [
          -4.40677966101695,
          4.40677966101695
         ],
         "type": "linear",
         "visible": false
        }
       }
      },
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA2kAAAH0CAYAAACn72USAAAgAElEQVR4Xu3dd9zVdd348TcyFQTEkXvPHGWWZppGag7U1MSZuDcoYLjCxC2ucOSqG8edmgNNw0xw5SKJVDJTc6GUG1cqsn+P8/UHt+NC5nn7/XKe1z939+W5vu/Peb6/Dx++OOc6NJsyZcqU8EWAAAECBAgQIECAAAECpRBoJtJKsQeHIECAAAECBAgQIECAQCEg0twIBAgQIECAAAECBAgQKJGASCvRMhyFAAECBAgQIECAAAECIs09QIAAAQIECBAgQIAAgRIJiLQSLcNRCBAgQIAAAQIECBAgINLcAwQIECBAgAABAgQIECiRgEgr0TIchQABAgQIECBAgAABAiLNPUCAAAECBAgQIECAAIESCYi0Ei3DUQgQIECAAAECBAgQICDS3AMECBAgQIAAAQIECBAokYBIK9EyHIUAAQIECBAgQIAAAQIizT1AgAABAgQIECBAgACBEgmItBItw1EIECBAgAABAgQIECAg0twDBAgQIECAAAECBAgQKJGASCvRMhyFAAECBAgQIECAAAECIs09QIAAAQIECBAgQIAAgRIJiLQSLcNRCBAgQIAAAQIECBAgINLcAwQIECBAgAABAgQIECiRgEgr0TIchQABAgQIECBAgAABAiLNPUCAAAECBAgQIECAAIESCYi0Ei3DUQgQIECAAAECBAgQICDS3AMECBAgQIAAAQIECBAokYBIK9EyHIUAAQIECBAgQIAAAQIizT1AgAABAgQIECBAgACBEgmItBItw1EIECBAgAABAgQIECAg0twDBAgQIECAAAECBAgQKJGASCvRMhyFAAECBAgQIECAAAECIs09QIAAAQIECBAgQIAAgRIJiLQSLcNRCBAgQIAAAQIECBAgINLcAwQIECBAgAABAgQIECiRgEgr0TIchQABAgQIECBAgAABAiLNPUCAAAECBAgQIECAAIESCYi0Ei3DUQgQIECAAAECBAgQICDS3AMECBAgQIAAAQIECBAokYBIK9EyHIUAAQIECBAgQIAAAQIizT1AgAABAgQIECBAgACBEgmItBItw1EIECBAgAABAgQIECAg0twDBAgQIECAAAECBAgQKJGASCvRMhyFAAECBAgQIECAAAECIs09QIAAAQIECBAgQIAAgRIJiLQSLcNRCBAgQIAAAQIECBAgINLcAwQIECBAgAABAgQIECiRgEgr0TIchQABAgQIECBAgAABAiLNPUCAAAECBAgQIECAAIESCYi0Ei3DUQgQIECAAAECBAgQIFDZSHtlzFjbI0CAAAECBAgQIECAwHQFllx4/krqiLRKrs2hCRAgQIAAAQIECBCYkYBIm5HQXP7nXkmby6AuR4AAAQIECBAgQGAeExBpyQsVacngxhEgQIAAAQIECBComIBIS16YSEsGN44AAQIECBAgQIBAxQREWvLCRFoyuHEECBAgQIAAAQIEKiYg0pIXJtKSwY0jQIAAAQIECBAgUDEBkZa8MJGWDG4cAQIECBAgQIAAgYoJiLTkhYm0ZHDjCBAgQIAAAQIECFRMQKQlL0ykJYMbR4AAAQIECBAgQKBiAiIteWEiLRncOAIECBAgQIAAAQIVExBpyQsTacngxhEgQIAAAQIECBComIBIS16YSEsGN44AAQIECBAgQIBAxQREWvLCRFoyuHEECBAgQIAAAQIEKiYg0pIXJtKSwY0jQIAAAQIECBAgUDEBkZa8MJGWDG4cAQIECBAgQIAAgYoJiLTkhYm0ZHDjCBAgQIAAAQIECFRMQKQlL0ykJYMbR4AAAQIECBAgQKBiAiIteWEiLRncOAIECBAgQIAAAQIVExBpyQsTacngxhEgQIAAAQIECBComIBIS16YSEsGN44AAQIECBAgQIBAxQREWvLCRFoyuHEECBAgQIAAAQIE5kBgypQp8dHYj6PtAvPP8CpjPx4XrVq2jObN55vhY7/sASJtjvhm/YdF2qyb+Ym5L1D7l02zZs2KC7/3/gfxj6eei402+GaTg8a8/V5MmDghFl9skbl/EFckQIAAAQIECJRY4O77h8ddf34kVlp+6Xh+1L/j2CP3jQ7t233hxLWIO2PAFbHUEovGf159Mzb41pqx7ZabFI/rf8GVxTU+/bXGqivEBWccPd1nLtKSbwqRlgxu3GcEnnjqubjo17+L8RMmxrJLLx59Du8W/3ntzbjkihtjwGk/a1Lr+luGxJtj3onuB+xKkwABAgQIECDQMAIvvPSfOLj3qXH1xafEEl9bJK763R/iX8+/FKf9vPsXDM4YMLB4Be2ow/eKDz74KHbc+6giwpZZ8mvR76zLYtcdfxQLtmtb/NzNg++OZZdeIn7adRuRVpa7SaSVZRONd44PPhwb+/Y4Mc49uXcRaLV/0Swwf5tYZ81Vp0Xap19hm5HQ1MfOys9M75qzeo3pPX5WrzOj5+ifEyBAgAABAo0rcP3vh8RDjzw+7RWvF0b9Ow4+6rS47ZoBMX+b1tNgJk6aFFvv0j3O/MURsd431ii+XwuzFZZdMrb90Sevpi3cqcO0xx/Q8+Q45oh9YpUVlxVpZbm9RFpZNlGNc1z6qwHRrNl8cfBhR8zxge9/+NH4w5D74+x+PYtrTZo0uXi/9DPPvRQX/vq6WHmFZeL2oQ/Gj7f+wbRXzU499zfFnxjtsE3n2GnbHxY/97833B5/fezJeGn0q7HN5hvHXfc/Eofvv0vxp0N33vNwvPPef+PxJ56Jg7rtFDtvv3nxtsprB90Rg+98ICZMnBiH7dc1Om/8nXj2hZej9ipdLRTvuPuh4vqH7ts1Ro1+NX5xxsXFrIkTJ0WLFs2Lf+ktufiiMfzRJ+PG24bGqJdfiX123z66bLFx8bjpfX+O0VyAAAECBAgQaFiB2n8HTYkpccJRBxYGb7/7fuy6/zHFu4/WXH2laS61/y45sNcpcfl5fWOF5ZYqvv+r/7khRv/nteK/YT799cprb8aRx58dN/xP/2m/etIUsLc7Jt92Ii0ZvITjevU4aKZO9fRTT8bfH3+seGy/U8+KAw/54kvrM3Wh//+g6wb9qfgXzR4/2fozP1aLtO7HnBkn9jko1v76KlH7052LzjwmvrbYwlF7ZeqW2++NDz78KLrtum3xcxf95vro2GHBIp6efPqF2LFL5yLOvv/ddePsi66Oy395QkyZMjl69z0vTvjZgVF7z/XDw0cWf7JUm39Qr1Pj3FN6xxtvvRM9jz87Tj720Fh9leVjn+4nxsAL+sVCHReMj8eNL2b97/W3x6tvvFWc7dXX34pDf3Z6XHFhv+KfdT+mf1x45jExbvz4Jr//6T+xmhUnjyVAgAABAvOCwMz+98a88Fw//Rw222Lr2Hb7HefK0+p22Amx/rfWmvaH1xMnToytd+1R/IHzjl0++cPr2te9D/41Tv/lwPjdr8+c9orZb2+8Pa763eAYctPFn4mx3//x3vj3K2/M8NdIRNpcWeHMX0SkzbzVvPjIu4bcEft32zV26rrbDJ/e3I602itgtVfOmoq0M8+/Ylr8nHDGJbHbTlvGmqutWJzxltvvidpbJffapcu0SKv9C2vM2+/GmHfeK+Kr9naAzTdZP4b++ZE46ZhDPgmsG26Pli1bxG47bhnvvPt+8Ypd7Xfbbh58T/Q8ZI9o0aJFXHDZtXHZeX2Lxx93yoWx927bFcFW+6q9Wnfhr38XF599fLRrO3/xC7dXXz84frDRt4t/ft9DI+KAvXaM8eMnNPn9TTb81gyNPYAAAQIECMyrAiJtzjfbdb+j44ff/07xTp/aV+0Pr3+082FFYNXeeTT1a8i9w4o/qL5x4FnFH2TXvmrvFvrNb2/5QqT17ntu7Nl1m2lvi5zeKUXanO9vlq4g0maJa557cC3SDjtw7/jXS2/M1HObm293rP0LZMi9f4lzTu71mdm1ePr0B4ec2P/S4m2Ka6+xcvG4piJtg/XWirfGvBvvvPd+rLv26kWkbbHpBnH/sEfjuJ77FT930213xfv//bD4l1iPY/sXr7jV3rI46A93x167dolWrVrFldfdNu3tlz8/7aL4adcuxStvr78xJrof2z/OOKFH8TbM2tfQ+/4SDw0fWbz1cupX7Xfr/vb4P5v8fqeO7WfK2IMIECBAgAABAk0JHHHcWcUHhkz9b5tx48bHtnscGaf37R7fWXfNaT/yt5FPxbEnXxD/c/6Jxe/9175qv/tf+wPni/ofO+1xtU/U3nnfPjH4uguidauWX4ou0pLvSZGWDF6ycbMaaXPz+O++99/Y+/BfxAVnHh3LLb1E8Yuw7RdsW8TS3Iq0y64aFJf/sm80n6959Ok3IHbebrMYN35CjHj8n8UvyNY+7eig3qdGnx7dphtpKy6/dPT6+Tmx3VabxNabbTSN4LU33oqj+51fvLLWdoE28cxzo2K1lZeP198c0+T3p/4VA3PT0LUIECBAgACBxhG49Mqb4l/PvRTnnXpU8aRHv/J67NejX/z2klOLXwuZ+vX2O+/FrgccG2ef1Cu+udaqxbdr71Jq3bpV9Dpkz2mPq70raNiIv0/7HbcvkxRpyfeZSEsGL9m4rzLSahS1V6IGXHpN1F5lqv2FjLXfGau9BfHSK26a9i+g2itpXbffPNZaY+Xik4mefPr54i2FtVfPju+1f/E7aRt+Z51446234913/xvfXHu1uPHWobHZJusXbzvstFCHeO7F0bH5JhvEgd12LCKt9kpax/YLFm8BGDX6leKDRtq0bh0Dr7112itpfU//Vey58zbF30Fy/mXXFn9yNfVr1x23LD4k5I9DH4xaCNb+fpLanFOOO7T4wJLpfb9k63ccAgQIECBAoEICL//7tdj/yJOmfQT/4Dvvj6eeHRV9uncrnkXtA0MefOSx4r9fan8XWu0j+Hsf9tMYP2FC7HP4idG/35HFR/BP/ar9d1Xtd/hr/800oy+RNiOhufzPRdpcBq3Y5b7qSKtx1T4m9uOx46JduwXmml7tPdoP//Xv8cCwR6NP973j43Hjigj89Fft99pqv1s2o6/atZr6mvrKWC36Joyf8IXzT+/7M5rnnxMgQIAAAQIEpidw74Mj4tY77ot1114tXntjTPHp1Qv9/1+pqP1qR+0Pj3//v+dF7T9fzjx/YCy6SKfiUquutOxn3hFU+2+j7fboGTddcXaTfxn25+eLtOR7UqQlg5dsXBkirV4ktVfpapF27JH71muE6xIgQIAAAQIE0gUmT55cfPJ07a8N+vRX7Q+WP6z9IfSn/uD7w4/GFo/7/K9d1P7qo9qvnszsp0+LtOQ1i7Rk8JKNm5cjrWTUjkOAAAECBAgQqKyASEtenUhLBi/ZOJFWsoU4DgECBAgQIECghAIiLXkpIi0ZvGTjRFrJFuI4BAgQIECAAIESCoi05KWItGTwko0TaSVbiOMQIECAAAECBEooINKSlyLSksFLNk6klWwhjkOAAAECBAgQKKGASEteikhLBi/ZOJFWsoU4DgECBAgQIECghAIiLXkpIi0ZvGTjRFrJFuI4BAgQIECAAIESCoi05KWItGTwko0TaSVbiOMQIECAAAECBEooINKSlyLSksFLNk6klWwhjkOAAAECBAgQKKGASEteikhLBi/ZOJFWsoU4DgECBAgQIECghAIiLXkpIi0ZvGTjRFrJFuI4BAgQIECAAIESCoi05KWItGTwko0TaSVbiOMQIECAAAECBEooINKSlyLSksFLNk6klWwhjkOAAAECBAgQKKGASEteikhLBi/ZOJFWsoU4DgECBAgQIECghAIiLXkpIi0ZvGTjRFrJFuI4BAgQIECAAIESCoi05KWItGTwko0TaSVbiOMQIECAAAECBEooINKSlyLSksFLNk6klWwhjkOAAAECBAgQKKGASEteikhLBi/ZOJFWsoU4DgECBAgQIECghAIiLXkpIi0ZvGTjRFrJFuI4BAgQIECAAIESCoi05KWItGTwko0TaSVbiOMQIECAAAECBEooINKSlyLSksFLNk6klWwhjkOAAAECBAgQKKGASEteikhLBi/ZOJFWsoU4DgECBAgQIECghAIiLXkpIi0ZvGTjvr/BOvHC88/FIYcfGSecdEbJTuc4BAgQIECAAAECZRAQaclbEGnJ4CUa9+QTI+NHnTcsTrTaGl+Pex4YUaLTOQoBAgQIECBAgEBZBERa8iZEWjJ4ycatt9ZK8dprr3olrWR7cRwCBAgQIECAQJkERFryNkRaMnjJxp19xslx/nn9499vfliykzkOAQIECBAgQIBAWQREWvImRFoyeMnGibSSLcRxCBAgQIAAAQIlFBBpyUsRacngJRsn0kq2EMchQIAAAQIECJRQQKQlL0WkJYOXbJxIK9lCHIcAAQIECBAgUEIBkZa8FJGWDF6ycSKtZAtxHAIECBAgQIBACQVEWvJSRFoyeMnGibSSLcRxCBAgQIAAAQIlFBBpyUsRacngJRsn0kq2EMchQIAAAQIECJRQQKQlL0WkJYOXbJxIK9lCHIcAAQIECBAgUEIBkZa8FJGWDF6ycSKtZAtxHAIECBAgQIBACQVEWvJSRFoyeMnGibSSLcRxCBAgQIAAAQIlFBBpyUsRacngJRsn0kq2EMchQIAAAQIECJRQQKQlL0WkJYOXbJxIK9lCHIcAAQIECBAgUEIBkZa8FJGWDF6ycSKtZAtxHAIECBAgQIBACQVEWvJSRFoyeMnGibSSLcRxCBAgQIAAAQIlFBBpyUsRacngJRsn0kq2EMchQIAAAQIECJRQQKQlL0WkJYOXbJxIK9lCHIcAAQIECBAgUEIBkZa8FJGWDF6ycSKtZAtxHAIECBAgQIBACQVEWvJSRFoyeMnGibSSLcRxCBAgQIAAAQIlFBBpyUsRacngJRsn0kq2EMchQIAAAQIECJRQQKQlL0WkJYOXbJxIK9lCHIcAAQIECBAgUEIBkZa8FJGWDF6ycSKtZAtxHAIECBAgQIBACQVEWvJSRFoyeMnGibSSLcRxCBAgQIAAAQIlFBBpyUsRacngJRsn0kq2EMchQIAAAQIECJRQQKQlL0WkJYOXbJxIK9lCHIcAAQIECBAgUEIBkZa8FJGWDF6ycVWItCH3DovnXhwdh+23S6F38cAbYo1VV4zOG3+7ZJqOQ4AAAQIECBCYNwVEWvJeRVoyeMnGVSHS7rxnWLz48n/ikH12LvQu+s31sc6aq8QmG37rSzWnTJkSzZo1i6n/d2boZ+WxM3M9jyFAgAABAgQIzAsCIi15iyItGbxk4+oVaf94YmQM/dPtscJKK8cOO33yCtjsfk0v0jq2XzAe+Mtjcfj+n1z//Muviy07fzf++tg/46+PPRkvjX41ttl847jr/keKx/xgo2/HtYPuiMF3PhATJk6Mw/brGp03/k48+8LLcf0tQ2KB+dvEHXc/FDtt+8M4dN+us3tcP0eAAAECBAgQmOcERFrySkVaMnjJxtUibcC5Z8Yuu/90rp7sriF/irfHvFVc8+ER/4jlll9xtq9fi7SbB98d66y5anGNx554Orrtum20ad067rz34TjhqAOL75/Y/9LYYZsfxEOPjIyOHRaMFi2ax5NPvxA7dukcd97zcBzXc794ePjIWO8ba8SUmBIH9To1zj2ld7zx1jvR8/iz4+RjD43VV1k+9ul+Ygy8oF8s3KnDbJ/ZDxIgQIAAAQIE5iUBkZa8TZGWDF6ycUP+dHsctO8esePOu87Vk83tSBs24u/x4603Lc446A93x486bxitW7WKIfcNazLS1v/WWjHm7XdjzDvvFVF2/e+HRL+jD4533n0/nnnupXhzzDtx8+B7ouche0SLFi3igsuujcvO61tc/7hTLoy9d9uuCDZfBAgQIECAAAECESIt+S4QacngJRtXi7QjDt0vnn7x9bl6snq+3bH2tsZ1114t2i4wf/z+j/fFKccdWpy9T78BscdPtipeSdtgvbXirTHvxjvvvR/rrr16EWk9Dtgtehzbv3hlbcnFFy1ib69du0SrVq3iyutui7P79Syu8/PTLoqfdu0Sa6y6wlw1cTECBAgQIECAQFUFRFry5kRaMnjJxtUr0ubm0/z876RNjbQVll0yep9wXgy84MSYOHFS8TbFfsccPN1I+/53140Rj/8zjjlin/jgg4/ioN6nRp8e3UTa3FyWaxEgQIAAAQLzpIBIS16rSEsGL9m4KkRa7SP4X3jp/z7dcWqk1aLrrAuvKn7PbKklFouxH48rPgzkkb/9Izb8zjrxxltvx7vv/je+ufZqceOtQ6NPj72LV9JqHzhS+521UaNfKT5QpPa7bQOvvXXaK2l9T/9V7LnzNl5JK9m96jgECBAgQIDAVycg0pLtRVoyeMnGVSHSPk/2+Y/Jf+/9D6L9gm2Lj9uf3tenf+aDD8dGu7bzl2wTjkOAAAECBAgQKK+ASEvejUhLBi/ZuCpGWskIHYcAAQIECBAgMM8LiLTkFYu0ZPCSjRNpJVuI4xAgQIAAAQIESigg0pKXItKSwUs2TqSVbCGOQ4AAAQIECBAooYBIS16KSEsGL9k4kVayhTgOAQIECBAgQKCEAiIteSkiLRm8ZONEWskW4jgECBAgQIAAgRIKiLTkpYi0ZPCSjRNpJVuI4xAgQIAAAQIESigg0pKXItKSwUs2TqSVbCGOQ4AAAQIECBAooYBIS16KSEsGL9k4kVayhTgOAQIECBAgQKCEAiIteSkiLRm8ZOO67b5T3HfP0Dj3/Euj6257lux0jkOAAAECBAgQIFAGAZGWvAWRlgxe53FnnPKLuOG6/43WbdrMcNKE8ePjtddeLR63xtfXirvuHz7Dn/EAAgQIECBAgACBxhMQack7F2nJ4HUet/SibeOwHr1j5VVXm+Gk999/L07qe0xMnjw5ttl2h/j1ldfO8Gc8gAABAgQIECBAoPEERFryzkVaMnidx9Ui7YZb7ojvbbzJTE265uqB0a/v0fHsy2/N1OM9iAABAgQIECBAoPEERFryzkVaMnidx81qpPngkDovxOUJECBAgAABAvOAgEhLXqJISwav8ziRVmdglydAgAABAgQINKCASEteukhLBq/zOJFWZ2CXJ0CAAAECBAg0oIBIS166SEsGr/M4kVZnYJcnQIAAAQIECDSggEhLXrpISwav8ziRVmdglydAgAABAgQINKCASEteukhLBq/zOJFWZ2CXJ0CAAAECBAg0oIBIS166SEsGr/M4kVZnYJcnQIAAAQIECDSggEhLXrpISwav8ziRVmdglydAgAABAgQINKCASEteukhLBq/zOJFWZ2CXJ0CAAAECBAg0oIBIS166SEsGr/M4kVZnYJcnQIAAAQIECDSggEhLXrpISwav8ziRVmdglydAgAABAgQINKCASEteukhLBq/zOJFWZ2CXJ0CAAAECBAg0oIBIS166SEsGr/M4kVZnYJcnQIAAAQIECDSggEhLXrpISwav8ziRVmdglydAgAABAgQINKCASEteukhLBq/zOJFWZ2CXJ0CAAAECBAg0oIBIS166SEsGr/M4kVZnYJcnQIAAAQIECDSggEhLXrpISwav8ziRVmdglydAgAABAgQINKCASEteukhLBq/zOJFWZ2CXJ0CAAAECBAg0oIBIS166SEsGr/M4kVZnYJcnQIAAAQIECDSggEhLXrpISwav8ziRVmdglydAgAABAgQINKCASEteukhLBq/zOJFWZ2CXJ0CAAAECBAg0oIBIS166SEsGr/M4kVZnYJcnQIAAAQIECDSggEhLXrpISwav8ziRVmdglydAgAABAgQINKCASEteukhLBq/zOJFWZ2CXJ0CAAAECBAg0oIBIS166SEsGr/M4kVZnYJcnQIAAAQIECDSggEhLXrpISwav8ziRVmdglydAgAABAgQINKCASEteukhLBq/zOJFWZ2CXJ0CAAAECBAg0oIBIS166SEsGr/M4kVZnYJcnQIAAAQIECDSggEhLXrpISwav8ziRVmdglydAgAABAgQINKCASEteukhLBq/zOJFWZ2CXJ0CAAAECBAg0oIBIS166SEsGr/M4kVZnYJcnQIAAAQIECDSggEhLXrpISwav8ziRVmdglydAgAABAgQINKCASEteukhLBq/zOJFWZ2CXJ0CAAAECBAg0oIBIS166SEsGr/M4kVZnYJcnQIAAAQIECDSggEhLXrpISwav8ziRVmdglydAgAABAgQINKCASEteukhLBq/zOJFWZ2CXJ0CAAAECBAg0oIBIS166SEsGr/M4kVZnYJcnQIAAAQIECDSggEhLXrpISwav87iqRNqkSZOjefP55khjzNvvxYSJE2LxxRYprlO75j0PDC/+90orLBMrLrfUl35/job7YQIECBAgQIBAAwmItORli7Rk8DqPK3ukPfr3p2PApdfE2I/HxWorLxfHHrFvtGu3wGypXH/LkHhzzDvR/YBdi5+fPHlyDBvxRDww7NFYbpklYvedtvrS78/WUD9EgAABAgQIEGhAAZGWvHSRlgxe53FljrS3330/uh99Zpzxix6x3NJLxJXX3RbLLbNkdN7423NVpRZvk6dMnhZpUy8+ve/P1eEuRoAAAQIECBCYBwVEWvJSRVoyeJ3Hze1Ie2LkYzF0yB2x0sqrxo933HmOTj/k3mHxr+dfnvbK15QpU6JZs2bFNe99cET8+uqb48OPxkbPQ/acFm59+g2I7663dtx6x33RvHnz6HXInrHOmqvEqef+Jv71/EuxwzadY6dtf/iZc81KpF076I4YfOcDMWHixDhsv67ReePvzNFz9MMECBAgQIAAgXlRQKQlb1WkJYPXedzsRNpB++4RO+78yVsGP/819M4/xjtvv118+y+PPhXLLLvcbD+D2itniy3SKbbZYuPPXGPU6FfjyOPOinNO7hUtmjePnj8/J84/4+hYfpklotthJ0SrVi3jnJN6xVtvvxvN55svVlhuqagF3i233xsffPhRdK+gNF4AACAASURBVNt129mOtIeHj4z1vrFGTIkpcVCvU+PcU3rHogsvNNvP0Q8SIECAAAECBOZFAZGWvFWRlgxe53GzGmm14/TqcdB0TzU3I+3SK2+KVVZcNjbbZP247MpB8dDwx2Orzb4XLVu0iFqo9enerTjH+ZdfFyssu2Rsv9WmRaTVvr/211f5whlvuf2e+ODDsbHXLl1mO9Leeff9eOa5l4rfbbt58D3R85A94htrrlrnLbk8AQIECBAgQKBaAiIteV8iLRm8zuNmJ9K+7Ehz8+2Otah6a8y7cWC3nYpPYvzLiL/HP556Pjp0aBevv/l2HHnQ7sVRBl57a7Rq2SJ+2rVLEWm1V9hqr8B9/mtOI6326ZA9ju0fO3bpHEsuvmgM+sPdsdeuXWLdtVev85ZcngABAgQIECBQLQGRlrwvkZYMXudxczvS5uZxX39jTBzU+9S49NyfxxJfWyRqv6P24kuvxMYbrhv9+l8al53XN9q0bhWHH31mHLpv11j/W2sWkXbuyb1j0UW++BbEOY20u+8fHiMe/2ccc8Q+8cEHHxVn69Ojm0ibm0t3LQIECBAgQGCeEBBpyWsUacngdR5X5kirPfXaB4RcPPCG6NC+XSy1xGKxy4+3iDVXXyl+d8udcd2gPxXfr70Fsvbx+bUPFalF2nmnHhWLdOr4Gbl+Z10WTz79fIwfPyE2WG+tOL7X/vHe+x/E0f0GxGtvjCl+j61Tx/bR/8Qji+s09f3WrVsVr6R1bL9gdOywYIwa/Uocvv8uIq3O96jLEyBAgAABAtUTEGnJOxNpyeB1Hlf2SKs9/YmTJhVxtcD8bT6jMW7c+IhmzaJ1q5bTvv/pT4CsF13t99ratZ2/Xpd3XQIECBAgQIBA5QVEWvIKRVoyeJ3H1SLt+pv/GBt9f9M6T3J5AgQIECBAgACBRhEQacmbFmlfDj527Ecx8rG/JW9l9sb98pwz4sH774slllwqHnns6eLvFfNFgAABAgQIECBAYE4FRNqcCs7iz4u0Lwc78rD946YbrptF1a/+4Q88MjJWXOmLH1v/1Z/MCQgQIECAAAECBKomINKSNybSvhz8kP33ik4LLxyn9f9l8mZmfdx5Z58e5511Wmy40SZx061/mvUL+AkCBAgQIECAAAECTQiItOTbQqTNXKSdftaA5M3M3rgqfHDI7D0zP0WAAAECBAgQIPBVCYi0ZHmRJtKSbznjCBAgQIAAAQIEKiYg0pIXJtJEWvItZxwBAgQIECBAgEDFBERa8sJEmkhLvuWMI0CAAAECBAgQqJiASEtemEgTacm3nHEECBAgQIAAAQIVExBpyQsTaSIt+ZYzjgABAgQIECBAoGICIi15YSJNpCXfcsYRIECAAAECBAhUTECkJS9MpIm05FvOOAIECBAgQIAAgYoJiLTkhYk0kZZ8yxlHgAABAgQIECBQMQGRlrwwkSbSkm854wgQIECAAAECBComINKSFybSRFryLWccAQIECBAgQIBAxQREWvLCRJpIS77ljCNAgAABAgQIEKiYgEhLXphIE2nJt5xxBAgQIECAAAECFRMQackLE2kiLfmWM44AAQIECBAgQKBiAiIteWEiTaQl33LGESBAgAABAgQIVExApCUvTKSJtORbzjgCBAgQIECAAIGKCYi05IWJNJGWfMsZR4AAAQIECBAgUDEBkZa8MJEm0pJvOeMIECBAgAABAgQqJiDSkhcm0kRa8i1nHAECBAgQIECAQMUERFrywkSaSEu+5YwjQIAAAQIECBComIBIS16YSBNpybeccQQIECBAgAABAhUTEGnJCxNpIi35ljOOAAECBAgQIECgYgIiLXlhIk2kJd9yxhEgQIAAAQIECFRMQKQlL0ykibTkW844AgQIECBAgACBigmItOSFiTSRlnzLGUeAAAECBAgQIFAxAZGWvDCRJtKSbznjCBAgQIAAAQIEKiYg0pIXJtJEWvItZxwBAgQIECBAgEDFBERa8sJEmkhLvuWMI0CAAAECBAgQqJiASEtemEgTacm3nHEECBAgQIAAAQIVExBpyQsTaSIt+ZYzjgABAgQIECBAoGICIi15YSJNpCXfcsYRIECAAAECBAhUTECkJS9MpIm05FvOOAIECBAgQIAAgYoJiLTkhYk0kZZ8yxlHgAABAgQIECBQMQGRlrwwkSbSkm854wgQIECAAAECBComINKSFybSRFryLWccAQIECBAgQIBAxQREWvLCRJpIS77ljCNAgAABAgQIEKiYgEhLXphIE2nJt5xxBAgQIECAAAECFRMQackLE2kiLfmWM44AAQIECBAgQKBiAiIteWEiTaQl33LGESBAgAABAgQIVExApCUvTKSJtORbzjgCBAgQIECAAIGKCYi05IWJNJGWfMsZR4AAAQIECBAgUDEBkZa8MJEm0pJvOeMIECBAgAABAgQqJiDSkhcm0kRa8i1nHAECBAgQIECAQMUERFrywkSaSEu+5YwjQIAAAQIECBComIBIS16YSBNpybeccQQIECBAgAABAhUTEGnJCxNpIi35lpvhuClTphSPadas2Qwf6wEECBAgQIAAAQL1FxBp9Tf+zASRJtKSb7kvHXfHXQ/FwGtvLR6zV9cusf3Wm37lx7t96IOx3jdWj8UXW+QrP4sDECBAgAABAgS+CgGRlqwu0kRa8i033XGjXn4lfnHmJXHhmcfE5MmT46Dep0a/ow+ONVdfabaPWHtV7vOvyE2aNDmaN5+vuGZtznzzffK/p/fV7bAT4tyTe8eiiyw00+doau5M/7AHEiBAgAABAgRKJiDSkhci0kTanNxyz/7r6Rh82y2xzDLLxc677jEnl4prbroj2i7QJnbYpnNxnetvGRL//eDDOGCvHaNPvwHx3fXWjlvvuC+aN28evQ7ZM9ZZc5W4dtAdMfjOB2LCxIlx2H5do/PG34lnX3i5+NkF5m8Td9z9UOy07Q/j0H27Trvm7UMfiEU6dYz27dvFwgt1iE2/t1488JfH4vD9dykec/7l18WWnb8bbdq0iV+ccXG8+vpbsejCC8WySy8eZ/7iiOIxszp3jmD8MAECBAgQIEDgKxYQackLEGnzXqRt2nnzWOxrX0u5k/58z13x+uuvFbOG/vmR+Pqaa8/23LMvujp++P3vxHrfWKO4Ri2c/vzQ36LvUQdE7dWsVq1axjkn9Yq33n43ms83X6yw3FLx8PCRxeOnxJQ4qNepce4pveONt96JnsefHScfe2isvsrysU/3E2PgBf1i/IQJcdp5v4kBp/eJ+x4cEX8b+VQcc8Q+MfzRJ+POex+OE446sJh7Yv9LY4dtfhDrrr161F5122qXw2PwtedHy5Ytpr3qNitzF+7UYbZN/CABAgQIECBAoAwCIi15CyJt3oq0G393TTz80J/T7qK5GWn9L7gyfrDRt2OD9dYqzj/yyX/FbX/6cxFPtUjr071brP31VT7z3N559/145rmX4s0x78TNg++JnofsES1atIgLLrs2Ljuvb/HY4065MPbebbuYNHlyXDLwhiLS7n3gr/Hci6OLV9ge+ds/Ysh9w5qMtNrbFn+082ExdNAlsz23Foq+CBAgQIAAAQJVFhBpydsTafNWpCXfPjE33+74vzfcHm1at4quP96ieBp/HPpgjP7P63HwPj8pIu2ck3vFYot0mvYUx7z9XvQ4tn/s2KVzLLn4ojHoD3fHXrt2iVatWsWV190WZ/frWTz256ddFD/t2qV4Va32+No1V1lx2fhZ972KDwOpvaL2+z/eF6ccd2jx+NpbK/f4yVbFK2njxo2Pbfc48jORNqtz11h1hey1mEeAAAECBAgQmKsCIm2ucs74YiJt+kYjH380frrLj6NV69Yx6LY7Y/kVZv8DLGa8CY94YdS/o+/pF8dF/Y8pQqsWVEcevEd8c61Vi0j7/Id33H3/8Bjx+D+Ltyx+8MFHxQeN9OnRbbqRtvSSi8XhR58Zp/XtHm0XmD8W6rBg8aEio//zWvQ+4bwYeMGJMXHipOLtkf2OObiItIkTJ8bWu/aIay8/vfhnC7ZrG4/87YlZmivS3NsECBAgQIBA1QVEWvIGqxhpJ/3i2LjmqoHRomWLump9PHZsjBs3rphxwkmnxyGHf/LKjK/6CdReDbv6+sHFgN1/slXstuOWxf+uRdp5px5VfODH1K+xH48rQq5j+wWjY4cFY9ToV4oP/2jTunXxMf5TX0nre/qvYs+dt4kVllsy9utxUiy3zBJR+9n/vPpGXHLO8cWHh5x14VXF77cttcRixT+rfQjJd9Zdsxg18Jpbiw8sqf1O3Ol9u8fSS35tluaKtPrdL65MgAABAgQI5AiItBznaVOqGGnLLt4+juh1dKy11jp11br/vnviqisuL2b8+oprY5vtdqjrPBf/RGDc+AlR+3usW7VsOY3kyz7S/oMPx0a7tvPPkK8Wf+3aLlB82mPt66LfXF+87XHLH25Y/P/vvf9BtF+wbZN/ifbH48ZF61atPvPPZnbuDA/mAQQIECBAgACBkguItOQFVTHSll+yY1x1zaDYtPNmddfas+uPo2OnheJXl11Z91kG1Ffg6WdHxTm/urr4vbcJEyYWf+9Z7UNJWrduVd/Brk6AAAECBAgQqLiASEteoEj7cvBD9t8rOi28cJx+1oDkzRhXL4FaoNX+Auupf6F1vea4LgECBAgQIEBgXhEQacmbFGkiLfmWM44AAQIECBAgQKBiAiIteWEiTaQl33LGESBAgAABAgQIVExApCUvTKSJtORbzjgCBAgQIECAAIGKCYi05IWJNJGWfMsZR4AAAQIECBAgUDEBkZa8MJEm0pJvOeMIECBAgAABAgQqJiDSkhcm0kRa8i1nHAECBAgQIECAQMUERFrywkSaSEu+5YwjQIAAAQIECBComIBIS16YSBNpybeccQQIECBAgAABAhUTEGnJCxNpIi35ljOOAAECBAgQIECgYgIiLXlhIk2kJd9yxhEgQIAAAQIECFRMQKQlL0ykibTkW844AgQIECBAgACBigmItOSFiTSRlnzLGUeAAAECBAgQIFAxAZGWvDCRJtKSbznjCBAgQIAAAQIEKiYg0pIXJtJEWvItZxwBAgQIECBAgEDFBERa8sJEmkhLvuWMI0CAAAECBAgQqJiASEtemEgTacm3nHEECBAgQIAAAQIVExBpyQsTaSIt+ZYzjgABAgQIECBAoGICIi15YSJNpCXfcsYRIECAAAECBAhUTECkJS9MpIm05FvOOAIECBAgQIAAgYoJiLTkhYk0kZZ8yxlHgAABAgQIECBQMQGRlrwwkSbSkm854wgQIECAAAECBComINKSFybSRFryLWccAQIECBAgQIBAxQREWvLCRJpIS77ljCNAgAABAgQIEKiYgEhLXphIE2nJt5xxBAgQIECAAAECFRMQackLE2kiLfmWM44AAQIECBAgQKBiAiIteWEiTaQl33LGESBAgAABAgQIVExApCUvTKSJtORbzjgCBAgQIECAAIGKCYi05IWJNJGWfMsZR4AAAQIECBAgUDEBkZa8MJEm0pJvOeMIECBAgAABAgQqJiDSkhcm0kRa8i1nHAECBAgQIECAQMUERFrywkSaSEu+5YwjQIAAAQIECBComIBIS16YSBNpybeccQQIECBAgAABAhUTEGnJCxNpIi35ljOOAAECBAgQIECgYgIiLXlhIk2kJd9yxhEgQIAAAQIECFRMQKQlL0ykibTkW844AgQIECBAgACBigmItOSFiTSRlnzLGUeAAAECBAgQIFAxAZGWvDCRJtKSbznjCBAgQIAAAQIEKiYg0pIXJtJEWvItZxwBAgQIECBAgEDFBERa8sJEmkhLvuWMI0CAAAECBAgQqJiASEtemEgTacm3nHEECBAgQIAAAQIVExBpyQsTaSIt+ZYzjgABAgQIECBAoGICIi15YSJNpCXfcsYRIECAAAECBAhUTECkJS9MpIm05FvOOAIECBAgQIAAgYoJiLTkhYk0kZZ8yxlHgAABAgQIECBQMQGRlrwwkSbSkm854wgQIECAAAECBComINKSFybSRFryLWccAQIECBAgQIBAxQREWvLCRJpIS77ljCNAgAABAgQIEKiYgEhLXphIE2nJt1xdxr33/gfxj6eei402+GZdru+iBAgQIECAAIFGFhBpydsXaSIt+Zabo3H3PvjXOP2XA6NjhwWjZYsWsfkPNoh9d98+/vX8y3HJFTfGgNN+NkfX98MECBAgQIAAAQJfFBBpyXeFSBNpybfcHI27+/7h8ewLL8ch++wc77z7fhx5/NnRp3u3aNWq1bRImzJlSjRr1uwLc2b1+9M76PSuM0dPzA8TIECAAAECBEosINKSlyPSRFq9b7n/ufzi+PjjsXHQIT2iZatWczTu05FWu9AJZ1wSP/z+t2PJxReLC399Xay8wjJx+9AH48db/yC6H7BrMWv4o0/GjbcNjVEvvxL77L59dNli4+L7Dz3yeFz/+yHx6utvxffW/0Ycuu/O8a/nXo4H/vJYHL7/LsVjzr/8utiy83dj9VVWaPLxbVq3nqPn44cJECBAgAABAlUQEGnJWxJpXw5+4D57xL+e+Wd869vrJ29m3hg36sXnY/hfhhVPpnvPn8VxfU+eoydWi7TaWx63/OH34t+vvB6/u/nO+J8LTowxb78X3Y85M07sc1Cs/fVV4oCeJ8dFZx4TkyZPjkN/dnpccWG/T85wTP+48MxjotNC7WOHvXrHJeccH19bdOG4dtAd8ZPtNot/PPV83Hnvw3HCUQcWjz+x/6WxwzY/iG+utVqTj19g/jZz9Hz8MAECBAgQIECgCgIiLXlLIu3LwceO/SiOP7pn8lbmnXH1iLRBf7g7vrXO6sWrYANO7xNrrrZiPPPcS3Hm+VdMi7HaK2y77bRlvPram3H19YPjBxt9u0C976ERccBeO8YmG34rfnnpNfHY35+OLlt8P7b4wQbRaaEO8cjf/hFD7hv2hUhbd+3Vm3z8vLMpz4QAAQIECBAgMH0BkZZ8d4i0ZPAGHFevtzued/FvY8Xll4odtulcRNqnPzik9grYzttvHq+9/lY8NHxk8ZipX8suvXh06tg+ar9b9tgTz8SQe4fFsL/+vXhF7qXRr8bv/3hfnHLcocXD+/QbEHv8ZKuoRVpTj1+kU8cG3KinTIAAAQIECDSagEhL3rhISwY3bo4EPv07aaNfeT169z03fnvJqTFq9KtNRtqiC3eMo/udHxeffXy0XaBNPPPcqFht5eXj/f9+GA//dWRsvdlGRXwd2OuUOOqwvaJd2/mj9wnnxcALToyJEyfFPt1PjH7HHBwrLrd0k49fY9UV5uj5+GECBAgQIECAQBUERFrylkRaMrhxcyTw+Q8OOWPAwFhz9ZVi9VWWj0uvuCnOO/Wo4vq1V9K6br95rLXGyvHHoQ/GZVcNig7t2xVvaay9SvbRRx/H2RddHa+89mYs2G6B4vfYah80Ugu2sy68Kh4ePjKWWmKxGPvxuDhsv66x7FKLN/n4OXoyfpgAAQIECBAgUBEBkZa8KJGWDG7cVyIwbvyEmDB+QrRrt8Bn5tcirEXz5tGyZYvPfL/2l2O3X7DtFz7Kf3qP/0qelKEECBAgQIAAgSQBkZYEPXWMSEsGN44AAQIECBAgQIBAxQREWvLCRFoyuHEECBAgQIAAAQIEKiYg0pIXJtKSwY0jQIAAAQIECBAgUDEBkZa8MJGWDG4cAQIECBAgQIAAgYoJiLTkhYm0ZHDjCBAgQIAAAQIECFRMQKQlL0ykJYMbR4AAAQIECBAgQKBiAiIteWEiLRncOAIECBAgQIAAAQIVExBpyQsTacngxhEgQIAAAQIECBComIBIS16YSEsGN44AAQIECBAgQIBAxQREWvLCRFoyuHEECBAgQIAAAQIEKiYg0pIXJtKSwY0jQIAAAQIECBAgUDEBkZa8MJGWDG4cAQIECBAgQIAAgYoJiLTkhYm0ZHDjCBAgQIAAAQIECFRMQKQlL0ykJYMbR4AAAQIECBAgQKBiAiIteWEiLRncOAIECBAgQIAAAQIVExBpyQsTacngxhEgQIAAAQIECBComIBIS16YSEsGN44AAQIECBAgQIBAxQREWvLCRFoyuHEECBAgQIAAAQIEKiYg0pIXJtKSwY0jQIAAAQIECBAgUDEBkZa8MJGWDG4cAQIECBAgQIAAgYoJiLTkhYm0ZHDjCBAgQIAAAQIECFRMQKQlL0ykJYMbR4AAAQIECBAgQKBiAiIteWEiLRncOAIECBAgQIAAAQIVExBpyQsTacngxhEgQIAAAQIECBComIBIS16YSEsGN44AAQIECBAgQIBAxQREWvLCRFoyuHEECBAgQIAAAQIEKiYg0pIXJtKSwY0jQIAAAQIECBAgUDEBkZa8MJGWDG4cAQIECBAgQIAAgYoJiLTkhYm0ZHDjCBAgQIAAAQIECFRMQKQlL0ykJYMbR4AAAQIECBAgQKBiAiIteWEiLRncOAIECBAgQIAAAQIVExBpyQsTacngxhEgQIAAAQIECBComIBIS16YSEsGN44AAQIECBAgQIBAxQREWvLCRFoyuHEECBAgQIAAAQIEKiYg0pIXJtKSwY0jQIAAAQIECBAgUDEBkZa8MJGWDG4cAQIECBAgQIAAgYoJiLTkhYm0ZHDjCBAgQIAAAQIECFRMQKQlL0ykJYMbR4AAAQIECBAgQKBiAiIteWEiLRncOAIECBAgQIAAAQIVExBpyQsTacngxhEgQIAAAQIECBComIBIS16YSEsGN44AAQIECBAgQIBAxQREWvLCRFoyuHEECBAgQIAAAQIEKiYg0pIXJtKSwY0jQIAAAQIECBAgUDEBkZa8MJGWDG4cAQIECBAgQIAAgYoJiLTkhYm0ZHDjCBAgQIAAAQIECFRMQKQlL0ykJYMbR4AAAQIECBAgQKBiAiIteWEiLRncOAIECBAgQIAAAQIVExBpyQsTacngxhEgQIAAAQIECBComIBIS16YSEsGN44AAQIECBAgQIBAxQREWvLCRFoyuHEEkgUmTZoc9zwwvJi60grLxIrLLZV8AuMIECBAgACBqguItOQNirRkcOMIzKTAkHuHxXMvjo7D9tul+ImLB94Qa6y6YnTe+NszeYVPHjZ58uQYNuKJeGDYo7HcMkvE7jttNUs/P/XBv/ntLfG1RReO7bbcpPjWn+55OF4a/WocvPdPZut6fogAAQIECBCojoBIS96VSEsGN47ATArcec+wePHl/8Qh++xc/MRFv7k+1llzldhkw2/N5BU++7DrbxkSk6dMbjLSpkyZEs2aNfvS61521aC468+PxHWXnx7RrFl0O/SE2PR768XB+8w40mbm+rP1pPwQAQIECBAgkCIg0lKY/29I1SLtiZGPRZcfbRKH9Tgqju3bL1nLOAJfLvDoiOFx3z13xWprfD26bLfDHHF9WaTd++CI+PXVN8eHH42NnofsWby6dmL/S6PHgbvFfM2axYQJE6NNm9bFq2/H9dyvOEdTkTb80SfjxtuGxqiXX4l9dt8+umyxcfHYPv0GxHfXWztuveO+aN68efQ6ZM8YNuLv8cehD8bh++9SfO/SK2+KzTfZoIi0ps7z7AsvFzMXmL9N3HH3Q7HTtj+MQ/ftOkcmfpgAAQIECBD4agREWrJ71SJt6802ir+PfKxQ6rL9DtG2bbtkMeMITF/gT7f/Id5//73iAY//88VYdLGvzTZXLdJuHnx3rLPmqsU1Hnvi6ei267ax7NJLxJHHnRXnnNwrWjRvHj1/fk6cf8bRcfMf7o7NNl0/nn3+5XjtjTHx4603jRtvvSt6H/bTJiPtldfejEN/dnpcceEnf9jR/Zj+ceGZx8TCnTpEt8NOiFatWsY5J/WKt95+N5rPN18Mue8vsVDH9jH4zvuLf7b9lpvGq6+/FVtu9r0mz1ME5PFnx8nHHhqrr7J87NP9xBh4Qb/i+r4IECBAgACBagmItOR9ibRkcOPmaYG5HWm1V69qsVX7GvSHu+NHnTeM198YE6NGvxp9uncrvn/+5dfFCssuWfzvSZMmxcgnn40XRv27eIXthZf+HTtvt3mTkVZ76+LV1w+OH2z0ye+43ffQiDhgrx2Lt1PWIq12/bW/vsq0fdXe7rj2GivHQ8NHxkIdFyx+P+4f/3wuOi3Uvsnz1D6k5ILLro3LzutbXOO4Uy6MvXfbrgg2XwQIECBAgEC1BERa8r6qFmne7ph8gxg3SwL1fLtjLcbWXXu1qL0C9vqbb8eRB+1enG3gtbdGq5YtYt21V4/BQx6I/37wYbRo0SIWXXih+PY3vx4brLdWk5E29L6/FMG1wzadpz3HZZdePDp1bF9EWu2VusUW6fSFSFtz9ZWiRYvm8fg//lVEWocO7Zo+zzprxJXX3RZn9+tZXOPnp10UP+3aJdZYdYVZMvVgAgQIECBA4KsXEGnJO6hapNV4ll+yY1x1zaDYtPNmyVrGEcgT+PzvpE2NtIU7dYx+/S8tXqFq07pVHH70mcXveq25+oqxw169Y+/dto2FOnaIAZdeE1dedFIstcRiTUbaa2+8FUf3Oz8uPvv4aLtAm3jmuVGx2srLFx8gUou0c0/uHYsustAXIu1763+j+F4t8GqRtvGG6zZ5nrZt5xdpebeLSQQIECBAoK4CIq2uvF+8uEhLBjeOwEwK1D6C/4WX/u/THadGWu3tiL+75c64btCfokP7drHVZt8rPrFxalz97PBuscjCHWPvw38Rd9xwUXz44dg4ut+A4vfUar9LVnulrP+JR0bHDgsWHwRSextj7TqdFuoQpxx3aCzYrm0RaeedelQs0qnjZyLtG2uuEt/99jqfibTaB4c0dZ6nnx1VvMo39ZW0vqf/KvbceRuvpM3k/j2MAAECBAiUSUCkJW9DpCWDG0dgNgU+/zH248aNLz4Kv3WrltOu+OnHzOzHkjxZzQAADk9JREFU3o8bPyEmjJ8Q7dot0OR1pn6zqet9+ntNnWc2n6ofI0CAAAECBEomINKSFyLSksGNI0CAAAECBAgQIFAxAZGWvDCRlgxuHAECBAgQIECAAIGKCYi05IWJtGRw4wgQIECAAAECBAhUTECkJS9MpCWDG0eAAAECBAgQIECgYgIiLXlhIi0Z3DgCBAgQIECAAAECFRMQackLE2nJ4MYRIECAAAECBAgQqJiASEtemEhLBjeOAAECBAgQIECAQMUERFrywkRaMrhxBAgQIECAAAECBComINKSFybSksGNI0CAAAECBAgQIFAxAZGWvDCRlgxuHAECBAgQIECAAIGKCYi05IWJtGRw4wgQIECAAAECBAhUTECkJS9MpCWDG0eAAAECBAgQIECgYgIiLXlhIi0Z3DgCBAgQIECAAAECFRMQackLE2nJ4MYRIECAAAECBAgQqJiASEtemEhLBjeOAAECBAgQIECAQMUERFrywkRaMrhxBAgQIECAAAECBComINKSFybSksGNI0CAAAECBAgQIFAxAZGWvDCRlgxuHAECBAgQIECAAIGKCYi05IWJtGRw4wgQIECAAAECBAhUTECkJS9MpCWDG0eAAAECBAgQIECgYgIiLXlhIi0Z3DgCBAgQIECAAAECFRMQackLE2nJ4MYRIECAAAECBAgQqJiASEtemEhLBjeOAAECBAgQIECAQMUERFrywkRaMrhxBAgQIECAAAECBComINKSFybSksGNI0CAAAECBAgQIFAxAZGWvDCRlgxuHAECBAgQIECAAIGKCYi05IWJtGRw4wgQIECAAAECBAhUTECkJS9MpCWDG0eAAAECBAgQIECgYgIiLXlhIi0Z3DgCBAgQIECAAAECFRMQackLE2nJ4MYRIECAAAECBAgQqJiASEtemEhLBjeOAAECBAgQIECAQMUERFrywkRaMrhxBAgQIECAAAECBComINKSFybSksGNI0CAAAECBAgQIFAxAZGWvDCRlgxuHAECBAgQIECAAIGKCYi05IWJtGRw4wgQIECAAAECBAhUTECkJS9MpCWDG0eAAAECBAgQIECgYgIiLXlhIi0Z3DgCBAgQIECAAAECFRMQackLE2nJ4MYRIECAAAECBAgQqJiASEtemEhLBjeOAAECBAgQIECAQMUERFrywkRaMrhxBAgQIECAAAECBComINKSFybSksGNI0CAAAECBAgQIFAxAZGWvDCRlgxuHAECBAgQIECAAIGKCYi05IWJtGRw4wgQIECAAAECBAhUTECkJS9MpCWDG0eAAAECBAgQIECgYgIiLXlhIi0Z3DgCBAgQIECAAAECFRMQackLE2nJ4MYRIECAAAECBAgQqJiASEtemEhLBjeOAAECBAgQIECAQMUERFrywkRaMrhxBAgQIECAAAECBComINKSFybSksGNI0CAAAECBAgQIFAxAZGWvDCRlgxuHAECBAgQIECAwFwRmDRpcjRvPt+XXqv2mHseGF48ZqUVlokVl1tqrsxutIuItOSNi7RkcOMIECBAgAABAgTmSODRvz8dAy69JsZ+PC5WW3m5OPaIfaNduwWavObkyZNj2Ign4oFhj8ZyyywRu++01RzNbtQfFmnJmxdpyeDGESBAgAABAgQIzLbA2+++H92PPjPO+EWPWG7pJeLK626L5ZZZMjpv/O0vveb1twyJyVMmz3SkffpVulrozTffl79iN9tPqCI/KNKSFyXSksGNI0CAAAECBAg0mMBlF59fPOODDztyjp/5kHuHxb+efzm6H7Brca0pU6ZEs2bNiv997aA7YvCdD8SEiRPjsP26RueNvzNtXlORNr3H1x57+9AHYpFOHaN9+3ax8EIdoseBu8W9D46IX199c3z40djoecieMwzDOX6yJbqASEteRhUjbdnF28dGG28aiy+xRLKWcQQIECBAgAABArMi8MzTT8XIx/5W/MgvTj5jjkOt9srZYot0im222PgLx3h4+MhY7xtrxJSYEgf1OjXOPaV3LLrwQsXjmoq0ph4/ceKkOO2838SA0/vEfQ+OiL+NfCqOOWKfGDX61TjyuLPinJN7RYvmzaPnz8+J8884OpZfpjH+e1SkzcpdPxceW8VIu2XQDXH/fXfNhWfvEgQIECBAgAABAvUUmNuRdumVN8UqKy4bm22yflx25aB4aPjjsdVm34s9frJ1vPPu+/HMcy/Fm2PeiZsH3xM9D9kjvrHmqtONtKYe36JFi7hk4A1FpN37wF/juRdHx6H7do0bbx1ahFqf7t2K651/+XWxwrJLxvZbbVpPvtJcW6Qlr6KKkZZMZBwBAgQIECBAgMAcCMzNtzvecvs98daYd+PAbjtF7ffG/jLi7/GPp56PnbffPHoc2z927NI5llx80Rj0h7tjr127xLprr95kpI15+70mH//NtVYrvj/6P68XMfiz7nvF4ostEr+75c54/c2348iDdi+uN/DaW6NVyxbx065d5kCmOj8q0pJ3JdKSwY0jQIAAAQIECBCYbYHX3xgTB/U+NS499+exxNcWidrvqL340iux8orLxIjH/1m8NfGDDz4qHtOnR7fpRtrd9w9v8vErr7BMHH70mXFa3+7RdoH5Y6EOCxa/8/bkMy9Ev/6XxmXn9Y02rVsVj6m9wrb+t9ac7edSpR8UacnbEmnJ4MYRIECAAAECBAjMkUDtAzwuHnhDdGjfLpZaYrHY5cdbxIrLL128Atax/YLRscOCMWr0K3H4/rvEisstHUf3GxCvvTEmWrVqGZ06to/+Jx4ZrVu3avLxa6y6QuzX46Ti4/prH/H/n1ffiEvOOb74EJHaq2nXDfpTMbf2Fsvax/lP/dCSOXpCFfhhkZa8JJGWDG4cAQIECBAgQIDAHAtMnDQpxo+fEAvM3+Yz1/rgw7HRru38M339zz/+6usHR7u2C8RO2/6wuMZFv7m+eNvjlj/csPj/x40bH9GsWbRu1XKmZ8wLDxRpyVsUacngxhEgQIAAAQIECJRW4OlnR8U5v7q6eEvjhAkTY9FFFooTjjqweOWtkb9EWvL2RVoyuHEECBAgQIAAAQKlF6gFWu0vsG7evLH/EuupixJpybesSEsGN44AAQIECBAgQIBAxQREWvLCRFoyuHEECBAgQIAAAQIEKiYg0pIXJtKSwY0jQIAAAQIECBAgUDEBkZa8MJGWDG4cAQIECBAgQIAAgYoJiLTkhYm0ZHDjCBAgQIAAAQIECFRMQKQlL0ykJYMbR4AAAQIECBAgQKBiAiIteWEiLRncOAIECBAgQIAAAQIVExBpyQsTacngxhEgQIAAAQIECBComIBIS16YSEsGN44AAQIECBAgQIBAxQREWvLCRFoyuHEECBAgQIAAAQIEKiYg0pIXJtKSwY0jQIAAAQIECBAgUDEBkZa8MJGWDG4cAQIECBAgQIAAgYoJiLTkhYm0ZHDjCBAgQIAAAQIECFRMQKQlL0ykJYMbR4AAAQIECBAgQKBiAiIteWEiLRncOAIECBAgQIAAAQIVExBpyQsTacngxhEgQIAAAQIECBComIBIS16YSEsGN44AAQIECBAgQIBAxQREWvLCRFoyuHEECBAgQIAAAQIEKiYg0pIXJtKSwY0jQIAAAQIECBAgUDEBkZa8MJGWDG4cAQIECBAgQIAAgYoJiLSKLcxxCRAgQIAAAQIECBAgUEaBZlOmTJlSxoM5EwECBAgQIECAAAECBBpRQKQ14tY9ZwIECBAgQIAAAQIESisg0kq7GgcjQIAAAQIECBAgQKARBURaI27dcyZAgAABAgQIECBAoLQCIq20q3EwAgQIECBAgAABAgQaUUCkNeLWPWcCBAgQIECAAAECBEorINJKuxoHI0CAAAECBAgQIECgEQVEWiNu3XMmQIAAAQIECBAgQKC0AiKttKtxMAIECBAgQIAAAQIEGlFApDXi1j1nAgQIECBAgAABAgRKKyDSSrsaByNAgAABAgQIECBAoBEFRFojbt1zJkCAAAECBAgQIECgtAIirbSrcTACBAgQIECAAAECBBpRQKQ14tY9ZwIECBAgQIAAAQIESisg0kq7GgcjQIAAAQIECBAgQKARBURaI27dcyZAgAABAgQIECBAoLQCIq20q3EwAgQIECBAgAABAgQaUUCkNeLWPWcCBAgQIECAAAECBEorINJKuxoHI0CAAAECBAgQIECgEQVEWiNu3XMmQIAAAQIECBAgQKC0AiKttKtxMAIECBAgQIAAAQIEGlFApDXi1j1nAgQIECBAgAABAgRKKyDSSrsaByNAgAABAgQIECBAoBEFRFojbt1zJkCAAAECBAgQIECgtAIirbSrcTACBAgQIECAAAECBBpRQKQ14tY9ZwIECBAgQIAAAQIESisg0kq7GgcjQIAAAQIECBAgQKARBURaI27dcyZAgAABAgQIECBAoLQCIq20q3EwAgQIECBAgAABAgQaUUCkNeLWPWcCBAgQIECAAAECBEorINJKuxoHI0CAAAECBAgQIECgEQVEWiNu3XMmQIAAAQIECBAgQKC0AiKttKtxMAIECBAgQIAAAQIEGlFApDXi1j1nAgQIECBAgAABAgRKKyDSSrsaByNAgAABAgQIECBAoBEFRFojbt1zJkCAAAECBAgQIECgtAIirbSrcTACBAgQIECAAAECBBpRQKQ14tY9ZwIECBAgQIAAAQIESisg0kq7GgcjQIAAAQIECBAgQKARBURaI27dcyZAgAABAgQIECBAoLQCIq20q3EwAgQIECBAgAABAgQaUUCkNeLWPWcCBAgQIECAAAECBEorINJKuxoHI0CAAAECBAgQIECgEQVEWiNu3XMmQIAAAQIECBAgQKC0AiKttKtxMAIECBAgQIAAAQIEGlHg/wFL0aUUnjyFVQAAAABJRU5ErkJggg==",
      "text/html": [
       "<div>\n",
       "        \n",
       "        \n",
       "            <div id=\"7b2ec838-45f0-45f9-87a9-d59d5cc12b3b\" class=\"plotly-graph-div\" style=\"height:500px; width:500px;\"></div>\n",
       "            <script type=\"text/javascript\">\n",
       "                require([\"plotly\"], function(Plotly) {\n",
       "                    window.PLOTLYENV=window.PLOTLYENV || {};\n",
       "                    \n",
       "                if (document.getElementById(\"7b2ec838-45f0-45f9-87a9-d59d5cc12b3b\")) {\n",
       "                    Plotly.newPlot(\n",
       "                        '7b2ec838-45f0-45f9-87a9-d59d5cc12b3b',\n",
       "                        [{\"line\": {\"color\": \"black\", \"shape\": \"spline\", \"smoothing\": 1.3, \"width\": 1.25}, \"mode\": \"lines\", \"showlegend\": false, \"type\": \"scatter\", \"x\": [null, null, 0, 0, null, 0, 0.03139195169925202, null, 0.03139195169925202, 0.03139195169925202, null, 0.03139195169925202, 0.08159484122097149, null, 0.08159484122097149, 0.08159484122097149, null, 0.08159484122097149, 0.09617900440577909, null, 0.09617900440577909, 0.09617900440577909, null, 0.09617900440577909, 0.0961878725442622, null, 0.0961878725442622, 0.0961878725442622, null, 0.0961878725442622, 0.10959571362620046, null, 0.0961878725442622, 0.0961878725442622, null, 0.0961878725442622, 0.12185690868725763, null, 0.09617900440577909, 0.09617900440577909, null, 0.09617900440577909, 0.1057597283305971, null, 0.08159484122097149, 0.08159484122097149, null, 0.08159484122097149, 0.11501597595635671, null, 0.03139195169925202, 0.03139195169925202, null, 0.03139195169925202, 0.12464469563427283, null, 0, 0, null, 0, 0.1930734042426634, null, 0, 0, null, 0, 0.7165441919466353, null], \"y\": [null, null, -2.6, -0.08125000000000004, null, -0.08125000000000004, -0.08125000000000004, null, -0.08125000000000004, 1.1375, null, 1.1375, 1.1375, null, 1.1375, 2.275, null, 2.275, 2.275, null, 2.275, 3.25, null, 3.25, 3.25, null, 3.25, 3.9000000000000004, null, 3.9000000000000004, 3.9000000000000004, null, 3.25, 2.6, null, 2.6, 2.6, null, 2.275, 1.3, null, 1.3, 1.3, null, 1.1375, 0.0, null, 0.0, 0.0, null, -0.08125000000000004, -1.3, null, -1.3, -1.3, null, -2.6, -2.6, null, -2.6, -2.6, null, -2.6, -3.9000000000000004, null, -3.9000000000000004, -3.9000000000000004, null]}, {\"hoverinfo\": \"text\", \"marker\": {\"color\": \"black\", \"size\": 3, \"symbol\": \"circle\"}, \"mode\": \"markers\", \"showlegend\": false, \"text\": [\"root\", \"edge.3\", \"edge.2\", \"edge.1\", \"edge.0\", \"Chimpanzee\", \"Human\", \"Gorilla\", \"Orangutan\", \"Rhesus\", \"HowlerMon\", \"Galago\"], \"type\": \"scatter\", \"x\": [0, 0.03139195169925202, 0.08159484122097149, 0.09617900440577909, 0.0961878725442622, 0.10959571362620046, 0.12185690868725763, 0.1057597283305971, 0.11501597595635671, 0.12464469563427283, 0.1930734042426634, 0.7165441919466353], \"y\": [-2.6, -0.08125000000000004, 1.1375, 2.275, 3.25, 3.9000000000000004, 2.6, 1.3, 0.0, -1.3, -2.6, -3.9000000000000004]}],\n",
       "                        {\"annotations\": [{\"font\": {\"family\": \"Inconsolata, monospace\", \"size\": 12}, \"showarrow\": false, \"text\": \"Chimpanzee\", \"textangle\": 0, \"x\": 0.12750931842486635, \"xanchor\": \"left\", \"xref\": \"x\", \"y\": 3.9000000000000004, \"yref\": \"y\"}, {\"font\": {\"family\": \"Inconsolata, monospace\", \"size\": 12}, \"showarrow\": false, \"text\": \"Human\", \"textangle\": 0, \"x\": 0.1397705134859235, \"xanchor\": \"left\", \"xref\": \"x\", \"y\": 2.6, \"yref\": \"y\"}, {\"font\": {\"family\": \"Inconsolata, monospace\", \"size\": 12}, \"showarrow\": false, \"text\": \"Gorilla\", \"textangle\": 0, \"x\": 0.12367333312926299, \"xanchor\": \"left\", \"xref\": \"x\", \"y\": 1.3, \"yref\": \"y\"}, {\"font\": {\"family\": \"Inconsolata, monospace\", \"size\": 12}, \"showarrow\": false, \"text\": \"Orangutan\", \"textangle\": 0, \"x\": 0.1329295807550226, \"xanchor\": \"left\", \"xref\": \"x\", \"y\": 0.0, \"yref\": \"y\"}, {\"font\": {\"family\": \"Inconsolata, monospace\", \"size\": 12}, \"showarrow\": false, \"text\": \"Rhesus\", \"textangle\": 0, \"x\": 0.1425583004329387, \"xanchor\": \"left\", \"xref\": \"x\", \"y\": -1.3, \"yref\": \"y\"}, {\"font\": {\"family\": \"Inconsolata, monospace\", \"size\": 12}, \"showarrow\": false, \"text\": \"HowlerMon\", \"textangle\": 0, \"x\": 0.21098700904132928, \"xanchor\": \"left\", \"xref\": \"x\", \"y\": -2.6, \"yref\": \"y\"}, {\"font\": {\"family\": \"Inconsolata, monospace\", \"size\": 12}, \"showarrow\": false, \"text\": \"Galago\", \"textangle\": 0, \"x\": 0.7344577967453012, \"xanchor\": \"left\", \"xref\": \"x\", \"y\": -3.9000000000000004, \"yref\": \"y\"}, {\"ax\": 0, \"ay\": 0, \"showarrow\": false, \"text\": \"0.07\", \"x\": 0.8240258207386306, \"xref\": \"x\", \"y\": 3.9000000000000004, \"yref\": \"y\"}], \"autosize\": true, \"font\": {\"family\": \"Balto\", \"size\": 14}, \"height\": 500, \"hovermode\": \"closest\", \"margin\": {\"b\": 50, \"l\": 50, \"pad\": 4, \"r\": 50, \"t\": 50}, \"shapes\": [{\"line\": {\"color\": \"black\", \"width\": 1.25}, \"type\": \"line\", \"x0\": 0.7165441919466353, \"x1\": 0.7881986111412989, \"y0\": 3.9000000000000004, \"y1\": 3.9000000000000004}], \"showlegend\": false, \"template\": {\"data\": {\"bar\": [{\"error_x\": {\"color\": \"#2a3f5f\"}, \"error_y\": {\"color\": \"#2a3f5f\"}, \"marker\": {\"line\": {\"color\": \"#E5ECF6\", \"width\": 0.5}}, \"type\": \"bar\"}], \"barpolar\": [{\"marker\": {\"line\": {\"color\": \"#E5ECF6\", \"width\": 0.5}}, \"type\": \"barpolar\"}], \"carpet\": [{\"aaxis\": {\"endlinecolor\": \"#2a3f5f\", \"gridcolor\": \"white\", \"linecolor\": \"white\", \"minorgridcolor\": \"white\", \"startlinecolor\": \"#2a3f5f\"}, \"baxis\": {\"endlinecolor\": \"#2a3f5f\", \"gridcolor\": \"white\", \"linecolor\": \"white\", \"minorgridcolor\": \"white\", \"startlinecolor\": \"#2a3f5f\"}, \"type\": \"carpet\"}], \"choropleth\": [{\"colorbar\": {\"outlinewidth\": 0, \"ticks\": \"\"}, \"type\": \"choropleth\"}], \"contour\": [{\"colorbar\": {\"outlinewidth\": 0, \"ticks\": \"\"}, \"colorscale\": [[0.0, \"#0d0887\"], [0.1111111111111111, \"#46039f\"], [0.2222222222222222, \"#7201a8\"], [0.3333333333333333, \"#9c179e\"], [0.4444444444444444, \"#bd3786\"], [0.5555555555555556, \"#d8576b\"], [0.6666666666666666, \"#ed7953\"], [0.7777777777777778, \"#fb9f3a\"], [0.8888888888888888, \"#fdca26\"], [1.0, \"#f0f921\"]], \"type\": \"contour\"}], \"contourcarpet\": [{\"colorbar\": {\"outlinewidth\": 0, \"ticks\": \"\"}, \"type\": \"contourcarpet\"}], \"heatmap\": [{\"colorbar\": {\"outlinewidth\": 0, \"ticks\": \"\"}, \"colorscale\": [[0.0, \"#0d0887\"], [0.1111111111111111, \"#46039f\"], [0.2222222222222222, \"#7201a8\"], [0.3333333333333333, \"#9c179e\"], [0.4444444444444444, \"#bd3786\"], [0.5555555555555556, \"#d8576b\"], [0.6666666666666666, \"#ed7953\"], [0.7777777777777778, \"#fb9f3a\"], [0.8888888888888888, \"#fdca26\"], [1.0, \"#f0f921\"]], \"type\": \"heatmap\"}], \"heatmapgl\": [{\"colorbar\": {\"outlinewidth\": 0, \"ticks\": \"\"}, \"colorscale\": [[0.0, \"#0d0887\"], [0.1111111111111111, \"#46039f\"], [0.2222222222222222, \"#7201a8\"], [0.3333333333333333, \"#9c179e\"], [0.4444444444444444, \"#bd3786\"], [0.5555555555555556, \"#d8576b\"], [0.6666666666666666, \"#ed7953\"], [0.7777777777777778, \"#fb9f3a\"], [0.8888888888888888, \"#fdca26\"], [1.0, \"#f0f921\"]], \"type\": \"heatmapgl\"}], \"histogram\": [{\"marker\": {\"colorbar\": {\"outlinewidth\": 0, \"ticks\": \"\"}}, \"type\": \"histogram\"}], \"histogram2d\": [{\"colorbar\": {\"outlinewidth\": 0, \"ticks\": \"\"}, \"colorscale\": [[0.0, \"#0d0887\"], [0.1111111111111111, \"#46039f\"], [0.2222222222222222, \"#7201a8\"], [0.3333333333333333, \"#9c179e\"], [0.4444444444444444, \"#bd3786\"], [0.5555555555555556, \"#d8576b\"], [0.6666666666666666, \"#ed7953\"], [0.7777777777777778, \"#fb9f3a\"], [0.8888888888888888, \"#fdca26\"], [1.0, \"#f0f921\"]], \"type\": \"histogram2d\"}], \"histogram2dcontour\": [{\"colorbar\": {\"outlinewidth\": 0, \"ticks\": \"\"}, \"colorscale\": [[0.0, \"#0d0887\"], [0.1111111111111111, \"#46039f\"], [0.2222222222222222, \"#7201a8\"], [0.3333333333333333, \"#9c179e\"], [0.4444444444444444, \"#bd3786\"], [0.5555555555555556, \"#d8576b\"], [0.6666666666666666, \"#ed7953\"], [0.7777777777777778, \"#fb9f3a\"], [0.8888888888888888, \"#fdca26\"], [1.0, \"#f0f921\"]], \"type\": \"histogram2dcontour\"}], \"mesh3d\": [{\"colorbar\": {\"outlinewidth\": 0, \"ticks\": \"\"}, \"type\": \"mesh3d\"}], \"parcoords\": [{\"line\": {\"colorbar\": {\"outlinewidth\": 0, \"ticks\": \"\"}}, \"type\": \"parcoords\"}], \"scatter\": [{\"marker\": {\"colorbar\": {\"outlinewidth\": 0, \"ticks\": \"\"}}, \"type\": \"scatter\"}], \"scatter3d\": [{\"line\": {\"colorbar\": {\"outlinewidth\": 0, \"ticks\": \"\"}}, \"marker\": {\"colorbar\": {\"outlinewidth\": 0, \"ticks\": \"\"}}, \"type\": \"scatter3d\"}], \"scattercarpet\": [{\"marker\": {\"colorbar\": {\"outlinewidth\": 0, \"ticks\": \"\"}}, \"type\": \"scattercarpet\"}], \"scattergeo\": [{\"marker\": {\"colorbar\": {\"outlinewidth\": 0, \"ticks\": \"\"}}, \"type\": \"scattergeo\"}], \"scattergl\": [{\"marker\": {\"colorbar\": {\"outlinewidth\": 0, \"ticks\": \"\"}}, \"type\": \"scattergl\"}], \"scattermapbox\": [{\"marker\": {\"colorbar\": {\"outlinewidth\": 0, \"ticks\": \"\"}}, \"type\": \"scattermapbox\"}], \"scatterpolar\": [{\"marker\": {\"colorbar\": {\"outlinewidth\": 0, \"ticks\": \"\"}}, \"type\": \"scatterpolar\"}], \"scatterpolargl\": [{\"marker\": {\"colorbar\": {\"outlinewidth\": 0, \"ticks\": \"\"}}, \"type\": \"scatterpolargl\"}], \"scatterternary\": [{\"marker\": {\"colorbar\": {\"outlinewidth\": 0, \"ticks\": \"\"}}, \"type\": \"scatterternary\"}], \"surface\": [{\"colorbar\": {\"outlinewidth\": 0, \"ticks\": \"\"}, \"colorscale\": [[0.0, \"#0d0887\"], [0.1111111111111111, \"#46039f\"], [0.2222222222222222, \"#7201a8\"], [0.3333333333333333, \"#9c179e\"], [0.4444444444444444, \"#bd3786\"], [0.5555555555555556, \"#d8576b\"], [0.6666666666666666, \"#ed7953\"], [0.7777777777777778, \"#fb9f3a\"], [0.8888888888888888, \"#fdca26\"], [1.0, \"#f0f921\"]], \"type\": \"surface\"}], \"table\": [{\"cells\": {\"fill\": {\"color\": \"#EBF0F8\"}, \"line\": {\"color\": \"white\"}}, \"header\": {\"fill\": {\"color\": \"#C8D4E3\"}, \"line\": {\"color\": \"white\"}}, \"type\": \"table\"}]}, \"layout\": {\"annotationdefaults\": {\"arrowcolor\": \"#2a3f5f\", \"arrowhead\": 0, \"arrowwidth\": 1}, \"colorscale\": {\"diverging\": [[0, \"#8e0152\"], [0.1, \"#c51b7d\"], [0.2, \"#de77ae\"], [0.3, \"#f1b6da\"], [0.4, \"#fde0ef\"], [0.5, \"#f7f7f7\"], [0.6, \"#e6f5d0\"], [0.7, \"#b8e186\"], [0.8, \"#7fbc41\"], [0.9, \"#4d9221\"], [1, \"#276419\"]], \"sequential\": [[0.0, \"#0d0887\"], [0.1111111111111111, \"#46039f\"], [0.2222222222222222, \"#7201a8\"], [0.3333333333333333, \"#9c179e\"], [0.4444444444444444, \"#bd3786\"], [0.5555555555555556, \"#d8576b\"], [0.6666666666666666, \"#ed7953\"], [0.7777777777777778, \"#fb9f3a\"], [0.8888888888888888, \"#fdca26\"], [1.0, \"#f0f921\"]], \"sequentialminus\": [[0.0, \"#0d0887\"], [0.1111111111111111, \"#46039f\"], [0.2222222222222222, \"#7201a8\"], [0.3333333333333333, \"#9c179e\"], [0.4444444444444444, \"#bd3786\"], [0.5555555555555556, \"#d8576b\"], [0.6666666666666666, \"#ed7953\"], [0.7777777777777778, \"#fb9f3a\"], [0.8888888888888888, \"#fdca26\"], [1.0, \"#f0f921\"]]}, \"colorway\": [\"#636efa\", \"#EF553B\", \"#00cc96\", \"#ab63fa\", \"#FFA15A\", \"#19d3f3\", \"#FF6692\", \"#B6E880\", \"#FF97FF\", \"#FECB52\"], \"font\": {\"color\": \"#2a3f5f\"}, \"geo\": {\"bgcolor\": \"white\", \"lakecolor\": \"white\", \"landcolor\": \"#E5ECF6\", \"showlakes\": true, \"showland\": true, \"subunitcolor\": \"white\"}, \"hoverlabel\": {\"align\": \"left\"}, \"hovermode\": \"closest\", \"mapbox\": {\"style\": \"light\"}, \"paper_bgcolor\": \"white\", \"plot_bgcolor\": \"#E5ECF6\", \"polar\": {\"angularaxis\": {\"gridcolor\": \"white\", \"linecolor\": \"white\", \"ticks\": \"\"}, \"bgcolor\": \"#E5ECF6\", \"radialaxis\": {\"gridcolor\": \"white\", \"linecolor\": \"white\", \"ticks\": \"\"}}, \"scene\": {\"xaxis\": {\"backgroundcolor\": \"#E5ECF6\", \"gridcolor\": \"white\", \"gridwidth\": 2, \"linecolor\": \"white\", \"showbackground\": true, \"ticks\": \"\", \"zerolinecolor\": \"white\"}, \"yaxis\": {\"backgroundcolor\": \"#E5ECF6\", \"gridcolor\": \"white\", \"gridwidth\": 2, \"linecolor\": \"white\", \"showbackground\": true, \"ticks\": \"\", \"zerolinecolor\": \"white\"}, \"zaxis\": {\"backgroundcolor\": \"#E5ECF6\", \"gridcolor\": \"white\", \"gridwidth\": 2, \"linecolor\": \"white\", \"showbackground\": true, \"ticks\": \"\", \"zerolinecolor\": \"white\"}}, \"shapedefaults\": {\"line\": {\"color\": \"#2a3f5f\"}}, \"ternary\": {\"aaxis\": {\"gridcolor\": \"white\", \"linecolor\": \"white\", \"ticks\": \"\"}, \"baxis\": {\"gridcolor\": \"white\", \"linecolor\": \"white\", \"ticks\": \"\"}, \"bgcolor\": \"#E5ECF6\", \"caxis\": {\"gridcolor\": \"white\", \"linecolor\": \"white\", \"ticks\": \"\"}}, \"title\": {\"x\": 0.05}, \"xaxis\": {\"automargin\": true, \"gridcolor\": \"white\", \"linecolor\": \"white\", \"ticks\": \"\", \"zerolinecolor\": \"white\", \"zerolinewidth\": 2}, \"yaxis\": {\"automargin\": true, \"gridcolor\": \"white\", \"linecolor\": \"white\", \"ticks\": \"\", \"zerolinecolor\": \"white\", \"zerolinewidth\": 2}}}, \"width\": 500, \"xaxis\": {\"visible\": false}, \"yaxis\": {\"visible\": false}},\n",
       "                        {\"responsive\": true}\n",
       "                    ).then(function(){\n",
       "                            \n",
       "var gd = document.getElementById('7b2ec838-45f0-45f9-87a9-d59d5cc12b3b');\n",
       "var x = new MutationObserver(function (mutations, observer) {{\n",
       "        var display = window.getComputedStyle(gd).display;\n",
       "        if (!display || display === 'none') {{\n",
       "            console.log([gd, 'removed!']);\n",
       "            Plotly.purge(gd);\n",
       "            observer.disconnect();\n",
       "        }}\n",
       "}});\n",
       "\n",
       "// Listen for the removal of the full notebook cells\n",
       "var notebookContainer = gd.closest('#notebook-container');\n",
       "if (notebookContainer) {{\n",
       "    x.observe(notebookContainer, {childList: true});\n",
       "}}\n",
       "\n",
       "// Listen for the clearing of the current output cell\n",
       "var outputEl = gd.closest('.output');\n",
       "if (outputEl) {{\n",
       "    x.observe(outputEl, {childList: true});\n",
       "}}\n",
       "\n",
       "                        })\n",
       "                };\n",
       "                });\n",
       "            </script>\n",
       "        </div>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "tree = result.tree\n",
    "fig = tree.get_figure()\n",
    "fig.scale_bar = \"top right\"\n",
    "fig.show(width=500, height=500)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Getting parameter estimates\n",
    "\n",
    "For a discrete-time model, aside from the root motif probabilities, everything is edge specific. But note that the `tabular_result` has different keys from the continuous-time case, as demonstrated below."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "2x tabular_result('edge motif motif2 params': Table, 'motif params': Table)"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "tabulator = evo.tabulate_stats()\n",
    "stats = tabulator(result)\n",
    "stats"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<style>\n",
       "tr:last-child {border-bottom: 1px solid #000;} tr > th {text-align: center !important;} tr > td {text-align: left !important;}\n",
       "</style>\n",
       "<caption style=\"color: rgb(250, 250, 250); background: rgba(30, 140, 200, 1); align=top;\"><span style=\"font-weight: bold;\">edge motif motif2 params</span><span></span></caption>\n",
       "<thead style=\"background: rgba(161, 195, 209, 0.75); font-weight: bold; text-align: center;\">\n",
       "<th>edge</th>\n",
       "<th>motif</th>\n",
       "<th>motif2</th>\n",
       "<th>psubs</th>\n",
       "</thead>\n",
       "<tbody>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">Galago</td>\n",
       "<td>T</td>\n",
       "<td>T</td>\n",
       "<td style=\"font-family: monospace !important;\">0.8751</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">Galago</td>\n",
       "<td>T</td>\n",
       "<td>C</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0649</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">Galago</td>\n",
       "<td>T</td>\n",
       "<td>A</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0409</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">Galago</td>\n",
       "<td>T</td>\n",
       "<td>G</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0192</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">Galago</td>\n",
       "<td>C</td>\n",
       "<td>T</td>\n",
       "<td style=\"font-family: monospace !important;\">0.1126</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">...</td>\n",
       "<td>...</td>\n",
       "<td>...</td>\n",
       "<td>...</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">edge.3</td>\n",
       "<td>A</td>\n",
       "<td>G</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0055</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">edge.3</td>\n",
       "<td>G</td>\n",
       "<td>T</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0000</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">edge.3</td>\n",
       "<td>G</td>\n",
       "<td>C</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0011</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">edge.3</td>\n",
       "<td>G</td>\n",
       "<td>A</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0039</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">edge.3</td>\n",
       "<td>G</td>\n",
       "<td>G</td>\n",
       "<td style=\"font-family: monospace !important;\">0.9950</td>\n",
       "</tr>\n",
       "</tbody>\n",
       "</table>\n",
       "<p>\n",
       "176 rows x 4 columns</p>"
      ],
      "text/plain": [
       "edge motif motif2 params\n",
       "===================================\n",
       "  edge    motif    motif2     psubs\n",
       "-----------------------------------\n",
       "Galago        T         T    0.8751\n",
       "Galago        T         C    0.0649\n",
       "Galago        T         A    0.0409\n",
       "Galago        T         G    0.0192\n",
       "Galago        C         T    0.1126\n",
       "   ...      ...       ...       ...\n",
       "edge.3        A         G    0.0055\n",
       "edge.3        G         T    0.0000\n",
       "edge.3        G         C    0.0011\n",
       "edge.3        G         A    0.0039\n",
       "edge.3        G         G    0.9950\n",
       "-----------------------------------\n",
       "\n",
       "176 rows x 4 columns"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "stats['edge motif motif2 params']"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python [conda env:c3dev] *",
   "language": "python",
   "name": "conda-env-c3dev-py"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.3"
  },
  "widgets": {
   "application/vnd.jupyter.widget-state+json": {
    "state": {},
    "version_major": 2,
    "version_minor": 0
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}
