Skip to content

Instantly share code, notes, and snippets.

@akelleh
Created May 5, 2019 00:53
Show Gist options
  • Save akelleh/dd51c4a9c03ee96f6ca98bb1c785c611 to your computer and use it in GitHub Desktop.
Save akelleh/dd51c4a9c03ee96f6ca98bb1c785c611 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [],
"source": [
"%matplotlib inline"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"import pandas as pd\n",
"import dowhy.api\n",
"import seaborn as sns\n",
"import matplotlib.pyplot as pp"
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {},
"outputs": [],
"source": [
"sns.set(rc={'figure.figsize':(7.5, 5)})\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Generating Data"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Here, there is no causal relationship between title length and click-through rate. For each population, the click-through rate is independent of title length. "
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {},
"outputs": [],
"source": [
"N = 10000\n",
"n = 1000\n",
"\n",
"# normally, ctr is independent of title_length\n",
"title_length = np.random.choice(range(25), size=N)\n",
"click_through_rate = np.random.beta(5, 100, size=N)\n",
"\n",
"# but one quirky person prefers titles in a much narrower range, and is especially talented\n",
"title_length_2 = np.random.normal(12, 3 ,size=n).astype(int)\n",
"click_through_rate_2 = np.random.beta(10, 100, size=n)\n",
"\n",
"# let's concatenate these together\n",
"all_title_lengths = np.array(list(title_length) + list(title_length_2))\n",
"all_click_rates = np.array(list(click_through_rate) + list(click_through_rate_2))\n",
"\n",
"df = pd.DataFrame({'click_through_rate': all_click_rates, 'title_length': all_title_lengths, 'author': [1]*N + [0]*n})"
]
},
{
"cell_type": "code",
"execution_count": 14,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x7fe33572b2e8>"
]
},
"execution_count": 14,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAeoAAAFFCAYAAAAjJ64BAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzs3XmUXXWZ6P3vHs58as6pSlVlJiEhmAFIQ6RDIEpMEwNBAjYX20j3C1zXUrGzLnjpi9ymaRu9rig40Poq16ZFWA1OibwqoCiKBoEEhGAmQqZKpeaqM097n73fPw5VJCSp2mesquT5rMWi6uwz/Gpnn/3s/RueR7Ft20YIIYQQE5I63g0QQgghxOlJoBZCCCEmMAnUQgghxAQmgVoIIYSYwCRQCyGEEBOYBGohhBBiApNALYQQQkxgEqiFEEKICUwCtRBCCDGBSaAWQgghJjAJ1EIIIcQEJoFaCCGEmMD08W5AtQ0NJbCs0uuQNDUFGRiIl6FFZx/Zd8WTfVc82Xelkf1XvOF9p6oKDQ2Bgl9/1gVqy7LLEqiH30sUR/Zd8WTfFU/2XWlk/xWvlH0nXd9CCCHEBCaBWgghhJjAJFALIYQQE5gEaiGEEGICk0AthBBCTGASqIUQQogJTAK1EEIIMYFJoBZCCCEmMAnUQgghxAQmgVqIs0jOssa7CUKIAkmgFuIskUwbvHU0QsbIjXdThBAFkEAtxFkgk81xsCtGNpujP5walzZYto1hykWCEIWSQC3EGc7MWRzqjqGpUBNwMRBNk86aVW/HYDTNW0cjEqyFKJAEaiHOYJZlc7g7Ri5n4fXoKIqCrqn0DFX3rtowLXoGU9gWHOmJSxUmIQoggVqIM5Rt2xwbiJNMG/h971a09Xk0IrEMqUz17qr7Iyls2ybg10mkDXqGklX7bCEmOwnUQpyh+sIpBqIZgn7XCY8rioLLpdI9WJ1gmcnm6AunCLxzsVDjd9E7lCISz1Tl84WY7CRQC3EGGoql6RpIUut3oSjKSdt9Hp1oIksibVS8Ld2DSXRNHWmHoigEfDodvXEyWRmvFmIsFQ3UTz31FGvXrmX16tU89thjJ23fvXs3GzZsYM2aNdx9992YZr4rrre3l9tuu41rr72WG2+8kaNHjwIQjUa57bbbuOqqq/jYxz5GX19fJZsvxKSUSBt09MYJ+vVTBulhXrdG90AS267ceHEibRCOZ/B79RMe1zUVXVM40huTtd1CjKFigbqnp4cHHniAxx9/nK1bt/LEE0+wf//+E55z5513cs899/DMM89g2zZPPvkkAJ/73OdYtWoVW7ZsYf369WzevBmABx98kGXLlvHLX/6SG264gX/7t3+rVPOFmJQy2RyHumL4PDqaOvrX2+PWSKQN4qnK3FXbtk1XfwKvRzvldq9HJ5PN0TWQqOjFghCTXcUC9bZt21i+fDn19fX4/X7WrFnD008/PbK9s7OTdDrN0qVLAbjuuut4+umnGRwcZM+ePdx4440AbNiwgX/8x38E4Pnnn+fqq68GYN26dfz+97/HMCrfdSfEZGCYFge7o2gauHRnX22vW6N7sDJ31ZFElkTGxOM6daAGCPh0BqIZhmIyXi3E6ehjP6U4vb29hEKhkd+bm5t54403Trs9FArR09NDR0cHbW1t3H///bz00ku0tbVxzz33nPQaXdcJBoMMDg7S0tLiuF1NTcFS/7Tj2lxTtvc628i+K96p9l3OstnfMURNjZegz13Q+4XjaVxeNw213nI1kVzO4lg4TXtL3ZgXDXW1FrFUlmlBLwGfa9TnlkqOu9LI/iteKfuuYoH6VFfox4+XnW67aZrs2rWLz3zmM9x999388Ic/5K677uLRRx895eeoY3TvvdfAQHnWcIZCNfT1xUp+n7OR7LvinWrf2bZNZ1+coXiWGr+LoUxhvUyGafHmW72cO70edZQx7UL0h1P0DySoDTi7aMgaOXb85RjzptWja5Xp6JPjrjSy/4o3vO9UVSnqZrFiXd8tLS309/eP/N7b20tzc/Npt/f19dHc3EwoFCIQCLBq1Sog38U9fCfe3Nw88hrTNInH49TX11fqTxBiUugNJxmMZgj6irvudukqhmGVbbmUYeboHkwSLODu2O3SsCw42hvHkvFqIU5QsUB96aWX8uKLLzI4OEgqleLZZ59l5cqVI9vb29vxeDzs2LEDgC1btrBy5UpmzJhBS0sLv/vd7wD47W9/y/nnnw/A5ZdfzpYtWwD4xS9+wbJly3C5KttVJsRENhhL0z2QoiZw6mVYTvm8+Rng5ZiB3RdOoygKqlpYewK+/JKx/sj45CIXYqKq6B31pk2b2LhxI9deey3r1q1j8eLF3HrrrezcuROAzZs388UvfpGrrrqKVCrFxo0bAfjmN7/Jww8/zLp16/j+97/P/fffD8BnP/tZ/vznP/PhD3+Yxx9/nP/9v/93pZovxIQXTxkcdbAMywldUzFyFuF4tqT3SWVM+iMp/N7TTyAbTdDvors/WbGZ6EJMRop9lq2LkDHq8Sf7rnjD+y6dNdl/NILHrTme4T2WXM4ibVjMn178OPHBrijprInPU/z0F8O0yBoW86bV4R5lxnih5Lgrjey/4k3YMWohRGUYZr4alktXyxakATRNxbJsBmPpol4fTxnEktmSgjTkx8xVFTp6pXiHECCBWohJJZezONITw7ZsPO7y3W0OC3h1egdTGGZhY9WWbXOsP4G3TG3yefLFO6qVj1yIiUwCtRCTSEdvjGTGxOetzMrK4QlgA9HCJnSF4xnSGbOsXdU1fhd94RRhKd4hznISqIWYJAwzR384XfQyLKf8Xp2+oTRZw1nBDDNn0T2QPKGUZjkoikLQp9PREyedrV5JTiEmGgnUQkwSkUQWRaHkGd5jUdX80qr+iLOx6sFompxlVyRRiaapuFwKh7vjmDkp3iHOThKohZgEbNumP5wi4KlO3gC/V6M/khqzDGXWyNEzmCJQoa54AK9bxzBzdPVL8Q5xdpJALcQkkMqYZE0LvYyzvEejKAqaptAbHn0yV284iaYWntykUAGfzlAsI+PV4qwkgVqISWAolqlYDuzT8XvywTGVOfX4cCpjMhjN4CsyuUkhFEXB59XpCxe3dEyIyUwCtRATnJmzGIxl8J2mrnOlKIqCS1fpHTp5Brht23QNJHDrasXHzIe5dJW0YcrEMnHWkUAtxAQXT2axbbtqAfF4Po9OOJEhmT4xOMZTBrGUgbfE5CaF0lCJJktLcyrEZCOBWogJrj+SLlsikWJ4dJXuwcTI75Zl09mfwF/lO3wAr0dlIJKWSWXirCKBWogJLJUxSWbLm0ikUF6PTjxljBTKCMczZA0Ll179NmmaimFapDLO1ngLcSaQQC3EBBZJZNCU8f+aetwa3QMJDNOiayBJoAoTyE5H0xSiCZn9Lc4e438GEEKckmXZDETSVZ9Edioel0Yyk6OzP45t22hVnoF+PJ9bZzCaloId4qwhgVqICSqeNshZdsXXKDvl86gMRtP4K5jcxAlVVTAtm+Rplo0JcaaRQC3EBDUQTuN2TZyvqEvXaKz1TogLB5euMlRkOU4hJpuJcxYQQozIGjliqSyecZxENpF53RrheJacJfm/xZlPArUQE1A0kUVVlHFZOz0ZKIoCNiRS0v0tznwSqIWYYGzbpi+SwjsBJpFNZG6XykBUur/FmU8CtRATTDJjYphW1XN7TzZul0osmcUwpftbnNnkTCDEBDMYTeOqUpWsyUxRFBQgnpKUouLMJmcDISYQM2cRiWfHNWXoZOJxawxK97c4w0mgFmICiSWzWBUqwPGnv3Sz7c3usr/veHK7NBKZHBlDUoqKM9f4Zi4QQpygv0KZyPYeGeLZV44CMKe1lqlN/rJ/xnhRyV/geOp8490UISqionfUTz31FGvXrmX16tU89thjJ23fvXs3GzZsYM2aNdx9992YZn6pxZYtW1ixYgXr169n/fr1PPDAAwAcPXqUj33sY6xfv56Pf/zjdHZ2VrL5QlRVKmOSyphlL3YRiWfY+odDTG304fNoPLu944yqPuV1awxEMmfU3yTE8SoWqHt6enjggQd4/PHH2bp1K0888QT79+8/4Tl33nkn99xzD8888wy2bfPkk08CsHPnTu666y62bt3K1q1b2bRpEwBf+9rX+PCHP8zWrVv50Ic+NBLAhTgTROIZtDJn/bIsm5/8/iCWbXP9Fedw+dI2DnXFeOtopKyfM550XSVj5EhnpftbnJkqFqi3bdvG8uXLqa+vx+/3s2bNGp5++umR7Z2dnaTTaZYuXQrAddddN7J9586dbNmyhWuuuYY77riDSCR/UrEsi3g8DkAqlcLr9Vaq+UJUlWXZDETT+DzlHY16/s+ddPTGWff+mTTWerlofoimWg+/3n70jMrqpar5SmNCnIkqNkbd29tLKBQa+b25uZk33njjtNtDoRA9PT0jP992220sXryYr371q9x333185Stf4bOf/Sw33ngjjz76KIZh8MQTTxTcrqamYAl/1YlCoZqyvdfZRvbdiSLxDMEaH/VBz5jPbagPOHrPvUeG+OMb3Vxy/lQuu3DGyOPrL5/L9576C3s6oqxY0l50mwth2zYv/aWbudPrmVKBseTaGotU1qSpKThqLnI57koj+694pey7igXqU40XHT+TdbTtDz300Mhjt9xyC1deeSUA//N//k/uu+8+rrzySp555hk+/elP87Of/aygGbIDA/GylMcLhWro64uV/D5nI9l3Jzt4LErWzDFkjp4Ss6E+wFA4Meb7xVMGj/5iF011Xj5wQesJr2lv9DJzag2/3HaIc1qDeN2Vn1P66r4+/r9th2mo8XDLuvPK3nMA+bSrRzp1Al7XKbfLcVca2X/FG953qqoUdbNYsa7vlpYW+vv7R37v7e2lubn5tNv7+vpobm4mFovxyCOPjDxu2za6rjM4OMiBAwdGgvaaNWvo6+tjaGioUn+CEFWRGS7AUaa107Zts+WFg2SMHNdfcc5Jk9MUReFDfzWNZMbkD29UfrnWUCzDsy930NLgI5LI8uPfHahILWldU4nEpftbnHkqFqgvvfRSXnzxRQYHB0mlUjz77LOsXLlyZHt7ezsej4cdO3YA+ZneK1euxO/38/DDD/P6668D8IMf/IDVq1fT0NCAx+Nh+/btAOzYsYNAIEBjY2Ol/gQhqmK4AEe5/HFnNweORVlz8QyaG07dzdzaFGDxOU28tKuHcKxywc2ybLb+4SCKonDjB+eydvkMDhyL8ptXy79iw+vWGIxmKnIRIMR4qlifV0tLC5s2bWLjxo0YhsH111/P4sWLufXWW7n99ttZtGgRmzdv5vOf/zyJRIKFCxeyceNGNE3jwQcf5N577yWdTjNr1iy+/OUvoygK3/zmN/nXf/1X0uk0gUCAb3zjG5VqvhBVYdk2/eFU2dZOd/TG+e1rnSyc1cCF504Z9bkfuLCdXYeGeO7VTjZcPqcsn/9eL+3q4UhPnPUrZlEX9HDhuSG6B5Jse7ObqY0+3jenqWyfpaoKlm2TSBvU+N1le18hxptin2WLD2WMevzJvntXPGVw4FiU2sCpx1Xfa7Qx6lTG5Ds/24WqKtx69XmOxp5/+1onL7zexT+sXcC05vJNtAToHUrx3ad2MXdaHR9ddc7IXJJczuL7z+yjayDJ369dQGsZk6+kMyZ+r4sZLSdP3JHjrjSy/4o3YceohRBjyxfgKL3b27ZtfvbHQ8RSBtddPsfxBLG/ft9Ugj4Xz75S3iQouZzFlhcO4nFrrHv/zBMmfGqayg2rzsHv0XjyN/tJpI2yfa7HrRFNZDFzZ87SMyEkUAsxTspZgOOVPX3sPRLmgxe10z7F2fItyOfKXnVBG0f7Euw+XL6Jmb9/vYvuwSTr3j+TgO/k3oKgz8VHPzCXeMrgR88fKNuabkXJd3/HU+UL/kKMNwnUQoyTaCILlF6Ao2sgya9e6WDetDqWL2wp+PVL5k6hpcHHr7cfLcud6NG+OH/Y2cWSc5pYMLPhtM9rmxLg6r+exeHuGL96Jw95OUhFLXGmkUAtxDiwbZv+SBpvieuJM0aOH//ubfxenfUrZhUV9FVV4cq/mkY4nuXl3b0ltccwc2x94SA1fjdrLpk+5vMXn9PEJQtbeHl3L39+q3/M5zvhcWnEUwZZqaglzhASqIUYB+lsPje1Sy/+K2jbNr948TBDsQwfWTkH/2kSfThxTlsdc9treeH1LpIljBk/t6OTgWiG9StmOR4nX71sGrNaa/j5i4fp7IsX/dnHU1CISfe3OENIoBZiHAzFMmglDk2/vn+AnQcGWbmkjVlTS0/tuPqvppM1c/z+9a6iXn/gWJSXd/dyyXnNzG6tdfw6VVW4/vI51PhdPPnbt4knSw+wXo/GQES6v8WZQQK1EFWWsywGo2l8JaTu7Aun+OVLR5g5tYbLFreWpV2heh8Xnhti+54++gsMcumMyc/+cJApdV4+cNG0gj/b781PLktnczz5/H5yJY6Vu3SVdNYknR09JasQk4EEaiGqLJEysW1GLR4xGsO0+PHvDuDSVK5bObvo9zmVK5a2oesKv95e2OSup1/uIJYyWH/Z7KK786c2+rnmr2dxtDfBL186UtR7HE9TVKLJbMnvI8R4k0AtRJX1R1K43cV/9X71Sge9QynWXza77Bm4Aj4XKxa1sq8jzMGuqKPX7D48xBtvD3DZ4taCloadyvmzG/nrRVN5dV8/O/b2lfReXo/KQCRd1vXhQowHCdRCVFHGyBFPm3hcxQ1Qv/5WH9v39vH+81uYN62uzK3LW76whbqAm1+9cnTMIBdPGfx822Fam/xctqQ8XfCrLmhnbnstv3zpCEd6is+EpWkqhmmRysjsbzG5SaAWoooiiUzRX7qhWIb/+tVe2qYE+MCFlasjresqH7ione7BJG+8PXDa59m2zf+37TAZI8e1l81GU8tzOlFVhetWzqE+6OaHv337nfXmxdE0hWhCKmqJyU0CtRBVYts2g5E03iILcPxqewc2sOHyOWhaZb+675vdSNuUAL95tfO065Ff3z/Avo58NrRQ/amrdBXL69H52w/MxchZPPmb/ZhmcZPLfG6dwWhaKmqJSU0CtRBVks7mMEwbvYggOxhNs+dwmMuWtNFQ46lA6040XLM6ljR48S89J20PxzM8/fIRZrYEuaSIbGhOhOp9fOSy2RwbSPLzFw8XNdasqgqmZZPMyOxvMXlJoBaiSmKp4utOv7y7F1VVWLGkcl3e7zWjpYbzZjaw7c1uYsfNnrZtm61/OATA+hWzS06BOpr5Mxq4fGkbr789UHTWNJeuMhSTNdVi8pJALUSVDEWzeIqY7Z3OmPz5rX7eN7uRumDl76aP98GL2slZNs+/dmzksZd29XK4O8aai2dQX4W7+5VLWpk/o55nX+ngWP+pS3yOxuvWCMezJa/NFmK8SKAWogrSWZOMaaIXscb41bf6yZpWxbqYR9NY6+Xi85p57a1+ugeT9IVTPLfjKOdOr2Pp3KaqtEFRFNavmIVLU4tasqUoCtgQK0PGMyHGgwRqIaognjRQKbyLOGdZvLy7l5lTa2ht8legZWO7bHErPo/Gs690sPWFg3hcGusuLa4ASLG8bp2Fsxv5y8HBooptuF0qfeFkBVomROVJoBaiCgZjGTxF1J3efShMNJFl+fnVv5se5vPorFzSxqGuGMcGknz4/TMJnqLGdKUtmdtE1rTYcyRc8GvdLpVoPItR5OxxIcaTBGohKiyTzZHJmgWn1rRtmz/t6qax1sO5FUpu4tSy+SHapvhZNj/EebNOX2O6kmY0B2mo8fDn/YWXw1QUBUWBSFzWVIvJRwK1EBUWTxtQRLd3R2+cY/1JLlnYUtVu5lPRNJX/58Pnsfb9M8etDYqisGRuE4e6YoRjhQfcgM9N10BS6lSLSUcCtRAVNhhN4/UU/lX7064evG6NJedUZ9LWWMb7YgEY2Revj5Ix7XQ0VUHTFLoHZaxaTC4SqIWooIyRI5XN4dILG58eTnBy0fwQ7iLzgp+J6oIeZrfW8Pr+/qISoPg8GkOxDPGUzAAXk4cEaiEqKJEyUIrIXjmc4OTi85rL36hJbuncKYTjWQ73xAt+raIo+DwanX0JSSsqJo2KBuqnnnqKtWvXsnr1ah577LGTtu/evZsNGzawZs0a7r77bkwzn+Zvy5YtrFixgvXr17N+/XoeeOABAHp7e7ntttu49tprufHGGzl6tLCauUJU21AsXXCSk+MTnJS7jOWZYMHMetwuldeLmFQG4HZpZIwcg5KtTEwSFQvUPT09PPDAAzz++ONs3bqVJ554gv3795/wnDvvvJN77rmHZ555Btu2efLJJwHYuXMnd911F1u3bmXr1q1s2rQJgM997nOsWrWKLVu2sH79ejZv3lyp5gtRMsPMkUibBXddj2eCk8nApWucP6uRXYeGip4YFvBpdA8kMUyZWCYmvooF6m3btrF8+XLq6+vx+/2sWbOGp59+emR7Z2cn6XSapUuXAnDdddeNbN+5cydbtmzhmmuu4Y477iASiTA4OMiePXu48cYbAdiwYQP/+I//WKnmC1GyeMooeK73REhwMhksnTcFw7TYdWioqNdrqoqqKPQMycQyMfFVLFD39vYSCoVGfm9ubqanp+e020Oh0Mj2UCjEZz7zGbZu3Uprayv33XcfHR0dtLW1cf/993PNNddw++2343JVP+mCEE4NRgtPcjIREpxMBtNCARprPUV3fwP4vBoD0QzJtEwsExObXqk3PtWMzOOXd4y2/aGHHhp57JZbbuHKK6/kpptuYteuXXzmM5/h7rvv5oc//CF33XUXjz76aEHtamoKFvT80YRCNWV7r7PNmb7vDDOHqz/JlAKKaNi2zfa9ewnV+7j4fW2nrbTVUB8oVzMnteXva+UX2w6RU1Sm1Dmrh/3efecPmCQMm+ntQVR1/JefTXRn+ve2kkrZd44DdXd3N3v37mXFihX09vbS2to66vNbWlrYvn37yO+9vb00NzefsL2//92r4b6+Ppqbm4nFYvz4xz/m5ptvBvInL13XCYVCBAIBVq1aBcC6dev4whe+4LT5IwYG4mWZ7RkK1dDXFyv5fc5GZ8O+C8fShCNJLNN5HeQjPTGO9MS4avkMIpFTd8k21AcYChdeQepMNK+9FoAXXu3gigvGLv95un0XSxhotkVjrbfsbTyTnA3f20oZ3neqqhR1s+io6/v555/nxhtv5F/+5V8YGBhg7dq1/PrXvx71NZdeeikvvvgig4ODpFIpnn32WVauXDmyvb29HY/Hw44dO4D8TO+VK1fi9/t5+OGHef311wH4wQ9+wOrVq5kxYwYtLS387ne/A+C3v/0t559/fsF/sBDVMBTLFtztXekEJ5ZlMxTLkLPOjHzXdQE3c9pqeX3/QFFrqof5vRpdA0nJAy4mLEeB+pvf/CZPPvkktbW1NDc38/jjj/P1r3991Ne0tLSwadMmNm7cyLXXXsu6detYvHgxt956Kzt37gRg8+bNfPGLX+Sqq64ilUqxceNGNE3jwQcf5N577+Wqq67iL3/5C3feeedIOx5++GHWrVvH97//fe6///4S/3whys/MWcRTBu4CcntXI8FJLGXQWOMhmTpzZjovndtEJJHlUHfxd3qapoINvTKxTExQjrq+Lcs6odv6vPPOc5RO8Oqrr+bqq68+4bHvfve7Iz8vWLCAH/3oRye9btmyZfz0pz896fE5c+YUPCYtRLUl0iYohaXcrHSCk2TapM7vZnpzDQeORUlnTbzuik1RqZr5MxrwuI7w+v4BZrfWFv0+fp/GQCRNQ40Xv3fy7xdxZnF0ye/z+Th27NjIiWf79u14PM4nyQhxNhmMpgu6m650ghPTtLCB9lAAVVVoDwXIGtYZkZnLpau8b3Z+TXUmW3xPgaIoeNwax/rjWCV0owtRCY7OJv/jf/wP/uEf/oEjR47wt3/7t3zqU5/ijjvuqHTbhJh0zJxFImXgdjkP1JVMcGLbNom0yfRQYCTfuM+j09zgI3GG5LteMrcJM2ex69BgSe/jcWsk0yZhKYUpJhhHfTwXXnghTz75JK+99hqWZbFkyRIaGxsr3TYhJp1k2sSynXd7VzrBSTxlEqr3URs4sQcsVO8jEjfIGrlJX/SjPRSgqdbLn/cPcMG5obFfMIqAT6erP0mt342uSSkEMTE4OhJvueUWamtrufzyy1m1ahWNjY189KMfrXTbhJh0wvEMbpfzselKJjjJZHO4dY2WxpPXGGuqyrTmAKlMrqQZ0xOBoigsnddER2+cwWhp+bs1TcW2bXqHUmVqnRClG/WO+vbbb+fgwYN0dHScMCnMNE1UVa42hTiembOIJLIEfc4mI9m2zZ92ddNY6+HcaXVlbYtl2WSMHPOm1aOd5rsa8LoI1fsYjKYJ+id3lr9Fc5r4zaud/Hn/AB+4cOw11aMJ+HT6Iykaajz4PDKxTIy/UY/Cz33uc3R2dnLPPfdwzz33jDyuaRrz5s2reOOEmEySGRMK6Pbu6I1zrD/JVctnFDRD3Il40qAtFBgz0DQ3+AjHMximhauACXDFME2LjJkj4C3/RUHtO2uq39jfzxVL20rKMqYoCm6XyrH+BHPaasv+byNEoUb9Fk+bNo1p06bx9NNPn3QHnUzKmkMhjheJZ9B15yf1SiU4SaVNagJumhxk2tI1lWmhAAe7YtQGXBULSmbOIpkx8XtdJJImAX/571SXzp3Cj393gEPdMea0Fb9UC8Dr1okkskTiGeprJGOZGF+Ovi2/+c1v+PrXv04ymcS2bSzLIhwO89prr1W6fUJMCjnLIhLPEnDY7T2c4OSvF00t62Qu07Sw7HzRCqdBtzbgob4mSyJp4HfY/kLkchbJtMmsqTX4vS4OdcVIpsyyf9b86fV43Rp/3t9fcqAGCHh1jg0kCcrEMjHOHB19X/7yl/nkJz9Ja2sr//zP/8xll102Um5SCDE829t2HBwrkeBkZClW87tLsZxqa/Jjkw+q5ZSzLOJpk5ktNdQGPOiayqzWGlwulVS6vBnS9HfWVO85PEQ66zzH+mnfT1OxLJv+sEwsE+PLccKTtWvXsnTpUjweD/feey9/+tOfKt3aCqJtAAAgAElEQVQ2ISaNSDzreIy3UglOEqdZiuWES9doneInni49wA3LWRaJpMnM5iB1x1UR0zWVWVNrUFWFVKZ8nwewZN4UzJxddJ3q9wr4dHrDqbK3U4hCODqzuN1ustksM2bMYPfu3aiqSjabrXTbhJgULMsmHM/gcdiFXYkEJxkjh+s0S7Gcagh6qPG5yhKULMsmnjSY1hw85RivS9eY3Zov+1eOu99hbU1+QvVe/vxW8XWqj6coCi5dpWsgMemXsYnJy1Gg/uAHP8htt93GZZddxiOPPMJnPvMZ6uvrK902ISaFZCbf7e1kpnElEpxYlk0mm2NGS/C0S7GcUBSF9ilBzFxp6UUtyyaWNJgWqhm1dKTbpTGntQ7Lyl9olIOiKCyZO4WjfQn6I6WtqR7m8+jEkwbRhNyciPHh6Ft9/fXXc//99zN16lT+/d//nWXLlo1ZPUuIs0U0kXE82agSCU7iSYO2prGXYjnhcWu0NgaIp4q7yx0O0u1TAjTVjT1b2uPWmN1ai2naZMsUrBfNaURR4PX95bmrBvB7dTr7E5hlHsMXwglHZ5e/+7u/o62tDYCFCxfyiU98gqamytTMFWIysex8jWevg9rTlUhwkkwb1ARcNDoIik411nnxebSCi1zYtk0sZdDa5GdKvfMueJ9HZ3ZrLVnDKktN6Bq/m7ntdbzx9kDZCo/ouoqZs0mWcQxfCKccBer29nZeffVVrDOk4LwQ5ZLKmOQsHHV7Dyc4uWRhS1nWK5s5C9vOd1erZVz/rCoK00JBMlnn6UVt2yaaMGip9xEqIEgP83t1ZrfVks6YmGUI1kvmNhFLGhzoipb8XsNcukIkIQU7RPU56it7++23uemmm9B1Hbfbjf3OMpRXX3210u0TYkKLxLM4XQn15sFBXLpalgQntm2TSJnMbq2pSFGN4QpbfZE0NWOkF7Vtm1jCoLnBR0ujv+iLkIDXxazWWg52RQkoOloJa5fPnV6Pz6Px+v4B5raXp/fC49KIJLK023ZZL4yEGIujQP3YY4+ddtvLL7/MxRdfXLYGCTFZWLZNOJbG62Bs2LZt9nVEmNNWW5bAGk+aTKnzFrUUy6lQg49wIjtqha3hO+kpdV6mlhCkh9X43cxsqeFQT4ygTy96cpyuqbxvdhOv7usjlTHLMn6vqgq2ZZPO5PB7JQe4qB7HXd+n+g/gi1/8YkUbKMRElc6Y5Cxns717hlJEE1nOnV76aonMO4FzagXKYh5PU1WmhYKkRukCjyVNmuo8tE5xngltLHVBDzOag8STZkljzEvmNpGzbP5ysLQ61cdTFZV4WmZ/i+oqOS+erC0UZ6toIouqOQtO+zrCAMwrsRt2eCnW9ObSlmI5FfS5mFLrJZk6eWJZNJGlPuimrcxj5AANNV6mNQeJJY2ig3Vrk5/meh+v7x8oW7vcbpVwTAK1qK6Sv+lSWUacjWzbZtDhbG+AfR0R2qYESi4nGU1maWsKVLXrtaXRj6JywiSvWCJLXcDDtFD5g/Swplov7VMCRJNGUTcEiqKwZF4Tnf0J+sqUBtSlq2SyJoZZ3vSnQoxGMs0LUYRUJoeZsxzd1caTBsf6E5w7vbS7acO08Li0si7FckLXVNqnBEbSi8YSBjV+N9ObgyWVk3RiSr2P1iYfsURxwXrRnKZ31lSX767aBpIZCdSieiRQC1GEaDKD5jBIvXU03+1d6vh0KpOjPRQYlxnHtQE3DUEPg9EMfq/O9JbKB+lhzfV+Qg3FBeugz8W8afk11bkyral26yrRuCzTEtUjY9RCFMguIMkJwL6jEWoDbloais/DnctZ6JpC3TjVRlYUhdYmPy0NPmZOranK+Pjxnz210U9TnZdEERnTls6dQjxlsPdweSaVuV0a4UQWS859okpK/rZdcskl5WiHEJNGOpvDyNmO1vmapsWBY1HOnV5X0nyOZNqkpdHv+C6+EtwujfZQcFxqMyuKQnODHxu74JuDedPq8Hl0Xt7VU5a2qKoCtk1aKmqJKnE0I+XAgQN873vfY2Bg4IQvybe//W3+6Z/+6bSve+qpp/jWt76FYRjcfPPNfOxjHzth++7du/n85z9PPB5n2bJl/Mu//Au6rrNlyxY2b948kqb0iiuuYNOmTSOv27VrFx/96Ed58803C/pjhSiHWDKL5jDoHuyOYZgW504rvtvbsvIJhuoC5SuJORm5dJWgz03WyOFx2JsBoGkqi+Y0smNfH6svasPvLW1CH4CiKsTTRlneS4ixOArUd9xxBxdddBGrV692fFfQ09PDAw88wE9+8hPcbjc33ngjl1xyCXPnzh15zp133skXvvAFli5dyv/6X/+LJ598kptuuomdO3dy1113sW7dupPeN5VKcd9992EYhsM/UYjysW2bwWgaj9vZXeVbHWFcer7+crGSaZPmBt+43MlONE21Xg53xwoK1AAXzQ/xyp5eXnijmzUXTy+5HV6XRjiWpbm+smvZhQCHgdowDO6+++6C3njbtm0sX758pBzmmjVrePrpp/n0pz8NQGdnJ+l0mqVLlwJw3XXX8fWvf30kUB8+fJjvfOc7nHvuudxzzz3U1eVnzH7pS1/i5ptv5rXXXiuoPUJAPjd3z2ASTVVQ3/lP0xQ0Vc0/pigoSr57U+Hdn1Ul3/2aLxxh4/WMHTRt22bf0Xw2Ml0vLsgOryFuGKex6YnG79VRFEbSGDsVqvdx8cKpbN/dw8XnNdNQU1pGN11XiY6RtU2IcnEUqNva2ujo6GD6dOdXor29vYRCoZHfm5ubeeONN067PRQK0dPTM/LzbbfdxuLFi/nqV7/Kfffdx1e+8hWee+450uk0f/M3f+O4He/V1BQs+rXvFQoVf5d0thuvfXesL46i67jdGpZtY9uQtWzsnI3NKMUgjhsWbWj0U+Mbuxu6szdONJFl7aWzaKgPFNXeWDJL29Q62lre3V9n+3GXsRUi8QyBArud/2b5LF7d08sf3+zh41edV3I7VJeOv8Z7Vl1Ene3HXilK2XejBupPfvKTAPT19XH99dezaNEidP3dl3z7298+7WtPNeHj+Cvg0bY/9NBDI4/dcsstXHnllfT19fGtb32LRx55ZLQmj2lgIF6W0nehUA19fbGS3+dsNJ777mBHGFUBO1d8N7JlwlBm7KGX7bu7AGhr9DEUThT8Obadr+0cqnGP7C857sA2TPoHEmQDhQXqhvoAFy9s5o87u7no3CZam4q7eBqWzpq8fdhkZkttSe8zWcixV7zhfaeqSlE3i6MG6jVr1hTdsJaWFrZv3z7ye29vL83NzSds7+9/t7B7X18fzc3NxGIxfvzjH3PzzTcD+ZOVrus8//zzhMPhEyakrV+/nscee4xgsHx3yeLMZZgWmaxJTZUmZe3rCNM+JUDQV9yEo1TGpKHGg0e6Vk/g9+roan5YoNC13H+9aCqv7uvn19uP8ncfOrekmfgel0Y0YRTVDiEKMWqg/shHPgLAsWPHTnhcURS83tG7ey699FK+8Y1vMDg4iM/n49lnn+Vf//VfR7a3t7fj8XjYsWMHF110EVu2bGHlypX4/X4efvhhLrjgApYsWcIPfvADVq9ezQ033MANN9ww8vr58+ezdevWgv9gcfZKZ02gOifUfDayJKsuaCvq9bZtY+ZsptQVv/b6TKUqCo21XgaiaQIFXgR53Torl7TyzMsdvH0sWlIJTEVRsO38cSWzv0UlORqj/m//7b/R29tLMBhEURRisRiaptHQ0MDXvvY1LrzwwpNe09LSwqZNm9i4cSOGYXD99dezePFibr31Vm6//XYWLVrE5s2b+fznP08ikWDhwoVs3LgRTdN48MEHuffee0mn08yaNYsvf/nLZf/DxdknljTQqnRzWmo2skw2R23AXZbyjGeiuqCH3ki6qNdeND/ES7t6eG77Uea01pZ0N6yp+eNKArWoJMV2kD3gn/7pn7jkkku49tprAXjmmWf44x//yI033sg///M/88Mf/rDiDS0XGaMef+O17/YcHsSlq44SlZTqief20z2Y5PbrFxXVvRpNZJnbXndSAJDjLs+2bfYeGULXVcfL1hrqAyNzBd48MMhPfn+A9StmsWTulKLbYeYscjmYP6P08qUTnRx7xSt1jNrREb5nz56RIA35ses333yThQsXynpmMSlkjRxZh9nESmWaFge6is9GljFyBLwuuUsbhaIoTKnzkS6yOMb5sxtobfLz29eOnVAVrFC6ppIxc2QMKdIhKsfRWcs0Tfbt2zfy+759+7Asi0wmg2lKGj0x8aWzOahSbuaD3VEM02JekdnIMlmLlkZJpDGWGr+76HzbiqJw5bJpRBNZXt7TW1I7FNsmlZYbFlE5jjOTffzjH2fevHlYlsXhw4fZvHkzX//617nyyisr3UYhShZNZHEVmXSkUPs6IkVnIzPMHF63RqCK9aYnK49bw+/Ri046Mru1lrnttfzhjS4umDel6PkAbpdGJJGl/ixaTy2qy9GRefnll/PMM8+wfft2NE3jwgsvpK6ujkWLFsnSKDHh2bZNNJl1XO2q1M96qyPMOUVmI0ulc8yYWlPSsqGzSVOdl6O98aKzg33womn8vz/bxQtvdPGhvyoutajbpRJLyjItUTmOAvWzzz478rNlWbz00ksAfOhDH6pMq4Qoo6xhYVoWqlr5u9TuwRTRpMEVRcz2NnMWLl2l9iwvvlGIoM+NTeEpRYe1NPpZMreJV3b3cvGCZuqLSC2qKAqWDamsWXC2NCGccHTmevTRR0d+NgyDvXv3cvHFF0ugFpNCKmOgVKl08L6O/LKsedMKX5+bTJtMCwVR5W7aMZeuUuN3kzFMvO7iLsSuuKCdvxwc5LevdfKRlXOKeg9Ny6d7lUAtKqHgQA1w5MgRvvKVr1SkQUKUWzRp4HJVZ3z6raNh2kOBghNx5CwLTVWoC8rddKGGK2p5i9x1dQE3F5/XwrY3u1l+/lRamwqfyOdxaYRjGaY2lpaWVIhTKersNWPGDA4cOFDutghRdpZtE0tmceuVH5+OJbMc608WleQklc7R0uBHU6WUZaECvncrahVrxaKp+Dwaz+04WtTrdU3Np6iVZVqiAgoeo7ZtmzfffPOE4hxCTFSZbI6cTVUm+bx1NALAudML6/YeTsBTzPioAE1VqQ96iKWyRc/c9np0LlvcyrOvHOXtzgjnFJNaVFFIpg3JzS7KruCub0VRaGxs5Etf+lLFGiVEuaQyZhXHpyPUBdw01xeWnzuVzjGl3us4w5Y4WX2Nh8FYBl8J1zrLFjTz8u5enttxlDlttQVPTnPrKuF49qwqeymqo6BAbZomtm3jcsmECTE5RBJZ3O7KB0DDtDhwLMoF85oKOsHbto2FRaOc3Evi9+romvLOWH9x/966prLqwnZ++vuD7DwwyOJzmgp6vdulEk8ZJbVBiFNxdDQNDAxw6623snTpUhYvXszGjRvp6empdNuEKIll2cRTBu4qJDo52BXFzFnMK3B8OpkxaazxFr0OWOTlK2p5ik4pOux9sxvzqUVf7Sw4taiSHygnVWIbhHgvR2ew++67jyVLlrBt2za2bdvGsmXLuPfeeyvcNCFKM5w2tBrJQ97qiOAuMBuZbdvkrHzSDlG6uoCHXInDHIqi8MGLphFJZHmliNSiqqYQT2ZLa4QQ7+EoUB86dIhPf/rT1NbW0tDQwO23386RI0cq3TYhSpLMGFUJ0rZts+9omDnttQWNM6ezOer87qLX/4oTed0abl0tqcgGwJy2Ws5pq+WFN7pIZQqrZeB1aQzFMyXNQBfivRwX5chkMiO/p1IpSXEoJrxIPIu7CuunuwdTxJIG5xZYhCNrWoQKnHgmTk9RFKbUevM9KSX64LJppLM5/rCzq6DXaZqKkbPIGqVdLAhxPEeX8mvXruXmm2/muuuuA+AnP/kJa9asqWjDxJkrV4Z64GN/hkUyYxL0Vf5utZhsZOmsSY3PhV+Kb5RVTcDNscFkye8ztdHP4nOaeHlXPrVoXdD5dHLFVkikDTxVyC0vzg6OzhKf+tSnmDp1Ki+88AKWZXHddddx/fXXV7pt4gw0FEvTE80QCrorurY5nc1h21RpfDrMtAKzkWWyFtNChVfXEqPzuDQCHq3oilrHW3VB2zupRY9x7WWzHb/O7VaJxLM01srcA1EejgL1Jz7xCf7zP/+TDRs2VLo94gw2FEvT0ROnps5HXFMqWnwikTKoxgqZWDLLsYEkqy5sd/yarJHD79WllGWFNNZ66SihotawuqCHixc28+KbPSw/v4WpDmuEu/X8Mi0zZ8naeFEWjo6iWCxGMll6d5I4e4VjaY70xAn4dQIeF33hyh5PkUS2KhmiislGls7maGn0yzyPCgn68heA5ZjQtWJRK163xnPbnacWHV6mVY6xciHA4R21z+dj1apVzJ8/H7//3avKb3/72xVrmDhzhOMZDvfGCfp1NFXF49boTpukMmbRKR9HY+YsUhmzKuUi93WEC8pGZpgWbpdGsMCiHcK54VKh6WzxFbWG+d5JLfqr7Uc5cCzKnLZaR6/TNIVoIiv/zqIsHB3FMh4tihWOZzjcEyPo1U/I1qRpCkOxTEUCdbXuZPLZyGJcMG+K47vjdDbHtFBASllWWGOtl8NdxVfUOt5fvZNa9Nfbj3Lr1ec5+rf2uDXC8QytTdJzIkrn6Cz5kY98hFwuRzwel/WBwrHI8UH6PWN1PrfOQDRNc4Ov7ON48WQWTav8yXE4G5nTbm/btsGGGr+Usqy0gDdfUcuy7JInLeq6yorFrfz8xcN0D6YclcHUVJWclSNj5GSdvCiZoyPo8ccf50tf+hKGYQD5E46iKOzevbuijROT12hBGvLVrGzbJprMlj3PdSRRnQpG+zrCuHWVmQ6zkRmm9U5OaplgVGmaqtJQ4yGSyJZlCdyCmfX84k+H2XNkyHm9ahsS6dK734VwdAT93//7f3niiSc477zzKt0ecQaIJjIc7o4R8J06SA/zeXT6htI0BD1l6x40zBwZw6TWXdm7Vtu2eetovhyi08CbyVq0hSTBSbXUBT0MRDNjP9GBgNfF9OYgew+HWXWBsxn+nneWaTXJMi1RIkdnmLq6uqKC9FNPPcXatWtZvXo1jz322Enbd+/ezYYNG1izZg133303pplP17dlyxZWrFjB+vXrWb9+PQ888AAAO3bsYMOGDaxfv55PfOITdHZ2FtwmUVnRRIZDXTH8YwRpyE/6yRgmyQLTNI4mlc0Ble/27h5M5rORFTDb21ZsWZJVRcdX1CqHBTMa6A2nGIymHT3fpask0vllWkKUYtQzaTgcJhwOs3TpUh555BH6+/tHHguHw6O+cU9PDw888ACPP/44W7du5YknnmD//v0nPOfOO+/knnvu4ZlnnsG2bZ588kkAdu7cyV133cXWrVvZunUrmzZtGnn+v/3bv7F161auvvpqvvCFL5Tyt4syiyYyHOrOB2mnd5kuXWUg4uzE50Q8mcWlVz5Q7+vIL8ua6zAbWS5n4dLUqnTJi7xyVdQatmBGPkXsniOjn/uGKUp+eKfQfOFCvNeoZ9Ply5fz/ve/f2SMesWKFSxfvnzk8dFs27aN5cuXU19fj9/vZ82aNTz99NMj2zs7O0mn0yxduhSA6667bmT7zp072bJlC9dccw133HEHkUiEbDbLZz/7WRYsWADA/Pnz6eoqLA+vqJyRIF3gGKzXrRFJZDDM0k+mtm3n609XIb/3vuFsZF5ny2/SRq6sXfzCmbqAh3Ld0NbXeJja6GfP4SHHr9E1lZhU0xIlGrUfbs+ePQBYloX6njRPY91R9/b2EgqFRn5vbm7mjTfeOO32UCg0UuM6FApx2223sXjxYr761a9y33338ZWvfIX169ePtOeb3/wmV155pZO/8QRNTcGCX3M6IUkBCUAknmaoP0l7Sx26w9rPDfWBkZ81dwbV4ybkMPPT6aSzJoGBVEF5mYsRiWfoGkjy4Utnn/B3jEaJp5k1o7EsM77luHPOtm2imdzI0kCn/16nc8H8Zp5+8RCqS6cuMPZxVpuzSGdzTJkSPCMu0uTYK14p+87RgNmGDRv46U9/esJjN910E7/4xS9O+5pTLeM6/kAdbftDDz008tgtt9xyQkDOZrPcddddmKbJf//v/91J808wMBDHKkNRiFCohr6+WMnvM9nFklkOdkXxe3ViprMuvob6AEPhxMjvuZzFvrf7UMyGktYXRxJZwpEklsN2FGvH3j4Apof8J/wdp2PbNomUSSKWIp0obXKTHHeF02yL7t4409vqHf17jWZGsx8beOXNLi6aHxrz+QCxhEFHp16RnAHVJMde8Yb3naoqRd0sjnr784lPfIILL7yQvXv3cuGFF478t3TpUurqRh+ba2lpob+/f+T33t5empubT7u9r6+P5uZmYrEYjzzyyMjjtm2j6/kDPJFIcMstt2CaJt/61rdwuSTrz3gaCdKe0pYcaZqKmbNIpIyS2hNNZHA7vKMvxb6OMPVBN6F6Z7N5M0aOGr/rhIQvonpq/G7KNZ2rud5HQ42noO5vRaGsEybF2WfUM8dDDz3Ez372M5YtW8ZTTz018t8zzzwzMov7dDnAL730Ul588UUGBwdJpVI8++yzrFy5cmR7e3s7Ho+HHTt2APmZ3itXrsTv9/Pwww/z+uuvA/CDH/yA1atXA/nJZDNnzuRrX/sa7govvxGji6eMd4N0GYKj26XRHy5+Uplt28SSRsmFGMZimDkOdkU5d3q9467MrGFRV4V0puLUhitqZcqQsU5RFBbMqOdgd4x01lnwdbtUwrHyLBMTZ6dR+2KCwSDBYJDvf//7p33Oxz72sZO6xSF/x7xp0yY2btyIYRhcf/31LF68mFtvvZXbb7+dRYsWsXnzZj7/+c+TSCRYuHAhGzduRNM0HnzwQe69917S6TSzZs3iy1/+Mrt27eK5555j7ty5XHvttUB+3Pu73/1uibtAFCqeMjh4LFK2IA35lIvRRJZMNldUHd+MkSOXKz0L1VgOdsUwczbzCliWBeB3OOlMVEZjrZeYw8A6lgUzG3jxLz28dTTCojlNYz7fpavEk6ZU0xJFK3nQZLSUoldffTVXX331CY8dH1gXLFjAj370o5Net2zZspOC/8KFC9m7d2+JrRWlMswch7tj+MoYpIdpqspQPM3UxsIn/FRrCcy+jjBul8qsFmcTQ0zTwuvWK36nL0YX9LmJZdMjWRVLMS0UIOhzsedw2FGgHv68VMaU9LGiKCWfac+EmYzCue6BJAqUPUgD+Dz57u9iElREHS7Lsiwby7KLylk/ko2srW7MZC7D0tkc9TVych5vLl2lPugpW/f3/On17O+MYJrOjlVdV2SZlija5J6GKKoqkTYYjGeo9VemG1dVFSwbYkmD+gKWWFm2TSxljJn1K5bM8u8//QsZ492TtaooKEr+sxVl+HcFVcmfkFX13Z/z71FYNjLLtgmWo4STKFlLo5+DHUN4y7B6b/6Menbs6+PAO/MVxuJxaYTjWaY22VI5TRRMArVwxLJtOvsS+NxaRXtRvG6VvnCqoECdyeawrLG7NN86GiFj5Hj/+S24dBXbzv9d7/7/nZ8t+z2Pv/N/y2ZmS5AFMxsctcuybHRVweuRbu+JIOh3Ux/0kMwYJS+Vmt1ag8elsedI2FGgVlUFM2cxGEkzxWHtciGGSaAWjgzFMqSzOWoDlZ0U5XZpRBMGqYzp+GSazJiO7lL2H41QG3Bz5bJpVRmyyRg5agMeuYOaQFoafezryOB1lzZWrWkq86bVse9IGOv9ziYxBn0uOgcSeD06QZ9MLhTOFT3QOFxAw+8vLZuUKFzOsqpaF9wwc3T1Jwj4qnNnqGs4LnwAEI2PPT6dy1kc6Ioyt722avMqTNOmNijd3hOJ160zpc5HMlX6WPWCmfUkMyZHeuOOnq+qCn6PzuGeGFmjPPnHnUplTDr741U9b4jycRSo31v5av/+/dxwww2n3CYqy7ZtDnfH6B489fr1SugZSqIoStUSdvg8OkOxjKOqQ5ZlE08buMaY3HakN07WsJjbXtiyqmINnxD9kzwb1ZkoVO/Fxi45Q+Hc9jo0VWHvEefJT1y6igJ09JYnQ6ITGSO/9r9vKCWJVyYpR2fe//qv/+LnP/85AP/xH//BTTfdNJJ3W1RXIm0SSxr0DqWIxCufRCGZNhiIZvB7qzfOOlx1KJIYe5ZsOmti22OvPtjfGUFVFWa31ZarmaMyTKvgAiWiOly6RnOjr+RMeG6Xxpy2WvYcDhd0p+r36iTSRlUuts2cxaGuGKoCXo9GTxUv8EX5OLrc/973vsff//3f85//+Z+oqsoTTzzB7NmzK9028R62bdM1kMDn0VBVhY7eOF63XlSCECcs26azP4G3whPITsXn0ekbStJYM3rFqWTaxEmOk7ePRpnREqxamcmsYdE6RSYNTVRNtV76w+mSk5AsmFHPW0cjdA+maG1yPgxY43fRF07h92jU1zhLRVuonGVxuDtGLmfh9+VP9dGEQTJtSAKeScZRPWqXy8XmzZvp7u7mk5/8JA0NDWNWzxLlF01kSWVM3C4NXVPRNYUjvbGi1h07EY5lSGdy41JDWddVsqZNIj16V10kkcHjHv1EG0lk6Q2nqtbtDWDDmMvFxPjRVJW2Jj/JVGldwefOqEdRYE8B3d+Q7wEK+nQ6+uIVSdYzvEojmTFHgjSAW1foGUqV/fNEZY16Jlm+fPlIN+SwT37yk0D+QNu9e3dlWydGWJZN12DyhJnQXo9OPGnQNZCgvcxl9AzTomsgib9KE8hOxa0rDETSp50ha+YsEukcNf7RA+L+oxEA5k2rTqDO5Sx0TRmXCxzhXG3QgzecJmvkis4cF/C6mN4cZO/hMKsuaC/otZqm4nZpHO6JMbe9rmzDJLZt0z2QJBzPnrRKw+vRiSayJNMmfrmQnDQc1aMGRlLv5XI5LMuSylVVFk5kyBrWSV+8gERYM1UAACAASURBVE/PjyF7XDTWlq8LrXcoP5Y1nhWfPG6NSCJD1vCf8kSazuZQGHuZzf7OCHUBN1PqKtPF+F4ZI0djjVey9k1wqqLQNiXAgWORklK8LpjRwLOvdDAYTRf8HfS4NBJJk86+BDNaynOxPRBJ0xdOUnuaQjAuPZ+rYOZUqS09WTg6C7/00ksjk8cOHDjAFVdcwWuvvVbRhol35SyL7oHkKSd0DXehHS1jF1oybTAQSVd1Atmp5DOEKYRPU8M5nsqOuX41l7M4eCzK3Pa6qgXOnJVPriEmvqDPRY3fXdJ3Z8GMfMKTPUeKGw4M+HUi8Qx9kdK7pIdiaTr749T43ac93r3vXABXKz++KJ2jQP1//s//4Ytf/CIA8+bN4zvf+c7I76LyhqKZUSe9aKqK153vQnOypGk0lm1zrD+BZxwmkJ2Kz6vRP5Q65VKWaMIYs3v5SG+crGkxt0rd3rZto5DPWy4mh6mNfgyz+NwE9TUepjb6Cx6nPl7Q76KrP0m8hJno8ZRBR2+coN816gWsoijoWv6uWkwOjgK1YRicf/75I7+ff/75ZLOSYL4azJxFz1CKgG/08SS3SyOXy08gKSWpQTiWIZk2KzaTvFCaqmK+s1b6eIZpkcmaYxYH2X80gqYqzG6tTjdf1rCo8bvGdchAFMbn0Wmq9ZJKF5+EZP6Meo72Jogniwu0qqrg9+oc7i4uGUoqY3KoK4rPozs69nwejXA847imthhfjs4mPp+P3//+9yO/v/jii5KRrEr6Iyls23b05Qv4SutCG55ANtZFQbV53Bp975mpmj/BOEgb2hlhRkuwamUms0aOutOMDYqJK1TvI2cVnwRlwcx89/fejuJXw7h0FVWFIz2FJUPJvpPQxOVSx0z8Myx/V63QX4budlF5js7Id999N5/61KfQ9fzTVVXlG9/4RkUbJvKpO/uG0gUFzqDflZ+t7XEVnE+4dygJNo5LOFaLx6URTWRJZ0287vy+iCUNNG30QB2JZ+gLp1k6d0o1mpmnKLJGdRJyuzRaGn30DKWoKaI6XHO9j4YaD3sOD3HR/FDR7fB5dGKJbH4lRyg45vPNnMWh7hiKQsGrDHwencFohlCdf8L0oIlTcxQBlixZwvPPP8++ffvQNI3Zs2fjdstdQ6X1hdP5MotOMnq84/h8wvPa6xzfSSbTJv3RdMVKWJZKUxWGYhlam4YTN2TwjJHfe39nFKBq49OmaeFxaVW7exfl1VjrpT+SJpezCr5YVRSFBTPqeWl37wkXlMUI+l30R9L4vToNoyRDsax8OmHDtIrqBVMUBV1VGYimaZsSKLq9ovJGPRq3bt0K5NOGPvroo7z00kts27aNxx57jP/4j/+oSgPPVplsjv5IqqiZ14XmEx6eQOZ1TYwJZKfi8+gMRPKZpLJGjmzOHvNk+tbR6i7LShs56mvkAnay0jWVqU1+kkWOVS+Y2YBl2SPr9oulKApBv87R3tOv5LBsm6N9cZIZs6ShKp9XYyCSqnqREFGYUf+FDx8+DMC+ffuq0hjxrt6hJLqmFB04/d58YoOeoSStTaNfLYfjGZIZ47TrLicCVVWwbJtYMpsfrx9jwpyZszjYFWXxOU1Vu/iwLZugd+LuQzG2+qCHvqHikqBMCwUIeHX2HAnzvjlNJbVDU1Xc7tMnQ+kZTDIUy1BXYnU2RVFQ37mrHus8IcbPqIH69ttvB5ClWFWWypgMxTNFjZUdr8bvoncohd+jUxf0nPI5/397dx4d1X0efPx7Z981WmYkIYHELpZgvGG84eAEg1ksvOQkIacyr2PS9rQ9pyRxHCfHx6mPm7g5PnGaze/bJm3axk5M6xiCY2OMXewEsGOwAxgwIEAs2vfZZ+72/jGSzCLEjGZGmkG/zwnH0tzRzE83V/eZ3/Y8ijqwgKwAshTZrUY6+2I4rKYrV8tqDyGP4bYsTdMxGiRsYltWQUsmQXFwsiWQdqCWJInZU7x8dLIHRdGuuCPhSqxmI+HopclQuvqidPRFs1Yb3mEz0tUXo6zIhtkkrt98NOLdec2aNSP+8NatW7PaGCGprSeSHL7OsCcoSRJOu2nE4h0dvVF07crDyNlysiVAtSZxhfTcwzKbkovKZEW74j7lxuaBbVljlH0pLqt4nFYMeTp1IKQumQTFTCyuYEuzTGndlGI+ONbFydYAsyZ7M27L0E4OmxG/10FfKE5zdxi33Zy1kaLBxELd/TEqRK86L414u3z88cd5/PHHefjhh1EUhccff5x169bR0tIylPNbyK5wTCYQTlyQ0zsTJqMBk2n44h3RuEJX/5X3aGeDpum8/qez/Gr7Mf7j1cOj3uttNhlIKOoVF9g1nhvbbVmKouPJcBhSyA+SJFFR4iQhp58EZWqlG6vZOOosZcNxOcy0dUXo7Itwpj2Iy2ZKa4FpKuw2I139MWQlNwV+hMyMeIdetGgRAOvXr+eLX/wiixYt4pprriEej7N582ZWrVo1Jo2cKHRdp7Urgi3LWyVsFhPB8IXFO/SB6joWc+Y99yuJxGT+5+2TNLUGqfY5OdcZ5lRrkGmjqA1tt5queH76gnG6+mNcN2tstmUN3swdWfpwJYy/5IprK8GonFbxCqPRwMzqIo6d6UO7Wc9KQB1MhnKuM4zLbs7J6NdgO3sCMcpLRI6MfJPS/+O9vb00NDQAYLVaWb9+PZ2dnVf8ua1bt7Jy5UqWLVvG888/f8nxI0eOcP/997N8+XK+/e1voyjJFY6bN2/mtttuo76+nvr6ep599lkAWlpa+NKXvsSKFSv467/+a8LhcMq/aCEIRWXCcTknexpdjuSq6d5gMm92fyhOOC5ntI0kFW3dEf516xHOtoe459ZaGpbPxuUw897h9lG/ZipFOIAxK2spKxoOmylr1Y+E/OAvcaBqWtpJUOpqvETiCmc7Qllri8lkoNhtTTmhyWg4bCY6+6IZpyHON7Ki0dwZorkzRH8omeN8tIltxktK/6+rqkp7+yc31q6urisOCbW3t/Pss8/ywgsvsGXLFl588UUaGxsveM4jjzzC448/zuuvv46u62zatAmAgwcP8s1vfpMtW7awZcsWNm7cCMA//MM/sG7dOrZt28b8+fP52c9+ltYvm8+SW6Qi2HOUeCC55cPMuc4QwUiC5jFYQPbRyW7+7dWP0XSd9SvrWDizDJPJwK0LJnH8XD9dOco13Njcj9dloXSMtmUlZA2ve/jFekLhspqN+Lx2IleoiX6x6VVFGA1SRrm/x4PBIKHr0BOMjXdTsiYUlTl+ro++UIJAJMGZ9hCNzf0caurhRHMf7T0RgpEECVnNKPVyrqUUqNevX8/atWv5xje+waOPPsq9997Lww8/POLP7N69m8WLF+P1enE4HCxfvpxt27YNHW9ubiYWi7Fw4UIA7rvvvqHjBw8eZPPmzdxzzz18/etfp7+/H1mWef/991m+fPklz78aBMIJ4rKS0znVweIdTW1BdE3PWQ9Q03TeeP8sv33nFJWlDjasmUvVeQkVbl0wCaNB4r0jHVl/7+S2rCAzqseuWpauUxCr5oX0lRXZB8r7pt7LtJqNTJvk4ePTfXl98x+O026io6fwe9WqptHaHeZESz9mkwGn3YTdasLtTC4UdNlNaLpOVyBKU2uQo2f7OHy6l9NtAbr6o4Rjcl6dg5TuLg888ADz58/n3XffxWg08uUvf5lZs2aN+DMdHR34fJ+k0vP7/Rw4cOCyx30+31Cv3efz8ZWvfIUFCxbwgx/8gCeffJJHH30Ul8s1lMb0/OcXOk3Tae0Oj8kcp8VsRJLIWZCOxhVeevskJ1sC3FDnY/mNky+ZU3M7LCyYXsr+xm6WXluV1QL2p9uTmZrGathb1TRMJint9I1CYRhMgtLcGcbjTP1vpm6Kl+Pn+mnriVJZWjhzvgaDhE5ynUeZ1z7ezRmV6MC0Qzyh4HEMvzpekiTMJuMF29E0TSeWUAlG5KE0DRazEbfDhNNuGdcc/infIevq6qirq0v5hYf7JHn+CRvp+E9/+tOhxx5++GE++9nP8o1vfGPE10tVaemV8+emyufLztafzt4oDqcVr2tshmpzpaUzxC9+/zH94ThfWDaLm+ZVXva5y26q5cPjXRw+08eyRTVZa8O5/a0YjRIL6yrGJHiGojLVlXb8/uxdV1eSretuIhrNuSspdSHrEkajhCXFfcaL5pt5Zc9pTneEmDt99Lm/x4PHrRFJKJSUOC/5kJ3P156u63T0RukIhiny2NPeWjccRdGIKyq9UZnaycUZdXAyOXc568KVl5ezd+/eoe87Ojrw+/0XHO/q6hr6vrOzE7/fTzAY5KWXXmL9+vVA8uSbTCZKSkoIhUKoqorRaBx6frq6u9OrTHM5Pp+bzs5gxq+jqBpHz/ZhMxvoVQo3jd+hUz38blcTNouRB1fMptrnordv+MV+xV4nNhNMm+ThnQ+buXZ6SdZWsh462U1NuZtIOEYkK684smA4QZHNmJVrIRXZuu4mokzOncMETW3BtJKMTPa7+POxThbPSf8+Nd4C4QTHjRKlnk86D/l87SVkleauEMGwjMthJqqoRKPxrL1+MJKgszM46kA9eO4MBmlUncWcLSG85ZZb2LNnDz09PUSjUbZv386SJUuGjldVVWG1Wtm3bx+QXOm9ZMkSHA4HP//5z9m/fz8Av/rVr1i2bBlms5kbbriBV1999YLnF7qeYCyZ1apAVwxrms6Oved46e2TlJfY2bB6LtUpVP0BWDy3nFBU5lBTdhbd9A5syxqrYW9d10GSrph8RSh8bocFp81EPI2c2HVTiunojdITKLzFWU6bifbu8CW5F/JRfyjO8XP9ROMKHpcl63vM80HOokN5eTkbN26koaGBtWvXsnr1ahYsWMCGDRs4ePAgAM888wzf+973uPvuu4lGozQ0NGA0GvnhD3/Id77zHe6++24OHTrEI488AsATTzzBpk2bWLlyJXv37uXv//7vc9X8MSErGh090YJdiBSNK/x6x3F2f9TG9bN9PDiw9SpV06s8+Lw23j3UnpVFN0PbssYobWhC1vA4zCnVChcKmyRJlHntxBOpB666KcnMZNlMfjJWjEYDqg79ocR4N+WyFFXjXGeQprYgVovhqi4vK+mFtiwxQ/k09N3WHaarP5ZWcMsXHb1RXnyrkf5wgpWLp3DdrNTn4Yq9zqFh8Q+OdfLK7tM0LJ9FbWX6CVDO9+sdx+nqj/G3980fkxXfwXCCyX4X3hFKEWZbPg8/5rtMz52iahw53YvLbkr5+vqX3x3CZDLw0Mo5o37f8aKoGrKsMXtKMQaDlFfXXiQmc6YjhDJQ4jPXf+/BSII5NSVX39C3MLKErNLZH83qiuexcqSpl1/8/giyovHgitlpBemLfWpaKQ6riXcPZ7ZVS1EGtmVVecZuW5YkXdWf4oULmYwGvC4LsUTqw9+zpxRzriNMKCLnsGW5YTIakFWd/nD25nozpek6Hb0Rjjf3Y5CS6VXztTRvNolAPU46+6IYJUPBzafs+aiN/955An+xnQ1r5jA5w9XOZpOB6+t8HDvbR3cGc3mn24Moqsb0MZqfVhQNm9k4ZrnEhfxQ7LahpJEPu64mOfx99GzhDX8DOKxG2noiaHkw8BpPqJxsDtDWE8VtN0+ovz0RqHNM0/XkEJKiEk+oROMKgXCC7kAMu62wLrTmrjA79p2jrsbLgytm43ZkZ1/hjbP9GA0Sf8ogrehQtazKsdk+EpNVvG5RhGOicViTBTFSnT7ze+0Uu60Fl6VskMlkQJY1AuHxm6uWFZWuvijHzvYhqyoep7ngOjiZKrxx1zygahq9wWTebFXVUDQNVdUH/puc21EHHtMBdGDwuhpYKWw1GwtqyEZRNLb84RQuu5l7bqnNasIUl8PM/Gkl/Lmxm09fWzWqymGN5/qprXCPWT1dXdNx2USgnmgMBolit43eUAxnCtMekiRRN8XLe0c6iCWUnOfWzwW7zURbT4TpNWPXq9Z0nUhMoScQoy8UR0LCYTdO2IWbE/O3zlA0rnLsTC/nOoO090boCcQJRhLEEiqqpmIwJNMIOu0m3A7zUNq65NcW3A5zTgpv5NLOP7fQ1R9jza21WUkkcLGb5pYjKxofHLtysZeL9QbjdAfiY7baW9OSVZFsYlvWhOR1WVDV1INWXU0xmqbTeK4/h63KHbPJQEJW6QnEcl7MIiEne89Hz/RysiVAOCYP3UMnapAG0aMeNZPBkLWh33x3tiPEnkNtXDerLGd7lCtKHNRWuvnTkQ4WzytP649y8AY4VoE6Lqt4nFYMBTQiImSP3WrCbDKgqlpK+Q+qfU6cNhMfn+lj/rTSMWhh9jmsJppaAwQCUVx2Mx6nBbslWXI202Howd5zV3+UYDiBJCU/BNutEzcwX0ycCWFEsqKy5Y+n8DgsLLtxck7fa/HccoIRmSNpJkA53txPsdt6QRalXFIUnSLXxPiQJlxKkiRKi2zE4qktKpMkiXlTSzhyupcTzYXZqzaZDHhdVlx2EwlFpbUrzInmPg419XCqNVnIIhJLr3xkXFbp6Ivw8eleTrb0E4sruBxmXA6zKBl7EXE2hBG99UEzPYE499xWm/Pc2TOriyj1pJcARVE0mlqDY5aNDEBHH5MCKkL+8jgsqHrqq7/vvK4Kv9fOS2+fzGh3w2goqsYru5s43NST8WtJA+trXAPTeC67CVlRaeuOpBS4NU0nEElwqiXA0TO9dPREsZoNeJwWbNbc74cuVCJQC5d1ui3Ie4c7uLHOx9QMk5GkQpIkbprrp6U7wtmOUEo/0zSwLWtGde7bB8k5NKfNJD7xT3A2iwmbJRmkUmExG/n8nTMwSBIvvtlILJFejevR0nWdzX84xQfHuvj9njNZf19JkrBcKXC3JAN3e0+Ej0/30NQaIKGoyXKTDnPBpk8eS+IMFZjTbUF+tf0YvcHcJiFIyCq/29VEsdvKZ66vzul7ne+aGaXYrUbeTXGrVuO5fkxGidqK1LZlxQfK2I3uX4JoQh3TTGRC/kpn+BvA67bywNLp9ATi/PadUzlfmKXrOtvfP8vhpl6unVlGNK6w51BuSwMPG7jVZODu6othtRrxOC1YLYW162W8ifG7AqKoGlt3NdETjPMfr31Mw4rZlORoXnbHvnP0BuM8ePfsMU0sYDYZuW6Wj10H2+gNxil2W0d8fmNzPzVpbMuKyyo1FR4cV1ixPdLIu9kkPt8K4LabAR1d11MOOrUVblbcNJlX3z3DWx8089kbcvcheM+hdt473MFNc/zctWgycVnl3UPt3Fjnx2Ufm4x6g4F7IiUnyQVxxykg7x5qpycYZ9kN1Siqzn9sO0p3f/bnu062BNj7cSc3zS2npnzs688umuPHkEIClJ5AjJ5APOX56cEbqstuGioaf7l/gzeX4f6JnoAAyeFsl91CIo1MZQA31Pm5YbaP3R+1ceBEd07advBENzv2nmNubTF3LZqMJEksvbYKRdX4w/7WnLynkDsiUBeIQDjBHw60Mnuyl5vnV9CwYjaqlgzWXX3RrL1PPKGydVcTpR4rd15XlbXXTYfbYWFebTEfHu8acU6tsTkAJBehpSIha7jsE3s/ppBdxR4rCTn9UpDLb5pMTYWbrbuaaO4avm77aJ1o6WfLriZqKtysvX3q0AfL0iIb1830se9YZ86nzoTsEnesArFj7zk0TeeuRcktUv5iOw+umI2uJ4N1R292gvX2988SiCSov23quA7xLp5XTkLR+PBY12Wf03iujxK3NeXh/4SsUiy2VQlZNDiEnG4RQqPBwOc+PQ23w8ymtxoJRrKTorO1O8x/v3WCsiIbn79z+iWLHpcsrMQgwc4Pm7PyfsLYEIG6AJxuC/LRqR5u+VTFBXO2Pm8yWEuSxH++fpT23khG79N4rp8Pj3dx87wKqjMstpGpylInNeUu/nSkY9hFN7Ki0dQWTCvJiQ6i2pWQVSajgSJnehW1BjlsZj5/5wxiCZUX3zqRVrGP4fQG47zwxnHsVhPrPjtz2HSlboeFRXPKOXiyh/aezO4XwtgRgTrPaZrOtvfOUOS0cNunKi45XjYQrI0Gif/cdoy2Uf7xxeIKW3c34fPa+PTCSZk2OysWzyunP5zgyOlLE6CcbguiqHrK89OyomGzmMSiFiHrit1W5FEG2fISB/fePpWWrjCv7Dmdds98UDgm8/wbx1A1nXXLZuJxXn7k6NZPVWCzGHnrA9GrLhQiUOe5vUc7ae+NsuzG6suubC4tsvHgitmYjRL/9fpRWrvTn/N6/U9nCUVl6m+biilPVjXPrPZS4rYOu1WrsTm5LasmxW1ZiYRGyRVWkAvCaDhtZoxpVNS6WF1NMZ++dhIHTnTz7ii2TyVklV/vaCQQTvDFz87E57WP+Hy71cQt8ys4fq6fM+3BUbVZGFv5cUcWhhWOyez8sJnaSjdzaopHfG6Jx8aDd9dhMRn5r9eP0ZLGApWjZ/rYf6Kb2z5VyaQyZ6bNzhqDQWLR3HKaO8OXJED5pFpWapewio5zjLakCBPLYEWtTJKJ3L6gkjk1xezYd47jaRTv0DSdl94+SWt3mPvvmJZyffib5ia3aL25r3nUvXhh7IhAncf+94Nm4rLKikVTUtoSVOy28uCK2dgsyWB9rvPK2b0iMYVXdjdRXmxnyTWV2Wh2Vi2cUYrNYuS983rV3YEYPcE4M6q9Kb2GpumYDBK2AqtYJhSOdCtqXUySJOpvq8XvtfPbt0/SlcK2S13XeWXPaY6f62fl4hpmTxn5w/z5zCYjdyycxNmOUFofDITxIQJ1nmrpCvPBsS4WzSnHXzzyUNb5vAPB2mEz8fz241dMxbntvTNE4yr1t03Ny1R+FnMyAcqR0730DWwpGaqWVZVa2tBYQsXrsoj9z0LO2K3JtLKqOvoFYRazkc9/ZgZGo8SLbx4nFh+5h77zzy38+XgXS66p5PrZvrTfb+HMUkrcVt7a15zzLGlCZvLvziyg6zqvvXcGp83EHQvT7+UWuZLB2mkz8fz2Y5edhzrS1MtHp3pYsrCSilJHps3OmRvn+AH408cdQHJ+usST+rYsVdXwOMX8tJA76VbUuhyvy8rnlk6nN5TgpbdPXjaA7v24gz/sb+XamWXcMcrFn0aDgaXXVdHRF+WjU5kX7BByRwTqPHTgRDfNnWE+c331sFssUuFxWnjw7tm4HWaef+M4p9suDNbhqMzv95ymstTBrcOsJs8nRU4L82pL+PBYF6GozOm2IDPTyEaGJGG/QspQQciUO82KWpdTU+5m5eIpnGgJsGPfuUuOf3y6l9feO8PM6iJW3VyT0UjR3NpiKkoc7PywGSWD0QAht0SgzjOxhMKOveeo8jm5ZkZmRebdDgsNK2ZT5LTwwo7jnGpNZvLSdZ1X3z1DXB4Y8i6ATF03zSsnLqts+cMpFFVnehrZyNwiG5kwBuzWwYpamQe862b5uLHOz7uH2tnf+EnSnzPtQX77zkkmlTq5/45pGAyZTedIksSd11fRF0rwwbHOTJst5Ii4e+WZd/a3Eo4p3H1TagvIrmQwWHtdFn69o5GTLQEOnerlyOlePn3tpLTmv8dTVZmTyX4XJ1oCmIwGalPMQZ6Qk/PTgjAWksPf6Sc/Gc5di6qprXDzyu7TnOsI0dkX5TdvNlLktPDFz87IWk6A6ZM81Fa4eWd/K3E5O23PJ/2hOAdPdtPeEynYuficBuqtW7eycuVKli1bxvPPP3/J8SNHjnD//fezfPlyvv3tb6MoFy6eOHz4MPPnzx/6PpFI8LWvfY01a9ZQX1/P7t27c9n8MdfZF+VPhzu4dmZZVrdJuezmgUpbVn7z5nFeffc0VWVObp439kPesqLSGxxdIZHF88oBmFrpTnmvt8hGJowlt92MTnaCgdFg4IFPT0+mGf3fEzz/xnFMRgPrls3K6jU92KuOxJQLdlcUOl3X+eBYJ89tPsTL75zi//3uME8//yH/9vsjbHvvDPsbu+jsixZE8M5Zmcv29naeffZZfvvb32KxWPjCF77ATTfdxIwZM4ae88gjj/DUU0+xcOFCvvWtb7Fp0ybWrVsHQDQa5cknn0SW5aHnb9myBU3T2Lp1K0ePHmXDhg288847ufoVxpSuJzOQWcwG7rw++8UwnDYzDctn8V+vH6M7EKP+9qkZD5ulKy6ryLJGcbGZ7u4E1jS3S82e7GXe1BIWTE9tSkBkIxPGmsVsxGkzk5DVrFx3DpuJL3xmJv/2+yMgwfoVdVcs/Toa1T4XdVO87P6ojRtm+wr+w20gnGDr7iZONAeorXSz9Noq+kNxWrojtHSF+fB4F386kpyiMJsMVJY6qCx1MqnMwaRSJyUea17tEslZoN69ezeLFy/G603udV2+fDnbtm3jb//2bwFobm4mFouxcOFCAO677z5+9KMfDQXqp59+mvXr1/Phhx8OvaamaUSjUVRVJRqNYrPlphbzePj4TB+nWoOsuGkyzhz9kThsZh5aVUckplDkGttV0LG4gqbD9KoiirwOWtoCaQdqg0Hi/jumpfz8REKjvKQwhvaFq0epx8qZjlDWPiD6i+38n1V1GCTpilnHMrH0uiqObunjDwfaWD5Q/KfQ6LrOwZM9bHvvDIqqs+KmKdxY50OSJCb7XcyflvyQr2k63YEYrQOBu6U7zL6jnbx3OBm8rWZjMngPBG6fd3xjTc4CdUdHBz7fJ3v7/H4/Bw4cuOxxn89He3ty2OXNN98kFouxYsWKC17z3nvv5eWXX+b2228nEAjwgx/8IO12lZZmXmzCFknQE+6l2Jud4emErLJj7zkqy5x89qapGMe4p5troZiM3WFlxmQvNosJTdMpK3VhNRlzmq7UYI5TM7m44HsHF/P5xr5G+NViLM6dt9hJIKHhcWRv73627jVXeo8b51Sw92g7yxfXUjzM9sexaMdoBSMJ/vvN4xw80UVtpYd1d83GV3z5baelJS5m1X7yvarpFcUMmwAAGjdJREFUtPeEOdse4mx7kLPtQf50pANV1Zk52cvyWy+tRpaOTK69nAXq4dLSnX/RXu54Z2cnzz33HL/85S8vOf6Tn/yEhQsX8utf/5qmpibWr1/PvHnzqKpKfai4uzuU8ZxEKJocju/ty04d2Z0fNtMbjNNway2BwNVV0SYcUbBYDNRWeAj2RwmSvGAtkk5LewC3MzdBVNN0YgmVUCBKeJRz4vnI53PT2SnyM4/GWJ47XVZp6whgs+bsFpsTN8/1se/jdn73TiP33Db1gmPFXmfW7nnZdqSpl9/vOU1cVvnsDdUsnluOQdLTbq/dJDGrys2sqmRQVVWNjr4o6NDZGRx1oB689gwGaVSdxZx1Z8rLy+nq+mRbQUdHB36//7LHOzs78fv97Ny5k76+Pr70pS9RX18PQH19PaFQiDfffJP77rsPSZKYOnUq11xzzQW99ELUG4yz62Ab82qLqU2xwEShCEZkHHYjUys9l+TkLhpIQJKrPMPJbGT5Nc8kTBwlHiuJLGzTGmtFLis3zvGz/0Q3nX3ZqXGfS9G4wm/fOcl/7zxBkcvChjVzuWV+RdbW3xiNBipLnbgc4zsql7NAfcstt7Bnzx56enqIRqNs376dJUuWDB2vqqrCarWyb98+ADZv3sySJUv43Oc+x44dO9iyZQtbtmwBkovIXC4XdXV17NixA4Cenh4++ugj5syZk6tfYUxsf/8sBoPEshsLc05oOLquEwjLFDkt1JR7hv0UajYZKPFYiWZpK8vFktnIxLYsYXw4bWYM0ugrao2n2z5Vgdlk4H/zvAzm8XN9PLf5EIdP9XLHwkk8tKquYLabpiunPeqNGzfS0NDA2rVrWb16NQsWLGDDhg0cPHgQgGeeeYbvfe973H333USjURoaGkZ8zccee4yDBw+yatUqHnzwQb761a9SW1ubq18h504093P0TB+3Lai8aoKKrusEIjKlHivVfteIn2yL3TaUDAoZjNQGkY1MGE8Gg0SJx0osUXj7kh02M7fMr+DjM32cu0KtgPEQT6hs3dXEr3c04rCa+PLqOdyxcNJVndRI0idYjbNszVH3hGU0ZfRl7VRV4/9uOYym6/z12nkZLVLIF5qmE4zIlJfYKS92XHbY+fy5wuPn+tB1PatbqOIJFZvFSE1FakU7ComYox69sT534ZjMieb+gvwQnpBVfvzSQcqKbDSsmI0kSXkxR32qNcDv/thEIJLglvkV3LFw0pjcO4ORBHNqSq6+OWphZO8d6aA7EGP5oslXRZBWVY1gRKbK56SixJny3LDPayeeyO5cXlxWKSrAm6NwdRmqqKUV3ly1xWzk9msmcbo9xImWwHg3h4Ss8tq7Z/iv149hMkr8n7vr+Mz11VfFvTMVhbUk8SoRjCR4588tzKwuYtbk1Goq5zNF0YjEFKZUuClOc3+222HGYABV07I6dHW1bckSCo9Bkij12Ojqj+G0F15AuX5WGe8eauOtfc1MnzS2o1ORmExXf4zu/hhdgRgfn+6jNxjnpjl+7ry+CrNpYk1riUCdJl3X+e07J+kOxDAbJcwmIxaTAbPJgMVsuPB7kwGz+cLvLSYjb+5rRtX0gk0qcD5ZUYklVKZO8uB2pN+LNRoM+Lx2OnqjuByZ38xkRcNuFdnIhPzgcVro6M3/1dPDMRoNfPraKjb/4RSHm3q5rTjzHBTnU1WN3mA8GZADnwTl7v7YBYtMjQYJf7GdhuWzqK28+qazUiECdZp0oKMnQmtPhISsIivaqLZh3PqpipTrKeereEJNVrKaVJRRD7bIZaWtJ4Ku6xlvp4onVCpK8re2tjCx2CxGLGYjiqLlNLlPrsyfWsLuj9r43w+aWfSpSSiKhqbryX8aaLqOrutomo6mM/D4J8cHj8kDQbl7ICh39cfoDcY5f4WUy26m1GNlTk0xZUU2SotslHpseF3WMU93nG9EoE6TQZL4q7XzL1hMpus6iqojKyoJRSMha0NfywP/kl+rJGQNs8nAdbN8V3in/BaNK+g6TK/yjLpm9iCr2UiR00o0LmecIEIDnHYx7C3kB0mSKC2y0tYdwVWAgdpgkLjzuip+82Yjj/7kjxm/ntGQnA4oL3Ywt7YkGZA9VkqLbBnfR65m4sxkgSRJmE0SZpOBidCXi8QUjAYDUye5szbEXFpk42RLHFsGKchVTcNkkLClmUNcEHLJ7bDQ0pWfGb1SMbO6iFU316AhEY/LGAwSBil530t+nfx+6GuDhCRx3tcSRoNEsdtKkdMy4XvHoyECtZCWaFzBbEqmBL0421gmnDZTxkOE8YSK12UT2ciEvGLNckWtsSZJEtfP9uXF9qyJqvDGYoRxk5BVdB1qyt1ZDdKQvBn4vXYi8Uz2pusFuWdVuPqVFmjyEyE/iEAtpERRNOKyxtRKT856BR5nstrQaBLSiGxkQj5zOSxISDnLbS9c3USgFq5I03TCMYWachf2HFYDMhkNlHpsREfRq47LKh6H+apOIygULpPRQJnXRjAij3dThAIk7mrCiAZzd1f7XHicGaz0SlGJx4o6ih51QtZENjIhr5UXO3A7zISjo5/eESYmEagLUCSmEEuMzR97MKJQXmyntGhs9nzbLCZcdjPxUczniW1ZQj4zGCQm+92YjAZiGazFGI3kFtLCS2UqJIlAXWASsprc8oCBUI6H0UKRZKnK8jFOIFLmtROXUw/Ug9nIJlpaQaHwmIwGaivcqFoyq99Y0HWdYEQhHFXG7D2F7BKBuoDouk4srjLZ72J6VTJlZyCUyMkClWhMxWYxUe13Yhjj7U4uuxmz0YCaYg8gnlDTzjEuCOPFajFSW+kmGldzXrBjMEiXeCxMm+QhGlPFgrYCJAJ1AQlFFXzFdpw2MyajgcnlLspL7QTCiawOayV7szo1Fa5xWZxlkCR8XjuRWGqf/kU2MqHQOG1mqv0uQhElZ4FT13UCYZkSj5VJZS7cDgv+EjvBiJgjLzQiUBeIhKxiMRnxF9uHHjNIEuXFTqZWJj8pj2Ze92KKoiHLGrWVnnEdSi5yWdDRr3gTE9nIhEJV4rZRXmwnGM7+FJau6wTDMmVFNiaVfTIq5i924LSZiMREsC4kIlAXAF3XiSZUJvudw/ZwPU4rM6qL0CGjFaWqlixXWVvpyek2rFSYTUaK3dYLqugMJ5mNzCqykQkFqbzEgddtzWqwHgzSpRcFaUh+uJ/sd6HryQ/lQmEQgboAhCIy5V77iBWq7FYTM6qKcNhMBMJy2sNpuq4TiihU+1248mQYudRjv+KQvshGJhQySZKo8jmxW7PTyx3qSXvtTCpzDvsB1mI2MtnvJBzL3bC7kF0iUOe5uKxiMZvwnTfkfTkmo4GaCjd+r41AWE55MRZAMCxTXmzPq9KbdqsRu9VE4jIrwAfLYopsZEIhMxoM1FQkaz1f7lpPxWDOgzKvncpSx4ijTB6nFZ/XnvOdI0J2iECdx3RdJ55IrvJOdVGXQZKoKHVSU+EmElNS+sMPRmS8buuYb8O6EmlgUVksMfwHjris4hbZyISrgNlkZGqlh7isjmph6GCQ9qcQpAdVlDiwWU1jvqdbSJ+4w+WxYETGX2zHYUt/vtjrsjKj2ouq6iMOqUViCg6riSrf8MNk4y0ZiBk2/7fIRiZcTexWEzXlbsJRJa1894Oru/1eOxUlqQVpGEzA4kJW9LRG34SxJwJ1noonkvuY/d7R93LtVhMzqr1YzUYC4Uv3W8cSCgaDxJRyd972So0GA6VFNqLDbNWSENuyhKuLx2mlqsxJKJLaOpPBOeny4vSC9KDBXAnBqJivzmf5eXee4DRNJy4nh7wzLbJuNhmYWumhrMienLceSLAgKxqqCrUV2S9ZmW3Fbhuqrl1wI5EVFZvIRiZchUqLbJQWXbmAR7InncBfbKd8FEF6ULHbRlmRVeQgz2P5fYeeoEJRmYoSR9a2SBkMEpPKnEwudxGOKsTiCtG4Qm2lG5tlfLdhpcJqNuJxWi7YJx5PaCIbmXBVkiSJylInLruZyGWC52CQLi9xZBSkB1WUODGZDGml7i0kCVklIRduVracBuqtW7eycuVKli1bxvPPP3/J8SNHjnD//fezfPlyvv3tb6MoF16Uhw8fZv78+UPfJxIJnnrqKdauXcuqVav44x//mMvmj4t4QsVuNVFWdOVV3ukqcduYUVWETnJuyjnCdq98U1ZkJ6F88kcmspEJV7PBKSmDUbqkAI+u6/SHZSpKHZQXZx6kYWDHSLmbeEIdVT34fDX4gQYkTEYj4ZhCIJIgEJYJx+SCCd45C9Tt7e08++yzvPDCC2zZsoUXX3yRxsbGC57zyCOP8Pjjj/P666+j6zqbNm0aOhaNRnnyySeR5U+Gf37+85/T29vLyy+/zA9/+EMee+yxgjjJqRoc8q72ZT7kfTkOm5nZU7wUu/NnG1YqnDYTFpOEomgiG5kwIZiMBqZWeFDV5FQVJO8R/WGZSaUOyouzuwDUbjVRVeYiGL06tmwlZJVgOLkgd0ZVEdMmeZhbW8Ls6mJqK92UeeyYjEZC0WTwDkbyN3jnLFDv3r2bxYsX4/V6cTgcLF++nG3btg0db25uJhaLsXDhQgDuu+++C44//fTTrF+//oLXfO2119iwYQOSJDFz5kz+/d//Pe9OaCZCUYWK0uwNeV/OWBfZyAZJkvAXO4jGFZGNTJgwrBYjNRVuojEFRdUIRpJB2l+cm62UJR4rXmdhz1cnkzfJqKrO9OoiKkqcQx0fgyRhtRjxOJJVAadN8jBvajJ411Qkg7fRYCAUVQhG5KHgPd5hJmcRoaOjA5/PN/S93+/nwIEDlz3u8/lob28H4M033yQWi7FixYoLXvP06dO8//77fOtb38JsNrNx40ZmzJiRVrtKS12j+XUuYIsk6An3Uux1Zvxag6JxBZfbzuwpxTnrTecTn8+d9s94i52EZQ1F05lWU4LbMTG3Zo3m3AlJhXjufIDbY+dkcz/zq4upKMv8HjaS4hInh091YzIasZovHLXK5j0vFxRFIxRNMHWymyr/6BfKappOQlaJyyrhqExcVin3ezK6N2dy7eUsUA/X0z2/B3S5452dnTz33HP88pe/vOS4qqq0tbXx0ksvcfToUR5++GFee+013O7UT0B3dyjjOZjQwNBQb184o9cZpGk64ajMjGov3d2hrLxmPvP53HR2Bkf1s0ZdpzcQIxyMEgvHs9yy/JfJuZvoCv3cVRRZMer6mPwORTYTjc39uO3moeBU7HVm7Z6XC5FYcv95td+F02ygrzc7bTUDZrMho3vz4LVnMEij6izmbOi7vLycrq6uoe87Ojrw+/2XPd7Z2Ynf72fnzp309fXxpS99ifr6egDq6+sJhUKUlZWxatUqJEmirq6OiooKTp06latfYcyEogqVpc5xL4RRCEo9NipK7Hm771sQcmUs7w9Om5nKUgehAhgC1zSdQCiBzWJk1mQv3qtwN0jO7na33HILe/bsoaenh2g0yvbt21myZMnQ8aqqKqxWK/v27QNg8+bNLFmyhM997nPs2LGDLVu2sGXLFgC2bNmCy+Vi6dKlvPrqqwCcPXuW1tZWpk6dmqtfYUzE4snMYCVFhbW4a7xYLcaczc8JgvCJsiI7LruJSCx/F5fFZZVQVGZSmZPaSg8W89W5wDSnPeqNGzfS0NDA2rVrWb16NQsWLGDDhg0cPHgQgGeeeYbvfe973H333USjURoaGkZ8za9//et0dHSwatUq/uqv/oqnnnoqrWHvfKNpOrKiUe13FuQCL0EQrl6flMSU8q4kpq7ryYQwusTMai9lXvtVfQ+V9Ktp2XQKsjVH3ROW0ZTMhoUCoQSTfM6c7JnOZ4U+VziexLkbPXHuRicQjnOqLciUSV76+yPj3RxkRSMSk/F5HZQXyDRYpnPUYlJ0nETjCk67Oa/KSgqCIFzM47RS7lXoC8U/KYupk0y2rw9+AUjJDpAkSclHpMEjEpIEkpTspRskCYNh8LH0esGRqAISTK/y4ppACY9EoB4lRdWID5SH++SiHLxghy5dpIGD0tATk9e2ompMqyy6qodrBEG4OvhLHLjcdjq7gkN7ijVdh+T/QNfRSU7nDT6uDvxX03R0dFQVFE1DVjTkhIqqJX/mfNLA6xqkTwK5wZD8bzyuUuSyMqnMmff1CbJNBOpRsJoNVPvddHcH0YYuWkDX0dDRtfMeA3RdG7iIBx/XqS5zYRWZtQRBKAAGScJpNxPJctphTdNRNX0owJ//taIkg7qi6iiahs9rp3iCJjoSgXoUzCYjPp8LizShpvcFQRCyymCQJkSCp0xNrPEDQRAEQSgwIlALgiAIQh4TgVoQBEEQ8pgI1IIgCIKQx0SgFgRBEIQ8JgK1IAiCIOQxEagFQRAEIY+JQC0IgiAIeUwEakEQBEHIYxMuM1k2s+CIjDqjJ87d6IlzN3ri3GVGnL/RyyQL24QrcykIgiAIhUQMfQuCIAhCHhOBWhAEQRDymAjUgiAIgpDHRKAWBEEQhDwmArUgCIIg5DERqAVBEAQhj4lALQiCIAh5TARqQRAEQchjIlALgiAIQh4TgTpNW7duZeXKlSxbtoznn39+vJtTUBoaGli1ahX19fXU19ezf//+8W5S3guFQqxevZpz584BsHv3btasWcNdd93Fs88+O86ty38Xn7/HHnuMu+66a+gafOONN8a5hfnpJz/5CatWrWLVqlV8//vfB8S1l6rhzl3G150upKytrU1funSp3tvbq4fDYX3NmjX68ePHx7tZBUHTNP3WW2/VZVke76YUjD//+c/66tWr9Xnz5ulnz57Vo9Gofscdd+hnzpzRZVnWH3roIX3nzp3j3cy8dfH503VdX716td7e3j7OLctvu3bt0j//+c/r8XhcTyQSekNDg75161Zx7aVguHO3ffv2jK870aNOw+7du1m8eDFerxeHw8Hy5cvZtm3beDerIJw8eRJJktiwYQP33HMPv/rVr8a7SXlv06ZNPPHEE/j9fgAOHDhATU0NkydPxmQysWbNGnH9jeDi8xeJRGhpaeHxxx9nzZo1/OhHP0LTtHFuZf7x+Xx885vfxGKxYDabmT59Ok1NTeLaS8Fw566lpSXj604E6jR0dHTg8/mGvvf7/bS3t49jiwpHIBDg5ptv5qc//Sm//OUv+c1vfsOuXbvGu1l57R//8R+54YYbhr4X1196Lj5/3d3dLF68mO9+97ts2rSJvXv38j//8z/j2ML8NHPmTBYuXAhAU1MTr776KpIkiWsvBcOdu9tvvz3j604E6jTowxQakyRR9i0V1157Ld///vdxOByUlJTwwAMP8Pbbb493swqKuP4yM3nyZH76059SWlqK3W7nL/7iL8Q1OILjx4/z0EMP8eijjzJlypRLjotr7/LOP3fTpk3L+LoTgToN5eXldHV1DX3f0dExNKwmjGzv3r3s2bNn6Htd1zGZJlw59IyI6y8zR48e5fXXXx/6XlyDl7dv3z7Wr1/P1772Ne69915x7aXh4nOXjetOBOo03HLLLezZs4eenh6i0Sjbt29nyZIl492sghAMBvn+979PPB4nFArx8ssvs2zZsvFuVkG55pprOHXqFKdPn0ZVVV555RVx/aVB13W++93v0t/fjyzLvPjii+IaHEZrayt/8zd/wzPPPMOqVasAce2larhzl43rTnycTEN5eTkbN26koaEBWZZ54IEHWLBgwXg3qyAsXbqU/fv3s3btWjRNY926dVx77bXj3ayCYrVaefrpp/m7v/s74vE4d9xxBytWrBjvZhWMuro6vvKVr/DFL34RRVG46667WL169Xg3K+/84he/IB6P8/TTTw899oUvfEFceym43LnL9LqT9OEmvgRBEARByAti6FsQBEEQ8pgI1IIgCIKQx0SgFgRBEIQ8JgK1IAiCIOQxEagFQRAEIY+JQC0IgiAIeUwEakEoYA899BA9PT1s2LCBxsbGCx4DuPPOOzl48OCoXz/Tnx9JNtspCFczkfBEEArYYGGTf/3Xf73ksXxXKO0UhPEmArUgFKjHHnsMgAcffJDGxkY2bdrECy+8MPTYv/zLv1zw/LfeeovnnnsOWZax2Ww8+uijaWWHu9zP//jHP6a5uZnOzk6am5spKSnh2Wefpby8nAMHDvCd73wHWZaZMmUKLS0tfPOb32Tz5s2XtPPFF1/kiSeeoKenh/r6ejZu3JiN0yQIhS/zUtmCIIyXWbNm6d3d3frSpUv1AwcOXPCYrutDj586dUpfvXq13tPTo+u6rh87dky/9dZb9XA4POLrp/LzP/rRj/TPfOYzejAY1HVd1//yL/9S/+d//mddlmV9yZIl+s6dO3Vd1/U9e/bos2fP1t99991h2/nkk0/quq7rHR0d+vz58/WWlpZsnipBKFiiRy0IE8CuXbvo6Ohg/fr1Q49JksSZM2eoq6vL6OcBFi1ahMvlAmDu3Ln09/dz7NgxAO644w4AFi9ezMyZMy/7HoP5j30+H2VlZXR3d1NZWZnW7ykIVyMRqAVhAtA0jZtvvpkf/vCHQ4+1tramXKpwpJ9/4403sNlsQ49LkoSu6xiNxktqaBuNxsu+x/ml/wZfQxAEsepbEAqa0WhEUZQrPrZ48WJ27drFiRMnAHj77be55557iMfjKb3PaH5++vTpWCwW3nnnHQAOHDjAsWPHkCTpsu0UBOFSokctCAVs2bJlrFu3jnA4fMljP/vZz4YemzlzJk8++SRf/epXhwrXP/fcczgcjpTeZzQ/bzKZ+PGPf8wTTzzBD37wA2praykrKxvqfQ/XTkEQLiXKXAqCkDP/9E//xJe//GXKyspobW2lvr6eHTt24PF4xrtpglAwRI9aECaw3/3ud/ziF78Y9tiaNWt4+OGHM3r9qqoq1q9fj8lkQtd1nnrqKRGkBSFNokctCIIgCHlMLCYTBEEQhDwmArUgCIIg5DERqAVBEAQhj4lALQiCIAh5TARqQRAEQchj/x+WZ0jeYwpNRgAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 540x360 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"sns.lineplot(data=df, x='title_length', y='click_through_rate')\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"You can see now the click-through rate depends on the title length! From the data-generating process, we know the dependence isn't causal: it was introduced because one quirky author happened to like 12-16 word titles, and happened to be really good at writing titles. If we were to control for author, this effect would go away, and we'd recover the underlying causal relationship! Let's use the do-sampler to do the adjustment."
]
},
{
"cell_type": "code",
"execution_count": 19,
"metadata": {
"scrolled": true
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"WARNING:dowhy.do_why:Causal Graph not provided. DoWhy will construct a graph based on data inputs.\n",
"INFO:dowhy.do_why:Model to find the causal effect of treatment ['title_length'] on outcome ['click_through_rate']\n",
"INFO:dowhy.causal_identifier:Common causes of treatment and outcome:['U', 'author']\n",
"WARNING:dowhy.causal_identifier:There are unobserved common causes. Causal effect cannot be identified.\n"
]
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"WARN: Do you want to continue by ignoring these unobserved confounders? [y/n] \n",
"Please respond with 'y' or 'n'\n",
"WARN: Do you want to continue by ignoring these unobserved confounders? [y/n] y\n"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"INFO:dowhy.causal_identifier:Instrumental variables for treatment and outcome:[]\n",
"INFO:dowhy.do_sampler:Using WeightingSampler for do sampling.\n",
"INFO:dowhy.do_sampler:Caution: do samplers assume iid data.\n"
]
}
],
"source": [
"causal_df = df.causal.do('title_length', \n",
" method='weighting', \n",
" variable_types={'title_length': 'd', \n",
" 'click_through_rate': 'c', \n",
" 'author': 'd'},\n",
" outcome='click_through_rate',\n",
" common_causes=['author'])\n"
]
},
{
"cell_type": "code",
"execution_count": 20,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x7fe3356597b8>"
]
},
"execution_count": 20,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAeoAAAFFCAYAAAAjJ64BAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzs3XmUVOWdP/733Wqv3qsautn3RaBFAogKYiRE1ACiBs0EmZzgzzmZmHgSjQkyUU8SjUOCk1Xna5Y5RmbAJIK4AJq4AoIgCEg30ND03rXv67117++PokuWrurabtPL53VOTuy6tTx1qarPfZ7n83weRlEUBYQQQgjpl9gr3QBCCCGEpEeBmhBCCOnHKFATQggh/RgFakIIIaQfo0BNCCGE9GMUqAkhhJB+jAI1IYQQ0o9RoCaEEEL6MQrUhBBCSD9GgZoQQgjpxyhQE0IIIf0YBWpCCCGkH+OvdAP6kscTgiwXvgdJZaUJLlewCC0amuj85Y/OXf7o3OWPzl1hus8fyzIoLzfm/PghFahlWSlKoO5+LpI/On/5o3OXPzp3+aNzV5hCzh8NfRNCCCH9GAVqQgghpB+jQE0IIYT0YxSoCSGEkH6MAjUhhBDSj1GgJoQQQvoxCtSEEEJIP0aBmhBCCOnHKFATQggh/RgF6jwkZAWKQlV6CCGEqI8CdR6cnjACYfFKN4MQQsgQQIE6DwlFQSQmXelmEEIIGQIoUOcpTIGaEEJIH6BAnSfqURNCCOkLqm5zuWPHDvz+97+HKIpYu3Ytvva1r110vL6+Ho899hiCwSDmzJmDJ554Aj6fD9/4xjdS9wkEAvB4PDh8+DD8fj++//3vo7W1FRUVFXj22WdhsVjUfAtpxaUEpIQMnqNrHUIIIepRLcrYbDZs2rQJmzdvxvbt27FlyxY0NjZedJ+HH34YGzZswK5du6AoCrZu3YrKykps374d27dvxyuvvILa2lo8+eSTAIBnn30Wc+bMwZtvvom77roLP/3pT9Vqfq8kSUFclK/Y6xNCCBkaVAvUe/fuxfz581FWVgaDwYClS5di586dqePt7e2IRqOoq6sDANxxxx0XHQeAv/3tb9Dr9bj99tsBAO+++27qv2+77Ta8//77EMUrk32dkBWIUuKKvDYhhJChQ7VAbbfbLxqWtlqtsNlsaY9bLJaLjicSCfz+97/H9773vR4fw/M8TCYT3G63Wm8hI44FInGapyaEEKIu1eaoeyoIwjBM1sc/+OADjB07FpMnT874Oiyb/bVGZaUp6/tm0uEMorLCBK1OC4vFXJTnHGrovOWPzl3+6Nzlj85dYQo5f6oF6urqahw8eDD1t91uh9Vqvei40+lM/e1wOC46/vbbb2PZsmUXPafVaoXT6cSwYcMgSRKCwSDKysqybpPLFYQsF6GiGMMgFIoiHIqiTM9ddIFBemexmOFwBK50MwYkOnf5o3OXPzp3hek+fyzL5NVhVG3oe8GCBdi3bx/cbjcikQh2796NhQsXpo7X1tZCq9Xi0KFDAIBt27ZddPzIkSOYM2fORc+5aNEibNu2DQDwxhtvYM6cORAEQa23kBHLALICSAlKKCOEEKIe1QJ1dXU1HnroIaxZswYrVqzAbbfdhpkzZ2LdunU4duwYAGDjxo146qmncMsttyASiWDNmjWpx7e2tmLYsGEXPed3vvMdHDlyBLfeeis2b96M//iP/1Cr+VlRFECUKFATQghRD6MMod0lijX0LTIMTjc5kUgAI6xGlJt1RWjd0EHDaPmjc5c/Onf5o3NXmH479D0U8DxDFcoIIYSoigJ1AXiOpZrfhBBCVEWBugA8xyASS9De1IQQQlRDgboA3cuyKKGMEEKIWihQF0pREKdATQghRCUUqAvEMAziItX8JoQQog4K1AVKzlNTQhkhhBB1UKAuEM+xCEUpUBNCCFEHBeoCcRyDmJiATJnfhBBCVECBukCpzG+REsoIIYQUHwXqYlAAkTbnIIQQogIK1EXAMEA0TvPUhBBCio8CdRHwPIMIJZQRQghRAQXqIhA4FiFaokUIIUQFFKiLgONYiJKMhEzz1IQQQoqLAnURUc1vQgghxUaBuoio5jchhJBio0BdJCxlfhNCCFEBBeoi4XmWan4TQggpOgrURcJzDMK0RIsQQkiRUaAuEo5lIUkKJKpQRgghpIgoUBeVQpnfhBBCiooCdREpoCVahBBCiosCdRFxHEMJZYQQQoqKAnURCRyLMAVqQgghRUSBuoh4CtSEEEKKjAJ1EbEsAzkhU+Y3IYSQoqFAXWwMg7hIgZoQQkhxUKAuNgUQpcSVbgUhhJBBggJ1kXFUoYwQQkgRUaAuMoFnEIlTj5oQQkhxUKAuMp5jEYmJUBTlSjeFEELIIECBusgYhoGsgDK/CSGEFAUFapXEqZQoIYSQIqBArQIGQFykeWpCCCGFUzVQ79ixA8uWLcOSJUvw0ksvXXa8vr4eq1atwtKlS7F+/XpIUjJb2m634/7778eKFSuwevVqtLW1AQA+/vhjzJs3D8uXL8fy5cvxwx/+UM3m541qfhNCCCkW1QK1zWbDpk2bsHnzZmzfvh1btmxBY2PjRfd5+OGHsWHDBuzatQuKomDr1q0AgEceeQSLFy/Gtm3bsHz5cmzcuBEAcOzYMXzjG9/A9u3bsX37djz11FNqNb8gVEqUEEJIsagWqPfu3Yv58+ejrKwMBoMBS5cuxc6dO1PH29vbEY1GUVdXBwC44447sHPnTrjdbjQ0NGD16tUAgFWrVuG73/0ugGSg3rNnD1asWIEHHngAnZ2dajW/IDzHIBJLQKbMb0IIIQVSLVDb7XZYLJbU31arFTabLe1xi8UCm82G1tZW1NTU4Gc/+xm+8pWv4MEHH4QgCAAAs9mMNWvWYNu2bVi0aBEeeughtZpfEIZhAEWBRAllhBBCCsSr9cQ9rSNmGKbX45Ik4cSJE/j2t7+N9evX4+WXX8ajjz6KF198EU8++WTqvvfccw9+8YtfIBAIwGw2Z9WmykpTHu/kch3OIEpLjTDphbT3YQUeJWUGmA2aorzmYGOxZPdvRi5H5y5/dO7yR+euMIWcP9UCdXV1NQ4ePJj62263w2q1XnTc6XSm/nY4HLBarbBYLDAajVi8eDEA4LbbbsNPfvITyLKM559/Hvfffz84jvv8DfDZvwWXKwhZLsJwNMPA5wtBjKUP1MGwiPZOFpUlusJfb5CxWMxwOAJXuhkDEp27/NG5yx+du8J0nz+WZfLqMKo29L1gwQLs27cPbrcbkUgEu3fvxsKFC1PHa2trodVqcejQIQDAtm3bsHDhQowaNQrV1dV47733AADvvPMOpk+fDpZl8dZbb2HXrl2p+8+aNQt6vV6tt1AQnmp+E0IIKQJGUbHW5Y4dO/D8889DFEXceeedWLduHdatW4cHH3wQM2bMQENDAx577DGEQiFMmzYNTz31FDQaDc6ePYsf//jH8Hg8MJlMePrppzFmzBicPn0aGzZsQCAQQEVFBZ555hkMHz486/YUq0ctMgxONzlh0KXvUUsJGbIMTBpZVvDrDTZ0dZ4/Onf5o3OXPzp3hSm0R61qoO5v+jJQK4qCYETC9LEVYC+Ymyf0pS8Enbv80bnLH527wvTboe+hrjtxThQp85sQQkj+KFCrSQHiEpUSJYQQkj8K1CpiWCBGNb8JIYQUgAK1iniOQYQyvwkhhBSAArWKBI5FiGp+E0IIKQAFahVxHAtRkpGQKaGMEEJIfihQq0wBIFLNb0IIIXmiQK0yBkCcAjUhhJA8UaBWGcsA0TjNUxNCCMkPBWqV8TyLcIQCNSGEkPxQoFaZwLEIU+Y3IYSQPFGgVhnLMkgkFEgJmqcmhBCSOwrUfUKhzG9CCCF5oUDdRyjzmxBCSD4oUPcBjmMQpXlqQggheaBA3Qd4SigjhBCSJwrUfYACNSGEkHxRoO4DLMtATsiUUEYIISRnFKj7iAKGAjUhhJCcUaDuIwwAUUpc6WYQQggZYLIO1F1dXXjvvfeQSCTQ2dmpZpsGJY5jEI7SPDUhhJDcZBWo3333XaxevRpPPPEEXC4Xli1bhrffflvttg0qAs8gHBOvdDMIIYQMMFkF6t/85jfYunUrSkpKYLVasXnzZvzqV79Su22DCs+xiMQTUBTlSjeFEELIAJJVoJZlGVarNfX31KlTwTCMao0ajBiGgaKAan4TQgjJSVaBWq/Xo6OjIxWcDx48CK1Wq2rDBisqJUoIISQXfDZ3+t73vodvfOMbcDgc+OpXv4pz587h17/+tdptG3QYALF4AkadcKWbQgghZIDIKlDPnj0bW7duxeHDhyHLMmbNmoWKigq12zbocByDaJwyvwkhhGQvq6Hvb37zmygpKcGiRYuwePFiVFRU4O6771a7bYMOlRIlhBCSq4w96gcffBBNTU1obW3F7bffnrpdkiSwLNVKyRXPMQhGJMiKApaS8QghhGQhY6B+5JFH0N7ejg0bNmDDhg2p2zmOw8SJE1Vv3GCTTMZTIEoytAJ3pZtDCCFkAMgYqEeMGIERI0Zg586dl/Wgw+Gwqg0bvBgK1IQQQrKWVTLZP//5T/zqV79COByGoiiQZRlerxeHDx9Wu32DkIKYmIBJT5nfhBBCepdVoH7mmWfw3e9+F//7v/+LdevW4e2334bRaFS7bYOSwLEIRyVUllzplhBCCBkIsi54smzZMtTV1UGr1eLxxx/HRx99pHbbBiWeZxGhzG9CCCFZyipQazQaxONxjBo1CvX19WBZFvF4XO22DUocyyAWlyDLVPObEEJI77Ia+v7iF7+I+++/H08//TRWr16NQ4cOoaysTO22DUoMwwDM+YQyDSWUEUIIySyrHvWdd96Jn/3sZxg2bBh+97vfYc6cOVntnrVjxw4sW7YMS5YswUsvvXTZ8fr6eqxatQpLly7F+vXrIUnJIWG73Y77778fK1aswOrVq9HW1gYA8Pv9uP/++3HLLbfga1/7GhwORy7vtd9QlGRCGSGEENKbrAL1v/zLv6CmpgYAMG3aNNx3332orKzM+BibzYZNmzZh8+bN2L59O7Zs2YLGxsaL7vPwww9jw4YN2LVrFxRFwdatWwEk128vXrwY27Ztw/Lly7Fx40YAwLPPPos5c+bgzTffxF133YWf/vSnOb/h/oDnWIQitDc1IYSQ3mUVqGtra/HJJ59AlrPf+Wnv3r2YP38+ysrKYDAYsHTpUuzcuTN1vL29HdFoFHV1dQCAO+64Azt37oTb7UZDQwNWr14NAFi1ahW++93vAgDefffdVIW02267De+//z5EceAFPA3PIkCBmhBCSBaymqM+c+YM7r33XvA8D41GA0VRwDAMPvnkk7SPsdvtsFgsqb+tViuOHj2a9rjFYoHNZkNraytqamrws5/9DPv370dNTU2qKtqFj+F5HiaTCW63G9XV1Vm92cpKU1b3602HM4jSUmNBa6F9wRjKK4zguaFZitViMV/pJgxYdO7yR+cuf3TuClPI+csqUPc0v9ztwIEDmDt37mW3K8rlWc3MBfWt0x2XJAknTpzAt7/9baxfvx4vv/wyHn30Ubz44os9vn4uNcddrmBxsq0ZBj5fCGIs/0DtD4lo7/DCMAS3vLRYzHA4Ale6GQMSnbv80bnLH527wnSfP5Zl8uowZj303dP/AOCpp57q8THV1dVwOp2pv+12O6xWa9rjDocDVqsVFosFRqMRixcvBpAc4u7uiVut1tRjJElCMBgcsNnnLANE4pRQRgghJLOCx1176hkDwIIFC7Bv3z643W5EIhHs3r0bCxcuTB2vra2FVqvFoUOHAADbtm3DwoULMWrUKFRXV+O9994DALzzzjuYPn06AGDRokXYtm0bAOCNN97AnDlzIAgDs0cq8CyCEVqLTgghJLOshr4zYdJs11hdXY2HHnoIa9asgSiKuPPOOzFz5kysW7cODz74IGbMmIGNGzfiscceQygUwrRp07BmzRoAwG9+8xv8+Mc/xn/+53/CZDLh6aefBgB85zvfwaOPPopbb70VZrM5lQ0+EAk8i1BESs33E0IIIT1hlHRd4iytXLkSr7zySrHao6pizVGLDIPTTc6C55cDYRGTR5ZBM8R20qL5rvzRucsfnbv80bkrTJ/MUROVKKDCJ4QQQjJSbY6a9I7lgFCU1lMTQghJr+BAPW/evGK0Y0jS8FShjBBCSGZZJZOdPXsWf/zjH+FyuS7qQT/33HP44Q9/qFrjBjueYxEMi5BlBSxLCWWEEEIul1Wg/v73v49rrrkGS5YsoQzlIuo+lzExAb224AR8Qgghg1BW0UEURaxfv17ttgxRDKJxiQI1IYSQHmU1R11TU4PW1la12zIkCTyDYES60s0ghBDST2Xsxj3wwAMAkuU977zzTsyYMQM8//lDnnvuOXVbNwQIAlUoI4QQkl7GQL106dK+aseQxbEsJEmCKMkQeFrWTggh5GIZA/XKlSsBAB0dHRfdzjAMdDqdeq0aYhQkE8ooUBNCCLlUVhlM99xzD+x2O0wmExiGQSAQAMdxKC8vx3/9139h9uzZardzUGNZBuGYWND+1oQQQganrAL1ggULMG/ePKxYsQIAsGvXLuzZswerV6/Gj3/8Y7z88suqNnKw0/AsgmEJ1oG5YychhBAVZTXW2tDQkArSQHLu+vjx45g2bRpEkSprFUrgWYSiImQqx0oIIeQSWQVqSZJw6tSp1N+nTp2CLMuIxWKQJFpaVKjuwidx2qCDEELIJbKuTPb1r38dEydOhCzLaG5uxsaNG/GrX/0KN998s9ptHBoUBTFRhk6j9sso8AZjKDNpqcocIYQMAFkF6kWLFmHXrl04ePAgOI7D7NmzUVpaihkzZsBkyn1vTXI5nmcRisRRalQ3UkdiCbTYguBYBiVGraqvRQghpHBZBerdu3en/luWZezfvx8A8KUvfUmdVg1BAs8iEFZ/GsEfjkGBgnZnCEa9AI6lJWGEENKfZRWoX3zxxdR/i6KIkydPYu7cuRSoi4jnWISjcUgJGTynTvBUFAUefwxmvYBQVILLH4W1zKDKaxFCCCmOnAM1ALS0tOAXv/iFKg0a2hjExIRqgToaT0CUFeg5Fia9AJsrgjKjFhqBU+X10lEou50QQrKWV0QYNWoUzp49W+y2DHkcwyASU2/4OxCOo3vba5ZlwHEMbO6waq/XE6cvguYuf5++JiGEDGQ5z1ErioLjx49ftDkHKQ5eYBAMi6gq1Rf9uRVFgcsfhV7zee9Zr+XgDsZQUaqDUad+VbRITEKnM4QySYG2TKfayAEhhAwmOQ99MwyDiooKPP3006o1aqgSuGThE0VRir50KhpPQJIU6LWfB0eGYaDTcOhwhjC+thSsisu1ZFlBmyMIQWAhKwpCUUn1DHdCCBkMcgrUkiRBURQIAtWkVgPLMpAVBXFJhrbI88aBSLzHQKwVOPiDcXiDMVSY1dtoxe4JIxKTUGLUQKfh4fZFKVATQkgWshp7dLlcWLduHerq6jBz5kysWbMGNptN7bYNUQxi8eJWKFMUBW5/DFptz//cRj2PTlcYUkIu6ut2C0VF2D0RmA3JCzytwCEYiUOUqBIbIYT0JqtA/eSTT2LWrFnYu3cv9u7dizlz5uDxxx9XuWlDE8cCoUhx66fHxATiGbLJOY6FIitweiNFfV0AkBIyWroC0Gm51HA+wzAAwyAQpjrxhBDSm6wC9blz5/Dv//7vKCkpQXl5OR588EG0tLSo3bYhScNzCBQ5UAfDYq/zz0Y9D7s3gmi8uFnnXe4QErJy2RIwnYaFyx8t6msRkok3EFNt1IgQNWW9KUcsFkv9HYlEqE60SnieRUxMICEX7wfF5Y9Bp8k8580wDASeRacrXLR1zr5gDC5fDEb95akQAs8hEpOKfmFASE/snjDqz7nQ6QzROn4y4GSVTLZs2TKsXbsWd9xxBwDg73//O5YuXapqw4a6WFyGQVf48qVYPIGYJKHE0Hvill7LwxeMIxgRYc7i/pmIUgJtjhCMej7tRR3HsPCH49BpaKkfUY/dG0aHK4zRNWVo6fRCr+NVWQJJiFqy+oX81re+hWHDhuGDDz6ALMu44447cOedd6rdtn4pEI7jw89sGFGp3hedARCNSzDoCg9ggUgcLLIf/TDoOLQ7Q5g0QgDL5jdqoigK2hwhsAwyrpXWaVm4fFFYSvU0QkNUYfeG0ekMo8SQ/Dyb9QLanSHoNDxMelq9QgaGrCLBfffdh//5n//BqlWr1G5Pv+fwRvHyP05j6dwRmDdtmCqvwfMMghERFSWFL5dy+2PQ9jLsfSGB5xAJxeH2R1FVlt/FiNsfQyAcR0kvy684joUYlRCJJYpyUULIhZzeCDqdYZgNn190siwDg5ZHsy2AibWlfV4+l5B8ZDW2GggEEA73banJ/mrscDNGDTPjQL0dsqzOXJeG5xCMiAXPpcXEBGJxCQKf2xC6SS+gyx3Oa/lUNC6hwxXMurfCcQz8oVjvdyQkB05vBO3O0EVBupvAs2AZoMUWVO07TEgxZdWN0ev1WLx4MSZPngyD4fPdlp577jnVGtZfMQyDW64dg+dfOYb6Zg+mj60o+muwLANJViBKckFX/MGICOQw7H3h67MMA5snjBEWc9aPk2UFbfZg8ocwy2FzvYaH2x+FtdyQ91A7IRfKFKS76bU8AqE4Ol0h1FpMfdxCQnKTVaAeqvPR6dRNtqCiRIsPj3Zi2phyVeZXGSXZIy4kULv9UejSFDnpjV7Hwe1LViszZFkH3O4NI3y++li2ui9KwjGJ5gxJwZy+3oN0N5NBgMMXhV7Hq1qVj5BCZRWoV65ciUQigWAwSEsbALAMg7lTrNh5oBWN7X5MHFFa/NfggHBUyjv7OiYmEIknUGLIL/gxDAPt+Trg47KoAx6OirB5Inm9nsCz8ASiFKhJQZy+CNod2QVpIPkZNxt4tNmD0Ak85UmQfiur7tbmzZtx9dVXY/78+bj22mtT/z+UTR1TjlKjBh8e7VTl+TU8i2AknvfjQxERTIHXVFoNh3BMgj+YeQ5ZSshosQeh13B5jS7oNBy8wTgVoyB5c/miOQXpbhzLQqfh0NzlhyjR54/0T1kF6j/84Q/YsmUL6uvrUV9fj4aGBtTX1/f6uB07dmDZsmVYsmQJXnrppcuO19fXY9WqVVi6dCnWr18PSUoWv9i2bRuuv/56LF++HMuXL8emTZsy3n4lcCyDa68ahlZ7EM22QNGfn+dYhGOJvJNdPIEotJrC12EbdDw6eqkDbnOHIRUwn84wDBQlOYJASK5c/ijaHIGcg3Q3jcBBAdDmCECmEUPSD2U11lNaWoqpU6fm9MQ2mw2bNm3C3//+d2g0GqxevRrz5s3DhAkTUvd5+OGH8ZOf/AR1dXX40Y9+hK1bt+Lee+/FsWPH8Oijj+K222676DnT3X6lXD2xEu9/2oE9Rzsxekn2SVfZ6O6ZxsQE9NrchuTiYgKhaG5zxenwHItINAGnN4JhlcbLjvtDMTh9UZQYCxu21grJNdXFaDMZOlz+KNrsAZgNmoKSEQ06Hv6QCLs73OPnnJArKWOXy+v1wuv1oq6uDn/+85/hdDpTt3m93oxPvHfvXsyfPx9lZWUwGAxYunQpdu7cmTre3t6OaDSKuro6AMAdd9yROn7s2DFs27YNX/nKV/D9738fPp8v4+1XisBzmD+tGo3tfnS6VFi+piCvEpuhqFjUBLfuOuCX7uolSgm0OkIw6tJXH8uWRmBpRy2SE3eRgnQ3s4GHzROBr5epHkL6Wsau2vz5888PSyaHg55++unUMYZhMg5/2+12WCyW1N9WqxVHjx5Ne9xisaS2zrRYLLj//vsxc+ZM/PKXv8STTz6JX/ziF2lvz1ZlZXGWYXQ4gygtNcKkF3DzvDHYe7wLBxocWHvrtKI8fzetToSg08CSwxIpAHCHRVRXmXMqdNIbXUREVFZSy7UURcGZNi/KSvUw6XPvBZeXXd5r4TUxCDotLOVU3jGTXD8Pg5HDE4E/lsDImnJwOQTpnj53Fyoxy/DFRNTW6HMeyRrs6HNXmELOX8ZPYkNDAwBAlmWw7MWd79561D1lh1/Y68p0/Le//W3qtm9+85u4+eabM96eLZerSAUOGAY+XwhiLDnce81kC/Yc60JjswuVpcVb5pGQZXh9YZiE7OeaRSmB9k4fzEYNil2jxukOgpVlmA0auP1RtNqDKDVp4InltttXeZkRHm/osttFKYHT55yAVFasJg86FosZDkfxcyIGEncgilZbEGaDAL8/+xGYdJ+7S8XiCRw63onxtSUZS+DmQpQSiEsyjFkudexv6HNXmO7zx7JMXh3GrD6FPZUOvffeezM+prq6Gk6nM/W33W6H1WpNe9zhcMBqtSIQCODPf/5z6nZFUcDzfNrb+4P506rBcwz2HO8q6vNyLAtJknPKRg0WeYvMC+k0HDqcYURiEtqdIZgMxT3/As8hGpMuG2InpFt3kDYZeNUK5Gg1XPKC11HYTluyrMAfjqOp04+GZg/OdvgRF+mzTXKXMVDfd999mD17Nk6ePInZs2en/ldXV4fS0sxrhxcsWIB9+/bB7XYjEolg9+7dWLhwYep4bW0ttFotDh06BCCZ0b1w4UIYDAa88MIL+PTTTwEAf/nLX7BkyZK0t/cHRr2A2ZMsONrogi+U/5KqdGI5fLm9gTg0RRzyvpBG4BATE2jq9EPgGHBscXobF2IYBv5w8c8hGbikhIxQVITDE0kF6Xw+e7kEXZNBgC8Yg8MXyfk1IjEJXa4Q6pvdONfpR0yUYDIIYMHATXuwkzxk7BL99re/hdfrxY9+9CM89dRTnz+I51Pzy+Fw+KKyot2qq6vx0EMPYc2aNRBFEXfeeSdmzpyJdevW4cEHH8SMGTOwceNGPPbYYwiFQpg2bRrWrFkDjuPw7LPP4vHHH0c0GsWYMWPwzDPPpL29v7h2ejUONjjw0fEuLJ03qmjPy7AMollW7RIlGcFoXNXCISY9j2hcgk6rzmvotBycvgiqSnW0o9YQJMsKYmICcTGBYEREICKe74UyYBnkFaQVRcGuA61odYSw5kuTss7dMBkEdDnDMGiFXr9ToiQjEI7D6YsiFk+AZZNlSi/s9euDn/SQAAAgAElEQVR1HJy+KCpL9TnX3ydDG6MUWGps5cqVeOWVV4rVHlUVa45aZBicbnJeVlpz+4dN+KzJg+/cNaNoc1FxMQGB4zC2pqTX+3oDUbTYg/1+iVNvc4X+UBwTasuoUlQPBtNcoaIoiEsyYmIC4aiIUEREOPb56BHPMRB4tuB54v0nbNh1oBUAMHtSFW5bMCbrx4qSjLgoY8KIUmgvqRMgKwrCUQkuXxT+UByAAp2WzxiEAyER1ZV6WMsu79z0Z4Ppc3clFDpHXfAvIZUU/dx1Vw3Dp40u7D9hx02za4vynALPIhQVIStKr2U8PcH4ZT8mA1H3jloUqAcXWVYQiUuIxCQEw2Lycy0rUMCA45LV+Ez6wpf6Xaix3YfdH7diyqgyDLeY8M6hNkwZXY4JtdmV/RV4Nll5zxbAuJoScCyLaFyCNxiD2x+DlJCT7TZk126jnofdHUGFWVe0RDUy+BX8SaHhyc9VlekxdXQ5Pq63Fy0himEYKDJ6TUKREjKCYRGaHDLE+6vuHbWoStTgoSgKWu1BnGn3o8sVRkxKFvIxGzUoMQow6gQIfH4laNNx+qL427tnYS3TY8UNY3HLtWNRVarDjj3nEI1lX59Ar+URjSfQZg+isc2LU61euHwxaDUsSowa6LTZX1ywbHK5q5fWapMcDPxf9X7muhnDEBMTOHjSXrTnVKAgJmbO/A5FJYAZHBdOyR21qKToYOIJxOALxlBiFGAyCNAKnKrbmkZiEv7vH6fBcQy++sUJ0AgcBJ7F8hvGIhgRsfvj1pyez6TnEYwkR7ZKjBoY9fkltAHJKmh2dxgJmWqLk+xQoC6ymiojxtWU4KPPbEUr8i/wLEK9bNDh8UehGUQJKgLPwBOgDNnBIBqX0O4Mwljk5XzpyLKCv713Ft5gHHcvHo8ykzZ1rLbKiOtmDMORRhdOtWauBXEhhmFg1AsFbTvbjeNYSLICX5BWN5DsDJ5f9n7k+pnDEYpKONLo7P3OWRB4FsFI+t6llJARjAyOYe9utKPW4CDLCtrsQQgcq8pyvp7s/rgVZzv8uPXa0RhVfXk1qIWzamAt0+O1vc2I5DAEXkwGHQ+bO1ycAkxk0Mv7m9O901VPS7OGutHVJoywGrH3WFdRhrd4jkUsnkgbtMJRCbIyOIa9uzEMA9COWgOe3RtGOCZB10flOD855cCBejvmTavG1ROrerwPz7FYfsMYhKNSKhu8r/EcCymhwB+iuWrSu6wC9aVbVDY2NuKuu+7q8RhJBpnrZwyHLxTH8bPuoj1vuoQyTyAKjTB4gnQ3jcDCRQUiBqxgRITdHYHZ0DdlM891BfDGvhaMrynBkjkjMt53eKUR188chqNnXDjZ4umT9l1Kp+XQ5YlQ0iTpVVaB+v/+7//w+uuvAwD+9Kc/4d5778Xy5ctVbdhAN3FEKarL9dhzrKsoS9gYBgj3MEwnJWT4w+KAWZaVy7n4fEctGv4eaKSEjFZbAHpdcTO50/EEYnj5nTMoN2uxatG4rBLVbpg5HMMqkkPgV2LkRuBZxMUEAipUMySDS1aB+o9//CN+//vf4+6778auXbuwZcsWrF27VuWmDWwMw+C6GcPh9EVxsiX7pJV0BJ5FMHx5He9wTAIGyLC3wxvBsy8fxdFGR1b37x7+DlBJ0QFFURR0OkOQ5WT9drXFxAS2/KMRiqJg9c0Tsh5m5zgWX7l+LCLxBHbub1G5lT3TaTjYPBGqR9GLoT6Xn9V+1IIgYOPGjejq6sIDDzyA8vLyXnfPIsC0MeUoN2vx4dHOgr+I3YVPLn0ebyAGge//QVpKyHjl/bMIhEW8daAl6/Oh1XA0/D3AeIMxuIMxGPTqB2lFUfDK+01w+CJYdeM4VJbktnvdsAoDFs4ajuNNbtSf6/shcI3AIRpPqLqZzkAnJWQ0tvsQig7dc5TTftQA8MADDwDofT9qklwPfN2MYXhtbzOaOpOVjQp5Lvl8ycXuYe6ELMMfisOozz5RJxqXwLMs+D5eyvXO4XZ0uSOYMqoMDS1etDtCGGHtvZSeRuAQCMURiycK3l87JiZgd4dhMggoNxdvO1LyuZiY3HXKpCtuhbF0/vlJO061evHleSMxvia7amOXum7GMJxs8eL1j5oxapipz7ei1GpY2NxhmPTCgBgZ62udzhDCURHtjhAmjCjttULjYJTx17qhoQH19fUX/f9nn32GY8eOUZDO0szxlTAbBHx4tLMIz8ZctJNWMttbyenL/eKuU/jvHScQjffdnFxTpx/7jtswe1IVVtwwFjoNh48bsi8IwzBMQcPfCVmG3RPGqVYP/OE4Wm1BuGmNdtHJSnIpFscx4PqgPOaxMy7sOdaF2ZOq8IUp1t4fkAbHslh+/VjE4gm8sS/70Z5i0QocwjGpxxyUoc4XjMEdiKHUpEE0nhiyu49l9W3av39/Knns7NmzuPHGG3H48GFVGzZY8ByL+dOrca4rgDZ7sKDn4lggfMHwjzcYy2kXHoc3gk5XGE5fFK+839Qn8z6RmIRtHzShskSLL31hJDQChy9MG4YT5zwIZTnc172jVq4/oN2lGk+2eGHzRGDUCTDqBRgNfDJYD9EvvVqc3ghCUQn6PliK1e4I4tU95zCq2oRb5o0quCdqLddjUV0N6ps9OHFFhsCTvWryOVGS0eYMwXh+dMao59DlCkOUht6e3ln9yv/85z9PbXM5ceJE/Pd///dF216SzK6ZZIFey+HDY10FPY/AcwicTyiTz1c20uUwHFzfnPwBun7mcJxu8+HdI+0Ftac3iqLg9X3NCEUkrFw4LlXV6bqZNUjICj45nV1BGJ5jEU8oiOZQPz0Sk3C2w4+WrgA0PAuzQUhlAnNs8u9We4CCdZGEoyK63GGYcpiGyZc/FMeWf56B2SDg7sXji9Z7X3DVMNRUGfHGR809Jm7my+WP4rW95zIu1dRpeAQiItUNOE9RFHS4QmAUpKbpOJYFyzDocg29C5qsPuGiKGL69Ompv6dPn454nDJxs6UROMydWo1TrV7YPPl/yHiOQTSeQEKWEY5JUHIc9q4/58FIqwmLr67B1ROr8OHRLnzWVLx13pc6esaFE+c8uPHqGtRUGVO3V1cYMHa4GYdOOrLu1XNMchisN6Iko90RxKk2L8REAiUmTY/z8SzLwGzQoNUeoGS1AkkJGS32IHQadet3A4AoJbDln42Iiwms/uLEy7aaLQTLMlhx/RjERRmv72sueAjcH4rjtb3n8LtXjuOTU0688VFzxs1ANDwLhzdS0GsOFr5QHN5gDIZLLvz0Og7uQGzIJd9lFaj1ej3ef//91N/79u2jimQ5mjvVCg3PYk8BverkciUFsbgMbyCWU0KYyxeFzRPBtDHlYBgGt8wfhREWI17dcw5dKgy5uf1RvPlRC0ZVm7DgqmGXHf/CFCv8oXjW9Zb12sw7asmyApcvilOtXngDcZQYBOg0mXt33cG6zR6E00c/kPnqcocgSXJR6mBnoigKXv3wHDpdYaxcOA7Wcn3RX6OqTI+bZtfiZKs372JFoWhy049f/+0YjjS6MGeyFffePBHReAL7PrOlfZxOw8EXil2xsqb9hSglExKNPWxzyzAM9FoO7Y7QkFqyldU41fr16/Gtb30LPJ+8O8uy+PWvf61qwwYbvZbHNZMt+OiEDTfW1aAix2Ukn2MQiYnJq80c9mvuHvaeMrocQHI4+a7F4/HCa/XY8o9GrLt9atF6J7Ks4JUPmsAwDFbeMLbHXtakkWUoMWrwcYM91aZMLtxRy6S/uJ3BSDIjNCYmYNRzOdWUTgZrAe2OEKAkf6gHI0VREIpKiMYkmI2aohXI8QVjcPmSu2Kp7YOjnfjsnAc3za7F5FFlOT1WlhW4/BGIYqLX9z5vWjXqWzx4c38Lxgw3w2zQZPUa0biEjz6zJTfkSciYOb4Si2bVoMyc3BRk2phyfHTChrlTrTDqLz9fDMOA51g4vRGM7KFG+VCgKAraHSEwDNLu160ROPhDcbj90UH7fb1UVr9os2bNwrvvvovf/e53eP7557Fz586LhsJJduZPrwbLMNh3PP1VdW8EgYHbH4OsIKdhxvpmD2otRpQaP//RMRs0uPumCQhGRLz87tmibbv3/qcdaHeEcOu1o1F6wc5FF2JZBtdMtqCpMwBnlsN9As/AG/h8+DsWT6C5K4CzHT4wjIISo5DXxg/dwbrDFcq6LQOFoigIRkScaffhTLsPXe4wTrZ40NwV6HFdfi7iYgJtjhCMenWWYimKgi53GAfq7Xj5nTN493AHZoyrwHUzLh+h6U0wLGJ4pQnReKLX98yyDJZfNxZSQsFre3sfAhelBPYe78Kv/3YM73/aifG1pXhg+XQsv35sKkgDwOKrayElZHx4LP0KEL2WgycYK9p+9gONJxiDPxTvtRNi0gvodIXTllUebDKeje3bt2P58uX405/+dNHte/fuBQD867/+q3otG4TMBg3qJlThSKMTi2fX5tQj7qY5n1CWS7a3JxBDpyvcY/3j2iojbr9uDLZ90IS3Pm7Dl+eNyrlNF2q1B/HB0U7MGFeBq8ZVZLzv7IlVeP9IBw6edGT1uskdtWKwluvhDkTh8ETAcSxKjNn1eDJhWQYmvYB2ZwgyFFjLBvbUTneAtrmTm2JoNRxKTZrUsUhMxNn2GDQCD2u5DiVGTU4XObKioM0RzNjzyVUiIaPTFUazLYAWWxCt9mAqgbDEIOCaSRYsnTsy54uCcFRCiVGDkdUmeL0hOL1RmHsZAags1eGL19Ri14FWHD3jwqwJl2/wkUjIOHzaifc/7UQwImJ8bQkWX117UT7Gpc9ZN6EKBxscmD+tuseLWIZhwHEMXP5o2ucZrOJiAh3nL/x6w7LJ89TpDmP0EBh9yHhGmpubAQCnTp3qk8YMBddMseDQKQeON7kwd2p1zo9nWQYMkFe299Q0Q8wzx1eiyx3GR5/ZUF2ux9WTLDm3C0j2cF95/yxKjRrcMr/3wGvUC5g2phyfNrqweHZtr0OSyeI7SM1rG/VCUZOXunvWnc4woADW8uIGaykhIxKToOG5gou3pNMdoLvOB2i9hrvsQoZhGOi0PHTa80tgHCEwjhCqynQoN+uyGhZ3+6IIhkWUmPK/SBKlZI+8xRZEsy2ANnsotUNcZYkOU0eXY/QwM0ZZTSg1afLqtXcv5am1mMAwDCxlengDMYiS3OvF7typVtQ3e7BzfyvGDi9JnUdZVnC8yY33jnTAE4hhpNWEVYvGYfSw3gPGwlnDcfSMC+9/2onbrxvT430MWh4ufwRVpTrV5/37C0VR0O4Mng/A2V346bUcvIEYKkt0l02HDTYZA/WDDz4IALQUq4iGVRgwvNKAw6ec+MIUa14/Prn+OJ4450FNpeGiYbhL3XzNCNg9Ebz+UQuqyvQYmUXVsEvt3N8CXyiO+748uddErm5fmGLFsbNuHDvjwpwsilYY9BwUpXi9uEt9PgyeTLArNFgne68JeINRuPyxZKY+kiVhy8xamA0a6DV8wRcciqIgEBHR5QojGk9Ap2EvmuZIR+BZCDybnMP1xmD3RFBq1KKqTAeDtuch7UhMQocrBFOOu2JFYhJa7cFUYO50hlPJgcMqDJg9qQqjq80YWW0qyg+vLCsIRxMYX1uaCso8x2KE1YSmTj94LnMlMIZh8JXrxuD5V09gx95zuPfmiTjZ4sU7h9vh8EYxrEKPe26eiAm1JVl/j0tNWsyZbMGBBjsWXDUMlaWX56owDAMWLNz+KIZVXvletZSQ4Qsmk9zUWiPvCcTgD4tZfWa7MQwDg45DmyOISSPKVF9xcCVlPOu33357xgfv2LGjqI0ZKuomVuHNj1rQ5Q5juMpfRG8whg5nCF+8pjbj/ViWwapF4/DCa/V4+Z0z+OZtU3MaUv6syY1Pz7hww6zhGJXDUFStxYjhlQZ83GDHNZMtvf7g5TMHnSuWZVBqTAZr5fwweK4XVKIkwx+Ow+mNIiZK4DkGpgvmcqWEDJc/CrsnCpYBzAYBpUYNDDohp16UrCgIhOKweSKIxiXoNFxeiV0sy8Bo4KEoCsIxEWfaY9AKPKrLdTBfMCyekGW02ILQ5rgU62SLF3999wwSsgKWZVBTZcC1V1VjVLUJIy0mVfarDkYkDK80XBb0zQYNykxaBMPiZct/LlVRosPN14zAm/tb8Ju/H4cnEENlSXKHru4VFLm6buZwfHLaiXcPt2PVjeN7vI9ex8Hpi6KyVJ/TNFexyLKCcEyCJxCFNxhHSYkefl8YtRYTys3aouYkxOLny87msQZf4DlEQnE4/ZEBP12VScYzs2HDBgBAZ2cnnnvuOTzxxBM4c+YMNm7ciCeffLJPGjgYzRhbgbc+bsWR007VA3VDL8PeF9JreXz1pgn44+v12PpOI9Z+eUpWS8B8oThe39eMmiojFs4anlP7GIbBF6ZY8eqec2i2BTEmi+HDvsAwyWDd5QpDkZNrv3v7cZIVBeGoBLc/Cl8wDjA4Hzgvv+DhOTY1KqAoCiJxCb7z2x3qNDzKzBqYdAJ0Wr7H2sbdAbrLHUFMlNK+Tj7vW6/lodcmh41bLxkWd3ojEKVETr1ppy+KVz44C2u5HkvmjECtxaR68IlEJZj0fNqs4OGVBpwK+5BIyL0Otc6ZYsHpNi8c3ihuv24MZo2vLKj3ZtILmDetGh8e7cR1rjCGVV4eYLqf3xOIFn0KJh1ZURCNSfAG4/AEopBlBTzPwqTnUWbSQorH0WYPIhgRUVNlLMqoVneuA88zeV+Em/QCbO4wSo3aAbPdb64yBuq5c+cCANauXYt77rkHc+fOxaxZsxCLxbBt2zbceuutfdLIwUan5TFldDmOnXVjyZyRqm6QcaLZg2EV+qyXg1nL9Vhxw1hsfecMXt/XjK9cPyZjgJJlBds+aEJCVnDHwrF5fdmmj63AWwdbcbDB3m8CNZAMWiVGDeyeZCZ4umAdFxPwh+JweCMQEwoEnoHJkH0mNMMw0Gl46M7HWUmSYfdEYJPDYFgGpUYNSo1a6LU8ZDlZFrXLncx41WuLE6B7IvAcBJ67aFhcUZBTjz0mJrD1n43gORZ33zQhp6HNfImSDBnASKsp7QYOAs9heJUBrfZgr21iGAb33Dwx9d/FsGB6NQ422PHO4fbUc1/KoONh90RQUaJTbapHUZIV//yhOFz+6PkLl+TF2qUXIxzLwmwU4A/FEY5JGF1tLngo3O2PIhTpPdchU3Enlk0G+S5XOKs8gYEoq399j8eDNWvWAAC0Wi3Wrl0LhyO7PYVJz+omVCEaT6ChCHtVp+MPxdFmD2XVm77QlNHlWFRXg0/PuLD/RObNM/Z91oXmrgC+PG9U3mvDBZ5F3YQq1Dd74A/1r4p3DMPAbBRg80TQ5Q6nlurIsgJ/OI6mDj8aWjzocoehEViUGAXo08ztZivZixFgNmpg0PIIRUQ02wJoaHbj09MOtNgC4NjkRURf7PfcPSxuNggoMWa/w1N3gRKXP4pVi8b1SZBWzo9qjK4293puyk1amHRCVgVGGIYp6nCvTstjwVXDcLrNhxZboMf7sCwDBbhoSWKxxMQEnL4ITrX60Njmg9MXhVZgk585XfokTYZhYDIIYBjgdJsXLl807yV+0biEzl5yHbrLEL/wWn3GGt8GHQ9vMDZo967PKlAnEgnYbJ+v/XU6nbTReYHGDjej1KjB4dPqXfA0tJwf9h6TeZlUTxbOGo7Jo8rw1sFWnO3w93ifTlcI7xzuwJTRZaibUFlQW+dMsUJRgE9O9b8LwGTPWoD9/KYmXe4QGprdONfpR1xKwGwQYDIIquwYxbLJDG2zQUj9iCYDdN/PW+YarPYe70J9swc3XzMCY4fnv8VrLgJhEdUV+qyS0RiGQa3FCEmSr0iVq7lTrTDqePzzk/a0v6cGHQebJ1yUGgeilNx9qrHNi5OtXnS5wuBYwGwUYNTzOX1+tQIHk15AmyOIFnswla2freROayEIHJtxGuFAvR2HTjrQ6QrjncMdGZ+zO7GsWPUg+pOs/mXWrl2LFStW4JFHHsEPfvADrFy5Et/85jfVbtugxjAM6iZWoakzoMoVM5DM9raW6VHVQ2ZpbxiGwYobxqKqVIe/vXcGnkvaKEoJ/P29Jhh1PG67NvPweDbKzVpMHFGKT045kcjxS98XGIZBiUGAyxeByxeDTsujxKiBVsP16R7C3ADJbD3T4cM/P2nH9DHlmD8992WI+YjEJBh1Qk5zujoNj+pKwxWpHa0RONwwqwYttmDai2GOZZGQFXiD2fcUZVlBLJ5AIByH0xdBqy2AU61eNLR40O4IQlYUlJy/uCxk2o1lGZSaNAiGRJxu8+a0oYjbF0U4JmVMIjzXFcDuj1sxeVQZZk+qwv4Ttow7EAo8B0mS4fINvtr9Wf0r3XnnnfjTn/6EadOm4aqrrsIf/vCHXjPCSe9mne+FfnrGVfTnDoZFtNiCmDomt2HvC2kFDl+9aQIUBdjyj8aLqgC99XEbXP4oll8/Nq/CLT2ZM8WKYERUdTqgEMlhcA2M+sKXUw1m3kAMf3+vCZZSPW6/rvCLuGwke8WZ56XTqSrRQyvwV6TK1exJVSg1anrpVfOwucOX9foTsoxoXEoGZG8yIDe0ePBZkwunWr041xlAlyuMUEwEyyaz3c1GTdHXZhsNPDiWQWO7N6vtaFPL+zJkeftCcfz13TOoLNFhxfVjsWTOSJgNGry65xwkKf2FvPF8Ytlgq+yW9eXUlClTsHbtWnz961/HpEmT1GzTkFFm0mJcTQmOnC7+VEJq2DvH+elLVZTosOrGcXD4Itj2YRMURcGpVi8OnnRg/vRqjKsp3pDmhNoSlJu1OFCfeV6c9F+iJOPld89AlhXcddP4PinYkaxjLmKk1ZjX67EsgxFWI6Kx3suLFhvPsVhUV4NOVzjtBSrPsZASCpy+CJzeCJptybyIE+c8ON3mQ1NnAF2eZEAWuOTFpNkowGxM9pp1Gl61ZLRumvND4e3OZAGbdEPhspzM8tYI6Ye8JUnGy+80QkrIuHvxeGg1yQJBty0YDacvivc+TT8EzrIMeJ5Fhys0qKZn+36ii1ykbkIVfKE4mjp7TijJ14lmD6pKdbCU5bv5x+fG15RiyZwRaGj2YvfHrXh1zzlUlyd3GSomhmEwZ7IFrfYgbCrs6EXUpSgK3tjXfH53q7GozHvjmdwEwiKs5QaUGNMX9OmNUSegqkyPUKTvd66aOb4SVaU6vPNJe9q5coOOQ5c7jC5PGNFYIhmQDUIqwc+kTwZkNfIkssWeX6EQjIg43epFOHr5dILTF0EkJqUtiKQoCt74qBkdzjBW3jD2ouV1E2pLUTehEnuPd6HDGUrbDr2Whz8cH1SJZRSor7Apo8qg03A4ctpZtOcMRUU0dwUwdXR+BRl6Mm9aNWaOr8T+E3bExQRWLhynylV63cQq8ByLjxuoVz3QHDrpwKdnXFg4azgmjcxtd6t8RWISDFoB1UVYa2wt1yd3acswtKoGlmVw49U1cPqiOHa252kw/nxNe5NegFbDXdGA3JtkYhqTzCb3fj4UHo5K6HKHMyb6HTrpwJFG1/lk1stHA7/0hZEw6gS8uudcxlwWg5ZHuyOUc5Jbf9V//7WHCJ5nMWNcBeqbPUXbh/ZkixeKgoLmpy/FMAxuvXY0rhpXgduvG6PKXsBA8mp4xrgKHDvrRnSI78s7kLTag9h5oBUTR5RiUV1Nn7ymlJCRSCjJeeki5AzwHIvaKiNCOSRF5SJTZvnU0eUYXmnAe0c6+mUyZa40AgeTQUCHM4TmrgBiYgJtjswV7VpsAezcn/kzpNPyuHXBaNg9EXx4rCvt6ws8C0lWBs0+8xSo+4G6iRYkzhf6L4YT5zyoMGtRXeRgKvAs7lg4DjPGFbYUqzdzplghSjKONBY3yc4TiGH/Cdug+CHsTwLhOF5+5wxKjRqsvGFsnySPda+XHmk1FXWDkxJjsrxoT8O2+ZIkGf5gHL6wmHbpEMMwWHx1LbzBOD4p4ujalcSyDEpMGoSjEk63ehGNS2krh/lDyc9Qmbn3z9DkkWWYMa4CH3zamXGKzKRPFoyJxgf+BT8F6n5geKUBwyr0RRn+DkclNHX6MTXPOsT9wfBKA0ZYjDh40l60hJBmWwAvvFaPXQda8coHTVdk3exglEjI+Ou7ZxETE7j7pvGq1OzuSTAioapUn3a/83wxDIPhlQbIMgpejyvLyTKvoiSndgELR9JnI4+vLcGoahM++LQzY3GPgcag56HTcjCnKWwiJZIJiHFJxt2LJ2T1GVo6dxT0Wg6v7jmX9rvMMAwEnkWHc+AnllGg7ifqJladL6ZRWBLVqdbzw94FZntfaV+YaoXbH0u7vjQXR0478eKuUzCcrwZ14pwHr+1rHvBf3lwkEjKOnHaisd1X1Pf91sE2tNqDuH3BmKLME2eje/OR6gp1pl80AofhVUaEwvn1xLqz0EMRCdWVBkwaVYZSkxblZi20Gi7tMjCGYXDT7FoEI+KgW/nAc2zajsOu/a1od4Sw4vqxWU+pGXQ8ls0fjU5XGHuPpx8C12t5BMNiqo7+QKVqoN6xYweWLVuGJUuW4KWXXrrseH19PVatWoWlS5di/fr1kKTkF2Pbtm24/vrrsXz5cixfvhybNm0CAHR0dOBrX/savvzlL+Pf/u3fEAqlz/wbaK4aWwmOZQruVZ9o9qDMpMHwHgr9DyRTR5fDqOPxcUP+lcoURcE/DrXh1T3nMLrahG/cOgU3zxmBhbOG48hpJ3YfaB0SwbrDGcILr9Xj1T3nsPmt0/jN349j3/GunApU9OToGRcO1Nsxf3o1rhqXe/W7fCQSMkQxOS+t5k5q5WYtDDoh52HTmJiAPyyiRK/B5FFlsJYZUu1kz1dCi8SktJ+7UdVmTBhRij3Huq54joYoJXCy1YvX9p7Dr7YexqeNxV9G+skpBw6dcuC6GcNyzqmZOqYcU0eX470jHXB40xRB/0gAACAASURBVM9FG3Q82uxBeAPR1LaqA41qn3SbzYZNmzZh8+bN2L59O7Zs2YLGxsaL7vPwww9jw4YN2LVrFxRFwdatWwEAx44dw6OPPort27dj+/bteOihhwAATzzxBO69917s3LkTV111FX73u9+p1fw+Z9DxmDyqDMfOuvLOVIzGJJzt8Bc12/tK4TkWV0+y4FSrN6/KbXExgZffOYM9x7pwzWQL7l0yMbWBwKK6GsybZsX+ejveO5K5LOFAJkky/nGoDX94vR6hqIS7F4/HHQvHwqQX8NbBNjz78qfY/mET2jMsdUmn0xXGa3vPYfQwM26+ZoQKrb9ccs9tCSOsxqz3O89Xd1CNi9mVF00kkvPQLMNgQm0pRlabe1zTbdQJqDDrEImlH9pefHUtovEE9n1mS3sftfhDcRw66cD/vn0a//m/R7DlH404ftaNYETE9g/P4f/tOFGUUS4AaLMH8eZHLRhfU4LFV+e31POW+aOgETIPgfM8C52GQ4stiNOtPgTC8QF3gc49/vjjj6vxxG+99RZYlsWtt94KQRDgdrvR2NiY2pGrvb0df/3rX/HII48AAEpLS/HSSy9h5cqV+O1vf4umpia88MILOHToEObNmweO4/DEE0/g5z//OViWxciRI/HMM8/gvvvuy7pNkUgcxfj3kRkGbm+46BsiaAUOn5xyorpCD0ua7fkyOdHsQX2zB1+aO1K1HZWKQa/TIJpFsk5FiQ4H6m3gWDanwir+UBx/2X0azbYAls4diRuvrgF7Qe+LYRiMrymBPyxi/wk7BJ7FSKspr/fS17I9d632IDa/dRonW72om1iF1TdNwLBKI6zlBlw9sQpTRpVBVoDjTW4cbHDgdKsXLMugslTba081HJXwl90nwXMsvv6lSUVN5uqJoiiIizKCURFVpTpUl+e3NazRqEU4h7W1As9CgQJvIJb2PSqKglBEQkJJ7q1eU2XsdatFvZaHwxuFwPc8HGw2CHB4Izh6xoXZE6tULRqjKAo6XWEcOunAWwfb8PbBNpxu80GWFVw1rhI3zq7Bsvmj8aX5Y2HSsjjd5sOBejvaHMGsa6r3JBgW8eLuU9BreXztS5Pyfo8aIbn/+oF6O7QaLu33mGUZaDUcZFmGwxtFKJqcPumLTW2Azz97DMPAYMj9t1m1y1K73Q6LxZL622q14ujRo2mPWyyW1MYfFosF999/P2bOnIlf/vKXePLJJ/GDH/wAJpMJPM9fdv9sVVYW58e4wxlEaakx7w9pOrNLDHj9oxYcb/LgurqROT++sb0JZWYtrppg7fc96vKy3n9sy8uMmDG+Cp82OrHixolZbUTRYgvgj280ICYmsO4rV2Hq2PQZ6l+/ZRoU1OPtg20oL9VjwYy+WVZUqEznLiYm8MaeJnxwpB1lZi3+v5UzMGX05cPS5WVGTBlnQTQm4f9v776Do7zv/IG/n7L77D7btU0FFSQECh2bOJQAITElBizHIRfHyRHOCZffTMaZcW5u4kyS8Y3vkvN5XC7x+XwpnvGNywXOmZiQcwHHAROMz4ABcaFIFAEqoF62PfuU7++PlTaACtuEdsXnNeMZa1e7++jLs8/n+bbP5/Dpq/jT8Tb87kAz3j3cgrvmFGPZvBL4R5lzNgyG7e+dwGBUxcNfWohpxRNTbCOu6ojFdRiMQeAAb5EEt0OC2y5ltYfY70+vDGJRkQ0nL/QkLvQ3BJNwTIWqGagtdiNYJKdVKEWwmNHaMQjXGEla6lfOwBMvH8Khxi58YdWMtI75ZuKqjsbLvTh5vgcnL3SjPxwHB6CyxIkNy6djTrUXxaOUdF1xRwWWzp+G/cdbseejS/jF707irtnFWL+0Cu40FvVpuoGXdzdCiev4f1+ej1J/dtflFYtkNLYMYO/Hbfjk7NHP2xtFFQ2doTg8PI8yn/2WLIJM99y71oQd3WhDC9f+w4/3/PPPP5987Jvf/CbuvvvuZM97rPdLRXd3KDerfTkO/f1hqEpuAzWA5LaD5tbetMoCKnEdp5p7sLjOj77+/M7q5XHb0NuX2nDrgpoiNJztwoFjl7Fghm/c3z3Z3IM39jfDZhHxN5+vQ8BjuennbFxSjnA0jv/+QxM0VZvwrWfZGq/tLrQPYNeBZvSF4vhknR+fvXMaJJNw0zaYU+nG7AoXLl0N4dDpDrx/tBV7P25BdakTn6zzo3aaO7n39b2PW3DmUi82LquE03Lz906VbhhQ4ol90QwMskWE2yZBtphhGdp7a8Q19PRkPm/r9zvQ2Zl+BkC7mce51j44bWZwHAdV0xGJ6XDaEmtBJI6hrze9duAMhkhYgRJVRu3VmXhgwQwfDhxvw6KaoqxXtw+E42hq6Ufj5T5caB+Ephswm3jUlLrwmUWlmDHNBZvlL9ezG68h1553C2uKMGuaE/sb2nHo1FUcOd2BpXOCWDav+KajCQDw1oeXcKFtAPevrIbVxOXkHFq7uAxnL/fhlbdP4evrZ6UUGxhjaGnrR/PlXvhcVvjdlgnrYQ+fezzPZdRhnLBAHQwGcfjw4eTPHR0dCAQC1z3f1fWXhVOdnZ0IBAIYHBzEb37zG2zduhVAojFFUURRURFCoRB0XYcgCMnfn2oWzvBh//F2HD/bhZULUu/hNbX0QTdYwa/2vlFVsQM+lwWHTneMGagZY9jf0I69R9swLWDDl1fPgC3F0Q5B4PGlz9TgtT2NeGP/BZhFftSMSPksFtfw7uEWfNzYhSKHhK+vn4XK4vTu3jmOQ2WxA5XFDgxG4jja1IWPz3Ri+3vn4LSZcedMH+xWE/7UcAV3zPThjpn+m7/pOIaHs+OaAYBB5Hm47WY4ZDOs0sTnpk6H3WqC12VB70AcDAxmk4DqUifs1tRrc9+I5zmU+Ww43zYAl3304LByQQlOnOvG+8fbsWl5VcrvbRgMnX1RXO4IoaUzjMsdoWT1O7fdjDtm+lBb7kJV0JHx6IRVErH2k+X4ZF0A733civ0N7TjS2IlVC0txx0zfmNMnx5q6cGgosOdyAaJDNmPdXeX43YFmHDrdibs+cfPYwHEcZIsIxhh6BxT0DMQQKLKiyGHJq/MPmMBAvWzZMjz33HPo6emB1WrF7t278Y//+I/J58vKyiBJEo4cOYI777wTb7zxBlauXAlZlvGrX/0KixYtwoIFC/DKK69gzZo1MJlMWLx4Md58801s2rQp+ftTjcchoarEgeNnu7FifknKF4KTF3tht5oKZq41VRzHYXFdAG//7yW0doVR5rt+2FfTDOz6oBknzvdgXnURNi2rSrt0n0nk8cDdtXj5nUa8vvc8vnJ3bU6LjUykxst9+J+DFxGKqlg6J4jPLCrNulfgkM1YuaAUn55XkizAMlwLuNRnw/pPVWT0vrpuIDpc+GKobGjAY4UsmWA2jb19Jx8UF9kQVXR4HBKKHJacZEJzyGa4bImEIKNVoHPbJdw5y49DpzuwbG4xvGOUq1VUHa1DAflyRwitnWEoQ1vAbBYR5UE7Fs/yo6bMBb/bktN29jgkfHFVNZbMDmDP4Ra89eElfHSyA3cvLsPMcvd1n9XWFcb/HLyI6SUOfC6NBYiGwTAQjsNmNY07vbBghhd/vtCDPxxpQe00FzyO1EYhOI6DTRZhGAxXe6Lo7IuhpEiG2y7lTZU8jk3g8rddu3bh5z//OVRVxebNm7Ft2zZs27YN3/nOdzBv3jycPn0aP/zhDxEOhzF79mz88z//M8xmMw4fPowf//jHiMViqKqqwpNPPgmHw4HW1lY8+uij6O7uRklJCZ555hm4XK6UjydXQ98qx6HpQhdkS+6HvgHgxLlu/Hb/BWxZNxNVJTcPGHFVx1O/PoZFtX58fklmF9FbKZ2hbyAxrP/sjuP4RKUH9SumJx8PR1Vs/+NZtHSEsXpRKT6dxo3NaKKKhv986wx6Qwq+tnZmXt70DLddJKZh96HLaDjXDb/bgnuXV6Esy7m+8XT3x3DqYi8WzPDCkcFiGN0wEI5oKPHaYLOaksPZt1KmQ98TSYkn5ottVtOo7RGKqnjuNycws9yNL66qBmOJ2tQtHSFc7gzh8tUQOvqiyUWyQY8V0wJ2lA/957abcxKYU/nOJirr9ePdw4kSuBVBO9YsnoYyvx3hqIpf/v4UOADbNn0i5WunYTAMRuIIuGV09EXhtI0/itEfUvDCzj+jzGfD19bOzOhv13UD4ZgOs4lHSZGcnPLIRrZD3xMaqPNNoQRqVTPwzPbjmFXhxn3XBKaxnGzuwet7z2PL+lmoSnPI81aLKRrsDiu0eHopGt/88CKONnbhkb+aD9liQkdvFP/1bhPCMQ33rajC7KrcDKOFIipeeus0IoqGLetmoTgH+9Hjqo5TF3txtKkL/aE4HDYTXLIZTtv1/7nkRK3r8S4KHrcNHxxrwZsfXkRU0bF8fjFWzC/Ju6G6ayUutirKg3YUOW5NRa3R5GOgBoArPWF09sXGzNz13sct+FPDFcwsd6OtK4xQNPHdMYs8yvy2ZFAu80/ctrV0bq51w8DRxi7sO9aGcEzDnCoPBqMq2rrC+Jt76lDiTW3V/nCQLg84UOS04Ep3eChYj3+jeORMJ/7n4EVsXFaZ1RSNqhlDRV9ElPrso456pCpv56hJ5kwij7nVRTh+thvrP1V+0y/fyeZe2CwiKvKwB3ituKpDZ4AgcIhpRlpD1ItnBXD4dCeONnUh4LHiN/vOwywK2Pr5WSj1ZbZdZzR22YSvrZuJl946g1f2NGLr5+vgG2PIcTyMMbR1hXG0qQv/d6EHcdVAkVNCRdCOgYiarD+s33DjyPOJIeFkAB8aHnXazJAtInYeaMbxpi4UF8n46toqFBfld2IbxhhCERUlPnlSg3Q+87ms6BlQoI3xnVg2pxgNZ7vR0RvB9BJHosfstyerfeUbgeexuC6AeTWJkpQH/+8qNN3AfSumpxykdcNAKKIlgzQABDwywjENUUVL5kQYzR0zffjzhR7sPnQZNWWutBblXssk8jCJZkQVDW3dIcwouzUV4UZDPeoMTHSPGgBau8J48fensGFpJe6cNfZdoarpeOrXxzG/xosNSysn7HiypesGIjENNWVuOFwWHD15Je293v/59hlc6Y4grukIemQ88LkZE7ZfvLs/hpfeOg2B57D183VwpzjfFYlpOHG+G0ebutDRG4Uo8JhT5cHCmT5UBOwjdj5EFA0D4TgGwnH0h9XE/0fiyccGIup156wgcFi1oBRL5wYnNDNXrgyEE3ufS7wjt/vcavnaowaAvpCCS1cG4bSPfj4zxia1/dKdrrrWYCSOzr5Yyus+kkF6lBEYRdVxtqUfkokf90a/d1DBf+z8MyqDdnzl7tqs2k7TDDAAtdMyD9TUo56iSr0yAh4rjjZ1jRuoz7UOQNWMvF7tbRgMoaiKymJnYtuN3QKLWURc1dNKdPCpTwSw44/nMKvCjS+smD6hiSC8Lgu+tnYm/vPtM3h5dyO2fn7WmHOzjDGcbx/AscauZC+51Cdjw9JKzJnuGXNEhOM42Cwm2CymMXsajCXabjhoz6z0QkBhVP8KRVS4bGYU50GQzncumxk2qwkxRRt1T28ht59DNqe8rmF4LcNY0ySSSUB5wIYL7YPjzld7HBI+e0cZ3vkosY7jZls78x0F6jzFcRwWzvBh96HL6OiNjpms/uTFXlglMW/nphNpH1WUeG3JpAg8z6HUZ8O51v60gm1dpQffqp+NgNt6Sy5cwSIZD65JrAZ/ZXcjvr6+7rp5qv5wHMeaunCsqQv94TgsZgF3zvJjUa0PwRwNSXMcl7zQlQHwuK0527s8kSLRxPDktIANfAEHmVuF4xLfiaaWPkjmye09T5bhnnRF0A7PONMkTpuEgEdD503mq+/6RAAnm3vxzkeXUV3qzGgRZL7I/7Gz29i8miLwPIejYxTq0DQDjZf7UFfhzsu5KiDRq/I6pBEpUW0WEQ45/aIHQc+t7Z1N89vxwOdmoGdAwWt7GhGJqTjZ3INXdzfip//dgH3H2lDklHD/ymp8968WYP2nKnIWpAtVVNEgCDwqix0FMTyfL6ySCJ/LinC08Osnp0vXUwvSw4IeGbIkjlu4hOM43Lu8Cppu4LU9TWhq6Su4HN/DqEedx2wWE2aVu3HiXDfuvrNsRHKCc20DiKtGWlVnEot7NMhWYcIvoon9oSaU+uwjgivHcSguknG2tR+SKb97ENNLnPjS6hrseO8cnt5+HIwBTtmEFQtKsHCGL+X9mreDuKqDMaCq1JHXK9HzVcBjRe+gAl03skqVWkh03UAopqGy2JFyKlKe51AecKCppQ+abox5rnldFty/shrvfHQZ//XuWRQXyVixoAR1Fe68vubciAJ1nltY68Wpi71obOkfMQ996mIvLGYB00tSH/ZWVB0WSUA4mkhKP1HzvEpcB89zqAg6xuztyxYTXHYJ4aia1daHW2FmuRubP1ONUxd7Mbfai5pSZ96OYkwWVTMQVw3UlLlSSiVJRhIFHiU+GS1XQ2MuLLvVGGPQtEQVsVyf88kgHUw9SA+TzAKmBey4eGX8+eq6Sg9qp7lw4nwP/tTQjv/+4zn43RZ8en4J5lQVFcT3OL+vjgQ1pS44ZBOONXVdF6g13cCZS32oq3Sn1TOOqwYqgjaYRQHNVwag6xqsOQ6SmmZA1QzMmOa6aaGCoMeKxpAy6ataU1FX6UFdHi/am0yJrGMaqktd426dITfntkvo7o9BUfVbfsNjGAyqbiQCMwPAAIDBZmeIKhoMAzCb+JxUTNN1A+GohqpiR8a5zN12CWGXip5BZcx96EAiVfDCWh/m13hxsrkH+xva8dv3L2Df0TYsn1+C+dVFeT2CQd+oPMfzHBbM8OLAiSsYjMSTCyIutA9AUfW0V3tzSMwPm0QBtdMShRgGwnE45MzzFl/LMBjCioaaUldKyRcsZhFepwV9oThsVjodC5FhsGSvKNcV5W5HPJfIA362pR/mMUphZosxBk03oGoMus4SFwYwCDwPWRLhcUiwmEWYRR5mUUAw6MQVm4hwTEN3XwwD4UTtbYskZDTFoekGItHEcHe2BUeKvTIiioZYXLvpNYfnOcyt9mLO9CKcudSH/Q3t2HWgGe8fa8OyecVYNMOXdgriW4GujAVg4Qwf/tRwBcfPduPT80sAAKeaeyENFQdIVVzVIQ8FaQAwiQKmlzhxpSeCzr4o7LKY1bw1YwwDERXlAXtaF2y/O5HwYSKG1sj1cj1yMbyqv+yaVf0ke7LFhCKXhIGQCjkHN7BxVUdcNYY6yIlc67IkwGM3J68JZhM/btAVeB5OOZGER1F1DIbj6OyPIhLTYBJ5WMxCSueWNpRToarEAecYZT7TIfCJevJNLX0wCanN7XMch7pKD2ZVuHG2dQD7j7fhrQ8vYf/xdiydG8SdM/0Tuv0zXRSoC0CR04KKoB3HmrqwfF4xDMZw+lIfZpa70rqbjasGSm7I4jW8VcoqCWjpCEHKYt56MKIh6LHC60wvA5XZJMDvsaCzNwaHjXpkE4GxRBpPjuPAGGASuZQvrOMZDKsIeKzwuUffPkgyF/TI6Bvsg24Yad9AM8YQ1wwocR0cEoG/2GtJ9pJNWfbUJZMAyW1FkcuCqKKhZyCG/lAcBmOwSsKYhWGSQbo4N0F6mMUsosxvx6WrIbhukg/8WhzHoXaaCzPKnGi+Moj9De3Yc6gFBxqu4FNzgvhknR9iHuxcoEBdIBbV+rDzT8241BGCphmIxXXMTmO1NwAYYLCNMR/tcSS+xM3tA4joo1fzGc9wcotMtyb5XFZ09ysZXZTI+IZzbXtdFhQXyYjFdfQNKugLKYmgbeIgmdIP2oPhODwO6bbfjjZRTKKAYq+Mtq4InLabfycYY1BUHXGVgUMiHW7QbYXNar7pWpFM8dcl7TEwGImjqz8xNC7wHKySmBwl0zQDESX3QXqYxy4hElXRF4rDPs589Wg4jsP0EiemlzhxuSOE/Q3t+OPHrfjgxBUsrvPjU7ODOT/edFCgLhCfqPTgrf+9hGNNXRB4DmaRR3Vp6pXDNN2AWeDHXZxilUTMmOZGS+cg+kPxm1aqGRZVNEhmAWX+zJNbiAKPYJEVbV3hCUsLejsaTiJR4pXhH0oUY7fysFtNKPElcif3DsQwEFbBwFLeCRCOarDLJpT57ZTQZAIVOSzo7lfGzOLHGEMsrkPTEvPMTtmEUq8FsuXW1/QWBR4ehwUeR6KX3R9W0N0fg24wiAIPTTcmLEgDQ1s+vTaEY3pWC/HKA3Y8eHct2rvD2N/QjgMnrsAsCpg73ZvjI04dBeoCYTYJmDu9CCfO90AUeNSWu9O6S1biOnyum2f0Mok8KouduNoTQUdvFHarOO6cT1zVYTCgMpj9vlmPQ0Jnb3TM4gQkPcPVf8ZKInHtnKOmGwhF4ugeUDAQjoMDB4vEjzqEGVU0mE08ygNjb70jucHzHMr8Npxv+0sWP8NI9Jw1zQDHcXDZzXDZJcjSrQ/OY7FKIqySiIBbRjimondQgcchTXh2MFHgURG042xLH0SBy2p0rsRrw1+tnoHegRikCapKlioK1AVk4QwfPm7sgqoZmJ3mam/DYONuX7gWz3Eo8dpglURc7hiE2SSMeneq64k5sJoyd04WXgg8jxKvjEt5tIe0UCmqDlU1UF3qSmlhnyjwcDsscDssiKs6BqMquvtHru5V4jqQoxszkhq71QSXTUJvKAae4yHwgMsmwe2QYJUmPnFRNnieSyvXdy5YJRGlfjtaO3JzHXHIZkx2PjMK1AWkzG+Dz2VBfziOGdNSX+1tGAwCz42a7H88brsEySTg4tVBRKLadatPbyy0kStOuwSpN5Z2wY5UGEai5KIg8FN6K1gkpoEDUFOW2Z5ms0mA1yTA67QgFtcwEImjpz+GiKJD4DjUlDnzakXs7aDEK0My8XDYzLCaRRrJuIkiRyKR0mBYhU0u/O964f8FtxGO47BxWSVCUXXMVZWjUVQdboeU0VyiVRJRU+pCS2cIA2EVjqGTfjCqosSX+y05PMeh1CfjfPtAToOBrhsIRVUUe2WEIhpCETXtBSeFIBRRIZkFVBU70jpHxmIxi7CYRfhdVsTiOniOy0myC5Ies0lAcYq1nMk1RU5ifRNy03+rUaAuMBXB9KtkqZqR1dBTYt7agY6eCK72RgEAXpcEv2tituTYrSbYLSYocT0nQUHVdEQVHZXFTrjtErxOA5euDmIwrE6Z7WCMMQyEEyvvpwXsOR+W5jiOMo6RgiIKPCqDTpxt7YMo8AU9CpG/kxskJ4YTXGR7keWHVlRWFTvgsUso9Y4stJEridWbiW1E2YopGuIqw4wyV7L3Lw5VdnLIJgxG1Kw/Y7IZBsNgWIXPZUFFMc0dEzJMtogo9dkwGFWhG4VRx300dIs8xcVVA07ZlLOLt8suZZ3yLxU2iwluu4SIomZ8kxGJJsotzihzjOiZCzyPiqADlztC6A8rBbslbHhIv8Rngz+FVf2E3G68TgsMxtDVF4Om6zlL9nMrUaCe4hRVR9BTmFmjAh4rmloUWMzppb0czsJlt5pQMc7q5ESpPDu4TqB3UEl533i+uHFInxAyEsdxCLhl+FxWRBUNvQMKeoeS/Ujm8XNL5AsK1FMcB8BWoIUSrJKIImd6+Y4NI5F72utMDM/fbF6K5zlMC9jBc0DPYBwOWSyIYK3EdWg6Q02ZCzZLYf77EnIrXZtFrdgrIxRV0TMQS+YNsEpC3uZvoEA9hana9UU4CpHfLaN3sDelgh3JYWCvLZmFKxU8x6HUbwfHhdHdn8g3ns/BOhrTwPGJbVKpVCgjhFxPFHi47RLcdglxVcfAUOrTSESFwCPvtsDRt3wKU+IGir2FnYdZMgnwu6zo6o+Nu50q22Fgfmg7B8dx6OyPwpmjsp/ZYoyBMcBgDIbB0B9WYBIFVBbbC/oGjJB8YTYJ8LkSxYSiio6BsILugUTqU7PI50WdagrUU5gBNiXqA3tdVnSPUwYzFteg68h6GJjjOJR4ZXAc0NkbnZCeta4b0A02FHivD8Q8x43IgMQhMTwvCjxEQUDAI8PkocIlhOQax3GQLSJki4iAJ5H6dLgq2GRPH1KgnqJSKcJRKEwij6DHivbuMBw3rM4eb2V3JjiOQ3FRIlhf6Un0rLMdAtN1AzHFgM6MZDrWRODlYBJ4iCIPnuPA8xwEnkv+P88nevrX3iz4/Q50dg5m+2cSQsZxbepTVTOg6ZO7tYsC9RSlxHV4XZa8GL7NBY9TQkdfFJpuQBR4MJZIBypbTKjM8d7hRLBOVAJr747AkUGwNgyGaFyDrjOYRB5+jwUO2UxJQwgpMKah+t2Tia4aU5RhMDhvYSL8iSbwPIq9Mlo7wrDLHAYiKrwOCaW+m6/szlTAI4PnOLR2hVMK1sMlB1XNgMBz8DgscNsTwXmq3DARQm49CtRTUKZFOPKd2y6hozeK/lAcpb70VnZnyjf0GS2dg3DI5hHBmjE2VKkqUQ/YbTfD47BAlvJr1SghpHBNrSs5AZBIcuK0ZVaEI5/xHIdpfhs0nd3SBB9elwUcD1y+GoJdFiHwPOKqnkxx6pDNKPVaIFvypx4wIWTqoEA9BamaAdcUredst07O31XksIBDIliD52CTBJQH7LBbzZM+f0UImdooUE8xuSrCQUbyOCywDCVCmAqr6QkhhYGu5lNMXDVgt+auCAe5Ht0AEUJuNbqaTzGKqsMzRYe9CSHkdkSBeoop5CIchBBCRprQQL1r1y7cc889WLNmDV599dURz586dQpf/OIXsW7dOvzgBz+ApmnXPX/y5EnMnTs3+XNbWxsWLVqE+vp61NfX4xvf+MZEHn7BUTUdFqmwi3AQQgi53oQF6qtXr+LZZ5/Fa6+9hp07d2L79u04e/bsdb/z93//9/jRj36Ed955B4wx7NixI/lcNBrF448/DlVVk4+dOHECmzZtws6dO7Fz5068+OKLE3X4BUmJG/A4qC4xg7jbUQAADHZJREFUIYRMJRMWqD/44AMsWbIEbrcbsixj3bp1ePvtt5PPt7a2IhaLYeHChQCA+++//7rnn3jiCWzduvW69zxx4gQaGxtx//33Y8uWLThz5sxEHX5BmipFOAghhPzFhC1h7ejogN/vT/4cCATQ0NAw5vN+vx9Xr14FAPzhD39ALBbD+vXrr3tPSZJw33334YEHHsC+ffvw7W9/G2+++SbM5tQWT3m99mz+pKS2rhBcLltGQTESU6GoOiSTADmLSk830nQDktWMaaXugkhX6fc7JvsQCha1Xeao7TJHbZedbNpvwgI1YzcW7MN1AWSs5zs7O/HCCy/gpZdeGvH8ww8/nPz/VatW4emnn8b58+dRV1eX0jF1d4dgGCM/N20ch/7+MFQl9UAbUzQomgGnbIbXZcHljhAGBQ7mHO3HDUdVeF0WdHWFcvJ+E4kqQGWO2i5z1HaZo7bLznD78TyXUYdxwoa+g8Egurq6kj93dHQgEAiM+XxnZycCgQD27t2Lvr4+fPWrX0V9fT0AoL6+HqFQCC+//DJ6e3uTr2GMQRTze1+rEtcxEI7DbBJQW+bC9BInnLIZ04sdiMV16DkqnzbVinAQQghJmLBAvWzZMhw8eBA9PT2IRqPYvXs3Vq5cmXy+rKwMkiThyJEjAIA33ngDK1euxJe+9CW8++67yQVjALBz507Y7XYcOnQIr7/+OgDgo48+gmEYqK6unqg/IStxVUd/KA5R4FFd6kJ1qeu6oW7ZYkJF0IFQVBt1dCEdhsHAT8EiHIQQQiZw6DsYDOKRRx7Bli1boKoqNm/ejPnz52Pbtm34zne+g3nz5uGpp57CD3/4Q4TDYcyePRtbtmwZ9z1/8IMf4NFHH8XOnTshSRKefvpp8Hx+bQWPqzpiig6rRURNmQs2y9glDt12CXGvjis9iZrHmc4tK6oO1xQswkEIIQTgWLbduQKSqzlqlePQdKHruh6yqhmIKhokk4gSnwyHNbXAyxhDS0cIfeE4HHJmi8sGwnFML3HCUSBD3zTflTlqu8xR22WO2i472c5R01hpljTNQFTRYTLxqAg64LSZ0+rZchyHUr8NimogGtNgtaT3T0JFOAghZGrLr3HjAqLpDANhFapmoDxgw8xyN9z2zIafBZ5HZXHiLktR9bReS0U4CCFkaqNuWAZ4AALPocRrSwRnPvu5YZMooKrEibMt/RB4LuXAq6g6gh5r1p9PCCEkP1E3LAM+j4y6Sg+KnJacBOlhVklEZbEdkaiW1lw6FeEghJCpiwJ1BkSBhzBBq82dNgmlPhsGI/GbbttSNR1WKsJBCCFTGgXqPOR1WeBzWTEYUcf9PSrCQQghUx8F6jzEcYn5b4dsQiSqjfl7jIGKcBBCyBRHgTpP8TyH8oADgsAjFh8ZrHXdgChykHKUK5wQQkh+okCdx0SBR1WxA7qeSKhyrVhch8chFUSlLEIIIZmjQJ3nJLOAymIHojENuvGXYK0boCIchBByG6BAXQDsVhOmBewIRRIFPBJFOEBFOAgh5DZAV/oCUeS0QFF1dPbFYBI5uKkIByGE3BaoR11AgkUyXDYzwjENLjsNexNCyO2AetQFhOc4lPltEAUqwkEIIbcLutoXGFHgUeZPv0waIYSQwkRD34QQQkgeo0BNCCGE5DEK1IQQQkgeo0BNCCGE5DEK1IQQQkgeo0BNCCGE5DEK1IQQQkgeo0BNCCGE5DEK1IQQQkgeu60yk/F87opY5PK9bkfUfpmjtssctV3mqO2yw/Ncxm3IMcZYjo+HEEIIITlCQ9+EEEJIHqNATQghhOQxCtSEEEJIHqNATQghhOQxCtSEEEJIHqNATQghhOQxCtSEEEJIHqNATQghhOQxCtSEEEJIHqNAnaZdu3bhnnvuwZo1a/Dqq69O9uEUlC1btmDDhg2or69HfX09jh8/PtmHlPdCoRA2btyIlpYWAMAHH3yATZs2Ye3atXj22Wcn+ejy241t9/3vfx9r165Nnn979uyZ5CPMX//2b/+GDRs2YMOGDXjyyScB0LmXqtHaLutzj5GUXblyha1evZr19vaycDjMNm3axJqamib7sAqCYRhs+fLlTFXVyT6UgnHs2DG2ceNGNmfOHHb58mUWjUbZqlWr2KVLl5iqquyhhx5ie/funezDzEs3th1jjG3cuJFdvXp1ko8s/x04cIB9+ctfZoqisHg8zrZs2cJ27dpF514KRmu73bt3Z33uUY86DR988AGWLFkCt9sNWZaxbt06vP3225N9WAXh/Pnz4DgO27Ztw7333otXXnllsg8p7+3YsQOPPfYYAoEAAKChoQGVlZUoLy+HKIrYtGkTnX9juLHtIpEI2tra8KMf/QibNm3Cz372MxiGMclHmZ/8fj8effRRmM1mmEwm1NTUoLm5mc69FIzWdm1tbVmfexSo09DR0QG/35/8ORAI4OrVq5N4RIVjYGAAS5cuxfPPP4+XXnoJv/71r3HgwIHJPqy89uMf/xiLFy9O/kznX+pubLvu7m4sWbIEP/nJT7Bjxw4cPnwYr7/++iQeYf6qra3FwoULAQDNzc148803wXEcnXspGK3tVqxYkfW5R4E6DWyUQmMcR6XfUrFo0SI8+eSTkGUZRUVF2Lx5M/bt2zfZh1VQ6PzLXHl5OZ5//nl4vV5YrVb89V//NZ1/N9HU1ISHHnoI3/ve91BRUTHieTr3xnZt21VXV2d97lGgTkMwGERXV1fy546OjuTQGhnf4cOHcfDgweTPjDGI4m1VDj1rdP5l7syZM3jnnXeSP9P5N74jR45g69at+Lu/+zt84QtfoHMvDTe2XS7OPQrUaVi2bBkOHjyInp4eRKNR7N69GytXrpzswyoIg4ODePLJJ6EoCkKhEH77299izZo1k31YBWXBggW4cOECLl68CF3X8fvf/57OvxQxxvCTn/wE/f39UFUV27dvp/NvDO3t7fj2t7+Np556Chs2bABA516qRmu7XJx7dEuZhmAwiEceeQRbtmyBqqrYvHkz5s+fP9mHVRBWr16N48eP47777oNhGHjwwQexaNGiyT6sgiJJEp544gk8/PDDUBQFq1atwvr16yf7sApCXV0d/vZv/xZf+cpXoGka1q5di40bN072YeWlF198EYqi4Iknnkg+9sADD9C5l4Kx2i7bc49jo018EUIIISQv0NA3IYQQkscoUBNCCCF5jAI1IYQQkscoUBNCCCF5jAI1IYQQkscoUBNCCCF5jAI1IQXsoYceQk9PD7Zt24azZ89e9xgAfPazn8WJEycyfv9sXz+eXB4nIVMZJTwhpIANFzb55S9/OeKxfFcox0nIZKNATUiB+v73vw8A+PrXv46zZ89ix44deO2115KP/eIXv7ju99977z288MILUFUVFosF3/ve99LKDjfW65977jm0trais7MTra2tKCoqwrPPPotgMIiGhgb8wz/8A1RVRUVFBdra2vDoo4/ijTfeGHGc27dvx2OPPYaenh7U19fjkUceyUUzEVL4si+VTQiZLDNnzmTd3d1s9erVrKGh4brHGGPJxy9cuMA2btzIenp6GGOMNTY2suXLl7NwODzu+6fy+p/97Gfsc5/7HBscHGSMMfatb32L/fSnP2WqqrKVK1eyvXv3MsYYO3jwIJs1axb78MMPRz3Oxx9/nDHGWEdHB5s7dy5ra2vLZVMRUrCoR03IbeDAgQPo6OjA1q1bk49xHIdLly6hrq4uq9cDwF133QW73Q4AmD17Nvr7+9HY2AgAWLVqFQBgyZIlqK2tHfMzhvMf+/1++Hw+dHd3o6SkJK2/k5CpiAI1IbcBwzCwdOlS/Ou//mvysfb29pRLFY73+j179sBisSQf5zgOjDEIgjCihrYgCGN+xrWl/4bfgxBCq74JKWiCIEDTtJs+tmTJEhw4cADnzp0DAOzbtw/33nsvFEVJ6XMyeX1NTQ3MZjPef/99AEBDQwMaGxvBcdyYx0kIGYl61IQUsDVr1uDBBx9EOBwe8di///u/Jx+rra3F448/ju9+97vJwvUvvPACZFlO6XMyeb0oinjuuefw2GOP4ZlnnkFVVRV8Pl+y9z3acRJCRqIyl4SQCfMv//Iv+MY3vgGfz4f29nbU19fj3XffhdPpnOxDI6RgUI+akNvY7373O7z44oujPrdp0yZ885vfzOr9y8rKsHXrVoiiCMYY/umf/omCNCFpoh41IYQQksdoMRkhhBCSxyhQE0IIIXmMAjUhhBCSxyhQE0IIIXmMAjUhhBCSx/4/z489PFV0wFgAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 540x360 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"sns.lineplot(data=causal_df, x='title_length', y='click_through_rate')"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"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.6.6"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment