{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Specifying a non-stationary model with time-heterogeneous parameters\n",
    "\n",
    "A model is considered time-heterogeneous if rate matrix terms differ between edges on a tree.\n",
    "\n",
    "In this example, we load the and display our sample primate tree, highlighting the edge sets that will be assigned different rate matrices."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "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,
          1,
          null,
          0,
          0,
          null,
          0,
          1,
          null,
          0,
          0,
          null,
          0,
          0.2,
          null,
          0.2,
          0.2,
          null,
          0.2,
          1,
          null,
          0.2,
          0.2,
          null,
          0.2,
          0.4,
          null
         ],
         "y": [
          null,
          null,
          2.6,
          3.9000000000000004,
          null,
          3.9000000000000004,
          3.9000000000000004,
          null,
          2.6,
          2.6,
          null,
          2.6,
          2.6,
          null,
          2.6,
          0.08125000000000004,
          null,
          0.08125000000000004,
          0.08125000000000004,
          null,
          0.08125000000000004,
          1.3,
          null,
          1.3,
          1.3,
          null,
          0.08125000000000004,
          -1.1375,
          null,
          -1.1375,
          -1.1375,
          null
         ]
        },
        {
         "line": {
          "color": "red",
          "width": 1.25
         },
         "mode": "lines",
         "type": "scatter",
         "x": [
          0.4,
          0.4,
          null,
          0.4,
          1,
          null,
          0.4,
          0.4,
          null,
          0.4,
          0.6000000000000001,
          null,
          0.6000000000000001,
          0.6000000000000001,
          null,
          0.6000000000000001,
          1,
          null,
          0.6000000000000001,
          0.6000000000000001,
          null,
          0.6000000000000001,
          0.8,
          null,
          0.8,
          0.8,
          null,
          0.8,
          1,
          null,
          0.8,
          0.8,
          null,
          0.8,
          1,
          null
         ],
         "y": [
          -1.1375,
          0,
          null,
          0,
          0,
          null,
          -1.1375,
          -2.275,
          null,
          -2.275,
          -2.275,
          null,
          -2.275,
          -1.3,
          null,
          -1.3,
          -1.3,
          null,
          -2.275,
          -3.25,
          null,
          -3.25,
          -3.25,
          null,
          -3.25,
          -2.6,
          null,
          -2.6,
          -2.6,
          null,
          -3.25,
          -3.9000000000000004,
          null,
          -3.9000000000000004,
          -3.9000000000000004,
          null
         ]
        },
        {
         "hoverinfo": "text",
         "marker": {
          "color": "black",
          "size": 3,
          "symbol": "circle"
         },
         "mode": "markers",
         "showlegend": false,
         "text": [
          "root",
          "Galago",
          "HowlerMon",
          "edge.3",
          "Rhesus",
          "edge.2",
          "Orangutan",
          "edge.1",
          "Gorilla",
          "edge.0",
          "Chimpanzee",
          "Human"
         ],
         "type": "scatter",
         "x": [
          0,
          1,
          1,
          0.2,
          1,
          0.4,
          1,
          0.6000000000000001,
          1,
          0.8,
          1,
          1
         ],
         "y": [
          2.6,
          3.9000000000000004,
          2.6,
          0.08125000000000004,
          1.3,
          -1.1375,
          0,
          -2.275,
          -1.3,
          -3.25,
          -2.6,
          -3.9000000000000004
         ]
        }
       ],
       "layout": {
        "annotations": [
         {
          "font": {
           "family": "Inconsolata, monospace",
           "size": 12
          },
          "showarrow": false,
          "text": "Galago",
          "textangle": 0,
          "x": 1.025,
          "xanchor": "left",
          "xref": "x",
          "y": 3.9000000000000004,
          "yref": "y"
         },
         {
          "font": {
           "family": "Inconsolata, monospace",
           "size": 12
          },
          "showarrow": false,
          "text": "HowlerMon",
          "textangle": 0,
          "x": 1.025,
          "xanchor": "left",
          "xref": "x",
          "y": 2.6,
          "yref": "y"
         },
         {
          "font": {
           "family": "Inconsolata, monospace",
           "size": 12
          },
          "showarrow": false,
          "text": "Rhesus",
          "textangle": 0,
          "x": 1.025,
          "xanchor": "left",
          "xref": "x",
          "y": 1.3,
          "yref": "y"
         },
         {
          "font": {
           "family": "Inconsolata, monospace",
           "size": 12
          },
          "showarrow": false,
          "text": "Orangutan",
          "textangle": 0,
          "x": 1.025,
          "xanchor": "left",
          "xref": "x",
          "y": 0,
          "yref": "y"
         },
         {
          "font": {
           "family": "Inconsolata, monospace",
           "size": 12
          },
          "showarrow": false,
          "text": "Gorilla",
          "textangle": 0,
          "x": 1.025,
          "xanchor": "left",
          "xref": "x",
          "y": -1.3,
          "yref": "y"
         },
         {
          "font": {
           "family": "Inconsolata, monospace",
           "size": 12
          },
          "showarrow": false,
          "text": "Chimpanzee",
          "textangle": 0,
          "x": 1.025,
          "xanchor": "left",
          "xref": "x",
          "y": -2.6,
          "yref": "y"
         },
         {
          "font": {
           "family": "Inconsolata, monospace",
           "size": 12
          },
          "showarrow": false,
          "text": "Human",
          "textangle": 0,
          "x": 1.025,
          "xanchor": "left",
          "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
        },
        "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.07531948881789136,
          1.234584664536741
         ],
         "type": "linear",
         "visible": false
        },
        "yaxis": {
         "autorange": true,
         "range": [
          -4.40677966101695,
          4.40677966101695
         ],
         "type": "linear",
         "visible": false
        }
       }
      },
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABB0AAAH0CAYAAABxdlVzAAAgAElEQVR4Xu3debxe87nw/ytzZA6h5qnGEqpaqhQpDk1UUVO1QhU1BDHEVCpFDTXUVKX1qDqnVDWGoooaKiKVqqm0lKrh1JiYGtLMz+tejlTYiWzZ97XWznrnn985ce91fdf7Wr8/ns9Z+06HmTNnzgx/CBAgQIAAAQIECBAgQIAAAQJtLNBBdGhjUZcjQIAAAQIECBAgQIAAAQIECgHRwYNAgAABAgQIECBAgAABAgQINEVAdGgKq4sSIECAAAECBAgQIECAAAECooNngAABAgQIECBAgAABAgQIEGiKgOjQFFYXJUCAAAECBAgQIECAAAECBEQHzwABAgQIECBAgAABAgQIECDQFAHRoSmsLkqAAAECBAgQIECAAAECBAiIDp4BAgQIECBAgAABAgQIECBAoCkCokNTWF2UAAECBAgQIECAAAECBAgQEB08AwQIECBAgAABAgQIECBAgEBTBESHprC6KAECBAgQIECAAAECBAgQICA6eAYIECBAgAABAgQIECBAgACBpgiIDk1hdVECBAgQIECAAAECBAgQIEBAdPAMECBAgAABAgQIECBAgAABAk0REB2awuqiBAgQIECAAAECBAgQIECAgOjgGSBAgAABAgQIECBAgAABAgSaIiA6NIXVRQkQIECAAAECBAgQIECAAAHRwTNAgAABAgQIECBAgAABAgQINEVAdGgKq4sSIECAAAECBAgQIECAAAECooNngAABAgQIECBAgAABAgQIEGiKgOjQFFYXJUCAAAECBAgQIECAAAECBEQHzwABAgQIECBAgAABAgQIECDQFAHRoSmsLkqAAAECBAgQIECAAAECBAiIDp4BAgQIECBAgAABAgQIECBAoCkCokNTWF2UAAECBAgQIECAAAECBAgQEB08AwQIECBAgAABAgQIECBAgEBTBESHprC6KAECBAgQIECAAAECBAgQICA6eAYIECBAgAABAgQIECBAgACBpgiIDk1hdVECBAgQIECAAAECBAgQIEBAdPAMECBAgAABAgQIECBAgAABAk0REB2awuqiBAgQIECAAAECBAgQIECAgOjgGSBAgAABAgQIECBAgAABAgSaIiA6NIXVRQkQIECAAAECBAgQIECAAAHRwTNAgAABAgQIECBAgAABAgQINEVAdGgKq4sSIECAAAECBAgQIECAAAECooNngAABAgQIECBAgAABAgQIEGiKgOjQFFYXJUCAAAECBAgQIECAAAECBEQHzwABAgQIECBAgAABAgQIECDQFAHRoSmsLkqAAAECBAgQIECAAAECBAiIDp4BAgQIECBAgAABAgQIECBAoCkCokNTWF2UAAECBAgQIECAAAECBAgQEB08AwQIECBAgAABAgQIECBAgEBTBESHprC6KAECBAgQIECAAAECBAgQICA6eAYIECBAgAABAgQIECBAgACBpgiIDk1hdVECBAgQIECAAAECBAgQIECg3UaH5ydMsj0CBAgQIECAAAECBAgQIFALgSUXWahd3qfo0C7X5tAECBAgQIAAAQIECBAgUCcB0SF52950SAY3jgABAgQIECBAgAABAgRKExAdkulFh2Rw4wgQIECAAAECBAgQIECgNAHRIZledEgGN44AAQIECBAgQIAAAQIEShMQHZLpRYdkcOMIECBAgAABAgQIECBAoDQB0SGZXnRIBjeOAAECBAgQIECAAAECBEoTEB2S6UWHZHDjCBAgQIAAAQIECBAgQKA0AdEhmV50SAY3jgABAgQIECBAgAABAgRKExAdkulFh2Rw4wgQIECAAAECBAgQIECgNAHRIZledEgGN44AAQIECBAgQIAAAQIEShMQHZLpRYdkcOMIECBAgAABAgQIECBAoDQB0SGZXnRIBjeOAAECBAgQIECAAAECBEoTEB2S6UWHZHDjCBAgQIAAAQIECBAgQKA0AdEhmV50SAY3jgABAgQIECBAgAABAgRKExAdkulFh2Rw4wgQIECAAAECBAgQIECgNAHRIZledEgGN44AAQIECBAgQIAAAQIEShMQHZLpRYdkcOMIECBAgAABAgQIECBAoDQB0SGZXnRIBjeOAAECBAgQIECAAAECBEoTEB2S6UWHZHDjCBAgQIAAAQIECBAgQKA0AdEhmV50SAY3jgABAgQIECBAgAABAgRKExAdkulFh2Rw4wgQIECAAAECBAgQIECgNAHRIZledEgGN44AAQIECBAgQIAAAQIlC0yfPiM6deo411M0PnP76HHFZz6+wjKx4nJLlXzqthkvOrSN4zxfRXSYZyofJECAAAECBAgQIECAQLsWuP/hx+LsC38ek/49OVZdabk46qBvRK9ePVq8pxkzZsTY+/4co8feH8sts0R8dfut2vW9v3t40SF5jaJDMrhxBAgQIECAAAECBAgQKEHg1dffjGFHnBqnfOfAWG7pJeLSK34dyy2zZAza6NNzPc2V19wSM2bOmOfo8N63KBrhomPHub9RkU0hOiSLiw7J4MYRIECAAAECBAgQIECgBYGxY+6KsWNGx8C114ktthzc5ka33DE2/vb3Z2PYXjsX1545c2Z06NCh+J8vH3VT3HDz6Jg6bVrsv+eOMWijz8ya31J0mNPnG5+98dbRMWDhftGnT69YpH/fOHDvXeKOu++Ln1x2dbz19qQYvu/XPjR0tPnNv+eCokMzdVu49oIaHa4Z9cu4687fJWsaR4AAAQIECBAgQIAAgY8mcP21o2LSpEnFDz/x7Pjo0aPlX3v4aFeP4s2GxQYsHIO32OgDl7hn3EOx7tqrx8yYGfscclKceeKhsegi/YvPtRQdWvr8tGnT43tnXRxnnzwi7rz7vvjTQ3+NIw/aI55+7oU4+OjvxxknHBKdO3WK4d8+I8455YhYfpklPuqtzNfPiQ7zxdf6H15Qo8Oyi/eJDTfaJBZfopwHufWb8BMECBAgQIAAAQIECNRZoNnR4cJLfxUrr7hsbLbxenHRpaNizLgHY6vNPhe7fuWL8drrb8bjTz4Tr0x4La6+4fYYvu+usfYaq8wxOrT0+c6dO8ePLvllER3uGP3HePIfz8V+39gxrrru1iI8jBg2tLjeOT++IlZYdsnYZqtNSlm36JDMvqBGh+WX7Bc/+/mo2GTQZsmixhEgQIAAAQIECBAgQKD1As3+9Yprbrw9xk94PfYeun00vnfhD/c9HI/89e+xwzabx4FHnRbbDRkUSy6+aIy6/rbYbechsc7A1VqMDhNefaPFz39yzVWLv3/uny8VcePwYbvF4osNiF9cc3O89MqrcfA+Xy2ud8nl10XXLp3j6zsOaT1SG/yE6NAGiK25hOjQGi2fJUCAAAECBAgQIECAQPsUeOnlCbHPoSfFhWd+O5b42IBofMfDP555PlZacZm478G/FL8KMXHi28VnRhw4dI7R4ba7xrX4+ZVWWCYOOOLU+N6xw6Jnj4Wif9/exXdGPPr4UzHytAvjorOOje7duhafabwBsd6n1igFUnRIZhcdksGNI0CAAAECBAgQIECAQEkCjS90vOCSX0bfPr1iqSUWi52+vEWsuPzSxRsK/fr0jn59e8fTzz0fB3xzp1hxuaXjiJFnx4svT4iuXbvEwv36xGnHHxzdunVt8fOrr7JC7Hngd4t/XrPxT3L+84WX40dnHFN8qWTjbYcrRv22mNv4lY7GP7/57pdYZlOIDsniokMyuHEECBAgQIAAAQIECBAoUWDa9OkxZcrU6LFQ99lOMfGtSdGr50LzfLL3f/6yK2+IXj17xPZbf6G4xvkXX1n8msWWX9ig+N8nT54S0aFDdOvaZZ5nNOODokMzVOdyTdEhGdw4AgQIECBAgAABAgQILIACjz3xdJzxw8uKX6GYOnVaLDqgfxx32N7FmxFV+iM6JG9DdEgGN44AAQIECBAgQIAAAQILsEAjOHTs2DE6depYybsUHZLXIjokgxtHgAABAgQIECBAgAABAqUJiA7J9KJDMrhxBAgQIECAAAECBAgQIFCagOiQTC86JIMbR4AAAQIECBAgQIAAAQKlCYgOyfSiQzK4cQQIECBAgAABAgQIECBQmoDokEwvOiSDG0eAAAECBAgQIECAAAECpQmIDsn0okMyuHEECBAgQIAAAQIECBAgUJqA6JBMLzokgxtHgAABAgQIECBAgAABAqUJiA7J9KJDMrhxBAgQIECAAAECBAgQIFCagOiQTC86JIMbR4AAAQIECBAgQIAAAQKlCYgOyfSiQzK4cQQIECBAgAABAgQIECBQmoDokEwvOiSDG0eAAAECBAgQIECAAAECpQmIDsn0okMyuHEECBAgQIAAAQIECBAgUJqA6JBMLzokgxtHgAABAgQIECBAgAABAqUJiA7J9KJDMrhxBAgQIECAAAECBAgQIFCagOiQTC86JIMbR4AAAQIECBAgQIAAAQKlCYgOyfSiQzK4cQQIECBAgAABAgQIECBQmoDokEwvOiSDG0eAAAECBAgQIECAAAECpQmIDsn0okMyuHEECBAgQIAAAQIECBAgUJqA6JBMLzokgxtHgAABAgQIECBAgAABAqUJiA7J9KJDMrhxBAgQIECAAAECBAgQIFCagOiQTC86JIMbR4AAAQIECBAgQIAAAQKlCYgOyfSiQzK4cQQIECBAgAABAgQIECBQmoDokEwvOiSDG0eAAAECBAgQIECAAAECpQmIDsn0okMyuHEECBAgQIAAAQIECBAgUJqA6JBMLzokgxtHgAABAgQIECBAgAABAqUJiA7J9KJDMrhxBAgQIECAAAECBAgQIFCagOiQTC86JIMbR4AAAQIECBAgQIAAAQKlCYgOyfSiQzK4cQQIECBAgAABAgQIECBQmoDokEwvOiSDG0eAAAECBAgQIECAAAECpQmIDsn0okMyuHEECBAgQIAAAQIECBAgUJqA6JBMLzokgxtHgAABAgQIECBAgAABAqUJiA7J9KJDMrhxBAgQIECAAAECBAgQIFCagOiQTC86JIMbR4AAAQIECBAgQIAAAQKlCYgOyfSiQzK4cQQIECBAgAABAgQIECBQmoDokEwvOiSDG0eAAAECBAgQIECAAAECpQmIDsn0okMyuHEECBAgQIAAAQIECBAgUJqA6JBMLzokgxtHgAABAgQIECBAgAABAqUJiA7J9KJDMrhxBAgQIECAAAECBAgQIFCagOiQTC86JIMbR4AAAQIECBAgQIAAAQKlCYgOyfSiQzK4cQQIECBAgAABAgQIECBQmoDokEwvOiSDG0eAAAECBAgQIECAAAECpQmIDsn0okMyuHEECBAgQIAAAQIECBAgUJqA6JBMLzokgxtHgAABAgQIECBAgAABAqUJiA7J9KJDMrhxBAgQIECAAAECBAgQIFCagOiQTC86JIMbR4AAAQIECBAgQIAAAQKlCYgOyfSiQzK4cQQIECBAgAABAgQIECBQmoDokEwvOiSDG0eAAAECBAgQIECAAAECpQmIDsn0okMyuHEECBAgQIAAAQIECBBoxwLTp8+I20ePK+7g4yssEysut1S7uhvRIXldokMyuHEECBAgQIAAAQIECBAoQeCWO8bGk/94Lvbfc6di+gWX/DJWX2XFGLTRp1t1mhkzZsTY+/4co8feH8sts0R8dfutWvXz73744v+5Jj626CLxpS03Lv7qt7ffE88890J8a/evfKTrzesPiQ7zKtVGnxMd2gjSZQgQIECAAAECBAgQIFBhgZtvHxv/ePafse8eOxSnPP/iK2OtNVaOjTf41Ec69ZXX3BIzZs5oMTrMnDkzOnToMNfrXvSzUfG7398bV/z45IgOHWLofsfFJp9bN761x4dHh3m5/pyGiw4fad0f/YcWxOjw54ceiCH/tXHsf+BhcdSxIz86jp8kQIAAAQIECBAgQIBAksDYMXfF2DGjY+Da68QWWw5u86lziw533H1f/OSyq+OttyfF8H2/Vrz9cPxpF8aBe+8SHTt0iKlTp0X37t2KtyOOHr5ncbaWosO4+x+Nq359azz97POxx1e3iSFbbFR8dsTIs+Oz6w6M6266Mzp16hSH7Pu1GHvfw/GbW++OA765U/F3F176q9h84/WL6NDSeZ546tliZo+FusdNt42J7bf+Quz3jR1b7SQ6tJps/n5gQYwOX9xsw3j4oQcKmCHbbBs9e/aaPyQ/TYAAAQIECBAgQIAAgSYLXH/tqJg0aVIx5Ylnx0ePHj3adGIjOlx9w22x1hqrFNd94M+PxdCdt45ll14iDj76+3HGCYdE506dYvi3z4hzTjkirr7+tthsk/Xiib8/Gy++PCG+/MVN4qrrfheH7v/1FqPD8y++EvsdfnL89Lx3/g+/w448Lc479chYZOG+MXT/46Jr1y5xxncPifGvvh6dOnaMW+78Q/Tv1yduuPmu4r9ts+Um8cJL42PLzT7X4nmKIHLM6XHCUfvFaisvH3sMOz4uOXdkcf3W/BEdWqPVBp8VHdoA0SUIECBAgAABAgQIECAwnwIZ0aHxdkEjHjT+jLr+tvivQRvESy9PiKefeyFGDBta/P05P74iVlh2yeJ/nj59ejz06BPx1NP/W7wB8dQz/xs7fGnzFqND41clLrvyhth0w3e+I+LOMffFXrttV/z6RiM6NK4/8BMrz1Jq/HrFwNVXijHjHor+/XoX3y/xyF+ejIX792nxPI0vrTz3osvjorOOLa5x9Innxe67fKkIEK35Izq0RqsNPrsgRge/XtEGD4ZLECBAgAABAgQIECCQKpD96xWNuLDOwFWj8YbCS6+8Ggfv89Xifi+5/Lro2qVzrDNwtbjhltHxr4lvRefOnWPRRfrHpz/5iVh/3TVbjA633vmHIiBsO3jQLLdll148Fu7Xp4gOjTcpFhuw8AeiwxqrfTw6d+4UDz7ytyI69O3bq+XzrLV6XHrFr+P0kcOLa3z7e+fH13ccEquvskKr9iQ6tIpr/j+8IEaHhsryS/aLn/18VGwyaLP5R3IFAgQIECBAgAABAgQItHOB93+nw7vRYZGF+8XI0y4s3iDo3q1rHHDEqcV3Jayx2oqx7W6Hxu67bB39+/WNsy/8eVx6/ndjqSUWazE6vPjy+Dhi5DlxwenHRM8e3ePxJ5+OVVdavvhCyUZ0OPOEQ2PRAf0/EB0+t97axd81gkUjOmy0wTotnqdnz4VEh/b4DIoO7XFrzkyAAAECBAgQIECAAIHWCTT+ycynnvnPv17xbnRo/PrDL665Oa4Y9dvo26dXbLXZ54p/keLdWHD4AUNjwCL9YvcDvhM3/fL8eOutSXHEyLOL73lofBdD402G044/OPr17V18MWTj1yYa11m4f9848ej9onevnkV0OOukw2LAwv1miw5rr7FyfPbTa80WHRpfJNnSeR574uniLYx333Q49uQfxtd2GOxNh9Y9BvmfFh3yzU0kQIAAAQIECBAgQIBA2QLv/2cnJ0+eUvzTld26dpl1tPd+Zl7/mcrJU6bG1ClTo1ev/3wRZks/+2F/19J52sLMr1e0hWIrriE6tALLRwkQIECAAAECBAgQIECgXQuIDsnrEx2SwY0jQIAAAQIECBAgQIAAgdIERIdketEhGdw4AgQIECBAgAABAgQIEChNQHRIphcdksGNI0CAAAECBAgQIECAAIHSBESHZHrRIRncOAIECBAgQIAAAQIECBAoTUB0SKYXHZLBjSNAgAABAgQIECBAgACB0gREh2R60SEZ3DgCBAgQIECAAAECBAgQKE1AdEimFx2SwY0jQIAAAQIECBAgQIAAgdIERIdketEhGdw4AgQIECBAgAABAgQIEChNQHRIphcdksGNI0CAAAECBAgQIECAAIHSBESHZHrRIRncOAIECBAgQIAAAQIECBAoTUB0SKYXHZLBjSNAgAABAgQIECBAgACB0gREh2R60SEZ3DgCBAgQIECAAAECBAgQKE1AdEimFx2SwY0jQIAAAQIECBAgQIAAgdIERIdketEhGdw4AgQIECBAgAABAgQIEChNQHRIphcdksGNI0CAAAECBAgQIECAAIHSBESHZHrRIRncOAIECBAgQIAAAQIECBAoTUB0SKYXHZLBjSNAgAABAgQIECBAgACB0gREh2R60SEZ3DgCBAgQIECAAAECBAgQKE1AdEimFx2SwY0jQIAAAQIECBAgQIAAgdIERIdketEhGdw4AgQIECBAgAABAgQIEChNQHRIphcdksGNI0CAAAECBAgQIECAAIHSBESHZHrRIRncOAIECBAgQIAAAQIECBAoTUB0SKYXHZLBjSNAgAABAgQIECBAgACB0gREh2R60SEZ3DgCBAgQIECAAAECBAgQKE1AdEimFx2SwY0jQIAAAQIECBAgQIAAgdIERIdketEhGdw4AgQIECBAgAABAgQIEChNQHRIphcdksGNI0CAAAECBAgQIECAAIHSBESHZHrRIRncOAIECBAgQIAAAQIECBAoTUB0SKYXHZLBjSNAgAABAgQIECBAgACB0gREh2R60SEZ3DgCBAgQIECAAAECBAgQKE1AdEimFx2SwY0jQIAAAQIECBAgQIAAgdIERIdketEhGdw4AgQIECBAgAABAgQIEChNQHRIphcdksGNI0CAAAECBAgQIECAAIHSBESHZHrRIRncOAIECBAgQIAAAQIECBAoTUB0SKYXHZLBjSNAgAABAgQIECBAgACB0gREh2R60SEZ3DgCBAgQIECAAAECBAgQKE1AdEimFx2SwY0jQIAAAQIECBAgQIAAgdIERIdketEhGdw4AgQIECBAgAABAgQIEChNQHRIphcdksGNI0CAAAECBAgQIECAAIHSBESHZHrRIRncOAIECBAgQIAAAQIECBAoTUB0SKYXHZLBjSNAgAABAgQIECBAgACB0gREh2R60SEZ3DgCBAgQIECAAAECBAgQKE1AdEimFx2SwY0jQIAAAQIECBAgQIAAgdIERIdketEhGdw4AgQIECBAgAABAgQIEChNQHRIphcdksGNI0CAAAECBAgQIECAAIHSBESHZHrRIRncOAIECBAgQIAAAQIECBAoTUB0SKYXHZLBjSNAgAABAgQIECBAgACB0gREh2R60SEZ3DgCBAgQIECAAAECBAgQKE1AdEimFx2SwY0jQIAAAQIECBAgQIBAhQTeeHNiPPLXJ2PD9T9ZoVM17yiiQ/NsW7yy6JAMbhwBAgQIECBAgAABAgRKELjj7j/GyT+4JPr17R1dOneOzTddP77x1W3ib39/Nn7006vi7O8dXsKp8keKDsnmokMyuHEECBAgQIAAAQIECBAoQeC2u8bFE089G/vusUO89vqbcfAxp8eIYUOja9eus6LDzJkzo0OHDh84XWv/fk63N6frZHKIDpnaESE6JIMbR4AAAQIECBAgQIAAgRYExo65K8aOGR0D114ntthycJsbvTc6NC5+3Ck/ii98/tOx5OKLxXk/uSJWWmGZuPHWu+PLX9w0hu21czF/3P2PxlW/vjWefvb52OOr28SQLTYq/n7MvQ/GldfeEi+8ND4+t97asd83doi/PflsjP7DA3HAN3cqPnPOj6+ILQd9NlZbeYUWP9+9W7c2v8d5uaDoMC9KbfgZ0aENMV2KwBwEDjlwHzYECBAgQIAAAQIE5ipw/bWjYtKkScVnnnj2lejRo2ebijWiQ+NXLLb8wufif59/KX5x9c3x/849Pia8+kYMO/LUOH7EPjHwEyvHXsNPiPNPPTKmz5gR+x1+cvz0vJHFOYYdeVqcd+qRsXD/PrHtbofGj844Jj626CJx+aib4itf2iwe+evf4+Y77onjDtu7+Pzxp10Y2w7eND655qotfr7HQt3b9P7m9WKiw7xKtdHnRIc2gnQZAnMQuPy/fxrf/c7RMXjrbRgRIECAAAECBAgQmKPA7NFhfPTo0aNNtRrRYdT1t8Wn1lqteEvh7JNHxBqrrhiPP/lMnHrOT2fFhcYbELtsv2W88OIrcdmVN8SmG366OMedY+6LvXbbLjbe4FPxgwt/Hg88/FgM2eLzscWm68fC/fvGvX96JG65c+wHosM6A1dr8fNtenOtuJjo0Aqstvio6NAWiq5BYM4CjehwwXk/iLvHPYyJAAECBAgQIECAwBwFMn+94qwL/idWXH6p2HbwoCI6vPeLJBtvKOywzebx4kvjY8y4h4rPvPtn2aUXj4X79YnGdzM88OfH45Y7xsbYPz5cvDHxzHMvxLW/uTNOPHq/4uMjRp4du35lq2hEh5Y+P2DhfqU8DaJDMrvokAxuXO0ERIfardwNEyBAgAABAgQqKfDe73R47vmX4tBjz4z/+dFJ8fRzL7QYHRZdpF8cMfKcuOD0Y6Jnj+7x+JNPx6orLR9v/uutuOePD8UXN9uwiAl7H3JiHLb/btGr50Jx6HFnxSXnHh/Tpk2PPYYdHyOP/FasuNzSLX5+9VVWKMVJdEhmFx2SwY2rnYDoULuVu2ECBAgQIECAQCUF3v9FkqecfUmssdrHY7WVl48Lf/qrOOukw4pzN9502HGbzWPN1VeK39x6d1z0s1HRt0+v4lcoGm8xvP32v+P08y+L5198JXr36lF8D0TjiycbAeL75/0s7hn3UCy1xGIx6d+TY/89d4xll1q8xc+XhSQ6JMuLDsngxtVOQHSo3crdMAECBAgQIEBggRKYPGVqTJ0yNXr1mv07JhpRoXOnTtGlS+fZ7veNNydGn949P/BPb87p89lYokOyuOiQDG5c7QREh9qt3A0TIECAAAECBAhUWEB0SF6O6JAMblztBESH2q3cDRMgQIAAAQIECFRYQHRIXo7okAxuXO0ERIfardwNEyBAgAABAgQIVFhAdEhejuiQDG5c7QREh9qt3A0TIECAAAECBAhUWEB0SF6O6JAMblztBESH2q3cDRMgQIAAAQIECFRYQHRIXo7okAxuXO0ERIfardwNEyBAgAABAgQIVFhAdPPsf+YAACAASURBVEhejuiQDG5c7QREh9qt3A0TIECAAAECBAhUWEB0SF6O6JAMblztBESH2q3cDRMgQIAAAQIECFRYQHRIXo7okAxuXO0ERIfardwNEyBAgAABAgQIVFhAdEhejuiQDG5c7QREh9qt3A0TIECAAAECBAhUWEB0SF6O6JAMblztBESH2q3cDRMgQIAAAQIECFRYQHRIXo7okAxuXO0ERIfardwNEyBAgAABAgQIVFhAdEhejuiQDG5c7QREh9qt3A0TIECAAAECBAhUWEB0SF6O6JAMblztBESH2q3cDRMgQIAAAQIECFRYQHRIXo7okAxuXO0ERIfardwNEyBAgAABAgQIVFhAdEhejuiQDG5c7QREh9qt3A0TIECAAAECBAhUWEB0SF6O6JAMblztBESH2q3cDRMgQIAAAQIECFRYQHRIXo7okAxuXO0ERIfardwNEyBAgAABAgQIVFhAdEhejuiQDG5c7QREh9qt3A0TIECAAAECBAhUWEB0SF6O6JAMblztBESH2q3cDRMgQIAAAQIECFRYQHRIXo7okAxuXO0ERIfardwNEyBAgAABAgQIVFhAdEhejuiQDG5c7QREh9qt3A0TIECAAAECBAhUWEB0SF6O6JAMblztBESH2q3cDRMgQIAAAQIECFRYQHRIXo7okAxuXO0ERIfardwNEyBAgAABAgQIVFhAdEhejuiQDG5c7QREh9qt3A0TIECAAAECBAhUWEB0SF6O6JAMblztBESH2q3cDRMgQIAAAQIECFRYQHRIXo7okAxuXO0ERIfardwNEyBAgAABAgQIVFhAdEhejuiQDG5c7QREh9qt3A0TIECAAAECBAhUWEB0SF6O6JAMblztBESH2q3cDRMgQIAAAQIECFRYQHRIXo7okAxuXO0ERIfardwNEyBAgAABAgQIVFhAdEhejuiQDG5c7QREh9qt3A0TIECAAAECBAhUWEB0SF6O6JAMblztBESH2q3cDRMgQIAAAQIECFRYQHRIXo7okAxuXO0ERIfardwNEyBAgAABAgQIVFhAdEhejuiQDG5c7QREh9qt3A0TIECAAAECBAhUWEB0SF6O6JAMblztBESH2q3cDRMgQIAAAQIECFRYQHRIXo7okAxuXO0ERIfardwNEyBAgAABAgQIVFhAdEhejuiQDG5c7QREh9qt3A0TIECAAAECBAhUWEB0SF6O6JAMblztBESH2q3cDRMgQIAAAQIECFRYQHRIXo7okAxuXO0ERIfardwNEyBAgAABAgQIVFhAdEhejuiQDG5c7QREh9qt3A0TIECAAAECBAhUWEB0SF6O6JAMblztBESH2q3cDRMgQIAAAQIECFRYQHRIXo7okAxuXO0ERIfardwNEyBAgAABAgQIVFhAdEhejuiQDG5c7QREh9qt3A0TIECAAAECBAhUWEB0SF6O6JAMblztBESH2q3cDRMgQIAAAQIECFRYQHRIXo7okAxuXO0ERIfardwNEyBAgAABAgQIVFhAdEhejuiQDG5c7QREh9qt3A0TIECAAAECBAhUWEB0SF6O6JAMblztBESH2q3cDRMgQIAAAQIECFRYQHRIXo7okAxuXO0ERIfardwNEyBAgAABAgQIVFhAdEhejuiQDG5c7QQuOO+suPTii2LcQ4/X7t7dMAECBAgQIECAAIGqCYgOyRtZUKPDsov3iYMOOSLWXHOtZFHjCPxH4MUXno9vH3Vo8Rff3Gf/OOHkM/AQIECAAAECBAgQKFVg5syZxfwOHTqUeo6yhosOyfILanT47neOip//7JLo3KVzsqhxBP4jMG3atHhr4sTiLzbe9Atxxa9uwEOAAAECBAgQIECgNIGbfjcmLrn8umL+bjsOiW2+uElpZ3l38I233h3rrr1aLL7YgJSziA4pzP8ZsqBGh2RG4wjMUWCrzTaMxx/7S/zg3Iti26/sRIoAAQIECBAgQIBAKQJPP/t8fOfUH8V5px4ZM2bMiH0OPSlGHvGtWGO1j3/k8zTemnj/GxPTp8+ITp06FtdszOnY8Z3/eU5/hu5/XJx5wqGx6ID+83yOlubO6w+LDvMq1UafEx3aCNJlCMxBwBdJejQIECBAgAABAgTmRWDsmLti7JjRMXDtdWKLLQfPy4+06jM//9VN0bNH99h28KDi56685pb418S3Yq/dtosRI8+Oz647MK676c7o1KlTHLLv12KtNVaOy0fdFDfcPDqmTpsW+++5Ywza6DPxxFPPFj/bY6HucdNtY2L7rb8Q+31jx1nXvPHW0TFg4X7Rp0+vWKR/39jkc+vG6D88EAd8853/A9w5P74ithz02ejevXt855QL4oWXxseii/SPZZdePE79zkHFZ1o7tzUQokNrtNrgs6JDGyC6BIG5CIgOHo8qCXT580PR88c/rNKRnIUAAQIECBD4P4GVrh0VL0yaVPxvTzw7Pnr06NGmNqeff1l84fOfiXXXXr24biME/H7Mn+LYw/aKxtsGXbt2iTO+e0iMf/X16NSxY6yw3FJxz7iHis/PjJmxzyEnxZknHhovj38thh9zepxw1H6x2srLxx7Djo9Lzh0ZU6ZOje+ddXGcffKIuPPu++JPD/01jjxojxh3/6Nx8x33xHGH7V3MPf60C2PbwZvGOgNXi8ZbEVvtdEDccPk50aVL51lvRbRm7iIL922Vk+jQKq75/7DoMP+GrkBgbgKig+ejSgKLrb9WzOjbL6atulqVjuUsBAgQIECAQEQ0Ozqcdu6lsemGn471112z8H7o0b/Fr3/7+yIGNKLDiGFDY+AnVp5tF6+9/mY8/uQz8cqE1+LqG26P4fvuGp07d45zL7o8Ljrr2OKzR594Xuy+y5di+owZ8aNLfllEhztG/zGe/MdzxRsQ9/7pkbjlzrEtRofGr0n81w77x62jfvSR5zbCR2v+iA6t0WqDz4oObYDoEgTmIiA6eDyqJLDoJuvH218bGm/tc0CVjuUsBAgQIECAQEQ0+9cr/vuXN0b3bl1jxy9vUXj/5ta747l/vhTf2uMrRXQ444RDYrEBC8/axYRX34gDjzotthsyKJZcfNEYdf1tsdvOQ6Jr165x6RW/jtNHDi8+++3vnR9f33FI8dZD4/ONa6684rJx+LDdii+HbLzxcO1v7owTj96v+HzjVzl2/cpWxZsOkydPia13PXi26NDauauvskKrnh/RoVVc8/9h0WH+DV2BwNwERAfPR5UERIcqbcNZCBAgQIBArsBTT/9vHHvyBXH+aUcW4aARCA7+1q7xyTVXKaLD+7/M8ba7xsV9D/6l+BWJiRPfLr54csSBQ+cYHZZecrE44IhT43vHDouePRaK/n17F18y+dw/X4xDjzsrLjn3+Jg2bXrx6xgjj/xWER0a/9rbF3c+MC7/8cnFf+vdq2fc+6c/t2qu6JD7HLV6mujQajI/QKBVAqJDq7h8uMkCokOTgV2eAAECBAhUXKDxtsJlV77zz7h/9StbxS7bbVn8z43ocNZJhxVfAPnun0n/nlyEiX59eke/vr3j6eeeL74Msnu3bsU/u/numw7HnvzD+NoOg2OF5ZaMPQ/8biy3zBLR+Nl/vvBy/OiMY4ovk/z+eT8rvh9iqSUWK/5b40spP7POGsWoS35+XfEFlo3vlDj52GGx9JIfa9Vc0aHiD53oUPEFOV67FxAd2v0KF6gbEB0WqHW6GQIECBAg8JEEJk+ZGh06RHTt0mXWz8/tn6Cc+Nak6NVzoQ+d1YgZvXr2KP41i8af8y++svg1iy2/sEHxv7/x5sTo07vnB/6JzcZ/+/fkydGta9fZ/tu8zv3Qg73vA369orVi8/l50WE+Af04gQ8REB08IlUSEB2qtA1nIUCAAAECC5bAY088HWf88LLieyOmTp0Wiw7oX3x5ZLduXSt1o6JD8jpEh2Rw42onIDrUbuWVvmHRodLrcTgCBAgQILBACDSCQ8eOHaNTp46VvB/RIXktokMyuHG1ExAdarfySt+w6FDp9TgcAQIECBAgkCAgOiQgv3eE6JAMblztBESH2q280jcsOlR6PQ5HgAABAgQIJAiIDgnIokMysnG1FhAdar3+yt286FC5lTgQAQIECBAgkCwgOiSDe9MhGdy42gmIDrVbeaVvWHSo9HocjgABAgQIEEgQEB0SkL3pkIxsXK0FRIdar79yNy86VG4lDkSAAAECBAgkC4gOyeDedEgGN652AqJD7VZe6RsWHSq9HocjQIAAAQIEEgREhwRkbzokIxtXawHRodbrr9zNiw6VW4kDESBAgAABAskCokMyuDcdksGNq52A6FC7lVf6hkWHSq/H4QgQIECAAIEEAdEhAdmbDsnIxtVaQHSo9ford/OiQ+VW4kAECBAgQIBAsoDokAzuTYdkcONqJyA61G7llb5h0aHS63E4AgQIECBAIEFAdEhA9qZDMrJxtRYQHWq9/srdvOhQuZU4EAECBAgQIJAsIDokg3vTIRncuNoJiA61W3mlb1h0qPR6HI4AAQIECBBIEBAdEpC96ZCMbFytBUSHWq+/cjcvOlRuJQ5EgAABAgQIJAuIDsng3nRIBjeudgKiQ+1WXukbFh0qvR6HI0CAAAECBBIERIcEZG86JCMbV2sB0aHW66/czYsOlVuJAxEgQIAAAQLJAqJDMrg3HZLBjaudgOhQu5VX+oZFh0qvx+EIECBAgACBBAHRIQHZmw7JyMbVWkB0qPX6K3fzokPlVuJABAgQIECAQLKA6JAM7k2HZHDjaicgOtRu5ZW+YdGh0utxOAIECBAgQCBBQHRIQPamQzKycbUWEB1qvf7K3bzoULmVOBABAgQIECCQLCA6JIN70yEZ3LjaCYgOtVt5pW9YdKj0ehyOAAECBAgQSBAQHRKQvemQjGxcrQVEh1qvv3I3LzpUbiUORIAAAQIECCQLiA7J4N50SAY3rnYCokPtVl7pGxYdKr0ehyNAgAABAgQSBESHBGRvOiQjG1drAdGh1uuv3M2LDpVbiQMRIECAAAECyQKiQzK4Nx2SwY2rnYDoULuVV/qGRYdKr8fhCBAgQIAAgQQB0SEB2ZsOycjG1VpAdKj1+it386JD5VbiQAQIECBAgECygOiQDO5Nh2Rw42onIDrUbuWVvmHRodLrcTgCBAgQIEAgQUB0SED2pkMysnG1FhAdar3+yt286FC5lTgQAQIECBAgkCwgOiSDe9MhGdy42gmIDrVbeaVvWHSo9HocjgABAgQIEEgQEB0SkL3pkIxsXK0FRIdar79yNy86VG4lDkSAAAECBAgkC4gOyeDedEgGN652AqJD7VZe6RsWHSq9HocjQIAAAQIEEgREhwRkbzokIxtXawHRodbrr9zNiw6VW4kDESBAgAABAskCokMyuDcdksGNq52A6FC7lVf6hkWHSq/H4QgQIECAAIEEAdEhAdmbDsnIxtVaQHSo9ford/OiQ+VW4kAECBAgQIBAsoDokAzuTYdkcONqJyA61G7llb5h0aHS63E4AgQIECBAIEFAdEhA9qZDMrJxtRYQHWq9/srdvOhQuZU4EAECBAgQIJAsIDokg3vTIRncuNoJiA61W3mlb1h0qPR6HI4AAQIECBBIEBAdEpC96ZCMbFytBUSHWq+/cjcvOlRuJQ5EgAABAgQIJAuIDsng3nRIBjeudgKiQ+1WXukbFh0qvR6HI0CAAAECBBIERIcEZG86JCMbV2sB0aHW66/czYsOlVuJAxEgQIAAAQLJAqJDMrg3HZLBjaudgOhQu5VX+oZFh0qvx+EIECBAgACBBAHRIQHZmw7JyMbVWkB0qPX6K3fzokPlVuJABAgQIECAQLKA6JAM7k2HZHDjaicgOtRu5ZW+YdGh0utxOAIECBAgQCBBQHRIQPamQzKycbUWEB1qvf7K3bzoULmVOBABAgQIECCQLCA6JIN70yEZ3LjaCYgOtVt5ZW941FVXxFkH7xdLL71M/L/b74levXpX9qwORoAAAQIECBBoloDo0CzZOVxXdEgGN652Ao3ocNr3RsaECeNrd+9uuFoCM2fOnHWgSy67MrYc/KVqHdBpCBAgQIAAAQIJAqJDAvJ7R4gOyeDG1VLgnrvvquV9u+lqCfzyiv+Oq678eXTt1CluvmtcrLLq6tU6oNMQIECAAAECBBIERIcEZNEhGdk4AgQIVETgXxusE9132TW6HDyiIidyDAIECBAgQIBAroDokOsd3nRIBjeOAAECJQr4IskS8Y0mQIAAAQIVEpg+fUZ06tRxvk404dU3Yuq0qbH4YgOK6zSuefvoccX//PEVlokVl1tqrn8/X8Pn44dFh/nA+yg/Kjp8FDU/Q4AAgfYpIDq0z705NQECBAgQaCuB+x9+LM6+8Ocx6d+TY9WVloujDvpG9OrV4yNd/sprbolXJrwWw/baufj5GTNmxNj7/hyjx94fyy2zRHx1+63m+vcfaWgb/JDo0AaIrbmE6NAaLZ8lQIBA+xYQHdr3/pyeAAECBAjMj8Crr78Zw444NU75zoGx3NJLxKVX/DqWW2bJGLTRp+fnsh/42UaMmDFzxqzo8O4H5vT3bTp8Hi4mOswDUlt+RHRoS03XIkCAQLUFRIdq78fpCBAgQKDeAmPH3BVjx4yOgWuvE1tsObjNMW65Y2z87e/PznozofEvW3Xo0KGYc8fd98VPLrs63np7Ugzf92uzQsSIkWfHZ9cdGNfddGd06tQpDtn3a7HWGivHSWdeHH/7+zOx7eBBsf3WX5jtrK2JDpePuiluuHl0TJ02Lfbfc8cYtNFn2vy+339B0aHpxLMPEB2SwY0jQIBAiQKiQ4n4Rn9AoN/w/SOmTyNDgAABAv8nsNK1o+KFSZOK/+2JZ8dHjx4f7dce5gTaeLNhsQELx+AtNprtI08/90IcfPT344wTDonOnTrF8G+fEeecckQsv8wSMXT/46Jr1y5xxncPifGvvh6dOnaMFZZbKhrB4pob74iJb70dQ3fe+iNHh3vGPRTrrr16zIyZsc8hJ8WZJx4aiy7Sv6nPhOjQVN4PXlx0SAY3jgABAiUKiA4l4hs9m0DXe++JAUM2j7d3/lpEx3f+r2z+ECBAoO4CzY4OF176q1h5xWVjs43Xi4suHRVjxj0YW232uejSuXM0wsOIYUOLFZzz4ytihWWXjG222qSIDo2/H/iJlT+wnmtuvD0mvjUpdttpyEeODq+9/mY8/uQzxXdDXH3D7TF8311j7TVWaeqjIDo0lVd0SOY1jgABApUSEB0qtY5aH+bd6PD8i29GdO5caws3T4AAgXcFmv3rFY1IMH7C67H30O2Lf2niD/c9HI/89e/Rt2+veOmVV+Pgfb5aHOWSy6+Lrl06x9d3HFJEh8YbEI03JN7/Z36jQ+NfvzjwqNNiuyGDYsnFF41R198Wu+08JNYZuFpTHwrRoam8okMyr3EECBColIDoUKl11PowokOt1+/mCRAoSeCllyfEPoeeFBee+e1Y4mMDovEdD/945vnYaIN1YuRpF8ZFZx0b3bt1jQOOODX2+8aOsd6n1iiiw5knHBqLDvjgrzzMb3S47a5xcd+Df4kjD9ojJk58uzjbiAOHig5zeD46zGz8Uks7/OPXK9rh0hyZAAECH1FAdPiIcH6szQVEhzYndUECBAjMk0DjCyMvuOSX0bdPr1hqicVipy9vEWus9vH4xTU3xxWjflv8feNXLhr/3GXjSyYb0eGskw6LAQv3m+36I79/UTz62N9jypSpsf66a8Yxh3wz3nhzYhwx8ux48eUJxfdALNyvT5x2/MHFdVr6+27duhZvOvTr0zv69e0dTz/3fBzwzZ1EB9Fhnp5lHyJAgACBCgqIDhVcSk2PJDrUdPFumwCBSghMmz69iAU9Fuo+23kmT54S0aFDdOvaZdbfv/dfuGjW4RvfC9Gr50LNuvwHruvXK9Ko3xnkTYdkcOMIECBQooDoUCK+0bMJiA4eCAIECBAoS0B0SJYXHZLBjSNAgECJAqJDifhGiw6eAQIECBCohIDokLwG0SEZ3DgCBAiUKCA6lIhvtOjgGSBAgACBSgiIDslrEB2SwY0jQIBAiQKiQ4n4RosOngECBAgQqISA6JC8BtEhGdw4AgQIlCggOpSIb7To4BkgQIAAgUoIiA7JaxAdksGNI0CAQIkCokOJ+EaLDp4BAgQIEKiEgOiQvAbRIRncOAIECJQoIDqUiG+06OAZIECAAIFKCIgOyWsQHZLBjSNAgECJAqJDifhGiw6eAQIECBCohIDokLwG0SEZ3DgCBAiUKCA6lIhvtOjgGSBAgACBSgiIDslrEB2SwY0jQIBAiQKiQ4n4RosOngECBAgQqISA6JC8BtEhGdw4AgQIlCggOpSIb7To4BkgQIAAgUoIiA7JaxAdksGNI0CAQIkCokOJ+EaLDp4BAgQIEKiEgOiQvAbRIRncOAIECJQoIDqUiG+06OAZIECAAIFKCIgOyWsQHZLBjSNAgECJAqJDifhGiw6eAQIECBCohIDokLwG0SEZ3DgCBAiUKCA6lIhvtOjgGSBAgACBSgiIDslrEB2SwY0jQIBAiQKiQ4n4RosOngECBAgQqISA6JC8BtEhGdw4AgQIlCggOpSIb7To4BkgQIAAgUoIiA7JaxAdksGNI0CAQIkCokOJ+EaLDp4BAgQIEKiEgOiQvAbRIRncOAIECJQoIDqUiG+06OAZIECAAIFKCIgOyWsQHZLBjSNAgECJAqJDifhGiw6eAQIECBCohIDokLwG0SEZ3DgCBAiUKCA6lIhvtOjgGSBAgACBSgiIDslrEB2SwY0jQIBAiQKiQ4n4RosOngECBAgQqISA6JC8BtEhGdw4AgQIlCggOpSIb7To4BkgQIAAgUoIiA7JaxAdksGNI0CAQIkCokOJ+EaLDp4BAgQIEKiEgOiQvAbRIRncOAIECJQoIDqUiG+06OAZIECAAIFKCIgOyWsQHZLBjSNAgECJAqJDifhGiw6eAQIECBCohIDokLwG0SEZ3DgCBAiUKCA6lIhvtOjgGSBAgACBSgiIDslrEB2SwY0jQIBAiQKiQ4n4RosOngECBAgQqISA6JC8BtEhGdw4AgQIlCggOpSIb7To4BkgQIAAgUoIiA7JaxAdksGNI0CAQIkCokOJ+EaLDp4BAgQIEKiEgOiQvAbRIRncOAIECJQoIDqUiG+06OAZIECAAIFKCIgOyWsQHZLBjSNAgECJAqJDifhGiw6eAQIECBCohIDokLwG0SEZ3DgCBAiUKCA6lIhvtOjgGSBAgACBSgiIDslrEB2SwY0jQIBAiQKiQ4n4RosOngECBAgQqISA6JC8BtEhGdw4AgQIlCggOpSIb7To4BkgQIAAgUoIiA7JaxAdksGNI0CAQIkCokOJ+EaLDp4BAgQIEKiEgOiQvAbRIRncOAIECJQoIDqUiG+06OAZIECAAIFKCIgOyWsQHZLBjSNAgECJAqJDifhGzyYw/sbrYvzuX42V//lqdOrWnQ4BAgQIEEgTEB3SqN8ZJDokgxtHgACBEgUWW3+tmLrOujF5081KPIXRdReYNHVKrHnMiHh10qTY6PObxJXX3FR3EvdPgAABAokCokMituiQjG0cAQIEShbo+d+XRu9Tvhszu/u/LJe8ilqPf2batPj48/8sDJZYYsm4789P1trDzRMgQIBAroDokOvtTYdkb+MIECBAgACBiJP32T0evfqq2OKUM2OPvfdDQoAAAQIE0gREhzTqdwb59YpkcOMIECBAgACB6HrvPTFgyObx/ItvRnTuTIQAAQIECKQJiA5p1KJDMrVxBAgQIECAwP8JiA4eBQIECBAoS0B0SJb3pkMyuHEECBAgQICANx08AwQIECBQmoDokEwvOiSDG0eAAAECBAiIDp4BAgQIEChNQHRIphcdksGNI0CAAAECBEQHzwABAgQIlCYgOiTTiw7J4MYRIECAAAECooNngAABAgRKExAdkulFh2Rw4wgQIECAAAHRwTNAgAABAqUJiA7J9KJDMrhxBAgQIECAgOjgGSBAgACB0gREh2R60SEZ3DgCBAgQIEBAdPAMECBAoGSBmTNnRocOHYpTvPHmxHjkr0/Ghut/ssVTTXj1jZg6bWosvtiAkk/dNuNFh7ZxnOeriA7zTOWDBAgQIECAQBsJdL33nhgwZPN4/sU3Izp3bqOrugwBAgQIfJjAn//6ZJz/k1/ElKnTYtmlF48RBwyNf774Svzop1fF2d87vMUfv/KaW+KVCa/FsL12/rDLt4v/Ljokr0l0SAY3jgABAgQIEPCmg2eAAAECJQhMfGtSfOPA4+PMEw4tgsPPfnF99Fioe6y1xiqzosN734D4sCO++9nW/Mycrtnaa8zp8/NyHdHhwzbbxv9ddGhjUJcjQIAAAQIEPlTAmw4fSuQDBAjUUGDsmLti7JjRMXDtdWKLLQe3ucBd99wf199yV5w+cnhx7enTZ0SnTh3j8SefifN+ckWstMIyceOtd8eXv7jprLcaTjrz4vjb35+JbQcPiu23/kLxc//9yxvjjw88Gs8890IM3nyj+N1d98YB39wpevfqGTfffk+89sa/4sE/Px77DN0+dthm8+LXOC4fdVPccPPomDptWuy/544xaKPPxBNPPRuNtyga4eOm28YU19/vGzvG08+9EN855YJi1rRp06Nz505x6ncOiiUXXzTG3f9oXPXrW+PpZ5+PPb66TQzZYqPic3P6+5YQRYc2f7TmfkHRIRncOAIECBAgQMCbDp6BSgl0mPzv6Hv4QZU6k8PUU2Cla0fFC5MmFTf/xLPjo0ePHm0KccWo38bMmBm7fuWLs123ER2GHXlqHD9inxj4iZVjr+EnxPmnHhkfW2yRaLw5cM2Nd8TEt96OoTtvXfzc+RdfGf369i5iwKOPPRXbDRlUxIbPf3adOP38y+LHPzguZs6cEYcee1Ycd/jesfoqK8Q94x6KdddevZi/MDvD5AAAFVBJREFUzyEnxZknHhovj38thh9zepxw1H6x2srLxx7Djo9Lzh0Z/fv1jn9PnvJO4Ljyxnjh5fHF2V54aXzsd/jJ8dPzRhb/bdiRp8V5px4Zk6dMafHvF1m4b4t+okObPlYffjHR4cONfIIAAQIECBBoWwFvOrStp6vNn0D/Pb8eXf7ySEz5zHrzdyE/TWA+BZodHRpvKDTebGgpOpx6zk9n/T/mjzvlR7HL9lvGGquuWNzRNTfeHo1fzdhtpyGzosN6n1ozJrz6ekx47Y0iJlx57S2x+cbrxa2/vze+e+S+7wSDX94YXbp0jl222zJee/3N4o2KxndDXH3D7TF8312jc+fOce5Fl8dFZx1bfP7oE8+L3Xf5UhEgGn8ab1Oc95NfxAWnHxO9ei4Uv/v9vXHZlTfEpht+uvjvd465L/babbuYMmVqi3+/8QafEh3m85lskx8XHdqE0UUIECBAgACBVgiIDq3A8tGmC/Tfe/eY0bdvvHHGuU2fZQCBuQk0+9crbrljbNxyxx/ijBMOme0YjRjw3i+SPP60C4tfixi4+kpzjA7rr7tmjJ/werz2xpuxzsDViuiwxSbrx11j74+jh+9Z/Nyvfv27ePNfbxW/rnHgUacVb0Q0fkVi1PW3xW47D4muXbvGpVf8etave3z7e+fH13ccUrwZ8dLLE2LYUafFKccdWPzaR+PPrXf+IcaMe6j4VY93/zS+m+JPD/6lxb9fuF8f0aEK/7+c6FCFLTgDAQIECBCol4DoUK99V/1uRYeqb8j52krg9Tf+Fbsf8J0499QjYrmll4gx9z4YfXr3LP4f/20VHS762aj48Q+OjU4dO8WIkWfHDl/aLCZPmRr3PfiXOPKgPWLixLdjn0NPihEHDp1jdFhx+aXjkG+fEV/aauP44mYbzrr9F18eH0eMPKd486Fnj+7x+JNPx6orLR8vvTKhxb9/958Efb+fX69oqydqHq8jOswjlI8RIECAAAECbSYgOrQZpQu1gYDo0AaILtFuBBpvCpx94c+j8RZAzx4LFd+50PiVhwt/+qs466TDivtovOmw4zabx5qrrxQjv39RPPrY34tfYWi83XDMId8svtNhg8+sFS+PfzVef/1f8cmBq8ZV190am228XvFrDgv37xtP/uO52Hzj9WPvodsV0aHxpkO/Pr2L74J4+rnniy+e7N6tW1xy+XWz3nQ49uQfxtd2GBx/f/p/45yLLo8lPjZgluvO221ZfGnkb269Oxpho2+fXsWcE4/er/gCyzn9fUuLER2SH1fRIRncOAIECBAgQMAXSXoGKiUgOlRqHQ6TIDBt+vT496TJ0atX231RZeMLJ+/548Mxeuz9MWLY7vHvyZOLqPHeP43vhWh8N8OH/Wlcq6U/77650IgYU6dM/cD55/T377+W6PBhG2jj/y46tDGoyxEgQIAAAQIfKuBNhw8l8oFEAdEhEduoBVqg8RZFIzocdfA3Kn2fokPyekSHZHDjCBAgQIAAAW86eAYqJSA6VGodDkOg6QKiQ9OJZx8gOiSDG0eAAAECBAiIDp6BSgmIDpVah8MQaLqA6NB0YtEhmdg4AgQIECBA4H0Cfr3CI1ElAdGhSttwFgLNFxAdmm882wRvOiSDG0eAAAECBAh408EzUCkB0aFS63AYAk0XEB2aTjz7ANEhGdw4AgQIECBAQHTwDFRKQHSo1DochkDTBUSHphOLDsnExhEgQIAAAQLvE/DrFR6JKgmIDlXahrMQaL6A6NB849kmeNMhGdw4AgQIECBAwJsOnoFKCYgOlVqHwxBouoDo0HTi2QeIDsngxhEgQIAAAQKig2egUgKiQ6XW4TAEmi4gOjSdWHRIJjaOAAECBAgQeJ+AX6/wSFRJQHSo0jachUDzBUSH5hvPNsGbDsngxhEgQIAAAQLedPAMVEpAdKjUOhyGQNMFRIemE88+QHRIBjeOAAECBAgQEB08A5USEB0qtQ6HIdB0AdGh6cSiQzKxcQQIECBAgMD7BPx6hUeiSgKiQ5W24SwEmi8gOjTfeLYJ3nRIBjeOAAECBAgQ8KaDZ6BSAqJDpdbhMASaLiA6NJ149gGiQzK4cQQIECBAgIDo4BmolIDoUKl1OAyBpguIDk0nFh2SiY0jQIAAAQIE3ifg1ys8ElUSEB2qtA1nIdB8AdGh+cazTfCmQzK4cQQIECBAgIA3HTwDlRIQHSq1Doch0HQB0aHpxLMPEB2SwY0jQIAAAQIERAfPQKUERIdKrcNhCDRdQHRoOrHokExsHAECBAgQIPA+Ab9e4ZGokoDoUKVtOAuB5guIDs03nm2CNx2SwY0jQIAAAQIEvOngGaiUgOhQqXU4DIGmC4gOTSeefYDokAxuHAECBAgQICA6eAYqJSA6VGodDkOg6QKiQ9OJRYdkYuMIECBAgACB9wn49QqPRJUERIcqbcNZCDRfQHRovvFsE7zpkAxuHAECBAgQIOBNB89ApQT67b17zOzbN94449xKncthCBBojoDo0BzXOV5VdEgGN44AAQIECBCYFR1e+uMjEZ07ESFQmsAtv78jhg7fP/p17RoXX3V9bLDh50s7i8EECOQIiA45zrOmiA7J4MYRIECAAAECRXToe/Th0eXhB2kQKFXg6Ig49f9OcMxxJ8YBBx9W6nkMJ0Cg+QKiQ/ONZ5sgOiSDG0eAAAECBAgUAp2ee4YEgdIFvOlQ+gocgEC6gOiQTC46JIMbR4AAAQIECBAgUCkB3+lQqXU4DIGmC4gOTSeefYDokAxuHAECBAgQIECAQKUE/OsVlVqHwxBouoDo0HRi0SGZ2DgCBAgQIECAAIEKC4gOFV6OoxFogoDo0ATUuV3Smw7J4MYRIECAAAECBAhUSkB0qNQ6HIZA0wVEh6YTzz5AdEgGN44AAQIECBAgQKBSAqJDpdbhMASaLiA6NJ1YdEgmNo4AAQIECBAgQKDCAqJDhZfjaASaICA6NAF1bpf0pkMyuHEECBAgQIAAAQKVEhAdKrUOhyHQdAHRoenEsw8QHZLBjSNAgAABAgQIEKiUgOhQqXU4DIGmC4gOTScWHZKJjSNAgAABAgQIEKiwgOhQ4eU4GoEmCIgOTUCd2yW96ZAMbhwBAgQIECBAgEClBESHSq3DYQg0XUB0aDrx7ANEh2Rw4wgQIECAAAECBColIDpUah0OQ6DpAqJD04lFh2Ri4wgQIECAAAECBCosIDpUeDmORqAJAqJDE1DndklvOiSDG0eAAAECBAgQIFApAdGhUutwGAJNFxAdmk48+wDRIRncOAIECBAgQIAAgUoJiA6VWofDEGi6gOjQdGLRIZnYOAIECBAgQIAAgQoLiA4VXo6jEWiCgOjQBNS5XdKbDsngxhEgQIAAAQIECFRKQHSo1DochkDTBUSHphPPPkB0SAY3jgABAgQIECBAoFICokOl1uEwBJouIDo0nVh0SCY2jgABAgQIECBAoMICokOFl+NoBJogIDo0AXVul/SmQzK4cQQIECBAgAABApUSEB0qtQ6HIdB0AdGh6cSzDxAdksGNI0CAAAECBAgQqJSA6FCpdTgMgaYLiA5NJxYdkomNI0CAAAECBAgQqLCA6FDh5TgagSYIiA5NQJ3bJb3pkAxuHAECBAgQIECAQKUERIdKrcNhCDRdQHRoOvHsA0SHZHDjCBAgQIAAAQIEKiUgOlRqHQ5DoOkCokPTiUWHZGLjCBAgQIAAAQIEKiwgOlR4OY5GoAkCokMTUOd2SW86JIMbR4AAAQIECBAgUCkB0aFS63CYJgrccsfYePIfz8X+e+5UTLngkl/G6qusGIM2+nQTp1bv0qJD8k5Eh2Rw4wgQIECAAAECBColIDpUah0O00SBm28fG/949p+x7x47FFPOv/jKWGuNlWPjDT4116kzZ86MDh06xLv/33k5Yms+Oy/Xa8vPiA5tqTkP1xId5gHJRwgQIECAAAECBBZYAdFhgV1tu7uxsWPuirFjRsfAtdeJLbYc3Obnn1N06Nend4z+wwNxwDffeQPinB9fEVsO+mz88YG/xB8feDSeee6FGLz5RvG7u+4tPrPphp+Oy0fdFDfcPDqmTpsW+++5Ywza6DPxxFPPxpXX3BI9FuoeN902Jrbf+gux3zd2bPP7mN8Lig7zK9jKnxcdWgnm4wQIECBAgAABAguUQP89vx5d/vJITPnMegvUfbmZ9iew0rWj4oVJk4qDP/Hs+OjRo0eb3kQjOlx9w22x1hqrFNd94M+PxdCdt47u3brFzXfcE8cdtnfx98efdmFsO3jTGHPvQ9Gvb+/o3LlTPPrYU7HdkEFx8+33xNHD94x7xj0U6669esyMmbHPISfFmSceGi+Pfy2GH3N6nHDUfrHaysvHHsOOj0vOHRmLLNy3Te9jfi8mOsyvYCt/XnRoJZiPEyBAgAABAgQILFACHSb/O/oeftACdU9upn0KZESHsfc9HF/+4iYF0Kjrb4v/GrRBdOvaNW65c2yL0WG9T60ZE159PSa89kYRGa689pYYecS34rXX34zHn3wmXpnwWlx9w+0xfN9do3PnznHuRZfHRWcdW1z/6BPPi913+VIRIKr0R3RI3obokAxuHAECBAgQIECAAAECBFoQyP71isavUawzcNXo2WOhuPY3d8aJR+9XnGrEyLNj169sVbzpsP66a8b4Ca/Ha2+8GesMXK2IDgfutUsceNRpxZsPSy6+aBEvdtt5SHTt2jUuveLXcfrI4cV1vv298+PrOw6J1VdZoVL7Fh2S1yE6JIMbR4AAAQIECBAgQIAAgRIE3v+dDu9GhxWWXTIOPe6suOTc42PatOnFr0WMPPJbc4wOn//sOnHfg3+JIw/aIyZOfDv2OfSkGHHgUNGhyTvtMLPx9Zzt8I/o0A6X5sgECBAgQIAAAQIECBBopUDjn8x86pn//OsV70aHRkT4/nk/K76nYaklFotJ/55cfDnkvX96JDb4zFrx8vhX4/XX/xWfHLhqXHXdrTHiwN2LNx0aX0DZ+M6Hp597vviCycZ3Q1xy+XWz3nQ49uQfxtd2GOxNh1buaU4fFx3aCNJlCBAgQIAAAQIECBAgQKD5Au//Zy3feHNi9Onds/jnMef0570/M/GtSdGr50LNP2gbT/DrFW0M+mGX86bDhwn57wQIECBAgAABAgQIECCwoAiIDsmbFB2SwY0jQIAAAQIECBAgQIAAgdIERIdketEhGdw4AgQIECBAgAABAgQIEChNQHRIphcdksGNI0CAAAECBAgQIECAAIHSBESHZHrRIRncOAIECBAgQIAAAQIECBAoTUB0SKYXHZLBjSNAgAABAgQIECBAgACB0gREh2R60SEZ3DgCBAgQIECAAAECBAgQKE1AdEimFx2SwY0jQIAAAQIECBAgQIAAgdIERIdketEhGdw4AgQIECBAgAABAgQIEChNQHRIphcdksGNI0CAAAECBAgQIECAAIHSBESHZHrRIRncOAIECBAgQIAAAQIECBAoTUB0SKYXHZLBjSNAgAABAgQIECBAgACB0gREh2R60SEZ3DgCBAgQIECAAAECBAgQKE1AdEimFx2SwY0jQIAAAQIECBAgQIAAgdIERIdketEhGdw4AgQIECBAgAABAgQIEChNQHRIphcdksGNI0CAAAECBAgQIECAAIHSBESHZHrRIRncOAIECBAgQIAAAQIECBAoTUB0SKYXHZLBjSNAgAABAgQIECBAgACB0gREh2R60SEZ3DgCBAgQIECAAAECBAgQKE1AdEimFx2SwY0jQIAAAQIECBAgQIAAgdIERIdketEhGdw4AgQIECBAgAABAgQIEChNQHRIphcdksGNI0CAAAECBAgQIECAAIHSBESH0ugNJkCAAAECBAgQIECAAAECBKoo0GHmzJkzq3gwZyJAgAABAgQIECBAgAABAgTat4Do0L735/QECBAgQIAAAQIECBAgQKCyAqJDZVfjYAQIECBAgAABAgQIECBAoH0LiA7te39OT4AAAQIECBAgQIAAAQIEKisgOlR2NQ5GgAABAgQIECBAgAABAgTat4Do0L735/QECBAgQIAAAQIECBAgQKCyAqJDZVfjYAQIECBAgAABAgQIECBAoH0LiA7te39OT4AAAQIECBAgQIAAAQIEKisgOlR2NQ5GgAABAgQIECBAgAABAgTat4Do0L735/QECBAgQIAAAQIECBAgQKCyAqJDZVfjYAQIECBAgAABAgQIECBAoH0LiA7te39OT4AAAQIECBAgQIAAAQIEKisgOlR2NQ5GgAABAgQIECBAgAABAgTat4Do0L735/QECBAgQIAAAQIECBAgQKCyAqJDZVfjYAQIECBAgAABAgQIECBAoH0LiA7te39OT4AAAQIECBAgQIAAAQIEKisgOlR2NQ5GgAABAgQIECBAgAABAgTat4Do0L735/QECBAgQIAAAQIECBAgQKCyAqJDZVfjYAQIECBAgAABAgQIECBAoH0LiA7te39OT4AAAQIECBAgQIAAAQIEKisgOlR2NQ5GgAABAgQIECBAgAABAgTat4Do0L735/QECBAgQIAAAQIECBAgQKCyAqJDZVfjYAQIECBAgAABAgQIECBAoH0LiA7te39OT4AAAQIECBAgQIAAAQIEKisgOlR2NQ5GgAABAgQIECBAgAABAgTat4Do0L735/QECBAgQIAAAQIECBAgQKCyAqJDZVfjYAQIECBAgAABAgQIECBAoH0LiA7te39OT4AAAQIECBAgQIAAAQIEKisgOlR2NQ5GgAABAgQIECBAgAABAgTat4Do0L735/QECBAgQIAAAQIECBAgQKCyAqJDZVfjYAQIECBAgAABAgQIECBAoH0LiA7te39OT4AAAQIECBAgQIAAAQIEKisgOlR2NQ5GgAABAgQIECBAgAABAgTat4Do0L735/QECBAgQIAAAQIECBAgQKCyAqJDZVfjYAQIECBAgAABAgQIECBAoH0LiA7te39OT4AAAQIECBAgQIAAAQIEKisgOlR2NQ5GgAABAgQIECBAgAABAgTat4Do0L735/QECBAgQIAAAQIECBAgQKCyAqJDZVfjYAQIECBAgAABAgQIECBAoH0LiA7te39OT4AAAQIECBAgQIAAAQIEKisgOlR2NQ5GgAABAgQIECBAgAABAgTat4Do0L735/QECBAgQIAAAQIECBAgQKCyAqJDZVfjYAQIECBAgAABAgQIECBAoH0LiA7te39OT4AAAQIECBAgQIAAAQIEKisgOlR2NQ5GgAABAgQIECBAgAABAgTat8D/B9BhGCOaKKDsAAAAAElFTkSuQmCC",
      "text/html": [
       "<div>\n",
       "        \n",
       "        \n",
       "            <div id=\"9f64383e-adf5-4124-b282-2e321ff7f54e\" 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(\"9f64383e-adf5-4124-b282-2e321ff7f54e\")) {\n",
       "                    Plotly.newPlot(\n",
       "                        '9f64383e-adf5-4124-b282-2e321ff7f54e',\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, 1, null, 0, 0, null, 0, 1, null, 0, 0, null, 0, 0.2, null, 0.2, 0.2, null, 0.2, 1.0, null, 0.2, 0.2, null, 0.2, 0.4, null], \"y\": [null, null, 2.6, 3.9000000000000004, null, 3.9000000000000004, 3.9000000000000004, null, 2.6, 2.6, null, 2.6, 2.6, null, 2.6, 0.08125000000000004, null, 0.08125000000000004, 0.08125000000000004, null, 0.08125000000000004, 1.3, null, 1.3, 1.3, null, 0.08125000000000004, -1.1375, null, -1.1375, -1.1375, null]}, {\"line\": {\"color\": \"red\", \"width\": 1.25}, \"mode\": \"lines\", \"type\": \"scatter\", \"x\": [0.4, 0.4, null, 0.4, 1.0, null, 0.4, 0.4, null, 0.4, 0.6000000000000001, null, 0.6000000000000001, 0.6000000000000001, null, 0.6000000000000001, 1.0, null, 0.6000000000000001, 0.6000000000000001, null, 0.6000000000000001, 0.8, null, 0.8, 0.8, null, 0.8, 1.0, null, 0.8, 0.8, null, 0.8, 1.0, null], \"y\": [-1.1375, 0.0, null, 0.0, 0.0, null, -1.1375, -2.275, null, -2.275, -2.275, null, -2.275, -1.3, null, -1.3, -1.3, null, -2.275, -3.25, null, -3.25, -3.25, null, -3.25, -2.6, null, -2.6, -2.6, null, -3.25, -3.9000000000000004, null, -3.9000000000000004, -3.9000000000000004, null]}, {\"hoverinfo\": \"text\", \"marker\": {\"color\": \"black\", \"size\": 3, \"symbol\": \"circle\"}, \"mode\": \"markers\", \"showlegend\": false, \"text\": [\"root\", \"Galago\", \"HowlerMon\", \"edge.3\", \"Rhesus\", \"edge.2\", \"Orangutan\", \"edge.1\", \"Gorilla\", \"edge.0\", \"Chimpanzee\", \"Human\"], \"type\": \"scatter\", \"x\": [0, 1, 1, 0.2, 1.0, 0.4, 1.0, 0.6000000000000001, 1.0, 0.8, 1.0, 1.0], \"y\": [2.6, 3.9000000000000004, 2.6, 0.08125000000000004, 1.3, -1.1375, 0.0, -2.275, -1.3, -3.25, -2.6, -3.9000000000000004]}],\n",
       "                        {\"annotations\": [{\"font\": {\"family\": \"Inconsolata, monospace\", \"size\": 12}, \"showarrow\": false, \"text\": \"Galago\", \"textangle\": 0, \"x\": 1.025, \"xanchor\": \"left\", \"xref\": \"x\", \"y\": 3.9000000000000004, \"yref\": \"y\"}, {\"font\": {\"family\": \"Inconsolata, monospace\", \"size\": 12}, \"showarrow\": false, \"text\": \"HowlerMon\", \"textangle\": 0, \"x\": 1.025, \"xanchor\": \"left\", \"xref\": \"x\", \"y\": 2.6, \"yref\": \"y\"}, {\"font\": {\"family\": \"Inconsolata, monospace\", \"size\": 12}, \"showarrow\": false, \"text\": \"Rhesus\", \"textangle\": 0, \"x\": 1.025, \"xanchor\": \"left\", \"xref\": \"x\", \"y\": 1.3, \"yref\": \"y\"}, {\"font\": {\"family\": \"Inconsolata, monospace\", \"size\": 12}, \"showarrow\": false, \"text\": \"Orangutan\", \"textangle\": 0, \"x\": 1.025, \"xanchor\": \"left\", \"xref\": \"x\", \"y\": 0.0, \"yref\": \"y\"}, {\"font\": {\"family\": \"Inconsolata, monospace\", \"size\": 12}, \"showarrow\": false, \"text\": \"Gorilla\", \"textangle\": 0, \"x\": 1.025, \"xanchor\": \"left\", \"xref\": \"x\", \"y\": -1.3, \"yref\": \"y\"}, {\"font\": {\"family\": \"Inconsolata, monospace\", \"size\": 12}, \"showarrow\": false, \"text\": \"Chimpanzee\", \"textangle\": 0, \"x\": 1.025, \"xanchor\": \"left\", \"xref\": \"x\", \"y\": -2.6, \"yref\": \"y\"}, {\"font\": {\"family\": \"Inconsolata, monospace\", \"size\": 12}, \"showarrow\": false, \"text\": \"Human\", \"textangle\": 0, \"x\": 1.025, \"xanchor\": \"left\", \"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}, \"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('9f64383e-adf5-4124-b282-2e321ff7f54e');\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": [
    "from cogent3 import load_tree\n",
    "from cogent3.app import io, evo\n",
    "\n",
    "tree = load_tree(\"../data/primate_brca1.tree\")\n",
    "fig = tree.get_figure(contemporaneous=True)\n",
    "fig.style_edges(\"Human\", tip2=\"Orangutan\",\n",
    "                outgroup=\"Galago\",\n",
    "                line=dict(color=\"red\"))\n",
    "fig.show(width=500, height=500)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Specify a distinct rate matrix for the great apes\n",
    "\n",
    "As illustrated above, this is done by specifying the names of two \"ingroup\" tips (their last common ancestor is the ancestor of all edges of interest) and an \"outgroup\" tip. We can define the clade by Orangutan and one of Human, Chimpanzee or Gorilla. Any tip not a member of the clade outside suffices as an outgroup.\n",
    "\n",
    "To construct a `model` with this information, we create a dictionary that we assign to the `time_het` argument as follows."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "time_het = evo.model(\"GN\", tree=tree,\n",
    "                     time_het=[dict(tip_names=[\"Human\", \"Orangutan\"],\n",
    "                                     outgroup_name=\"Galago\")])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The outcome of this setting to `model` is that the black edges will share one rate matrix, the red edges another. We fit the model to the sample data set."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "reader = io.load_aligned(format=\"fasta\")\n",
    "aln = reader(\"../data/primate_brca1.fasta\")\n",
    "result = time_het(aln)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Examining the fitted likelihood function, the rate parameters differ between the edge sets."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<h4>GN</h4>\n",
       "<p>log-likelihood = -6982.2248</p>\n",
       "<p>number of free parameters = 36</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 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>parent</th>\n",
       "<th>length</th>\n",
       "<th>A&gt;C</th>\n",
       "<th>A&gt;G</th>\n",
       "<th>A&gt;T</th>\n",
       "<th>C&gt;A</th>\n",
       "<th>C&gt;G</th>\n",
       "<th>C&gt;T</th>\n",
       "</thead>\n",
       "<tbody>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">Galago</td>\n",
       "<td>root</td>\n",
       "<td style=\"font-family: monospace !important;\">0.1730</td>\n",
       "<td style=\"font-family: monospace !important;\">1.0064</td>\n",
       "<td style=\"font-family: monospace !important;\">3.7857</td>\n",
       "<td style=\"font-family: monospace !important;\">1.0385</td>\n",
       "<td style=\"font-family: monospace !important;\">1.7266</td>\n",
       "<td style=\"font-family: monospace !important;\">2.3267</td>\n",
       "<td style=\"font-family: monospace !important;\">6.3442</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">HowlerMon</td>\n",
       "<td>root</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0451</td>\n",
       "<td style=\"font-family: monospace !important;\">1.0064</td>\n",
       "<td style=\"font-family: monospace !important;\">3.7857</td>\n",
       "<td style=\"font-family: monospace !important;\">1.0385</td>\n",
       "<td style=\"font-family: monospace !important;\">1.7266</td>\n",
       "<td style=\"font-family: monospace !important;\">2.3267</td>\n",
       "<td style=\"font-family: monospace !important;\">6.3442</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">Rhesus</td>\n",
       "<td>edge.3</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0215</td>\n",
       "<td style=\"font-family: monospace !important;\">1.0064</td>\n",
       "<td style=\"font-family: monospace !important;\">3.7857</td>\n",
       "<td style=\"font-family: monospace !important;\">1.0385</td>\n",
       "<td style=\"font-family: monospace !important;\">1.7266</td>\n",
       "<td style=\"font-family: monospace !important;\">2.3267</td>\n",
       "<td style=\"font-family: monospace !important;\">6.3442</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">Orangutan</td>\n",
       "<td>edge.2</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0077</td>\n",
       "<td style=\"font-family: monospace !important;\">0.2103</td>\n",
       "<td style=\"font-family: monospace !important;\">3.5008</td>\n",
       "<td style=\"font-family: monospace !important;\">0.2093</td>\n",
       "<td style=\"font-family: monospace !important;\">0.9013</td>\n",
       "<td style=\"font-family: monospace !important;\">1.3534</td>\n",
       "<td style=\"font-family: monospace !important;\">4.9566</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">Gorilla</td>\n",
       "<td>edge.1</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0025</td>\n",
       "<td style=\"font-family: monospace !important;\">0.2103</td>\n",
       "<td style=\"font-family: monospace !important;\">3.5008</td>\n",
       "<td style=\"font-family: monospace !important;\">0.2093</td>\n",
       "<td style=\"font-family: monospace !important;\">0.9013</td>\n",
       "<td style=\"font-family: monospace !important;\">1.3534</td>\n",
       "<td style=\"font-family: monospace !important;\">4.9566</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">Human</td>\n",
       "<td>edge.0</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0061</td>\n",
       "<td style=\"font-family: monospace !important;\">0.2103</td>\n",
       "<td style=\"font-family: monospace !important;\">3.5008</td>\n",
       "<td style=\"font-family: monospace !important;\">0.2093</td>\n",
       "<td style=\"font-family: monospace !important;\">0.9013</td>\n",
       "<td style=\"font-family: monospace !important;\">1.3534</td>\n",
       "<td style=\"font-family: monospace !important;\">4.9566</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">Chimpanzee</td>\n",
       "<td>edge.0</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0028</td>\n",
       "<td style=\"font-family: monospace !important;\">0.2103</td>\n",
       "<td style=\"font-family: monospace !important;\">3.5008</td>\n",
       "<td style=\"font-family: monospace !important;\">0.2093</td>\n",
       "<td style=\"font-family: monospace !important;\">0.9013</td>\n",
       "<td style=\"font-family: monospace !important;\">1.3534</td>\n",
       "<td style=\"font-family: monospace !important;\">4.9566</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">edge.0</td>\n",
       "<td>edge.1</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0000</td>\n",
       "<td style=\"font-family: monospace !important;\">0.2103</td>\n",
       "<td style=\"font-family: monospace !important;\">3.5008</td>\n",
       "<td style=\"font-family: monospace !important;\">0.2093</td>\n",
       "<td style=\"font-family: monospace !important;\">0.9013</td>\n",
       "<td style=\"font-family: monospace !important;\">1.3534</td>\n",
       "<td style=\"font-family: monospace !important;\">4.9566</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">edge.1</td>\n",
       "<td>edge.2</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0034</td>\n",
       "<td style=\"font-family: monospace !important;\">0.2103</td>\n",
       "<td style=\"font-family: monospace !important;\">3.5008</td>\n",
       "<td style=\"font-family: monospace !important;\">0.2093</td>\n",
       "<td style=\"font-family: monospace !important;\">0.9013</td>\n",
       "<td style=\"font-family: monospace !important;\">1.3534</td>\n",
       "<td style=\"font-family: monospace !important;\">4.9566</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">edge.2</td>\n",
       "<td>edge.3</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0121</td>\n",
       "<td style=\"font-family: monospace !important;\">1.0064</td>\n",
       "<td style=\"font-family: monospace !important;\">3.7857</td>\n",
       "<td style=\"font-family: monospace !important;\">1.0385</td>\n",
       "<td style=\"font-family: monospace !important;\">1.7266</td>\n",
       "<td style=\"font-family: monospace !important;\">2.3267</td>\n",
       "<td style=\"font-family: monospace !important;\">6.3442</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">edge.3</td>\n",
       "<td>root</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0076</td>\n",
       "<td style=\"font-family: monospace !important;\">1.0064</td>\n",
       "<td style=\"font-family: monospace !important;\">3.7857</td>\n",
       "<td style=\"font-family: monospace !important;\">1.0385</td>\n",
       "<td style=\"font-family: monospace !important;\">1.7266</td>\n",
       "<td style=\"font-family: monospace !important;\">2.3267</td>\n",
       "<td style=\"font-family: monospace !important;\">6.3442</td>\n",
       "</tr>\n",
       "</tbody>\n",
       "</table>\n",
       "<table>\n",
       "<thead style=\"background: rgba(161, 195, 209, 0.75); font-weight: bold; text-align: center;\">\n",
       "<th>edge</th>\n",
       "<th>G&gt;A</th>\n",
       "<th>G&gt;C</th>\n",
       "<th>G&gt;T</th>\n",
       "<th>T&gt;A</th>\n",
       "<th>T&gt;C</th>\n",
       "</thead>\n",
       "<tbody>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">Galago</td>\n",
       "<td style=\"font-family: monospace !important;\">9.0145</td>\n",
       "<td style=\"font-family: monospace !important;\">1.2462</td>\n",
       "<td style=\"font-family: monospace !important;\">0.7647</td>\n",
       "<td style=\"font-family: monospace !important;\">1.4577</td>\n",
       "<td style=\"font-family: monospace !important;\">3.5709</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">HowlerMon</td>\n",
       "<td style=\"font-family: monospace !important;\">9.0145</td>\n",
       "<td style=\"font-family: monospace !important;\">1.2462</td>\n",
       "<td style=\"font-family: monospace !important;\">0.7647</td>\n",
       "<td style=\"font-family: monospace !important;\">1.4577</td>\n",
       "<td style=\"font-family: monospace !important;\">3.5709</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">Rhesus</td>\n",
       "<td style=\"font-family: monospace !important;\">9.0145</td>\n",
       "<td style=\"font-family: monospace !important;\">1.2462</td>\n",
       "<td style=\"font-family: monospace !important;\">0.7647</td>\n",
       "<td style=\"font-family: monospace !important;\">1.4577</td>\n",
       "<td style=\"font-family: monospace !important;\">3.5709</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">Orangutan</td>\n",
       "<td style=\"font-family: monospace !important;\">4.6227</td>\n",
       "<td style=\"font-family: monospace !important;\">1.1149</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0000</td>\n",
       "<td style=\"font-family: monospace !important;\">0.3394</td>\n",
       "<td style=\"font-family: monospace !important;\">2.9912</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">Gorilla</td>\n",
       "<td style=\"font-family: monospace !important;\">4.6227</td>\n",
       "<td style=\"font-family: monospace !important;\">1.1149</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0000</td>\n",
       "<td style=\"font-family: monospace !important;\">0.3394</td>\n",
       "<td style=\"font-family: monospace !important;\">2.9912</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">Human</td>\n",
       "<td style=\"font-family: monospace !important;\">4.6227</td>\n",
       "<td style=\"font-family: monospace !important;\">1.1149</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0000</td>\n",
       "<td style=\"font-family: monospace !important;\">0.3394</td>\n",
       "<td style=\"font-family: monospace !important;\">2.9912</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">Chimpanzee</td>\n",
       "<td style=\"font-family: monospace !important;\">4.6227</td>\n",
       "<td style=\"font-family: monospace !important;\">1.1149</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0000</td>\n",
       "<td style=\"font-family: monospace !important;\">0.3394</td>\n",
       "<td style=\"font-family: monospace !important;\">2.9912</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">edge.0</td>\n",
       "<td style=\"font-family: monospace !important;\">4.6227</td>\n",
       "<td style=\"font-family: monospace !important;\">1.1149</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0000</td>\n",
       "<td style=\"font-family: monospace !important;\">0.3394</td>\n",
       "<td style=\"font-family: monospace !important;\">2.9912</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">edge.1</td>\n",
       "<td style=\"font-family: monospace !important;\">4.6227</td>\n",
       "<td style=\"font-family: monospace !important;\">1.1149</td>\n",
       "<td style=\"font-family: monospace !important;\">0.0000</td>\n",
       "<td style=\"font-family: monospace !important;\">0.3394</td>\n",
       "<td style=\"font-family: monospace !important;\">2.9912</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">edge.2</td>\n",
       "<td style=\"font-family: monospace !important;\">9.0145</td>\n",
       "<td style=\"font-family: monospace !important;\">1.2462</td>\n",
       "<td style=\"font-family: monospace !important;\">0.7647</td>\n",
       "<td style=\"font-family: monospace !important;\">1.4577</td>\n",
       "<td style=\"font-family: monospace !important;\">3.5709</td>\n",
       "</tr>\n",
       "<tr>\n",
       "<td style=\"background: rgba(161, 195, 209, 0.25); font-weight: 600;\">edge.3</td>\n",
       "<td style=\"font-family: monospace !important;\">9.0145</td>\n",
       "<td style=\"font-family: monospace !important;\">1.2462</td>\n",
       "<td style=\"font-family: monospace !important;\">0.7647</td>\n",
       "<td style=\"font-family: monospace !important;\">1.4577</td>\n",
       "<td style=\"font-family: monospace !important;\">3.5709</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.3754</td>\n",
       "<td style=\"font-family: monospace !important;\">0.1768</td>\n",
       "<td style=\"font-family: monospace !important;\">0.2080</td>\n",
       "<td style=\"font-family: monospace !important;\">0.2398</td>\n",
       "</tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/plain": [
       "<cogent3.evolve.parameter_controller.AlignmentLikelihoodFunction at 0x11762a5f8>"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "result.lf"
   ]
  }
 ],
 "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
}
