{ "cells": [ { "cell_type": "markdown", "id": "5d0da7f5", "metadata": {}, "source": [ "# Multiverse integration\n", "\n", "Comet includes multiverse integration methods based on the weighting schemes introduced in:\n", "\n", "> Cantone, G. G., & Tomaselli, V. (2025). Characterisation and calibration of multiversal methods. \n", "> Advances in Data Analysis and Classification, 19(4), 989-1021. https://doi.org/10.1007/s11634-024-00610-9\n", "\n", "For illustration and validation purposes, we use the \"unidentified multiverse\" example from the paper, which was introduced in a previous paper by Del Giudice and Gangestad (2021). https://doi.org/10.6084/m9.figshare.12089736" ] }, { "cell_type": "code", "execution_count": 1, "id": "9999bfd8", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Removed 16 out of 256 universes:\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
UniverseDecision 1Value 1Decision 2Value 2Decision 3Value 3Decision 4Value 4Decision 5Value 5Decision 6Value 6Decision 7Value 7Decision 8Value 8
0Universe_1X10X20X30X41ANTE0MEDI0COLL0CONF0
1Universe_2X10X20X30X41ANTE0MEDI0COLL0CONF1
2Universe_3X10X20X30X41ANTE0MEDI0COLL1CONF0
3Universe_4X10X20X30X41ANTE0MEDI0COLL1CONF1
4Universe_5X10X20X30X41ANTE0MEDI1COLL0CONF0
......................................................
235Universe_236X11X21X31X41ANTE1MEDI0COLL1CONF1
236Universe_237X11X21X31X41ANTE1MEDI1COLL0CONF0
237Universe_238X11X21X31X41ANTE1MEDI1COLL0CONF1
238Universe_239X11X21X31X41ANTE1MEDI1COLL1CONF0
239Universe_240X11X21X31X41ANTE1MEDI1COLL1CONF1
\n", "

240 rows × 17 columns

\n", "
" ], "text/plain": [ " Universe Decision 1 Value 1 Decision 2 Value 2 Decision 3 Value 3 \\\n", "0 Universe_1 X1 0 X2 0 X3 0 \n", "1 Universe_2 X1 0 X2 0 X3 0 \n", "2 Universe_3 X1 0 X2 0 X3 0 \n", "3 Universe_4 X1 0 X2 0 X3 0 \n", "4 Universe_5 X1 0 X2 0 X3 0 \n", ".. ... ... ... ... ... ... ... \n", "235 Universe_236 X1 1 X2 1 X3 1 \n", "236 Universe_237 X1 1 X2 1 X3 1 \n", "237 Universe_238 X1 1 X2 1 X3 1 \n", "238 Universe_239 X1 1 X2 1 X3 1 \n", "239 Universe_240 X1 1 X2 1 X3 1 \n", "\n", " Decision 4 Value 4 Decision 5 Value 5 Decision 6 Value 6 Decision 7 \\\n", "0 X4 1 ANTE 0 MEDI 0 COLL \n", "1 X4 1 ANTE 0 MEDI 0 COLL \n", "2 X4 1 ANTE 0 MEDI 0 COLL \n", "3 X4 1 ANTE 0 MEDI 0 COLL \n", "4 X4 1 ANTE 0 MEDI 1 COLL \n", ".. ... ... ... ... ... ... ... \n", "235 X4 1 ANTE 1 MEDI 0 COLL \n", "236 X4 1 ANTE 1 MEDI 1 COLL \n", "237 X4 1 ANTE 1 MEDI 1 COLL \n", "238 X4 1 ANTE 1 MEDI 1 COLL \n", "239 X4 1 ANTE 1 MEDI 1 COLL \n", "\n", " Value 7 Decision 8 Value 8 \n", "0 0 CONF 0 \n", "1 0 CONF 1 \n", "2 1 CONF 0 \n", "3 1 CONF 1 \n", "4 0 CONF 0 \n", ".. ... ... ... \n", "235 1 CONF 1 \n", "236 0 CONF 0 \n", "237 0 CONF 1 \n", "238 1 CONF 0 \n", "239 1 CONF 1 \n", "\n", "[240 rows x 17 columns]" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "Starting multiverse analysis for all universes...\n" ] }, { "data": { "application/vnd.jupyter.widget-view+json": { "model_id": "5f0c6f7531ef48968761e0b826a304cb", "version_major": 2, "version_minor": 0 }, "text/plain": [ "Performing multiverse analysis:: 0%| | 0/240 [00:00" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "mverse.multiverse_plot(measure=\"beta\", n_bins=20, sig_col=\"p_value\", baseline=0, figsize=(7,9))" ] }, { "cell_type": "markdown", "id": "7b238452", "metadata": {}, "source": [ "Comet currently contains two methods for integrating over the multiverse:\n", "\n", "`uniform`\n", "- All universes are weighted equally (the estimate is just the plain mean/median across universes). This corresponds to the density function as plotted in the multiverse plot above.\n", "\n", "`mli` (Maximum Local Influence)\n", "- Each universe is compared with its nearest neighbours (specifications differing by a single decision). Universes whose estimate barely changes when one decision is flipped are weighted more, down-weighting results that hinge on a single arbitrary choice.\n", "\n", "As an example, we can use the `.integrate()` method to simply get an integrated measure over all universes. For example, using uniform weighting and mean aggregation:" ] }, { "cell_type": "code", "execution_count": 3, "id": "177e98c0", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Mean beta: -0.0066\n" ] } ], "source": [ "estimate, weights = mverse.integrate(measure=\"beta\", method=\"uniform\", type=\"mean\")\n", "print(f\"Mean beta: {estimate:.4f}\")" ] }, { "cell_type": "markdown", "id": "33dbe1c3", "metadata": {}, "source": [ "We can also use the `compare_integration` and `plot_integration` convenience functions to get a summary for both integration methods." ] }, { "cell_type": "code", "execution_count": 4, "id": "d3c42976", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
Schememedianmeangini_werr_medianerr_mean
0Uniform-0.0210-0.00660.0000.2210.207
1MLI-0.0328-0.01440.2430.2330.214
\n", "
" ], "text/plain": [ " Scheme median mean gini_w err_median err_mean\n", "0 Uniform -0.0210 -0.0066 0.000 0.221 0.207\n", "1 MLI -0.0328 -0.0144 0.243 0.233 0.214" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "table, weights = mverse.compare_integration(measure=\"beta\", true_value=0.2)\n", "table" ] }, { "cell_type": "markdown", "id": "91f5e652", "metadata": {}, "source": [ "The columns are, for each weighting scheme:\n", "\n", "- `median` / `mean`: the integrated estimate (weighted median/mean of `beta`)\n", "- `gini_w`: how concentrated the weights are (0 = all equal, → 1 = dominated by a few universes).\n", "\n", "If `true_value` is passed, additional `err_*` columns give the absolute error of each estimate against it. Usually, this is not available in real multiverses, but we can include it here as we have the ground truth from the simulation." ] }, { "cell_type": "code", "execution_count": 5, "id": "daad2d7a", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAArIAAAEiCAYAAAAF9zFeAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjcsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvTLEjVAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAfMRJREFUeJzt3XdYFNf79/H3LlW6KIoo2LB3UbFXrMTee+9GE02ixtgTNdGoMZbYS+zd2LvYRVHsvYEKdpAiLOzO8wdP9he+GqO4y7Jwv65rr+DM7Dn3jAQ/zJ45R6UoioIQQgghhBBmRm3qAoQQQgghhEgJCbJCCCGEEMIsSZAVQgghhBBmSYKsEEIIIYQwSxJkhRBCCCGEWZIgK4QQQgghzJIEWSGEEEIIYZYkyAohhBBCCLNkaeoCPodOp+PJkyc4OjqiUqlMXY4QQgghhPhMiqIQFRWFh4cHavWH77madZB98uQJnp6epi5DCCGEEEIYWGhoKLly5frgMWYdZB0dHYGkE3VycjJxNUIIIYQQqUej0fDrr78CMGzYMKytrU1ckWG8efMGT09Pfc77ELMOsn8PJ3BycpIgK4QQQogMRaPRYGtrCyRlofQSZP/2McNG5WEvIYQQQghhliTICiGEEEIIsyRBVgghhBBCmCWzHiP7sbRaLQkJCaYuQ/wPa2vr/5xWQwghhBDvZ2VlRf/+/fVfZ0TpOsgqikJ4eDgRERGmLkW8h1qtJm/evOlucLoQQgiRGlQqFdmyZTN1GSaVroPs3yE2W7Zs2NnZyaIJacjfi1mEhYXh5eUlfzdCCCGE+GTpNshqtVp9iM2SJYupyxHv4ebmxpMnT0hMTMywH4n8G0VRCIsO42HEQxQUrC2ssbawxtbSljwuebC2kLvYQgiR0Wm1Wo4dOwZAtWrVsLCwMHFFqS/dBtm/x8Ta2dmZuBLxb/4eUqDVajN8kH0W84ydt3ay5+4erj+/zt3Xd4lNiH3vsVZqK4q4FaG0e2nKuJfhi4Jf4O3qncoVCyGEMDWtVktAQAAAlStXliCb2saNG8f48eOTbStUqBA3btwwWB/ykXXaldH/bl69fcXi84vZcmMLpx+dRkFJtl+tUuPp5Iml2hKNVoNGqyFaE01MQgyXnl7i0tNLrLi4gq/3fo1PDh/aFW9Hm2Jt8HL2MtEZCSGEEKnL5HdkixUrxoEDB/R/trQ0eUlCGFVoZCjTT01n4fmFxCTE6LeXzVGWxgUb45vTF29Xb3K75H5nCIGiKDyMfMjF8ItcfHqRYyHHOHz/MEFhQQSFBfHd/u9oVbQVP1T/gZLZS6b2qQkhhBCpyuSp0dLSEnd3d1OXkS7lyZOHr776iq+++gpIevitc+fOnDx5EisrK5nNIZWFRYXx/aHvWXlpJYm6RABKZS9FX5++NC7UmFxOuf6zDZVKRR6XPORxyUPTwk0BeB7znE3XN7H2yloCHgaw4doGNlzbQPPCzRldfTRlcpQx6nkJIYQQpmLySTxv376Nh4cH+fLlo2PHjoSEhJi6JJOrWbOmPnz+07Jly3Bxcfnods6ePUufPn30f54xYwZhYWEEBwdz69YtA1QqPoZWp2Xu2bkUnlOYZcHLSNQlUitPLfZ03MOFvhfoX77/R4XYf+Nm70a/cv040u0Il/pdok2xNqhQseXGFsouKEvf7X2JiIsw3AkJIYQQaYRJg6yvry/Lli1jz549zJs3j/v371OtWjWioqLee3x8fDxv3rxJ9hL/zs3NLdnDbnfv3sXHx4cCBQqkeN45jUZjqPIyhODwYCovqczAXQN5E/+Gch7lONXzFIe6HqK+d32DjxMukb0E61qt48qAK7Qv3h6ABecXUGROETZc3YCiKP/RghBCCGE+TBpkGzZsSOvWrSlZsiT169dn165dREREsH79+vceP3nyZJydnfUvT0/PVK447ejWrRvNmjVj2rRp5MiRgyxZsjBw4MBkK5jlyZOHmTNn6r/etGkTK1asQKVS0a1bNwBCQkJo2rQpDg4OODk50aZNG54+fapvY9y4cZQuXZpFixaRN29ebG1tgaSPuOfPn88XX3yBnZ0dRYoU4dSpU9y5c4eaNWtib29P5cqVuXv3bqpdk7REURRmnZlF+YXlCXwciKO1I783/J3TPU9TMVdFo/df1K0oq1uuJqBbAIWyFCI8Opw2G9vQZG0TnsU8M3r/QgghRGow+RjZf3JxcaFgwYLcuXPnvftHjhzJ0KFD9X9+8+bNJ4VZRVGIjX3/lEbGZKzFGA4fPkyOHDk4fPgwd+7coW3btpQuXZrevXu/c+zZs2fp0qULTk5O/Pbbb2TKlAmdTqcPsQEBASQmJjJw4EDatm3LkSNH9O+9c+cOmzZtYvPmzcmm9pg4cSLTp09n+vTpDB8+nA4dOpAvXz5GjhyJl5cXPXr0YNCgQezevdvg556WRcVH0Wt7L9ZfTfqFrHnh5sxuNBsPR49Ur6V67upc7HeRyccnM+nYJHbc2kHpP0qzrtU6quWulur1CCGEMBxLS0t69eql/zojSlNnHR0dzd27d+ncufN799vY2GBjY5Pi9mNjY3FwcEjx+1MqOjoae3t7g7ebOXNmZs+ejYWFBYULF8bf35+DBw++N8i6ublhY2NDpkyZ9A/X7d+/n8uXL3P//n39LwQrVqygWLFinD17lvLlywNJwwlWrFiBm5tbsja7d+9OmzZtABg+fDiVKlVi9OjR1K9fH4AhQ4bQvXt3g593Wnb12VVarm/JzZc3sVRbMq3uNAb7DjbpVGM2ljaMqzmO1kVb02ZjG649v0at5bWYVGcS31T+BrXK5EPlhRBCpIBarSZnzpymLsOkTPov2DfffENAQAAPHjzg5MmTNG/eHAsLC9q3b2/KssxGsWLFkt0hzZEjB8+effzHxtevX8fT0zPZXe2iRYvi4uLC9evX9dty5879TogFKFny/6Z3yp49OwAlSpRIti0uLi7DjGU+cO8AFRdX5ObLm+R0zMnRbkcZUnFImpkvt1i2YgT2CqRTyU5oFS3DDwyn2dpmvInPGH8/Qggh0h+T3pF99OgR7du35+XLl7i5uVG1alVOnz793tBkCHZ2dkRHRxul7f/q91M4OTkRGRn5zvaIiAicnZ31f/7f1bBUKhU6nS5lRX7Av91N/mf/f4e1920zRk1pzfqr6+m0uRMJugRq5qnJ+lbrcbM3zvfx57C3tmdFsxVU86rG4N2D2X5rO9WWVmNXh13kdMrYv9ULIYS50Wq1nD59GoCKFSvKyl6pbe3atanan0qlMspH/IZWqFAh9u3b98728+fPU7BgQYP1U6RIEUJDQwkNDdXflb127RoREREULVrUYP2kd3MC5/Dl7i9RUGhdtDV/Nv8TG8uUD4ExNpVKRR+fPpTNUZYvVn/BpaeXqLi4Irs67KJE9hL/3YAQQog0QavV6heVKl++fIYMsjI4Lg3q378/t27dYvDgwVy6dImbN28yffp01qxZw7BhwwzWj5+fHyVKlKBjx46cP3+ewMBAunTpQo0aNShXrpzB+kmvFEVh3JFxDNo9CAWFAeUGsKblmjQdYv+pnEc5Tvc6TeGshXn05hFVl1blwL0D//1GIYQQIo2QIJsG5cuXj6NHj3Ljxg38/Pzw9fVl/fr1bNiwgQYNGhisH5VKxbZt28icOTPVq1fHz8+PfPnysW7dOoP1kZ79dOwnxgeMB2B8zfHMbjQbC7V5/TacxyUPJ3ucpHru6ryJf0PDVQ3ZfH2zqcsSQgghPopKMeMZ0t+8eYOzszORkZE4OTkl2xcXF8f9+/eTzX0q0hZz/juafmo6w/Yl3R2fVncawyob7k65KcQnxtNlaxfWX12PhcqClS1W0q54O1OXJYQQ4gM0Gg2TJ08GkqYotba2NnFFhvGhfPe/5I6sEJ9o3tl5+hA7sdZEsw+xkDRF1+oWq+laqitaRUvHzR1ZHrzc1GUJIYQQHyRBVohPsOLiCgbsGgDAyKojGVVtlIkrMhwLtQVLmi6hd9ne6BQd3bZ1Y0HQAlOXJYQQQvwrCbJCfKT9d/fTY1sPAIb4DuGn2j+lmTliDUWtUjP/i/l8WeFLAPru6MvSC0tNXJUQQgjxfmlqZS8h0qorz67QakMrtIqWTiU7MaP+jHQXYv+mUqn4rcFvWKgsmHlmJr2298LOyo62xduaujQhhBD/YGlpSdeuXfVfZ0QZ86yF+ATh0eH4r/bnTfwbqueuzqLGi9JtiP2bSqViev3pxCbEsuD8Ajpt6YSdlR2NCzU2dWlCCCH+P7VaTZ48eUxdhknJ0AIhPiA2IZbGaxoTEhlCwSwF2dJ2i9nME/u5VCoVc/3n0rFERxJ1ibTa0ErmmRVCCJGmSJAV4l8oikKXLV049+QcWTJlYWeHnbhmcjV1WanKQm3BsmbLaF64ORqthqZrm3LuyTlTlyWEEIKklb0CAwMJDAxEq9WauhyTkCArxL/45cQvbLq+CSu1FVvbbcXb1dvUJZmEpdqSNS3XUDdfXWITYvFf7c+91/dMXZYQQmR4Wq2W3bt3s3v3bgmyQoj/c+DeAb4/9D0Avzf8napeVU1ckWnZWNqwqc0myriX4VnMMxqsbMDzmOemLksIIUQGJ0E2DerWrRsqlYp+/fq9s2/gwIGoVCq6deumP7ZZs2b/2laePHmYOXOmcQpNpx5GPKTdxnboFB3dS3enj08fU5eUJjjaOLKzw05yO+fm9qvbNF7TmNiEWFOXJYQQIgOTIJtGeXp6snbtWt6+favfFhcXx+rVq/Hy8jJhZelbXGIcLde35OXbl5TNUZY5jeak+xkKPkUOxxzs6bQH10yunHl8hnYb25GoSzR1WUIIITIoCbJpVNmyZfH09GTz5s36bZs3b8bLy4syZcqYsLL0bcjuIQSFBeGayZVNbTaRySqTqUtKcwpnLcz29tuxtbRl+63tDNw5EEVRTF2WEEKIDEjmkU3DevTowdKlS+nYsSMAS5YsoXv37hw5csS0haVTG65uYMH5BahQsablGvK45DFaX1qtlhs3bhAUFMSDBw948uQJT548ISwsjISEBNRqtf6VLVs2PD098fT0xMvLi7Jly1K4cGHUatP9HlrZszJrWq6h5fqWLDi/AE9nT36o/oPJ6hFCCJExZaggqyiKScb02VnZpejj6U6dOjFy5EgePnwIwIkTJ1i7dq0EWSN4EPGA3tt7AzCi6gjq5a9n0PZ1Oh2BgYH89ddfnDx5kqCgIKKjo1PcnrOzM76+vlSqVAl/f3/KlSuX6kMgmhVuxu8Nf2fgroGMPjyaXE656Fa6W6rWIIQQImPLUEE2NiEWh8kOqd5v9Mho7K3tP/l9bm5u+Pv7s2zZMhRFwd/fn6xZsxqhwowtUZdIh00diIyPpGKuioyvOd4g7ep0Og4dOsTGjRvZtm0b4eHhyfbb29vj4+NDoUKF8PDwwMPDA3d3d2xsbFAUBUVRSEhIIDw8nNDQUEJDQ7l79y7nz58nMjKSffv2sW/fPsaPH4+XlxctW7akVatWVKpUKdVC7YDyAwiNDGXKiSn0+qsX7g7uNPBukCp9CyFERmdpaUn79u31X2dEGfOszUiPHj0YNGgQAHPmzDFxNenT+CPjOfXoFE42TqxusRorC6vPau/Vq1csW7aMefPmcefOHf12Jycn/P39qVu3LhUqVKBw4cJYWFh8cvuJiYlcvnyZ06dPc+jQIXbv3k1ISAgzZsxgxowZFC9enEGDBtGpUyfs7T/9F6hPNanOJB5HPebPS3/SekNrjnU/Rmn30kbvVwghMjq1Wk3BggVNXYZJZagga2dlR/TIlH+c+zn9plSDBg3QaDSoVCrq169vwKoEwJEHR/jp2E8ALPhiAXkz501xW/fv32fSpEmsXLmSuLg4ICm8tm/fnhYtWlCzZk2sra0/u2ZLS0vKlClDmTJl6N+/P2/fvmXv3r1s2rSJLVu2cOXKFfr168eIESPo2bMnw4YNI0eOHJ/d779RqVQsarKIJ1FPOHj/IP6r/TnT6wy5nHIZrU8hhBACMliQValUKfqI35QsLCy4fv26/uv3iYyMJDg4ONm2LFmy4OnpaezyzFpkXCRdtnRBQaFnmZ60Ld42Re2Ehoby008/sXjxYhITk6aiKlWqFAMHDqRDhw5GvyuaKVMmmjVrRrNmzYiIiGDZsmXMnj2bu3fv8uuvvzJ37lwGDx7Md999h6urcZbYtbawZmObjVRZUoVrz6/hv9qfY92P4WTjZJT+hBBCJD04fPnyZQBKlCiRok/5zJ1Mv2UGnJyccHL690Bw5MgR/R26v1/jxxtmnGd6NmTPEELfhJI/c35mNpj5ye+PjIxk6NCheHt7M3/+fBITE6lbty7Hjh3jwoUL9O7dO1U+2v8nFxcXvvrqK27dusXOnTupVKkSb9++5eeffyZv3rxMnDgx2dzEBu3b1oVdHXaR3T47l55eou3GtjLHrBBCGJFWq2Xbtm1s27Ytwy5Rq1LMeALIN2/e4OzsTGRk5DtBLy4ujvv375M3b15sbW1NVKH4EFP+HW29sZXm65qjQsWx7seo4lXlo9+rKAqrVq3im2++4enTpwDUqFGDCRMmUL16dWOVnCKKorBr1y5GjRrFxYsXAcibNy+///47/v7+Runz3JNz1FhWg9iEWPqU7cMfX/whi0oIIYQRaDQaJk+eDMDIkSMNMnwtLfhQvvtfckdWZDjPYp7RZ3vSsrPfVfnuk0LstWvXqFmzJp07d+bp06cULFiQPXv2cOTIkTQXYiFpOI2/vz/nz59nzZo15MqVi/v37/PFF1/QokULQkJCDN5nOY9yrG6xGhUqFpxfwNSTUw3ehxBCCAESZEUGoygK/Xb043nsc0pkK/HRU21ptVqmTZtGmTJlOHr0KJkyZWLSpElcunTJLB7CU6vVtGvXjuvXr/Ptt99iaWnJli1bKFq0KPPnzzf4ylxNCzfVD9cYfmA4G65uMGj7QgghBEiQFRnMyksr2XJjC1ZqK1Y0X4GNpc1/vufBgwfUrl2bb7/9Fo1Gg7+/P9evX2fkyJHY2Pz3+9MSBwcHfvnlFy5cuEDVqlWJiYmhX79+fPHFF4SFhRm0r8G+gxlcYTAAnbd05mToSYO2L4QQQkiQFRlGeHQ4Q/YMAWBsjbEfNdfpqlWrKFmyJEePHsXe3p6FCxeyfft2cufObeRqjat48eIEBAQwffp0bGxs2LVrFyVKlGDTpk0G7Wd6/ek0KdSEeG08Tdc25e6ruwZtXwghRMYmQVZkGIN2DeJ13GvKuJfhuyrfffDY+Ph4+vfvT6dOnYiKiqJKlSpcvHiRXr16pZsHl9RqNV9//TVBQUGULl2aly9f0qpVKwYNGkR8fLxB+rBQW7C6xWp8cvjwIvYFjVY34mXsS4O0LYQQQqT7IGvGkzKke6n5d7Pp2iY2Xd+EpdqSJU2XfHD1rpCQEKpXr84ffyQ9bT927FgCAgLInz9/qtWbmooVK8aZM2cYMWIEkLSCXLVq1Xjw4IFB2re3tmdHhx14OXtx6+Utmq9rTnyiYYKyEEJkZJaWlrRq1YpWrVpl2CVq022QtbJKCiqxsbEmrkT8G41GA/z7Qg+G8urtKwbuGgjA8CrDPzik4ODBg5QtW5bAwEBcXV3ZtWsX48aNS/eTTFtbWzN58mR27tyJq6srZ8+epWzZsuzcudMg7bs7uLOrwy6cbJw4FnKM7tu6o1N0BmlbCCEyKrVaTbFixShWrBhqdbqNdB+UbueRBQgLCyMiIoJs2bJhZ2eXbj4STg90Oh1PnjzBysoKLy8vo/7ddNvajeUXl1MkaxEu9L3wrw94LVy4kP79+6PVavHx8WHjxo3kyZPHaHWlVQ8fPqRNmzYEBgaiUqn48ccfGTlypEH+jg7eO0iDVQ1I1CUyqtoofqz9owEqFkIIkZ58yjyy6TrIKopCeHg4ERERqV+c+E9qtZq8efMadQLnvXf20mBVA1SoONHjBJU8K71zjE6nY8SIEUydmjTfaadOnVi4cGGGXkhDo9Hw9ddfM3fuXADatWvH4sWLsbOz++y2l15YSo+/egCwsPFCepXt9dltCiFERqTT6fTL2BcpUiTd3JX9lCCbZgZUTJkyhZEjRzJkyBBmzpxpkDZVKhU5cuQgW7ZsJCQkGKRNYTjW1tZG/Z/ubcJbBuwaACRNBfW+EBsbG0unTp3YsmULAOPHj2f06NEZ/u69tbU1c+bMoWTJkgwaNIi1a9dy+/Zttm7dSq5cuT6r7e5lunPv9T1+PPYjfXf0Jbt9dhoXamygyoUQIuNITExk48aNQPpa2etTpIkge/bsWebPn0/JkiWN0r6FhUW6H+Mo3jXp2CTuvb5HTsecTKw18Z39r1+/xt/fn1OnTmFtbc3SpUvp0KGDCSpNu/r27UvhwoVp2bIlQUFBVKhQgV27dlG6dOnPandCrQk8jnrM0uCltNnYhoNdDlLZs7JhihZCCJFhmPwedHR0NB07dmThwoVkzpzZ1OWIdOLGixv8fOJnAGY1nIWjjWOy/WFhYdSoUYNTp06ROXNmDh48KCH2X9SoUYOzZ89SrFgxwsLCqF69OgcPHvysNlUqFQsaL8C/gD9xiXF8sfoLrj2/ZqCKhRBCZBQmD7IDBw7E398fPz+//zw2Pj6eN2/eJHsJ8b8URaH/zv4k6BLwL+BP88LNk+2/d+8eVatW5fLly+TIkYOjR49StWpVE1VrHvLmzcvx48epUaMGUVFRNGzYkFWrVn1Wm5ZqS9a3Xk/FXBV5Hfea+ivrExoZaqCKhRBCZAQmDbJr167l/PnzTJ48+aOOnzx5Ms7OzvqXp6enkSsU5mjV5VUceXCETJaZ+L3h78nGu165coWqVaty79498uXLx/HjxylevLgJqzUfLi4u7N27lzZt2pCQkECnTp2YOnXqZ80HbGdlx472OyictTCP3jyi3sp6PIt5ZsCqhRBCpGcmC7KhoaEMGTKEVatWffTT4SNHjiQyMlL/Cg2VuzciuddvXzN071AARlcfTd7MefX7rly5Qq1atQgLC6N48eIcP36cfPnymapUs2RjY8OaNWv4+uuvAfjuu+/46quv0Gq1KW4zi10W9nbai6eTJzde3KD+yvpExEUYqGIhhBDpmcmCbFBQEM+ePaNs2bJYWlpiaWlJQEAAs2bNwtLS8r3/MNrY2ODk5JTsJcQ/jTw4kuexzymStQjDKg/Tb798+TK1atXixYsXlC1bloCAAHLkyGHCSs2XWq1m+vTp/PrrrwDMmjWLdu3aERcXl+I2vZy9ONjlINntsxMcHkyjVY2I1kQbqmQhhBDplMnmkY2KiuLhw4fJtnXv3p3ChQszfPjwj/q491PmGRPp3+lHp6m8uDIKCgHdAqieuzqQFGJr166tD7EHDhyQBwsNZO3atXTt2hWNRkO1atXYtm3bZ13bS08vUXNZTV7HvaZO3jrs6LADW8uMO5+vEEJ8iFar5fLlywCUKFEi3czQ9Cn5zmR3ZB0dHSlevHiyl729PVmyZJExi+KTJeoS6bejHwoKXUt11YfYK1eu6EOsj4+PhFgDa9euHXv27MHJyYljx45RrVo1njx5kuL2SmYvye6Ou3GwduDg/YM0X9ecuMSU3+kVQoj0zMLCgtKlS1O6dOl0E2I/lclnLRDCEGYHzubi04tkts3M1LpJK3TduXMHPz8/fYjdv3+/hFgjqFWrFsePH8fDw4OrV69StWpV7t69m+L2fHP5sr39djJZZmLPnT0SZoUQQvyrdL1ErcgYHr15RJE5RYjWROuXPA0NDaVatWo8fPiQkiVLcvjwYVxdXU1darr24MED/Pz8uHv3Lu7u7uzdu/ezFjk58uAI/qv9iU2IpX7++mxtt1WGGQghxD/odDru3LkDgLe3d4ZcojZ9nLHI0L7a8xXRmmgq5apEjzI9ePbsGXXr1uXhw4cUKFCAffv2SYhNBXny5OH48eOULFmS8PBwatSowcmTJ1PcXs08NdnVYRd2VnbsvbuXpmub8jbhrQErFkII85aYmMiaNWtYs2YNiYmJpi7HJCTICrO2584eNl3fhIXKgj+++IOoN1HUr1+fmzdv4unpyYEDB8iePbupy8ww3N3dCQgIoEqVKkRERODn58fevXtT3F6NPDX0YXbf3X3UX1mfyLhIA1YshBDCnEmQFWYrQZvAV3u+AmCw72AKuRSiWbNmBAcHky1bNg4cOICXl5dpi8yAXFxc2LdvHw0bNuTt27c0btyY9evXp7i9GnlqsLfTXpxsnDgWcoyay2vyNPqp4QoWQghhtiTICrM15+wcbr68iZudG6OrjaZLly4cOXIER0dH9uzZQ8GCBU1dYoZlZ2fH1q1badu2LQkJCbRr144FCxakuL2qXlUJ6BZANvtsBIcHU3VpVR5EPDBcwUIIIcySBFlhlp7HPGfckXEATKoziYk/TGT9+vVYWVmxZcsWypQpY9oCBdbW1qxatYp+/fqhKAp9+/ZlypQpKV7StrR7aY53P05u59zceXWHKkuqEBwebNiihRBCmBUJssIs/XDoByLjIymboyyvD71mxowZACxdupQ6deqYuDrxNwsLC+bOncuoUaOApGWmv/vuuxSH2QJZCnCixwmKuRXjSdQTqi6pyl83/zJkyUIIIcyIBFlhdoLDg1l4fiEATa2b8t233wHwyy+/0LFjR1OWJt5DpVLx448/6pe0nTZtGr17937vMtQfI6dTTo51P0advHWISYih2dpm/Hry1xSHYyGEEOZLgqwwK4qiMGTPEBQU6rrXZXL/yQB8+eWXfPPNNyauTnzI0KFDWbJkCWq1msWLF9O2bVvi4+NT1FbmTJnZ3XE3fcr2QUHhm/3f0HdHX+ITU9aeEEKYIwsLCxo2bEjDhg0z7MpesiCCMCvrr66n7ca22FrYYr/Unpf3XuLv78+2bdsy7P/E5mbLli20a9cOjUaDn58fW7ZswcHBIUVtKYrCzNMzGbZvGAoK5T3Ks6H1BnK75DZw1UIIIVKLLIgg0qXYhFi+3f8tAE6XnHh57yUlS5ZkzZo1EmLNSPPmzdm1axf29vYcOHAAPz8/Xr58maK2VCoVX1f6mh0ddpDZNjNnn5ylzPwy7Lq9y8BVCyGESIskyAqzMe3kNEIiQ7CNt+XZ1me4u7uzY8cOHB0dTV2a+ER16tTh0KFDuLq6cubMGWrUqMGTJ09S3F6jAo043/c85T3K8zruNf6r/Rl5YCQarcaAVQshRNqi0+l48OABDx48QKfTmbock5AgK8xCaGQoU45PASDurzgyWWZi+/bteHp6mrgykVIVKlTg2LFj5MyZk6tXr1KlShX9muEpkcclD8e6H2Ng+YEATDkxhYqLKnL56WVDlSyEEGlKYmIiy5cvZ/ny5bJErRBp2fADw3mb+BYeAldhxYoVlCtXztRlic9UtGhRjh8/jre3Nw8ePKBq1apcunQpxe3ZWNowu9Fs1rdaj2smVy6EX6DcwnJMOT6FRF3G/CEvhBDpmQRZkeYdDznOmitrQAF2w+jRo2nVqpWpyxIGkidPHo4fP06pUqV4+vQp1atX58SJE5/VZutirbk64CqNCzZGo9Uw8uBIqiypwrkn5wxUtRBCiLRAgqxI07Q6Lf3/6p/0hyBoVrEZ48aNM2lNwvCyZ8/OkSNHqFKlCpGRkdStW5c9e/Z8VpvuDu5sa7eNZU2X4WTjRODjQCosrEDf7X15EfvCQJULIYQwJQmyIk1beHYhV15egTgo9LgQK1asQK2Wb9v0yMXFhX379tGwYUPevn1LkyZNWLly5We1qVKp6Fq6K9cHXqdjiY4oKCw4v4CCvxfk9zO/y7yzQghh5iQRiDQrMi6SoTuHAmAXaMeujbtkhoJ0zs7Ojq1bt9K+fXsSEhLo3LkzEydO/OxVuzwcPVjZYiVHux2lZPaSvI57zeA9gynwewEWBC2Q2Q2EEMJMSZAVaVa7ee14a/EWXsD6b9eTL18+U5ckUoG1tTUrV67ku++Slh4eM2YMPXr0QKP5/LBZLXc1gvoEMc9/HjkdcxL6JpS+O/pSaHYh5p+bT2xC7Gf3IYQQIvXIyl4iTdp2fBvN9jUDC+ik6sSfY/40dUnCBObPn8/AgQPRarXUrl2bjRs3kjlzZoO0HZcYx4KgBUw6NomnMU8ByGybmZ5lejKg/ADyZs5rkH7SA52iI0YTQ0xCDG8T3gJJwzYALNWWOFo74mjjiFol90aESE1arZbTp08DULFixXSzONCn5DsJsiLNiYiIIOe3OYnNFYvbazfCp4fLuNgMbM+ePbRu3Zro6Gi8vb3566+/KFKkiMHaj02IZUHQAn4P/J17r+8BoEJFfe/6dCjegWaFm+Fokz6HtLyJf8PdV3cJiQzh0ZtHhL4J5XHUY57FPON5zHNexL7g5duXH3WnWoUKRxtHMttmxtPZEy9nLzydPMnrkpcS2UtQPFtxnGzk57QQ4r9JkBVmS6fTUblrZc54nwEdnOlyhgr5K5i6LGFily5dokmTJjx8+BBHR0fWrFmDv7+/QfvQ6rTsvrObWWdmsf/efv32TJaZaFyoMc0LN6duvrpkscti0H6NLTIuktuvbnPn1R396+8/P4t59kltqVCRySoTAIqioKCQoE1Aq2g/6v15XfJSJkcZqntVp0aeGpTIVgILdfq4gySEMBwJssJsTZoyiVGPR0FW6JSvE392liEFIsnz589p1aoVR48eRaVSMWnSJIYPH67/iNuQbr28xZrLa1h9ZTW3Xt7Sb1ehopxHOernr0+13NUo71GezJkMM9QhpRRF4Xnsc+6+usvd13f1//07rP7XVGNudm7kcclDLqdceDp5ktMpJ9nts+Nm74abnRtZ7LLgaO2IvbU9mSwzvXO9FUUhLjGON/FviIyP5EXsC0IjQwl9E0pIZAi3X93m8tPLPI56/E7fLrYu1MxTk2aFmvFFwS/M7pcEIUxNp9MRFhYGQI4cOdLNp5cSZIVZOnz4MHVG1UGpr+CodiT021CcbZ1NXZZIQzQaDUOGDOGPP/4AoHnz5ixduhRnZ+N8nyiKwvmw86y7uo49d/Zw+dm7y90WcC1A+ZzlKeZWjAKuBSiYpSDert7YW9sbpIYYTQzh0eGERYcl/TcqjJDIkKTQ+vou917fI1oT/cE23B3c8Xb1xtvVmwKuBfRfe7t6p9rH/S9jX3Lp6SXOPD5DwMMAjoccT1a3hcqCarmr0apIK9qXaI9rJtdUqUsIc6bRaJg8eTIAI0eOxNra2sQVGYYEWWF2Hj9+TOkqpXnR/gXYwsLGC+lVtpepyxJp1B9//MGQIUPQaDTkz5+fjRs3Urp0aaP3+/jNY/bd3ceB+wc48+gMd1/f/ddjnW2cye6QHXcHd7LZZ8Peyh47KzvsrOywsbABkh6i0ik6EnWJRGuiiU6IJloTzeu3r/Xh9b9CKiTdKc7llIv8rvnJnznpVSBLUmDNnzl/mhzjm6hL5ELYBXbe3snWG1u5+PSifp+NhQ3NCjejR5ke1MlbR4YfCPEvJMhKkBVpQEJCAjVr1uSk60koB6Wzl+Zcn3Pyj5f4oLNnz9K6dWsePnyIjY0Ns2fPpmfPnkYZavBvXsa+5NyTcwSFBXHz5U1uvbzFrZe3ePX2lUH7sbOyI4dDDtwd3HF3cE8KrZnz64Nrbpfc2FraGrTP1Hb/9X223NjC8ovLufT0kn57Xpe8fFXxK7qX7p4mA7kQpiRBNoVB9t69e2liTk8JsunD119/zcy1M6EvoIJj3Y9R1auqqcsSZuDVq1d06dKFnTt3AtCqVSv++OMPsmQx7VjLiLgIwqPDCY8O52n0U57HPic2IVb/ikuMQ4UKtUqNWqXGUm2Jg7WD/uVs64y7g7s+vDpYO6RqQDclRVG4EH6BpReWsuryKl7HvQaS7nD39enLkIpD8HD0MHGVQqQNEmRTGGTVajU1atSgZ8+etGrVCltb09wJkCBr/rZt20azZs2gO5Ab2hVvx5qWa0xdljAjOp2OX375hdGjR5OYmIiHhwfLli2jbt26pi5NkBRMo6Ojefr0KeHh4bx8+ZLo6Giio6OJiYlBo9GgUqn0L1tbW5ydnXF2dsba3ppTsadYfX81dyOShnHYWNjQv1x/RlQdQXaH7CY+OyFMS4JsCoNscHAwS5cuZc2aNWg0Gtq2bUvPnj2pUCF1p0mSIGveQkNDKVWqFK89XkPrpGmObg66iaezp6lLE2bo3LlzdOzYkVu3kmYZGDJkCD/99BP29oZ56Ep82PPnz7l06RLXr1/nzp073L17l7t37/Lw4UNiYz9zxTQV2JawRVVNxVu3pAUZrFXWtM7dmp8a/URut9wGOAMhzI8E2c8cI5uYmMhff/3FsmXL2LNnDwULFqRHjx507twZNze3lDb70STImq/ExERq1arF8TPHsf7aGo2dhvE1xzOmxhhTlybMWExMDN9++y3z5s0DIG/evPzxxx/Uq1fPxJWlL2FhYZw5c4bAwECCgoK4dOkS4eHhH3yPg4MD7u7uZMmSBQcHB/3LysoqaU5aRUGn0xEXF8ebN2+IjIwkMjKSp0+f8vr16/9rKD9QC8j1///8FtxvutPAtQEVK1SkevXqFC5cOMMMxRAZmwRZAz3sFR8fz9y5cxk5ciQajQZra2vatGnDzz//TI4cOT63+X8lQdZ8jR07lgkTJmBd1xpNFQ1ezl5cH3gdOys7U5cm0oHdu3fTt29fQkNDAejUqRMzZswga9asJq7MPD19+pT9+/ezf/9+Dh8+rL+u/6RSqcifPz/FihWjQIEC5M+fn/z585M3b17c3d1xcHBIcf+xsbE8fvyY0NBQbt26xbXr1zgWfowr2a+QmCUx6aAXwF7gdtJ8mrVr16Z27do0bNjQqP8OCWFKWq2WY8eOAVCtWjVZovZTnTt3jiVLlrB27Vrs7e3p2rUrPXv25NGjR4wfP543b94QGBiY0ub/kwRZ83TkyBFq166N4qhgPcwajaJhXat1tCnWxtSliXQkKiqKH374gd9//x1FUXB1dWXcuHH069cPKysrU5eXpr19+5bjx4+zb98+9u/fz8WLF5PtV6vVFCtWDF9fX8qXL0/p0qUpVqxYqg/j0Oq0TD8ynZ9O/URkYmRSbXfU6LbrIPL/jqtQoQJNmzalSZMmFCtWTO7WCpHGGT3ITp8+naVLl3Lz5k0aNWpEr169aNSoUbIVJR49ekSePHlITEz813bmzZvHvHnzePDgAQDFihVjzJgxNGzY8KPqkCBrfl68eEGpUqV48uQJeb/Ny337+1TPXZ0jXY/IPy7CKM6cOUPv3r25fDlpMYPChQszbdo0GjVqJN9z/5+iKFy6dIn9+/ezb98+jh07RlxcXLJjypYtS926dalbty6+vr6fdYfV0CLjIvnx6I/8duY3EnQJ2KptqZpQlYi9EZwLPJfs2KJFi9KhQwfat2+fJmbfEUK8y+hBtkCBAvTo0YNu3br960c2Go2GNWvW0LVr139tZ/v27VhYWFCgQAEURWH58uVMnTqVCxcuUKxYsf+sQ4KseVEUhcaNG7Nz5068qnkRUicEFSrO9z1PaffSpi5PpGOJiYksWrSI0aNH8+JF0pKttWvXZty4cVSrVs3E1ZlGWFiYfrjA/v37efr0abL9uXLlom7dutSrV486deqkynMPn+va82v03dGX4yHHASiboyyTK07mwekHbNu2jQMHDqDRaPTHV6xYkR49etCuXTscHWWOWmF+FEXh+fPnALi5uaWbX84/Kd8pKXD//n1Fq9W+s12n0ykPHz5MSZN6mTNnVhYtWvRRx0ZGRiqAEhkZ+Vl9itQxY8YMBVCsba2VwjMKK4xD6fNXH1OXJTKQiIgI5dtvv1Wsra0VQAGUWrVqKYcPHzZ1aUYXExOj7NmzRxk6dKhSvHhx/fn//bKzs1P8/f2VmTNnKteuXVN0Op2pS04RrU6rLAxaqLhMcVEYh2I5wVIZf2S8oknUKK9fv1aWLFmi+Pn5KWq1Wn/u9vb2Ss+ePZVTp06Z7XmLjCk+Pl4ZN26cMm7cOCU+Pt7U5RjMp+S7FN2RtbCwICwsjGzZsiXb/vLlS7Jly4ZWq/3UJtFqtWzYsIGuXbty4cIFihYt+s4x8fHxxMfH6//85s0bPD095Y6sGQgKCqJSpUokJCTQcVpHVkWvwtnGmdtf3sbNPu3f6RHpy4MHD5g8eTJLly4lISEBgEqVKjFo0CBatWqVLp781el0XLx4UT/O9dixY8nuRqpUKnx8fKhXrx5169alUqVK2NjYmLBiw3oa/ZSBuway6fomAHxy+LCi+QqKuiX92xIWFsaqVatYtGgRN2/e1L/Px8eHIUOG0KZNm3R1PUT6JLMWkLI7siqVSnn69Ok72x88eKDY2dl9UluXLl1S7O3tFQsLC8XZ2VnZuXPnvx47duzYd+4iIHdk07w3b94o3t7eCqB80fILJdvUbArjUKafnG7q0kQG9/DhQ6V///7J7tBmy5ZNGTVqlHL//n1Tl/dJdDqdcvfuXWXJkiVK+/btFTc3t3d+Vnp6eio9e/ZU1q1bp7x48cLUJRudTqdTVl1apb87azPRRplxakayu646nU45evSo0rlzZ8XGxkZ/rbJnz66MGzdOef78uQnPQIgPkzuyn3hHdujQoQD89ttv9O7dGzu7/5sqSavVcubMGSwsLDhx4sTHNolGoyEkJITIyEg2btzIokWLCAgIkDuy6UjXrl1ZsWIFnp6eNJndhDkX5lAoSyEu9b+EtUX6+O1RmLewsDAWLlzI/PnzefLkiX57xYoVadeuHa1bt8bDI20ti6ooCrdu3SIgIICAgACOHj3Ko0ePkh3j4OBArVq19HddCxYsmG7G0H2KJ1FP6PVXL3bf2Q1AowKNWNp0Kdnsk3+q+Pz5cxYsWMDcuXP13weZMmWid+/eDBs2DC8vr1SvXYgPkTuyn/iwV61atQAICAigUqVKyS6YtbU1efLk4ZtvvqFAgQIpLB38/PzInz8/8+fP/89j5WGvtG/Dhg20adMGtVrNn7v+pGtgVxJ1iezuuJsG3g1MXZ4QySQkJLBt2zbmzZvH4cOH+fvHo0qlokKFCvqn9itWrJjq/2DExMQQHBxMUFAQJ06c4OjRo+8sQmBlZUX58uWpVasW9evXx9fXN938w/a5FEVh7tm5DNs3jHhtPNnts/Nn8z+pm//dpYwTEhLYuHGj/uFjAEtLSzp27MioUaM+6984IQxJgmwKZy3o3r07v/32m1HCY+3atfHy8mLZsmX/eawE2bTtyZMnlChRglevXjF69GjOFz7Pzts78S/gz44OO0xdnhAfFBYWxsaNG1m7di0nT55Mts/e3p5y5cpRtmxZ/St//vwGGVP59u1b7t+/z/Xr1/Wvixcvcv36dXQ6XbJjbWxsqFixIjVq1KBGjRpUrFgx2Sdl4l2Xn16m/ab2XH1+FYBR1UYxvuZ4LNTvTiSvKAoHDhxgypQpHDp0CEh6RqRTp0788MMPeHt7p2rtQvwvCbIGWtkrpUaOHEnDhg3x8vIiKiqK1atX8/PPP7N3717q1n33t+T/JUE27VIUhYYNG7J37158fHwYu3IsTdY1wVJtydUBVymYpaCpSxTioz169IgDBw6wf/9+Dhw4wLNnz945RqVS4eHhQZ48ecidOzeurq44Ozvj7OycbGonRVFITEwkMjKSiIgIIiIieP78OaGhoYSGhuqnB3ufHDlyUK5cOSpUqECNGjWoUKGCPJCUAm8T3jJ071D+CPoDgLr56rK65Wqy2v37ym9nzpxh4sSJ7Ny5E0gKtF26dGHChAnkypXrX98nhDFJkP2EINuiRQuWLVuGk5MTLVq0+OCxmzdv/qhCe/bsycGDBwkLC8PZ2ZmSJUsyfPjwjwqxIEE2LZs7dy4DBw7E1taW02dP0/ZIW26+vMmwSsOYVm+aqcsTIsV0Oh3Xrl0jKCiIoKAgzp8/T3BwMDExMQbrw9HRkSJFilCkSBEKFy5M8eLF8fHxkaVWDWz15dX03t6b2IRYvJy92NRmE+U8yn3wPYGBgYwfP55du3YBYGtry9dff83w4cNxdnZOjbKF0NNqtRw8eBCAOnXqyBK1H9K9e3dmzZqFo6Mj3bt3/+CxS5cu/fhqP4ME2bTp5s2blClThrdv3zJr1iwSyyUydN9Qstln49agWzjbyg97kb4oisKLFy948OAB9+/fJyQkhIiICCIjI4mMjCQ6OhpA/6CVpaUlzs7OuLi44OLigqurK7ly5cLLywtPT09cXFwy5ENZpnD56WVarG/BnVd3sLawZk6jOfQq2+s/33f69GmGDx/O0aNHAciaNStjxoyhb9++6eaumBCmYjZDCz6XBNm0JyEhgSpVqnD27Fnq1q3Lik0rKDynMJHxkSxqvIieZXuaukQhhEgmMi6Srlu7su3mNgB6lO7B7EazyWSV6YPvUxSF7du3M3z4cG7cuAGAt7c3U6ZMoUWLFvLLiBAp9Cn5Tp2SDt6+fUtsbKz+zw8fPmTmzJns27cvJc2JdGTSpEmcPXsWFxcXli5dyg+HfyAyPhKfHD50K93N1OUJIcQ7nG2d2dx2M5NqT0KtUrMkeAlVl1blQcSDD75PpVLRpEkTLl++zB9//EH27Nm5c+cOrVq1okqVKpw6dSp1TkBkWIqi6Mfam/F9yc+SoiDbtGlTVqxYAUBERAQVKlTg119/pWnTpsybN8+gBQrzERgYyMSJEwGYN28eT3jCkgtLAJjVcNZ7nwoWQoi0QK1SM7LaSPZ22ktWu6ycDztP2fllOXjv4H++19LSkr59+3L79m3GjBmDnZ0dp06donLlynTv3v29DwcKYQgJCQn89ttv/Pbbb/pVCjOaFAXZ8+fPU61aNQA2btyIu7s7Dx8+ZMWKFcyaNcugBQrzEBsbS+fOndFqtbRv3542bdsweM9gFBQ6l+xMZc/Kpi5RCCH+k18+P4L6BFHeozyv417TYFUD5p/773nNIekhvfHjx3P79m169OgBwLJlyyhUqBDz5s1L0fLtQogPS1GQjY2N1U8ns2/fPlq0aIFaraZixYo8fPjQoAUK8/Ddd99x69YtcubMyZw5c1h5aSWnH53GwdqBKX5TTF2eEEJ8NC9nL452P0qHEh1I1CXSb2c/vtrzFYm6xI96v4eHB4sXL+bUqVOUKVOGiIgIBgwYQMWKFTl79qyRqxciY0lRkPX29mbr1q2Ehoayd+9e6tWrB8CzZ8/koasMaM+ePcyZMwdImrHC0s6S4QeGA/BDtR/wcExbS3sKIcR/sbW0ZWXzlfxY60cAfjvzG03WNCEyLvKj2/g7uP7+++84Oztz7tw5fH196devH69evTJW6UJkKCkKsmPGjOGbb74hT548+Pr6UqlSJSDp7myZMmUMWqBI216+fKn/CG3w4MHUrVuXH4/+SHh0ON6u3nxV8SvTFiiEECmkUqkYVX0UG1pvIJNlJnbf2U3lJZW59/reR7dhYWHBoEGDuHnzJp07d0ZRFObPn0/BggVZtmxZhn1ARwhDSVGQbdWqFSEhIZw7d449e/bot9epU4cZM2YYrDiRtimKQv/+/QkLC6Nw4cJMmTKFWy9vMeN00vfAzPozsbGUVYeEEOatVdFWHOt+DA9HD649v4bvIl+OPTz2SW1kz56dFStWEBAQQLFixXj58iXdu3enQYMGMiRPiM+QoiAL4O7uTpkyZVCr/6+JChUqULhwYYMUJtK+1atXs2HDBiwtLVm5ciWZMmXi671fk6BLoFGBRvgX9Dd1iUIIYRA+Hj4E9grEJ4cPL2JfUGdFHZYFL/vkdqpXr86FCxeYMmUKNjY27Nu3j+LFizN37lx0Op3hCxcinUtRkI2JiWH06NFUrlwZb29v8uXLl+wl0r+QkBAGDhwIwNixY/Hx8WHnrZ3sur0LK7UVM+rLnXkhRPqS0yknR7sfpVXRViToEui+rTtjDo/55OEBVlZWDB8+nIsXL1KlShWio6MZOHAgtWrV4vbt20aqXqRHarWacuXKUa5cuWQ3FjOSFK3s1b59ewICAujcuTM5cuR4Z/WSIUOGGKzAD5GVvUxDp9Ph5+fH4cOHqVixIseOHUOLluLzinPn1R2+q/wdP9f92dRlCiGEUegUHWMOj+GnYz8B0LlkZxY1WYS1xacvTavT6Zg7dy4jRowgJiYGW1tbfvzxR7766issLGTubZExGX2JWhcXF3bu3EmVKlVSXKQhSJA1jZkzZ/L1119jZ2fHxYsX8fb25ufjPzPi4AjcHdy5NegWjjaOpi5TCCGMatH5RfTb0Q+toqV23tpsbrMZZ1vnFLX14MEDevfuzYEDB4CkoXpLly6laNGihixZCLNg9CVqM2fOjKura4qKE+bt6tWrjBgxAoDp06fj7e3Nk6gnTDyatKLXL36/SIgVQmQIvcr2YkeHHThYO3Do/iGqLq1KaGRoitrKkycP+/btY9GiRTg7OxMYGIiPjw+zZs2SsbPiXymKQkxMDDExMRl2BowUBdmJEycyZswYYmNjDV2PSMM0Gg2dOnUiPj6eRo0a0adPHwC+2/8dMQkxVMpViY4lO5q4SiGESD0NvBtwtNtRcjjk4MqzK1RcXJHg8OAUtaVSqejZsydXr16lYcOGxMXFMWTIEBo0aMDjx48NW7hIFxISEpg2bRrTpk2TJWo/xa+//srevXvJnj07JUqUoGzZssleIn0aP348wcHBZMmShcWLF6NSqTjy4AirLq9ChYpZDWehVmXMweZCiIyrTI4ynO51mmJuxXgS9YRqS6ux987eFLeXM2dOdu7cyZw5c8iUKRP79++nRIkSrFu3zoBVC5E+WKbkTc2aNTNwGSKtO3HiBFOmJC01u2DBAtzd3UnQJjBwV9LMBf3K9aOcRzlTliiEECbj5ezF8R7HabGuBYcfHMZ/tT8LGi+gR5keKWpPpVIxYMAA6tSpQ6dOnTh37hzt2rVj+/btzJ49GxcXF8OegBBmKkUPe6UV8rBX6oiKiqJ06dLcu3ePrl27smzZMgCmnpjKdwe+w83OjZuDbpI5U2bTFiqEECam0Wro+VdPVl5aCcDEWhMZVW3UO7P7fIqEhAQmTpzITz/9hE6nw9PTk+XLl1OrVi1DlS3MlEajYfLkyQCMHDkSa+tPnzkjLTL6w14AERERLFq0iJEjR+rXjD5//ryM40mHhg4dyr1798idOze//fYbAKGRoYwPGA/AL3V/kRArhBCAtYU1K5qt4Puq3wMw+vBovtz9JVqdNsVtWllZMWHCBE6cOEH+/PkJDQ2lTp06DB8+PMOOixTibykKspcuXaJgwYL8/PPPTJs2jYiICAA2b97MyJEjDVmfMLHt27ezaNEiVCoVy5cvx9k5aWqZr/d+TUxCDFU8q9ClVBcTVymEEGmHSqXipzo/MavBLFSomHN2Du02tSMuMe6z2q1YsSLBwcH07t0bRVH45ZdfqF69Og8ePDBM4UKYoRQF2aFDh9KtWzdu376Nra2tfnujRo04evSowYoTpvXs2TN69eoFwLBhw6hRowYAe+/sZdP1TVioLJjrP1ce8BJCiPf40vdL1rZai7WFNRuvbaThqoZExkV+VpsODg4sWLCADRs24OzszOnTpylTpgybN282UNVCmJcUJZCzZ8/St2/fd7bnzJmT8PDwzy5KmJ6iKPTp04dnz55RokQJfvzxRwDiEuMYtHsQAF9W+JKS2UuaskwhhEjT2hRrw+6Ou3G0duTIgyPUWFaDsKiwz263VatWBAcH4+vrS0REBC1btmTgwIHExX3eXV9hXtRqNaVKlaJUqVIZdonaFJ21jY0Nb968eWf7rVu3cHNz++yihOktW7aMbdu2YWVlxZ9//omNjQ2Q9IDXnVd3yOGQg/G1xpu4SiGESPtq561NQLcAsttn5+LTi1ReUplbL299drt58uTh2LFjDB8+HIC5c+fi6+vLjRs3PrttYR4sLS1p1qwZzZo1w9IyRRNRmb0UBdkmTZowYcIE/SBzlUpFSEgIw4cPp2XLlgYtUKS++/fvM3jwYAB+/PFHSpUqlbT99X0mHZ8EwPT603GykZkihBDiY5TJUYaTPU/i7erNg4gHVFlShbOPz352u1ZWVkyZMoU9e/bg5ubGpUuX8PHxYfny5QaoWoi0L8ULIkRHR+Pm5sbbt2+pUaMG3t7eODo68tNPPxm6RpGKtFotXbp0ITo6mmrVqjFs2DAgaajBl7u/JC4xjtp5a9O2WFsTVyqEEOYlX+Z8nOhxAp8cPryIfUGt5bU+a+GEf6pfvz4XL16kdu3axMbG0q1bN3r06MHbt28N0r5ImxRFQaPRoNFoMuwStZ81j+yJEye4ePEi0dHRlC1bFj8/P0PW9p9kHlnD+/nnnxkxYgSOjo5cunSJPHnyALDtxjaarWuGldqKS/0vUThrYdMWKoQQZioqPoqW61uy/95+LNWWLG26lE4lOxmkba1Wy+TJkxk7diw6nY5SpUqxceNGvL29DdK+SFtkHtkU3JHV6XQsWbKEL774gr59+zJv3jyOHz/OkydPMuxvA+nFxYsXGT16NAC//fabPsTGJsQyZM8QAIZVGiYhVgghPoOjjSM7OuygQ4kOJOoS6bylM9NOTjNI2xYWFvzwww/s27cPNzc3Ll68iI+PD1u3bjVI+0KkNZ8UZBVFoUmTJvTq1YvHjx9TokQJihUrxsOHD+nWrRvNmzc3Vp3CyOLi4ujUqRMJCQk0a9aMbt266ff9ePRHHkY+xMvZix+q/2C6IoUQIp2wtrDmz+Z/8nXFrwH4dv+3DNs7DJ2iM0j7derU4cKFC1SuXJk3b97QvHlzhg8fTmJiokHaFyKt+KQgu2zZMo4ePcrBgwe5cOECa9asYe3atVy8eJEDBw5w6NAhVqxYYaxahRH98MMPXLlyhWzZsrFgwQL9coqXnl5i6smpAPzW4Dfsre1NWaYQQqQbapWaX+v9yi9+vwAw/fR0umzpgkarMUj7OXPm5MiRI3z9dVJY/uWXX/Dz85NpMkW68klBds2aNXz//ffvXd+5du3ajBgxglWrVhmsOJE6jhw5wvTp0wFYvHixfgo1rU5L7+29SdQl0rxwc5oVbmbCKoUQIv1RqVR8W+VbljdbjoXKglWXV9F4TWOiNdEGad/Kyorp06ezYcMGHB0dCQgIoEyZMrJ4kUg3PinIXrp0iQYNGvzr/oYNG3Lx4sXPLkqknsjISLp27YqiKPTu3ZsvvvhCv2/O2TkEPg7EycaJ2Y1mm7BKIYRI37qU6sL29tuxs7Jj39191Fpei2cxzwzWfqtWrTh79izFihUjPDyc2rVrM3PmTHm2RZi9Twqyr169Inv27P+6P3v27Lx+/fqzixKpZ/DgwYSEhJAvXz79XVmAkMgQvj/4PQA/+/2Mh6OHqUoUQogMoWGBhhzqcogsmbJw7sk5qi6pyv3X9w3WfqFChThz5gydOnVCq9Xy9ddf07VrV5miS5i1TwqyWq32gytHWFhYfNJA8smTJ1O+fHkcHR3Jli0bzZo14+bNm59SkvgMmzZtYsWKFajVav78808cHByApIf6Bu4aSExCDFU8q9DHp4+JKxVCiIzBN5cvJ3qcILdzbm6/uk3lJZW5GG64Tzrt7e1ZsWIFM2fOxMLCgj///JNq1aoRGhpqsD5E6lGr1RQtWpSiRYtm2CVqP2keWbVaTcOGDfXLlf6v+Ph49uzZg1ar/aj2GjRoQLt27ShfvjyJiYl8//33XLlyhWvXrmFv/98PFck8sikXFhZGiRIlePnyJd9//32yhSw2XN1Am41tsFJbEdwvmKJuRU1YqRBCZDxPop7QYGUDLj+7jJONE9vabaNmnpoG7ePQoUO0adOGly9f4ubmxsaNG6levbpB+xAiJT4l331SkO3evftHHbd06dKPbTKZ58+fky1bNgICAj7qfyYJsimjKAr+/v7s3r2bMmXKcPr0af0kyq/fvqbInCI8jXnK2BpjGVdznGmLFUKIDCoiLoIma5pwLOQY1hbWrG6xmpZFDbsM/IMHD2jevDnBwcFYWloyY8YMBg4cqJ+5RghTMFqQNbY7d+5QoEABLl++TPHixd/ZHx8fT3x8vP7Pb968wdPTU4LsJ5o3bx4DBgzAxsaG8+fPU7To/91x7f1XbxZdWEThrIUJ7huMjeX7774LIYQwvrcJb+mwuQNbb2xFhYq5/nPpV66fQfuIjY2lV69erFmzBki6aTV37lxsbW0N2o8QH8uoK3sZi06n46uvvqJKlSrvDbGQNKbW2dlZ//L09EzlKs3fjRs3GDZsGABTpkxJFmIDHgSw6MIiABY2XighVgghTCyTVSY2tt5In7J9UFDov7M/446MM+hsA3Z2dqxatYpp06ahVqtZunQpNWrU4NGjRwbrQxiHRqNh/PjxjB8/Ho3GMPMPm5s0E2QHDhzIlStXWLt27b8eM3LkSCIjI/UvGZz+aTQaDR07duTt27fUrVuXwYMH6/fFJcbRZ0fSQ119ffpS1auqqcoUQgjxDxZqC/744g/GVB8DwPiA8fTf2R+t7uOeR/kYKpWKYcOGsWfPHlxdXQkMDKRcuXKcOHHCYH0IYQxpIsgOGjSIHTt2cPjwYXLlyvWvx9nY2ODk5JTsJT7e+PHjOX/+PK6urixbtizZE46Tjk3i1stb5HDIwRS/KSasUgghxP9SqVSMrzWeuY3mokLF/KD5tNnYhrjEOIP2U7duXc6ePUuJEiV4+vQptWrVYv78+QbtQwhDMmmQVRSFQYMGsWXLFg4dOkTevHlNWU66duzYMSZPngzA/Pnz8fD4v3lhrzy7wpTjSeH194a/42LrYooShRBC/If+5fuzvvV6rC2s2Xx9M/VX1iciLsKgfeTLl49Tp07Rpk0bEhIS6NevH3369MmwH12LtM2kQXbgwIGsXLmS1atX4+joSHh4OOHh4TI5s4FFRkbSuXNnFEWhW7dutGrVSr8vQZtAt63dSNAl0KRQE1oUaWHCSoUQQvyXVkVbsbfTXpxsnDj68CjVl1bnSdQTg/Zhb2/P2rVrmTJlCiqVioULF1KnTh2ePTPcamNCGIJJg+y8efOIjIykZs2a5MiRQ/9at26dKctKdwYPHszDhw/Jly8fs2bNSrZvyvEpBIUFkdk2M3/4/yFTrgghhBmomacmAd0CcHdw5/Kzy1ReXJlbL28ZtA+VSsXw4cPZuXMnTk5OHD9+nPLlyxMcHGzQfoT4HCYfWvC+V7du3UxZVrqyfv36ZKt3OTo66vcFhwcz4egEAGY3mk0OxxymKlMIIcQnKu1empM9TuLt6s3DyIdUWVKFwMeBBu+nYcOGnDlzhoIFCxISEkKVKlXYuHGjwfsRIiXSxMNewjgePXpE3759ARg1ahSVK1fW79NoNXTd2pVEXSLNCzenffH2pipTCCFECuXNnJcTPU7gk8OHF7EvqL28Nnvv7DV4P4ULF+b06dPUr1+f2NhYWrduzdixY9HpdAbvS3w8tVpNgQIFKFCggCxRa45kZa9/p9PpqFu3LocOHaJ8+fKcOHECKysr/f7Rh0bz47EfyWqXlasDrpLNPpsJqxVCCPE5ouKjaLm+Jfvv7cdSbcmypsvoWLKjwftJTExk+PDhTJ8+HYDmzZuzYsUKHBwcDN6XyLjMckEEYVgzZ87k0KFD2NnZsXLlymQh9uzjs0w+njSDwTz/eRJihRDCzDnaOLKjww7aF29Poi6RTls6MePUDIP3Y2lpya+//srSpUuxtrZmy5YtVK5cmQcPHhi8LyE+hgTZdCg4OJiRI0cCMGPGDAoWLKjfF5cYR9etXdEqWtoVb0eroq3+rRkhhBBmxNrCmpUtVjLEdwgAQ/cNZfj+4QZdBexv3bp148iRI7i7u3P58mXKly9PQECAwfsR4r9IkE1noqOjadeuHRqNhiZNmtC7d+9k+8ceHsv1F9fJbp+d2Q1nm6hKIYQQxqBWqZlRfwZT6iTNDf7LyV/otq0bCdoEg/dVqVIlzp49i4+PDy9evMDPz08WT0hlGo2GSZMmMWnSpAw7z68E2XTmyy+/5ObNm+TMmZMlS5Ykm07rVOgppp2aBsCCxgvIYpfFVGUKIYQwEpVKxfCqw1nadCkWKgtWXFxBs3XNiNHEGLyvXLlycezYMdq3b09iYiL9+vVj4MCBJCQYPjiL90tISMjQ11uCbDqyevVq/dKzq1atIkuW/wuqsQmxdN3aFZ2io0upLjQp1MSElQohhDC2bqW7sbXdVjJZZmLX7V3UWVGHl7EvDd5PpkyZWLVqFZMnT0alUjF37lzq1avHixcvDN6XEP9Lgmw6cefOHfr16wfA6NGjqVGjRrL93x/8ntuvbuPh6MFvDX4zRYlCCCFS2RcFv+Bgl4Nkts3MmcdnqLq0KiGRIQbvR6VSMWLECLZt24ajoyNHjhyhfPnyXL582eB9CfFPEmTTAY1GQ/v27YmKiqJatWr88MMPyfYfvn+Y384khdfFTRbjYutigiqFEEKYQiXPShzvcZxcTrm48eIGlRdX5sqzK0bpq3Hjxpw+fZr8+fPz4MEDKlWqxNatW43SlxAgQTZd+P777zl37hyurq6sWrUKS0tL/b5Xb1/ReUtnAHqX7U0D7wamKlMIIYSJFHUryqmepyjqVpTHUY+ptrQaRx8eNU5fRYsSGBhInTp1iImJoXnz5vz4449GmT1BCAmyZm737t38+uuvACxZsgRPT0/9PkVR6LujL4+jHlMwS0Fm1Df8nIJCCCHMQy6nXBzrfozKnpWJiIug7p91WXtlrVH6cnV1Zffu3Xz55ZdA0pC3du3aERsba5T+RMYlQdaMhYWF0bVrVwAGDRpE06ZNk+1fFryMjdc2Yqm2ZFWLVdhb25uiTCGEEGmEayZXDnQ+QIsiLdBoNbTf1J6pJ6Ya5W6plZUVs2bNYuHChVhZWbF+/XqqVq1KSIjhx+hmVCqVity5c5M7d+5ksxRlJLJErZnS6XTUq1ePgwcPUqpUKU6fPo2tra1+/51Xdyj9R2liEmKYXGcyI6qOMGG1Qggh0hKtTsuwfcP0z08MLD+Q3xr8hoXawij9HT9+nBYtWvD8+XOyZcvG5s2bqVKlilH6EuZPlqjNAH7++WcOHjyInZ0d69atSxZiE7QJdNzckZiEGGrmqcm3lb81YaVCCCHSGgu1BTMbzGR6vemoUDHn7Bxarm9JbIJxPvqvWrUq586do3Tp0jx79oxatWqxZMkSo/QlMhYJsmYoICCA0aNHAzB79mwKFSqUbP8Ph34g8HEgLrYurGi2wmi/YQshhDBvX1f6mvWt12NjYcO2m9uovbw2z2OeG6UvLy8vjh8/TqtWrUhISKBnz5589dVXJCYmGqU/kTFIkDUzT548oW3btmi1Wjp16kS3bt2S7d99eze/nPwFSJpqy9PZ8z2tCCGEEElaFW3FgS4HcM3kypnHZ6i0uBK3X942Sl/29vasX7+eCRMmAPDbb7/RsGFDXr16ZZT+0juNRsPUqVOZOnWqLFEr0r6EhATatGnD06dPKVGiBPPnz082uPtJ1BO6bO0CJI13alGkhalKFUIIYUaqelXlZI+T5HXJy93Xd6m8pDKnQk8ZpS+VSsXo0aPZvHkz9vb2HDhwAF9fX65fv26U/tK72NjYDD0bhARZMzJ8+HBOnDiBk5MTmzZtws7OTr9Pq9PSYVMHXsS+oLR7aabVm2bCSoUQQpibQlkLcarnKcp5lONF7Atqr6jNlutbjNZf8+bNOXXqFHny5OHOnTv4+vqyY8cOo/Un0icJsmZi/fr1zJiRNA/s8uXLKVCgQLL9E49OJOBhAA7WDqxrtQ5bS9v3NSOEEEL8q+wO2TnS9Qj+BfyJS4yj5fqW/H7md6P1V6JECc6ePUuNGjWIioqiSZMm/Pzzz7J4gvhoEmTNwPXr1+nRoweQdFe2WbNmyfYfuHeACQFJ443mfzGfglkKpnaJQggh0gl7a3u2tttKX5++KCgM3jOYYXuHodVpjdJf1qxZ2b9/P/3790dRFEaMGEGnTp14+/atUfoT6YsE2TQuKiqKFi1aEBMTQ61atfjxxx+T7Q+NDKX9pvYoKPQq04sOJTqYqFIhhBDphaXaknn+85hcZzIA009Pp+X6lsRoYozSn5WVFXPnzmXu3LlYWlqyevVqqlevzuPHj43Sn0g/JMimYYqi0KtXL27cuIGHhwdr1qzB0tJSvz8+MZ7WG1rzIvYFZXOU5fdGxvv4RwghRMaiUqkYUXUEq1us1k/PVW1pNR6/MV647N+/P/v37ydLliycO3eO8uXLc+bMGaP1J8yfBNk07LfffmP9+vVYWlqyYcMGsmfPnmz/sH3DOPP4DJltM7Ox9UYZFyuEEMLg2pdoz6Guh3Czc+NC+AUqLKrA+bDzRuuvZs2anD17luLFixMWFkaNGjVk8YR/oVKp8PDwwMPDQ5aoNUfpeYna48ePU6tWLRITE5k1axZffvllsv2rLq2i05ZOAOzssJNGBRqZokwhhBAZxP3X9/lizRdce34NOys7VrdYTdPCTY3WX1RUFF26dGHr1q0A9OzZk9mzZydbyVKkT7JErZl79OgRrVu3JjExkfbt2zNo0KBk+4PDg+mzow8Ao6uPlhArhBDC6PJmzsvJHiepl78esQmxNF/XnGknpxlthgFHR0c2bdrEpEmTUKvVLF68mCpVqnD//n2j9CfMkwTZNCY2NpZmzZoRHh5O8eLFWbBgQbKPC57HPKfZ2mbEJsRSP399xtYYa8JqhRBCZCTOts7s7LCT/uX6o6Dw7f5v6bO9DwnaBKP0p1arGTlyJHv37iVr1qycP38eHx8fdu/ebZT+hPmRIJuGKIpCz549CQoKImvWrPz11184ODjo9ydoE2izsQ0PIx/i7erNmpZrsFBbmLBiIYQQGY2l2pI5jeYws/5M1Co1iy4sov7K+ryIfWG0Pv38/Dh//jy+vr68fv0af39/xo4di1ZrnCnBzEVCQgIzZ85k5syZJCQY55eJtE6CbBoyefJk1q5di6WlJRs3biRv3rzJ9g/dO5QjD47gYO3AtnbbyJwps4kqFUIIkZGpVCqGVBzCX+3+wsHagcMPDlN+YXkuP71stD49PT0JCAhgwIABKIrChAkT8Pf35+XLl0brM61TFIXIyEgiIyMz7CISEmTTiG3btjFq1CgAZs+eTY0aNZLtX3x+MbPPzgZgVYtVFHUrmuo1CiGEEP/kX9CfUz1PkS9zPh5EPKDS4kpsurbJaP3Z2NgwZ84cVqxYQaZMmdi7dy9ly5bl9OnTRutTpG0SZNOAoKAgOnRIWshg4MCB9O3bN9n+Iw+O0H9nfwAm1JxAk0JNUr1GIYQQ4n2KZyvO2d5n8cvnR0xCDK02tGLs4bHoFJ3R+uzcuTOnT5/G29ubkJAQqlWrxtSpU9HpjNenSJskyJpYaGgojRs3JjY2lrp16zJjxoxk+2++uEmLdS1I0CXQumhrRlUfZaJKhRBCiPdzzeTK7o67+cr3KwAmHJ1Ay/UtiYqPMlqfJUuW5Ny5c7Rt25bExES+++47/P39efbsmdH6FGmPSYPs0aNHady4sX4i37/nisso3rx5g7+/P2FhYRQvXpwNGzZgZWWl3/8i9gX+q/15HfeairkqsrzZctQq+d1DCCFE2mOptmRGgxksbboUawtrtt7YSqXFlbj76q7R+nR2dmbNmjUsWLAAW1tb9uzZQ+nSpTl06JDR+hRpi0lTUUxMDKVKlWLOnDmmLMMkEhISaNOmDZcvX8bd3Z2dO3fi7Oys3x+fGE/zdc25+/oueVzysLXtVjJZZTJhxUIIIcR/61a6GwHdAsjhkIOrz69SYVEF9t/db7T+VCoVvXv35uzZsxQtWpSwsDD8/PwYM2YMiYmJRutXpA0mDbINGzbkxx9/pHnz5qYsI9UpisLAgQPZu3cvdnZ27NixAy8vL/1+naKjx189OB5yHCcbJ3Z22El2h+wfaFEIIYRIOyrmqsi5PueokLMCr96+ov7K+kwMmGjUcbPFixfn7Nmz9OzZE0VRmDhxInXq1OHRo0dG69PUVCoVbm5uuLm5ZdglauVzahMYM2YMCxcuRK1Ws2bNGnx8fJLt/27/d6y+vBoLlQUbW2+UGQqEEEKYHQ9HDwK6BdCrTC8UFMYcGcMXq7/gZazxpsuys7Nj0aJFrF69GgcHB44ePUqpUqXYsGGD0fo0JSsrKwYMGMCAAQOSDU3MSMwqyMbHx/PmzZtkL3Mze/ZsfvzxRwDmzZtHkybJZyD49eSv/HrqVwCWNF1C3fx1U71GIYQQwhBsLW1Z2GQhS5suxdbSlt13dlN2QVnOPj5r1H7bt2/PhQsX8PHx4dWrV7Rp04ZOnToRERFh1H5F6jOrIDt58mScnZ31L09PT1OX9EnWr1/P4MGDAZgwYQJ9+vRJtn/lpZV8s/8bAH7x+4Uupbqkeo1CCCGEoXUr3Y3TPU/j7epNSGQIVZdWZd7ZeUadxN/b25uTJ08yatQo1Go1q1atokSJEhw8eNBofYrUZ1ZBduTIkfoVLCIjIwkNDTV1SR/twIEDdOrUCUVRGDBgAD/88EOy/Xvu7KH7tu4ADK04lG8qf2OKMoUQQgijKOVeinO9z9G8cHM0Wg0Ddg2g85bOxGhijNantbU1P/74I8ePH8fb25tHjx7h5+fHV199xdu3b43Wb2pJSEhg7ty5zJ07V5aoNQc2NjY4OTkle5mDY8eO0bRpUxISEmjVqhWzZs1KNij76MOjtFjXgkRdIh1LdGRqvakZdtC2EEKI9MvZ1plNbTYxre40LFQWrLq8Ct9Fvlx/ft2o/VaqVIng4GD69esHwG+//UbZsmU5d+6cUfs1NkVReP78Oc+fP5clak0hOjqa4OBggoODAbh//z7BwcGEhISYsiyDCgwMxN/fn9jYWOrXr8/KlSuxsLDQ7z/z6Az+q/15m/iWRgUasaTpEpkrVgghRLqlUqkYVnkYh7oewt3BnavPr+KzwIc/zv1h1DBmb2/PvHnz2LVrF+7u7ty4cYNKlSoxfvx4NBqN0foVxmXSxHTu3DnKlClDmTJlABg6dChlypRhzJgxpizLYIKDg6lfvz5RUVHUrFmTzZs3Y2Njo99/IewCDVY1IFoTTe28tdnUZhPWFtYmrFgIIYRIHdVzV+dC3wvUy1+Pt4lv6b+zP83XNedF7Auj9tuwYUOuXLlC69atSUxMZNy4cZQtW5ZTp04ZtV9hHCYNsjVr1kRRlHdey5YtM2VZBnHt2jXq1q1LREQElSpVYvv27djZ2en3X312lXor6xERF0EVzypsa7cNW0tbE1YshBBCpC53B3d2d9zN9HrTsbawZtvNbZScV5ID9w4Ytd8sWbKwbt061qxZg5ubG1evXqVKlSoMHjyYqCjjLasrDE8+wzaCS5cuUbNmTV68eIGPjw+7d+/GwcFBv//y08vUWl6LF7EvKOdRjp0dduJg7fCBFoUQQoj0Sa1S83WlrznT6wxFshYhLDqMun/W5Zt93xCfGG+0flUqFe3ateP69et07doVRVH4/fffKVasGDt27DBav8KwJMga2Pnz56lVqxbPnz+nbNmy7N27N9nSsxfCLlBreS2exz6njHsZ9nTcg7Ot8wdaFEIIIdK/0u6lOdfnHP3L9Qfg11O/UmlxJW68uGHUfrNkycKyZcvYt28fefPmJTQ0lMaNG9OuXTuePn1q1L7F55Mga0CBgYHUqVOHV69eUaFCBQ4ePEiWLFn+b//jQGqvqM3Lty+pkLMCB7scJItdlg+0KIQQQmQcdlZ2zPWfy7Z228hql5UL4RcoO78sM0/PRKvTGrXvunXrcuXKFb799lvUajXr1q2jSJEiLFq0CJ3OeEvrfg6VSqWfWz+jznakUsx4voY3b97g7OxMZGSkyafiOnbsGP7+/kRFRVG5cmV2796drKYTISdouKohUZooKntWZnfH3TjZmMf0YUIIIURqC4sKo+vWruy/tx+ASrkqsbjJYoq4FTF63+fPn6dXr15cuHABAB8fH2bNmkXlypWN3rf4tHwnd2QNYNu2bdSrV4+oqChq1KjB3r17k134Hbd2UPfPukRpoqiRuwZ7O+2VECuEEEJ8QA7HHOzptIf5X8zH0dqRU49OUXp+aX46+hMJWuNO/l+2bFkCAwOZPn06Tk5OBAUFUaVKFTp16sTjx4+N2rf4NBJkP9OiRYto0aIFcXFxfPHFF+zatSvZg11LLyyl2dpmvE18i38Bf3Z13CUPdgkhhBAfQa1S08enD9cGXsO/gD8arYYfDv9A+YXlOR923qh9W1pa8vXXX3P79m169eqFSqVi1apVFCpUiEmTJhEXF2fU/sXHkSCbQoqiMGnSJHr37o1Op6NHjx5s2bJFP8WWoihMOT6FHn/1QKto6VqqK1vabsHOyu4/WhZCCCHEP+VyysX29ttZ2XwlWTJl4eLTi1RYWIGRB0YSl2jcQJktWzYWLlzI2bNnqVy5MjExMYwaNYpixYqxbds2k66olZCQwMKFC1m4cKEsUSs+XkJCAv369WPUqFEAjBw5kkWLFmFpaQlAoi6RL3d/yciDIwH4rvJ3LG26FCsLK5PVLIQQQpgzlUpFx5IduTbwGm2KtUGraJlyYgql/ijFwXsHjd6/j48Px48fZ+XKlXh4eHDv3j2aNWtGrVq1TLaYgqIoPHnyhCdPnsgSteLjvHr1ivr167NgwQJUKhUzZ85k0qRJ+qcFI+IiaLSqEXPOzgHg13q/8nPdnzPs04RCCCGEIWWzz8a6VuvY0nYL7g7u3Hp5C78//Wi7sS2P3jwyat8qlYqOHTty8+ZNvv/+e2xsbAgICKBy5co0adKES5cuGbV/8S4Jsp/g1q1bVKxYkcOHD+Pg4MBff/3FkCFD9PvvvLpDpcWV2H9vP3ZWdmxpu4WhlYaasGIhhBAifWpWuBnXB17nywpfolapWX91PYVnF+aXE78YdSEFAAcHB3766Sf9+FkLCwu2b99O6dKl6dChA3fu3DFq/+L/SJD9SLt378bX15fbt2/j5eXFiRMn+OKLL/T7D90/hO8iX268uEEup1wc736cZoWbma5gIYQQIp1zsXVhVsNZBPUJorJnZWISYhh+YDjF5hZj8/XNRv+43dPTk4ULF3Lt2jXatm2LoiisWbOGwoUL07dvX5nhIBVIkP0PWq2WsWPH4u/vT0REBBUrViQwMJCSJUsCoFN0TDo2ibp/1uXV21dUyFmBwF6BlMlRxsSVCyGEEBlDaffSHOt+jKVNl+Lu4M7d13dpub4lNZfXJOhJkNH7L1iwIGvXruX8+fM0bNgQrVbLggUL8Pb2ZtCgQTx48MDoNWRUEmQ/4MWLFzRq1IgJEyagKAr9+/fnyJEjZM+eHYDXb1/TdG1TRh0ahU7R0aN0D450PUIOxxwmrlwIIYTIWNQqNd1Kd+P2l7f5odoP2FracvThUcotLEe7je24/fK20WsoU6YMu3bt4ujRo1StWpW4uDjmzJmDt7c3nTp14vLly0avIaORIPsvjh49StmyZdm3bx+ZMmXizz//ZO7cudjY2ABJy836LPBhx60d2FjYsKjxIhY3XUwmq0wmrlwIIYTIuBysHZhYeyI3B92kY4mOqFCx7uo6iswpQt/tfXn8xvgf91erVo2jR49y8OBB6tati1arZdWqVZQsWZLGjRtz4sQJg/VlZ2enn/ozI5Ilav+HRqNh7Nix/PzzzyiKQoECBdi0aRMlSpQAkqbWmnRsEhMCJqBVtOTLnI+NrTfKUAIhhBAiDboYfpFRh0ax8/ZOAGwsbOhVthcjqo4gl1OuVKkhKCiIn3/+mY0bN+rH7VatWpURI0bQsGFD1Gq5r/hPn5LvJMj+w40bN+jYsSPnzyetFtKjRw9mzpyJo6MjAPde36PT5k6cepQ0X1zbYm2Z5z+PzJkyf3bfQgghhDCe4yHH+f7g9xwLOQaAtYU1PUr3YHjV4eRxyZMqNdy6dYupU6eyfPly/QIG+fPnp3///nTv3h1XV9dUqSOtkyD7iRITE5k+fTpjx44lLi4OV1dXFi5cSIsWLYCkB7r+OPcHww8MJ1oTjZONE3MbzaVDiQ4yP6wQQghhJhRF4ciDI4wPGE/AwwAALFQWtC3elm8rf0tp99KpUsfjx4+ZOXMmCxcuJDIyEgBbW1vat2/PwIED8fHxSZU60ioJsp8gKCiIXr16ERwcDEC9evVYunQpHh4eAFx7fo3e23tzMvQkANW8qrGi+YpU++1NCCGEEIYX8CCAiUcncvD+/60K5pfPj2GVhlEvfz3UKuN/3B8TE8Pq1auZM2cOFy9e1G+vUKECAwcOpE2bNtja2v7r+xMSEli1ahUAHTt2xMoqfawg+in5LsMOynjz5g3Dhg2jQoUKBAcHkzlzZpYuXcqePXvw8PAgLjGO8UfGU2Z+GU6GnsTB2oHfG/7O4a6HJcQKIYQQZq5Gnhoc6HKAoD5BtCveDrVKzYF7B2i4qiFF5hRh1plZvIl/Y9Qa7O3t6d27NxcuXODEiRN06NABKysrAgMD6dq1Kzlz5mTAgAGcPn36vXPiKorCw4cPefjwYYZdojbD3ZHVarUsW7aM77//nmfPngHQrl07Zs6cSfbs2VEUhc3XN/PN/m94EPEAAP8C/szzn4ens6exTkUIIYQQJnT/9X1mnp7J0uClRGmigKQZEDqW6Eivsr3wyeGTKsMJnz17xuLFi/njjz8ICQnRby9QoABdunShU6dO5MmTB0h6QH3y5MkAjBw5Emtra6PXlxpkaMG/OHr0KF999RUXLlwAkr4pZs6cSaNGjQC49PQSX+35isMPDgOQ0zEnv9b7lTbF2shYWCGEECIDiIqP4s9LfzI7cDbXX1zXby+VvRS9yvaiQ4kOuGYy/kNZWq2WQ4cOsWLFCjZv3kxsbKx+X/Xq1encuTNNmzZl7ty5gARZs/SxJ3r27FlGjx7N3r17AXBycmLs2LEMGjQIa2trbr+8zfiA8ay+vBoFBVtLW76t/C3DqwzH3to+tU5HCCGEEGnE3w+GLbqwiE3XNhGvjQfASm1FwwIN6VSiE18U/CJV5o+Piopiy5YtrFixgkOHDumHEVhZWZE7d26KFCnC3LlzyZUrdaYTMzYJsv/fxYsXGTt2LNu2bQPA0tKSnj17MmHCBLJly8bDiIdMPDqRZcHL0CpaAFoXbc0vdX+RcbBCCCGEAODV21esvryaxRcWExwerN/uZONE44KNaVa4GQ28G+Bg7WD0Wh49esSqVatYsWIF165d029Xq9VUrVqVFi1a0Lx5c7y8vIxei7Fk6CCrKAoHDx5k6tSp7Nu3D0j6y+3UqRNjxowhf/78XHp6iWknp7HmyhoSdYlA0jjYCbUmUDZHWZOdjxBCCCHStqvPrrLq8ipWXV5FSOT/jWG1sbDBL58fzQo3o3HBxmR3yG70Wi5evMjQoUO5fv06YWFhyfb5+PjQpEkT6tWrR/ny5bGwsDB6PYaSIYOslZUVGzZsYMaMGfqptNRqNW3btmXMmDEULFSQQ/cP8eupX9lzZ4++jTp56zCx1kQqeVYy0VkIIYQQwtzoFB2nQk+x9cZWttzYwt3Xd/X7VKio7FmZZoWbUS9/PYpnK26U6bw0Gg3Tpk0DoHXr1uzatYvNmzdz7NixZLMYZM6cGT8/P+rXr0+9evXw9EzbD69nuCDbv39/1q5dy+vXr4GkdYd79uzJ119/jWN2R5YFL2N+0HzuvLoDgFqlplXRVnxb+VvKeZQz5SkIIYQQwswpisK159fYemMrW29u5dyTc8n2u9m5UStvLWrnqU2dfHXInzm/UR8if/bsGdu3b2f37t0cOHBAv+jC34oUKUL9+vWpVasWVapUIUuWLEarJSUyXJD9m5eXF3369KFH7x6ce32OVZdXseXGFjRaDQCO1o50LdWVryt9Tb7M+UxVthBCCCHSsdDIUP66+Rc7bu/g6MOjxCbEJtvv6eRJnXx1qJWnFhVzVaSAawGjBdvExETOnj3L3r172bt3L4GBgeh0umTHFC1alGrVqlG1alWqVauGl5eXSWdrynBBtmHDhvQb0A+7InZsvL6RDdc28OrtK/1xPjl86FeuH+2Kt0uVgdhCCCGEEAAarYbAx4Ecun+Ig/cPcir0FAm6hGTHZLbNTPmc5fHN6UuFnBXwzemLm72bUep59eoVhw4dYt++fRw7dowbN268c0yuXLmoVq0aFSpUwMfHhzJlyuDgkHr5KcMF2TZ/tmHf431ExEXo9+VwyEG74u3oWKIjPh4Ze81iIYQQQqQNsQmxnAg5wcH7Bzn68CgXwi8Qlxj3znF5XfJSzqMcJbKVoET2EpTIVoK8mfMmG2ubmJjI+vXrAWjTpg2WlpafXM/z5885ceIEx48f59ixY5w/f57ExMRkx6hUKgoXLoyPjw/lypXTh1t7e+NMUZrhgiwjANukMSiNCzamQ4kO1MxTEwu1+TyhJ4QQQoiMJ0GbwKWnlwh8HMiZx2cIfByYbCGGf7KzsqOYWzF9uC3gUoD96/bjjDOjRo4yyIIIMTExnDlzhhMnThAUFERQUBCPHj165ziVSkXevHkpVqxYslfhwoXJlOnz5tbNcEF2yJYhtC7bmoq5Kkp4FUIIIYRZi4yL5OyTs1wIu8DlZ5e5/Owy159f1y/K8L/UqPF29aZAlgIUcC2At6s33q7e5Mucj1xOuT570YanT58SFBTEuXPn9OH28ePH7z1WpVKRL18+ihYtire3N/nz59e/8uTJg5WV1X/2Z3ZBds6cOUydOpXw8HBKlSrF77//ToUKFf7zfZ+6RK0QQgghhDlK1CVy59UdLj+9rA+3N57f4PbL22jRfvC9bnZueDl74ensiZeTl/5rD0cP3B3ccXdw/+RniJ49e8a1a9e4evVqstfLly//9T0WFhZ4eXmRP39+8ubNi6enJ7ly5dL/N1euXDg4OJhXkF23bh1dunThjz/+wNfXl5kzZ7JhwwZu3rxJtmzZPvheCbJCCCGEyKg0Gg0/Tf6JN7yhfvv6PHzzkNuvbnPn1R1uv7rNw4iHxCTEfFRb9lb2uDu4k90he1K4tXcni10WXDO5vveV2TYzVhbJ764qisKzZ8+4evUqN27c4O7du8leb9++/c86XFxcyJEjB9evXzePIOvr60v58uWZPXs2ADqdDk9PT7788ktGjBjxwfdKkBVCCCFERqXRaJg8eTIAI0eOfGeMrKIovI57TWhkKCGRIYREhhD65v++Do8OJyw67J3pwT6Wo7WjPti62LrgaOOIg7UDjtZJ//3n1/bW9iTGJBLxLIJXYa94Ef6CF2EveBb2jLDQMMJCw4h6HQX/mBnsY/Ldpz/eZkAajYagoCBGjhyp36ZWq/Hz8+PUqVMmrEwIIYQQwrypVCp90CzlXupfj4vWRBMeHZ7s9TT6Ka/evuJV3Kuk//7j9fcsUVGaKKI0UTyMfJiC4gCP///y+edmFVYaKzSTNB/VjEmD7IsXL9BqtWTPnnw94uzZs793XrP4+Hji4/9voPPfK1W8efPGuIUKIYQQQqQxGo2GuLikqbvevHnzWbMWZLPMRjaXbJR0Kfmfx2p1WiLjInkd95rXb18n/TfuNTGaGGISYoiKj9J/HZ0QTXR8dNLXmmiiNUlfx2hi0Gg1xCfGk6j7v+m+FBQ08Ukh9mMGDZg0yH6qyZMnM378+He2p/U1g4UQQgghjGnKlCmmLsHgoqKikq3g+j4mDbJZs2bFwsKCp0+fJtv+9OlT3N3d3zl+5MiRDB06VP9nnU7Hq1evyJIli0mXUkstb968wdPTk9DQUBkT/AFynT6OXKePI9fp48m1+jhynT6OXKePkx6vk6IoREVF4eHh8Z/HmjTIWltb4+Pjw8GDB2nWrBmQFE4PHjzIoEGD3jnexsYGGxubZNtcXFxSodK0xcnJKd18sxqTXKePI9fp48h1+nhyrT6OXKePI9fp46S36/Rfd2L/ZvKhBUOHDqVr166UK1eOChUqMHPmTGJiYujevbupSxNCCCGEEGmYyYNs27Ztef78OWPGjCE8PJzSpUuzZ8+edx4AE0IIIYQQ4p9MHmQBBg0a9N6hBCI5Gxsbxo4d+87wCpGcXKePI9fp48h1+nhyrT6OXKePI9fp42T062TyBRGEEEIIIYRICbWpCxBCCCGEECIlJMgKIYQQQgizJEFWCCGEEEKYJQmyacyrV6/o2LEjTk5OuLi40LNnT6Kjoz/4ngULFlCzZk2cnJxQqVRERES8c0yePHlQqVTJXua8CoixrlNK2k3LUnI+cXFxDBw4kCxZsuDg4EDLli3fWbTkf7+XVCoVa9euNeapGNScOXPIkycPtra2+Pr6EhgY+MHjN2zYQOHChbG1taVEiRLs2rUr2X5FURgzZgw5cuQgU6ZM+Pn5cfv2bWOeQqow9HXq1q3bO983DRo0MOYppIpPuU5Xr16lZcuW+p/JM2fO/Ow2zYWhr9O4cePe+X4qXLiwEc8gdXzKdVq4cCHVqlUjc+bMZM6cGT8/v3eOT68/n/QUkaY0aNBAKVWqlHL69Gnl2LFjire3t9K+ffsPvmfGjBnK5MmTlcmTJyuA8vr163eOyZ07tzJhwgQlLCxM/4qOjjbSWRifsa5TStpNy1JyPv369VM8PT2VgwcPKufOnVMqVqyoVK5cOdkxgLJ06dJk309v37415qkYzNq1axVra2tlyZIlytWrV5XevXsrLi4uytOnT997/IkTJxQLCwvll19+Ua5du6b88MMPipWVlXL58mX9MVOmTFGcnZ2VrVu3KhcvXlSaNGmi5M2b12yuyfsY4zp17dpVadCgQbLvm1evXqXWKRnFp16nwMBA5ZtvvlHWrFmjuLu7KzNmzPjsNs2BMa7T2LFjlWLFiiX7fnr+/LmRz8S4PvU6dejQQZkzZ45y4cIF5fr160q3bt0UZ2dn5dGjR/pj0uPPp3+SIJuGXLt2TQGUs2fP6rft3r1bUalUyuPHj//z/YcPH/5gkH3fDwJzZKzr9LntpjUpOZ+IiAjFyspK2bBhg37b9evXFUA5deqUfhugbNmyxWi1G1OFChWUgQMH6v+s1WoVDw8PZfLkye89vk2bNoq/v3+ybb6+vkrfvn0VRVEUnU6nuLu7K1OnTtXvj4iIUGxsbJQ1a9YY4QxSh6Gvk6IkBdmmTZsapV5T+dTr9E//9nP5c9pMq4xxncaOHauUKlXKgFWa3uf+3ScmJiqOjo7K8uXLFUVJvz+f/kmGFqQhp06dwsXFhXLlyum3+fn5oVarOXPmzGe3P2XKFLJkyUKZMmWYOnUqiYmJn92mKRjrOhn7+qe2lJxPUFAQCQkJ+Pn56bcVLlwYLy8vTp06lezYgQMHkjVrVipUqMCSJUtQzGAmP41GQ1BQULLzU6vV+Pn5vXN+fzt16lSy4wHq16+vP/7+/fuEh4cnO8bZ2RlfX99/bTOtM8Z1+tuRI0fIli0bhQoVon///rx8+dLwJ5BKUnKdTNGmqRnznG7fvo2Hhwf58uWjY8eOhISEfG65JmOI6xQbG0tCQgKurq5A+vz59L/SxIIIIkl4eDjZsmVLts3S0hJXV1fCw8M/q+3BgwdTtmxZXF1dOXnyJCNHjiQsLIzp06d/VrumYKzrZMzrbwopOZ/w8HCsra1xcXFJtj179uzJ3jNhwgRq166NnZ0d+/btY8CAAURHRzN48GCDn4chvXjxAq1W+87KgdmzZ+fGjRvvfU94ePh7j//7evz93w8dY26McZ0AGjRoQIsWLcibNy93797l+++/p2HDhpw6dQoLCwvDn4iRpeQ6maJNUzPWOfn6+rJs2TIKFSpEWFgY48ePp1q1aly5cgVHR8fPLTvVGeI6DR8+HA8PD31wTY8/n/6XBNlUMGLECH7++ecPHnP9+nWj1jB06FD91yVLlsTa2pq+ffsyefLkNLMaSFq4TuYgLVyn0aNH678uU6YMMTExTJ06Nc0HWWFa7dq1039dokQJSpYsSf78+Tly5Ah16tQxYWXCHDVs2FD/dcmSJfH19SV37tysX7+enj17mrAy05gyZQpr167lyJEj2NramrqcVCNBNhUMGzaMbt26ffCYfPny4e7uzrNnz5JtT0xM5NWrV7i7uxu0Jl9fXxITE3nw4AGFChUyaNspZerrlJrX/3MY8zq5u7uj0WiIiIhIdlf26dOnH7wGvr6+TJw4kfj4+DTzi9H7ZM2aFQsLi3dmYfjQ+bm7u3/w+L//+/TpU3LkyJHsmNKlSxuw+tRjjOv0Pvny5SNr1qzcuXPHLINsSq6TKdo0tdQ6JxcXFwoWLMidO3cM1mZq+pzrNG3aNKZMmcKBAwcoWbKkfnt6/Pn0v2SMbCpwc3OjcOHCH3xZW1tTqVIlIiIiCAoK0r/30KFD6HQ6fH19DVpTcHAwarX6nY+eTcnU1yk1r//nMOZ18vHxwcrKioMHD+q33bx5k5CQECpVqvSvNQUHB5M5c+Y0HWIBrK2t8fHxSXZ+Op2OgwcP/uv5VapUKdnxAPv379cfnzdvXtzd3ZMd8+bNG86cOfPBa5aWGeM6vc+jR494+fJlsn9gzUlKrpMp2jS11Dqn6Oho7t69m+G+n3755RcmTpzInj17kj0TAenz59M7TP20mUiuQYMGSpkyZZQzZ84ox48fVwoUKJBsuqRHjx4phQoVUs6cOaPfFhYWply4cEFZuHChAihHjx5VLly4oLx8+VJRFEU5efKkMmPGDCU4OFi5e/eusnLlSsXNzU3p0qVLqp+foRjjOn1Mu+YmJdepX79+ipeXl3Lo0CHl3LlzSqVKlZRKlSrp9//111/KwoULlcuXLyu3b99W5s6dq9jZ2SljxoxJ1XNLqbVr1yo2NjbKsmXLlGvXril9+vRRXFxclPDwcEVRFKVz587KiBEj9MefOHFCsbS0VKZNm6Zcv35dGTt27Hun33JxcVG2bdumXLp0SWnatKnZT29j6OsUFRWlfPPNN8qpU6eU+/fvKwcOHFDKli2rFChQQImLizPJORrCp16n+Ph45cKFC8qFCxeUHDlyKN98841y4cIF5fbt2x/dpjkyxnUaNmyYcuTIEeX+/fvKiRMnFD8/PyVr1qzKs2fPUv38DOVTr9OUKVMUa2trZePGjcmmIYuKikp2THr7+fRPEmTTmJcvXyrt27dXHBwcFCcnJ6V79+7JviHv37+vAMrhw4f128aOHasA77yWLl2qKIqiBAUFKb6+voqzs7Nia2urFClSRJk0aZJZ/+NhjOv0Me2am5Rcp7dv3yoDBgxQMmfOrNjZ2SnNmzdXwsLC9Pt3796tlC5dWnFwcFDs7e2VUqVKKX/88Yei1WpT89Q+y++//654eXkp1tbWSoUKFZTTp0/r99WoUUPp2rVrsuPXr1+vFCxYULG2tlaKFSum7Ny5M9l+nU6njB49WsmePbtiY2Oj1KlTR7l582ZqnIpRGfI6xcbGKvXq1VPc3NwUKysrJXfu3Erv3r3NOpz97VOu09//z/3vq0aNGh/dprky9HVq27atkiNHDsXa2lrJmTOn0rZtW+XOnTupeEbG8SnXKXfu3O+9TmPHjtUfk15/Pv1NpShmMGeOEEIIIYQQ/0PGyAohhBBCCLMkQVYIIYQQQpglCbJCCCGEEMIsSZAVQgghhBBmSYKsEEIIIYQwSxJkhRBCCCGEWZIgK4QQQgghzJIEWSGEEEIIYZYkyAohRCqrWbMmX331lanLEEIIsydBVgghzMyyZctwcXExdRlCCGFyEmSFEEIIIYRZkiArhBAmkJiYyKBBg3B2diZr1qyMHj0aRVEAiI+P55tvviFnzpzY29vj6+vLkSNHADhy5Ajdu3cnMjISlUqFSqVi3LhxAPz555+UK1cOR0dH3N3d6dChA8+ePTPRGQohhPFJkBVCCBNYvnw5lpaWBAYG8ttvvzF9+nQWLVoEwKBBgzh16hRr167l0qVLtG7dmgYNGnD79m0qV67MzJkzcXJyIiwsjLCwML755hsAEhISmDhxIhcvXmTr1q08ePCAbt26mfAshRDCuFTK37cAhBBCpIqaNWvy7Nkzrl69ikqlAmDEiBH89ddf7Nmzh3z58hESEoKHh4f+PX5+flSoUIFJkyaxbNkyvvrqKyIiIj7Yz7lz5yhfvjxRUVE4ODgY85SEEMIk5I6sEEKYQMWKFfUhFqBSpUrcvn2by5cvo9VqKViwIA4ODvpXQEAAd+/e/WCbQUFBNG7cGC8vLxwdHalRowYAISEhRj0XIYQwFUtTFyCEEOL/REdHY2FhQVBQEBYWFsn2feiuakxMDPXr16d+/fqsWrUKNzc3QkJCqF+/PhqNxthlCyGESUiQFUIIEzhz5kyyP58+fZoCBQpQpkwZtFotz549o1q1au99r7W1NVqtNtm2Gzdu8PLlS6ZMmYKnpyeQNLRACCHSMxlaIIQQJhASEsLQoUO5efMma9as4ffff2fIkCEULFiQjh070qVLFzZv3sz9+/cJDAxk8uTJ7Ny5E4A8efIQHR3NwYMHefHiBbGxsXh5eWFtbc3vv//OvXv3+Ouvv5g4caKJz1IIIYxLgqwQQphAly5dePv2LRUqVGDgwIEMGTKEPn36ALB06VK6dOnCsGHDKFSoEM2aNePs2bN4eXkBULlyZfr160fbtm1xc3Pjl19+wc3NjWXLlrFhwwaKFi3KlClTmDZtmilPUQghjE5mLRBCCCGEEGZJ7sgKIYQQQgizJEFWCCGEEEKYJQmyQgghhBDCLEmQFUIIIYQQZkmCrBBCCCGEMEsSZIUQQgghhFmSICuEEEIIIcySBFkhhBBCCGGWJMgKIYQQQgizJEFWCCGEEEKYJQmyQgghhBDCLEmQFUIIIYQQZun/ASnsCvQzJIoeAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "mverse.plot_integration(measure=\"beta\", weights=weights, true_value=0.2)" ] } ], "metadata": { "kernelspec": { "display_name": "comet", "language": "python", "name": "python3" }, "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.13.0" } }, "nbformat": 4, "nbformat_minor": 5 }