{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# DoWhy example on ihdp (Infant Health and Development Program) dataset" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "# importing required libraries\n", "import dowhy\n", "from dowhy import CausalModel\n", "import pandas as pd\n", "import numpy as np" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Loading Data" ] }, { "cell_type": "code", "execution_count": 2, "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", " \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", "
treatmenty_factualy_cfactualmu0mu1x1x2x3x4x5...x16x17x18x19x20x21x22x23x24x25
0True5.5999164.3187803.2682566.854457-0.528603-0.3434551.1285540.161703-0.316603...1111000000
1False6.8758567.8564956.6360597.562718-1.736945-1.8020020.3838282.244320-0.629189...1111000000
2False2.9962736.6339521.5705366.121617-0.807451-0.202946-0.360898-0.8796060.808706...1011000000
3False1.3662065.6972391.2447385.8891250.3900830.596582-1.850350-0.879606-0.004017...1011000000
4False1.9635386.2025821.6850486.191994-1.045229-0.6027100.0114650.1617030.683672...1111000000
\n", "

5 rows × 30 columns

\n", "
" ], "text/plain": [ " treatment y_factual y_cfactual mu0 mu1 x1 x2 \\\n", "0 True 5.599916 4.318780 3.268256 6.854457 -0.528603 -0.343455 \n", "1 False 6.875856 7.856495 6.636059 7.562718 -1.736945 -1.802002 \n", "2 False 2.996273 6.633952 1.570536 6.121617 -0.807451 -0.202946 \n", "3 False 1.366206 5.697239 1.244738 5.889125 0.390083 0.596582 \n", "4 False 1.963538 6.202582 1.685048 6.191994 -1.045229 -0.602710 \n", "\n", " x3 x4 x5 ... x16 x17 x18 x19 x20 x21 x22 x23 \\\n", "0 1.128554 0.161703 -0.316603 ... 1 1 1 1 0 0 0 0 \n", "1 0.383828 2.244320 -0.629189 ... 1 1 1 1 0 0 0 0 \n", "2 -0.360898 -0.879606 0.808706 ... 1 0 1 1 0 0 0 0 \n", "3 -1.850350 -0.879606 -0.004017 ... 1 0 1 1 0 0 0 0 \n", "4 0.011465 0.161703 0.683672 ... 1 1 1 1 0 0 0 0 \n", "\n", " x24 x25 \n", "0 0 0 \n", "1 0 0 \n", "2 0 0 \n", "3 0 0 \n", "4 0 0 \n", "\n", "[5 rows x 30 columns]" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "data= pd.read_csv(\"https://raw.githubusercontent.com/AMLab-Amsterdam/CEVAE/master/datasets/IHDP/csv/ihdp_npci_1.csv\", header = None)\n", "col = [\"treatment\", \"y_factual\", \"y_cfactual\", \"mu0\", \"mu1\" ,]\n", "for i in range(1,26):\n", " col.append(\"x\"+str(i))\n", "data.columns = col\n", "data = data.astype({\"treatment\":'bool'}, copy=False)\n", "data.head()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 1.Model" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAwAAAABPCAIAAADTBTX9AAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nOy9d1gb17Y3PBK992666WA6tsG9494T18ROnO6c+CY3557yJs5xie3Ejh077r2DARsbbAIYbAwGAUIgEB1RJCGEEEIVldHM98e+R49epBlGIue+93m+8/vDT6LZWmzt2WvttVcloSgK/Rv/xr/xb/wb/8a/8W/8/wnmOM/YbHZdXd3g4CAEQSiKkkgkd3f3lJSU0NDQKf5ViUTS0NDQ2dmpUqlUKpWlpaWNjU1sbGx8fLyVldVUKGs0mra2NhqNJhaLYRg2Nzc3MzMLCgpKSUlxd3ef4rQHBwfr6+vZbDaKogiCkMlkNze3pKSksLAwEok0FcoymYxGo7W3tyuVSu2CREdHJyQkWFtbT4UygiDt7e00Gm1sbEyj0ZiZmZHJ5ICAgNTUVE9Pz6lQhiBoaGiovr6exWIhCAKIu7i4JCYmRkREkMnkqVAeHx+n0WhtbW0KhQIsiLW1dWRkZGJioq2t7VQooyja2dnZ0NAwOjoKXiKZTJ42bVpKSoqPj89UKINpNzU1tba2jo+Pg2lbWFgEBQUlJye7ublNkTgMwwwGg06ni8VisNokEsnX1zc5Odnf33+KxFEU7e7ubmhoGBkZAZciEonk5uaWmJgYFhY2xbcJQRCHw6FSqRwOB0EQwJj29vYzZsyIjY21sLCYInGhUEilUplMplqt1u6WqKio+Ph4Ozu7KRKXy+WNjY26W9HCwiIkJCQ5OdnFxWWKxNVqNXinEokELAuZTPbz80tOTvbz85sicQRBwDsVCATad+ru7p6UlDR9+vQpiiwIglgsFpVK5XK5WvZ3dHScMWNGTEyMuTneyUIEAoGASqX29fVp36mNjU1UVFRCQoKNjc0UiUul0sbGxo6ODu07tbS0DA0NTUpKcnZ2niJxtVrd3NxMp9NlMpn2DJo2bVpycvLUJQyCIJ2dnTQabXR0FBzKEAR5eHgkJydP/VyGIKi/v7+hoYHL5aIoCt6pk5NTfHx8VFTU1N8pn89vaGjo6+uDYVj7TmNiYuLj46d4zEEQJJFIaDRaZ2en9gy1tLScPn16UlKSk5MT1rdI+hagjo6OoqIiBEGCg4OTkpICAgK0j4aGhhoaGrq6uhAEmT9/flJSkrFTzMnJEQqFrq6uSUlJkZGRlpaW4JFcLqfT6U1NTTKZzN/ff8OGDcYud3FxcUtLi6WlZWxsbEJCgnYfazQaJpNZV1fH4/Gsra03bNjg7e1tFOWenp7nz5+r1erAwMDk5OTAwECt7ODz+VQqtaOjA0GQ2bNnz5o1yyjKcrk8JydnZGTEyckpKSkpOjpaq/8pFIrm5mYajSaTyby8vDZv3qxdK4J49epVQ0ODubl5TExMYmKiq6sr+BxBkN7e3rq6Oi6Xa2lpuW7dumnTphlFeWBgID8/H4ZhPz+/1NTUwMBA7QEpEAgaGhra2tpgGE5NTZ07d65RlJVKZU5ODo/Hc3BwSExMjI2N1TKGUqlkMBgNDQ1isdjT03Pz5s3G8kxVVRWFQjE3N4+MjExKStIqxCiK9vf3A03O3Nx8zZo1QUFBRlGWSqU5OTmjo6OOjo4JCQkxMTFaGQ3DcE9PD5VKHR4eBsQDAwONIg7DcH5+fn9/v7W1NbgkaPkZRVEWi9XQ0NDf34+i6IIFCxISEowiDkFQeXk5jUYzNzcPDw9PTEz08vLSPuLz+TQaraOjA4bhuLi4JUuWGEu8ubn55cuXEAQFBAQkJSXp8o5YLKbT6S0tLePj435+fhs2bDBWE2Kz2fn5+Wq12t3dPTk5efr06VoKCoWitbUVXIScnJw2bdqEIwQNQiQS5ebmjo2NOTk5JSYmRkdHa/ebWq3u7u6mUqkjIyOWlpZr1641loPUavXjx4/ZbLaNjU1cXNyMGTMcHR3BI7AbGxoaBgYGIAhavHhxXFycUcQhCCopKWlubrawsIiIiEhISNC96vB4PHBOaDSahISEhQsXGku8sbGxvLycTCYHBgYmJiYGBARo36lIJGpqamppaVEqlQEBAevWrTNWjPf39z979gyGYU9PT3CoaymMj48zGIzGxkaJROLi4rJp0yYHBwejiAuFwtzcXLFY7OLikpCQoCtv1Wp1V1cXlUoVCARWVlbr1683VllRqVS5ublcLtfOzi42NnbGjBna6SEIAt4pi8UikUjLli2LiooyijiKor///ntra6uFhQW4B+re57lcLo1G6+7u1mg0KSkpxkpdCILq6+vfvHlDJpODg4MTExN1L1RCobCpqYnBYKhUquDg4DVr1piZmRlFnMlkFhYWajQab2/v5OTk4OBg7TuVy+UtLS1NTU1SqdTNzW3Tpk3G3lgEAkFubq5MJgN378jISO07ValU4K4rEAhsbW3Xr1+vK9wA/i8FSCAQXL9+ffr06WvXrsX/kSiKvnz5sq6ubvv27UQEOoqi9+7dGx0d3bVr16TXpoGBgYcPH86YMWPFihWTUoYgqKGhoaioaNWqVfHx8fgjFQrFgwcPZDLZvn37iJiaRCLR1atXAwICiChkFRUVlZWVW7dunT59OpFpZ2VlcbncnTt3TmqaGhwcvH//fnh4+Nq1a4lQbmlpefr06dKlS1NTU/FHKpXK7Ozs0dHRffv2EbGsSKXSq1evAv1jUoWspqamrKxs/fr10dHRRKb9+PFjJpO5Y8eOSTXU4eHhu3fvBgQEbN68mQjlzs7OnJyc+fPnZ2Rk4I9Uq9W5ubmDg4P79u0jIl4RBLl9+7ZEItmxY4dWxcSCSqXKy8sbHBzcs2cPQeNBYWFha2vr5s2bg4OD8UeiKFpaWkqlUrdu3RoSEkKEOJVKLS4uXrJkyaT7BIKgxsbGoqKiBQsWENTyBwYGHjx4EBcXt2LFikkNSCwWKysra/r06evXrydCXCKRXLt2zd3dfcuWLZMyskgkunfvnqWl5Z49e4jIbo1Gc/PmTYVCsWPHjkntAUBl5/P5e/fu1Sox+MjPz+/q6tq6davu3dIgEAT5/fff6XT6u+++S1BvplAoZWVlmZmZRFTh+vr6kpKSpUuXpqSkECHe29ublZWVlJS0dOnSSQ1IfX19jx49ioqKWr16NRHiY2Nj169f9/Hx2bRp06SyRSgU3rt3z97efvfu3UTMkzAMX79+HUGQ7du3T/qaFApFdnb22NjY3r177e3tiUw+Nze3r6/v3XffndRuhyBIYWFhW1vbjh07CBr5KisrKysrV61aRUQVrqmpKS8vJ7gBIAjq6urKzc1NS0tbuHDhpO+0p6cnNzc3Pj5++fLlRIgLBIKbN28Ci8ak15uRkZH79++7uLjs3LmTiHlSpVJdu3bNzMxs+/btk76m8fHxhw8fSqXSDz/88P+yIKL/RGNj47Fjx5RKJUoYGo3mt99+Ky8vxx8mlUp/+OGHrq4u4pRRFK2srDx9+rRGo8Efdv/+/ZycHKMoj46Ofv/99xwOB39Ya2vr4cOH5XI5ccoIgly5cqWoqAh/mEKhOHz4MIPBIE4ZRdG6urqffvoJhmH8YXl5effu3TOKskgkOnjwYF9fH/6wnp6eH374QSKRGEX89u3b+fn5+GNUKtXx48cbGhqMokyn03/88cdJd2xBQcGNGzeMoiyTyQ4dOtTZ2Yk/bHR09ODBgwMDA0YRVygUx48fb2xsxB+mVqt//vnn2tpao4gjCHLt2rXnz59POvLWrVtPnjwxijiKooWFhVevXp10WElJyYULFyZl3gkAIkihUOAPYzAYR44ckclkRhEfHBw8ePDg8PAw/jA+n3/w4MFJhcMEyOXyo0ePtrS04A9TKpXHjx+n0WhGEddoNBcvXiwpKZl05NWrVwsKCowijqJofn7+rVu3Jh32/Pnzq1evIghiFHEgtVQqFf6wpqamY8eOjY+PG0WcxWIdPHgQOPhwwOVyDx48yOPxjCIukUgOHz7c0dGBP2x8fPzo0aPNzc1GEddoNOfOnZv06ERR9OLFi7///rtRxFEUzcnJuX///qTDnj59evPmTWPf6du3b0+dOjXpSVRfX//TTz8ZpVGgKNrb2/vDDz+MjY3hD2Oz2QcPHgT+euIQiUT/+Mc/enp6tJ/8twLU399/5swZo2hpcffuXSqVivUUQZAffvhBLBabQLm3t/fcuXM4A548efLmzRsTKMMw/MMPP0ilUqwBXC73p59+MoEyiqK5ublVVVU4A44cOQI8uMaCw+GcPHkSZ0BRUVFxcbEJlBEEOXz4MM62EwgER48eNZZVAJ4/f15aWooz4MSJE0NDQyZQ5vP5P/74I86AV69ePX361ATKKIoeO3YM57xUKBTff/+9seytxW+//dbd3Y0z4Pjx46atCYqiBQUFL1++xBmAz7P4aGlpuXbtGs6AN2/e5OXlmUZcIBAcOXIEZ5sNDAz88ssvphFXq9UHDx7EudLI5fKDBw9Oelpj4fTp08AXiYUjR44YK7K1yMvLw5d1169fN/YY1qKhoeHOnTs4A8rKyp49e2YacR6Pd+zYMZwBPT09+HIeB0ql8vvvv8dRmiUSyQ8//DDpaY2Fn3/+eXBwEOspgiD/+Mc/Jj2tsfDw4UMKhYIz4OLFi5NqYFioqanJysrCGVBUVGSCagXA4XB+/vlnnAHt7e2XLl0yjfj4+Pj333+Pw4ZjY2OHDh0y9n4FgCCIrmD/bwXo6NGjppED+PHHH7F22N27d5lMpsmUi4uLsbYIn8+/ePGiyZQlEsnp06exnuL8IiI4ceIE1tGYl5dnrO1HF2/fvsU628RisclaLIqi4+PjOHsa5xcRwS+//IKlbhYVFdXV1ZlMuampCUvFwf9FkwLYpbCenjt3zjQtFgBBkEOHDmE9zc/Pb2pqMpk4iqK//PIL1q2js7PzwYMHUyH++PFjOp1u8JFcLjf55gDQ0dHx8OFDrKeHDx+eiqQSi8U4XH/69GmRSGQycY1Gc/jwYayn2dnZbW1tJhNHUfSnn37C0t6am5sfP348FeIPHz7EOmslEonJSicAnU7Hmd6hQ4dMu1kBjI6Onj17FuvpyZMncS66kwKG4SNHjmA9vXv3rq45wQTgeF3q6uqImHJxcOvWLSy7vlAoxFk0IsCZHr5wI4Lh4eELFy5gPT1x4oSxxkJdqFQq7bWZDEEQjUZLS0vD96QqFIpdu3bV19cbfLply5Znz54ZfDQ0NIQfvtDU1OTl5RUdHT02Nqb/dOnSpW/evDH4xezs7F27duFQ5nK5kZGRAQEBAoFA/6m9vb2dnZ1IJNJ/1NnZGRMTgx8uIJfLZTLZhx9+aPDpjh07cnNzDT5iMpn4MTEFBQUkEolEIsEwrP909uzZNBrN4BezsrJ2796NQxmCIDabnZCQMHPmTIVCMeGRtbW1l5fX0NCQwW/5+/vjO+YPHjwIpl1aWqr/dNeuXdnZ2Qa/2NLSgh+CUFxc7OTk5Ovr29vbq/90xowZnZ2dBr/46NEj/B0CQVB5ebmHh0d8fLxQKJzwyMLCYvr06Qb/qEwmAylvBmnqMktLS8u0adOSkpIm7DQSiZSZmfnq1SuDFLq6umbMmIE1Zy19HK7cs2dPVlaWwa8XFhZu3bp1UuIQNm+uW7eupKTE4NcfPXr03nvvESGOxZ7h4eFsNtvg1+vq6ubNm4clqXSJY7Gng4ODvb29QVEjEons7e2xAkR0iWNxKJlMnj9/PoVCMUhhYGAgMjLS4CNd+jgc+v7772MxUXFx8bp16yYlDmHz6ZYtWwoLCw1+PSsr6/333ydCHItV4+Lienp6DH69oqJi+fLlWDEfusSxWNXFxcXc3Fwmk+l/nc/nu7u7Y8XV6m5vLD41MzNLS0tramoySIHH4+HE22npDw0NYfHprl27sE6KioqKzMzMSYmPjY1h8en27dufPHli8OvZ2dk471SXIBafpqSktLa2Gvx6cXExTjCfLnEsPvXw8IBhWKlU6n+dw+H4+/tjZcDoEsfiUwsLi9jY2I6ODgiCyBAE1dXVzZkzR5+Wu7v74cOHSSQS8G3PmzcP6yeFhoYalFkikUg/7noC8R9//LG+vv7zzz/HijnFCnJEEMRg3K6W8pEjR/z9/UNCQrCEWmpqKp1O1/+cQqEYDKTXXZBz587Z29tfu3bt73//u/5IX19fg1qXSqXCmowu8dTU1BkzZhgU0xD2gigUCoMxm7qUr169KhAIamtrDXJdenp6Q0OD/uf19fWT7pAdO3agKFpdXW0wV8jNzc2ghEJRFCupVUtcIpFYWVm5ubl5eHgYHGlra2tQWRSJRAZT/XWn/ejRo/z8/PXr17948UJ/5Ny5c2tra/U/ZzAYycnJWGSbm5u1zHLz5s2HDx/u3bu3oKBgwvjk5GSDEkSpVOq/R4PMiMOVTk5O+icogIWFxQQdAovTnz17ZpA3SSQSljYskUj0X5PBlcnPz8diTy8vL319FIKg+vr69PR0IsRx2HPWrFmNjY36xBsbG2fOnEmEOARBWByanp5u8H4CUhexiOsuOw6Huru7SyQSfeIQBFlZWU3QIbDeKRafkslkrDBVfcGCRTwwMBCLVZ2dneVyuT7x5uZm/fuPwWXHYdWUlJSWlhZ94lQqVT/vweDRg8On8+bNMygEhoeH9ctP6K6Mln3s7e2x+BTrpIAMyXmDxB0cHLD4FGTgGySuVqsnBA5jEcfhUzs7O4OCt7u7OzY2FmvyusuOw6fx8fFAR5mA2tpa/cXEUiew+HTBggU1NTUQUIBUKpXBM7W2tpZEIuXk5MTGxtrY2JhQN2J0dBQrNUZL/MGDB62tratXr66urjaWPj5lHx+fq1evrlixQn9PA7i5uRncfFKp1KCaorsg3377LYqiH3/88eHDh4nPTSwWYylAWuKWlpYPHz7csGEDlunL0tJSpVIR/6O605bJZKdPn7548aJBYYS1IEKh0KC1Q5dyWFiYTCYzqOXgQKlUYulzWuIUCiU7O/vDDz/EsjI6ODhIpVLif1R32gKBICMj4/r16waVaWdnZ4NqqEAg0N/YWrKpqakEmQU1VIZ0bGxMf5MYZEZ8riTOsFic/re//Q2LN42SBgZXRiAQYLGnu7v76OioPh0EQfQFukHiOOzp7u4+MjKiT9ygsDJIvLOzE4tDyWSyRqMxSFyfgwwuOz6HTv2d4vDp1Ilfu3YNi1VdXV2xbnT6fxdrw2CxqmnvlODRgyVv8YVAbGysln0aGxv/cD7V5c0/hE+xiOPwqaOjI5ZGjkNfd9n/cD7VJY7Dp3Z2dv99ZKAo+vDhQ1D4aAICAwN/+eUXCIIeP37c1dUFQZCVlZVBpxoMwwaj2NRq9W+//WbwK1riWm0OK8QPK64F63Mt5UOHDgUEBEyfPh3LD1paWmowcaOgoMBgzprugqAo+sEHH0AQdOXKFf2RCIL8+uuvpk37559/DgkJiYuLw4rANXlBHj9+3NjY6OvrGxoa2tvbqz+SQqFUV1frf/769WuDKVoTFuTGjRtqtdrgHKYy7U8++cTV1TUwMBArLWsqC8Jisfz9/ZcsWWIwjrK1tdVgUHl3d3dhYSEW2Z9++knLLE1NTb6+vjNmzBAKhRPGy+Vyg5sHRVH9zWOQGfG5EmsHEiSOouipU6eweJM4cRRjZZqbm7HY8/z58wY30p07d/TjrgwSR7HZs7y83GAAU3Nzc1lZGRHi3d3dWBw6NjZmMKMKhmF9SWhw2fE5dOrvFMXm06kTLyoqwmLVs2fPGgz0uXbtmn5Cn8Flx2HV58+fG5TYtbW1+gLN4NGDw6dDQ0MGo4kVCoV+nK/uymjZ58aNG1h8iiAIViwO/rLr8uYfwqdYxHH4FIv4xYsX9QObsE58LD7Nzc0FNYcnoKKiQj+F1iBxHD7t7e0FkaMQiqJisXgq0cQoimZnZ2OltBw/fnwqAW4ikQgrmPz58+fGJpROANbcFAqFyVkJAIWFhVgZGURS2XEwPj6OpVNWVFTgZ59NihMnThicGwzDUwmvRlH01atXb9++Nfjo5MmTU4loU6vVWBGajY2NRJKHcfDLL79ghShOMdQXRdEbN25gpeaePHkSR5UkgpqaGqwbxZ07dwxKFuIYGRnBqixAo9HwE9AmBYIgWKmOIpHo8uXLUyGO4kqkqYTMA1y+fBkrNP7EiRNTCd9G/1l6zeCja9euTZrhjw8ul3v37l2Dj6qrq01LttUChmEsJsXZS8SBlWWGIMipU6emSPy3337DyieY4umG4p4UFy9e1NfGjMLAwEB2drbBR6WlpVNJPUFRVKVSYSlAXC7X2FIs+jhx4oTBz2EYxsljIIgzZ86AfAIyBEEODg5WVlbd3d365iYiGBkZYbPZWHW4169ff+fOHdMoQxB04cKFd9991+CjFStWgDqwplEuLy9PSkoyaCG0srJydXVlMBimURaLxS0tLfpOUICtW7deu3bNNMoQBF24cAErgnXu3LllZWVYkR+ToqamJioqyqDP2MzMLCAgACsEflLI5fKqqqrZs2cbfLpjx47Lly+bRhmCoMuXL2/ZssXgo/j4+Pr6eqO8Y7pobGz08/PDCnaZPXs2lmuVCPr6+pRKJVY3knfeeefSpUsmE5fJZCUlJQbDtiAI2rp1KyjoYhpxFEXPnz//zjvvGHyakJBAoVDEYrFpxCEIunr1KlYEpaOjo7m5OVbMOxGUlJSkpKRg+QWSkpKKi4tNJt7V1YUTGr958+arV6+aTFwsFlMoFKx0gW3btoFLrGnEEQS5dOkSFh/NmjXr1atXxnq3dXHp0iWsmqVubm4qlcpgqgFBFBYWYhU4JZFI0dHRWKkGRNDa2uro6IgVn7p+/fqbN2+aTFwoFDIYDKyTYvv27efPnzeZuEajuXbt2oYNGww+Xbx48fPnz8fHx02mf/78eawN4+3tLRQKsVIZiODx48fz5883+MjMzCwkJOTt27cmE29sbPT29v7v2FOgECEI8uOPP5qQr87n8/FT9lEUff78uQl1QYBtEL+yBZ/P/+GHHyatnKaP6urq69ev44/5+eef29vbjaUsFAr/8pe/4JdPLCsrI1KlSh+//fYbfmW8sbGx77//3tgacSiKUiiUSdXqs2fPmpCbLZFIvvvuO/zyiW/fvjXtFnj79m38u6lUKv3uu+9MyG1uaWmZdEEePnxoWtWlnp6eSYsqVVVV4ZfbwcLY2NikC85isY4cOWJCXQO1Wn38+HF8QSGXy7/77rtJy9MZxN27dyetDnfixAkTGBNF0fLycvxqNyiK3rlzB79mFRY6OjpwiiYAvHr16vbt2yYQFwgE3333Hb5UYTKZx44dM8FwqFKpjhw5gl/PE/CRaQVvrl+/XllZiTMAVNPBL4uFheLiYpyiCdoJvH792gTiDAYDv+4aiqIlJSWmCfPh4WH8CkYoinZ0dBw7dswEj4FCoTh06JDByBYtRCLRX//6V9NK9F26dAn/JNJoNAcPHjToxp0UhYWFubm5+GMuX76M5VLAR2Njo67hyuzgwYMQBJFIpIyMjKysrP7+/qioKIKRU6WlpZcuXQoODg4PD8dpUBcWFjY2Nvbw4UPd1k74YLFYx48fF4vFMTExOI0RbG1t4+PjT5065ezsTLDDl0qlunjx4rVr19LS0sLDw3Fy3dPT058+fdre3h4TE0NwQSoqKm7evNnb29vR0REUFIRVRz84OFitVt+6dSsqKopga08ul3vmzJnm5ubm5mZ7e3usZlXW1tYpKSm//PKLjY0NwTrrMAzfunWrqqqqra2tt7c3KioK61WmpaWVlJQ0NjbGxsYS7I5ZU1Nz8OBBMpkcFBSElQ8IQZC/v7+FhcWVK1ciIiII1p7n8/lnzpwRCoVkMjkiIgKrUYmlpeWsWbPOnj0LWsASoYwgyL179yorK83NzZ2cnHBaxsbGxvb39z958iQmJoZgE18EQbKyskDIXl9fX3h4ONbM/f39bWxsLl68iJPGqI+ysrIrV64kJCQkJibibFpHR8fo6OjTp087OTkRb43HYDAuX778wQcf4Le+srCwyMjIuHTpkkQiIdgWBoIgLpd79uzZuXPn6qdiTUBGRkZBQUFLS0tMTAzBfSiVSi9duuTk5LRx40asMUqlEnSnMTc3Lysri4mJIdh9T6PRPHjwoLe394svvsAZNjg4WF5erlQqX758OX36dOJ9rIqLi4uKiv7zP/8TR3jCMFxRUcFms0GuOA6vTUBTU9O1a9c++eQT/NZX7e3t3d3d169fd3FxIdhoBYIgNpv966+/Llu2DKdrJJfL/fXXXxkMRm9v7/DwMHFhK5FILly44O3tvWbNGqwxCoXi9u3bL1++5HA4IDWJYL85GIbv3LkzNDT0ySef4Azr7Oy8evVqU1NTd3d3eHg48T5Wz58/Ly8v/+abb3D2WENDw82bN7u6uuh0uq+vL/Fm3lQq9datW59//jlW2iyY+eXLl5lMZmNjo5WVFfEGhf39/WfPnl2zZg1OX46RkZErV660tLQwGAyxWExcoxgbGzt//nxQUBBOFywYhnNzc8vLy5ubmzkcTmxsLMF+cyqV6tatWyKRCEQdAUxshtrR0ZGfn+/n57d27VosFlUoFIWFhd3d3UqlcuPGjf7+/o8ePRobG0tLS8vIyMBSKWQy2a1btzQaTWZmJo5YrK+vr6io8PDwGB0djY6ORhAEvP5Vq1bh9OUBLZMSExMXLlyINYGhoaFnz56JxeKdO3dKJJKrV6/K5fKwsLAlS5bgtKZjMpk5OTne3t5r167FmoBKpXrx4kVHR0d6evqcOXMUCkVOTs6rV69gGN63bx+WeXZ8fPz27dtKpXLp0qU4E2hqaiorK3N2dt65c6eFhUVdXd39+/e5XO6CBQs+/PBDrHev1VSWLl2KNWZkZCQ/P18oFG7bts3Pz08gENy9e7e6utrT0/PPf/4zlv40MDCQlZXl5ua2fv16rBQ/GIZB677k5ORFixa9ePHiyZMn06ZNCwsLW716NZaKo1Kp7ty5I5FI8BtAMhiM0tLSjo6OEyOt/7AAACAASURBVCdO2NvbDw4O5uTkQBC0fPnyiIgIrG+9fv26trY2MjJyxYoVWEJQKBTm5+dTqdSvv/46KCgIQZCioiIGgxEaGrp69WqcxO/bt28jCIK/kfh8/rNnzwQCwaZNm8ARMjQ0lJ+fL5FINm/ejKXRwjD88OHD4eHhWbNmzZ49G0uOyGSygoKCgYGBefPmzZw5E3BxSkrKokWLsOYD8OLFCwaDERUVtWzZMqxlgWEYpAtERETgHDb6aGhoKC0t9ff3X716Nc55X1tbW1VVJRKJ0tLSVq5cSYSyWq3+xz/+YWlp6ePjs3r1ahwdrqurC5S73b17N1Y/VJFI9PjxYz6f/+6774LEZpFIdPv2bdC3Mjw8HIs4j8d79uzZ2NgYTmMvFEXLy8upVKqvr+/mzZutrKzUavXJkycVCkVmZiaOtieVSp89e8ZmsxctWqRfc0GL3t7ewsJCGIZXr16t0WiysrIiIiJYLFZsbOySJUuweF+tVhcXF7e1tUVHR+OsuUgkevbsGZfLlclklpaW+/bt6+3tffXqVWBg4OrVq7HOexRFa2pqqqurPTw8tm3bhjWHjo6OM2fO8Pn85cuXb9u2zc7Orre3Nycnx83Nbc2aNTiHd3t7e0lJCZlM3r17N9a+4vP5p0+fbm5uTktL27Nnj5+fn1AovH37trm5+YoVK3Aapw8ODhYUFEgkEvzGXoWFhTdv3rS0tDxw4EBKSopKpbp7965IJJo3bx7OyxKLxc+ePeNwOCtWrMAq9IUgyNOnTx8+fGhra7tp06bMzEwSiZSfn9/d3R0fH79o0SKsA06lUv3+++8dHR3x8fFLly7FmkNZWdmdO3eUSuXy5cu3bNlia2v79u3bysrKkJCQVatWYd1+URQFzaR9fHzeeecdrDl0dXWdP3+ew+FkZGRs27bN09Ozu7v78ePHHh4ea9ascXNzw5oVg8F4+fKlhYXF7t27sfaVSCT67bffqqurY2Nj33///YiIiJGRkbt371pZWWVmZuJ0sGaz2YWFhTKZbNu2bRMUfQPd4CEIGh4efvbsGbC4mpmZWVtbk8lkuVwOjPaWlpaZmZlAZ8zOzpZKpXv27CGRSFQqtbKyEkGQsLCwjIwMg75wlUpVVFTU398PQRBYRJCQpq1WlJCQMHfuXCDo3759W1RU9NVXX6nV6ufPn4+Njbm4uMydOxdr+9Lp9NevX0MQBMof2dvbKxQKGIbBtN3d3desWaPLMDU1NSUlJdOmTROJRGZmZikpKampqQbZVSAQ5OfnA9eSmZmZlZUVKL2FoiiCIObm5suXL9efVVtb2/nz5/v7+xMTE7/44guDLA3DcHFxcXd3N4lEMjc3RxDE3t4eLIhGo0EQJDY2duHChRNuuiKR6MaNG+Xl5S4uLgcOHMDqAtva2goSW7QLolKpVCoVWBAXF5e1a9dOOBUQBHn9+vXNmzcFAsGuXbs2b95scK8LhcKnT58C86mZmZmlpaWlpaW24iqJRFq6dKmuOiKXy7WhYHK5PDAwcMmSJQYNGxqNpqysrLW1lUwm6y8IiqIRERFLliwZGxs7d+7cRx99BDY0ULk6Ozvt7e2XL1+OdRp1dnaWlJQgCAIWxM7OTq1Wq1QqjUZDIpEcHBzWrl0LmkJ/9tln2m+BM0atVqenp6elpRnUQmAYLikp6erq0r5HBwcHuVwOXiIEQfb29hwOR7/chUqlevLkycDAQEJCwqJFi7BMGjU1NXV1dZDO3gavEuTMW1tbr1y5coK8rq2tLSkpyczMxLl/a5eltLSUw+F4e3uTyWRra2sSiaSNTCeTyQKB4G9/+5uxXaABhoaGCgsLgTBRKpVOTk7W1tZaDx0EQUlJSenp6SQSCfD49u3b8QlyOJwLFy588cUX3t7eEomkoKAApMuS/wm1Wm1mZgbkSXBw8LJly7CU1/b29sLCQhcXl3feeUdf7KrV6qdPn5aXlwMTI3inMpkMQRDwTt3d3XF0u76+vsLCQqVSuWDBAt1XUF1d3dDQ8Nlnn719+xaU3dK+U6VSqVarwW60tbVdtWoVlm7H5XKLioqEQmFISEhmZqaFhcWjR49gGN6+fTvYn21tbWVlZdqtDg427Ts1MzNbvHgx1oVBKBSWlJSw2WxQ1bCxsXHLli26t1YOh/PixQtATZ/9IQhKTU2dOXOmQU6RyWTnz5+vqKjw9fXdv3+/fgSMSCQqKCgAdRDIZLKZmZmtra1EIiGTycAZFBYWhnWvU6vVWVlZjx49srOze++995YsWTJh0yqVyqKiooGBAQiCtO8UzBy8Uy8vr1WrVhk8g2EYbmxsvH79OovFiouL+/vf/z7Bfo8gSGVlJaiaCBjZwcFBoVBo36mdnd3q1asNGpVRFG1qarp48SKbzZ41a9a+ffv0zXgtLS2vXr3SynNbW1sEQRQKBZi5ubn50qVLsYwLfX19p06d6u7ujomJ2bdvn75az2KxXrx4Adzi4J1aWFhooykgCJo1axZW1+TR0dErV65UVla6uLh8+umns2bNmvDqhUJhQUEBqIOg+061Uj0yMnLx4sUGJYxMJrt3715hYaGtre3WrVtXrlw5wdyuUChevHgBQo503ykEQaAmhbe396pVqww6WwwrQBMAlhjLWdPX13fr1q2PP/5Yy6tMJrOqqmp0dJREIrm7u8+dO1e/YBQAgiBisRjHtDM+Pn7p0qXQ0FBw9RSJRJWVlUBdsLa2Tk9Pj46OxjozJBKJtbU1vs0TQZCHDx/y+fxdu3aBk0+tVpPJ5ISEhJSUFCx1WKlUwjBM3OD5+++/A99/fHy8n59feno6lhkZRVGg5xGkDHyFbDY7PDw8ISFh4cKFOKZsqVQKRBURymq1+ty5c6WlpVFRUZGRkfPnzw8LC8MaDA7jSb1XXV1d9+7dW7lypZeX18uXLyUSibm5eUZGRnx8PI6NVCgUOjs76w8AvXjT0tJ0+5NLpdLi4mIWi0UikWbOnJmamorjJZHJZObm5vreq56enuzs7G+++UZ38yAIArQQBEESExPnzp2LrxCAysK6Y7q7u8vKyj766COD45ubm1++fEkmk5csWYJfKxzSeZWnT5/+5JNPcJwjL1++pFAoW7duxfdGIQhy5syZAwcOQP88JnX5fWhoKDs7+8svv8SfFT4UCsXdu3d3796tVCqxlIaKioq2traPP/4Yi8irV68aGxv3799vcPHBr9i/fz++VRx05O7o6IiMjFy5ciXODvnpp5/27dunFVBA3cT3u3E4nMLCQqlUGhwcnJmZOeHV5OfnazQag544UPAThz1RFKXRaG/fvoVh2NfXd9myZWBiDAbj0aNHH3zwAZaYhSAIhLviBCoA4tXV1Wq12tnZeenSpXK5PC8vb/bs2TglcAEmZX+NRtPQ0EChUOh0+ujo6NatW3Fqket/VyKR4JwRGo2GRqNRKBQGg8Hj8ebNm/f5558TdItAECQUCp2cnLDeqVQqrampaWtrg2G4paXFysrqo48+IthoHZrsDJJKpXV1dS0tLWKxmE6nOzk5ffPNNzgWxwmQy+XgumLwKYIg7e3tNTU1XC63qanJ2tr6wIEDiYmJBIkDXRznnXZ3d1dXV/P5/KamJpVK9e677+LUIp8AjUYjlUqxLLIQBHE4nKqqKhaL1dTUpFAoli1btnfvXoL+bgj7yNAHIQVoUqhUqvPnz8fFxS1evHjCo+Hh4aqqqoGBAaC0RkREpKSkYLlOsFBTU1NeXr5v3z5dP6hcLq+pqWEwGOCW4+fnl5qaSjDOYwJkMtnNmzetra137dplaWmp0Wiamprq6+vHx8dJJJKzs3NqampERATxF2AQIpEoNze3r6/Py8sLhmESiWRlZZWYmBgXF4cjmAiisrKyuLh4ZGQEhIy4uLgsXrzY2HXGolxeXg6sNQqFws/PLzMzE8eYOSkKCwtra2v379/v7u6u0WgaGxspFIpEIvHz81u5cqWxc378+LFcLt+xY8eEzxEEoVAoIFIvICAAy+CEBYFAcO7cua+++soglzIYDHC9Dg0NXbp0KcEAIAiC6urqenp6sBIbwbTLysqoVKqdnd3WrVtxIpAABgcHX79+vW3bNpwxwKje0dHx3nvvYVkUfv/9d29vbyxTIgRBQJJiJZUQQWlpqY+PT0xMDP6wpqamioqK/fv3T/gcRdHr168HBATgmPchCLp48eK+ffuwdFPAg8PDw2vWrJl0JjU1NWKxeNmyZfjDAEZHR0tLS8fGxuzt7SdYmrW4e/duYGCgwSrzOBgfH6+qqmIwGGq1OiEhYf78+doDFYbha9euWVtb47QfITJtLpcLQVBiYmJ6erq5ublQKLx3756Dg8OuXbumIvSYTGZ5efnIyAiCIFKp1M7Obv369ZMq98SJv3r1is/nA+Lm5uarVq1KS0ubIlmNRtPZ2dnY2Dg8PAxBkIODQ2xsbEtLi0Ag2LZtG37026RAUZTJZFIoFIFAAGzbTk5OfX19Dg4Ou3fvnvop0NfXV1dXB2KfyWSypaXlyMiIh4fHzp07CQab4oDFYlGpVGA8gyDIzc0NZOmvW7fOtGNXFzwer6GhQRsITyaTRSKRlZWV1iv9L8IfowABAFfi559/jmUaQRCks7Ozvr4emDfNzMwCAwPj4uKIRGCp1eobN24AHcWgZsfhcOrq6lgsFrDg2dnZxcXFEY8yhiBoZGQkOzvbzs5u586dugJ0bGystra2s7MTaFqWlpaRkZEzZswgbqfRhfaECw0NzczMbGlpaW5uBlVfzczMgCEHx/+NDzab/eTJE7VaHRcXJxKJ2Gy2TCYLCwtbvny5Uce/PiorK6uqqmJiYhITE8vLywUCAVBn09PTiYdzajE+Pn7lyhVra+sPPvhAu9SDg4PgCIEgKDIycs6cOQTfHWi1+O2332JJkP7+fq3BaebMmUlJSUTEulKp/Omnn9577z0cDmxvb3/58qVSqQwMDFy+fDmRCO7ff/9dqVSuXbsWfxgIF2Cz2X5+fhs3bsRZimPHjv3Xf/3XpH9XpVI9evSIx+N9+OGH+pvh5MmTX3/9NT6F27dvx8bGTupQw8KZM2e+/PJLIteylpaWmpoa3Q5BfD7/119//eCDD3A8/QDnz5//5JNP9N9ve3s76H61Y8cOIpyrVCrPnj37zTff4A8DrqixsTEQ5YB1qYVh+NSpU+vWrcMJU5uA4eHh0tJSPp9vbm6+aNEi/fCy5ubmvLy8Dz/8kGC6gxYymayysrK9vR2CIDc3t8WLF2vNxmq1Ojs7e2xsbO/evaadxwMDA2/evAE5gCEhIampqYWFhSBMauonWWdnZ1VVFag+HBMTEx4eXlhYKBaLt2/fbrJqolKpurq6mpqa+Hy+QqEgk8mRkZGpqane3t4gSm98fPz999/HMUHhYGxsjE6nt7e3A9cSmUwGa+Lm5gZOgZiYmFWrVpnQZQEQb2lpaW9vB6YgBEGCgoJSU1O9vLymTlwikbS3t7e0tACfLwRB/v7+ycnJAQEBVCq1trbWzs5u8+bNpulVMAz39PS0tLRwOBzwiZeXV3Jysr+/f0FBQXd3d2ho6Pr164mb8UzGH6kAQRAkFosvXLiQkZGBVYNEFyiK9vX10el0oLUAq3tERERUVBRWxHtPT8/9+/e3bNmC01ZQO5Ouri4GgyEUCm1sbMbHx11cXGJiYmJiYvAz0fr7+7OyssDWMTgAhuGOjo7W1lahUAgMy66urtHR0bGxscQtARAEdXV1vXjxwszMbNu2bcDsodFoACuCmmYkEsnb2zsyMjI8PNwoygiCvHr1ikajOTo6rl+/XiAQVFZWAr+Gr69venq6r68vcWq6oNPpJSUl2jg4JpNZXV0NIjB8fX0XLFhglOrGYrHy8vLS0tL0SwRpr3dmZmYJCQlz5szBf2ujo6M///zzJ598gn8XgWG4pqamoaEBQRA3N7e5c+fiH6goil65ciU1NXVS03F/f39JSYlEIrGxsVm4cCH+OZeXl+fl5YUVHT8BLBbr+fPnfD5/1qxZixcv1hdnVCpVKpVilc2YAKlUmpWVpVQqP/jgA+2mksvl9+/fx+rsq4ujR49+9NFHxBNSdPHrr78Sd6JRKJSenh4QD1RZWfn27dsDBw4QSeH57bffPv/8c+3/ajSaoqIikDJm1GFw6tSpPXv2YKlKfX19JSUlcrncw8Nj1apVOMZ8CIJAN/VPP/10Uu6QSqUVFRWgtrKfn9+yZcsM3lvUavXdu3dtbW2xqjEZ/EptbS2NRtNoNHZ2dnPmzNEXoYWFhVQq9aOPPiKeGAgwQQ4sXLgQ9DHIyspCEGT79u0mG4wVCkV9fT2dTgcmcxBd6uDg0NXVVVRUZG1tvWPHDqPOYJVK1d3d3dbWNjQ0BEIh7ezs4uPjJ2Qo0+n04uJiNze37du3E5e9QJdqbm4G1iMymezk5BQXFxcdHa31bCoUivz8/N7eXiI2SH3iuueOra1tVFRUUlKSdgWkUmlubi6Xy129ejVWeSGDAMclk8nU9uADimB8fLz2fqtSqZ49e9bd3Z2cnGyw4SMWEATp6+trbW3t7+/XaDRmZmbgVcbExGhPIg6HU1JSIpVKV61ahd89/Y/FH6wAAYD8gk8//ZRguIkWCoWivb29o6NjeHgYuMwgCAoODo6MjAwMDNSaCp48edLd3b1//36j1ILh4eHm5uauri4QSUcmkz08PKZPnx4aGqovvygUysuXL1euXEnE3TsyMkKn0zs6OrRtgGxtbadPnx4RETFpPqpIJHr06JFIJMrIyNCNZQHg8XhtbW3d3d0gDAvEX4eFhUVFRU3qHIEgSCwW5+fnj4yMgABwZ2fnoaGhuro6JpOpUCisra2Tk5PT0tKMfU2gjriNjc3u3bu1UoPL5VZWVg4ODkIQZGdnN3v2bILsXVZWVltbu3fvXoO/CEGQ5ubmmpoapVJJIpFiYmLS0tIMGlpQFL1x40ZISMiCBQuI/F2JRALiBiAIcnR0TEtLw5rw9evX/fz8li9fToSsQqGorKzs6ekZHx8H7kKDs7148eLs2bNxXE4TgKLo27dvKRSKubn54sWLJ0z1xIkT3377LUFSEARxOJxbt255e3u///77ZDL50aNHM2fOJGLHVqvVP/74owkB0RKJJDc3F6cBtT7Ky8tFIhEIiSMeXnDu3DmQjs7j8R4/fiyTyTZu3GisSKVQKMBNNuHzxsbGN2/ewDAMXJ9EbCQcDufGjRtff/011mBdL7CVldXcuXNx0oggCKqrqystLZ0QD2AQMpkM7HAEQaysrIAeb9D2+fbt21evXq1bt44gz4KAGAqFolAowGE2Z84c7UnJZrOzsrJcXFzeffddEywEPB6vtraWy+WCY36CjCouLm5qaoqJiQHpUfikxsbGenp6mEwmkEvgJhwSEhITExMcHKz/dW3i54wZM5YtW4ZDH0GQgYGBrq6u3t5eUAYPOAdCQ0Pj4+MNvhomk5mfn29pafnOO+/gvzsYhgcGBnp6egBxCIKAdSAmJiY6OtqgTtzV1VVQUGBlZfXOO+9Mqm4ODg4ymczOzs7x8XEEQchksp2dXUxMTFRUlEFhxeFw8vLyQOzapFJidHS0s7Ozvb0dZK6AlQkODo6KigoKCpqwpFp/iL+//5YtWwgWKfgD8S9RgCAIGh4evnTp0saNG41ScvUBlMeOjo6+vj7dLoO2trY0Gm3lypWZmZkmExcKhUwmU1fthSAIvK2QkJDg4OAXL17Q6fT169dPanCaAKVSyeFwGAxGT0+PNgUDUI6JiTEYp0ylUisqKiwsLLZt24azgxUKRWdnJ9ARNRoNmUwmkUienp5gwlgBNKOjoyUlJSAsVxugoFar6XR6Y2OjVCqFYTg8PDw9PZ34XQ0cMAqFYs+ePRM0SD6fX1VV1dfXB1pMz5gxIykpCee0UCqVN27csLW13blzJ75zislk1tTU8Pl8CIJcXV31TTggsd+ogxb6Z6gjjjJUUFAglUpxYncMoqOjo7y8XKFQODg4zJs3b0IU+cmTJzdv3ky8AgeARqOprq5ubW2FYXjRokVgZ+bk5CQmJuIk92JND2TLNzY2/sd//AfBb5kWEF1YWBgREUG8LBAEQWKx+NNPP12zZo1Ry37u3LnZs2dXVFTY2dnt2LGDeKaCFkql8uTJk3/961/B/0ql0pcvXzKZTBKJFBcXN3/+fOKWeQqFQqFQ9u/fr3+UgkwRYJeNiYmZNWvWpGSlUunly5ejoqJwhJ5AIHj79m1PTw+ZTLaxsQH3EJyDvK2trbCwMDk5eeHChfh/ncvlUigUkMBrbm4eFxeXmpo6gakHBgZyc3NBBhbx26lSqQSCCGg8wDg6waUFVBMqlbp8+XKDgT4CgYDJZPb09ABDFEgvcnd3j4iIIFJdTGs7WbVq1YQaHCiKcrnc7u7u7u5uMENwqAcGBoaHhwcHB0+aZ4NzxqtUKjabzWQyWSyWtquBQqHw8fEJCQmZNIRDN6BC32ckk8nAAcfhcHQrv/v6+oaEhISFhU3KHSDywdPT85133pkwE6VS2d/f39fX19/fr1QqIQgikUjAsh4REREeHo4fGsHj8XJycpRK5apVq4j7hf9w/KsUIIDs7Gw+n//RRx/94Zod0F0KCgooFMr8+fOBpRoU3AsJCQkICDDZfSgSibq7u3t6eng8HgRBKIq2trbyeLwNGzakpaUFBgaaHKo2NjbW1dXV3d0NugWBD83Nzf38/IKCggIDAx0dHUUi0ZMnTzgcTlpamkF/h0FoNbnBwUFtKwwQYhUSEhIdHa2dM6hyIZPJAgMDMzMzdX9LR0dHTU2NSCQCs4qNjU1MTJw0vkcsFmdnZ4+Pj+/YscOg+gU8hloXu62tbXp6usGbKIfDuXnz5vLly7Hq/U8Al8utqqoCXmRra+ukpKT4+HhLS8uOjo4HDx786U9/Mi1IC0sZqqure/PmzZ/+9CcTssF1XRuenp4LFy709vZGEOTo0aMG812JAGhCVCoVQZDMzMznz58TV2J0UVVVdfz48YMHDxIP7jEhIPrMmTN/+tOfiI9vaGh49uzZt99+e/Xq1Q0bNhCJ8FAqlaA6y5dffmmUiX4CTp8+vWvXLpFIBFJZURRdsWKFsVcgCILy8/OVSqU23Umj0dDp9Lq6OsChcXFx6enpxLWEFy9e0Gi0/fv3T2BJjUbT0tJSV1cnk8lIJJKrqytOkqkuWCxWTk5OYGDghg0bDMoZhUJBo9EaGxtBdWkfH59Zs2ZhGQBARVBnZ+cdO3ZMKnuBTGhsbNS214iOjtZ15ehCqVTm5+d3dnZu3bo1PDxco9GwWCxg19HKOhKJ5ObmFhISEhoaaqx/lslkPnv2jEQibd++XalU9vX19fX1aa/EwBEB1JGwsDCjDFpisfjJkyeDg4MrV64ExZlYLNbAwACQsWDNraysQkJCpk+f7u/vb5RgAfdPoVC4du3ayMhIDocDjgDdjiX29vbgYjxt2jSjgtlFIlF+fv7Q0NDSpUtjY2P5fD6Xy51gKVCpVF5eXiEhIZGRkUZdM6hU6ps3b2xtbU27n/yx+NcqQBAEDQ8P37p1a/78+VOPzzcIGIZv3rwpkUg+++wzFEV7e3uZTObAwADwc4Hta2ZmNm3atGnTpvn6+pp22Mjl8mvXrtHp9IiICGtra6BNk8lkFEU9PDz8/f0DAwN9fHxMK5Sir75YW1sPDAzw+XwQbYNf1RcLgJmZTGZ/fz+4W4AFsbOzA8zW0tICw3B0dLR+wTQgoahUKkgPIZFIUVFRGRkZWBYmmUz28OFDoVD47rvv4p9VIKK8o6MDbDx/f//ExERdKw4wQX/88ccES0IDqNXqzs5OCoUCjEMODg59fX3r1q0jGGeDBZFI9ObNm56eHgiCzMzM/Pz8aDTagQMHTFOtAEZGRsrLy5lMJgRBQUFBDAbjyy+/NC2qBkCpVBYXF9+6dWvGjBk7duww1g5UWVlpY2MzNjZGJFteC2MDoo0KALp//76VldWmTZsgCEIQ5NChQ3/9619xLlF9fX35+fkkEmnz5s15eXn4FZlxoFAoLl682NvbGxwcjOO+nBQgYS08PDw5ObmhoYFOp6vVatCXKiMjw9gbFI/Hu3Pnjm6NbBaL1dDQ0N/fDxI+4uLiUlJSiJ8lAoHgwYMHjo6O27dv12V8uVxOp9PpdDqQQiBHNTExEf/6Chzinp6eW7ZswVJ9NBpNe3v7BHkyZ84cHD6Sy+U0Gu3+/fujo6Px8fH29vbaajfTpk0LDQ0NDQ0l2FdgAhQKBZvNBikXTU1Nzs7OiYmJ5ubmZDLZx8cnKCgoKCjIZAYXCoWDg4O1tbUvX75UKBTJycnApq5VF0JDQ02Lp1ar1Vwul8ViVVRU0Ol0c3Pz+Ph4GxsbjUYDbtEhISEhISEmqxTDw8NcLvfVq1dVVVUQBCUnJwNtz8LCwt/fPzg4ODg4eCp5agqF4smTJ0wmMy0tbSr3kz8W/3IFCODFixf19fUHDhwwTaBMiqGhoatXr6alpRnMWVWr1UCHHRwcFAgEQAEH/mBgJgUmwUntRmq1Oi8vr6+vTzedeHh4mM1ms1gsDocDIvWAJRAwqpeXl5+fn5eXl7E2MARBuFwug8HIy8sbHBwEli1zc3Mg8nx8fHx8fPz9/b28vIy1dQH33ODgIJfLFQqFvb29nZ2doN1peHg44KIJS6FWqxkMBo1GA/kXFhYWUVFRM2bMmKAPqVSqvLy83t7ejRs3ErRqDg4Otra2stls4OoGYX0hISF379719/fHaoo5KQQCQV1dXV5enkgkAkHfSUlJISEhpiVEAAAhTqFQ8vLyUlJSvLy84uLiTIig0gWTySwrK8vJyVmwYMGMGTPmzJljcrKeQqE4e/ZscHAwi8Wys7NbsWIFwdzUX3/99fPPPzczMyOSLa8L4gHRo6OjL1680C9VoA+JRHLq1KlNmzbphnCC9Ezd0pQAWvu/h4fH9u3bEd9zAwAAIABJREFUwXGojQEijt7e3vLycolEAsMwi8U6ffq0UV+fABaL9fe//z0oKMjV1RU/8mZSoCiak5PD4/GWLFnS3NwMtAcIgvz9/ZOSkox1nkIQJJPJ7t+/r9Fo3nvvPRsbG6lU2tjY2NnZCYpGgrKrWJYYfQAXamho6IYNGyb8QIFAALKfQCFQc3PzmJiYpKSkCeYr7aHO4XBAQjs4j4RCYXt7u6+v7+7du8Gd09hfKhaL2Ww2UHTEYjGJRCKTydqYAXD9WLVq1aSl0vUxPj7O5XLBaTI8PAzUMhBUCsNwf38/DMMZGRlbtmwx4T6sUqlGRkaAuQUU0tP2K5VIJIODgyQSacmSJevWrTNhRykUisHBQSD5R0dHddtmg6JBEASlp6dv2LBhiinDE9Dd3V1YWKhSqaZeR+APx/+QAgRBkEgkunjxYmpqqgl7jiAqKysrKip27dpFPN8SQZChoSE2m83hcLhcLgzD5ubmarXa3Nxco9HY2Nh4eXl5enr6+Ph4enoC/VepVN69e1cqlW7fvh0nrQOGYS6Xy2azh4eHORwOiqIgAF7LhF5eXt7e3j4+Pr6+vpNq1lQq9dWrVzY2Nps3b3ZzcwNzBkwIaowCKQPDsLOzs5+fn7e3t5eXF3G7AoIgb9++ff36tUgkiouLk0gkWrKgqC5oGjVt2jQfHx9HR8e2tjZQGwNwvqOjI0h2sLKyQhDkxYsXDQ0NmZmZBJ1ZWoyOjtJotI6ODhiG2Wx2V1fXnj17Fi5caEKmPUBXV9f9+/d37do1ODjY1dWlVqvHx8fJZHJiYmJqaqpRQfS6uH79uqurq7u7O4h8BzmGc+fONS1/QaVSHT9+fOvWrQwGg8ViQRDk6uq6YMECY9OGr169unbtWk9Pz/Hx8dLS0tbWViJBtadPn/7qq690J4OTLa8LtVp99OjRv//975MKeqAyTqqQtbW15efn79+/X/8WW1hY6OHhobUi83g8oN3qR2wQVIBkMtmbN2+AJTI4OHjJkiV2dnZHjhw5cOCAsXG7wFUEysGJRCIajXb48OEpVrsBruqCgoLU1FRXV9eAgICYmBjiTbj0oVarHzx40NnZGRYWJhKJQNa0vb19QkJCXFycsdezrq6ux48fh4WFrV+/HhTxB4mxwIZtbW0NEmNByq3uuQt8KNqqjGq1GkhXX19fkObS0NBQUlISEhKyceNG/H0llUo5HM7w8PDg4ODIyIi2YDFQdIAY9PPzCwgI0L5QkHirVCrxT2KFQjE6OgqkqzYWWwsymayV29qLYldX1/Pnz2EYxo+4h2F4eHiYx+Nxudzh4WGtZgbOHdBlAczZ19dXG4759u3b6upqR0fHjRs34sdoCoVCHo83MjIyPDwM0twgCNKV4T4+PkCGa++uoPeUjY3Nhg0bjM3+w4dGowHlRsPDw1evXm2ae+Rfjf85BQigvLycQqF88MEHJpe6wYdGowE9Wfbu3Tt1a5NSqRQIBEKhEDCDQqHQMoNMJuvs7BQKhStXroyJiXFxcZk2bRp+NqwuEATh8Xg8Ho/D4fB4PFCiBrAB+Nfe3t7DwwP0NfTw8AAakkQiycvLAxnRWIUGwNWHx+MNDg6OjY1p9S1AFoIgT09Pb29vT09PPz8//SUSCoW5ubljY2OZmZm6gcByuZzFYvF4PBaLBepfgQh/CIJApLNKpZJKpaB4qJWVlYeHh0Kh4HA4qamp+uUxCUIqlf7yyy+dnZ2gBDtIKwWpEMQvhUql8ty5c4mJiVrNG5xSoLorIBsWFhYbG2tUdYCysjKQ6giuYiwWq7a2VlvWwsfHB0QlE7Q5jY+PHzt27OuvvwY6B4fDef36NRDrbm5uBKM6FArFhQsXQDVnAG2ckEKhyMjIyMjImDAfjUZz8eJF3bxxAIPZ8voAAdH79++nUCj6OYwQBIlEIuDpP3LkCL5zLTs7WywW46TiHzly5C9/+QuNRqNSqdbW1ps2bTJo7cdRgGAYbmpq0s23SkpK0i7Is2fPvLy8CHrqmUxmZWUll8u1tra2traOiYlJSUkBRaEOHDhg7O1cqVS2tbU1NTUB4zQEQb29vdOnT//ss8+mYrPk8XigTgyorJGWlpaSkmLsPp+Avr4+kBAUFBQE8jCAlcLd3d3BwWF8fBz8BO2hrtFoQD8NPz+/adOm4ajU2kpjoFqBRqMBugIwtIhEIkANrC2CIE5OTsC+7unp6enpiZ+xVV5e3tLS4ufnt379erFYPDw8zOfzgZYAxujKSUdHRy8vL2C/x6eszQwHMcgIgggEgpGREe2cgX9KGzIBQZCvr6+np6eXl5ePjw++u2psbOzp06csFmvmzJmLFy9WKBR8Pn9oaIjP5/P5fNBMBmgVMAyD3HLQvNnb29vd3d3LywtH5wDBplwuNyUlhXiwKUH8LwlwJoL/aQUI+uddRK1Wv//++/8irVAqlYJOdXv37v2XFlOSSqW3b9+m0+lJSUlmZmagICn0zxYkgKPs7e3d3Nzc3Nzc3d3d3d0JepfHx8eFQqGu7gXp3Jw6Ozu7u7tdXV1Xrlzp7e3t4uLi7+9PxG6JIMjw8DC4NgEWBVykZSEIglxdXT08PEDwUFxcHBFbLpgtmKrwnwC3qJGRkdHRUQcHh7S0tLi4uPj4+JCQEKO84IODgzdv3ly6dGlqaqpu7S9gTvPy8gJJB/gX97Kysqqqqq+//lp/mEaj6enpodPpg4OD4GZsZ2cXFRUVERGBf98CTan279+vH1gmFAoZDEZnZ6dEIlEoFHZ2dgkJCcnJyfipcMePH//8888n/NGRkZHq6ure3l4iAR937tyZN2+evnMEQRAQ0Do2NhYeHq7tCNvV1dXT04PVe3lCtrzBMaWlpUePHgVnrcF2JdOnT+/r6wsMDKRQKBMWqqGhISIiwszM7MKFC3PmzMHqNARBkFgsvnz5MoVC+fjjj/EDCCYoQCDJnEajIQii0WjS0tIMOqQGBwfz8/M//fRTgzRRFO3p6aHRaKDGrpmZWVRUVEpKiu6F5+nTp+Pj40QK86hUKmAv4fF4YA/b2NhERETEx8c7OTlVVVWVlJTgt/LWB4Ig/f397e3twLkDSrZ6eXlJJBIul7t161acPjb4EIlEDAajoaGBSqV2dHSYmZm5u7t7eHh4e3v7+fmBZoseHh6+vr4+Pj4uLi4+Pj5ETlOhUKjVQl6+fNnf3x8VFRUYGAg0BjAG2Mi9vb19fX0JmuUkEsnIyAhQEYAiQqPRFAoFSE0HIs7Z2RmoIOCGSfCAGB0d5fP5QL8ZGRnp7+9vaWkBXXHc3NxAKzFLS0tPT09A1tfXl7gjaXR0VCAQjI6Ogv9obm4GtfESEhLs7e0Bcd0rsaenp8mxPqBlJ5lMnkp9Jhzi9fX1NjY2WPeT/234f6AAAbDZ7JycnFmzZhm8OP5Rf+LRo0cBAQEgmvJfBxA80dbWZjD2RSqVCgSC4eFhwDzaSx5weAMTKARBDg4OunoSERUB+Hf4fH5CQgKKonK5HNhRtQOAsHbTgaurK75hDEXRER00NTXRaDQXFxdwVwYOLxKJ5Ojo6O7uDrjdzc1t0r1OpVKBRcHS0lJbTF2tVvv6+rq6urq6ulpYWDg4OLi4uABJCrpy6lIoLCxsamrSz4IZGRkBtaO0GaqWlpbBwcHh4eEBAQG6REQi0fnz51esWDFpVUO5XA5oartsQhA0bdq0yMjIkJAQXX8BDMPAsYvT1huCIKlUSqfTW1paQEFYc3Pz8PDw6OjoCXdxlUp17Nixzz77DMt3CVJ+amtr5XI5SPmZNWuWrmVFo9H8/PPPf/7zn3EmQ6PR3rx5o9FogoODlUrl7Nmz8Z1Tur3lKysrSSSSNrRcJpOlp6fT6XRra+snT54YLJUUFxfX2dl57NgxXdMUBEHgT/v4+MycOfOrr77COi0aGhrKysqsra23bdt2+/btTz75BN9ffPbs2f3792sLaZJIpKSkpHnz5uH32ALmJa2ir5sEAFpBBwUFYTmhUBS9du1aeHi4wZ5ZoKYLg8HQdaaEhISkpKRMKIchEolu3rwZExMzaYgoKFvX1dUlEAiAYx2CoKCgIFBqBWxXUNonMzNz0t2uvb2AgiNMJhPci0gkElAO3NzcQGuzvXv3+vr64l+HYBgGgk4gEAgEApD7qZVyurYKZ2fn1tbW0dHRXbt2TVotRSgUamkKBAIQjKjtYAqEEqQjRXk8Xn9/v7OzM/5JrFAotLc18C/o0KmlrP2xzs7O4PpKo9H4fH5kZOSkpYq1l1gttPdYbfwN6BRmb2/v4uKComh9fT1oLmFs8MCkGBoaevz4sVQqXbRoEb5P3AQA1wSbzZ45c+b/ngBnIvh/pgABFBYW1tXVffHFF1PJgsFHfX19UVHRmjVriBedMw0IguTl5TGZzDVr1ujXrZ8U2rsL4HOtnqQPKysrFx3Y2dm9efOGzWbHxMToGzM1Gg24W2j/BW1yIQjSvnrwFVtbW1dXVxcXF2dnZ/Cv9kzi8Xi5ubkWFhY7d+4Ex49YLObxeGC2o6OjODS1ipebm5tYLH7+/Lm2gqJKperp6WlvbweRTFKpVCqVjo+PW1lZ2djYODs7g+beut79hoaG4ODg9evXA4IGvTMwDPf29ra3t2s1LRRFra2tQQoGlUpls9lffPGFsX4KrVEHhmFwktnZ2YHg9MbGRqFQ+MEHHxAkBcNwd3c3OBSBTglBEKhv6e3tfejQoU8//RSnqa0WIyMjoGIy+N/AwMCUlJTm5mZ3d3ciArSvr+///J//k5ycbG1tvXr1avz4xOrq6pcvXxYVFYH4G23kTV1d3Z///OfKysotW7bcu3dP/4uzZs1ydnYuKiqa8PkXX3xx/vx5S0vL77777v9r77vj2jiztUe9I3rvvRrbYJtqx2AbG9txd9w3Ltmss9mUjXezu0n2bsneTW6yyWaTzSZO4ps4zXZijCsuFONKEwZsgREgBBISSKDeR+X74/wyny7NQhKpev7gB2I0MxrNvO/zPuec52CmO0Kh8KOPPvrTn/4ENe1CodDxrpZKpVVVVVP1veru7r59+/a5c+cWL148o3qrjz/+OCsrS6fTcblciIqSyeSMjAxHG9ypoNVq33jjjd27d0Pmh0KhAKsL6BUDdBwCrNM4l9rt9uPHj0ul0scee2zcOet0Otgh5spht9vZbHZycnJiYuKkwyY4Jebn5xcWFmKKBeSFYCEki8Uil8vh4bVarUwmE4ImmZmZSUlJmOvs4OBgRUWFn5/fzp07LRYLzOIYv4EGPhj5QL5JNwkKCoJVHDz7Ex9StVp98uTJ4eHhNWvWUCgURwriOB85DmVsNhsWh7Dimup7GRsbq6ioGBsby8rKio2NHcc8JsJoNMJwh2Hi0gvQ1tZWV1eHomhRURGNRgOLTqVSCVQM+aZgHrQ3eAW0fxip4Oek3NFsNldVVfX29oaFhW3YsMH9jmDjdn7hwoW+vr7w8PD169d7ducIgty9e/fKlSvQFsPjetK3gO+YACEIYjAYPvzwQ39//+3bt7vZbXQaVFZWdnd3b9++3f22bdPDbrefPXv2/v37eXl5D+yi7BpMJtO4hQX0zxoYGOjt7WWxWFjlKlxPNpvt4+MDP6FP6qTdj/V6PYxBSqUSfqrVascNNBpNS0sLFLYwGAyMJ8HA4evrO/Hx1uv1GPECoCiqVqtbW1spFEpubi6sy8GqwNfX19fXl8FgaDQaGBAddxUWFpaQkODv79/U1HTlypX8/HwCgeA4rmHrKjKZHBAQ4DiohYSEoCgK3iFCoVAmkzU2Ni5cuDAoKCgyMhKsMlzIGNPr9WBKKZPJxGLxvXv3Fi9eDK5OcXFxM8ootFqtfD6fy+UODQ2hKFpTU5OXlzd//vy0tLRxOtY0GBgYaGlpGRoaunz5cllZWVJSUk5OzvSRFNBLdDpdbW0tn8+32WwLFizAUoXee++9/fv3Y4rXiRMn9u3bp9Pp0tPTL126hBEmm8325z//+ejRo59++qnRaLRarXDnKBQKJpP58ssvP//88xERERQKBcIZDAbjxo0b69evJ5FIGRkZK1asABtrDoeza9cui8UC/e03bdo08Rq+/vrrv/71r+GCQF1eU1OTWq0mEonR0dGLFi2qqKj4xS9+Mf2FMhgMnZ2d9+7d02g0QqFQKpWWlpZCFr/z6cBKpfLq1avHjh0rLCx0nOPBaNgZ/gpob2+vqKjYvHkzi8Xi8/n9/f2YUITD4eh0elJSUlJS0rjo4aQh8qGhIWgqPm/ePBqNZjabYZkBlJ1KpWJhbgaDkZiYGBMTQyaTHWPWMJIgCKJQKNra2hgMRm5uLlxtKpXq5+cHEzlwkUmZDUClUsEvcrkc8/eDDcAUIC8vD2QPWHQB/P39J4bPIK8IA4xO4LyHAQq7hEIhDodbsGBBREQEtlv4ZVLmYTQaNRoNjDYYoQFW53j0u3fvqtXq0NDQjIwMf39/0K4AbDbb5foMBEG4XO7Vq1ehM+CMPEKdAYfDqa+vBz3J4zvXarWnTp2SSqWZmZnLly+fvYl7tvHdEyDA4ODgl19+mZWVVV5ePkuHsNlsZ8+e7e7u3rVrlztpgE6itbX12rVrwcHBW7du/RaaumEYGRk5fvw4iUTasWMHm8222+1KpRLGI8efUJ6KvQvTqCEwx2KxgDOBDkSn0+l0OmQvjY2NHT16NCQkpKSkBBuP4CfoxnBHAaex2WzYUSgUio+PD+yZxWLZ7fbTp08/9NBDWOeKiYqx49INBlOoVYacDAaDkZKSEhoampSUlJqamp6ejmndKpXKUTDHhnUALKNbWlosFktxcbHRaNTpdKASUSgUuBRhYWHR0dExMTHOp2JADnJaWlpUVJREIhkaGsKGaYvFwmQywW1hUg/+cbDb7W+88UZaWhrkdkA8AkEQKpWamJiYnJw8/fzK5XI7OzvnzJlz584daEuEIEhoaOjcuXMTExMdR6txGTMWi+X69evt7e12uz0xMfHgwYNRUVEffPABFKV/+umnoAAJBAK9Xr9jx46+vr7e3l6Iwvj7+4eEhGD0msViwW0/ODjIYrFUKpXZbIb6FxA7AwICSkpKIC03MzOztrb22Wef7e/vRxDklVdewaJ4HA4HLE/gz1u3bnE4HAKBANNqRkaGY5+43t7enp6ecUbJJpOpt7e3s7NTIpFAzg2dTocaJQKB8O677z6wiwjEU/h8PpQ9wtchEolwONyhQ4dmtPBVKpUikWhwcHBoaAhcpnx9fXNzc8E1IyIiwtfX12AwQDqIQqEAbRWkBSwAjSAIm80GOQQUkebm5tOnTzOZzLCwMKPRqNfr9Xo9NDqg0Wiw7JlIAkwmE5vNdlwqBAcHy+XykydP0mi01atXazQaRzYDsjT2jAPglODZ9/kGwA9gDIFFDo/Hq6ystNlsK1eutFgsSqVSq9XqdDqtVgutJMYpSdj+aTQanCQQDgDGvSQSSWVlpVQqzcnJiYmJ0Wg0Wq0WqBLwSMgKwLoIOP5JpVLhJFksFkZoYAyxWq11dXUcDickJGTTpk3usJyJgL5Xcrl8zpw5Hk9AHhwcvHjxolqtzsvLm1j04D4c65Gdacf0Pcf3hQABmpubL168uG7dujlz5szSIVAU/eKLLwYGBlzo+ecC+vv7L126hCDItm3bXPO/cg1arfbEiROjo6OusT3HgQ8GFBhSYagiEolQrx4ZGZmeng7TPIlEgjxTKpXKZrOBM4GiQ6fTgfeYzWaM1gDduXHjBpfLXbRokY+PD8izIORQqVT4E2pJYPUJoNFoUGfb0tJSUVFRUFAgk8kEAgF4ycMDjxXnY9IXk8mEX1gsFovFYjKZvr6+CoXio48+WrhwoZ+fHyx/Yb6x2+3AtyC6hyU/0el04DGxsbFJSUnJyckTM1fOnj0rEAieeOKJcfPN8PBwf3+/QCCAvCIAgUCAAhlwLhi3kDp8+HBGRoajl6PBYIDkD6x6BcphkpKSxmlOb7311p49exwz7iGEB8tZoHqZmZlNTU3j/KO/+uqrTZs24fH4119//be//a3dbo+Kitq3b9/DDz9cU1MDQtGcOXPi4+NTUlIWLVqUlJTkghMStFLifoPGxsaRkRFI/4qOji4oKDhw4IDBYPjtb39bUVGxevXqtWvXQjMKKpXa09Pz6quvTrrovHz5Mo1GgwijyWSCPkcUCiUxMRHCi+O2//vf//70009jCbYSiQTc/SF5CLuXQkJCwFQdQk6Q9BMXFzdVeSNYxYC1qUKh0Gg0sBJQKBTQi5RGowkEApVKVVJS4qijTIzI0Gg0jUajUqkUCgWc29DQkFwux0ItEIwOCAjIzs728/MDecbv/wKPx0+Vj+J4zqBigpe01WoF2gSPG+wnKCgIAmHY8+v4O/bYOj7Cer2ey+VKpdLw8PC5c+fSaDTHp9jPzw9yk9VqtUajMRgMWq1WrVbDaKNUKiEgDsMOiUSCZRs8WV1dXUKhEKorICgGg4yvry/87lq4h8vl1tbWoii6du1alzPHJ4XBYKiqqhIKhUFBQevWrfNsjrBCobhy5YpEIomIiHj44YfdMSebFFA1JhaLoarX47zqu8L3iwABzp8/39raeuDAAedl5JkCzE5GR0fHTRKzBIh5j42Nbd26dbZjcI4wGo1g0Lx79+7ZyLK6du3a7du3Dx48OI4HYAPixFES+caCEn5BEEShUNTX1xcVFcGIAGUa2PgI8xnUqeFwONCrNRoNGP/b7fb6+vqoqCjITiWTyUajEcryDQYDiqJWqxVWmVCXi8fjQTmAUi8SiQRTkVwuLykpgTpb+AjYSpHBYEB2No1Gg1eUSqVUKhUKhSKRaFyUEIo1CARCX19fWVlZamoqrFxhVvPx8RnHisDucmBgQCKRjIyMYJ+UQCDY7XY/P7/e3t7o6Ojpm0qq1ere3t7+/n6JRIJ5E6Ao2tbW9vTTT8fFxU16h4N39quvvgpmLZDnm5ycXFpaGhER8Zvf/IbL5ba1tdlsNqFQKBAI5s2bt3nz5pKSklmqax0eHq6rq6utrb127dqqVavCwsJOnDjR2tqKIMjixYv//e9/v/fee88++2xCQgLWVQP6KHG5XDCIMxgMN2/eXL9+fUpKyjRNlHQ6nVAoFIvFx48fDwgIcKxaCg0NhaY008h+Go3mn//8Z2lpKTitg9uewWDQ6/XAyex2O0zzoFuQSCT4xd/fn06n4/H4rq6uixcvLly4MDAw0DEfFnJrIA3OYrEAZyIQCLAAIJFIUO8NMUS1Wt3d3d3S0hIUFJSeng5GxnBvk8lkaF5BIBDgRTqdDmoukUik0+lQr4QgCKgvSqWyt7e3ubnZ19d3zpw58PjY7XZIh8IWNgiCUKlUJpPJYDAYDIafnx8UKMGh4YOAzqfRaPR6vVarbW5uHh0dTU1NjYqKgqcej8eTyWRYMmHMBnZLo9GAwdBoNAaDAYuoiSSmt7f3woULdru9vLzcgwRFKpVWVlZqNJp58+Y99NBDHozpgMUaCJmTRnXdAbjAd3Z2+vj4bNmyxeODvMViuXLlCpfL9fPz27Rp07e5hv928H0kQAiCoCh64sSJkZGRxx57zLPyoyN0Ot2nn35qs9l27949e0fBAKxraGjooYcemqXGIFMdF/xk9+/f71mLTwRBDAbD//7v/6alpT2wmeI0sNvtH374YUJCQklJicViAf0JVoEKhQJ+gVfodLper6dQKCiKQoo0iqKQIrps2TK9Xk+lUk0mE+YqRqPRYDFqtVqtVivEL4xGo9lsNpvNUqlUq9WaTKbR0dHBwcHQ0FBQ7GHo9/X1hT34+/szGAzMTwymGZAWsF4r2M4hlXt0dLStrY1KpYaFhYGeb7FYLBYLmUyGvAHgVUwmEwrv4TpQKBTQ5AE2m02r1VZXV4vF4vT0dCzLEliUv78/+DlBHGTcmgzyapVKJZPJxHJLYQ/BwcHQvCU0NPTdd9+FEJjNZisqKmKz2VevXoW6/YyMjOTk5Lt37+7bt2/Hjh2zV6YwDiaT6cSJE6+88orBYMjOzjabzZ2dnQwGg8vlPv/88wQCQSgUZmVlweyYnJycmpqKnRsW0VMoFOBiN86+C0EQOp0eGhqq0+nUavVjjz0G122cQDIyMgK9aMDl0mw2k8lkcLoaHByMjIxkMBjgdwXMBkEQMplMIBBIJBKQDxwORyKRIFxrt9vhNtbpdG1tbSQSKSwsDHPSAq4D1ITBYPj4+ICFNMZgIMOJQqHg8Xjg6FKplMvl+vr65uTkgEYCtzeVSoU/4bgUCgVSf8BmHcRUuG8tFguDwTAYDDKZrKOjIzQ0FKKcWq2WyWSO+8lgMGBdQSQSzWYzEBf4yGQymcVi0el0uHXJZDKTyRwdHa2urmaxWNu2bfOU35tcLj99+rRSqUxISFi1apWnmkuqVKqzZ8/KZLLAwMD169d7cBaw2Wy3b99uamrC4/ElJSXj7DrdhNlsrq6u7u7uJpPJa9asccET/IGAUJfZbC4rK3O+480PDt9TAgRQKBRHjx6Frs4eT1/HoFarP//8cxRFH3nkEdc6hc0UYJZAp9M3btz4rZklqFSqo0ePstnsnTt3etx+qa6ujsvlPvHEE+6snM6fP280Gl3zLOju7v76668PHTo0LqBgMBh0Op3ZbFar1ZB5AATLZDIBrwJGBT97enr4fP6SJUtYLBaKohaLZXR0VK/Xj4yMwAIXy3UAozaYcjBOQ6PRYCYODAzU6XR0Or27u5vL5WZnZ5NIJKPRCBEZhUKBoqjJZILcIxKJRKVSYSVNJBL1ej1wOwRBCAQCnU4nEongpJCfn08ikWCyodPpQBNBMIA1PZlMplAooGyRSKTAwMCrV68++uijkBMaEBAAk/3w8DD0ZQTRpbi4GI/Hv/DCC8C3aDRadHR0UlJSW1vbk08+eejQoe/KwrW6uvrpp58eHR2FSkYEQaAN3+HDh+Pj44OCgqC/DXR/g9n96tWrhYUs8kurAAAgAElEQVSFJpMJyASNRiOTycCecTgcMBjIQGptbY2KisJkQrjgIH5gUgRM50QiEXrm9Pb2KhSKzMxM4NbwRgKBAHoG5NzgcDjgCvBfBEFMJhOBQIBcYwRBsrKyQBCCYi5QIo1GI3AdpVIJ3yzEeuD2AHqEIAgOhxsaGgKmPmfOHBwOx2Kx4BwsFgtwL7AkBbpGoVBMJhOFQoEmxHDbkMlkuEM6Ojr4fH5ubu7ixYuZTCYUY7vzfdlstpqamra2tvj4+HXr1nkk8dFsNl+6dInP59Pp9A0bNniKhRuNxitXrvD5fAaDsWbNGs+qMlwut6GhQafTzZ8/37NZOOBr2tnZabFYli5d6kKt8QPR1dV15coVlUqVn5//Ywp1TYXvNQECjI2NHT9+HI/H79mzZ6Yu9c7DbDafOnWqr69v48aNLnR+dgFSqfTkyZOT2vnPHoaGho4ePTpv3rypvO/c2fOHH3741FNPuRNSvHbtmkAg2LNnjwvvVSgUr7/++jPPPOPOotNgMHzwwQdUKnX//v0TJ34URcHtGlKFFAoFxLCg0g0sT2DCBscjmAKhOCU2NhbmMxRFyWQyTA8gEkD/EIvFAtMeHo+HKQ3idyiKQvaVWCyGPCH4F4QS6HQ6eMcRCAQCgQDKAZFIBE9bBEH6+/uhTS8mHZFIJDKZDMKVUCicP3++Xq8/duwYfEYo29br9bdu3Zo4MXz22WeZmZlz586d/jI6udkDN0ZR9Gc/+9mtW7eg0R68SKfTMzMzoTgOFBRIp4XMG5gjQa4AqopdMVBKmEymWCxOTk6mUql0Ot1sNhOJRJ1OB+zTZDIBPcWCoSA2yGQyIFXASCBjF+KzTCbTaDQGBATAQaHYikAgsFgsnU4XGBgIbVjy8/PZbDaNRiORSCqVis1moygK8SmdTkehUOBLBI0WYyTYz3v37rW2tubn509qtuQ8TCZTZWXlwMDAypUrPZVqKRKJKisrzWZzaWmpp9xGZqOICUttJpFI69atm2nb4OnB5XJv375tNBrnzp07S7wHRdHS0tLZmJ4gV1Wv10/aG/tHjB8AAQKMjY0dO3YMRdF9+/Z5PI6DAfpYtbe3l5SUzJ5D47gjQkPHqKioLVu2eEranR7gFLJ161bPDgEmk+nNN99cs2aNYyfLmeLmzZs8Hm/v3r0uvBdrquVmngqfzz969OjKlSs9dQ+0t7dfuHDhscceCwwMBAlKoVBA4ifU14AoBRtTKJSxsTEqlapWq/F4vMlkgiiJTqdraWlJTU0lk8l4PB6iezCpm81maIEkl8uJRCIYJcB7wYqaxWKBMmG3281mM1RR2Ww2tVpNp9NVKhVmZ4IgCOSpJCcnT/wgmZmZXC4XvLv+8Ic/QOP0n/3sZwsXLly7du2vfvUrBEHsdjtsRiAQVq9eLRQK8/LyGhoaWltbz5w5AxsrlcqgoKDW1lbwDoB9TupdBO3tdDqd4+kBjSCRSEBEQLlBEMRoNDKZTDKZDFFFCH3icDgKhQKRTUgjYzAYFAqFxWKBcAJ1UkBDw8PDIe4JZd5ARisqKtavXw+BBqhxIxKJEC4BZWWq710mkx05cmTu3Lkusxa73X7x4sX29vbS0tJpnLKdwfDwcEVFhV6v3759e0REhDu7AmB8IigoaPo0NecxODhYVVUF3X48pUBAFk5LS4vNZisvL/csgeByuXV1dUqlMi8vr6SkxIPJQ5Dfw+fziUTibPAeqIS9deuW3W6Pi4tbsWLF7IVZvrf4wRAgwLdDg5AJjWlm70AYBALBmTNncDjcunXrvoVEacgRgXolz17JL774gslkPvzwwy7voaWlpbm5eaq+BA/EZ599FhMTU1xc7PIJAC5cuHDv3r39+/d7xODLZDK99957sbGx69atc3knNpvtnXfeyc3NzcnJ0ev1EOMD3QIifSBQIf/Xbba6ujotLQ3mPGASEHjVaDS3b99esWJFbW3t8ePHYWMqlRoVFcXj8SY9ga+++kqr1WIJsMuXL3/kkUdOnDiBIEhaWtqFCxf+8pe/HDlyBDbjcDjFxcUSiWTBggUVFRV//etf9+7dCxu//vrrERERer0+MzNzYGAA9rl+/fpJD/rkk09evHgR+mwjCEKj0fbv379r1y7IjIHMYjabfe7cucLCwumbpjU2NorF4g0bNjh5wVtaWmpqap5++umZNiQ3mUwff/wxgiCPPvqoaz13URT9+uuvhUKh+1JNU1PT1atXIyIiNm/e7HIDYEfweLyqqiqbzbZ27VqPyDMikejChQtarRZSfDxSxGS1Wq9fv97W1obD4RYvXvxAR2znYbPZWlpaGhoaoBUGBJE9tXOdTnflyhWBQEChUJYvX+5xCx+j0Qi1t2azGTK+fzp6z0T8wAgQQC6Xf/7553g8/pFHHpnV3Mz29vbq6mpIgHe+0ak7MJvNZ86cGRwcDAoKevjhh2f7oBqN5siRI35+fjt27PDgY3Djxg0Oh/Pkk0+6nD7S2tra2NjoMgf6+uuv7Xb7li1bXHs7BqPR+Mknn9hstj179ngkW+vOnTvnzp3btm2bOwUsZ86ckUql+/fvd5KaW63Wv/71r88///zEFd6RI0e2bt1aWVl58ODBsrKyRx99tKysLDMz886dO5Mu6Ht6eoqKisCxsLe3VyQSnT59eu/evZmZmY888sjLL7+cmJi4c+fO0tLSoqIiFosVGxsrEAiKioqOHTv2zDPP5Ofnw8Zms5nJZCqVytDQ0DfffBP2OZUV9erVq19++eV58+ZB0X5ERERxcfHEjqfj2tpPxPDw8FdffQUylTM4deoUgiDOsyUMoNk8+uijruUUgoeFQqHYsmWLO2shlUp18uRJuVy+cOFCj5iyDg0NnTt3TqfTJSUlrVq1yv0RQygUVlVVKRSKiIiIdevWeSQH2WAwVFdX8/l8u92em5vrwWgUlBlyuVyTyTR//nzPUgc+n19TU6PT6VgsVmlpaWxsrKf2DJBKpZcvXx4dHaVQKKWlpZPquz9B/CAJEMBoNJ4+fZrP5z/88MMPbCXjDlQq1enTpwcGBvLz87+1Riez7evgCJFIdOLEidjY2A0bNnhqsBCLxYcPH55pT0dHdHZ2Xrp06ZlnnnHtlO7cuXPhwoXf/va37kcVNRrNiRMnHtga3UmA9sbn85966imX004HBgY+//zzX/ziF1CC9EBoNJo333zzpZdeGncx79y5o9VqwYgPu1Bvv/327du3jx49+n1YGr777rs1NTUnT550fHHSlu/TEyCj0fi3v/3tT3/6kzOkHC7Xxo0bZxrMbW5urqqq2rBhg2tZfdBmBEXRXbt2OfnNTgoOh3Pt2jUSibRt2zb3l4gKhQLCUkFBQWvXrnWfpgwODl66dEkul3uQ94yNjV28eHF0dJRIJC5fvtyDE7xUKq2urpbJZEQicdmyZR60gYDknubmZjweHx4evnLlSs8WI2u12vr6+t7eXgRBgoKCli9f7qmivB8NfsAECGCz2c6dO9fd3Z2WlrZq1apZrViBPro4HG737t3fgnsQAJw9pVKpB4Pik6Kjo6Oqqio1NXXt2rUeUXRRFP3Pf/6TkJCwevVq1/YgEAiOHj367LPPujYuSCSSw4cP79q1yyOpTs60RnceSqXygw8+iIiI2LFjh2t7MJlM77zzTn5+fkFBgTPbCwSCqqqqcaIaxNSeeuqpcRu/9NJL165d+9e//jXbHfSmgVQqfeGFF8Cwx5EpGo3Gzz777MCBA44by2Sy6urq7du3T7W3V1999fHHH3fGyKSzs/Ps2bNPPvnkjDQ/kUj02WefZWdnj/OhdhK9vb0XL14kk8m7d+92ORUDlmpCoXDRokXujxUGg+HChQu9vb0+Pj6bN292f+4cGBi4fPmyXC4PDw/fsGGDm0VnAIFAcPnyZYPB4OPjU15e7sEyXkhq1ul0AQEBa9eu9aAYPzo6Wl9fLxKJEAQpKiqCPtOe2rnFYmlvb29ra9NqtTabbfHixR5vffpjwg+eAGHo6empqqoymUzbt2+fvq2jmxgdHT158qRGoykpKfk2DRJu377d0NBAJBJXrFgxS2Z0CILcv3//7NmzERERGzdunGnew6RoamqqqanZv3+/a77per3+jTfe2LhxY3p6ugtvt9vtn3zyCXBWj7C67u7uU6dOQVqr+8NWa2vrhQsXNm3a5HJF6+nTp8HhyRnqf+PGDalUunHjRscX33jjjXFm0IDbt28/99xz6enpTzzxxLdsBCISiY4cOXL8+PHf/e53u3btGned79+/LxKJxmmxx44de+ihh6aqZ/7oo48KCgoeeJHtdvvRo0dpNNrWrVudP9vR0dFPP/00ODh427ZtM12A2Wy2S5cucbnctLS08vJy1+4oFEUvX77M4/F8fX03btzo5lStUqkuXrwoFouZTOaaNWvcdKO12+2tra23bt1CUTQhIaGsrMz9UcVisTQ2NkLj9Pj4+OXLl3sqe1epVNbV1Q0MDOBwuNzc3Ly8PE+tqC0WS0tLC7TfCQkJWbFihQdbh1qt1vb29sbGRnBAyM/PB0NLT+3/R4wfDwEC6PX6ioqKkZGR+Pj48vJyj2T8TQq73V5XV3fnzh0qlbpq1arpUy89CLPZXFNT093djcPhioqKZondDw0NnTp1CkXRZcuWuV+iD7aWUql03759ro3OlZWVfD7/iSeecG307O/vP3bsWH5+PtZ3zE20t7dfvnw5Ojp648aNbobY7Hb7hQsX7ty5s337dteUKolE8sknnzz88MPOcMTLly9rNBpHs6WmpiaDwbBkyZJJz62ysvLIkSNisXjHjh1lZWUZGRmzp0GKRKKampovv/xSoVDs2bPn0UcfnVSDqaqqSk5OHnetpqJxCIJcunSJRCKVlJRMf3SlUnn48OH169c7Hz3R6XRffPEFpIjNdA7WaDSnTp0Si8XuuGBwudz6+nqTyVReXu7mokgmk125ckUoFLLZ7LVr17pZJgbtdTs7O/F4/Lx58xYvXux+EH9kZOTmzZtCodBoNObk5HgqBQeTTKxWK4lEWrp0qafyb+x2+717927dumU0GolEYm5ubm5urqcYlVKpvH37dm9vLxhmZmdnL1q06NspIv4x4cdGgDD09/efP3/ebDZDd+vZo8Mmk6mqqqq/v5/BYJSVlc2GKeeksFqtN2/e5HA4CIKkp6cvWbLEI4KNI2w2W3V19d27d1kslvsfTa1Wf/zxxyQSadOmTS6oQXK5HBowbdiwwbVB5ObNm9evX1+8eLGTMaMHYmBgoLKy0m637927181lt81mO3PmDI/H2717t2tr7srKyp6enqeffvqBM83ly5fVavXmzZuxV1555ZXf/e5307xFIpF89dVXNTU1XV1d8+fPz8vLy8zMnFHP80mhUqmgF1hLS8vNmzf9/f1LSko2bdo0PSF4//339+3b5zjWQ6uZSa0TWlpa7t+/v2vXrunPpL6+vqmpyZmrB0BR9Pjx46Ojo48++uhM+wP09PTU1dVBkr5riT4DAwOXLl3ySFicz+fX1taOjY1FRES4H+iBTF65XM5ms8vLy92vZlWr1fX19Xw+H0GQiIiIZcuWeaQbg8ViaW1t5XA4YB2Zn5+fnZ3tqZL7e/fuNTY2gitmZmZmfn6+p0ZmPp9/48YNcNRkMBj5+flgEO+Rnf808aMlQBjAp0GtVmdlZa1cuXL2ODJWfWCz2YB1fWu3Jp/Pv3nzplwuRxBk+fLlrkWLpgH20Ugk0qJFi9yRnbDUddea6gkEgtOnT9NotD179rg2rEAiV1BQ0NatWz2ygoTZVyaT7dy50811M0yrY2NjrrWoGxkZOXz4cHl5+QO/IDB7xTjQtWvXcDicM8YBVqsVfMzv3r3b1dU1NDQUGBgYEhISEREREhISGhrq5+cH3dzgJ5TlQ4m+TqeDbvDDw8NisVgqlTKZzMzMzPT09Pnz5xcUFDhpJPOvf/1rXNLS559/vmzZson5HwKB4Ny5cxPTpR1hsVjAWtpJa1CXqSrm+JWQkOAaiZfL5efOnZNIJAkJCe4URthstjt37jQ2NqpUqszMzLKyMnfkGaw2SqfTZWVlrVixwk3d3WAwcDicnp4esAxdsGCBR2pcQOlpaGiAPiGZmZnFxcUeCRGgKNrR0QFpN1Ag5vzNPD30en1ra2t3d7dWq7VYLMnJyQUFBR6MnXnx4ydAGNra2q5fv261Wv38/JYtW+YRK7BJYbVar1271t7ejsPhAgMDS0pKZq+r6zjodLq6urq+vj4EQfz9/YuLiz1bTomiaF1dXWdnJw6HCw0NLSkpcTk1sqGh4fbt2wQCYcmSJTPNtB0aGoJC91WrVrmm/Pf09Jw7d45IJK5evdoj4UuDwXDixAmZTJabm7tkyRJ3uK9Op/vss8/0ev2mTZtcWENXVlb29fXt3bt3enWhpqZGIBDs378f/vzb3/723HPPucApFQqFI6eBnt7QaQQ2gN4LWq322rVrv/71r4ODgyMjI4ODg0NCQly7ShMJ0Guvvfab3/xm3GZSqfSDDz74wx/+MM1Rurq6vvrqq8cff9yZ5Fmr1Xry5EmBQLB58+YZ3TNglIyi6MqVK11I9pJIJBcuXFCr1WCN4bJrl0Qiqa6uHhsbw+Px+fn5ubm5Lt+lOp3u1q1b9+/ft9lsdDq9pKTEzToDqVR669atgYEBPB5PpVILCws9soobGBhoaGgYHh7G4XDApebOnet+EMpisXR2dnI4HI1GA81cFy1alJWV5f6eoVvc3bt3odsanU5fuHBhRkbGd9WO5kePnxABwqBQKGpqaoaGhhAEiY2NLSwsnD0zodHR0draWngC/f39CwsLPW7wMBW0Wm1DQwOXy0UQxMfHBx4kD+5fKpWCOm21WpOSksrKylwYmjGTVqPRWFRUVFRUNKP3wpI6ODh4y5YtLhSVYFa2Lu9hIjgcTkNDA4lEWr9+vWt53wAURc+cOXP//v2SkpL8/PwZvddsNh87dkyn0+3bt2+aNW5PT8+xY8cOHTpEo9GUSuVHH3303HPPuXzCD8TmzZt37tzpgq2OI4aHh69fv+7o8NTS0qLRaMb14lWpVG+99dYf/vCHqUQ+MCOwWCwPjI4hCIKi6MmTJ/l8/pYtW5w3cHLTKBnKPw0GQ1BQ0OrVq10LTlksloaGhqamJhRF/f393UlqhvGkt7cX2vcuWLDAzfIlPp9/+/bt0dFRBEECAgIKCgrcX4rodLo7d+50d3dDG5mAgIAlS5Z4JC1BLBZzOByJRALNcVNTUwsLC92xKgBotdq2traOjg6NRgONBbOysnJycn6CpszfCX6KBMgRAwMDN27cGBsbs9vtnnoIp8LY2NitW7cgnk0kElNSUubPn+/+I+QMJpIhzwaP+Xz+1atXoS1lTk7OkiVLZhpqxJiQwWAoLi6eERMaHh6uqqoaGxubO3euaykREomksrJSr9cvWrSooKDA/YwxSHEdHR1NTExcuXKlOyEGMCWPjIycaZGRRCL59NNPwc18qm3kcvlbb7118ODB0NDQGzduqNXq8vJyl0/1geezfPnyhoYGd4jmF198UVpa6ijYvPrqq88//7zjNnq9/n/+539+97vfTSVo9ff3f/HFF7t3736gwIaiaEVFBY/H27Fjh/MiB3SURBCkvLx8pk6+PB6vurraaDTGxMSsWrXKtUgKPI8ajYZMJi9cuNBlpqLRaBobG2HccHNXCIIolcqmpqb79+/b7XYcDpeSkpKfn++mDT2Kovfu3QM7K7vdzmKxcnNzPSKZDAwMcDgcKFbH4XBJSUkLFixwP/yk0Wja29vv3r2rVqsRBGEwGHPnzs3NzfV4BqcXzuCnToAcAQSlv78f/oyLi8vNzZ2l6JXVah0YGOByuRKJBEEQrVYbGBiYk5PzLSS1qVSq69ev9/b24nA4Op2em5ublZXlKcs7i8XS3NzM4XCghXV2dnZubu6MAu2wbr537x4ej8/JyZlRJeqNGzcaGxupVOqKFStcs1puamq6desWgiAFBQULFy50YQ/jAOt46K7sTlrY/fv3z507x2KxNm7cOKOwY1tbW1VVVW5u7vLlyyfdwGKx/Oc//0lMTFy1atUHH3zgqejDpHjjjTckEslrr73m8h7efPPNZ599Fvuzu7uby+U6FvYbDIa///3vzz///KQVZFDojiDInj17pv8uzGbzl19+KRQK9+7d62S4XKFQnDlzZnh4OCEhYd26dc6vASwWy40bN9ra2ux2e1JS0vLly13ITRGLxfX19VKpFEGQlJSUxYsXu0aeRkZGbty4IRKJYE1YVFQUFxfnwn4QBLFYLN3d3e3t7dAnOCAgID8/Py0tzZ0FBoqiPB6vvb0dGhIbjca4uDj3xRiz2Xzv3j1QYnA4nN1uj4mJycnJcTNTwmq19vX1dXR0iMViUM7YbPa8efPS09O/Dy6jXngJ0JTo7+9vaWmRSCTwPISFhWVnZyckJMxSOFYgENy5cwc6hyMIEhwcnJaWBj2rZ+NwAEjfu3Pnjs1mMxgMfn5+xcXFLo9342C1Wu/fv9/Y2KjRaKxWa1RUVElJyYyWUFwu99atW3q9nsFgrF271kmXMyz0QCaT16xZ4xoTgmbUZrN5poG5SWG322/evNnU1EQgEPLy8hYuXOjyivz8+fMjIyNgJef8GNrR0VFTU5OUlLRmzZpJN2hoaLh06dKhQ4fef//9jRs3zlKg1mKx5Ofnf/DBB042ih8Hg8Fw9OjRxx9/HHvl5ZdffuGFF7CLqdfrX3311eeee25SXaG7u/v48eM7d+6cXstRKpXHjh0zGo07duxwJohpNpsvXLjQ19cXFBS0bt0652NVcrm8pqYGJOHCwkIX+DEEBMViMY1Go9FoS5cudcECDTJtOzo6IDs4JCSkqKjINSs1SI65c+eOSqVCEIREImVnZ8+bN8+dgI5Op7t7925HR4fJZIKEm6ysrOzsbDezjJVK5b179yDbhkaj2e32hISEuXPnupkOoVAouFwulsQDwbKMjIy4uDhvudb3EF4C5CyGh4c7Ojp6enqsVisej8fj8XFxcSkpKTExMbNBiUZHR7lcbk9Pj8lkgiNSqdTExMTU1NSp3N7ch2Oc27PlDAiC9PT03L59Wy6X4/H4wMDA4uLiqKgoJ987Ojp67do1oVBot9uLi4ud1OGNRuOVK1c6OzuJROKGDRtcCG5ardb6+vr29naQo/Lz893/rrlcbkNDg1KpTE5Odrkssaurq7q62mQybdq0yXnC2tzcXF1dPX/+/BUrVky8gCqV6t13312yZEltbe3Pf/5zdxKYpj+Hp5566ubNmy7IAJWVlXPmzMG+x7q6OjKZXFhYCH8C+zl06NBE33CLxfL5558TCISdO3dOc+dIJJIzZ84YDIaf/exnD6zCs1gs1dXVnZ2dZDK5vLzc+bvr7t27169fN5vNfn5+K1eunKl5sUQiuXHjBp/Pp1Kp4eHhS5cudWHOxgqqaTSazWbLzMxcsGCBC7ITaDwwaCAIgsPh0tLS8vPzXaYRNpuNz+d3dHQMDQ3hcDibzcZgMObMmTNnzhx3Krb0ev39+/c7OzuhhhxBkKCgoOzs7OTkZHeeaLPZ3NPTw+Fw4Erq9Xp/f/+cnJzU1FRv2vIPAl4C5DrEYnFnZ6dQKFSr1UajkUql0un0uLi4zMzM2eAoJpOJx+N1d3ePjIzYbDYcDofD4cLCwpKTk+Pj4z1FUzBA1SiHw4Gm4r6+vvPmzUtLS/OIciuVSsHTzG634/H45OTkBQsWOCNim83ma9euQVJCaGjo0qVLnZmn9Xr9hQsXRCIRmUwuKipyrbc2Zo3v7++/bt06N3MXYIe1tbV6vT49Pd21UmTwFBCJRGFhYc43aero6Lh8+TKLxZrUvq++vr6+vt5oNO7fv98jLUQm4uc///nChQvH9bJwBq+//vqhQ4fgdxRF//GPf2D2RVqtFmrBJiYYtbW1nT59+sCBA9OEM9ra2mpra8PCwjZv3jw9JcWS1fR6/cqVK530yIZOHRCfyszMLCoqcn46t9lsXV1dTU1NIKvExsYWFxfPNBlleHi4ubm5r68Ph8Ph8fjY2FjXCqrlcnlbWxtkGdvtdhKJlJWVNXfuXJdbBatUqrt37/J4PIvFYjAYUBRNSUlxUzIBXtLZ2SkWi+EVAoGQmpo6Z84cd2i9SCTq7u7u6+szmUww/NLp9KysrPT0dI8Pv158O/ASIE9Co9EARxkbG4MnBEGQyMjIlJSUuLi42bClHhoa6unp6e/v1+l02ItUKjUuLi42NjY6OtpTvkdqtbqjo6OjowNFUVCMMzIy5s6d65HKKbFYfPPmTYgF4HC4+fPn5+fnP3BIxcrQYClcWlr6QAqFdR80m83x8fGrVq1y4fwFAsGVK1e0Wi2DwViyZIn7bUk6Oztra2tRFI2KiiotLXXBAWhsbOz8+fPQOnf9+vXOfKihoaGvvvqKQqFs27Zt3BH1ev0HH3zQ0dGxb98+TF/xIORyeWFh4bVr12aUzDQ0NHT9+vVt27bBn++8886WLVtAPhkdHX3nnXd+85vfjLtnNBrN4cOHU1JSpgr8Id/kmEOq+DQzLvCe5uZmk8nkZFQUlhDt7e0ajYbBYKxatcr5hBLH4iA8Hu9CzZFQKGxraxMIBHDyYWFhixYtmmk9lMViGRwc5HK5fD4f1nj+/v7p6ekZGRmuheZlMhkI2waDAY/HQ5pRVlZWamqqyyOVUqnk8XhcLlev11utVgKBQKPR0tPT09LSXPZ1NJvNIpGIy+UODQ3h8XiDwQBB/Pj4+PT0dG+J1o8GXgI0u7Db7UKhkMfjCQQCk8mEvUgmk6OiomJjY2NjYz3+OBmNxv7+/v7+/sHBQRRFsdf9/PxiY2OjoqLCw8PdJEaQigjSN5VKhTXlggUL3Je+oKyjtbUVKB2VSk1PT58zZ870iotAIKivr4clcmpqalFR0QPXZDwer6amxm2vGPYAABrfSURBVGQysdnsZcuWOR+PwwAWcPfu3cPhcBERESUlJW5mYoLZ/+joqMFgyMjIcKGBgFAovHjxolwud9IrT6PRHD9+XKlUPvTQQ7m5uY7/6u7u/vOf/5yenv7iiy/O+JM8CJ988kl9ff2RI0ecf8s777yzf/9+eFja29v7+/vXr1+PIIhYLP70009//etfj7ulL168yOVyH3vssUnvHJ1O9/XXX8tkspUrV07f8p3L5V69elWn0xUUFDwwTcdms0FXJpPJRCKRCgsL58yZ46SS4ejzSyAQ8vPz582b52SgEEJRWIgHKsBzcnJmlHFsNpt5PF5XV9fIyAhIs1CrlZ2d7VpICxJienp6oOIJQRA/P7+MjIy0tDTXJBOTydTb28vj8UDasdvtdrvdz88vJSUlNTXVtZbJFotFKBT29fX19/ebzWagUBQKJTExMS0tzYPdVb34HsJLgL4boCgqFArFYrFEIoExy2AwQMuY6Ojo+Pj4xMRED/YfBigUCoFAIBKJRCKR1WpFEMRms8E6zN/fPzo6OiIiIiIiwgWlymaz9fb23rlzB0R+HA7HZrNhEeamOIwxLYiyY+Vy04zsXV1dN2/eNBgMCIIEBQUVFhZOT26USmVtba1QKEQQJDIycunSpS7wGOhjpVQqcTgcrNddDgog33j1QsSNQqGUlZXN1ECPx+NduXIFRdGIiAhnLJo4HA6Hw6HRaBs3bnQ888OHD3/99devv/66a3HDqWC325ctW/biiy+O8++ZCiaT6f333wf/Q5PJ9NprrwEt6+rqOnfu3HPPPed4M/B4vK+//vrhhx+elNnw+fyzZ8+azeYdO3ZMJclYrdZbt25Bn5n8/PwHJq1jNhDQ+spJi2FIucPKTlNSUvLy8px86kFs5vF4mPNkWlraokWLnA/xAJnA4kRYInB8fLwL4SelUtnT08Pj8eRyOSZ+h4WFZWRkJCQkuBA3t9lsg4ODPB4PeAnsk0KhJCQkJCcnu1acBboOn88XiUSQpIwgiNFoDAsLi4+PT01NdeeZ9eKHCC8B+n4BiNHAwIBAIIA6T8gEJBAIIDOEhIRERUWFhIR4topybGxsaGhocHBwaGjIYrEAPSIQCBaLxdfXNzg4OCIiYkbuvSqVqre3l8vlYlmHdDo9LS3N/ajZ4OBgR0eHQCCA9R+0U0hNTZ10jh8ZGWlqaoK0bgiuTZ/sidUA63Q6l90d+Xx+Q0ODVqvVarVRUVHLli1zIaqFQaPR1NXVwTQZExNTWFg4o8gR1uTS19e3vLx8ei44MjJSUVFhMBgKCgry8vLgRZVKdfDgwcDAwKeeemqmrjbTgMfjbd++/datW3V1ddO0odBqtRQK5csvv1y2bFl4eDiCIK+++uqBAwcCAgKam5sbGxsdO13o9fpPP/3U19d369atE29UiHZFRUVt2bJlUgUUJD0+n28wGHJycqbRe+x2e1dXF3zLCIJkZGQUFBQ8UMqFfBdo8oAgCI1Gmz9/vjPOF2azGQSekZERAoFgs9lYLFZ2dnZGRoaTGqFSqezr63N8HhkMRmxsrAtpMTqdrq+vr6enZ3h4GPlGhvH3909KSkpKSnLhVlepVHw+n8/nQyY1gEAgREdHJycnx8XFuTDWSaVSgUDQ19cHRosAGo2WkJCQkJAQGRnpbZbuBeIlQD8gWCyWkZERkUg0PDw8NDRkt9uhOgy+QRqNFhYWFhERERYWFhwc7MGSS4PBAO0OQK8CZQWOiCAIiqIsFis8PDwsLCwyMnKa9atKperq6urs7NTr9SCwM5nMjIyMpKQkd7QunU7X2dnJ5XLBCQ1OLCUlJT09fVyOp9ls5nA4bW1tsKAMDAycP39+cnLyNErSjRs3jEYjmILk5+e7kEHJ4/Fu3rypVqtxOBz4oLjjtAkVQLBkn2k3RI1Gc+nSJaFQSCAQFixYsHDhwmkKVTgcTmNjo81mW716NVSZnT9//uOPP87IyNixY4fz/dKnx2OPPXbp0iW5XK5UKqea5D788MN//etfixcvPnjw4KVLl+Lj42k0WllZWV1d3fDw8Pbt27Etq6qqurq6Dhw4MI6zajSaioqKkZGRsrKySZuujI2NXbx4cXR0lEQilZaWTpXUhbUNhwZnKSkpxcXF0/fmHBsbg15ONpsNQZCgoKDc3NzExMTpvzJYBWEmYQiCgDbjZPdZjUbT19fX29sLBAWHw1mt1uDg4NTU1JSUlBkF3MViMZ/P7+vrgxsYXmQwGAkJCUlJSTP1SAMBpr+/HzKKsNfZbHZ8fHx8fDwQXOdhtVrFYjEsF6ETIjzLNpstODg4Li4uISFh9lz+vfgRwEuAfiQwGo1isVgsFg8NDcGix2q1QsWH3W632WxMJjMoKCgoKCg0NDQ4ONhTZQsWi0UqlUokEolEIpVKtVotrFDtdjs2v0K/zKCgoMDAwKCgIIxwqNXqrq4uHo+nVCqRbwavoKCgxMTEpKQk1yL6CIIYjcb79+/fv39fKpXCyZDJZBiyo6KisKM7lp8gCEKn0+fNmzdV/SpEKyBNGyxxnc/PwCCTyRoaGvr7++12O5FIzMrKmj9/vst6GLj09vb2Qm5ZTk7OokWLnEntgvgaqBdsNrusrGyqQnqDwXDmzJmhoaGAgIC1a9dSqdTDhw8PDg6GhIQUFBQ40zx1Ghw+fPjpp582Go2BgYFXrlyZyhnoqaeeevvttyE0LJPJkpOTTSbTc8895+/vn5+fD2JYb2/viRMn1qxZMy5O197eXl1dzWAwtm7dOjGs2d/fX11drdfrfXx8Vq1aNWn62sSLvHDhwqlEF5vN1tfX19bWhtWEBwQEzJs3LyUlZZpbZZy6Y7VaaTRaUlJSZmbmAydvhUIB2gmm68Bx4+Li4uPjY2NjnbxFh4eHBwcHBwcHocIUez08PDw+Pj4hIWFGOqharR4cHBQIBEKhEEVRjJSQyeTIyEjgOjPKoTYYDIODgyCNGwwG0JwQBCESiREREdHR0TExMd+Opb4XPzJ4CdBPBQaDQSaTSSQSmUwmk8nA8BSjRzDr+/n5AU3x9/cPCAjwSIWXzWaTyWSjo6MymWxkZARTpDF+Bn8GBwcHBQXBKCaXywUCAdS/wH9DQkKSkpISEhJcOyUURfl8Po/HGxwcBBt+BEEoFEp8fHxSUlJkZCQOh1MoFBwOB+p7kW8CdpmZmRPVKZgUOzs74c/IyMicnBwX6msg1xtiKFg5sWsf0Gg0Njc3t7W1gT9CVFRUbm6uM2ndGo2mtra2v78fh8NBffWkE4lCoaiqqhIKhUFBQXl5eZcvX4aQYkxMzIYNG1wrCDIajdu3b6+trTUYDG+++eYvf/nLSTd75JFHvvrqqxUrVojF4rt37xKJRIvFUlhYeOrUqbKysn379tntdh8fn927d2N3i8lkOn/+fE9PT3x8/MaNGx0ZrcViuX37dmtrq91uT0xMLC0tnaiIQE2iMzIblFtjTsdEIjE5OTkjI2MqkQ/yWqBsExgVHo8nk8nO0J2JXMdgMNBotPhvMM17sbMViUSg44rFYqjqotFomIgbHx/vvKuZ0WgUiUQSiWRwcBDyfjBe4uvrGxMTEx0dHRkZ6XyxxdjYmFgsHhwcFIvFYH5GJBJhlGAwGOHh4dHR0XFxcd4KLC88CC8B8uL/QyaTjY2NyeXysbGxsbExmJsRBMFuEpgG6HR6QEAAkCT46WZNmcViGR0dHXMAFmiDI2o0GqlUOjo6SiAQ6HQ6g8FgMBi+vr5JSUmxsbFxcXEuCFomkwlSGaDdD4BEIkVHR8fGxkZERAwODkJZDZwMUKL58+ePO5ZIJOJwOAMDA/BndHR0Tk7OTGvKUBS9e/cuhFcQBCEQCImJiVlZWTMNCgBArwLDEqvVGhkZWVRU9MA5cmBg4Pr162NjYwiCREVFLV68eOKUPDw8fP78eY1Go9VqDQZDQkKCyWTSaDTLli1z0hHHEXa7/a9//et///d/l5WVnT59Gl6ExHxsm7y8vLCwsEceeWTv3r2BgYFisZjJZKalpVksFg6Hk5ycXF9fLxAIXnzxxerqah6Pd/HiRbPZvHXrVsc+X6Ojo/X19SKRyGg0FhYWjmv3hqJoW1sbh8MBUhIfHz+pR47FYunp6bl3755EIgGphsFgZGVlZWZmTnr7YZQFYxsEAiEkJGSaUuppYjrh4eEJCQmJiYkPVEYVCgWwHLAoQ76JVpvNZji68wUWBoNBIpHARwDiRaPR4HHA4/GhoaHx8fFhYWGhoaHOSE06nU4sFo+MjAwPD0ulUovFAnQW3uvn5xcdHR0eHh4REeFtjOXFtwMvAfJixtDr9UCS5HI5EBes2H7cQhmPx/s6wM/Pz9fX1501nEKhwKiSTCYbGBiAM4GaDriZCQQCSFkxMTERERHYcX19fR9YmzOxOg9gNpuZTKbVajUajZB6hcPhqFRqUlJSamqqY26QUCjkcDjQ0sRms02jJE0DCKa0t7dDCyHI9U5JSUlJSXHBvK6/v5/D4WCmcBEREdnZ2fHx8dNMWkKh8Nq1azKZDIfDBQcHFxQUjJO4dDpdbW1tY2NjR0dHRkbG0qVLBQIBDodbvny52Wz+5S9/eeTIkalYYH9/f3t7e1dXF6z4Ozs7RSJReHg4g8HQ6XS+vr5Go9FisUD4cnR0NCsra3BwMCoqqqurKzc39/e///2777577NgxNpsdGBi4YsWKCxcuyGSygwcP5ubmCoVCMEuESN+dO3d0Oh2DwVi9ejWWs2K323k8XnNzM5AMEok0Z86cnJyccfMu5nQKtze0nZoo8FitVqgtguxp5JunICgoKDk5OSkpaWJtkaMYM2kRaFhYWFhY2FSpQjqdTiQSjYyMDA0NyeVyYBJYRmBgYCCQ+NDQ0OnlHJvNNjIyMjIyIhaLpVKpSqXC9gNXyc/PD1wzIiMjH/jYoigqlUqHh4chGq7RaBz3Bm2wgoODw8PDoZzCm4bsxXcOLwHyYhZht9uVSqVCoVAoFMpvgKk7E0EikVgO8PPz8/HxYbFYM6rMN5vNQqFQJBL19fUNDQ2ZTCatVms0GvV6vcFgoNPpbDbbx8eHyWTCBAOBAKPRCFQJQ0BAAHZco9EoEAgEAsHAwACKoqD2W61WqVSq1+sJBIKvr6+Pjw+JRIIcpvj4eEiPNRqNnZ2dd+/e1Wq1MBmwWKzU1NTk5OQZ1ctotVooe5bJZMg3wkB4eHhKSkpiYuKM7IIc2xXRaDSTyRQbGztN4AYq6UDGgAYphYWF2HRotVqrqqq++OILnU63YMGC69evV1dX43C4hISEv//979Cm1G63t7e3V1VVXblyRSQSxcbGzp07Nz09PTIyMjw8PCoqSiKRJCYmTjojWq3W4eFhrEqIy+U2NjbKZLLMzMzt27cPDQ0dPnwY4qp5eXlKpdJut7/00kvNzc0JCQkLFy5csGAB7FYikbS0tIAVMoFAiImJycvLw+rpDAZDd3c3j8eTSqVYRBgui2NbNMi0G1dHbTKZQAtJTk52lGf0er1QKJRIJCKRSKFQQDIy7JlCoUDmSkxMzKTEAqQXYEhyudwxZRj5puIBKg9CQkKmYjlarVYmk4F0Ojo6qtVqseQ87ExCQkJCQ0NDQ0PDw8OnkVHtdjvsBALomBk9PAVQo4ogCNSKgibkjVV58f2HlwB58T2CzWZTqVQqlUqtVmu+gUKhgNgETF02mw3L04RB3Gq1AmFis9lMJpNOpzOZTDabTafTx43CKIoODw8LhcKhoSGZTAbyO5SkwfobyBBUJGk0GmgPOe4M8Xg8TCFMJhMmPIPBoFarFQoFiUSCrHClUqnVaqlUKovF8vHx8ff3z8zMTEhIiIuLw1qajI2NgZKEIEhUVFRKSkpsbOyMqJ5IJOLxeBDqgk8BRimJiYnh4eFOJnNAZOfu3bvDw8PwuchkMkzn0dHR40gJiqLQIAW+ETabDYXcBAJBr9dnZmYC34WN2Wz2tm3bfH19Kysr09LSVq1aVVZWNtNkqUlhMplu3rxZUVFRUVHBYrEKCgqgdbFOpyOTyWlpaSiKnjp1qr+/H7O8cqw5BwrY2dkJaTQGgwFaJWDJNHK5fHBwcGBgAFwhkG9EHRaLBUnB4eHh8AoWbwKC6Eju6XQ6EItxcg7IWhifwMgEJuQQCAQ2mx0eHg5vn7TQDGNIsLqAQ2NpPQiC6HQ6PB4PRwdCP5WkhKm5APj6gNzA3Q7v8vf3Dw4ODg4ODgwMDAwM9JTFvBdefIfwEiAvfgwAqqTT6YA56fV6nU4HahOJRLJarTCUE4lEWLiTyWSLxYKxGQaDQaPR6HQ6tHXTarVyuZxKpeJwOBqNhsfjYXLy8fEJDAwMCAiAaYBAIMDEo/gG2DwEEyGRSFSpVFKpVCaTyeVySJeBeAqZTCaTyXQ6PSQkJDU1NSIiwsfHZ2xszNEunM1mx8XFxcXFRUREOFnrbjKZIMcWm/gRBAEHOchLjYmJeWCChcVi6e/v7+7uxhKbbDabn58fdJ1zTAxSq9UcDqerq6u9vd1qtX7yySdwUKBNeDw+JCTk0KFDjz/++OxldXR2dr744ov19fUpKSlSqVQkEplMJiqVSiaTt2zZ8vbbb4NuhPU0sNlsUIuUmJioUqmEQuHg4CB4GMInRRDE398/NjY2JiaGxWKNjIxIpVIIFSEIggVAIUsXYjqRkZGhoaGgwwElgvsBIlwYKYG7gkQiMZlMICXwXkzAUzoARFOdTjeRiyAIgt2EgHESoFarVU7AVBcQxE5/f3/I53PHs8oLL35Y8BIgL7xAlEqlXq83Go0qlcpisahUKpPJpNfrwfFFr9eTyWSbzabT6cxms1wuN5vNKpUK/IRQFAU6RSAQoL49KCjI19c3JCQEFs3h4eEUCsXX1xdmPgRBdDod2CRCo6Wenp6BgQGJRALkDJvkLBYLSEpEIpHJZIJARaFQmExmaGgoZDiB9xKVSgU5aio3HQjhgaoBHAuiIUQi0dfXF2IWISEh0yQYgdUv5iwH4waWM/7nP//53Llzer2ewWAEBwcbDAYGg/Huu++uWLHC09/V5FCpVHv27GlubiaTyVKpFCOgy5Yt27NnD51O1+l00KEPcz/H4XA+Pj40Go1MJqMoCvzGMUiEx+N9fHyCg4MhzR+Px6tUKmA2Go1m0tMAzgellIGBgUwm02QyOVJkYMmTvtdisTCZTBqNRqVSsTisj4+PwWDQaDQqlQpuG/gdAq+Ob3ccyZlMJiS9sdlsLP3OY9faCy9+LPASIC+88CQMBsPQ0JBIJIKyXsjaNplM0OwatCgEQQgEAoqiFAoFj8czGAzIFkdR1N/fH4fDMZlMPB5vNpshKRsmb4vFgqKozWbT6/UWiwUoGkRPCN+ASCSClOXj4xMUFOTj4+Pr6xsYGAgxMgj3gCAE6VZUKhVFUaPRCIVdOp2ORCKRyWQKhQIxNQRBwsLC/L9BQECAY1QRRdHBwcGVK1f29/eD9uPr6xsaGpqdnX3kyBFnSvNGR0fnzJmzcePGd955x/mLnJmZ+d57703sSHrmzJl9+/YZDAZwi4Hr/MQTT0BIFM6cQCDABiAUUSgUiKsCxZmotAG1JRKJcNFgY6vVCr3qIECGCYo4HI5EIoHIB/Khj48PmUxmMBgkEgmyi4DHgDxJIBDweDxs7/g7EFmLxQIpaywWC6t8hN+97ce98MJ9eAmQF158l1AoFBCw02g0arVap9Pp9XqlUmk0Gh1nRKvVCqQH8pZAgbBarVQqFeZyFEUh11upVBIIBKVSabVaDQYDsCjscGQyGcQPIpEI/IZGo2GNlmCCh+MiCALMDCQu4EagrICpHTAkEolUWVkJOwd/7aioqKamJiczkOrq6pRK5YYNG8a9/tlnn2VmZk7ljlhZWRkYGDhpS/avv/567969QA3hlT179rDZbKvVajKZbDYbXFgGgwGp6yiK+vj4gOoDChwwQgqFotfrwccP8l2AG/n6+tpsNvgJZWsTf9LpdIPBQKVSTSYTjUYDUYfNZgOJ8fHxgRx8b7G3F158t/BkPykvvPBipoBIh2vvtVgsWEAEaBOCIBC2M5vNIFFoNBoWi6XX6ykUitVqVavVeDxeoVAQiUQsaQlKqQ0GA8QB7XY7iqIgOEFcD3gAlmqNpaRYLBZHyQRqhbZt2zaR/UA47+233/7444+ffPLJRx99FF7//e9/39jY+Mwzz/zzn/9EEOTUqVP79u3bt28f9HJvbm7etWvXsWPH1q9ff/Xq1fPnz0Pnr1OnTk11TTZt2vSLX/wCuA68QqVSIVmHTCabTCZwQ1Cr1cCBEATBQodMJhO4DoPBgKwaeAU2YLPZ3sptL7z4McFLgLzw4ocKIpHoDn9yH0C2AgICQFUC3+RFixZN3JJIJJ4+fTosLCwhIWHVqlXY6y+//DKCIB999FF/f/9f/vKXEydO3Lx58+jRo//1X/+l1WpFItEf//hHtVp94MABBEFu3LgBm01zSjgcLiQkZGBg4MyZM2+99VZHR8fBgwenUpK88MKLnzK8ITAvvPDCLeTl5ZWXl7/00ks4HO7evXsHDhy4fv36pGXSb7/9NgSkHN/b1tb273//+8UXX0xMTGQymY2NjYWFhf/4xz+Ki4s5HI5CoSguLl6yZAmLxVq6dClsduPGjdLS0urq6omHOHXq1JdffnnixAn4E1JzZuST5IUXXvxE4CVAXnjhhSfx2muvXb9+/ZNPPhknTeXk5PzqV7+C4BcWOPPs+HPmzJkXXnjhypUrk3Y29cILL7xwhJcAeeGFFx7G0aNHX3nllT/+8Y9bt279dvJmxGLxc889p9Vq33//fdcaqHnhhRc/NXhz+rzwwgsPY8+ePXV1ddXV1XPnzv3nP/85jQuf+7h169bOnTvLy8s3bNhw9uxZL/vxwgsvnIRXAfLCCy9mCzKZ7KOPPvryyy8DAgJWrlxZVlaWkZExlVuj8xgbG6uvr7948eL169czMjIOHjxYUlLipFm2F1544QXAS4C88MKLWUd/f39VVVV1dXVXVxeLxZo7d25GRgY0Q42Ojg4JCZmKFanVapFIBL3Te3t7Ozo6eDyer69vQUHBqlWriouLvW46XnjhhWvwEiAvvPDiW4VCoWhra+vq6gLLbJFINDIyYjabiUQiOAxRKBTMTjAoKCgkJCQ6Ojo8PDw2NjY7Ozs5OdlJl0UvvPDCi2ngJUBeeOGFF1544cVPDv8PgAol9dgEwhEAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Create a causal model from the data and given common causes.\n", "model=CausalModel(\n", " data = data,\n", " treatment='treatment',\n", " outcome='y_factual',\n", " common_causes=[\"x\"+str(i) for i in range(1,26)]\n", " )\n", "model.view_model()\n", "from IPython.display import Image, display\n", "display(Image(filename=\"causal_model.png\"))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 2.Identify" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Estimand type: nonparametric-ate\n", "\n", "### Estimand : 1\n", "Estimand name: backdoor\n", "Estimand expression:\n", " d \n", "────────────(Expectation(y_factual|x18,x1,x21,x8,x24,x10,x15,x23,x22,x5,x16,x1\n", "d[treatment] \n", "\n", " \n", "4,x2,x7,x3,x19,x20,x6,x13,x11,x17,x25,x9,x4,x12))\n", " \n", "Estimand assumption 1, Unconfoundedness: If U→{treatment} and U→y_factual then P(y_factual|treatment,x18,x1,x21,x8,x24,x10,x15,x23,x22,x5,x16,x14,x2,x7,x3,x19,x20,x6,x13,x11,x17,x25,x9,x4,x12,U) = P(y_factual|treatment,x18,x1,x21,x8,x24,x10,x15,x23,x22,x5,x16,x14,x2,x7,x3,x19,x20,x6,x13,x11,x17,x25,x9,x4,x12)\n", "\n", "### Estimand : 2\n", "Estimand name: iv\n", "No such variable(s) found!\n", "\n", "### Estimand : 3\n", "Estimand name: frontdoor\n", "No such variable(s) found!\n", "\n" ] } ], "source": [ "#Identify the causal effect\n", "identified_estimand = model.identify_effect(proceed_when_unidentifiable=True, method_name=\"maximal-adjustment\")\n", "print(identified_estimand)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 3. Estimate (using different methods)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### 3.1 Using Linear Regression" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "*** Causal Estimate ***\n", "\n", "## Identified estimand\n", "Estimand type: nonparametric-ate\n", "\n", "### Estimand : 1\n", "Estimand name: backdoor\n", "Estimand expression:\n", " d \n", "────────────(Expectation(y_factual|x18,x1,x21,x8,x24,x10,x15,x23,x22,x5,x16,x1\n", "d[treatment] \n", "\n", " \n", "4,x2,x7,x3,x19,x20,x6,x13,x11,x17,x25,x9,x4,x12))\n", " \n", "Estimand assumption 1, Unconfoundedness: If U→{treatment} and U→y_factual then P(y_factual|treatment,x18,x1,x21,x8,x24,x10,x15,x23,x22,x5,x16,x14,x2,x7,x3,x19,x20,x6,x13,x11,x17,x25,x9,x4,x12,U) = P(y_factual|treatment,x18,x1,x21,x8,x24,x10,x15,x23,x22,x5,x16,x14,x2,x7,x3,x19,x20,x6,x13,x11,x17,x25,x9,x4,x12)\n", "\n", "## Realized estimand\n", "b: y_factual~treatment+x18+x1+x21+x8+x24+x10+x15+x23+x22+x5+x16+x14+x2+x7+x3+x19+x20+x6+x13+x11+x17+x25+x9+x4+x12\n", "Target units: ate\n", "\n", "## Estimate\n", "Mean value: 3.9286717508727214\n", "p-value: [1.58915682e-156]\n", "\n", "Causal Estimate is 3.9286717508727214\n", "ATE 4.021121012430829\n" ] } ], "source": [ "# Estimate the causal effect and compare it with Average Treatment Effect\n", "estimate = model.estimate_effect(identified_estimand,\n", " method_name=\"backdoor.linear_regression\", test_significance=True\n", ")\n", "\n", "print(estimate)\n", "\n", "print(\"Causal Estimate is \" + str(estimate.value))\n", "data_1 = data[data[\"treatment\"]==1]\n", "data_0 = data[data[\"treatment\"]==0]\n", "\n", "print(\"ATE\", np.mean(data_1[\"y_factual\"])- np.mean(data_0[\"y_factual\"]))\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### 3.2 Using Propensity Score Matching" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "Causal Estimate is 3.97913882321704\n", "ATE 4.021121012430829\n" ] } ], "source": [ "estimate = model.estimate_effect(identified_estimand,\n", " method_name=\"backdoor.propensity_score_matching\"\n", ")\n", "\n", "print(\"Causal Estimate is \" + str(estimate.value))\n", "\n", "print(\"ATE\", np.mean(data_1[\"y_factual\"])- np.mean(data_0[\"y_factual\"]))\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### 3.3 Using Propensity Score Stratification" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Causal Estimate is 3.4550471588628207\n", "ATE 4.021121012430829\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n" ] } ], "source": [ "estimate = model.estimate_effect(identified_estimand,\n", " method_name=\"backdoor.propensity_score_stratification\", method_params={'num_strata':50, 'clipping_threshold':5}\n", ")\n", "\n", "print(\"Causal Estimate is \" + str(estimate.value))\n", "print(\"ATE\", np.mean(data_1[\"y_factual\"])- np.mean(data_0[\"y_factual\"]))\n", "\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### 3.4 Using Propensity Score Weighting" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "Causal Estimate is 4.028748218389918\n", "ATE 4.021121012430829\n" ] } ], "source": [ "estimate = model.estimate_effect(identified_estimand,\n", " method_name=\"backdoor.propensity_score_weighting\"\n", ")\n", "\n", "print(\"Causal Estimate is \" + str(estimate.value))\n", "\n", "print(\"ATE\", np.mean(data_1[\"y_factual\"])- np.mean(data_0[\"y_factual\"]))\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 4. Refute\n", "##### Refute the obtained estimate using multiple robustness checks.\n", "##### 4.1 Adding a random common cause" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "Refute: Add a random common cause\n", "Estimated effect:4.028748218389918\n", "New effect:4.028354731290574\n", "p value:0.42000000000000004\n", "\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n" ] } ], "source": [ "refute_results=model.refute_estimate(identified_estimand, estimate,\n", " method_name=\"random_common_cause\")\n", "print(refute_results)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "##### 4.2 Using a placebo treatment" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "Refute: Use a Placebo Treatment\n", "Estimated effect:4.028748218389918\n", "New effect:0.0012777161602278353\n", "p value:0.48\n", "\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n" ] } ], "source": [ "res_placebo=model.refute_estimate(identified_estimand, estimate,\n", " method_name=\"placebo_treatment_refuter\", placebo_type=\"permute\")\n", "print(res_placebo)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### 4.3 Data Subset Refuter" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n", "/home/amit/py-envs/env3.8/lib/python3.8/site-packages/sklearn/utils/validation.py:993: DataConversionWarning: A column-vector y was passed when a 1d array was expected. Please change the shape of y to (n_samples, ), for example using ravel().\n", " y = column_or_1d(y, warn=True)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "Refute: Use a subset of data\n", "Estimated effect:4.028748218389918\n", "New effect:4.026016493307437\n", "p value:0.4\n", "\n" ] } ], "source": [ "res_subset=model.refute_estimate(identified_estimand, estimate,\n", " method_name=\"data_subset_refuter\", subset_fraction=0.9)\n", "print(res_subset)" ] } ], "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.8.5" }, "toc": { "base_numbering": 1, "nav_menu": {}, "number_sections": false, "sideBar": true, "skip_h1_title": true, "title_cell": "Table of Contents", "title_sidebar": "Contents", "toc_cell": false, "toc_position": {}, "toc_section_display": true, "toc_window_display": false } }, "nbformat": 4, "nbformat_minor": 4 }