{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# A Simple Example on Creating a Custom Refutation Using User-Defined Outcome Functions\n", "In this experiment, we define a linear dataset. In order to find the coefficients, we make use of the linear regression estimator. In order to test the effectiveness of the linear estimator, we now replace the outcome value with a dummy produced with the help of a linear expression based on the value of the confounders. This effectively means that the effect of the treatment on the outcome should be zero. This is exactly, what we should expect from the results of the refuter." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Insert Dependencies" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [], "source": [ "from dowhy import CausalModel\n", "import dowhy.datasets\n", "import pandas as pd\n", "import numpy as np\n", "import logging" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Create the Dataset\n", "You can change the values of the hyper params to see how the effects change, as each parameter changes\n", "Variable Guide:\n", "\n", "| Variable Name | Data Type | Interpretation |\n", "|-----------------|-----------|--------------------|\n", "| $Z_i$ | float | Insrument Variable |\n", "| $W_i$ | float | Confounder |\n", "| $v_0$ | float | Treatment |\n", "| $y$ | float | Outcome |\n", "\n" ] }, { "cell_type": "code", "execution_count": 10, "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", "
Z0W0W1v0y
00.0-1.0862650.566131-4.817235-51.000055
11.0-1.4461551.0789492.89188926.582560
21.00.3871381.05235512.759312134.107854
30.00.4401140.6280571.89223423.833082
41.0-0.116181-1.6179848.54340977.761051
\n", "
" ], "text/plain": [ " Z0 W0 W1 v0 y\n", "0 0.0 -1.086265 0.566131 -4.817235 -51.000055\n", "1 1.0 -1.446155 1.078949 2.891889 26.582560\n", "2 1.0 0.387138 1.052355 12.759312 134.107854\n", "3 0.0 0.440114 0.628057 1.892234 23.833082\n", "4 1.0 -0.116181 -1.617984 8.543409 77.761051" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Value of the coefficient [BETA]\n", "BETA = 10\n", "# Number of Common Causes\n", "NUM_COMMON_CAUSES = 2\n", "# Number of Instruments\n", "NUM_INSTRUMENTS = 1\n", "# Number of Samples\n", "NUM_SAMPLES = 100000\n", "# Treatment is Binary\n", "TREATMENT_IS_BINARY = False\n", "data = dowhy.datasets.linear_dataset(beta=BETA,\n", " num_common_causes=NUM_COMMON_CAUSES,\n", " num_instruments=NUM_INSTRUMENTS,\n", " num_samples=NUM_SAMPLES,\n", " treatment_is_binary=TREATMENT_IS_BINARY)\n", "data['df'].head()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Creating the Causal Model" ] }, { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "INFO:dowhy.causal_model:Model to find the causal effect of treatment ['v0'] on outcome ['y']\n" ] } ], "source": [ "model = CausalModel(\n", " data = data['df'],\n", " treatment = data['treatment_name'],\n", " outcome = data['outcome_name'],\n", " graph = data['gml_graph'],\n", " instruments = data['instrument_names'],\n", " logging_level = logging.INFO\n", ")" ] }, { "cell_type": "code", "execution_count": 18, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "WARNING:dowhy.causal_graph:Warning: Pygraphviz cannot be loaded. Check that graphviz and pygraphviz are installed.\n", "INFO:dowhy.causal_graph:Using Matplotlib for plotting\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAV0AAADnCAYAAAC9roUQAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nO2dd5hV1bnGf9/AMDPSEUGaYMMCoibW2Etiwa4xltiSa48m0cTr5HoTE6+iEUuMmhhTRGMhMVgiGqKxd41KsWsYQCz0MjB93vvHWgcOyAwzMOfsU77f8/Bs9j77rPPNnnXe2ftb33qXScJxHMfJDiVJB+A4jlNMuOg6juNkERddx3GcLOKi6ziOk0VcdB3HcbKIi67jOE4WcdF1HMfJIi66juM4WcRF13EcJ4u46DqO42QRF13HcZws4qLrOI6TRVx0HcdxsoiLruM4ThZx0XUcx8kinZMOwFkzwyon9gNOA0YBvYBFwBTgjqoxo+cmGZuTW3hfyS/MTcxzi2GVE3cGKoFDAAEVaS/XAAY8BoypGjP6texH6OQK3lfyExfdHGJY5cRzgOuAclpP/TQDtcDFVWNG/zYbsTm5hfeV/MVFN0dI+xJt0I63Lce/TEWH95X8xgfScoD4mNjeLxHx/OuGVU7cqeOjKg7M7HQzez7pONrKsMqJO9fMmHzDJ7eclkhfMbP+ZvasmS01s+vWp611+Oynzey/svmZmcAH0hLGzCq7DBj+/QGnXV+eOjb7tjPp3Hsg/Y//OenHeu31bcoGbs28R2+k/tMP6NRjI/p8/ezyik13rASOTSL+bGBmAraU9FHascuBLSR9O7HAkqHSsC6tnSCJpf/+O9Vv/YPGxV9QUt6NsoFb03OPE8u79Bu2vn3lLGAe0EP+mLxO+J1uwnQZuPWUhvmf9FdzUwlAY/UC1NRE/Rcfo+YmUscaF35G2ZCRzH34l3TpvxmDv38PvfY+hbkPXl3StGzRocMqJ26U6A/itIiZdcjNTaxSOIS1fG8XPvE7lr7+ML0PPJsh37+PgWfdRsXw3aj5+LUSYH37ylDgnXwS3I66/h2Fi27C9D/hiu1obqR+znQA6ma9TfnQUZT2GbzKsc69BqCGWuq/+Jhee55MSWkZXbfegy4bDWX5+y90IpQMFSVmtq+ZfWJmF5vZHDP7zMzOSHu9p5ndaWZzzWyGmV1mZiWrNmE3m9liM3vPzA5Ie+F0M/tPfJyebmYnp732HTN718wWmtkkMxua9prM7Hwz+xD40Mx+Y2ZjV4v7ITO7KP5/oJn9LcY43cwuTDuvwszumHHt0VWzbz+nrO7zD1u8Fg0LZrP0jYn0PeLHVAzbHutcSklpOd1G7EfP3b8JoMbFc85u6Xqk0i1mNjb+XNPN7JD42h2EfnaJmVWb2YFmVmZmN5rZp/HfjWZWlt7Waj+zzGyLVHtmdouZTYzX9xUz2zzt3K/H38diM7uZUI2R3lZ7rr+Z2Q2xfywxs6lmNrLFC5lBXHQTpqRLxYguA7eibuY0AOpmTaN88LaUDd52lWNlQ0bQMHcGnXttTEnZynRel36b0jBvVimhRrOY2RjoCQwCvgvcYma942u/jq9tBuwDnAqckfbeXYGPgb7Az4AJZtbHzLoCNwGHSOoOfA14C8DMjgR+AhwDbAQ8B9y7WkxHxba3ja99y8wsvr838A3gvih4fwcmx/gPAH5gZgfFdn4GbD7onNv/3v9bV9iyqf9q8SLUzphMp+4bUjZwq5ZOqZj/2E1ntOF6vB+vxy+BP5iZSToduBv4paRukp4A/gfYDdgB2B7YBbisxQC/zAnAz4HewEfAlQBm1heYENvqS/j97JF60zpc/28AewPD489+PDC/HXF2GC66ydOrfMhIamcFga2d9TZlQ0ZQPmTEKsfKN9mO5oZaSsq6rvJmK+tKc30NhKL4YqYB+IWkBkmPAtXAVmbWifDFrpS0VFIVYdDylLT3zgFujO8dTxCc0fG1ZmCkmVVI+kzS2/H4OcAYSe9KagSuAnZIv9uKry+QVEMQBQF7xdeOA16S9CmwM7CRpF9Iqpf0H+D2GDcEgbiyc/e+G3TusRHddzqixYvQXLOUTt36tPi6mpuonTFl07VcjxmSbpfUBIwDBgD9W2jyZMJ1nyNpLkFAT2nh3DXxgKRX4zW8myDeAIcCb0u6X1IDcCPwedr72nv9G4DuwNaEqq13JX3Wjjg7DBfd5FlUNmQkdZ+8Q1PNUpprllDaZxBlg7ahbvZ7NNUspWHeTMqGjKCktBzVLV/lzapbTkmXCgizkAqVJqB0tWOlhC9Sivnxy5diOdCNcJdUCsxIe20G4Y4yxezVcpQzgIGSlgHfInzBP4uPwVvHc4YCvzKzRWa2CFhAePxNb3dW6j+x/fuAE+Ohkwgik2prYKqt2N5PWCl0A2NbiwA692g5JVtS0Z2m6gUtvt5cswTUbGu5HivETVKqw3VrocmBa2hrYIsBfJl0IU39zlLtrn79ZqWd297r/yRwM3ALMMfMfmdmPdoRZ4fhops8U8oGbV3TXLec6smTKBu0DQAlZRvQqVsfqidPolO3PpT22pjSjYbSsOhzmtOEt37OdDpvOLiBMO2zUJkJDFvt2Kas+mVviXkEcU6/A9oEmJ22Pyj12J/2+qcAkiZJ+jrhbu89wh0ohC/02ZJ6pf2rkPRiWjurDzbdCxwX78Z2Bf6W1tb01drqLunQ+PpnwBDC77imcUnLM3vLh25P09L51H225rxvSUWPGsya1nI92sOna2jr0/j/ZaSVQZrZxu1oN/Uzp95r6fusw/WXdJOkrxLSDcOBH7cjng7DRTd5xpWUllnZgC1Y8uqDlA0ZseKF8sHbrnKstM8guvTfjMUv3Isa61n+/ovUz62iYstdSz+5+dRvm9mIlj4kzxkPXGZmg82sxMwOBA4H7l/bG+Mj8l+AK82sexS8i4A/p53WD7jQzErN7JvANsCjFmpSj4y53TpCyqI5vue3QGXqmlsYrPvmWmJ5k/BH4PfAJEmpp5NXgaVm9t9x0KyTmY00s53j638BKhe/MuHBxiVzbem/H2nxM0r7DKL7jocy7+FfUjtjCmpqQI31LHvnGRa/9FespJNZ5y4PrOV6tId7Cb+bjWIe9qdpbU0GRpjZDmZWDlzejnYnxvceY6H64EJC3j5Fu66/me1sZruaWSnhj0EtK3+XWcVFN2GqxoyeAzxWNmSkmpcvomzwtiteKxsygubliygfsnKQdaMjLqHusw+ZdeMJLHxmHH2PurS5uXrBx03VC7YHpsVR2yNt1dH5fOcXwIvA88BCwuDOyZKmtfH9FxC+aP+JbdwD/DHt9VeALQmCeCVwnKT5hO/HRYQ7twWEQadzASQ9AFxDGAhbAkwjlHOtjXuAA+OW2FYTcBghnzmdlcLcM57yc2DGoqf++Nrnf76koeuI/Vot1+r99bPp/pXDWPD4b5l1wwnM/u1/sfyDl6nYYpdm4FE11J21luvRHv4PeJ1wFz4VeCMeQ9IHhN/dE8CH8bPahKR5wDeBqwkDXlsCL6S93t7r34PwlLKQ8IQ0H7i2rfF0JD4NOAeIM9Kepv0z0iDkwfaZcfVhkwn5wnFpr10O3CBpyfrG6OQGHdFXqsaMfr1Dg3LaRSHdDeUt0QHqYsKXos1IzSkjk9fjyPudkoyQL5xJEN3FZvZEev2jk7+sa19hpfeCC27C+J1uDtFW5yhJzWqsK1n4rz9Q/dZjPVu6kzWzAYTHs1PjoSbCI9gT+TSjyPkybe4rzc1SU73VzXrnV1+M/98fZC1Ap0VcdHOMaEhSSahTbMkj9dHqaU/+ef4j108APgE2aU1E4wyhM4DfpB3+EfCbtJIgJ89oS1+Rmh/9fNxFR9V//lEJoQwukdpUZyUuujlKnB+/ptUAxqVWAzCzEwmDIGMk/WRtbcaymz0Js59SgzQPABdK+qTDfwgnK6T6SvXUJ6/t0m/ox136b/4iaX0lVhWk6sxKV6tndrKMi26eY2YPEKY77iWpzaPDZrYJ8Kv4Xgiifhjwoqce8hMLbmxjJX2p/tTMdgVeBh6RdHjWg3NW4ANp+U+qNvE5M+vX1jdJminpaKAroUi8F6Gkp9nMzjJr3T7QyS8kvULwcDgsPiE5CeGim+fER8XUtMsvotdAe96/XNJYQl/4OqFg/Dagzsz+1B4hd3KeKwjGMfeknL6c7OOiWwDEwZGUHeGd69iGJD0hqROwOfAUcDpByKeb2Q6tNuDkPDFttH3c/dDMKlo738kMLroFQjT0uBE4ycyOWtv5a2nrP5L2J8zi+TnB9+DNONvtpPbeTTu5QzTxSfk++grBCeCiW1hcRJiu+sBqFnfrRLT+uxzoBKT8BO8GGi2YVfdu8c1OzhKn555G8Db4UdLxFBsuugVEfHzcLO5Wxfrcjmi3WdLf42y3bQjz678PLDCzt9LsDp08QdKdwOPAtZ46yi4uugWGpMWsNIJ+MgPtvxft8XoTZkRtD7wbUw9HFJjRTqGTMmp/08x6tnqm02H4F6QAkTQZ+AHwNTM7N0OfsUjSjwgrSn8rHn4IaDKzy82sJdNrJ0eIKzKkKl9mr+Yp7GQIF90CRdKvgH8Dt2ZyAT5JTZL+ElMPowglST8j+MM+bWbDMvXZzvoTK18OItRr35RwOEWBi25hs2fcTjWz7pn+MElTJW1BWCLndoL/7PSYetjf76RyE0n/JNRmf8/M9k86nkLHRbeAkVTLyoG197MlepLmSzoL6AJ8Jx7+F2G228VxFQEntzgXqAf+5RNiMouLboEjaTphqvAA4vLWWfzsBkl/iqmHnQmmK2OBGjN7yMzas4Chk0Fi5UtKbNs9s9FpOy66RYCk+4G/EtaU2nNt52cohtcl9SOsc/UXQt3vbDNbama7JRGTsyqx8uUrcfeBJGMpZFx0i4eUyclzZtbyGt4ZRtIXkr4FlBHWLusGvBTzvmfFhQOdhIiLZ14CHG5mJyQdTyHiolskxMUP+8fdOUk/Pkqql3QzoQ/uRcgn3gbUm9m46AHrJICkawnLzd/r1Scdj4tuESFpDqGiAOCuJGNJEY12npdUBgwG/kFYXmiumX1iZqOSjbBo2TFup3fUzEYn4KJbZEh6lrCE+Ynra4zT0UiaLekQwkq3lwKDgMkx9XBi0nfnxcRqlS//TjKWQsNFtzi5FPiCDjLG6Wgk1Ui6htA/vxEP30Mw2rnJzHokF13xECtfTiAY41ycdDyFgotuERLLg7aMux1mjNPRxNTD47HkbFPgJcLg22Ize9vMtmy9BWd9kTQeeAQY66mejsFFt0iRtBRITQ9+KslY2oKkKklfI1Q7XAlsC3yQZrTjs90yRyoNNTkbMxsLHRfdIkbS28D5wO5mdl7S8bQFScskXUbou0fGww8RZrv9wsw2SC66wiRWvqTKDOf4H7j1w0W3yJF0K2GV2FvMbETS8bSVmHp4OKYetgLeAf4XWGZmz5vZkGQjLCwkzQP2A8oJK5Q464iLrgOwb9xOy8fHR0kfSBoB9ARuBvYAZsbUw35+Z9YxSHoa+BVwoZntm2w0+YuLroOkOiBVxfBhvoqUpCWSLiAsL3RSPPwkIfXwo1wdMMwzfggsB57yCSzrhouuA4CkmYQBk/5k2Rino4nLC90bUw/bAZ8C1wK1ZjbRzPq33oLTEqsZ48z1lULaj18wZwWSHiLUwyZmjNPRSJomaRCwIWF5+kOBz82sxsx2STa6/CSuKJyqfHkoyVjyERddZ3VOidvnCunxUdICSacBpcBZhAGhV2Le92wz65xshPlFrHz5PnCYmX1rbec7K3HRdVZBUjNh5QcIj48FNfVWUqOk29M8fpcBvwUazOxuM+uTbIT5g6SbgCnAfWa2SdLx5Asuus6XkDSfUAEA8OckY8kk0eO3GyFH+TBh8G2+mX2RyXXlCoyd4naGmXVJNJI8wUXXWSOSXiQMqJ2Qa8Y4HY2kuZKOJCwvdBFBhKfG1MNJPljUMnFF4VRN9FtJxpIveGdyWuN/CSP/DxTD42NcXuiGmHrYOx6+m7Cs/C35WMOcDSR9AhwDbGNmFyUdT67jouu0SCwPGh53i+rxUdJzUXwHAc8A5wFLzOwDM9s82ehyD0kPAH8DrjOz7ZKOJ5dx0XVaJZYHbRN3n08yliSQ9KmkfYEK4GcEd7aP3GhnjRwft1PMrGuikeQwLrrOWpH0HnAmsHO+GON0NJJqJf2C8J05OB5OGe1cZWYVyUWXG8TKl95xd6H/QVozLrpOm5D0e8Kd7i1mtm3S8SRFNNqZFFMPw4DJQCWw3MxeNbPBiQaYMJIWESpfSnFjnDXiouu0h/3j9m0z65ZoJDmApBmSdgC6AmMJdb+zUkY7yUaXHLHy5RrcGGeNuOg6bSaWB6Xu5Kb742NA0nJJPyZ8n46Nh5+M4vvfxTQAmULSpcBigjGOTzhJw0XXaReSZgOHEWatXZ1wODlFTD1MiKmH4UAV4RrVmdkkM9uo1QYKjwFxO99rnVfiF8JpN5ImAuOASwrFGKejkfShpE2BHsDvCQtszjGzpmIx2pFUA2wddx9OMpZcwkXXWVe+E7fPmdmGiUaSw0haKulMgsfvaYTvXMpo57xC87ZYHUnvA+cCo90YJ+Ci66wTsTwolaub54+PrRM9fu+MqYdRwELgFsKy8uPNrFeyEWYOSb8FXicY4xT9Mkr+RXHWGUkLgdSj8r1JxpJPSJoqqQ+hpvWvhEkFC81sQT6tU9dOdo/bmWZWmmgkCeOi66wXkl4DLgeON7OjEw4nr5C0SNLxQGfCqsy9CevUycxOKaSnB0mNrBxYm5ZkLElTML9UJzkk/RyYCUzwx8f2I6lJ0q0x9bAb0ExY5aLJzG4rlJpoSZ8DRwDDi9kYx0XX6ShSo9Qzi7EutaOQ9IqkTgR7yX8SVrlYamYfF4LRjqS/E1JR1xWrZ7GLrtMhxPKgLePuC0nGUghEj9+DCB6/lcBmrDTaOTLPJ6acHLdTi9GzwkXX6TAkfQScDuxkZt9LOJyCIHr8Xh1TD6mpxQ8SjHauMbPyBMNbJ6JlaI+4uyzJWJLARdfpUCSNA54Efl3MxjiZQNLTaR6/rwKXADVm9oaZDUo2uvYhaSmwK2BmdlPS8WQTF10nExwUt2+7r2rHEz1+dyWsaHwVsCPwSXw5b1b4kPQqcAVwQTEZ41i403ecjsXMNgY+I0wC2FDe0TJGzO8eCjySdvgy4JfRpCinMbN5wIZAn1j7XdD4na6TEWJ50EGE2tNrEg6noIlGOxPj7u+AD4D/A+rN7Ekz65tcdG0ilRpZkOcDhG3CRdfJGJL+CdwO/NiNcbLGEklbETx+byEMvs2NVQ+7JhvampFUB2wRdx9NMpZs4KLrZJqz49aNcbJI9Pj9HuE7fkI8/HIU3+/lmtGOpI+B/wIOLnRjHBddJ6PEXG7PuOvGOFkmph7Gx6qHbYA5wK8JRjsTzKxn6y1kD0l/AF4mGOMU7LJH/gVwMo6kJcBX4+74JGMpZiS9J6k/0B34M3A0sMjMluTQ7LA94naWmXVONJIM4aLrZAVJbwA/AY4zs2OSjqeYkVQt6RSCx++ZBBGeGlMPpyU5mBUtQ/vF3XeTiiOTuOg6WUPSGOBj4G9ujJM80eP39zH1sCNQB9xBmO32RzPbIKG45gKHAFuY2cVJxJBJXHSdbJN6jC16X9VcQtJbksoJJX4TgTOAZWY2KwmjHUn/IDitjc2h1EeH4KLrZBVJtQTzFgiDJk4OET1+DyN4/F5MWP05ZbRzdJZTD6fH7dR89JhoCRddJ+tImg58G/iKmV2QdDzOl4kev9fH1ENqcGsCIfVwvZmVZSEGASkv4ZpMf162cNF1EkHS3QS/2JvcGCe3kfRiFN9+BNvOHwK1ZjY100Y7kpYRK1/M7OZMfla2cNF1kuTQuHVjnDwgevzuSfD4/TkhP/9JTD0cmMHPfQP4KXB+IRjjuOg6iSGpiZXlQZ8Ww7z7QiB6/F4e736/EQ8/HsX3Z5mor5V0BfAF8FS+r5zsouskSiwPOoBgau3GOHmGpMej+A4B3iYsUtpgZs9nwGgnVWa4MJ//QLvoOokj6UngVtwYJ2+R9ImkkQSP3xsJg28po53dOugzGoBN4+5jHdFmErjoOrnC94BGgjFOn6SDcdYNSXWSfkjQltTMw5ei+P5gfb03JFURSskOMrMT1yvYhHDRdXKCWB7UO+7Od2Oc/CYa7TwQUw+bE1a2uIGwrPzfzaxH6y202vY44Dngnnw0xvGO7eQMkqqB7ePuX5OMxek4JP1H0hBgA+BPwGHAYjOrNbNR69jsPnE7K9dsKteGi66TU0iaAvwYOMaNcQoLSTWSvkPQndOAMmByTD18pz2DY/HJKOXP/FHHR5s5XHSdnEPSWOA9gjFO3j0+Oq0TUw93xtTDSGAp8AfCbLe72mq0I2kBoWRtmJn9OHMRdyw5sTDlsMqJ/Qh/+UYBvYBFwBTgjqoxo+cmGZuTDGbWheB6BdAlHxZYTBozEzBWUt4IUAoz6w6MI3j8QjBb30PSWu9izex2wqoToyRNhdzWlERFd1jlxJ2BSoKNm4CKtJdrACOUhoypGjP6texH6CSJmW0CzACmSNp+becXO/ksuiniAOq5QPqU3+OACa2tKB1/dgae/bu9S3sP/CE5rCmJpReGVU48B3gaOJJQ21ex2ikV8fiRwNPxfKeIkDSTsL7XKDfGKQ6ix+8tMfWwczx8PyH18Ov4BLQmNuq2wyF07jVgEjmuKYmIbvxhryOMZq4thpJ43nUuvMWHpPHAI7gxTtEh6fUovn2Apwi13HVm9v7qRjtDL33kuA0PPn+5mVWQ45qS1fSCmVWWVHQ/fMj3792e8EMz+7Yz6dx7IP2P//mK82bfdia99vo2DfNmsvzDl2mYN4ueX/sWvfY6eTmwT9WY0a9nLWgnceIjZ1Pc7SppeZLx5CrxEfs6ST9KOpZMEEvDLgGuSjt8wiY/mjDdOnd5iqgp7STrmpLtO91n1VC/s5qbygEaqxegpibqv/gYNYfvVGP1AhoXfkbZkJF07j2Q3vueQcXmqacMygk5YKeIiOtmpcqD5ubzvPsskPzIeIaIHr9j4t3vvvHwffXzZt1J0IZ1IeuaklXR3fiUsVWgzvVzppcA1M16m/KhoyjtM5j6OdNJHevcawCdu29It+0OoGLznSgpW5GaKQEOHVY5caNsxu0kTywP2odwNzM24XCchJH0jCTrsfvxo7r032woq2nZ4lf+xtwJV63yngWP38aCx29bvamsa0pWRbds0NYndRkwvLlu5jQA6mZNo3zwtpQN3pb0Y2VDRrTWjAilIE6RIelZgpnKRWa2V9LxOMnTe59TD05VLqTTdcR+1Ez/N8211QCouYll7z5L1+0OWFMzWdWUbK8rP6p8k+1KamdNo8cuR1E7622673QEnbpvyNI3H1txrMcuR7fWRsWy957/qdlhSQyq5eqjW67G1RrrE7OAZ83sY1bmerNNrl7z75jZ4UkHkSW22vjU6yaXDdxq9SoFOnfrQ9mQkSx773m673AwNf/5N50qelC28RZraqeCUM+bFbItur3Khoxk6RsTaapZSnPNEkr7DKJT197Me+QGmmqW0jBv5trudLFOpd2B7tkJ2clhsr5KbR7QJ/4rCtTY0OJyQd1G7s/SNx+j+w4Hs+ztp+g6cr/WmsqaMXq2RXdR2aCtaa5bTvXkSZQN2gaAkrIN6NStD9WTJ9GpWx9Ke23caiMbbLnrXZJOzUbATm5iZiOAacCDklp9NCoWCmFyRHsZVjnxLsIip19ig+G7s2DSrdTPraLmo9fove8ZrTW1KCMBroFsVy9MKSktqykbsAVLXn1wlTva8sHbfumYmhpRYz2SkJpRYz3NTQ21hOl8ThEj6W3g+8BRZnZc0vE4iTGFFlYKts5d2GDrPZj38Fi6DBhO55791nQa8f1Z05Rsi+44wMqGjKR5+SLKBq+sdS8bMoLm5YsoHzJyxbH5j/2amWOPYfk7z7DkxfHMHHsMy6b+q3zWTSdfZmZ7e+lQcSPpJsLd7l/dGKdoGUeY2rtGuo48gIa5VXRrPbVgsZ2skHXvhWGVEycQpuG1W/AlNTfMrXr7sz9esF3a4fOAP0mq7agYnfzBzEqB+rhb1MY4xZhegNY1pXHxHD69/VwGX3AXJWVrnDvRDDxYNWb0sRkOcwVJTAMeA6yTQJpZbZd+m6b8OA8iXLBbgRozu8PMWk8GOwVHFNnUXe4bScbiJMYaNUVqZslrD7LBNnu1JLjE943JZHCrk3XRjc4+FxOm37WH5cDFVWNGvx79OP8pqRMwnDAv+zTgMzOrMrOdOjZqJ5eRNBs4FhhpZt9POh4nu1SNGf2aGusvaW6sX1E+2Fxfy6zrj6e26i167XVyS29doSlZCTSSmLVjmulNOa2LfzPhr9HFVWNG/7alk+KaSxcDP007fApwn6TG9Y/YyXXMbALBj3VkHGgrKoo1vQBgZlXddjhkaJ+DzqsxszI6QFMyRWLWjvGH3Qd4kHABVh+BrInHHyQYUrR6cSQtkfQzoBNwVDx8F9BgZjeZ2YYtv9spEFJ5uWltXX3AyX/M7FxgaPVbjx1lZnvTQZqSKXJl5YiNWLPL+7j1cXmPVoB3AV+Jh6YCJxbjXVCxYGa9gIWEVScqWjO+LjSK8U7XzIYD7wP3S/pm6nimNKUjyAnRzTRm1hu4DLgo7fAxwEPRwcopIMxsD+B54EZJP0w6nmxRbKK72pJOJfnyB7YoFqaUtFDSxUApcGI8PAFoMrMrYj7YKRAkvQBcC/wgPm46hcnCuO2eL4ILRSK6KSQ1Srov+nHuSFi6+TJgsZk9a2Y+l79AkHQJUA08Y2ZF40VQLJjZNQSbz90lVScdT3soKtFNR9JbkrYE+gG3A3sBH5mZzOwbPtutIEjN+5wfV59wCgAz252wgsTVkl5OOp72UvQdUdJcSWcBXYDvxsOTCAvhXeKj4PmLpBpgm7j7YJKxOB1DXKr9RWCZpLxcRaboRTeFpAZJf4yph12BOcA1wDIze9jn9ucnkt4DzgcON7Nvru18J3eJT59L4m7epoxcdNeApFcl9QcGAuOBw4FZZlZtZnt66iG/kHQr8CbwF9DKmsQAABEMSURBVP/jmdfcH7fDJdW3emYO46LbCpI+k3QCYdbcBUBX4DlC6uHcWLLi5Ae7xO2saJLj5BFmdjShzPNcSR8mHc/64KLbBiTVSbqZcL32Ibha3QrUmdldZtaiUaeTG8Sp4APjrvsx5xFmNoBQ4vmmpERmkXUkLrrtIBrtPCupDBgK/IPgWv+Fmc02sx2TjdBpDUmfEaaIb21mRTNpIp8xs07Ap3G3IIysXHTXEUkzJR1CSDn8hHAX9UYsOft27CxOjiHpIeAvwPVmNnJt5zuJ81bcblwos0dddNcTScsljWGlxy8Ev4dGM7s5TkF2cosT4naqlwTmLmZ2ATASOELSF0nH01G46HYQaR6/Rlil9iVCqdICM3vHzLZONkInRZwy2jPuLvFqlNzDzLYBbgLGS/p70vF0JC66GUDSfyR9jbBM/FWEAv13Y+rhaJ8dlTySlgC7E6xAb0w4HCcNMysH3om7J7Z2bj7iX/4MIqla0v+wqsdvymjnSjPrllx0TpxCOga40Mz2TTgcZyWpCRA98snIpq246GYBSc2SHoqph60Jf8V/Aiw1sxfNbFiS8RUzkn5C8Fp9yvPvyWNm1xPcAHeTtDTpeDKBi26WkfS+pBEEY+VfEx5xp8fUw4GeX0yE1IKmCzz1kxxmthfwQ+AqSa8kHU+m8A6WEJIWS7oQ6AykVs57nDDb7b9jXsvJApLqCAucAhTUoE2+EFf8eBZYHFNyBYuLbsJIapJ0T0w9bA/MBq4mLCv/qJkNbL0FpyOIU0vPBg41s28lHU8xEZ/uUobk/ZOMJRu46OYQkqZIGgz0Be4EDgFmm1ld9BB1Moik3wGvAveZ2ZCk4ykiHo7b4fGpo6Bx0c1BJM2XdBphQCHl9ftizPue54YtGeVrcTvTr3PmiXabhwHn5LuRTVtx0c1h4vJCt6d5/FYDtwD1ZnavmfVNNsLCQ1ITKwfW3mntXGf9iDabfwFel3Rb0vFkCxfdPCF6/HYnCMJDhKmsc81srpltn2x0hUWccno4sIWZXZx0PIWImXUGZsXd3ZKMJdu46OYZkr6QdBRQRlhSvi/wVkw9nOJGOx2DpEeAe4CxZrZd0vEUIG/H7cbx6aJocNHNUyTVS7ohph72iYfvJBjt/MaXle8Qvh23U9wYp+Mws4sIJXqHF5KRTVtx0S0AosevAYOBp4BzCMvKf2hmw1t/t9MScQpq6o/XsiRjKRTiU8N1wH3xaaLocNEtICTNlrQ/UAH8FNgCeD/NaMdnu7WTOBV1FwAzuznhcPKa+LSQWrXjpCRjSRIX3QJEUq2kKwi/34Pj4QmE2W5j/FG5fUh6DfgFcL6Z7Zd0PPlI/INfHXd7FqKRTVtx0S1gosfvpJh62JSwIu6lhGXlXzOzTZKNMH+Q9DNgHvCkG+OsE78GjGBks2RtJxcyLrpFgqQqSV8BugHXEtabmuFGO+1iUNy6MU47iE8H5wNXFrKRTVvxjlNkSFom6RLC7/7YeDhltFNpZmXJRZfbSKon5MkBHksylnzBzDYEngQWSLos6XhyARfdIiWmHibE1MNWwHTCKhe1Zva4mRW88ci6IOlj4LvAN8ys4FY16Eji08C8uDuotXOLCRddB0kfSNqMsG7Y7cCBwOdm1mxmRTVbqC1I+iPwInCPG+O0yqNxO1xSbaKR5BAuus4KJC2RdBZheaFTCQMfL8W87/fi1E0nsFfczvTr8mXM7CTC6thnF4uRTVtx0XW+RFxe6K40j9/5hNHnBjP7q5n1STbC5JHUDPSLu+8nGUuuEZefuht4JdplOmm46DqtEj1++wJ9gPHAccB8M1toZqOSjS5ZJM0leB5vZmaXJB1PLmBmXQjjAwB7JhlLruKi67QJSQslnUDw+D2fsMbb5Jh6OLVYS6gk/QMYB1xT7H+EIu/G7QBJjYlGkqMU5RfFWXeix++taR6/jQTRaTKz282se7IRJsIZcTvZzCoSjSRBzOxSYDPgCEmfJx1PruKi66wz0eO3lJDbnAT8F7DEzKab2ZbJRpc94pTWbnF3eZKxJIWZ7QCMAe6V5It7toKLrrPeSJor6WDCskKXAsOAD2Lq4ZhimO0maRnwVQAz+03C4WSV+HTzZtw9ubVzHRddpwOR1CDpmph6SBnD/I0w2+3aQn/0lvQGwd3tHDPbP+l4skH8g5ryUuhdzEY2bcVF18kIkp6O4jsIeAX4EbDczN4q5AkF0d3tM+BfRVJal1rbbDdJixKNJE9w0XUyiqRPJe1G8Pi9klD3OzOmHr6ebHQZY2jczi/k1IqZfQM4E/g/N7JpOy66TlaIHr+XEfrcYfHwP6P4XhbrOwsCSQ0EK00IZkIFR/TmmATMlfS/SceTT7joOlklGu1MjKmHzYAPgCuAOjN7xsz6td5CfiCpCjgNOMDMCmpwKS5+mioJG9rauc6XcdF1EkPSdElbAV2Bm4G9gS/i3W/eG+1IuhN4FvhzgRnGp+7et5ZUk2gkeYiLrpM4kpZLuoDQH0+Ih1NGOxfm+bLyqSqOGYVgjGNmpxN+prMluefEOuCi6+QMMfUwPqYetiE8wv6KsKz8A2bWK9kI2080xukbd/+TZCzri5ltAfwJeNmNbNYdF10nJ5H0nqQBhCXQ/wwcBSw0s6X55nEgaT7wdWBInCqbd8Qa65RF495JxpLvuOg6OY2kpZJOIXj8nkmYbpsy2jk9X4x2JD0B/B4YY2bbJx3POpBKJQyK1RnOOpIXHdZxosfv72Pq4StALeFRt8nM7jCzbq23kDySzoz/fSufZueZ2eXAEIKRzacJh5P3uOg6eYekNyVVEDx+HyGUZi01s0/ywGhng7jNC2McM9sF+BlwjxvZdAwuuk7eEj1+Dwc6AxcRphynjHaOzcXZYLHEagcAM8vpwSgz602Ywg1wSpKxFBIuuk7eI6lJ0g0x9bBHPHw/wWjnBjMrTzC8LyFpMlAJnJmrU6FjrnxB3N0wVmE4HYCLrlNQSHoxim9/4HngB0CNmU0zs8HJRrcSSVcDswhToXPRGOeOuP2apAWtnei0DxddpyCRNEfSXgSP38uBEcCsmHr4RqLBrWSzuM0pYxwzO4yQTrhK0ktJx1NouOg6BU30+P15vPtNie2kKL6Xm1lpgrE1stK74F9JxZGOmQ0C/k6YmHJZwuEUJC66TtEg6fEovkOAqYRR+XozeyEpox1JM4FvA/uZWaKDVdHp7ZO4u4UbkmcGF12n6JD0iaRRQDlwA/A1VhrtfC2BeO4m3OneaWZJunY9HbfbxuWHnAzgousULZLqJF1E+B4cEw+/EMX3oiwb7aSqGKqSMMYxs3OA3YFzJb27tvOddcdF1yl6otHOAzH1sAWhquA6gtHOxGwY7cRH+VQVw8xMf146ZrYt8BvgJUm/zeZnFyMuuo6ThqSPJW1C8Pj9I3AowWinNtOeCZIWAvsDA8zsfzL5WSni9Om34+6+2fjMYsdF13HWQPT4/S7hO3IqUEbwTJCZfTdTJV6SniLcdf6fme2Qic9IEX+Gj+LuJpLqM/l5TsBF13FaIaYe7oqph5GE5cZ/T5jt9mcz65qBzzwv/vfNDBvjXEWYRHKUpFkZ/BwnDRddx2kjkt6W1JPg8fs34GSg2sy+yIDRTkpsM2KMY2Z7ApcSjGweysRnOGvGRddx2kn0+D2O4PH7PaAfK412ju+I1IOkWmAUgJn9YX3bS8fMNgKei7undmTbztpx0XWcdSR6/N4SUw+7xMPjCamHm82sbD3bnwr8GPhOR01djmVwc+LuRpKaOqJdp+246DpOByDptSi+GwJPAucDtWb2gZkNWY92xwLTCVOXN+yAUO+J2z0lzeuA9px24qLrOB2IpAWSDiB4/P4E2BKYGVMPB69js6l88bz1SV2Y2bHA8cAYSS+sazvO+uGi6zgZIHr8jol3v/vGw49F8b2yPbPOYgogdbf8zLrEY2abEjyGPwOyUgPsrBkXXcfJMJKeieI7AHiDcAfcYGavxEGttrTxCXACsJeZndbKqdutfiCauKeWf9/KjWySxfz65ybDKif2I6z9NQroBSwCpgB3VI0ZPTfJ2Jz1I9pJXgH8d9rhvSQ934b3PgYcDGwqqQpW9pVl7z57RdngkW907t7nY9L6ipm9BuwEbCdpWgf/OE47cdHNMYZVTtyZsJTLIYBYWa8JUAMY8BgwpmrM6NeyH6HTkZjZaMLimikuAa5raXmcmNNtBhh66SO7sJa+Uj93xgfzJ944qv7zD8+T9JtM/AxO+3DRzSGGVU48h2C0Uk7rqZ9mwhLkF1eNGe0GJQVAtHR8kpWrSfwTOF7S4jWcO6jnHic90muvk4azlr6i5mbU3Nhc0rnL+d5XcgMX3RwhTXA3WNu5aSzHhbegiLW9NwLnxEPNwFclvZU6x/tKfuOimwMMq5y4c/W0J59bMOnWLxXTq6GWnnueTM/djmP+pFtY/v4LWGkZPXc9lh67HA3hy7RP1ZjRr2c9cCdjxDTC8cB9aYev3+RHE8Zb5y5P0YLgqrGhpX4C3ldygqybJTtrpLLbyP1Lu43cf5WDSydPYtGzd9Fth4NY9PzdNC78lEHn/Ynm6oV8fm8lpX03oWKzr5YT8nrHJhK5kxFihcF4YLyZDQdeBS5qrlu+Y6fOXVpcUr6VfgIhFeF9JWG8ZCxh4sjzIaz2u6j//GMWPnE7Gx1xCZ279WHZtCfpuccJdCrvRmnfIXTf/iCqpz5BfN+hwyontqn0yMk/JH0gqdfgC+/ZuFPXXrsvfvn+krkPXLXKOQsev40Fj9/WWj8B7ys5gYtu8pxGGHleQXNtNXMfHEPPPU6gfOgommqraapeQJd+m644p7TfpjTMW7HAgGI7TgHTaYMepwLqus3e1Hz8b5rrggGZmptY/t7zdN1277X1E/C+kjguuskzirRSH0nMe+R6SvsOpceu4SlQ9TUAlJSttG4tKetKc11NareieuqT18bZTv6vQP/VffrBiUBF55796LLx5iz/4CUAamdMwUrL6NS979r6CbGvjergPuy0A8/pJs8q628tefl+GubNZMDpN2Jxmr11CZrcXLecTp27rPh/SdnKssxOXXvVsNI9anUyNVqaiXZbatOy/Hm5xuaduvZa4QjWddt9WPbus3Tb7gCWvfMMXbfdp039JJLxNd+clnHRTZ5Fqf/UzpjC4pf+wsYnXU1JebcVJ3Qq70anbn2onzOdik13BKBhznRK+26y4pyKzb5yvyT3Ri1ghlVOvIsws4wNtt6ThU/+gcYl81j+wUsMOGVsm/pJZNHqB5zs4emF5JkC1DRWL2Duw7+k9wFn0mXjzb90UteR+7P4xfE01VbTMH8W1ZMn0W27A1Mv18R2nMJmCuF3TacNelK2yXbMf/RGOvfsT2nf4Iezln4C3lcSx+90k2cc8IvqtybRvGwRC5/4HQuf+N0qJ3QdsS99Djyb+ZNuYfatZ2Cdy+i527GpMiAIj97jshy3k33GAb9I7XTddh/mP3I9vfY7Y8UJvfY8ubV+At5XEscnR+QAwyonTgCOZN2ePJqBB6vGjPbayyLA+0r+4+mF3GAMwUthXaiN73eKA+8reY6Lbg4Q3cIupv0rv6bm0/u0ziLB+0r+4+mFHMJdxpy24n0lf3HRzTGGVU7ciTA//lBa9tN9lOCn63ctRYz3lfzERTdHifPj17RyxDhfOcJJx/tKfuGi6ziOk0V8IM1xHCeLuOg6juNkERddx3GcLOKi6ziOk0VcdB3HcbKIi67jOE4WcdF1HMfJIi66juM4WcRF13EcJ4u46DqO42QRF13HcZws4qLrOI6TRVx0HcdxsoiLruM4ThZx0XUcx8ki/w87wIJxjTZPAgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "model.view_model()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "In the above figure, we have a causal graph that shows the relationships between the treatment, outcome, confounders and the instrument variable.\n", "- The Confounders $W_0$ and $W_1$ affect both the treatment and the outcome\n", "- The instrument variable $Z_0$ is able to effect the outcome $y$ through the treatment $x$" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Identify the Estimand" ] }, { "cell_type": "code", "execution_count": 20, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "INFO:dowhy.causal_identifier:Common causes of treatment and outcome:['W1', 'W0', 'Unobserved Confounders']\n", "WARNING:dowhy.causal_identifier:If this is observed data (not from a randomized experiment), there might always be missing confounders. Causal effect cannot be identified perfectly.\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "WARN: Do you want to continue by ignoring any unobserved confounders? (use proceed_when_unidentifiable=True to disable this prompt) [y/n] y\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "INFO:dowhy.causal_identifier:Instrumental variables for treatment and outcome:['Z0']\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "Estimand type: nonparametric-ate\n", "### Estimand : 1\n", "Estimand name: backdoor\n", "Estimand expression:\n", " d \n", "─────(Expectation(y|W1,W0))\n", "d[v₀] \n", "Estimand assumption 1, Unconfoundedness: If U→{v0} and U→y then P(y|v0,W1,W0,U) = P(y|v0,W1,W0)\n", "### Estimand : 2\n", "Estimand name: iv\n", "Estimand expression:\n", "Expectation(Derivative(y, [Z0])*Derivative([v0], [Z0])**(-1))\n", "Estimand assumption 1, As-if-random: If U→→y then ¬(U →→{Z0})\n", "Estimand assumption 2, Exclusion: If we remove {Z0}→{v0}, then ¬({Z0}→y)\n", "\n" ] } ], "source": [ "identified_estimand = model.identify_effect()\n", "print(identified_estimand)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Estimating the Effect" ] }, { "cell_type": "code", "execution_count": 22, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "INFO:dowhy.causal_estimator:INFO: Using Instrumental Variable Estimator\n", "INFO:dowhy.causal_estimator:Realized estimand: Wald Estimator\n", "Realized estimand type: nonparametric-ate\n", "Estimand expression:\n", " -1\n", "Expectation(Derivative(y, Z0))⋅Expectation(Derivative(v0, Z0)) \n", "Estimand assumption 1, As-if-random: If U→→y then ¬(U →→{Z0})\n", "Estimand assumption 2, Exclusion: If we remove {Z0}→{v0}, then ¬({Z0}→y)\n", "Estimand assumption 3, treatment_effect_homogeneity: Each unit's treatment ['v0'] is affected in the same way by common causes of ['v0'] and y\n", "Estimand assumption 4, outcome_effect_homogeneity: Each unit's outcome y is affected in the same way by common causes of ['v0'] and y\n", "\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "*** Causal Estimate ***\n", "\n", "## Target estimand\n", "Estimand type: nonparametric-ate\n", "### Estimand : 1\n", "Estimand name: backdoor\n", "Estimand expression:\n", " d \n", "─────(Expectation(y|W1,W0))\n", "d[v₀] \n", "Estimand assumption 1, Unconfoundedness: If U→{v0} and U→y then P(y|v0,W1,W0,U) = P(y|v0,W1,W0)\n", "### Estimand : 2\n", "Estimand name: iv\n", "Estimand expression:\n", "Expectation(Derivative(y, [Z0])*Derivative([v0], [Z0])**(-1))\n", "Estimand assumption 1, As-if-random: If U→→y then ¬(U →→{Z0})\n", "Estimand assumption 2, Exclusion: If we remove {Z0}→{v0}, then ¬({Z0}→y)\n", "\n", "## Realized estimand\n", "Realized estimand: Wald Estimator\n", "Realized estimand type: nonparametric-ate\n", "Estimand expression:\n", " -1\n", "Expectation(Derivative(y, Z0))⋅Expectation(Derivative(v0, Z0)) \n", "Estimand assumption 1, As-if-random: If U→→y then ¬(U →→{Z0})\n", "Estimand assumption 2, Exclusion: If we remove {Z0}→{v0}, then ¬({Z0}→y)\n", "Estimand assumption 3, treatment_effect_homogeneity: Each unit's treatment ['v0'] is affected in the same way by common causes of ['v0'] and y\n", "Estimand assumption 4, outcome_effect_homogeneity: Each unit's outcome y is affected in the same way by common causes of ['v0'] and y\n", "\n", "## Estimate\n", "Value: 10.006229589077737\n", "\n" ] } ], "source": [ "causal_estimate = model.estimate_effect( identified_estimand,\n", " method_name=\"iv.instrumental_variable\",\n", " method_params={'iv_instrument_name':'Z0'}\n", " )\n", "print(causal_estimate)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Refuting the Estimate" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Using a Randomly Generated Value" ] }, { "cell_type": "code", "execution_count": 24, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "INFO:dowhy.causal_refuters.dummy_outcome_refuter:Refutation over 100 simulated datasets of Random Data treatment\n", "INFO:dowhy.causal_estimator:INFO: Using Instrumental Variable Estimator\n", "INFO:dowhy.causal_estimator:Realized estimand: Wald Estimator\n", "Realized estimand type: nonparametric-ate\n", "Estimand expression:\n", " -1\n", "Expectation(Derivative(dummy_outcome, Z0))⋅Expectation(Derivative(v0, Z0)) \n", "Estimand assumption 1, As-if-random: If U→→y then ¬(U →→{Z0})\n", "Estimand assumption 2, Exclusion: If we remove {Z0}→{v0}, then ¬({Z0}→y)\n", "Estimand assumption 3, treatment_effect_homogeneity: Each unit's treatment ['v0'] is affected in the same way by common causes of ['v0'] and dummy_outcome\n", "Estimand assumption 4, outcome_effect_homogeneity: Each unit's outcome dummy_outcome is affected in the same way by common causes of ['v0'] and dummy_outcome\n", "\n", "INFO:dowhy.causal_refuters.dummy_outcome_refuter:Making use of Bootstrap as we have more than 100 examples.\n", " Note: The greater the number of examples, the more accurate are the confidence estimates\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "Refute: Use a Dummy Outcome\n", "Estimated effect:10.006229589077737\n", "New effect:2.0838738761455298e-06\n", "p value:0.99\n", "\n" ] } ], "source": [ "ref = model.refute_estimate(identified_estimand,\n", " causal_estimate,\n", " method_name=\"dummy_outcome_refuter\"\n", " )\n", "print(ref)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The result shows that the treatment does not lead to the outcome. The esitimated outcome $2x10^{-6}$ is a value that tends to zero, which matches our expectation. This shows that if we replace the outcome by randomly generated data, the estimator correctly predicts that the influence if treatment is Zero." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Using a Function that Generates the Outcome from the Confounders" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Let us define a simple function that generates the outcome as a linear function of the confounders." ] }, { "cell_type": "code", "execution_count": 26, "metadata": {}, "outputs": [], "source": [ "coefficients = np.array([1,2])\n", "bias = 3\n", "def linear_gen(df):\n", " y_new = np.dot(df[['W0','W1']].values,coefficients) + 3\n", " return y_new" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The basic expression is of the form\n", "$y_{new} = \\beta_0W_0 + \\beta_1W_1 + \\gamma_0$\n", "\n", "where,\n", "$\\beta_0=1$, $\\beta_1=2$ and $\\gamma_0=3$" ] }, { "cell_type": "code", "execution_count": 28, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "INFO:dowhy.causal_refuters.dummy_outcome_refuter:Refutation over 100 simulated datasets of Random Data treatment\n", "INFO:dowhy.causal_estimator:INFO: Using Instrumental Variable Estimator\n", "INFO:dowhy.causal_estimator:Realized estimand: Wald Estimator\n", "Realized estimand type: nonparametric-ate\n", "Estimand expression:\n", " -1\n", "Expectation(Derivative(dummy_outcome, Z0))⋅Expectation(Derivative(v0, Z0)) \n", "Estimand assumption 1, As-if-random: If U→→y then ¬(U →→{Z0})\n", "Estimand assumption 2, Exclusion: If we remove {Z0}→{v0}, then ¬({Z0}→y)\n", "Estimand assumption 3, treatment_effect_homogeneity: Each unit's treatment ['v0'] is affected in the same way by common causes of ['v0'] and dummy_outcome\n", "Estimand assumption 4, outcome_effect_homogeneity: Each unit's outcome dummy_outcome is affected in the same way by common causes of ['v0'] and dummy_outcome\n", "\n", "INFO:dowhy.causal_refuters.dummy_outcome_refuter:Making use of Bootstrap as we have more than 100 examples.\n", " Note: The greater the number of examples, the more accurate are the confidence estimates\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "Refute: Use a Dummy Outcome\n", "Estimated effect:10.006229589077737\n", "New effect:-5.262651999486695e-06\n", "p value:1.0\n", "\n" ] } ], "source": [ "ref = model.refute_estimate(identified_estimand,\n", " causal_estimate,\n", " method_name=\"dummy_outcome_refuter\",\n", " outcome_function=linear_gen\n", " )\n", "\n", "print(ref)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Like the previous experiment, we observe that the estimator shows that the effect of the treatment is Zero. The refuter confirms this as the value obtained through the refutation is $5x10^{-6}$ and has a p value of 1.0 across 100 simulations." ] } ], "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 }