Simulating_a_news_personalization_scenario_using_Contextual_Bandits.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Simulating a news personalization scenario using Contextual Bandits\n",
    "\n",
    "In the Contextual Bandit(CB) introduction tutorial, we learnt about CB and different CB algorithms. In this tutorial we will simulate the scenario of personalizing news content on a site, using CB, to users. The goal is to maximize user engagement quantified by measuring click through rate (CTR).\n",
    "\n",
    "Let's recall that in a CB setting, a data point has four components,\n",
    "\n",
    "- Context\n",
    "- Action\n",
    "- Probability of choosing action\n",
    "- Reward/cost for chosen action\n",
    "\n",
    "In our simulator, we will need to generate a context, get an action/decision for the given context and also simulate generating a reward.\n",
    "\n",
    "In our simulator, our goal is to maximize reward (click through rate/CTR) or minimize loss (-CTR)\n",
    "\n",
    "- We have two website visitors: 'Tom' and 'Anna'\n",
    "- Each of them may visit the website either in the morning or in the afternoon\n",
    "\n",
    "The **context** is therefore (user, time_of_day)\n",
    "\n",
    "We have the option of recommending a variety of articles to Tom and Anna. Therefore, **actions** are the different choices of articles: \"politics\", \"sports\", \"music\", \"food\", \"finance\", \"health\", \"cheese\"\n",
    "\n",
    "The **reward** is whether they click on the article or not: 'click' or 'no click'\n",
    "\n",
    "Let's first start with importing the necessary packages:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "from vowpalwabbit import pyvw\n",
    "import random\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Simulate reward\n",
    "\n",
    "In the real world, we will have to learn Tom and Anna's preferences for articles as we observe their interactions. Since this is a simulation, we will have to define Tom and Anna's preference profile. The reward that we provide to the learner will follow this preference profile. Our hope is to see if the learner can take better and better decisions as we see more samples which in turn means we are maximizing the reward.\n",
    "\n",
    "We will also modify the reward function in a few different ways and see if the CB learner picks up the changes. We will compare the CTR with and without learning.\n",
    "\n",
    "VW optimizes to minimize **cost which is negative of reward**. Therefore, we will always pass negative of reward as cost to VW."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "# VW tries to minimize loss/cost, therefore we will pass cost as -reward\n",
    "USER_LIKED_ARTICLE = -1.0\n",
    "USER_DISLIKED_ARTICLE = 0.0"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The reward function below specifies that Tom likes politics in the morning and music in the afternoon whereas Anna likes sports in the morning and politics in the afternoon. It looks dense but we are just simulating our hypothetical world in the format of the feedback the learner understands: cost. If the learner recommends an article that aligns with the reward function, we give a positive reward. In our simulated world this is a click."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "def get_cost(context,action):\n",
    "    if context['user'] == \"Tom\":\n",
    "        if context['time_of_day'] == \"morning\" and action == 'politics':\n",
    "            return USER_LIKED_ARTICLE\n",
    "        elif context['time_of_day'] == \"afternoon\" and action == 'music':\n",
    "            return USER_LIKED_ARTICLE\n",
    "        else:\n",
    "            return USER_DISLIKED_ARTICLE\n",
    "    elif context['user'] == \"Anna\":\n",
    "        if context['time_of_day'] == \"morning\" and action == 'sports':\n",
    "            return USER_LIKED_ARTICLE\n",
    "        elif context['time_of_day'] == \"afternoon\" and action == 'politics':\n",
    "            return USER_LIKED_ARTICLE\n",
    "        else:\n",
    "            return USER_DISLIKED_ARTICLE"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "\n",
    "## Understanding VW format\n",
    "\n",
    "There are some things we need to do to get our input into a format VW understands. This function handles converting from our context as a dictionary, list of articles and the cost if there is one into the text format VW understands.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "# This function modifies (context, action, cost, probability) to VW friendly format\n",
    "def to_vw_example_format(context, actions, cb_label = None):\n",
    "    if cb_label is not None:\n",
    "        chosen_action, cost, prob = cb_label\n",
    "    example_string = \"\"\n",
    "    example_string += \"shared |User user={} time_of_day={}\\n\".format(context[\"user\"], context[\"time_of_day\"])\n",
    "    for action in actions:\n",
    "        if cb_label is not None and action == chosen_action:\n",
    "            example_string += \"0:{}:{} \".format(cost, prob)\n",
    "        example_string += \"|Action article={} \\n\".format(action)\n",
    "    #Strip the last newline\n",
    "    return example_string[:-1]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "To understand what's going on here let's go through an example. Here, it's the morning and the user is Tom. There are four possible articles. So in the VW format there is one line that starts with shared, this is the shared context, followed by four lines each corresponding to an article."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "shared |User user=Tom time_of_day=morning\n",
      "|Action article=politics \n",
      "|Action article=sports \n",
      "|Action article=music \n",
      "|Action article=food \n"
     ]
    }
   ],
   "source": [
    "context = {\"user\":\"Tom\",\"time_of_day\":\"morning\"}\n",
    "actions = [\"politics\", \"sports\", \"music\", \"food\"]\n",
    "\n",
    "print(to_vw_example_format(context,actions))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Getting a decision\n",
    "\n",
    "When we call VW we get a _pmf_, [probability mass function](https://en.wikipedia.org/wiki/Probability_mass_function), as the output. Since we are incorporating exploration into our strategy, VW will give us a list of probabilities over the set of actions. This means that the probability at a given index in the list corresponds to the likelihood of picking that specific action. In order to arrive at a decision/action, we will have to sample from this list.\n",
    "\n",
    "So, given a list `[0.7, 0.1, 0.1, 0.1]`, we would choose the first item with a 70% chance. `sample_custom_pmf` takes such a list and gives us the index it chose and what the probability of choosing that index was."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "def sample_custom_pmf(pmf):\n",
    "    total = sum(pmf)\n",
    "    scale = 1/total\n",
    "    pmf = [x * scale for x in pmf]\n",
    "    draw = random.random()\n",
    "    sum_prob = 0.0\n",
    "    for index, prob in enumerate(pmf):\n",
    "        sum_prob += prob\n",
    "        if(sum_prob > draw):\n",
    "            return index, prob"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We have all of the information we need to choose an action for a specific user and context. To use VW to achieve this, we will do the following:\n",
    "\n",
    "1. We convert our context and actions into the text format we need\n",
    "2. We pass this example to vw and get the pmf out\n",
    "3. Now, we sample this pmf to get what article we will end up showing\n",
    "4. Finally we return the article chosen, and the probability of choosing it (we are going to need the probability when we learn form this example)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "def get_action(vw, context, actions):\n",
    "    vw_text_example = to_vw_example_format(context,actions)\n",
    "    pmf = vw.predict(vw_text_example)\n",
    "    chosen_action_index, prob = sample_custom_pmf(pmf)\n",
    "    return actions[chosen_action_index], prob"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "\n",
    "\n",
    "## Simulation set up\n",
    "\n",
    "Now that we have done all of the setup work and know how to interface with VW, let's simulate the world of Tom and Anna. The scenario is they go to a website and are shown an article. Remember that the reward function allows us to define the worlds reaction to what VW recommends.\n",
    "\n",
    "\n",
    "We will choose between Tom and Anna uniformly at random and also choose their time of visit uniformly at random. You can think of this as us tossing a coin to choose between Tom and Anna (Anna if heads and Tom if tails) and another coin toss for choosing time of day.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "users = ['Tom', 'Anna']\n",
    "times_of_day = ['morning', 'afternoon']\n",
    "actions = [\"politics\", \"sports\", \"music\", \"food\", \"finance\", \"health\", \"camping\"]\n",
    "\n",
    "def choose_user(users):\n",
    "    return random.choice(users)\n",
    "\n",
    "def choose_time_of_day(times_of_day):\n",
    "    return random.choice(times_of_day)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We will instantiate a CB learner in VW and then simulate Tom and Anna's website visits `num_iterations` number of times. In each visit, we:\n",
    "\n",
    "1. Decide between Tom and Anna\n",
    "2. Decide time of day\n",
    "3. Pass context i.e. (user, time of day) to learner to get action i.e. article recommendation and probability of choosing action\n",
    "4. Receive reward i.e. see if user clicked or not. Remember that cost is just negative reward.\n",
    "5. Format context, action, probability, reward in VW format\n",
    "6. Learn from the example\n",
    "    - VW _reduces_ a CB problem to a cost sensitive multiclass classification problem.\n",
    "\n",
    "This is the same for every one of our simulations, so we define the process in the `run_simulation` function. The cost function must be supplied as this is essentially us simulating how the world works.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "def run_simulation(vw, num_iterations, users, times_of_day, actions, cost_function, do_learn = True):\n",
    "    cost_sum = 0.\n",
    "    ctr = []\n",
    "\n",
    "    for i in range(1, num_iterations+1):\n",
    "        # 1. In each simulation choose a user\n",
    "        user = choose_user(users)\n",
    "        # 2. Choose time of day for a given user\n",
    "        time_of_day = choose_time_of_day(times_of_day)\n",
    "\n",
    "        # 3. Pass context to vw to get an action\n",
    "        context = {'user': user, 'time_of_day': time_of_day}\n",
    "        action, prob = get_action(vw, context, actions)\n",
    "\n",
    "        # 4. Get cost of the action we chose\n",
    "        cost = cost_function(context, action)\n",
    "        cost_sum += cost\n",
    "\n",
    "        if do_learn:\n",
    "            # 5. Inform VW of what happened so we can learn from it\n",
    "            vw_format = vw.parse(to_vw_example_format(context, actions, (action, cost, prob)),pyvw.vw.lContextualBandit)\n",
    "            # 6. Learn\n",
    "            vw.learn(vw_format)\n",
    "\n",
    "        # We negate this so that on the plot instead of minimizing cost, we are maximizing reward\n",
    "        ctr.append(-1*cost_sum/i)\n",
    "\n",
    "    return ctr"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "\n",
    "We want to be able to visualize what is occurring, so we are going to plot the click through rate over each iteration of the simulation. If VW is showing actions the get rewards the ctr will be higher. Below is a little utility function to make showing the plot easier.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "def plot_ctr(num_iterations, ctr):\n",
    "    plt.plot(range(1,num_iterations+1), ctr)\n",
    "    plt.xlabel('num_iterations', fontsize=14)\n",
    "    plt.ylabel('ctr', fontsize=14)\n",
    "    plt.ylim([0,1])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Scenario 1\n",
    "\n",
    "We will use the first reward function `get_cost` and assume that Tom and Anna do not change their preferences over time and see what happens to user engagement as we learn. We will also see what happens when there is no learning. We will use the \"no learning\" case as our baseline to compare to.\n",
    "\n",
    "### With learning\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAEOCAYAAACNY7BQAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO3deZxddX3/8dfnLrNktuwhZCFBghBBkEZAReCnIItK/NX+BKrVWltaq/7cWh9QW7R0QS0WsUUpVX5QqyJ1jRrLTlEUJMiaxIRJCCQhy2SbZNa7fX5/nDOTe27uJDMwc+7MnPfz8biPOds993Nubu77nu8553vM3RERERmQqnUBIiIyvigYREQkQsEgIiIRCgYREYlQMIiISISCQUREImILBjO7xcx2mtkzQ8w3M/uymbWb2VNmdlpctYmIyEFx7jHcClx4mPkXAUvCxxXAV2OoSUREKsQWDO7+ILDnMIssB/7DAw8DU81sbjzViYjIgEytCygzD9hcNr4lnLatckEzu4Jgr4KmpqbfOeGEE2IpUERksnjsscd2ufusavPGUzAMm7vfDNwMsGzZMl+1alWNKxIRmVjM7Pmh5o2ns5K2AgvKxueH00REJEbjKRhWAO8Nz046E+h090OakUREZGzF1pRkZt8GzgVmmtkW4DNAFsDdbwJWAhcD7UAP8P64ahMRkYNiCwZ3v/wI8x34UEzliIjIEMZTU5KIiIwDCgYREYlQMIiISISCQUREIhQMIiISoWAQEZEIBYOIiEQoGEREJELBICIiEQoGERGJUDCIiEiEgkFERCIUDCIiEqFgEBGRCAWDiIhEKBhERCRCwSAiIhEKBhERiVAwiIhIhIJBREQiFAwiIhKhYBARkQgFg4iIRCgYREQkQsEgIiIRCgYREYlQMIiISISCQUREIhQMIiISoWAQEZEIBYOIiEQoGEREJELBICIiEQoGEZEq3B13r3UZuDt9+WKstWRieyURGdf68kV27O9jQ0cXhaJTcmfJnBZyhRLplNGYTdNUn6G5PkM6ZaQMzOyw6yyWnF1d/eQKJfryRXpyRZ7f08PWvb3s682xrzvP7u4cTfVpcoUSuUKJ3d05GrIpGrNpUmakUkZLfYbpTXU01Wcolpzd3Tk6e3N09RfZ250jXyzRkyvSkyuQLzrHzmoiZUY6ZfQXSsyf1khnTx4zaG3I4ji5QomSQ75Yolhy0imjIZtmW2cv2/b10dVfwB2yaWPhjCYMKLmTTafIF0u4Q75UImVGa0OG1sYss5rr6ckVKZac+myKnlyR3lyRQqlEb75E2oL3xIF8MfjCNwMc6rNpmurStDZm6eorsKcnx+6ufrr6g23Kpo3pTXX05IrUZ1J09xe5+u1Lufz0haP+WYg1GMzsQuAGIA18zd0/VzF/IXAbMDVc5kp3XxlnjSITSV++CEDKjBf2dNNcn2VXVz+zW+qZ2VzPwPe2OziwvzfPmm376TjQz8Zd3Wze08PGji729uTZuq+XYmn4v0rrMymy6RRtjVm6cwVKJWdKXRAac9sa2NbZx4udvQz1Q9cMZjTV01yfpuhOJpWiIZtmRlMdvfki+3v7KbnjDnt7cuzpzlEoOWYwfUodzQ0Z2hqzzGiuoy6dYkpdmsa6NLmC8+K+XvLF0mCdT27ex4ymOjBj854e0ikjk0qRSRvZdIp0yiiUnJ37+5k7tYGT57XR2pClUAq+vLd39g2+l4VSUGvKIJtOUSiV2N9boONAP09t6WRqY5Z8qYRhTKlLM6UuTToVhIc7pFKGEQROY12GUviP05cv0tVfYOeBPprqMiyZ3cyZx06npSFLc32Gzt48nT156jIpCiWnqS7NCUe1vMRPzuHFFgxmlgZuBM4HtgCPmtkKd19TtthfA3e4+1fNbCmwElgUV40iY6VUcjp7g//U+/vybOzo5kBfnjmtDTTWpTl2ZjMld7Z19rH6xU62d/ZRKDn5QokNHV30hr+20ykjbca6HQfY15Onq78ABF+ylV/AKQvCIG3Bl16llMGc1gYWz2zi6KmNLD/1aOZNbeQVs5upC38Vt+/soqk+Q2/YlLG3J8+BvjzFUvBFViw5Xf0FGrIpMqkUOw/0kTLjxc4+XrtoGvOnzWNOW0MYIsFex+KZzcyb1khTXfqIexyVcoUS2bSN+HkyMnHuMZwOtLv7RgAzux1YDpQHgwOt4XAb8GKM9Ym8bD25AvmCs7cnx+oX9/Pb7ftZu+0Aq57fw76e/Eta5+yWeprqM7Q2ZCi6058vccqCqcxsqmNGcz0pg/5CiaOnNtKXLzJ1Spa93Xl2dfVjBn35Epm0kUkZUxvrOGFuCzOa6lk8s4nGuvRhX3vZoukvqeaxUpfRYdE4xBkM84DNZeNbgDMqlvkscJeZfQRoAs6rtiIzuwK4AmDhwtFvX5PJyd3ZuKub/nyJuowxp7WB5vpM1V+fhWKJdTsOsHN/P7liiZb6DJv39tDZm6e7P2ha2NXVT0dXP9s7+8gVS/TnS/SGTTvlFk6fwvknzmHRzCZyhRJtjdnBNvDO3jx7unN0HOgnkzaOam3gVUe3MXdqA9lUirpM6ohf3iKjbbwdfL4cuNXdv2hmrwO+YWYnuXupfCF3vxm4GWDZsmW1P21Axg1355Hn9vBQ+y6KJWdDRxdb9vayY38f/fkSB8Kml3KvnNPCxSfPZeu+HnYe6OfFfb2s39F12Ndpbcgws7memS31vO4VM6jPpKjPpJnRXEex5MxpbeDEuS0sndumL3aZcOIMhq3AgrLx+eG0ch8ALgRw91+ZWQMwE9gZS4UyIfx2+35u++XzPLVlHyWHXKHI9KY6GusyPLvjANs6+waXPbqtgXnTGll6dBst9RlOWdDGlPCA3+6uHHt7cty7difX37Oe5voMc1rrmdvWyBuXzOJVR7cyp7UBd9jd3c8rwwN9C6dPYUrdePtNJTJ64vx0PwosMbPFBIFwGfD7Fcu8ALwZuNXMTgQagI4Ya5Qa27qvl4ee3UVrY5b6bIqfr9/Fxl1dHDermWd3drH6xf3s6uonnTLmtNQzvTloZ+/uL7C7K8eSOS188i2v5NQFbcybOmVYv9b/6uIiW/b2cuzMJlIpHdQUiS0Y3L1gZh8G7iQ4FfUWd19tZtcAq9x9BfBJ4N/N7OMEB6L/0MfDFSbysrk7JYe12/YzvamOuW0NdPUX+OWG3fyyfRdPb+0E4MktnVVPmXxgXfD74JQFU/mzc47lnafNZ1pT3ajU1pBNc9zs5lFZl8hkEOv+cHhNwsqKaVeXDa8B3hBnTTL6XtzXyyfueIKOA/1s6Oge9vPSKWP5qUfT2pANLq6a3cwpC6bS1pilqT5o0xeRsaeGUhkVd67ezqd/8DS7unJV589srufEuS0cM2MKe7pz9OaKNGTTXHjSUZxz/Cyy6eA0xKZ6fSRFak3/C+UlyRdLPLxxN7c+tIl7fxs9N2Dp3FYuPvko3vXaBcxuaahRhSLyUikYZNjcnV1dOb50z3q++cgLkXl/ecErefcZC2lrzOqqVJEJTsEgVT32/F4u//eHyRVKQy6z7JhpfOTNSzjn+FkxViYiY03BkHC5Qomrvv803/vNFl7/ihm8a9kCPvadJw77nG/98Rm8/riZMVUoInFTMCSQu7Njfz9/86NnuHvNjsHpv9ywm19u2A3A/GmNfO+Dr2dOq44RiCSNgiEh+vJF/uDrj/Dopr2HzPu7d5zEyfPa+OPbVnH64mlcf+mp1GfUjYNIUikYJjl35xftu/iDr//6kHl/es6xXPHGY5kRXh+w6q+r9lkoIgmjYJjEOnvz/O5XHhq8yOwL73w1bVOynL5o+qhdNSwik4+CYRIolXUhceX3n+KOVVsOWeaHH3oDpy6YGmdZIjJBKRgmsG898gJf+8VGNh6m24krLzqBPzvnFTFWJSITnYJhAtqxv49rV67lh08ceoO7tsYsd338bDIpGzx2ICIyEgqGCeTG+9v5pzvXRaa9ck4Lt/7Ra3GHo6c21qgyEZlMFAwTxA8f33pIKPzyyjcpDERk1CkYJoDr717PDfc+C8B1/+cUXrNwKq+YpfsHiMjYUDCMc/ev2zkYCnd//GyWzGmpcUUiMtkpGGIycHzgjMXT+c6fvu6Iyz+4voP33nLworSff+p/sWD6lLEsUUQEgFStC0iCrft6B48PPPLcHhZd+VNWbdpTddmdB/r43M9+GwmFL116qkJBRGKjPYYxlC+W+ObDz/PZH685ZN7v3fQrAGY01fG9D76enz697ZCDy3PbGvj4+cfzjtfMi6VeERFQMIypt335F6zbcWBwfNPn3sqPntjKR28/2K317u4c5173wCHP/dSFr+TPzz0ujjJFRCIUDGPkztXbI6HwzN9eAMDyU+ex/NR5fPWBDaQMrv3ZbweX+czbl/KmE2ZzzIym2OsVERmgYBgDF1z/4GAo/N83L+HdZyykueIm9x88N+im4oS5rSycPoWF06eQTumWmCJSewqGYbp7zQ66+wtHbO9/ZOPuwVBoqc/wifOPP+zyui2miIw3CoZh+pP/WAXAO14zj/vX7WRWcz0nzWuLLNPZk+fSmx8eHP/iu06JtUYRkdGgYBiGb//6hcHhRVf+dHB40+feGlnulGvuGnKeiMhEoesYjuCpLfu46vtPV5130Q0/xz24F8L2zr7B6e3/cFEstYmIjAUFwxE8uL5jyHlrt+3nO49uBuCWh54D4P1vWEQmrbdVRCYufYMdwXV3rT9k2u8cM21w+M7V2+nLF3ly8z6mTsly9duWxlmeiMio0zGGYXru2osxO3g66TNbO3nbv/yC+9d1cOo1d9GXL3HSvNbIMiIiE5H2GA4jXywNDld+4ZefkdSXD5Y7Zb7uqSwiE5+C4TB27A8OKF910QlV56+55oLB4fpMin/43yfHUpeIyFhSMBzGpl09AJw4t7Xq/Cl1Ge75xNkAfO+Dr4+tLhGRsaRjDIfxnq8/AgS9nA7luNktumZBRCYV7TEMw1GHCQYRkckm1mAwswvNbJ2ZtZvZlUMs8y4zW2Nmq83sW3HWN5SWhmytSxARiU1sTUlmlgZuBM4HtgCPmtkKd19TtswS4CrgDe6+18xmx1VfpWe2dtbqpUVEairOPYbTgXZ33+juOeB2YHnFMn8C3OjuewHcfWeM9UXs68kDcNFJR9WqBBGRmogzGOYBm8vGt4TTyh0PHG9mD5nZw2Z2YbUVmdkVZrbKzFZ1dAzdZcXLcaAvCIaPvGnJmKxfRGS8Gm8HnzPAEuBc4HLg383skKvG3P1md1/m7stmzRqb+xl0dPUDMKulfkzWLyIyXsUZDFuBBWXj88Np5bYAK9w97+7PAesJgiJ2uw70kzKY3lRXi5cXEamZOIPhUWCJmS02szrgMmBFxTI/JNhbwMxmEjQtbYyxxkFfvq+dkqPbbYpI4sQWDO5eAD4M3AmsBe5w99Vmdo2ZXRIudiew28zWAPcDf+nuu+OqUUREYr7y2d1XAisrpl1dNuzAJ8JHzQz0kdRcrwvDRSR5xtvB53HhjH+8FwhuuiMikjQKhsOYOkUHnkUkeRQMFUolHxy+e832GlYiIlIbCoYKd63ZMTj81pPn1rASEZHaUDBU2BVe2AbwnjOPqWElIiK1oWCoMLM5uNL5Jx85S/dvFpFEUjBU2B/2kdTWqK62RSSZFAwV9vcGwdCqYBCRhFIwVLh/XdDTd4subhORhBpWMJhZ1sy2m9mrxrqgWnuoPeiBI6U+kkQkoYYVDO6eB/KAH2lZERGZ2EbSlPQvwFVmpjYWEZFJbCRf8m8EzgG2mtkzQHf5THe/pOqzJpCgDz8RkWQbSTDsAr43VoWMB/2FEgAfP+/4GlciIlI7ww4Gd3//WBYyHnT3FwCY1qRTVUUkuYZ9jMHM7qt2/2UzazWz+0a3rNro7i8CMKVOh1FEJLlGcvD5XKBaP9QNBMcfJryucI+huT5d40pERGrniD+Nzey0stFXm9mesvE0cAGwdbQLi9t5//w/bO8M7tzWpIvbRCTBhvMNuIrg+gUH7qoyvxf4yGgWVQvtO7sGh7v6CjWsRESktoYTDIsBAzYCryU4O2lADtjp7sUxqK1m1KuqiCTZEYPB3Z8HMLPPAa9195vK55vZn5nZPHf/mzGqccwViqXI+PlL59SoEhGR2hvJwef3AI9Xmf4b4L2jU05t7DzQHxlPq58kEUmwkQTDbKCjyvRdwIT+iX3v2h1HXkhEJCFGEgwvUP201LOBLaNTTm3UZQ6+DZ95+9IaViIiUnsjOS/z34DrzawOGLig7c3AtcDnR7uwON2xKsi1r777NC46eW6NqxERqa2RdInxRTObCXyZgxe65YAb3P0LY1FcXB57fi8ADVld2CYiMqIrudz9KjP7e2CgvWWtu3cd7jkiIjKxjPgSX3fvBh4dg1pid8M9z3L28TMHx133IRIRGXkwTCbX37Oe6+9ZPzi+YNqUGlYjIjI+jOSspEnr989YCMCSOS01rkREpPYUDARXPs9ta6h1GSIi44KCAcgXPXItg4hIkiX6GMOAHzw+4XsNFxEZNfqZLCIiEQoGERGJiDUYzOxCM1tnZu1mduVhlnunmbmZLYurtmbdtU1EBIjxGIOZpYEbgfMJOt171MxWuPuaiuVagI8Cj8RV26IZU3j1/KlxvZyIyLgW5x7D6UC7u2909xxwO7C8ynJ/R9ApX19cheWLTiatezCIiEC8wTAP2Fw2viWcNsjMTgMWuPtPD7ciM7vCzFaZ2aqOjmq3iBiZQqlENqXDLSIiMI4OPptZCvhn4JNHWtbdb3b3Ze6+bNasWS/7tQvaYxARGRRnMGwFFpSNzw+nDWgBTgIeMLNNwJnAirE6AO1+sMO83d059nTnxuJlREQmnDiD4VFgiZktDm/2cxmwYmCmu3e6+0x3X+Tui4CHgUvcfdVYFOMVHan+7JntY/EyIiITTmzB4O4F4MPAncBa4A53X21m15jZJXHVISIihxfryfvuvhJYWTHt6iGWPXdMaxnLlYuITGDj5uBz3LyyLUlERIAEB0Ol973umFqXICIyLiQ2GCr3Fwol7UGIiECSg6EiB/LFUm0KEREZZxIbDJXOO3FOrUsQERkXEhsMXtGYdOLc1hpVIiIyviQ2GCqlU+oSQ0QEEhwMlccYMgoGEREgwcFQKaVgEBEBFAyDtMcgIhJIbDBUNiXpGIOISCC5wVBxVpKCQUQkkNhgqKRgEBEJJDYYDmlKMgWDiAgkORgqxrXHICISSGwwVDLtMYiIAAkOBt2PQUSkuuQGQ60LEBEZpxIbDCIiUl1ig0EtSSIi1SU2GMrbktoas7WrQ0RknEluMJSZ0VxX6xJERMaNxAZDeZcYalYSETkoucFQFgYlJYOIyKDEBkO5YknBICIyILHBUB4F2mEQETkoucFQlgZ//MbFNaxERGR8ydS6gFq7ZvmreO/rFtW6DBGRcSO5ewy1LkBEZJxKbDAMUJ+qIiJRiQ0GHXAWEakuucEw0Jik+zCIiEQkNhgGKBZERKKSGwxqShIRqSrWYDCzC81snZm1m9mVVeZ/wszWmNlTZnavmR0zVrUM5IJakkREomILBjNLAzcCFwFLgcvNbGnFYo8Dy9z91cB3gS+MeV1qTBIRiYhzj+F0oN3dN7p7DrgdWF6+gLvf7+494ejDwPyxKkZnJYmIVBdnMMwDNpeNbwmnDeUDwM+qzTCzK8xslZmt6ujoeEnFDJyVpKYkEZGocXnw2czeAywD/qnafHe/2d2XufuyWbNmvbzXelnPFhGZfOLsK2krsKBsfH44LcLMzgM+DZzj7v1jVYyakkREqotzj+FRYImZLTazOuAyYEX5Amb2GuDfgEvcfedYFqOzkkREqostGNy9AHwYuBNYC9zh7qvN7BozuyRc7J+AZuC/zOwJM1sxxOpGjc5KEhGJirXbbXdfCaysmHZ12fB5MdYS10uJiEwo4/LgcxwGc0E7DCIiEYkNhgHKBRGRqMQHg4iIRCU2GHyw123tM4iIlEtsMAxQLIiIRCU2GFz9bouIVJXcYNAN3EREqkpsMAxQMIiIRCU2GNSQJCJSXWKDYYC6xBARiUpsMKhLDBGR6pIbDOFfHWMQEYlKbDCIiEh1iQ0GtSSJiFSX2GBg8J7PaksSESmX4GAIKBZERKISGwxqShIRqS65wRD+VUuSiEhUYoNhgC5wExGJSmwwqClJRKS65AbD4FlJNS5ERGScSWwwDFAuiIhEJTYY1JQkIlJd4oNBTUkiIlGJDYaDlAwiIuUSGwy657OISHXJDQY1JYmIVJXYYBigXBARiUp8MIiISFRig+FgU5L2GUREyiU2GAYoFkREohIbDDorSUSkusQGwwC1JImIRCU2GNQlhohIdckNhvCv9hhERKJiDQYzu9DM1plZu5ldWWV+vZl9J5z/iJktGvOadPhZRCQitmAwszRwI3ARsBS43MyWViz2AWCvux8HXA98fqzq6ekvjNWqRUQmtDj3GE4H2t19o7vngNuB5RXLLAduC4e/C7zZxuhCg/9evR2AukxiW9NERKrKxPha84DNZeNbgDOGWsbdC2bWCcwAdpUvZGZXAFeEo11mtu4l1jTzrM9H150AM0HbnADa5mR4Odt8zFAz4gyGUePuNwM3v9z1mNkqd182CiVNGNrmZNA2J8NYbXOc7ShbgQVl4/PDaVWXMbMM0AbsjqU6EREB4g2GR4ElZrbYzOqAy4AVFcusAN4XDv8ecJ+7rjgQEYlTbE1J4TGDDwN3AmngFndfbWbXAKvcfQXwdeAbZtYO7CEIj7H0spujJiBtczJom5NhTLbZ9INcRETK6VxNERGJUDCIiEhEYoPhSN1zTCRmdouZ7TSzZ8qmTTezu83s2fDvtHC6mdmXw+1+ysxOK3vO+8LlnzWz91V7rfHAzBaY2f1mtsbMVpvZR8Ppk3mbG8zs12b2ZLjNfxtOXxx2H9MedidTF04fsnsZM7sqnL7OzC6ozRYNn5mlzexxM/tJOD6pt9nMNpnZ02b2hJmtCqfF+9l298Q9CA5+bwCOBeqAJ4Glta7rZWzP2cBpwDNl074AXBkOXwl8Phy+GPgZwT2KzgQeCadPBzaGf6eFw9NqvW1DbO9c4LRwuAVYT9DNymTeZgOaw+Es8Ei4LXcAl4XTbwI+GA7/OXBTOHwZ8J1weGn4ea8HFof/D9K13r4jbPsngG8BPwnHJ/U2A5uAmRXTYv1sJ3WPYTjdc0wY7v4gwVlc5cq7F7kNeEfZ9P/wwMPAVDObC1wA3O3ue9x9L3A3cOHYVz9y7r7N3X8TDh8A1hJcNT+Zt9ndvSsczYYPB95E0H0MHLrN1bqXWQ7c7u797v4c0E7w/2FcMrP5wFuBr4XjxiTf5iHE+tlOajBU655jXo1qGStz3H1bOLwdmBMOD7XtE/I9CZsLXkPwC3pSb3PYpPIEsJPgP/oGYJ+7D/QIWV5/pHsZYKB7mQm1zcCXgE8BpXB8BpN/mx24y8wes6D7H4j5sz0hu8SQkXF3N7NJd16ymTUD3wM+5u77ray/xcm4ze5eBE41s6nAD4ATalzSmDKztwE73f0xMzu31vXE6Cx332pms4G7zey35TPj+GwndY9hON1zTHQ7wl1Kwr87w+lDbfuEek/MLEsQCt909++Hkyf1Ng9w933A/cDrCJoOBn7gldc/VPcyE2mb3wBcYmabCJp73wTcwOTeZtx9a/h3J8EPgNOJ+bOd1GAYTvccE1159yLvA35UNv294dkMZwKd4S7qncBbzGxaeMbDW8Jp407Ybvx1YK27/3PZrMm8zbPCPQXMrBE4n+DYyv0E3cfAodtcrXuZFcBl4Rk8i4ElwK/j2YqRcfer3H2+uy8i+D96n7u/m0m8zWbWZGYtA8MEn8lniPuzXesj8LV6EBzNX0/QTvvpWtfzMrfl28A2IE/QlvgBgrbVe4FngXuA6eGyRnDDpA3A08CysvX8EcGBuXbg/bXersNs71kE7bBPAU+Ej4sn+Ta/Gng83OZngKvD6ccSfMm1A/8F1IfTG8Lx9nD+sWXr+nT4XqwDLqr1tg1z+8/l4FlJk3abw217MnysHvhuivuzrS4xREQkIqlNSSIiMgQFg4iIRCgYREQkQsEgIiIRCgYREYlQMIgcgZktMjM3s3F1o/nxWpdMfAoGkSPbTNCj6xMAZnZu+IU8M64CzOwBM/vXw9UlMlrUV5LIEXjQR9H2sVi3mWXdPf9SnjuWdUmyaY9BxpXwl/FXzOwfzWyXBTcgus7MUuH8TWb2F1We869l45vM7Gozu9XMDpjZZjO71MymmtntZtYV3rzkLcOsabDJJuzN9f5wVkc4/dZwOTOzT5nZBjPrteBmK++psp7Lzew+M+sF/tTMZpjZt81sS/i81Wb2/rLn3QqcA3wofL6H6zqkKcnMzrbgJjV9ZrbDzK4Pu30Z1vsbLvO7Ftz0pdfM9pjZ/5jZQG+ekgAKBhmP3g0UgNcDHwY+Blw6wnV8jKBbhNMIbuxyG8HNXlYCpwIPAv9pZg0jXO9m4J3h8KsImnI+Go7/PUF3JB8iuDnMtcC/mdlbK9ZxLfCVcJkfEnTl8BvgbeE6bwif9+Zw+Y8CvwL+X/h6c4l2qQyAmc0juGnL4wRdkX8AuDx8vXJDvr9mdhRBh3W3AScS3ATqG8N4X2QyqXXfIHroUf4AHgB+VTHtbuBr4fAm4C+qPOdfy8Y3Ad8uG28m6Fvpy2XTFoXTlg2jpsiyBP32OGV32QKagF7gjRXP/RKwsmI9nxzGa94+sM3VtnGIuv6BoC+dVNkyfwj0A1OG+f6eFq7zmFp/FvSo3UPHGGQ8eqpi/EVg9ktdh7t3mVkPQSdjA3aEf0e63qEsJfjl/98W7Ss/SxBU5VaVj5hZmuB2jZcS3EylnuCWsw+MsIYTgYfdvVQ27Rfhuo7j4HtyuPf3SYJO2p4xs7vC4e+6e8cIa5EJTMEg41HlwVjnYLNniaBHyXLZYa4jXzEOo9ecOrCetwMvHKGW7orxvwA+SdBk9DTQBfwjoxdacHB7q9Uz+P66ezE89nImQVfNHwCuNbNz3P3JUaxHxjEFg0w0HQRt7ACExwhOIGhXj0su/Jsum7aGoMnmGHe/b4TrO0mvud8AAAFqSURBVAv4sbt/AwbvN3E8sK/iNdNVnltuLfAuM0uV7TWcFT53w3CLcXcnOKbxKzO7hqD750sJ9iYkARQMMtHcB/yRma0gCIlPE//n+HmCX9lvNbMfA73ufsDMrgOuC7/YHyQ4tnEmUHL3mw+zvvXApWZ2FrAL+AiwmGjYbQJOD8+K6gL2VFnPVwgOJH/FzG4g6Nv/cwTHJnqGs2HhzV7OI7ipyw6Cg9gLCIJPEkJnJclEcy1BOPwIuIugDT3OvQU8uPXiZwgO9u4ABk6V/RvgswRNQ6sJDuq+E3juCKv8e4IzqH5GECjdwDcrlrmO4Jf/GoJAXDhEXRcRfJk/AdxCcBOnvxrB5nUS3FLzJwQHsr8I/J27/+cI1iETnG7UIyIiEdpjEBGRCAWDJJ6Z3RReDV3tcVOt6xOJm5qSJPHMbDbQOsTs/e6+M856RGpNwSAiIhFqShIRkQgFg4iIRCgYREQkQsEgIiIR/x+6Z5qi4rK/RgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Instantiate learner in VW\n",
    "vw = pyvw.vw(\"--cb_explore_adf -q UA --quiet --epsilon 0.2\")\n",
    "\n",
    "num_iterations = 5000\n",
    "ctr = run_simulation(vw, num_iterations, users, times_of_day, actions, get_cost)\n",
    "\n",
    "plot_ctr(num_iterations, ctr)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Aside: interactions\n",
    "You'll notice in the arguments we supply to VW, **we include `-q UA`**. This is telling VW to create additional features which are the features in the (U)ser namespace and (A)ction namespaces multiplied together. This allows us to learn the interaction between when certain actions are good in certain times of days and for particular users. If we didn't do that, the learning wouldn't really work. We can see that in action below."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAEOCAYAAACNY7BQAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO3deXxcdb3/8ddnJnvSNEmTdN+ghVIUAUNBAeHBIgW94O+qV1Afbvwu159y1Z/L/YELKoKioiJXEbleLy5XEHErWARkVS9bWFpoC6Ut3fclaZNmmcx8fn+ck3ROOmmTtjlZ5v18PObROcuc+Zxpkvec7/ec7zF3R0REpFtiqAsQEZHhRcEgIiIRCgYREYlQMIiISISCQUREIhQMIiISEVswmNlPzWyrmb3Ux3Izs5vNbIWZLTazk+OqTURE9onziOF2YP4Bll8IzA4fVwA/iqEmERHpJbZgcPfHgZ0HWOUS4OceeBKoMrOJ8VQnIiLdCoa6gCyTgXVZ0+vDeZt6r2hmVxAcVVBeXv7GOXPmxFKgiMho8eyzz25397pcy4ZTMPSbu98G3AbQ0NDgjY2NQ1yRiMjIYmZr+lo2nM5K2gBMzZqeEs4TEZEYDadgWAB8IDw76TSg2d33a0YSEZHBFVtTkpndAZwN1JrZeuDLQCGAu98KLAQuAlYAe4EPx1WbiIjsE1swuPtlB1nuwMdjKkdERPownJqSRERkGFAwiIhIhIJBREQiFAwiIhKhYBARkQgFg4iIRCgYREQkQsEgIiIRCgYREYlQMIiISISCQUREIhQMIiISoWAQEZEIBYOIiEQoGEREJELBICIiEQoGERGJUDCIiEiEgkFERCIUDCIiEqFgEBGRCAWDiIhEKBhERCRCwSAiIhEKBhERiVAwiIhIhIJBREQiFAwiIhKhYBARkQgFg4iIRCgYREQkQsEgIiIRCgYREYlQMIiISISCQUREImINBjObb2avmNkKM7sqx/JpZvaImT1vZovN7KI46xMRkRiDwcySwA+BC4G5wGVmNrfXal8E7nL3k4BLgVviqk9ERAJxHjHMA1a4+yp37wTuBC7ptY4DleHzscDGGOsTERHiDYbJwLqs6fXhvGxfAd5vZuuBhcC/5tqQmV1hZo1m1rht27bBqFVEJG8Nt87ny4Db3X0KcBHwCzPbr0Z3v83dG9y9oa6uLvYiRURGsziDYQMwNWt6Sjgv2+XAXQDu/gRQAtTGUp2IiADxBsMzwGwzm2lmRQSdywt6rbMWOBfAzI4jCAa1FYmIxCi2YHD3LuBK4H5gGcHZR0vM7Fozuzhc7TPAP5vZIuAO4EPu7nHVKCIiUBDnm7n7QoJO5ex512Q9XwqcHmdNIiISNdw6n0VEZIgpGEREJELBICIiEQoGERGJUDCIiEiEgkFERCIUDCIiEqFgEBGRCAWDiIhEKBhERCRCwSAiIhEKBhERiVAwiIhIhIJBREQiFAwiIhKhYBARkQgFg4iIRCgYREQkQsEgIiIRCgYREYlQMIiISISCQUREIhQMIiISoWAQEZEIBYOIiEQoGEREJELBICIiEQoGERGJUDCIiEiEgkFERCIUDCIiEqFgEBGRCAWDiIhEKBhERCQi1mAws/lm9oqZrTCzq/pY55/MbKmZLTGzX8VZn4iIQEFcb2RmSeCHwPnAeuAZM1vg7kuz1pkNXA2c7u67zKw+rvpERCQQ5xHDPGCFu69y907gTuCSXuv8M/BDd98F4O5bY6xPRESINxgmA+uypteH87IdAxxjZn83syfNbH6uDZnZFWbWaGaN27ZtG6RyRUTy03DrfC4AZgNnA5cB/2FmVb1Xcvfb3L3B3Rvq6upiLlFEZHSLMxg2AFOzpqeE87KtBxa4e8rdXwOWEwSFiIjEJM5geAaYbWYzzawIuBRY0GudPxAcLWBmtQRNS6tirFFEJO/FFgzu3gVcCdwPLAPucvclZnatmV0crnY/sMPMlgKPAJ9z9x1x1SgiImDuPtQ1HJaGhgZvbGwc6jJEREYUM3vW3RtyLRtunc8iIjLEFAwiIhKhYBARkQgFg4iIRCgYREQkIq+DobMrw1fvWcKu1s6hLkVEZNjI62C476VN/NffV3P9wmVDXYqIyLCR18HQfQlHKp0Z2kJERIaRvA4Gs+DfzMi+xk9E5IjqVzCYWaGZbTaz4we7oKEw0q/+FhE5kvoVDO6eAlLAqPoLauEhw6jaKRGRwzSQpqR/B642s9huBzrYEmFTkpJBRGSfgfyRPxM4C9hgZi8BrdkL3f3inK8axowgGTJqShIR6TGQYNgO/HawChkK3UcMygURkX36HQzu/uHBLGQo7DsrSckgItKt330MZvZwrvsvm1mlmT18ZMuKizqfRUR6G0jn89lAUY75JQT9DyOOqSlJRGQ/B21KMrOTsyZPMLOdWdNJ4AJgw5EuLA4J02lJIiK99aePoZHgL6cDD+RY3gb865EsKm668llEZJ/+BMNMgsb4VcApBGcndesEtrp7ehBqG3TpMBF05bOIyD4HDQZ3XwNgZjcAp7j7rdnLzeyjZjbZ3b80SDUOmu5AUCyIiOwzkM7n9wPP55j/HPCBI1NOvNLdwaBkEBHpMZBgqAe25Zi/HRh/ZMqJV3ffgq5jEBHZZyDBsJbcp6W+BVh/ZMqJV2aAvc7pjLN1T/sgVSMiMjwMJBh+DHzPzP7ZzI4OH1cA3wFuG5zyBlemj6akLbvb+fIfX+oJDnfn3sUbOfrzC5l3/UPsaOmIu1QRkdgMZEiM75hZLXAz+y506wS+7+7fGoziBlvPWUm9up9P/fpDAJwxu47zjqvn2C/9mc6ufXd5+8o9SzlpahXX3ruUn31kHmcdUxdf0SIig2xAd3Bz96uBWuC08FHn7lcNRmFx6D5S6KuLYU97iq17OiKhAHDPoo1ce+9SAD7406cjtwa9Z9FGNjW3DUq9IiJxGPC9Fdy9FXhmEGqJXfdZSX11Pn/6rkX8wxsmHXQ7X16whOa2FH9avCky//kvnU91ea5RREYed++5sdFAbNvTQd2Y4p7pNTtaMYwp1aUkEgPfnogMvlFz051D0VcfQ7Z7Fm2MTL/31Gn86qm1kXm9p7ud9LUH+fa7TuDdDVMPr9AYvLplD9cvXMaMceWcM6eeqrJC5kyo5Ka/LOeWR1cC8Njnzmb6uPKDbqs9lWbNjr184fcv0rhmV5/rzaqvoKq0kMY1u6goLuDf5h/Le+dNI5V2igsSCg6RIZLfwRD2MWzMavrp6Mp9EXfjF89jY1Mb02vKe4Lg7SdM5N5eRwm9fe7uxbx17gTGlhUeoaqPvBvvf4UfPLIinNrG7f+zOud6Z3370cj0KTOqedcbp7CpuZ2b/vJqv97rxKlVpDPOixuaWbG1pWd+S0cX1/xxCdf8cUlk/ak1pRQkErz1+PEcXVvBlJpSfv/cBl7b3sqEsSW8/YRJXHD8+EM6mhGR3PI7GMIjhXU723hteyvTaspo7dg/GCZUllBbUUxtRdAk8tGzjubWx1Zy47vfwGXzpvG+nzwFwFF15Tz06bP4f79dzOVnHMUFNz0OwBuufYCnP38um5rb+dSvX+Duj76JcRXF+71PXDIZ57sPLs8Kg8BxEytZtml3ZN784ydw1YVzWLW9hevuXcaq7ftu3PfM6l08szr3EcG73jiFM2fXcsmJk/us46lVOxhfWcKM2nK2t3TwH39dxS+eWENZUZLtLZ1A8H8D8OPHVkVeO668iMY1u7h38SZm11fw3lOncd5x45laU9b/D0JEcrKRPk5QQ0ODNzY2HtJrf/q313o6kQH+5ayjeP+p0znzW49w8RsmsSBsRvrdx97MydOqc27D3fnfP2vkoZe30vjF83rCAyCVzjD7C/flfN0D//ctHDN+zCHVnb39J1ft4IxZtf3+xpxKZzjmi/ft13x23yfP5LiJlf3ahruzs7WTexZt5EePrWRCZQl3ffRNNO1NUVyQoKrsyPartHWmWbS+ibZUmqdf28n7T5vO5KpSUukMv3hiDbc8upLt4SnEp8yoprqsiG0tHXzozTO4+A2TdDQhkoOZPevuDTmX5XMw/OSvq7juT8t6pk+YMpYb/vEELrr5r3zinFnc/HDwjXrV1y865PbulzfvZv5Nf825bPUNbzukbQJ0pTOc/73HeS38Bl9SmOBDb57JrY8F/QHzj5/An5ds5tw59Tz08tac2/jEubO5/IyZVBQXkBzB7fmpdIa/r9jO/Us28/DLW9mye991JidNq+KcY4M+k6WbdjNvZg1nHVNPTT9PCmhuS7G9pYOjast7AsbdeXVrC5ua28lknKKCBCWFSTLunDBlLMUFyQNus6Wjiy272ykpTLKzpZMlG5spSCY4flIlzW0pCpMJXje58qDbETkcBwqGPG9KioZiV9r546Lg1hLTsjpZD6cTdM6ESlbf8DZ+07iOna2dNMyo5p0/egKAGVf9iU+cO5ubH9rXPn/7h0/h7GPrD7jNrXvamXf9Q5F57alMTygA/HnJZoCcoTCrvoL7P/WWER0G2QqTCc4+tj7yuXWlM/zyyTV87y+v8vza5T3z73h6XfgaY8LYEk6aWk1pYZJdezupKS8inXE2NbezZGMzyUSi50hkclUpRQUJ9nZ2kfHgbKu+lBUlmVZTRm1FMWt2trKluYPOdIaqskJKCpJs3t2/q+fLi5KcPquWdMbZ0dpJZ1eGo+srOO2oGi45cTIVxXn96zuidZ/l5+5s2xP8fOxu62J3e4pkwtjY1MbyLXtoT2XYsKuNVDpDWypNUUGChBnpjFNdVsi7G6Zy+qzaI15fXh8x3PrYSm647+We6ak1pT1t2vd98kwu/H7wTf9wvtnnsmZH634dudme/eJ5VJcVYcZ+zSB/fmkTH/3lcz3Tv/vYm7l30SZeWLeL59Y28enzj+FHj66kLZXm0lOm8urWFtydydVlfPr8Y5hZe/CzikaTTMbZtLud17a1Mqu+gle27OHx5dvY3tLBpuZ2XlzfTFsqTVVZIU17UxQkjGk1ZUwfV0ZZUQETxpZQUpjgf1buoK0zzZTqMszg1Jk11FeWUFyQ6Pklb+tM8/zaXbR0pHlpQzPrd+1l1vgxHFNfQUHS2NWaoi2VZnJ1KUfXVZA0KC1KMqu+gozDK5v3kDDDDFZvb2X9rjYWb2jCHarLiihMGss27aEtFfSDvfnocUypLmVGbTmnzhzHMeMrGFMyfE9yGC3aU2nSGWdPexepdIZdezupLitiY1MbG5vbWL6lhYriAsygaW+KHS2d7GjtIJ1xXt68h+a2FF3pDOVFBaQyGdpTmZzvk0wYhUlj4thSigsSFCQN9+AsymTC2LW3k89dcOwB+/EOZNgcMZjZfOD7BHd++4m739DHeu8E7iYY5vvQ/ur3Q7rXWEndoQDBL+JgmT6unLs/+ibedWtw5HD6rHGcMauO7z24nM50hjde95eedVdcfyEJM95xy99ZvL45sp2VX7+IZML26//4xLmzB632kSaRMCZXlTK5qhSACWNL+rxSvSudoSA5oGs+9/OOkw7tlxTglBk1B13H3XlubRMLXtjAw69s5cUNzexp7+pZPru+gtNn1XLqzBpOmVkT6fPKF+5OKh2Ma9aeSgNGKp1hd1uKMSWFtKXS7GjpoLgwyZKNzaS6nJaOFGt37mXrng7KipK0dKQpK0ySTBjNbSl2t6do7eiirTNNa+eBbz9TkDC6wr8tRckElaWF1I0pJpkIwnxCZQmJhLGnPUXSjJm15RQXJikpTFBTXowB9ZXFwfwhak6M7YjBzJLAcuB8gkH3ngEuc/elvdYbA/yJYNiNKw8WDIdzxPCDh1/lxgeW51z28tfm09yWCv+TSg5p+wfTHUzZTTqv+/L9tHR09fUSAN536jSu/1+vH5SaZOTZ2NTG06/tZOmm3Ty/dhcvrGsilXYSBqfPquXougpm1Vdw4tQq5kwYc9jhF4c1O1pJmJFIGFt2t9O0t5PdbV2MryxhQ1MbhUmjpaOLJ1ftZPnmPexNdVFVGjQFrt2596C/Q72ZwdTqMiaOLaGjK0NpYZL2rjQZh8qSAsaWFlJRXEBRQYKa8iIMo6aiiKQZ4yqKaNrbyfjKEiZVlTJjXDkZd7oyTnlRctie/DBcjhjmASvcfVVY1J3AJcDSXut9Dfgm8LnBLqivUAAoKUxSUji4aZ2rjX/hJ87k141rOWdOfU9fRLab3nMiF/fjamzJH5OqSnnHSZN7jlZ2tHTw/Nomnlm9kz+8sIG/r9jec2p2TXkRFxw/gbe9fiKnHlVD4REKCXenLZWmNPyd6Uxn6OzKsLm5nZ2tnbR0dLGxqY2yogJaO7v4+4rtNO1NsW7nXsqKC2jt6KKoIMG2PR0kzPr9h71uTDHHT6qkoriA3e1dpDMZXj95IhOrSqguKwoulDSjpChJdVkhO1s7qSwppKwoSWc6w/GTxlJWlKS4IDFs/4APhTiDYTKwLmt6PXBq9gpmdjIw1d3/ZGZ9BkM4qusVANOmTTvihU4bwnPhp40r43MXzAGCvo2fP7GamvIi5h8/YUR805OhN66imPPmjue8ueO5+qLjcHdWbmvhqdd28rdXt/P759dzx9Nrqa0o5uxj6ziqrpwTp1ZxdF0FdRXFOU+2SGec17a34h4MObl+115aO9K8uKGZx5dvY8OuNvZ0dFFckKAznTnoza/qxhQzqaqUN86ooa0zTVlRcFZXbUUxZkFnf1lRAal0hinVpVSVFTK2tJANTe1MGluCE5x0ML2mTFfID4Jhc1qDmSWA7wIfOti67n4b4VDfDQ0NR7wt7OoL5xzpTR6yD7xpxlCXICOcmTGrfgyz6sfwvlOn09yW4tFXtnLv4k08sGQzu7P6KEoLkxw3cQylRcmeTtHNze1saMo9MKRZ0Dcyd1IlM8eVs2tvisKkUVKYDDryw078ypICxlcGzTRlRUkmVx3aWFmz6g/v2h/pnziDYQOQPWjQlHBetzHA64BHw0O6CcACM7t4MDugc3nj9NwXs4mMBmNLC7nkxMk9Z7Ns3d3OSxuDIUpe3ryHpRt305bKUFGcxDCOn1TJ298wkWk1ZSTNKEgmmFRVwtjSQiaNLR01A0XKPnEGwzPAbDObSRAIlwLv7V7o7s0EQ3oDYGaPAp+NOxRg8DqbRYaj+soSzqks4Zw5I/IOvTIIYmu0dvcu4ErgfmAZcJe7LzGza83s4rjqEBGRA4u1j8HdFwILe827po91z46jJhERidJpLiIiEqFg6EWnMotIvhs2p6sOpX95y1H8+PFVPPyZszTWjIjkPQUDcNWFc/jMW4+lqEAHUCIiCgaCC4CKCtSGJCIC6mMQEZFeFAwiIhKRt8Ew0GF5RUTyRd4Gw3X39h7tW0REII+DoWlvaqhLEBEZlvI2GHQhm4hIbgoGERGJyN9gQMkgIpJL3gaDckFEJLf8DQYREckpb4Oh+4BhWk3ZkNYhIjLc5G0wdFMntIhIVN4GgykRRERyyttgEBGR3PI2GHS8ICKSW/4Gg5JBRCSn/A2GoS5ARGSYyttg8KEuQERkmMrbYOimIwcRkai8DwYREYlSMIiISETeB4P6GkREovI2GNS3ICKSW94GQzfXIYOISETeBkMivMLN1ZgkIhKRt8HQ3ZaUyQxtGSIiw03eBoNu7Skiklv+BoNyQUQkp7wNhm6u3mcRkYhYg8HM5pvZK2a2wsyuyrH802a21MwWm9lDZjZ9sGpJhEcMigURkajYgsHMksAPgQuBucBlZja312rPAw3ufgJwN/CtQasn7GPI6IhBRCQiziOGecAKd1/l7p3AncAl2Su4+yPuvjecfBKYMljFdPcxKBdERKLiDIbJwLqs6fXhvL5cDtyXa4GZXWFmjWbWuG3btkMqRp3PIiK5DcvOZzN7P9AAfDvXcne/zd0b3L2hrq7usN5LBwwiIlEFMb7XBmBq1vSUcF6EmZ0HfAE4y907Bq+c8MpnJYOISEScRwzPALPNbKaZFQGXAguyVzCzk4AfAxe7+9bBLCbR05SkZBARyRZbMLh7F3AlcD+wDLjL3ZeY2bVmdnG42reBCuA3ZvaCmS3oY3OHrbuPIaNcEBGJiLMpCXdfCCzsNe+arOfnxVWLhsQQEcltWHY+x0lXPouIROVtMJiufBYRySl/gyH8VwcMIiJR+RsM3TfqUTKIiETkbTB0UyyIiETlbTD0DImhZBARicjbYKguKwLghKljh7gSEZHhJW+DYVJVKQBfvfh1Q1yJiMjwkrfB0N3pXFyQtx+BiEhO+qsoIiIRCgYREYlQMIiISISCQUREIvI2GHT5gohIbnkbDN1072cRkai8DwYREYlSMIiISISCQUREIvI3GNT7LCKSU/4GQ8jU+ywiEpH3wSAiIlEKBhERiVAwiIhIRN4Gg6v3WUQkp7wNhm7qehYRicr7YBARkSgFg4iIRORtMLi6GEREcsrbYOim69tERKLyPhhERCRKwSAiIhEKBhERicjbYFDfs4hIbnkbDN1Ml7iJiETkfTCIiEhUrMFgZvPN7BUzW2FmV+VYXmxmvw6XP2VmM+KsT0REYgwGM0sCPwQuBOYCl5nZ3F6rXQ7scvdZwPeAbw5WPXs704O1aRGRES3OI4Z5wAp3X+XuncCdwCW91rkE+Fn4/G7gXBukW6y9uL4JgIKk+hhERLIVxPhek4F1WdPrgVP7Wsfdu8ysGRgHbM9eycyuAK4IJ1vM7JVDrKm27pvRbeeBWtA+5wHtc344nH2e3teCOIPhiHH324DbDnc7Ztbo7g1HoKQRQ/ucH7TP+WGw9jnOpqQNwNSs6SnhvJzrmFkBMBbYEUt1IiICxBsMzwCzzWymmRUBlwILeq2zAPhg+PxdwMPuGgdVRCROsTUlhX0GVwL3A0ngp+6+xMyuBRrdfQHwn8AvzGwFsJMgPAbTYTdHjUDa5/ygfc4Pg7LPpi/kIiKSTVc+i4hIhIJBREQi8jYYDjY8x0hiZj81s61m9lLWvBoze9DMXg3/rQ7nm5ndHO73YjM7Oes1HwzXf9XMPpjrvYYDM5tqZo+Y2VIzW2Jmnwznj+Z9LjGzp81sUbjPXw3nzwyHj1kRDidTFM7vc3gZM7s6nP+KmV0wNHvUf2aWNLPnzezecHpU77OZrTazF83sBTNrDOfF+7Pt7nn3IOj8XgkcBRQBi4C5Q13XYezPW4CTgZey5n0LuCp8fhXwzfD5RcB9gAGnAU+F82uAVeG/1eHz6qHetz72dyJwcvh8DLCcYJiV0bzPBlSEzwuBp8J9uQu4NJx/K/B/wucfA24Nn18K/Dp8Pjf8eS8GZoa/B8mh3r+D7PungV8B94bTo3qfgdVAba95sf5s5+sRQ3+G5xgx3P1xgrO4smUPL/Iz4B1Z83/ugSeBKjObCFwAPOjuO919F/AgMH/wqx84d9/k7s+Fz/cAywiumh/N++zu3hJOFoYPB84hGD4G9t/nXMPLXALc6e4d7v4asILg92FYMrMpwNuAn4TTxijf5z7E+rOdr8GQa3iOyUNUy2AZ7+6bwuebgfHh8772fUR+JmFzwUkE36BH9T6HTSovAFsJftFXAk3u3hWukl1/ZHgZoHt4mRG1z8BNwL8BmXB6HKN/nx14wMyetWD4H4j5Z3tEDokhA+Pubmaj7rxkM6sAfgt8yt13W9Z4i6Nxn909DZxoZlXA74E5Q1zSoDKztwNb3f1ZMzt7qOuJ0RnuvsHM6oEHzezl7IVx/Gzn6xFDf4bnGOm2hIeUhP9uDef3te8j6jMxs0KCUPhvd/9dOHtU73M3d28CHgHeRNB00P0FL7v+voaXGUn7fDpwsZmtJmjuPQf4PqN7n3H3DeG/Wwm+AMwj5p/tfA2G/gzPMdJlDy/yQeCPWfM/EJ7NcBrQHB6i3g+81cyqwzMe3hrOG3bCduP/BJa5+3ezFo3mfa4LjxQws1LgfIK+lUcIho+B/fc51/AyC4BLwzN4ZgKzgafj2YuBcfer3X2Ku88g+B192N3fxyjeZzMrN7Mx3c8JfiZfIu6f7aHugR+qB0Fv/nKCdtovDHU9h7kvdwCbgBRBW+LlBG2rDwGvAn8BasJ1jeCGSSuBF4GGrO18hKBjbgXw4aHerwPs7xkE7bCLgRfCx0WjfJ9PAJ4P9/kl4Jpw/lEEf+RWAL8BisP5JeH0inD5UVnb+kL4WbwCXDjU+9bP/T+bfWcljdp9DvdtUfhY0v23Ke6fbQ2JISIiEfnalCQiIn1QMIiISISCQUREIhQMIiISoWAQEZEIBYPIQZjZDDNzMxtWN5ofrnXJyKdgEDm4dQQjur4AYGZnh3+Qa+MqwMweNbMfHKgukSNFYyWJHIQHYxRtHoxtm1mhu6cO5bWDWZfkNx0xyLASfjO+xcy+bmbbLbgB0Y1mlgiXrzazz+Z4zQ+ypleb2TVmdruZ7TGzdWb2HjOrMrM7zawlvHnJW/tZU0+TTTia6yPhom3h/NvD9czM/s3MVppZmwU3W3l/ju1cZmYPm1kb8C9mNs7M7jCz9eHrlpjZh7NedztwFvDx8PUebmu/piQze4sFN6lpN7MtZva9cNiXfn2+4Tr/aMFNX9rMbKeZPWZm3aN5Sh5QMMhw9D6gC3gzcCXwKeA9A9zGpwiGRTiZ4MYuPyO42ctC4ETgceCXZlYywO2uA94ZPj+eoCnnk+H0dQTDkXyc4OYw3wB+bGZv67WNbwC3hOv8gWAoh+eAt4fb/H74unPD9T8JPAH8V/h+E4kOqQyAmU0muGnL8wRDkV8OXBa+X7Y+P18zm0AwYN3PgOMIbgL1i358LjKaDPXYIHrokf0AHgWe6DXvQeAn4fPVwGdzvOYHWdOrgTuypisIxla6OWvejHBeQz9qiqxLMG6Pk3WXLaAcaAPO7PXam4CFvbbzmX68553d+5xrH/uo63qCsXQSWet8COgAyvr5+Z4cbnP6UP8s6DF0D/UxyHC0uNf0RqD+ULfh7i1mtpdgkLFuW8J/B7rdvswl+Ob/Z4uOlV9IEFTZGrMnzCxJcLvG9xDcTKWY4Jazjw6whuOAJ909kzXvb+G2ZrHvMznQ57uIYJC2l8zsgfD53e6+bYC1yAimYJDhqHdnrLOv2TNDMKJktsJ+biPVaxqOXHNq93b+AYg/QoYAAAHCSURBVFh7kFpae01/FvgMQZPRi0AL8HWOXGjBvv3NVU/P5+vu6bDv5TSCoZovB75hZme5+6IjWI8MYwoGGWm2EbSxAxD2EcwhaFePS2f4bzJr3lKCJpvp7v7wALd3BnCPu/8Ceu43cQzQ1Os9kzlem20Z8E9mlsg6ajgjfO3K/hbj7k7Qp/GEmV1LMPzzewiOJiQPKBhkpHkY+IiZLSAIiS8Q/8/xGoJv2W8zs3uANnffY2Y3AjeGf9gfJ+jbOA3IuPttB9jecuA9ZnYGsB34V2Am0bBbDcwLz4pqAXbm2M4tBB3Jt5jZ9wnG9r+BoG9ib392LLzZy3kEN3XZQtCJPZUg+CRP6KwkGWm+QRAOfwQeIGhDj/NoAQ9uvfhlgs7eLUD3qbJfAr5C0DS0hKBT953AawfZ5HUEZ1DdRxAorcB/91rnRoJv/ksJAnFaH3VdSPDH/AXgpwQ3cfr8AHavmeCWmvcSdGR/B/iau/9yANuQEU436hERkQgdMYiISISCQfKemd0aXg2d63HrUNcnEjc1JUneM7N6oLKPxbvdfWuc9YgMNQWDiIhEqClJREQiFAwiIhKhYBARkQgFg4iIRPx/bIrPVJJbnHAAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Instantiate learner in VW but without -q\n",
    "vw = pyvw.vw(\"--cb_explore_adf --quiet --epsilon 0.2\")\n",
    "\n",
    "num_iterations = 5000\n",
    "ctr = run_simulation(vw, num_iterations, users, times_of_day, actions, get_cost)\n",
    "\n",
    "plot_ctr(num_iterations, ctr)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "\n",
    "### Without learning\n",
    "Let's do the same thing again (but with `-q`, but this time show the effect if we don't learn from what happens. The ctr never improves are we just hover around 0.2."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAEOCAYAAACNY7BQAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAfHklEQVR4nO3de5RcZZnv8e9Tl66+Jp1LJ4EkJFHCQHAQY8ugoCA4XBVmHVwKMw7KcMQ5Iy5d6pkF4xEd5Igo6ugRD3K8gMwIojNqhCigIMwoYBqBQAiBEBJyv3bS1+q6PeePvTup3alOukn37kv9PmvV6tqX2vW83dX1q73fXe82d0dERKRfYqwLEBGR8UXBICIiEQoGERGJUDCIiEiEgkFERCIUDCIiEhFbMJjZ981sh5k9N8hyM7NvmtlaM1tpZkvjqk1ERA6Ic4/hduC8Qyw/H1gc3q4C/m8MNYmIyACxBYO7PwrsOcQqFwM/9MDjQLOZHRVPdSIi0i811gWUmQtsLJveFM7bOnBFM7uKYK+ChoaGNx9//PGxFCgiMlk8+eSTu9y9pdKy8RQMQ+butwG3AbS2tnpbW9sYVyQiMrGY2YbBlo2ns5I2A/PLpueF80REJEbjKRiWAZeHZyedCuxz94MOI4mIyOiK7VCSmd0FnAnMNLNNwOeANIC73wosBy4A1gI9wBVx1SYiIgfEFgzuftlhljvw0ZjKERGRQYynQ0kiIjIOKBhERCRCwSAiIhEKBhERiVAwiIhIhIJBREQiFAwiIhKhYBARkQgFg4iIRCgYREQkQsEgIiIRCgYREYlQMIiISISCQUREIhQMIiISoWAQEZEIBYOIiEQoGEREJELBICIiEQoGERGJUDCIiEiEgkFERCIUDCIiEqFgEBGRCAWDiIhEKBhERCRCwSAiIhEKBhERiVAwiIhIhIJBREQiFAwiIhKhYBARkQgFg4iIRCgYREQkQsEgIiIRsQaDmZ1nZmvMbK2ZXVNh+TFm9rCZPWVmK83sgjjrExGRGIPBzJLALcD5wBLgMjNbMmC1/wXc4+5vAi4Fvh1XfSIiEohzj+EUYK27r3P3HHA3cPGAdRyYEt6fCmyJsT4RESHeYJgLbCyb3hTOK/d54ANmtglYDnys0obM7CozazOztp07d45GrSIiVWu8dT5fBtzu7vOAC4A7zeygGt39NndvdffWlpaW2IsUEZnM4gyGzcD8sul54bxyVwL3ALj7Y0AtMDOW6kREBIg3GFYAi81skZnVEHQuLxuwzqvA2QBmdgJBMOhYkYhIjGILBncvAFcD9wOrCc4+WmVm15vZReFqnwI+bGbPAHcBH3J3j6tGERGBVJxP5u7LCTqVy+ddV3b/eeC0OGsSEZGo8db5LCIiY0zBICIiEQoGERGJUDCIiEiEgkFERCIUDCIiEqFgEBGRCAWDiIhEKBhERCRCwSAiIhEKBhERiVAwiIhIhIJBREQiFAwiIhKhYBARkQgFg4iIRCgYREQkQsEgIiIRCgYREYlQMIiISISCQUREIhQMIiISoWAQEZEIBYOIiEQoGEREJELBICIiEQoGERGJUDCIiEiEgkFERCIUDCIiEqFgEBGRCAWDiIhEKBhERCRCwSAiIhGxBoOZnWdma8xsrZldM8g67zOz581slZn9KM76REQEUnE9kZklgVuAvwQ2ASvMbJm7P1+2zmLgWuA0d283s1lx1SciIoE49xhOAda6+zp3zwF3AxcPWOfDwC3u3g7g7jtirE9ERIg3GOYCG8umN4Xzyh0HHGdmvzezx83svEobMrOrzKzNzNp27tw5SuWKiFSn8db5nAIWA2cClwH/z8yaB67k7re5e6u7t7a0tMRcoojI5BZnMGwG5pdNzwvnldsELHP3vLu/ArxIEBQiIhKTOINhBbDYzBaZWQ1wKbBswDo/J9hbwMxmEhxaWhdjjSIiVS+2YHD3AnA1cD+wGrjH3VeZ2fVmdlG42v3AbjN7HngY+J/uvjuuGkVEBMzdx7qGI9La2uptbW1jXYaIyIRiZk+6e2ulZeOt81lERMaYgkFERCIUDCIiEqFgEBGRCAWDiIhEKBhERCRCwSAiIhEKBhERiVAwiIhIxJCCwczSZrbNzE4c7YJERGRsDSkY3D0P5IGJPX6GiIgc1nAOJf0f4Fozi+1yoCIiEr/hvMm/HTgD2GxmzwHd5Qvd/aKKjxIRkQllOMGwC/j30SpERETGhyEHg7tfMZqFiIjI+DDkPgYze6jS9ZfNbIqZPTSyZY2+57d0cMUP/kiuUBrrUkRExpXhdD6fCdRUmF9L0P8woVzzHyt5eM1OVm/tGOtSRETGlcMeSjKzpWWTJ5nZnrLpJHAusHmkCxMRkbExlD6GNoLvLzjwQIXlvcDHRrKoOOmLGSIiUUMJhkWAAeuAtxCcndQvB+xw9+Io1DaqbKwLEBEZpw4bDO6+AcDMvgS8xd1vLV9uZn9vZnPd/bOjVOOoctc+g4hIueF0Pn8AeKrC/D8Bl49MOTEy7TOIiFQynGCYBeysMH8XMHtkyomR9hRERCoaTjC8SuXTUt8BbBqZckREZKwNZ0iM7wBfN7MaoP8LbWcDNwI3jXRho06HkkREKhrOkBhfNbOZwDc58EW3HPANd//yaBQXBx1QEhGJGtYQ2u5+rZndACwJZ612966RL2v0aX9BRKSyYV9bwd27gRWjUIuIiIwDVX/NZ52cJCISVfXBICIiUQoGERGJqNpg0NmqIiKVVW0wHKBOBhGRclUbDNphEBGprGqDQUREKos1GMzsPDNbY2ZrzeyaQ6x3iZm5mbWOdk06XVVEJCq2YDCzJHALcD7BN6cvM7MlFdZrAj4OPBFXbSIickCcewynAGvdfZ2754C7gYsrrPcFgkH5sjHWJiIioTiDYS6wsWx6UzhvPzNbCsx39/sOtSEzu8rM2sysbefOSpeIODzT+aoiIhWNm85nM0sAXwM+dbh13f02d29199aWlpbRL05EpIrEGQybgfll0/PCef2agDcAvzOz9cCpwLLR7oBW37OISFScwbACWGxmi8KL/VwKLOtf6O773H2muy9094XA48BF7t42GsXoQJKISGWxBYO7F4CrgfuB1cA97r7KzK43s4viqmOgDbt7xuqpRUTGpWFfj+FIuPtyYPmAedcNsu6Zo1lLV18BgE//5Bne++Z5o/lUIiITyrjpfI6bvtgmIlJZ1QaDzlYVEamsioNBySAiUkn1BsNYFyAiMk5VbTCIiEhlCgYREYlQMIiISETVBkNJ56uKiFRUtcGgs5JERCqr3mAY6wJERMapqg0GERGpTMEgIiIRVRsM6mIQEamsaoNBREQqUzCIiEhE1QaDDiWJiFRWvcGgE1ZFRCqq2mAQEZHKqjYYdChJRKSyqg0GERGpTMEgIiIRVRsMOpIkIlJZ1QaDiIhUVr3BoN5nEZGKqjcYRESkoqoNBu0viIhUVrXBICIilSkYREQkomqDQX3PIiKVVW0wiIhIZVUbDNphEBGprHqDQceSREQqqtpgEBGRyqo2GLS/ICJSWazBYGbnmdkaM1trZtdUWP5JM3vezFaa2W/NbEGc9YmISIzBYGZJ4BbgfGAJcJmZLRmw2lNAq7ufBPwU+HIctZ35lYfjeBoRkQkhzj2GU4C17r7O3XPA3cDF5Su4+8Pu3hNOPg7MG61iyvue1+/uGXxFEZEqE2cwzAU2lk1vCucN5krgV5UWmNlVZtZmZm07d+58TcWYehlERCoal53PZvYBoBX4SqXl7n6bu7e6e2tLS0u8xYmITHKpGJ9rMzC/bHpeOC/CzN4FfAY4w937Rq0a7TCIiFQU5x7DCmCxmS0ysxrgUmBZ+Qpm9ibgO8BF7r4jxtpw9zifTkRk3IotGNy9AFwN3A+sBu5x91Vmdr2ZXRSu9hWgEfiJmT1tZssG2dyIy+ZLZPNFlj+7VSEhIlUtzkNJuPtyYPmAedeV3X9XnPWUe+rVdv76u08A0FCT5OnPnUM6OS67YERERlXVvvMN7GK46dcv7L/fnSvyiR8/HW9BIiLjRNUGw+LZjZHpZzbti0zft3Ir53z9kThLEhEZF6o2GBbNbDzsOi9u7+L+VdsolUa3z+Hzy1ax8Jr7WHjNffzt957g5Z1d+5f96dV2dnb2USo5uUJpVOsQEYGY+xjGk6F2MH/kzifJpBKsueH8EXvuFev3cPsf1tNYkyJXLPGzpw6ctfufL+3i7K8eek+lqTZFZ7YAwKKZDSye1cg/XXACs6ZkSCaMTCo5YrWKSPWp2mAYzL0fO5350+v52F1P8eiLwbeq+wol9nTnaK5Lk0gc2Rcg7vjDej63bNVB82+65M8598Q53PrIOm595OVDbqM/FABe2dXNK7u6eeD57ZF1jpvdyJcuOYmlx0w7onoPpydX4MkN7bxxfjNNmZSucyEyCdhEPzWztbXV29rahv247/7nOm64b/VB89d/6UIg2KNYdO3yg5b/+hNv5/g5U/ZPr9qyj289tJaObJ6vve9kZk+prfh87s7ND6zhlocPftO//uITufytCyPzCsUS+aJTm05gZrj7/p+PvrSLqXVp5k2r4+UdXTy0ZgfLn93Kxj29B2373ScdxXvfPI8zjmuJvGlvau/hih+s4KUdXVyydB73rtzC2xe38J43HsVZx88imTCWXHf//vXPO3EOHzptIS/t6OKeFRt5bss+Tl00g8fW7Y483xvmTuEdi1toXTiN+dPqeV1LI8myMO1vR65QYl9vnqbaFLXp176Hky+WeGTNTorutHfn6MkV6cwW6MkVeH1LI299/QzmTK3VGWZDUCw57T05dnX1sbsrR3N9ms5sgVTC2N2dY8PubgDcIZkwalIJDCg5ZFIJpjfUUCw5c6bWMqUuTV++RCpp1KaS9BWKvLyzi3zRqUsnWb+7m5pUgqbaFAkz9vXmWbezm/aeHEkzGmtT9OSKuMPUujQ1qQSFYolk0sjmikypSzO1Ls2U2jTJhLGtI0tfvkh9JkUqfL3li05vvkh7d46iO021KRprUsxozNBYm2JXZx/bO7J09hXoy5eor0mSK5RwHHdoqk2TThp9hRIJM0rupJOGe/C6m1pfw/T6NMlkUFux5PTmiuzuzpEwozadYEZjhhkNNftf5wbkS07Jne6+AiWH2lSCvb15OrMFcoUSNakEzXVpatNJCqUSqUSChkySQtEplErs6OyjozdPT67IuSfO4Y3zm1/T39vMnnT31orLFAwHLJ7VyIOfPGP/9HGf+RW54sHH9X/y92+ldcE08kXnwz9s45Fwz+LP507lZ//wNlIV3oRO/eJv2daRBeBtr5/Bjz586rBrHqreXJEXtnVw38qtfPe/Xtk/P5kwTjt25v49oZGy9JhmFs5oYN2ubl7c3klPrhhZvmBGPQDt3Tk6yvZ2AFIJ45RF01kwo4G29Xs4urmO17c00lyfpr0nx8IZDXTnCjy0egfZQpFFMxv504Z2GjMpzOCFbZ1DqrExk2JWU4aiO42ZFC1NGV43s5GjptZydHMdMxtrmNGYYXpDDblCiUwq+Buu29XFxj29LJrZQDIRvCl05wps78jS0pihpSnDvGn11NUkcXdyxRK7unLs7uojnUzQmEkxraGGrXt72d0dvNlOr6+h5JBIQF++RF1Nkt5ckabaFNl8iY3tPezq7KM3X2Rzey9bO7LMasrQmS2wo7OPXZ191KQSbNnbS9GDN9odnX3U1yQ5amotMxozNNQkmdGYYWb4xgSQK5ToyOZJmNGbL9LVV2D7viy7unO8tL2T7R1ZRrk77ZBq0wmaatPUJBN09RWoSycxg47ePLli8AaZLRRpqEnRnSsw8K3LjIPmAdSlk9TVJOnM5skXoyvUpIK/UV06WJ5JJ0lYMJZaZzZPT75IY02KkjvJhO1/P0gnEnRVqAGgKROs31coURjFX2gyYXzh4jfw139xzGt6vIKhgv5guP2Kt7By0z6+9uCLwIE9BoBXd/fQtmEPe3vyXH/v80Pe9u1XvIUz/2zW/un/fscKfrM6+CL3KQun8+OPnBrbIZet+3r54vIX+OUzWw5adsVpC/nshUvY2N7D8me3cdHJR3Pfyi38/KktbN3XSzqZ4PfXnEU6maC7r8Avn9nCb1Zv57JTjuHsE2aTDz8lDfzEv7Ozjz+8vItnNu6jM5vnhW2d7OjM4g7zp9eTNKMjm2fJ0VOoSyf5/dpd+0e4HeyfG6C5Ps2+3jzucPycJjqzBY6f08QJR03hlEXTmTUlw9S64BNkS2OGl3Z08cS63WzZl2Vfb55t+7Jk80VK7uzuyvHKru4R+8etSydJJ+2g4Bupbffmi9SkEtSlk8xtrqNQKrFgRgOphNHVV6ClMUNfscS2fVl2d/XRnQs+KR+ufc31aVoaMyyY0cDrZzXQ0phhSm2alqYMHdk8U2rTmAXBunhWE46TSBhegr5iEcMolILXwZ7wk/KOziwdvYXgU37JyeaLuDvHzW6iviZFZzbPUc111CQTdGSDv2djJsXsKZnD/l+USsHzl0pOZ7ZARzZPoeRMqw/2HrKF4v43/3TSSCUS1KQOfFDrKxTZ051jb0+elqYgNA/1nP3PV0mhWGJvb55SyUklEyQTRjJhNGaCI/TuTkdvgV3dwSf8XKFEyYO6ABrCvZtsvsSUuhRT69KkkgmKRWdvb47efHH/YJ8d2TyZVPAcU+vSzGqqxYwj2hNWMFTQHwwrP38Ov3h6C5/9+XNANBjKPbd5H0c317H0Cw9WXP74tWdz6o2/3T99+VsX8J43Hk06meCvbvk9AFe/81g+fe6fDbvWkdLVV2DFK3tork/zplHuexiujmzwD9ZcX0OhWGJbRxYzY3dXH53ZAotmNnB0c93+Q1EjoVhyOnrzvLK7m/buHHu6c7T35EglEuSLJfoKJeZMqWXx7Ea2d/SxtydHQyY4JDC9IU1vrsTWfb1s2ZulvSdHb65IS1OGmY01zJlaR75YojNboL0nx6ymDM31NfTkCuztyVMsOf3vN7liiabaNF3ZAomE0dKUYXZThil1aZpqU8ybVk82XySTSgyr7cWSs7cnaFeh5DTUpKjPJElYcBioNnyjUb9QdTpUMKjzGZgZ7mqfePSUQdd5w9ypQBAc33roJW5+4EUWzKjnzONauOjko5kztZZV/3wuT25o5/Lv/5EfPraBHz62Yf/j3/PGo/nUOceNbkMOozGT4p3Hzzr8imNgSm16//1UMsG8acHhp7nNdZH1RvJNLJkwpjXUMC38+49nr6UfJpmw4Bh3Y2YUKpLJTMEALAkD4d0nHT2k9a8+azFXn7X4oPkNmRTvOK6F7/ztm/nInU9Gln3z0pP1yUxEJoSqDwYDFsxo4PFrz6alaWQ+WZ174hzWf+lCdnf1cfeKjZx27EyFgohMGFUfDP3mTK18mumRmNGY4aPvPHbEtysiMpp0creIiEQoGEREJELBICIiEQoGERGJqNpgmODf6xMRGTVVGwz9dBqpiEhU1QeDiIhEKRhERCRCwSAiIhEKBhERiVAwiIhIRNUGg6PzVUVEKqnaYOink1VFRKKqPhhERCRKwSAiIhEKBhERiVAwiIhIhIJBREQiqjYYNLqqiEhlVRsM/TS4qohIVNUHg4iIRCkYREQkItZgMLPzzGyNma01s2sqLM+Y2Y/D5U+Y2cI46xMRkRiDwcySwC3A+cAS4DIzWzJgtSuBdnc/Fvg6cFNc9YmISCDOPYZTgLXuvs7dc8DdwMUD1rkYuCO8/1PgbNO1N0VEYpWK8bnmAhvLpjcBfzHYOu5eMLN9wAxgV/lKZnYVcFU42WVma15jTTMbbopuuwrMBLW5CqjN1eFI2rxgsAVxBsOIcffbgNuOdDtm1uburSNQ0oShNlcHtbk6jFab4zyUtBmYXzY9L5xXcR0zSwFTgd2xVCciIkC8wbACWGxmi8ysBrgUWDZgnWXAB8P77wUectd3lEVE4hTboaSwz+Bq4H4gCXzf3VeZ2fVAm7svA74H3Glma4E9BOExmo74cNQEpDZXB7W5OoxKm00fyEVEpJy++SwiIhEKBhERiajaYDjc8BwTiZl938x2mNlzZfOmm9mDZvZS+HNaON/M7Jthu1ea2dKyx3wwXP8lM/tgpecaD8xsvpk9bGbPm9kqM/t4OH8yt7nWzP5oZs+Ebf7ncP6icPiYteFwMjXh/EGHlzGza8P5a8zs3LFp0dCZWdLMnjKze8PpSd1mM1tvZs+a2dNm1hbOi/e17e5VdyPo/H4ZeB1QAzwDLBnruo6gPe8AlgLPlc37MnBNeP8a4Kbw/gXArwADTgWeCOdPB9aFP6eF96eNddsGae9RwNLwfhPwIsEwK5O5zQY0hvfTwBNhW+4BLg3n3wr8j/D+PwC3hvcvBX4c3l8Svt4zwKLw/yA51u07TNs/CfwIuDecntRtBtYDMwfMi/W1Xa17DEMZnmPCcPdHCc7iKlc+vMgdwF+Vzf+hBx4Hms3sKOBc4EF33+Pu7cCDwHmjX/3wuftWd/9TeL8TWE3wrfnJ3GZ3965wMh3eHDiLYPgYOLjNlYaXuRi429373P0VYC3B/8O4ZGbzgAuB74bTxiRv8yBifW1XazBUGp5j7hjVMlpmu/vW8P42YHZ4f7C2T8jfSXi44E0En6AndZvDQypPAzsI/tFfBva6eyFcpbz+yPAyQP/wMhOqzcC/AP8IlMLpGUz+NjvwgJk9acHwPxDza3tCDokhw+PubmaT7rxkM2sE/h34hLt3WNl4i5Oxze5eBE42s2bgZ8DxY1zSqDKzdwM73P1JMztzrOuJ0enuvtnMZgEPmtkL5QvjeG1X6x7DUIbnmOi2h7uUhD93hPMHa/uE+p2YWZogFP7N3f8jnD2p29zP3fcCDwNvJTh00P8Br7z+wYaXmUhtPg24yMzWExzuPQv4BpO7zbj75vDnDoIPAKcQ82u7WoNhKMNzTHTlw4t8EPhF2fzLw7MZTgX2hbuo9wPnmNm08IyHc8J540543Ph7wGp3/1rZosnc5pZwTwEzqwP+kqBv5WGC4WPg4DZXGl5mGXBpeAbPImAx8Md4WjE87n6tu89z94UE/6MPufvfMInbbGYNZtbUf5/gNfkccb+2x7oHfqxuBL35LxIcp/3MWNdzhG25C9gK5AmOJV5JcGz1t8BLwG+A6eG6RnDBpJeBZ4HWsu38HUHH3FrgirFu1yHaezrBcdiVwNPh7YJJ3uaTgKfCNj8HXBfOfx3Bm9xa4CdAJpxfG06vDZe/rmxbnwl/F2uA88e6bUNs/5kcOCtp0rY5bNsz4W1V/3tT3K9tDYkhIiIR1XooSUREBqFgEBGRCAWDiIhEKBhERCRCwSAiIhEKBpHDMLOFZuZmNq4uND9e65KJT8EgcngbCUZ0fRrAzM4M35BnxlWAmf3OzL51qLpERorGShI5DA/GKNo2Gts2s7S751/LY0ezLqlu2mOQcSX8ZPxtM/uime2y4AJEN5tZIly+3sw+XeEx3yqbXm9m15nZ7WbWaWYbzez9ZtZsZnebWVd48ZJzhljT/kM24WiuD4eLdobzbw/XMzP7RzN72cx6LbjYygcqbOcyM3vIzHqBj5jZDDO7y8w2hY9bZWZXlD3uduAM4KPh4z3c1kGHkszsHRZcpCZrZtvN7OvhsC9D+v2G6/w3Cy760mtme8zsETPrH81TqoCCQcajvwEKwNuAq4FPAO8f5jY+QTAswlKCC7vcQXCxl+XAycCjwL+aWe0wt7sRuCS8fyLBoZyPh9M3EAxH8lGCi8PcCHzHzC4csI0bgW+H6/ycYCiHPwHvDrf5jfBxZ4frfxx4DPhB+HxHER1SGQAzm0tw0ZanCIYivxK4LHy+coP+fs1sDsGAdXcAJxBcBOrOIfxeZDIZ67FBdNOt/Ab8DnhswLwHge+G99cDn67wmG+VTa8H7iqbbiQYW+mbZfMWhvNah1BTZF2CcXucsqtsAQ1AL/D2AY/9F2D5gO18agjPeXd/myu1cZC6/jfBWDqJsnU+BPQB9UP8/S4Nt7lgrF8Luo3dTX0MMh6tHDC9BZj1Wrfh7l1m1kMwyFi/7eHP4W53MEsIPvn/2qJj5acJgqpcW/mEmSUJLtf4foKLqWQILjn7u2HWcALwuLuXyub9V7itYznwOznU7/cZgkHanjOzB8L7P3X3ncOsRSYwBYOMRwM7Y50Dhz1LBCNKlksPcRv5AdMwcodT+7fzHuDVw9TSPWD608CnCA4ZPQt0AV9k5EILDrS3Uj37f7/uXgz7Xk4lGKr5SuBGMzvD3Z8ZwXpkHFMwyESzk+AYOwBhH8HxBMfV45ILfybL5j1PcMhmgbs/NMztnQ780t3vhP3XmzgO2DvgOZMVHltuNfA+M0uU7TWcHj725aEW4+5O0KfxmJldTzD88/sJ9iakCigYZKJ5CPg7M1tGEBKfIf7X8QaCT9kXmtkvgV537zSzm4Gbwzf2Rwn6Nk4FSu5+2yG29yLwfjM7HdgFfAxYRDTs1gOnhGdFdQF7Kmzn2wQdyd82s28QjO3/JYK+iZ6hNCy82Mu7CC7qsp2gE3s+QfBJldBZSTLR3EgQDr8AHiA4hh7n3gIeXHrxcwSdvduB/lNlPwt8nuDQ0CqCTt1LgFcOs8kbCM6g+hVBoHQD/zZgnZsJPvk/TxCIxwxS1/kEb+ZPA98nuIjTPw2jefsILql5L0FH9leBL7j7vw5jGzLB6UI9IiISoT0GERGJUDBI1TOzW8NvQ1e63TrW9YnETYeSpOqZ2SxgyiCLO9x9R5z1iIw1BYOIiEToUJKIiEQoGEREJELBICIiEQoGERGJ+P+emyWYwg+p8AAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Instantiate learner in VW\n",
    "vw = pyvw.vw(\"--cb_explore_adf -q UA --quiet --epsilon 0.2\")\n",
    "\n",
    "num_iterations = 5000\n",
    "ctr = run_simulation(vw, num_iterations, users, times_of_day, actions, get_cost, do_learn=False)\n",
    "\n",
    "plot_ctr(num_iterations, ctr)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Scenario 2\n",
    "\n",
    "In the real world people's preferences change over time. So now in the simulation we are going to incorporate two different cost functions, and swap over to the second one halfway through. Below is a a table of the new reward function we are going to use, `get_cost_1`:\n",
    "\n",
    "### Tom\n",
    "\n",
    "| | `get_cost` | `get_cost_new1` |\n",
    "|:---|:---:|:---:|\n",
    "| **Morning** | Politics | Politics |\n",
    "| **Afternoon** | Music | Sports |\n",
    "\n",
    "### Anna\n",
    "\n",
    "| | `get_cost` | `get_cost_new1`  |\n",
    "|:---|:---:|:---:|\n",
    "| **Morning** | Sports | Sports |\n",
    "| **Afternoon** | Politics | Sports |\n",
    "\n",
    "This reward function is still working with actions that the learner has seen previously.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "def get_cost_new1(context,action):\n",
    "    if context['user'] == \"Tom\":\n",
    "        if context['time_of_day'] == \"morning\" and action == 'politics':\n",
    "            return USER_LIKED_ARTICLE\n",
    "        elif context['time_of_day'] == \"afternoon\" and action == 'sports':\n",
    "            return USER_LIKED_ARTICLE\n",
    "        else:\n",
    "            return USER_DISLIKED_ARTICLE\n",
    "    elif context['user'] == \"Anna\":\n",
    "        if context['time_of_day'] == \"morning\" and action == 'sports':\n",
    "            return USER_LIKED_ARTICLE\n",
    "        elif context['time_of_day'] == \"afternoon\" and action == 'sports':\n",
    "            return USER_LIKED_ARTICLE\n",
    "        else:\n",
    "            return USER_DISLIKED_ARTICLE"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "To make it easy to show the effect of the cost function changing we are going to modify the `run_simulation` function. It is a little less readable now, but it supports accepting a list of cost functions and it will operate over each cost function in turn. This is perfect for what we need."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [],
   "source": [
    "def run_simulation_multiple_cost_functions(vw, num_iterations, users, times_of_day, actions, cost_functions, do_learn = True):\n",
    "    cost_sum = 0.\n",
    "    ctr = []\n",
    "\n",
    "    start_counter = 1\n",
    "    end_counter = start_counter + num_iterations\n",
    "    for cost_function in cost_functions:\n",
    "        for i in range(start_counter, end_counter):\n",
    "            # 1. in each simulation choose a user\n",
    "            user = choose_user(users)\n",
    "            # 2. choose time of day for a given user\n",
    "            time_of_day = choose_time_of_day(times_of_day)\n",
    "\n",
    "            # Construct context based on chosen user and time of day\n",
    "            context = {'user': user, 'time_of_day': time_of_day}\n",
    "\n",
    "            # 3. Use the get_action function we defined earlier\n",
    "            action, prob = get_action(vw, context, actions)\n",
    "\n",
    "            # 4. Get cost of the action we chose\n",
    "            cost = cost_function(context, action)\n",
    "            cost_sum += cost\n",
    "\n",
    "            if do_learn:\n",
    "                # 5. Inform VW of what happened so we can learn from it\n",
    "                vw_format = vw.parse(to_vw_example_format(context, actions, (action, cost, prob)),pyvw.vw.lContextualBandit)\n",
    "                # 6. Learn\n",
    "                vw.learn(vw_format)\n",
    "\n",
    "            # We negate this so that on the plot instead of minimizing cost, we are maximizing reward\n",
    "            ctr.append(-1*cost_sum/i)\n",
    "        start_counter = end_counter\n",
    "        end_counter = start_counter + num_iterations\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "def run_simulation_multiple_cost_functions(vw, num_iterations, users, times_of_day, actions, cost_functions, do_learn = True):\n",
    "    cost_sum = 0.\n",
    "    ctr = []\n",
    "\n",
    "    start_counter = 1\n",
    "    end_counter = start_counter + num_iterations\n",
    "    for cost_function in cost_functions:\n",
    "        for i in range(start_counter, end_counter):\n",
    "            # 1. in each simulation choose a user\n",
    "            user = choose_user(users)\n",
    "            # 2. choose time of day for a given user\n",
    "            time_of_day = choose_time_of_day(times_of_day)\n",
    "\n",
    "            # Construct context based on chosen user and time of day\n",
    "            context = {'user': user, 'time_of_day': time_of_day}\n",
    "\n",
    "            # 3. Use the get_action function we defined earlier\n",
    "            action, prob = get_action(vw, context, actions)\n",
    "\n",
    "            # 4. Get cost of the action we chose\n",
    "            cost = cost_function(context, action)\n",
    "            cost_sum += cost\n",
    "\n",
    "            if do_learn:\n",
    "                # 5. Inform VW of what happened so we can learn from it\n",
    "                vw_format = vw.parse(to_vw_example_format(context, actions, (action, cost, prob)),pyvw.vw.lContextualBandit)\n",
    "                # 6. Learn\n",
    "                vw.learn(vw_format)\n",
    "\n",
    "            # We negate this so that on the plot instead of minimizing cost, we are maximizing reward\n",
    "            ctr.append(-1*cost_sum/i)\n",
    "        start_counter = end_counter\n",
    "        end_counter = start_counter + num_iterations\n",
    "\n",
    "    return ctr"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### With learning\n",
    "Let us now switch to the second reward function after a few samples (running the first reward function). Recall that this reward function changes the preferences of the web users but it is still working with the same action space as before. We should see the learner pick up these changes and optimize towards the new preferences.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAEOCAYAAACNY7BQAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO3deZRdZZnv8e9zxkoNqSQkgYwQmiBEFIgRoRkVWyaFvtpXQF044MV2YOl1utCOjbbgPDQqRptGUUFFL0YaGpRBl1emMEoSAyFzQuahKjWe4bl/7F2VsyunkqqkzlC1f5+1auXs4ex69tmV8zv7ffd5t7k7IiIifRK1LkBEROqLgkFERCIUDCIiEqFgEBGRCAWDiIhEKBhERCSiasFgZjeb2RYze26Q5WZm3zGzFWb2rJnNr1ZtIiKyVzXPGG4Bzt/P8guAueHPVcD3q1CTiIgMULVgcPc/ATv2s8olwE888AgwwcymVac6ERHpk6p1ASVmAOtKpteH814auKKZXUVwVkFTU9OrjjvuuKoUKCIyVjzxxBPb3H1KuWX1FAxD5u4LgYUACxYs8MWLF9e4IhGR0cXM1gy2rJ6uStoAzCqZnhnOExGRKqqnYFgEXBFenXQqsNvd92lGEhGRyqpaU5KZ3QacA0w2s/XA54A0gLvfBNwNXAisADqBd1erNhER2atqweDulx9guQMfrFI5IiIyiHpqShIRkTqgYBARkQgFg4iIRCgYREQkQsEgIiIRCgYREYlQMIiISISCQUREIhQMIiISoWAQEZEIBYOIiEQoGEREJELBICIiEQoGERGJUDCIiEiEgkFERCIUDCIiEqFgEBGRCAWDiIhEKBhERCRCwSAiIhEKBhERiVAwiIhIhIJBREQiFAwiIhKhYBARkQgFg4iIRCgYREQkQsEgIiIRCgYREYlQMIiISISCQUREIhQMIiISoWAQEZEIBYOIiERUNRjM7HwzW25mK8zsmjLLZ5vZg2b2lJk9a2YXVrM+ERGBVLV+kZklge8C/wCsBx43s0XuvrRktU8Dv3T375vZPOBu4Khq1SgyXO5Oe0+etds76ejJs7m9h67ePJlUgiMPa+LkWRMws1qXKTIsVQsG4BRghbuvBDCz24FLgNJgcGB8+LgV2FjF+mSUKxadtTs62dOTpydfYNW2TgzIFYqs3dHJ1vYemrIpsukEubyzvaOHnlyRhnSCiU0ZUgljXDpJNp2kszfPjo5eUokEU1qybG7rpuhOe3eeQtFp687Rkyvy/OZ22rrzg9b0ypmt3Hj5fGYf1li9F0LkEFUzGGYA60qm1wOvGbDO54H7zOxqoAl4fbkNmdlVwFUAs2fPHvFCpTbcnafX7eKZdbvIF53xDWleMbOVmRPH0ZxNUSg6yYSRKzjLN7Xz3MbdrNiyh12dOdZs72DZS2109BYG3f6kpgx7eoI3dgMmN2fJphN09RbY1ZUjVyjivnf9loYU7eGbfjaVIJtK0JRNMS6TZHxDmkwqwUWvnMbsSU3MnDiO1nFpjmhtIJ1MkC8UeXjldr509zLe8K0/8k+vmsnVr5vL4eMbKvwqykgrFp2uXPB31ZsvkisW6ewpkEoamWSC3kKRhBnduQKdvQUSZuQKRfJFJ18o4uHzOnvz5AqOA/lCkc7eAr35Ir2FIu3dOZJmEG6nO1fo/z8wLp0klTTyBafoTmMmSTqZoCtX4Oxjp3DCjNYR3+dqBsNQXA7c4u5fN7PTgFvN7AR3L5au5O4LgYUACxYs8DLbkTpQLDp/29ROV67Auh2d9OaLrNi6h7auHLMmNZJNJXhk5Y7+/xhPr9sVeWMuJ5UwzCBXCFZMJoyJjWmmTxjHxSdNZ9608bQ2ZhiXTjKttYF80WnOppg5cRwN6WT/dtx9nyaeYtEpuLOnO8+ExjRmRiF8U2jKJIfdJDT38BbOnDuFr9+3nNseW8fvnnmJj5/3Mv7xpOm0NKSHtS05NN25Alvbe9jR0cuenjzb9vSwuytHd67Ajo4c2/b00NaVo+jQ0ZNnV1eOYtHZtqeH7R29Fa8vYVAM//azqQQN6SQJgz09+f6/9XJaGlIVCQbzA/1PHKlfFLzRf97dzwunrwVw9+tL1lkCnO/u68LplcCp7r5lsO0uWLDAFy9eXNHa5cA2t3Xz+6WbuevZjTy5dhe9+WLZ9dJJi/yhN2WStI5L09qY4bgjWjj16EmcfsxkcgWnrSvHsxt2s6Wtm217esiEn5KasilOnj2RedPGM3tSI5lU/V9ct2JLO//n13/liTU7SSaME2e28ndTmjnr2CkcP20801obaMrW2+e0oXF3esLj7Q5duUL/G3BjJsnMiePozhXJF4t09RbYtqeH7lyR1nFBOPbkg8fduQK7u3Ls6cmztb2HdNLY2t5DV65AruDs7syxcXcXW9p6yBeLpJMJevJFdnflSJgxoTHdH+7JhNGTK9LWlaMzV6BQHPx9Lp00Jjdn++tpzqaY0JgmYcZhzVmmNGdoyqYouNOQSpJOJWhMJykUnd5CkXTScIdsOkFjJoW7k0okSCaNdCL42wyWJckkE5hBMhFMZ1MJMqkE48IPLe6QSEQ/gBSKTm++SCppJM3Y05unWHQa0snIh53hMrMn3H1B2WVVDIYU8DxwLrABeBx4m7svKVnnHuAX7n6LmR0P3A/M8P0UqWCojpd2d/HnF7YxsTHD+p2d/N+nNzKxMU1PLmgyGWhaawO7OnO87+yjmT2pkaZsiiPGN3D8tPEUwr6AhMHRU5pJJuLROevuPLpqB/cu2cTDL25n466uSP/EiTNbecepR/Lm+TOr+pq4O7mC09mbZ832Tra099CbL1JwZ3dXDgPW7uhkS1s3W9p7+vtsunoLuHPAN95DlTBIJRO0jgvODKe2ZEknjd6805AO5ucKRdq78yQTRiI808umE4xvSNOUTdKUTTGpMcNhzVmasykmN2dobQyaaZoyqX3ejOOgLoIhLORC4FtAErjZ3f/NzK4DFrv7ovBKpB8CzQQd0Z909/v2t00Fw8jqzRdZ+lIbP390DZ29BRav3smmtu79PieVME6cNYFLTprOzInjeO3LpupKnCEoFoOgWLezkzXbO/jt0xtZv7OLedPGc+KsCcybPp7Tjj6MY6Y2H9T23Z2dnTme39zOmu0dtHfnWbejkw27utnVGXyi37iri/ae/AGb8DLJBIc1Z5g6voFsKsGRYdhD8Gm4JZvCHZLJoAN/UlOmv49m/c4uGjNJEmY0ZpJMbsmSTSX6+2+SZnT05smmgjfyxvBDRL5YDPqBwjNC/U2NrLoJhkpQMBy6QtHZuKuLO55Yz7fvf6HsOnOnNjNnchNTx2eZP3siJ8+eSK5QZFprg9rLR0ix6Nz2+Fp+/uhaVm3roDPsSD/qsEamtGSZPmEcJ0xvZVwmycqtHWxpD66USiYSGEE7ene+SG++wJb2Hjbt7u7fRp+mTJJpE8YxuTlDUybFjLBjP5NK9PfFTGnJ0pAOmj1aGtIU3ZnakiWVrP8mOxk6BYPQ1Vvga/ct5y8vbmf5pjaKDsce3syuzhxb2nv612sdl+adpx3JtAnjeMWMVl4+fbw+qdVAoeis2raHe5dsZunGNlZu62DT7i52duYAaEgnmNrSQDJh5IvBVTHj0kmSCSOTSjCttYGpLQ3MmtTI0ZObOGZqMy0NKVrHpXU8Bdh/MIzO3i45oGLRuW/pZn7z5HruW7q57DrPb94DwJSWLHMmN/H+s/+O1x43tZplyiCSCeOYqS0cM7Wlf567s21PL4WiM7k5o0/wUjEKhjHmqbU7+c79L/CXF7f3XynS5x9Pms57zzyarlyBLW09vLS7izfPn8mkpkyNqpXhMDOmtGRrXYbEgIKhju3uzPHwyu28/vippJIJ2rpz/OQvq1m8ZidLN7Zx9blzmTethUwyyZtu/HPkWmiAmRPHcdmrZ3HJSTOYNUnfvBWRoVEw1Al357bH1nHDPcu46+ozOe9bf+r/tuVgPnPnc5HpvlB452lH8oHXHqNv2YrIQVEw1FCh6Fy28GGeWhsMAdHnrK8+uN/n/Y+TZ7Bpd3fk+wONmSTnv/wIvvo/T4zN9wJEpDIUDFX26MrtXLrwkQOud+UZc/j0RcdjZjy2agf3L9vMW141k2MPbzngc0VEDoWCoYrau3NlQyGVMB79l3M5rLl8x+IpcyZxypxJlS5PRARQMFSFu/O+W5/ov2z0e2+fzxGtDSzd2MbbTpkdy6/ji0j9UjBUwQ/+tLI/FF4+fTwXvmIaAPNnT6xlWSIiZSkYKuy+JZu44Z6/AXDnB0/npFkTalyRiMj+KRgq5NdPrOdjv3qmf/rbl52kUBCRUUHfqa+Ae/76UiQUXnXkRN74yuk1rEhEZOh0xjDC2rtzvP9nT/ZP/+y9r+H0YybXsCIRkeFRMIywV3w+uH3EibMm8NsPnl7jakREhk9NSSMoX9g7aN2dH/j7GlYiInLwFAwj6F9/txSAG992ssa8F5FRS8EwQvKFIrc+sgaA1+meBiIyiikYRsifV2wD4BPnvYzGjLpuRGT00jvYIejNFzn20/cAkE4GTUfvPXNOLUsSETlkOmM4BJ9btKT/ca4QDJudTSVrVY6IyIhQMBykzW3d3PbY2si8ez9yVo2qEREZOWpKOkiv+dL9/Y9X33AR7q4rkURkTNAZw0HoLrnl5q1XngKgUBCRMUPBcBDe8aNHAfjhFQs4c+6UGlcjIjKyFAzD5O4sXrMTCAbHExEZaxQMw7R+ZxcAU1uyTGrK1LgaEZGRp2AYpoeWbwHg5ne9usaViIhUhoJhmD7z2+C7C8cd0VLjSkREKkPBcBBaGlKkknrpRGRs0rvbMFxx82MAtHfna1yJiEjlKBiGaN2OTv70/FYAfvyeU2pcjYhI5SgYhmjJxt39j88+Vt9dEJGxS8EwRP/80+A+zj/R2YKIjHFVDQYzO9/MlpvZCjO7ZpB13mpmS81siZn9vJr1DcWZcyfXugQRkYqqWjCYWRL4LnABMA+43MzmDVhnLnAtcLq7vxz4SKXr6s0X+fyiJezo6B10nT8s3VxaY6VLEhGpqWqeMZwCrHD3le7eC9wOXDJgnf8FfNfddwK4+5ZKF/XfSzZxy19W86W7l5Vd7u689yeLK12GiEjdqGYwzADWlUyvD+eVOhY41sz+n5k9Ymbnl9uQmV1lZovNbPHWrVsPqah8oQhAoehllz+2akf/4799oWw5IiJjSr11PqeAucA5wOXAD81swsCV3H2huy9w9wVTphzaFUIe5sFgDUR3Pr2h/3FDWndnE5Gxr5rBsAGYVTI9M5xXaj2wyN1z7r4KeJ4gKCrmnudeAqCtzJfWnly7k9seC05yHrn23EqWISJSN6oZDI8Dc81sjpllgMuARQPWuZPgbAEzm0zQtLSykkX9YVnQjfHi1j38df3uSJPSTQ+92P/4iNaGSpYhIlI3qhYM7p4HPgTcCywDfunuS8zsOjO7OFztXmC7mS0FHgQ+4e7bq1Hfqm0dvOnGP3Pydff1z+u7AEkXIolInFT1ns/ufjdw94B5ny157MBHw5+quPjE6Sx6ZmP/dGmT0lNrdwGw8ksXVqscEZGaq7fO56qb2JgedNmW9h5A310QkXip6hlDPfrxw2sGXdaSTTFft+8UkZiJ/RnDQH0D5HX1FmjvyTNv+vgaVyQiUl0KhgEefjHo6/7w7U8BsCG8x7OISFwoGAboLRRp785xXzg+0vknHFHjikREqkvBUMaqbR39jy98xbQaViIiUn2x73wuJ5VIMKUly9GTm2pdiohI1SkYymjvzrG1vYet4eWqIiJxEuumpCfX7iw7/9KFjwC6haeIxNOQgsHM0ma2ycxeXumCqml3Vy4y/fbXzI5MHzFe4yOJSPwMqSnJ3XNmlgPK37RglEon9ubi6hsuYmdHLz97dG3/vDedOL0WZYmI1NRwmpL+HbjWzMZMv0Q2Hex+NhX8O7Epwwkz9n6h7djDm2tSl4hILQ3nTf5M4Gxgg5k9B3SULnT3i8s+q471jYD0jbee1D/vghOm8dyGNiAIChGRuBnOGcM24NcEo6OuBbYP+Bl1+u69MKFkIL03z997t9F0MtZ98yISU0M+Y3D3d1eykFoohPf1TJSMnjq5OVurckRE6sKQPxKb2QPl7r9sZuPN7IGRLas6lm4MmoyKvrdPve8s4YxjJtekJhGRWhtOH8M5QLlG9waC/odR54v/tQyA5ze3c3pJEKy+4aJalSQiUnMHDAYzm18y+Uoz21EynQTOAzaMdGHVlEmpL0FEpM9QzhgWE3x/wYH7yizvAq4eyaKqLaNOZhGRfkMJhjkEV3auBF5NcHVSn15gi7sXKlBb1eiMQURkrwMGg7uvATCzG4BXu/tNpcvN7J/NbIa7f6ZCNVZcMqF7OouI9BnOR+V3AE+Vmf8kcMXIlFMbukRVRGSv4QTDVGBrmfnbgMNHppzaOPXow2pdgohI3RhOMKyl/GWpZwHrR6YcERGpteF8j+EHwDfNLAP0faHtXOB64MsjXZiIiNTGcIbE+LqZTQa+w94vuvUC33b3r1SiOBERqb5hDaHt7tea2ReBeeGsZe6+Z+TLEhGRWhn2vRXcvQN4vAK1iIhIHdA3u0REJGLM3I3tYLxiRitTWvQdBhGRUrE+Yyi6oy89i4hExTwYwEzJICJSKtbB4DpjEBHZR6yDIWhKUjKIiJSKeTCgYBARGaCqwWBm55vZcjNbYWbX7Ge9t5iZm9mCStZTdEe5ICISVbVgMLMk8F3gAoJvTl9uZvPKrNcCfBh4tNI1uc4YRET2Uc0zhlOAFe6+0t17gduBS8qs9wWCQfm6K12QLlcVEdlXNYNhBrCuZHp9OK+fmc0HZrn7f+1vQ2Z2lZktNrPFW7eWu0XE0KjzWURkX3XT+WxmCeAbwMcOtK67L3T3Be6+YMqUKQf9O4tFfY9BRGSgagbDBmBWyfTMcF6fFuAE4CEzWw2cCiyqZAe0vscgIrKvagbD48BcM5sT3uznMmBR30J33+3uk939KHc/CngEuNjdF1eqIF2uKiKyr6oFg7vngQ8B9wLLgF+6+xIzu87MLq5WHaV0uaqIyL6qOrqqu98N3D1g3mcHWfecStejsZJERPZVN53PtaEzBhGRgWIdDEWHpJJBRCQi5sGgMwYRkYFiHQwaEkNEZF+xDoaie61LEBGpO7EOBnTGICKyj1gHg/oYRET2FetgcNCQGCIiA8Q6GIIzBiWDiEipWAeDO2pKEhEZIPbBoM5nEZGoWAdD0R3FgohIVKyDIeh8VjSIiJSKdTDoclURkX3FOhhcw26LiOwjtsHg4XAYigURkagYB0Pwr/oYRESiYhsMfQPoKRdERKJiGwx946pqSAwRkajYBsPeMwYlg4hIqdgGQ18fg3JBRCRKwaDrkkREIuIbDGEvg/oYRESiYhsMRV2uKiJSVoyDQZerioiUE9tg6Otj6MkXa1uIiEidiW0wLNm4G4Cv3ru8xpWIiNSX2AZDW1e+1iWIiNSl2AaDiIiUp2AQEZGI2AbD5OYMAIePz9a4EhGR+hLjYAgC4RPnHVfjSkRE6ktsg6FPMvavgIhIVGzfFv3Aq4iIxFJVg8HMzjez5Wa2wsyuKbP8o2a21MyeNbP7zezIStWy99ae+uqziEipqgWDmSWB7wIXAPOAy81s3oDVngIWuPsrgTuAr1S+rkr/BhGR0aWaZwynACvcfaW79wK3A5eUruDuD7p7Zzj5CDCzUsWoKUlEpLxqBsMMYF3J9Ppw3mCuBO4pt8DMrjKzxWa2eOvWrQdVjCsZRETKqsvOZzN7B7AA+Gq55e6+0N0XuPuCKVOmHORv0a09RUTKSVXxd20AZpVMzwznRZjZ64FPAWe7e0+li1IsiIhEVfOM4XFgrpnNMbMMcBmwqHQFMzsZ+AFwsbtvqWQxakoSESmvasHg7nngQ8C9wDLgl+6+xMyuM7OLw9W+CjQDvzKzp81s0SCbO/R6wn/VkiQiElXNpiTc/W7g7gHzPlvy+PXVrAf0PQYRkYHqsvO5GtSUJCJSXnyDAd3zWUSknPgGQ3jGoFwQEYmKbTD00RmDiEhUbINBfQwiIuXFNxj2XrBa0zpEROpNfIOhr49BuSAiEhHbYOijXBARiYp9MIiISFRsg2FvU5LOGURESsU2GPooFkREomIbDK57uImIlBXfYNBVSSIiZcU3GMJ/FQwiIlGxDYY+GnZbRCQqtsHgGhNDRKSs+AZD3wOdMIiIRMQ2GPooF0REomIbDGpJEhEpL7bBQP8d3HTOICJSKrbBoDu4iYiUF9tg6KMTBhGRqNgGg7oYRETKi28w9Dcl6ZRBRKRUbIOhj5qSRESiYhsM+uaziEh58Q2G8F+dMIiIRMU3GJQMIiJlxTYY+qjzWUQkKrbBoDu4iYiUF9tgQHdwExEpK77BEFIuiIhExTYY1JAkIlJefIOhvylJ5wwiIqXiGwz9w27XuBARkTpT1WAws/PNbLmZrTCza8osz5rZL8Llj5rZURWvqdK/QERklKlaMJhZEvgucAEwD7jczOYNWO1KYKe7HwN8E/hyperp6ClUatMiIqNaNc8YTgFWuPtKd+8FbgcuGbDOJcCPw8d3AOdahToB/vj8FgAyqdi2pomIlJWq4u+aAawrmV4PvGawddw9b2a7gcOAbaUrmdlVwFXh5B4zW36QNU0+8cvRbcfAZNA+x4D2OR4OZZ+PHGxBNYNhxLj7QmDhoW7HzBa7+4IRKGnU0D7Hg/Y5Hiq1z9VsR9kAzCqZnhnOK7uOmaWAVmB7VaoTERGgusHwODDXzOaYWQa4DFg0YJ1FwDvDx/8EPOC6cYKISFVVrSkp7DP4EHAvkARudvclZnYdsNjdFwH/AdxqZiuAHQThUUmH3Bw1Cmmf40H7HA8V2WfTB3IRESmlazVFRCRCwSAiIhGxDYYDDc8xWpjZLDN70MyWmtkSM/twOH+Smf3ezF4I/50Yzjcz+06438+a2fySbb0zXP8FM3vnYL+zXphZ0syeMrO7wuk54VAqK8KhVTLh/EGHWjGza8P5y83svNrsydCY2QQzu8PM/mZmy8zstLF+nM3sf4d/18+Z2W1m1jDWjrOZ3WxmW8zsuZJ5I3ZczexVZvbX8DnfGdKXht09dj8End8vAkcDGeAZYF6t6zrIfZkGzA8ftwDPEww58hXgmnD+NcCXw8cXAvcQDBN1KvBoOH8SsDL8d2L4eGKt9+8A+/5R4OfAXeH0L4HLwsc3Ae8PH38AuCl8fBnwi/DxvPDYZ4E54d9Estb7tZ/9/THw3vBxBpgwlo8zwRdeVwHjSo7vu8bacQbOAuYDz5XMG7HjCjwWrmvhcy84YE21flFqdCBOA+4tmb4WuLbWdY3Qvv0W+AdgOTAtnDcNWB4+/gFwecn6y8PllwM/KJkfWa/efgi+B3M/8DrgrvCPfhuQGniMCa6EOy18nArXs4HHvXS9evsh+E7PKsILRgYev7F4nNk7EsKk8LjdBZw3Fo8zcNSAYBiR4xou+1vJ/Mh6g/3EtSmp3PAcM2pUy4gJT51PBh4FDnf3l8JFm4DDw8eD7ftoe02+BXwSKIbThwG73D0fTpfWHxlqBegbamU07fMcYCvwn2Hz2Y/MrIkxfJzdfQPwNWAt8BLBcXuCsX2c+4zUcZ0RPh44f7/iGgxjjpk1A78GPuLubaXLPPioMGauSzazNwJb3P2JWtdSRSmC5obvu/vJQAdBE0O/MXicJxIMrDkHmA40AefXtKgaqMVxjWswDGV4jlHDzNIEofAzd/9NOHuzmU0Ll08DtoTzB9v30fSanA5cbGarCUbpfR3wbWCCBUOpQLT+wYZaGU37vB5Y7+6PhtN3EATFWD7OrwdWuftWd88BvyE49mP5OPcZqeO6IXw8cP5+xTUYhjI8x6gQXmHwH8Ayd/9GyaLS4UXeSdD30Df/ivDqhlOB3eEp673AG8xsYvhJ7Q3hvLrj7te6+0x3P4rg2D3g7m8HHiQYSgX23edyQ60sAi4Lr2aZA8wl6KirO+6+CVhnZi8LZ50LLGUMH2eCJqRTzawx/Dvv2+cxe5xLjMhxDZe1mdmp4Wt4Rcm2BlfrTpcadvZcSHAFz4vAp2pdzyHsxxkEp5nPAk+HPxcStK3eD7wA/AGYFK5vBDdMehH4K7CgZFvvAVaEP++u9b4Ncf/PYe9VSUcT/IdfAfwKyIbzG8LpFeHyo0ue/6nwtVjOEK7WqPG+ngQsDo/1nQRXn4zp4wz8K/A34DngVoIri8bUcQZuI+hDyRGcGV45kscVWBC+fi8CNzLgAoZyPxoSQ0REIuLalCQiIoNQMIiISISCQUREIhQMIiISoWAQEZEIBYPIAZjZUWbmZlZXN5qv17pk9FMwiBzYOoLByJ4GMLNzwjfkydUqwMweMrMb91eXyEip2j2fRUYrdy8QDGQ24sws7cFwD8NWybok3nTGIHUl/GT8PTP7kpltC29g8jUzS4TLV5vZx8s858aS6dVm9lkzu8XM2s1snZldasGNbm43sz3hzUzeMMSa+ptswhFsHwwXbQ3n3xKuZ2b2STN70cy6wpujvKPMdi43swfMrAt4n5kdZsFNaNaHz1tiZu8ued4twNnAB8Pne7itfZqSzOwsC25S021mm83sm+GwL0N6fcN13mzBTWC6zGyHmf3RzPpG95QYUDBIPXo7kAf+HvgQ8BHg0mFu4yMEwyLMJ7ixy48JbupzN8HQEn8CfmpmDcPc7jrgLeHjlxM05Xw4nP4iwXAGHyS4Ocz1wA/M7KIB27ge+F64zp0EQzk8Cbwx3Oa3w+edG67/YeBh4D/D3zeN6BDLAJjZDIIbsTxFMPz6lQTj718/YNVBX18zO4JgYMIfA8cT3ETm1iG8LjKW1HqcEP3op/QHeAh4eMC83wM/Ch+vBj5e5jk3lkyvBm4rmW4mGE/qOyXzjgrnLRhCTZF1CcZncmByyTpNQBdw5oDnfgu4e8B2PjaE33l73z6X28dB6vo3grF1EiXrvAvoARqH+PrOD7d5ZK3/FvRTux/1MUg9enbA9EZg6sFuw933mFknwaBjfTaH/w53u4OZR/DJ/7/NrHQAsjRBUJVaXDphZkmCeytcSnATlRNT7iYAAAIRSURBVCzBrTsfGmYNxwOPuHuxZN6fw20dw97XZH+v7zMEg7Y9Z2b3hY/vcPetw6xFRjEFg9SjgZ2xzt5mzyLBCJOl0kPcRm7ANIxcc2rfdt5EMFz0/mrpGDD9ceBjBE1GfwX2AF9i5EILojd6GfT1dfdC2PdyKsHQzVcC15vZ2e7+zAjWI3VMwSCjzVaCNnYAwj6C4wja1aulN/w3WTJvKUGTzZHu/sAwt3cG8Dt3vxX677FxLLBrwO9MlnluqWXAW80sUXLWcEb43BeHWoy7O0GfxsNmdh2whOBsRsEQEwoGGW0eAN5jZosIQuJTVP/veA3Bp+yLzOx3QJe7t5vZ14CvhW/sfyLo2zgVKLr7wv1s73ngUjM7g+AG9lcT3M6yNOxWA6eEV0XtAXaU2c73CDqSv2dm3ya4b8ENBH0TnUPZsfDmL68nuPHLZoJO7FkEwScxoauSZLS5niAcfgvcR9CGXs2zBTy4Sf3nCDp7NxPc/ATgM8DnCZqGlhB06r4FWHWATX6R4AqqewgCpQP42YB1vkbwyX8pQSDOHqSuCwjezJ8Gbia4Ccy/DGP3dhPcPvMugo7srwNfcPefDmMbMsrpRj0iIhKhMwYREYlQMEjsmdlN4behy/3cVOv6RKpNTUkSe2Y2FRg/yOI2d99SzXpEak3BICIiEWpKEhGRCAWDiIhEKBhERCRCwSAiIhH/HzsKTKrfrPJVAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# use first reward function initially and then switch to second reward function\n",
    "\n",
    "# Instantiate learner in VW\n",
    "vw = pyvw.vw(\"--cb_explore_adf -q UA --quiet --epsilon 0.2\")\n",
    "\n",
    "num_iterations_per_cost_func = 5000\n",
    "cost_functions = [get_cost, get_cost_new1]\n",
    "total_iterations = num_iterations_per_cost_func * len(cost_functions)\n",
    "\n",
    "ctr = run_simulation_multiple_cost_functions(vw, num_iterations_per_cost_func, users, times_of_day, actions, cost_functions)\n",
    "\n",
    "plot_ctr(total_iterations, ctr)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "**Note:** The initial spike in CTR depends on the rewards received for the first few examples. When you run on your own, you may see something different initially because our simulator is designed to have randomness.\n",
    "\n",
    "### Without learning"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAEOCAYAAACNY7BQAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAd00lEQVR4nO3debhcdZ3n8fentruFbCRISKKJgktcWvE2Im60ogIqPNPOKDzN49K0ODPioyNqg7a2o3ajtt0uIzYwSuPQKrJ0a0QUW8DxcUTkIhgIEAlrboAsZL17Ld/545yEOpe6yb1wb92lPq/nqSdn+dWp76lzU586S/2OIgIzM7N9ctNdgJmZzSwOBjMzy3AwmJlZhoPBzMwyHAxmZpbhYDAzs4ymBYOkSyRtlXTnGPMl6euSNkpaJ+noZtVmZmZPaOYew6XAiQeYfxJwVPo4C/jnJtRkZmajNC0YIuJXwI4DNDkV+D+R+C2wUNKy5lRnZmb7FKa7gDrLgU11473ptEdHN5R0FsleBV1dXS9//vOf35QCzczmiltvvXV7RCxtNG8mBcO4RcTFwMUA3d3d0dPTM80VmZnNLpIeGmveTLoqaTOwsm58RTrNzMyaaCYFw1rgXenVSccCuyPiSYeRzMxsajXtUJKk7wPHA0sk9QJ/CxQBIuJC4FrgZGAjMAC8t1m1mZnZE5oWDBFx+kHmB/CBJpVjZmZjmEmHkszMbAZwMJiZWYaDwczMMhwMZmaW4WAwM7MMB4OZmWU4GMzMLMPBYGZmGQ4GMzPLcDCYmVmGg8HMzDIcDGZmluFgMDOzDAeDmZllOBjMzCzDwWBmZhkOBjMzy3AwmJlZhoPBzMwyHAxmZpbhYDAzswwHg5mZZTgYzMwsw8FgZmYZDgYzM8twMJiZWYaDwczMMhwMZmaW4WAwM7MMB4OZmWU4GMzMLMPBYGZmGQ4GMzPLcDCYmVmGg8HMzDKaGgySTpS0QdJGSec2mP9MSTdKuk3SOkknN7M+MzNrYjBIygMXACcBa4DTJa0Z1exvgCsi4mXAacA3m1WfmZklmrnHcAywMSLuj4gR4HLg1FFtApifDi8AHmlifWZmRnODYTmwqW68N51W7zPAGZJ6gWuBDzZakKSzJPVI6tm2bdtU1Gpm1rJm2snn04FLI2IFcDJwmaQn1RgRF0dEd0R0L126tOlFmpnNZc0Mhs3AyrrxFem0emcCVwBExE1AO7CkKdWZmRnQ3GC4BThK0mpJJZKTy2tHtXkYeAOApBeQBIOPFZmZNVHTgiEiKsDZwHXA3SRXH62X9FlJp6TNzgHeJ+kPwPeB90RENKtGMzODQjNfLCKuJTmpXD/t03XDdwGvamZNZmaWNdNOPpuZ2TRzMJiZWYaDwczMMhwMZmaW4WAwM7MMB4OZmWU4GMzMLMPBYGZmGQ4GMzPLcDCYmVmGg8HMzDIcDGZmluFgMDOzDAeDmZllOBjMzCzDwWBmZhkOBjMzy3AwmJlZhoPBzMwyHAxmZpbhYDAzswwHg5mZZTgYzMwsw8FgZmYZDgYzM8twMJiZWYaDwczMMhwMZmaW4WAwM7MMB4OZmWU4GMzMLMPBYGZmGQ4GMzPLcDCYmVlGU4NB0omSNkjaKOncMdq8Q9JdktZL+l4z6zMzMyg064Uk5YELgDcCvcAtktZGxF11bY4CzgNeFRE7JR3WrPrMzCzRzD2GY4CNEXF/RIwAlwOnjmrzPuCCiNgJEBFbm1ifmZnR3GBYDmyqG+9Np9V7LvBcSf9P0m8lndhoQZLOktQjqWfbtm1TVK6ZWWuaaSefC8BRwPHA6cD/lrRwdKOIuDgiuiOie+nSpU0u0cxsbmtmMGwGVtaNr0in1esF1kZEOSIeAP5IEhRmZtYkzQyGW4CjJK2WVAJOA9aOavNDkr0FJC0hObR0fxNrNDNreU0LhoioAGcD1wF3A1dExHpJn5V0StrsOuBxSXcBNwIfi4jHm1WjmZmBImK6a3hauru7o6enZ7rLMDObVSTdGhHdjebNtJPPZmY2zRwMZmaW4WAwM7MMB4OZmWU4GMzMLMPBYGZmGQ4GMzPLcDCYmVmGg8HMzDLGFQySipIek/TCqS7IzMym17iCISLKQBmY3f1nmJnZQU3kUNL/As6T1LTbgZqZWfNN5EP+NcDrgM2S7gT662dGxCkNn2VmZrPKRIJhO3D1VBViZmYzw7iDISLeO5WFmJnZzDDucwySbmh0/2VJ8yXdMLllmZnZdJnIyefjgVKD6e0k5x/MzGwOOOihJElH142+RNKOuvE88GZg82QXZmZm02M85xh6SH6/EMDPG8wfBD44mUWZmdn0GU8wrAYE3A/8KcnVSfuMAFsjojoFtZmZ2TQ4aDBExEMAkr4A/GlEXFg/X9J/lbQ8Ij41RTWamVkTTeTk8xnAbQ2m/x541+SUY2Zm020iwXAYsK3B9O3AMyanHDMzm24TCYaHaXxZ6muB3skpx8zMpttEusS4CPiKpBKw7wdtbwDOB7442YWZmdn0mEiXGP8oaQnwdZ74odsI8LWI+NJUFGdmZs03oS60I+I8SZ8H1qST7o6Ivskvy8zMpsuE760QEf3ALVNQi5mZzQC+57OZmWU4GMzMLMPBYGZmGQ4GMzPLcDCYmVmGg8HMzDIcDGZmltHUYJB0oqQNkjZKOvcA7d4uKSR1N7M+MzNrYjBIygMXACeR/HL6dElrGrQ7BPgQcHOzajMzsyc0c4/hGGBjRNwfESPA5cCpDdp9jqRTvqEm1mZmZqlmBsNyYFPdeG86bT9JRwMrI+InB1qQpLMk9Ujq2bat0S0izMzsqZoxJ58l5YB/As45WNuIuDgiuiOie+nSpVNfnJlZC2lmMGwGVtaNr0in7XMI8CLgl5IeBI4F1voEtJlZczUzGG4BjpK0Or3Zz2nA2n0zI2J3RCyJiFURsQr4LXBKRPQ0sUYzs5bXtGCIiApwNnAdcDdwRUSsl/RZSac0qw4zMzuwCd+P4emIiGuBa0dN+/QYbY9vRk1mZpY1Y04+m5nZzOBgMDOzDAeDmZllOBjMzCzDwWBmZhkOBjMzy3AwmJlZhoPBzMwyHAxmZpbhYDAzswwHg5mZZTgYzMwsw8FgZmYZDgYzM8twMJiZWYaDwczMMhwMZmaW4WAwM7MMB4OZmWU4GMzMLMPBYGZmGQ4GMzPLcDCYmVmGg8HMzDIcDGZmluFgMDOzDAeDmZllOBjMzCzDwWBmZhktHQyVao1zrvgDj+wanO5SzMxmjJYOhuvWb+Hq3/dy3BdumO5SzMxmjJYOhiCmuwQzsxmnpYOhkGvp1Tcza6ilPxlLBU13CWZmM05Tg0HSiZI2SNoo6dwG8z8i6S5J6yRdL+lZU1mP9xjMzJ6saZ+MkvLABcBJwBrgdElrRjW7DeiOiJcAVwFfmsqainkHg5nZaM38ZDwG2BgR90fECHA5cGp9g4i4MSIG0tHfAiumqpi+4QqX/uaBqVq8mdms1cxgWA5sqhvvTaeN5Uzgp41mSDpLUo+knm3btj2lYj5/zV1ct37L/vGRSu0pLcfMbK6ZkcdSJJ0BdAP/0Gh+RFwcEd0R0b106dKn9Bq7B8uZcf/Izcws0cxg2AysrBtfkU7LkHQC8EnglIgYnqpicspekXTGt2/mspseZOveoal6STOzWaGZwXALcJSk1ZJKwGnA2voGkl4GXEQSCluntJpRV6r27hzkUz9azzF/d/2UvqyZ2UzXtGCIiApwNnAdcDdwRUSsl/RZSaekzf4BmAdcKel2SWvHWNzTNnqPod7WPd5rMLPWVWjmi0XEtcC1o6Z9um74hGbVcqCftm3YspfD5rc3qxQzsxllRp58boYD7DDwyK5BfnT7Zv76qnXUasHASIUb7tky9hPMzOaQpu4xzCSNcmFeW4G+4Qp/ffUd+6f9oGdTps2DX3jLFFdmZja9WjcYGuwy9A1XDvq83p0DrFjU+bRfv1oLnvOJ5Kja/PYCf/PWNfyXl69oWJeZWTP5UBLw5y870O/ssl79xRvZtGPg4A1H2bJniH+/rZffbNzO925+eH8oAOwZqvDxq9ax+rxr+eov/si7Lvkdq879CWdeegt39O6e8GuZmT0drbvHUHcwaX5Hcf/wFe9/Je+46CYAPvj6I9k7VOE9x61i8bwSL/nMzwF4zZdu5C0vWcbnT30RnW15br5/B++65HdPqY4TX3g4P1v/2P7xr/7i3v3D19+zlevvSa7afenKhTyya5Cte5OfdrznuFV87M3P45YHd7DmiPks7CixrW+YjmKexV0lAIbKVXp3DvKcpV2ZPZF1vbt4ZNcQg+UKr1h9KIs6k/a1CIr5HKVCy35fMDNaORjq9hg6SnkA3vzCZ3DM6sVc88FXc+0dj/KRNz4384F6z+dO5Pmf+hkAP1n3KD9Z9+hTfv33HLeKT5z8gsyH8HClyrlX38G/37aZ773vFdxw91a+9eukP6fbN+3KPP/S3zzIpb95sOGy2wo5hht08VHIiUpt4jcnOmJBO686cgnHHXkozzq0i65SgVVLOhkaqdFRyj/lINnRP0JeYkFn8eCNJ6BSrfF4/wjD5RpLD2ljz1CZvuEK89uLzGsrsK9T3V0DZXYNJL+AL1drHDqvxKLOEjv6R9g9WKaYF9Va8reydc8w1QhygvZint6dA2zfO0LvzgF2DZbZO1RhUWeJjlKOpfPaWbagnVIhx/yOApJoL+QZqlTZvneYYj7HYYe0gWDJvDbyOdFezNNWyDE4koT5PY/t4bE9Q+weKBORhPbOgRFqAcW8WLaggxWLOljcVWLXQJn+4QpKaztsfjt5iZySQ6Yj1RoDwxW62goMlqsMl6vsHa5QzOUYqSZ/J/M7ihRyohZB/3CFHf1lImL/9K62AjmJaq2GJOa1FdgzVGZnf5mdAyN0lpIvJIWcCJL16mrLMzhSo3+4wu7BMlv3DlGtQVsxR/9whV0DZR7ZNUi5FgyXq0iiVMgxMFwhJyEl6z1cqdE/XKWtkGN+R5H2Yo62Qo5qLTikvcjCzmS75nMiL5HPi65SgfntBRZ0FpnfXmRBR5FSIYcQ7cUctYCutjxdpQJtxRzFfC7TqWatFvs/I8rVoJATuVwyYahcZXCkSkcpTyEnCnOwM05FzO67mHV3d0dPT8+En/exK//Albf2AvDdv3oFV/Zs4pw3PY+Viw98/uDxvmHecdFN3Let/4DtTnjBMzjj2Gdy5GHzKORy7OgfYc0R8ydcJyQfdD+98zG+/7uHOfnFy5jfUeS7v32Imx/Ysb/Na45awm0P72LFog7ueWzv/umvOvJQfv/QLgbLVQDyOVGtBW958TJWLu7kzs27+fXG7SyZV2J+R5FHdw1RyCX/QSUxv73A/dvHXtd8Trxo+QJevHw+z1k6j0WdJXI5MTRSZWCkQi1gW98w92/r496tfezoH2HpvDb2DlV4LP29yIKOIod2lRgYqdI3XKGzlKe9mN//gdRezLFyUSelQo5te4ep1IJ8TvQNV1jUWWRgpMreoQpdbXnyuRz3btnLo7ub81sUpR/uS+a1sXXPECOVGn0jFSbjv1VOsLirDSm5WGJhZ5FiPsdgucqW3UP0j1Sf/osc4LULdcFxIIekgXOwLx2lfI5iXgyWq3QU8xzSXmTVkk5KhTylfA5IQmBeW4FaBLVI6mgr5Oko5ilXa+weLDNSrTFcriFB/0iFnf1l+kcqVKpBtZY8xlP3aB3FPJVajdGrUU0npLnwpPk5QVdbgUPaCuRyorOUp1qLZL0KOQq5JIxK+RwBPN43wmC5Sq0W1CIYGKmidD272gp0FHO0F5P/Azlp/3pXasl701lK3p+I5Avmnz3/sAmvK4CkWyOiu9G8lt1jqP+BmwRfPe1l43reofPauP6c4/nR7ZvJ58SrnrOEH697hNOPeeYBu/E+fMFT/11EIZ/jbX9yBG/7kyP2TzulbrgZ+ocrbNiyl/u39dM3VGZ73wi1CLb3DfPHLX1c2dPbcC9lnyXz2nje4fN48fIF7B4s0zdU4T8dvZyOYp5Hdw/yyK4h5rUXWNxZon+kwsBwlWoklwr3DVX41b3b6Rsus+rQLiq1oFKt0Z4+d9/y9wwm35pXLu7k/a99Nh2lPI/tHqatmHxD3z1YZudAmVotKBWSQ2bLFrSj9Nv1zoEyu/pHkukLO6jWauQk+oerHL6gLfkPWQvKteDQrhJHLOxgYUdx/zfJfZJDeAOUq8m371oke4NthTyLu4qA2Lp3iFoNdgyMUK3VGBipMlSuUSrkWNJV4qXPXMjh89vHvBghItg9mGyHxV0lOop5JBgcqbKtb3j/h1m1FrQVkg+agZEqnaU8+ZxY2Fnc/004nxN7hsqUq4FI9joWpXtxI9Ua5WrQN5RcmJHLQQTsHapwSHuBRZ0lSoUctVqwd6hCNU3Ex/uG6Ruu0FHKp+ud1Fgq5IiIKb/IolytsWewzO50b27XYJlypUak20eCgeHki8hItcZIpZbuJebI55IP/7ySPaiutgLlao1KNdmLKOZzzGsrMFSpUq4EI9Uq/cPJl5Mg2eaFXLLXPlypUotgqFzb3z/bos4SK9o6yOeEJDqLyRGL4UpSz1C5xlC5yp6hMrUaFPJiQbrn1j9cZcueIXJK9mCGK1Pz5aBl9xjOvXodl9+SXIp6yydPYOkhbZNdWkup1YIte4fYsmeY4fTb45J5bZSrNZ4xv31S3t9mfKCYtQrvMTRQ/wHjUHj6crnkuPeyBR1T9hoOBbPmmHtnTcbJnzFmZo21bDAUck4GM7NGWjYY8g4GM7OGWjYYDnQFkZlZK2vZT0cfSjIza8zBYGZmGS0bDPlcy666mdkBteynYyHvPQYzs0ZaNxh8KMnMrKGWDQZfrmpm1ljLBoMvVzUza6xlPx29x2Bm1ljLBoPPMZiZNda6weBDSWZmDbXsp6N3GMzMGmvZYNjX7XZ7sWXfAjOzhlr2U3HfjeuufP9x01uImdkM07LBsM/8jpa9iZ2ZWUMtHwxmZpblYDAzswwHg5mZZTgYzMwso2WDYd9VSWZmltWywbCP8C/dzMzqtXwwmJlZVlODQdKJkjZI2ijp3Abz2yT9IJ1/s6RVzazPzMyaGAyS8sAFwEnAGuB0SWtGNTsT2BkRRwJfAb44VfUMlKtTtWgzs1mtmXsMxwAbI+L+iBgBLgdOHdXmVOA76fBVwBskTclJgNse2glAseBzDGZm9ZrZH8RyYFPdeC/wirHaRERF0m7gUGB7fSNJZwFnpaN9kjY8xZqWHPHF7LJbwBLwOrcAr3NreDrr/KyxZszKjoIi4mLg4qe7HEk9EdE9CSXNGl7n1uB1bg1Ttc7NPJS0GVhZN74indawjaQCsAB4vCnVmZkZ0NxguAU4StJqSSXgNGDtqDZrgXenw/8ZuCHCP0UzM2umph1KSs8ZnA1cB+SBSyJivaTPAj0RsRb4NnCZpI3ADpLwmEpP+3DULOR1bg1e59YwJessfyE3M7N6/uWzmZllOBjMzCyjZYPhYN1zzBaSVkq6UdJdktZL+lA6fbGk/5B0b/rvonS6JH09Xe91ko6uW9a70/b3Snr3WK85U0jKS7pN0jXp+Oq0K5WNadcqpXT6mF2tSDovnb5B0punZ03GR9JCSVdJukfS3ZJeOde3s6T/kf5d3ynp+5La59p2lnSJpK2S7qybNmnbVdLLJd2RPufr4/rRcES03IPk5Pd9wLOBEvAHYM101/UU12UZcHQ6fAjwR5IuR74EnJtOPxf4Yjp8MvBTQMCxwM3p9MXA/em/i9LhRdO9fgdZ948A3wOuScevAE5Lhy8E/ls6/N+BC9Ph04AfpMNr0m3fBqxO/yby071eB1jf7wB/lQ6XgIVzeTuT/OD1AaCjbvu+Z65tZ+C1wNHAnXXTJm27Ar9L2yp97kkHrWm635Rp2hCvBK6rGz8POG+665qkdfsR8EZgA7AsnbYM2JAOXwScXtd+Qzr/dOCiuumZdjPtQfI7mOuB1wPXpH/024HC6G1MciXcK9PhQtpOo7d7fbuZ9iD5Tc8DpBeMjN5+c3E780RPCIvT7XYN8Oa5uJ2BVaOCYVK2azrvnrrpmXZjPVr1UFKj7jmWT1MtkybddX4ZcDPwjIh4NJ31GPCMdHisdZ9t78lXgY8DtXT8UGBXRFTS8fr6M12tAPu6WplN67wa2Ab8S3r47FuSupjD2zkiNgNfBh4GHiXZbrcyt7fzPpO1XZenw6OnH1CrBsOcI2kecDXw4YjYUz8vkq8Kc+a6ZElvBbZGxK3TXUsTFUgON/xzRLwM6Cc5xLDfHNzOi0g61lwNHAF0ASdOa1HTYDq2a6sGw3i655g1JBVJQuG7EfFv6eQtkpal85cBW9PpY637bHpPXgWcIulBkl56Xw98DViopCsVyNY/Vlcrs2mde4HeiLg5Hb+KJCjm8nY+AXggIrZFRBn4N5JtP5e38z6TtV03p8Ojpx9QqwbDeLrnmBXSKwy+DdwdEf9UN6u+e5F3k5x72Df9XenVDccCu9Nd1uuAN0lalH5Te1M6bcaJiPMiYkVErCLZdjdExF8AN5J0pQJPXudGXa2sBU5Lr2ZZDRxFcqJuxomIx4BNkp6XTnoDcBdzeDuTHEI6VlJn+ne+b53n7HauMynbNZ23R9Kx6Xv4rrpljW26T7pM48mek0mu4LkP+OR01/M01uPVJLuZ64Db08fJJMdWrwfuBX4BLE7bi+SGSfcBdwDddcv6S2Bj+njvdK/bONf/eJ64KunZJP/hNwJXAm3p9PZ0fGM6/9l1z/9k+l5sYBxXa0zzur4U6Em39Q9Jrj6Z09sZ+J/APcCdwGUkVxbNqe0MfJ/kHEqZZM/wzMncrkB3+v7dB3yDURcwNHq4SwwzM8to1UNJZmY2BgeDmZllOBjMzCzDwWBmZhkOBjMzy3AwmB2EpFWSQtKMutH8TK3LZj8Hg9nBbSLpjOx2AEnHpx/IS5pVgKRfSvrGgeoymyxNu+ez2WwVEVWSjswmnaRiJN09TNhU1mWtzXsMNqOk34y/KenvJW1Pb2DyZUm5dP6Dkj7a4DnfqBt/UNKnJV0qaa+kTZLeqeRGN5dL6ktvZvKmcda0/5BN2oPtjemsben0S9N2kvRxSfdJGkxvjnJGg+WcLukGSYPA+yUdquQmNL3p89ZLem/d8y4FXgd8IH1+pMt60qEkSa9VcpOaIUlbJH0l7fZlXO9v2ubPldwEZlDSDkn/V9K+3j2tBTgYbCb6C6ACHAecDXwYeOcEl/Fhkm4Rjia5sct3SG7qcy1J1xK/Av5VUvsEl7sJeHs6/EKSQzkfSsc/T9KdwQdIbg5zPnCRpLeMWsb5wDfTNj8k6crh98Bb02V+LX3eG9L2HwJuAv4lfb1lZLtYBkDScpIbsdxG0v36mST9758/qumY76+kw0k6JvwO8AKSm8hcNo73xeaS6e4nxA8/6h/AL4GbRk37D+Bb6fCDwEcbPOcbdeMPAt+vG59H0p/U1+umrUqndY+jpkxbkv6ZAlhS16YLGAReM+q5XwWuHbWcc8bxmpfvW+dG6zhGXX9H0rdOrq7Ne4BhoHOc7+/R6TKfNd1/C35M38PnGGwmWjdq/BHgsKe6jIjokzRA0unYPlvSfye63LGsIfnm/zNJ9R2QFUmCql5P/YikPMm9Fd5JchOVNpJbd/5ygjW8APhtRNTqpv06XdaRPPGeHOj9/QNJp213Svp5OnxVRGybYC02izkYbCYafTI2eOKwZ42kh8l6xXEuozxqHCbvcOq+5byNpLvoA9XSP2r8o8A5JIeM7gD6gL9n8kILsjd6GfP9jYhqeu7lWJKum88Ezpf0uoj4wyTWYzOYg8Fmm20kx9gBSM8RPJ/kuHqzjKT/5uum3UVyyOZZEXHDBJf3auDHEXEZ7L/HxnOBXaNeM9/gufXuBt4hKVe31/Dq9Ln3jbeYiAiScxo3SfossJ5kb8bB0CIcDDbb3AD8paS1JCHxSZr/d/wQybfst0j6MTAYEXslfRn4cvrB/iuScxvHArWIuPgAy/sj8E5Jrya5gf0HSW5nWR92DwLHpFdF9QE7GiznmyQnkr8p6Wsk9y34Asm5iYHxrFh685cTSG78soXkJPZKkuCzFuGrkmy2OZ8kHH4E/JzkGHoz9xaI5Cb1f0tysncLyc1PAD4FfIbk0NB6kpO6bwceOMgiP09yBdVPSQKlH/juqDZfJvnmfxdJID5zjLpOIvkwvx24hOQmMJ+YwOrtJrl95jUkJ7L/EfhcRPzrBJZhs5xv1GNmZhneYzAzswwHg7U8SRemv4Zu9LhwuuszazYfSrKWJ+kwYP4Ys/dExNZm1mM23RwMZmaW4UNJZmaW4WAwM7MMB4OZmWU4GMzMLOP/AyVcqDX20G2IAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Do not learn\n",
    "# use first reward function initially and then switch to second reward function\n",
    "\n",
    "# Instantiate learner in VW\n",
    "vw = pyvw.vw(\"--cb_explore_adf -q UA --quiet --epsilon 0.2\")\n",
    "\n",
    "num_iterations_per_cost_func = 5000\n",
    "cost_functions = [get_cost, get_cost_new1]\n",
    "total_iterations = num_iterations_per_cost_func * len(cost_functions)\n",
    "\n",
    "ctr = run_simulation_multiple_cost_functions(vw, num_iterations_per_cost_func, users, times_of_day, actions, cost_functions, do_learn=False)\n",
    "plot_ctr(total_iterations, ctr)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Scenario 3\n",
    "In this scenario we are going to start rewarding actions that have never seen a reward previously when we change the cost function.\n",
    "\n",
    "### Tom\n",
    "\n",
    "| | `get_cost` | `get_cost_new2` |\n",
    "|:---|:---:|:---:|\n",
    "| **Morning** | Politics |  Politics|\n",
    "| **Afternoon** | Music |   Food |\n",
    "\n",
    "### Anna\n",
    "\n",
    "| | `get_cost` | `get_cost_new2` |\n",
    "|:---|:---:|:---:|\n",
    "| **Morning** | Sports | Food|\n",
    "| **Afternoon** | Politics |  Food |\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [],
   "source": [
    "def get_cost_new2(context,action):\n",
    "    if context['user'] == \"Tom\":\n",
    "        if context['time_of_day'] == \"morning\" and action == 'politics':\n",
    "            return USER_LIKED_ARTICLE\n",
    "        elif context['time_of_day'] == \"afternoon\" and action == 'food':\n",
    "            return USER_LIKED_ARTICLE\n",
    "        else:\n",
    "            return USER_DISLIKED_ARTICLE\n",
    "    elif context['user'] == \"Anna\":\n",
    "        if context['time_of_day'] == \"morning\" and action == 'food':\n",
    "            return USER_LIKED_ARTICLE\n",
    "        elif context['time_of_day'] == \"afternoon\" and action == 'food':\n",
    "            return USER_LIKED_ARTICLE\n",
    "        else:\n",
    "            return USER_DISLIKED_ARTICLE"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "\n",
    "### With learning\n",
    "Let us now switch to the third reward function after a few samples (running the first reward function). Recall that this reward function changes the preferences of the users and is working with a **different** action space than before. We should see the learner pick up these changes and optimize towards the new preferences\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAEOCAYAAACNY7BQAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO3deZhcdZ3v8fe3lq5e0p09ELKQKFEIoAIti6LkCiKLA/cOzkhGARWGe1Uc3B8YZ9TLOOKCOnpli4AwjoqOuESIF1RQRwVM2AIJBJKwJAGSztrptbbv/HFOd6o61Ul36D6nuuvzep56crY69T11OvWp8zunfsfcHRERkT6JuAsQEZHqomAQEZEyCgYRESmjYBARkTIKBhERKaNgEBGRMpEFg5ndYmZbzOyJQeabmX3LzNaa2UozOzaq2kREZI8ojxhuBc7Yx/wzgQXh41Lg+ghqEhGRASILBnf/A7B9H4ucC/y7Bx4AJpnZzGiqExGRPqm4CygxC9hQMr4xnPbSwAXN7FKCowqampqOO/zwwyMpUERkvHjooYe2uvv0SvOqKRiGzN2XAEsAWltbfcWKFTFXJCIytpjZ84PNq6arkjYBc0rGZ4fTREQkQtUUDEuBC8Ork04Edrn7Xs1IIiIyuiJrSjKzHwKLgGlmthH4HJAGcPcbgGXAWcBaoAt4f1S1iYjIHpEFg7sv3s98Bz4cUTkiIjKIampKEhGRKqBgEBGRMgoGEREpo2AQEZEyCgYRESmjYBARkTIKBhERKaNgEBGRMgoGEREpo2AQEZEyCgYRESmjYBARkTIKBhERKaNgEBGRMgoGEREpo2AQEZEyCgYRESmjYBARkTIKBhERKaNgEBGRMgoGEREpo2AQEZEyCgYRESmjYBARkTIKBhERKaNgEBGRMgoGEREpo2AQEZEyCgYRESmjYBARkTIKBhERKaNgEBGRMgoGEREpo2AQGQXZfJFcoRh3GSIHJBV3ASLVpjtboDObZ2dXjvaeHFvae2muT9G2u5fdPTma69OYwXNbu1jX1sGW3T1sbu9l6+5e6lIJCu7s7MoxtamOJRe2ctyhk+PeJNkPdydbKJJKJEgmjGLR6c0XyeaLZNLB92czKBYhkYBC0ckVnHyhSE+4XKFYpFCEnlyBZMJIJozefJFUwii605BOUpdKkE4m6M0X6ckVaO/OkUpa/+v25gsUilCfTpArOJlUgkwqgRnMnNhAfTqJu5MP68ukEtSnkyP+fkQaDGZ2BvBNIAnc5O5fGjB/LnAbMClc5gp3XxZljTL2FYvOru4c2zp7yRedJza1s7m9h60dvXT25jGMSU1pNm7vBoPNu3rIF52O3jw7OrPs6MpS9KG91sEt9cya3MARM5uZcthUurIF6pIJWhrS3LXyJRYveYDLT1vAe088lIkN6dHd8DGm74MRIJEwdnXngg/KRIJEIvgQ3tbZy5bdvbTt7qUumWBXd47N7T1s78ySLRSpSyboyRf63/e6VPAhnisUaUgnMTPcnZ5ckbpUgvaeHN3ZAulkgi27e9jdkyeZMPIFpztXIJUw6lIJurKFON+aIfvi/zqavzth7oiv19yH+D/glb6QWRJ4Gng7sBFYDix299UlyywBHnH3681sIbDM3efta72tra2+YsWK0Stcqo67s35rJy9s72Lj9i7+vG4bT2/eza7uPO7O7t482fzezTh1yQSNmSTd2QLZQpGDmusxgxkt9RgwvTnDtAkZpk2oY2JDmub6FE2ZFDMn1tPZW+CQSQ0016fYuKObpkySQ6c00VA3+Le1HZ1ZPn3HSn69ejMAh0ys58I3zWPxG+cysTH+kGjvyfHyruDDsSubJ51M0JMr8NKuHqZNyNCUSbKrK0d3rkAhDNuW+jQNdUl2defo7M2TKxR5ub2Hzt4CEzIpevPBh27RoaM3T9vuHl7eFXyQFx0yqQTZfBEzaO/JD7vmhMGM5nomN9X1r6s+naApk6KjN4978M0+lTC6c4X+8YZ0kt58keb6FA3pFLlCkRnNGVoa0uQLRdLJBJMa03RlC2TzRRrrkjTUpUgnjWyhiHvwd5dIGO6QMOsPkUwqQSadCKcF4wV3CkWnLpkgX3RSCaMzmydXcHrzBRrSSRrSSZrr0+SKRQoFp+BOfTpJwqA3VySVNHpyBXIF73+fs/kiCQuORjKpBG8+bBpHzGw5oP1vZg+5e2uleVEeMRwPrHX39WFRtwPnAqtLlnGgbysnAi9GWJ9Ugd58gZd29jB3SiMd2TxPvtjO01s6WP1iO+09OZ7f1skTm9rLnjO1qY4jZrbwhjn1JBPQUp9m5qQGWupTJMw4YmYLs6c0MKEuRSJsJsgWigd8CH5QS/2QlpvcVMd3LmzlkRd2sPSxF1m1qZ0v/eopvn7P05zy2umcedTBzJ7cyPTmDHOnNJJM2F7rKBadrvCDeUdnlp3dOdp29/LSrm66s8GHeL5YpC6ZJJ00JjamcYdkwtjdk+PFncEHc0++QE+uyM6uLLt78mQLxYrheSBawgDt6M3TWJckmy+STCSYkEkybUKGo2ZNZGpTHYmEkc0HH8LuzvTmDBMb68gXgg+7lobgQzvgJMyY0lTHjOZ6ZrRk6M4WaGlIV3yfZGRFGQyzgA0l4xuBEwYs83ngHjP7CNAEnFZpRWZ2KXApwNy5I38YJSOnszfP6pfaWd/WQTZfZHpzPT25Amu3dNCVLZArFEkmjA3bu9i0s5tntnRQqNCO01SXJJkwFhzUzPveNI9XTW9i5sQGZk9u4PCDmzEb+odFImHUJ0a+XXYwx8ydzDFzg/MMj23Yye3LX+DXqzf3H0kATJsQhFt7T55CsUh7d56dXdn+9uvBNKSTZNLBN+d8IQi8PsmEMaM5w0Et9TRlkkGAHtxMc32KTDoYP3hiPfXpJBMyKYrhN9aDW+rZ1pmlqzfPpMY6GuuSmAWBu6MrS77otNQHR1QJs30eNY2k0WhLl8qibEp6F3CGu18Sjl8AnODul5Us8/Gwpq+Z2UnAzcBR7j7o/ww1JUXDPWhKSCeDdtofPvgCE+pTPPz8ThozSa4443BmtNTTnS1wy5+e5Ru/fpr8UBvqCdrq505t5Jg5k2ioS7Jldy8zmjMsmNHM3CmNLDykZVx9U8wViqxr6+CZzR1s78zyx7Vb2bijm8mNaRJmTGxIM7kpTWNdiqlNdZjBlKYMExvSzGjOMHNSPQ3hB3ppKHb05ul7mzKp5Lh6z2RkVUtT0iZgTsn47HBaqYuBMwDc/X4zqwemAVsiqbDGFYqOu5NMGGu3dHDnypdYuXEnv3u6jf19f/jpwwN3ZWDx8XM4+bDpzJrcQCph7OjKkkokeMOcSdSlgisxcmEbby1JJxMcfnALhx8ctJxe9KZ5I7LeCRldaCivXJR/RcuBBWY2nyAQzgf+bsAyLwCnArea2RFAPdAWYY016dENOznv+j9XbMLpc+jURo6fN4UnX24nX3COmTuJs46eybQJGX7+6CZu/P165k5p5OhZE3nPCXM56dVTh9y8U2uhIFLtIgsGd8+b2WXA3QSXot7i7qvM7CpghbsvBT4BfMfMPkZwIvp9HlVbVw34xaObuPz2R/e5TEt9ioa6JDOa6/ngolfzutkTOWRiA4l9NEkcMbOFK888YqTLFZGYRHrcGf4mYdmAaZ8tGV4NvDnKmsazvks3r7zjce56/KW95jfXpzhkYgPHzJ3EhxYdxtypjTFUKSLVRg2S40ix6PxoxQYa0kn+8WeP7/UjnbOPnskHF72aw2ZM0BUeIjIoBcMYVyg6m9t7eM9ND/Ls1s6Ky9x4wXGcvvCgYV3SKSK1S8EwBnX25jnv+j/z1Mu7K84/8pAWLj91AacfeXDElYnIeKBgGCN6cgW+evcabv7jsxXn3/aB43nrgmk6KhCRV0zBMAbs7slx4hd/S2fJOYMLTjyUT53xWlrq4+9zR0TGFwVDldvdk+Poz98DBJ2w3fvJRTpxLCKjSsFQxa64YyW3Lw+6l7r81AV87O2vibkiEakFCoYq4+7Mv7L8FhR//5b5CgURiYyCoYrc+qdn+fwvV5dNe+ifTmPqhExMFYlILVIwVIl1bR39oTBnSgP3fWIRKfUhJCIxUDBUgYee38F51/8ZgGX/8BYWHnJgd2QSERkJCoYYZfNFXvNPv+ofv+ZvXq9QEJHYqa0iRhfd8pf+4Q8uejXvOm52jNWIiAR0xBCRB9Zv47O/eIJMKslbFkzjut+tA4IeTld+7nT9YllEqoaCIQKPbdjJ+Use6B9/fNOu/uH/+vT/UCiISFVRMIwyd+fca/+01/RJjWlWfOY0XXkkIlVHwTDKPvT9hwE47YgZ3HTRG2OuRkRk//R1dRR19ub51RMvA3Dte46NuRoRkaFRMIyiz/zscQBuvqiVTEod34nI2KBgGCX5QpGfP/oiAG87fEbM1YiIDJ2CYRQUi95/FdK3Fh+jq45EZEzRyedR8MZ//Q3bOrMAnHWUbq8pImOLjhhG2APrt/WHwpILjtPlqCIy5uhTawQVSpqQrv7rozn9SB0tiMjYo2AYQQ+/sKN/ePHxc2OsRETkwCkYRtBHb38UgBX/dFrMlYiIHDgFwwjZ3pll085uAKbpjmsiMoYpGEbIBTc/CMCFJx0acyUiIq+MgmEEbGnvYdWL7QBcde5RMVcjIvLKKBhGwM8e2QTA+988L95CRERGgH7g9grMu+KusvHP/dWRMVUiIjJyFAwHwN255p41cZchIjIqFAwHYP6Vy8rGb7qwlVNeOz2makRERpaCYZjuW7OlbPy5L50dUyUiIqMj0pPPZnaGma0xs7VmdsUgy/ytma02s1Vm9oMo6xuK9393ef/w7z+1KL5CRERGSWRHDGaWBK4F3g5sBJab2VJ3X12yzALgSuDN7r7DzKrqRgb/9Uxb//CzV5+l7rRFZFyK8ojheGCtu6939yxwO3DugGX+HrjW3XcAuPsWqoS7c8HNfwHgkpPnKxREZNyKMhhmARtKxjeG00q9BniNmf3JzB4wszMqrcjMLjWzFWa2oq2trdIiI670hPMVZx4eyWuKiMSh2n7glgIWAIuAxcB3zGzSwIXcfYm7t7p76/Tpo381UG++0D98/XuO1T0WRGRci/ITbhMwp2R8djit1EZgqbvn3P1Z4GmCoIjVb1bvadE68+iZMVYiIjL6ogyG5cACM5tvZnXA+cDSAcv8nOBoATObRtC0tD7CGiv68A8eBmDpZW+OuRIRkdEXWTC4ex64DLgbeBL4sbuvMrOrzOyccLG7gW1mthq4D/iUu2+Lqsb9OXrWxLhLEBEZdZH+wM3dlwHLBkz7bMmwAx8PH1UhKAleP2eSrkQSkZqgs6j70d6TB+CdOrcgIjVCwbAfL+0K7sp28MT6mCsREYmGgmE/Nrf3AjBTwSAiNULBsB+b23sAOKhFwSAitUHBsB9bwmCY3pyJuRIRkWgoGPZjc3svkxrT1KeTcZciIhIJBcN+rNy0i6lNdXGXISISGd2oZx92dmV5bMPOuMsQEYmUjhj24dxr/xR3CSIikRvSEYOZpQm6zD7V3VeNbknxW/7cdv7wdBvPb+sC4OaLWmOuSEQkOkMKBnfPmVkO8FGupyr8zQ33l42fesRBMVUiIhK94TQl/T/gSjMb1+clcoVi3CWIiMRqOB/ybwFOATaZ2RNAZ+lMdz+n4rPGmMc37Sob//2nFsVTiIhITIYTDFuBO0arkGpx18qX+oef+9LZMVYiIhKPIQeDu79/NAupFq+bHdxz4Y4PnhRzJSIi8RjyOQYzu7fS/ZfNrMXM7h3ZsuLz1bvXADBnSmPMlYiIxGM4J58XAZV+AlxPcP5hXNi4I+hme2JDOuZKRETisd+mJDM7tmT0dWa2vWQ8CbwD2DTShcUtk1LfSCJSm4ZyjmEFwe8XHLinwvxu4CMjWVTcXj9nrxYzEZGaMZRgmA8YsB54I8HVSX2ywBZ3L4xCbZF7Ifyl89bdvTFXIiISn/0Gg7s/D2BmXwLe6O43lM43s/9jZrPc/Z9HqcbI/M/rgr6RNu3sjrkSEZH4DOfk83uBRypMfxi4cGTKidf2ziwAHz1tQcyViIjEZzjBMANoqzB9KzDmOxPa1ZXrH37n62bGWImISLyGEwwvUPmy1LcCG0emnPh8+o7H+ocPm9EcYyUiIvEaTpcYNwLfMLM6oO8HbacCVwNfHunCovbntdviLkFEpCoMp0uMr5nZNOBb7PmhWxb4prt/ZTSKi9Lu3nzcJYiIVIVhdaHt7lea2ReAheGkJ929Y+TLEhGRuAz71p7u3unuy8PHuAmFdxwZnD+/6x9OjrkSEZF46Z7PoWy+yFGzWjjykIlxlyIiEqtxfTe24VjX1kldSjkpIqJgCL2wvSvuEkREqoK+IouISBkFA+DuAFxy8vyYKxERiZ+CgT0357npj8/GXImISPwUDMCal3cDcML8KTFXIiISv0iDwczOMLM1ZrbWzK7Yx3LnmZmbWeto1/Tc1k4u+fcVACw+fu5ov5yISNWLLBjMLAlcC5xJ8MvpxWa2sMJyzcDlwINR1PX05t39wxldrioiEukRw/HAWndf7+5Z4Hbg3ArL/QtBp3w9URS15A/r+4f1OwYRkWiDYRawoWR8Yzitn5kdC8xx97v2tSIzu9TMVpjZira2SreIGLqOks7zMqnkK1qXiMh4UDVfkc0sAXwd+MT+lnX3Je7e6u6t06dPf0Wv+9TLJU1J6ap5O0REYhPlJ+EmYE7J+OxwWp9m4Cjgd2b2HHAisDSKE9B90kkFg4hIlJ+Ey4EFZjY/vNnP+cDSvpnuvsvdp7n7PHefBzwAnOPuKyKsUUSk5kUWDO6eBy4D7gaeBH7s7qvM7CozOyeqOvalUPS4SxARiV2knei5+zJg2YBpnx1k2UVR1HTesbO54+HgltUt9epTUESk5hvV88UiMyfW84NLTmDBQc1xlyMiEruaD4ZcoUhTJsWbDpsWdykiIlVBwVBwXY0kIlKi5j8Rc4Ui6aTFXYaISNVQMBSKOmIQESlR85+IQVOSjhhERPooGHTEICJSpuYv3H/khZ1Mm1AXdxkiIlVDX5WBrR3ZuEsQEakaNR0M2zsVCCIiA9V0MPz+6S1xlyAiUnVqOhgSpquRREQGUjCIiEiZmg6GZELBICIyUE0Hg3JBRGRvNR0MoGQQERmoxoNBREQGqulgaKxLAnDO6w+JuRIRkepR08HQd1HSBScdGm8hIiJVpKaDoVB0QJetioiUqulg8CAXdNmqiEiJmg6GPUcMMRciIlJFajoYiq6mJBGRgRQMKBhERErVeDAE/+ocg4jIHjUdDB09eUDnGEREStV0MFzx05UArGvrjLkSEZHqUdPB0NeUpFMMIiJ71HQw9Dlh/pS4SxARqRoKBiCV1NsgItJHn4hASmefRUT6KRhQMIiIlFIwoN8xiIiUqulgyKSCzTddliQi0i/SYDCzM8xsjZmtNbMrKsz/uJmtNrOVZvZbMxvVGyWcesQMDpsxYTRfQkRkzIksGMwsCVwLnAksBBab2cIBiz0CtLr764CfAF8ZzZqKRUjqaEFEpEyURwzHA2vdfb27Z4HbgXNLF3D3+9y9Kxx9AJg9mgUV3PXjNhGRAaIMhlnAhpLxjeG0wVwM/KrSDDO71MxWmNmKtra2Ay7I3XXiWURkgFTcBVRiZu8FWoFTKs139yXAEoDW1lY/0Ne5f902Cn7ATxcRGZeiDIZNwJyS8dnhtDJmdhrwGeAUd+8dzYI6s4XRXL2IyJgUZVPScmCBmc03szrgfGBp6QJmdgxwI3COu2+JsDYREQlFFgzungcuA+4GngR+7O6rzOwqMzsnXOyrwATgP83sUTNbOsjqRERklER6jsHdlwHLBkz7bMnwaVHWIyIie6vpXz6LiMjeFAwiIlJGwSAiImUUDCIiUkbBICIiZWo+GF41vSnuEkREqkpNB0PC4OyjZ8ZdhohIVanpYFAvSSIie6vtYHBQ36oiIuVqOhgAdEMGEZFyCgYRESlT88Gg4wURkXI1GwyuG/SIiFRUw8EQ/KtTDCIi5Wo2GPqYGpNERMrUbDCoIUlEpLLaDYawLUlNSSIi5Wo2GPooF0REytVsMKgpSUSkspoNhj5qShIRKVezwaCfMYiIVFa7wUDfyWcdMoiIlKrZYBARkcpqNhjUlCQiUlnNBkMftSSJiJRTMOiXDCIiZWo2GNSUJCJSWc0GQx81JYmIlKvZYHD99llEpKLaDYa++zHEW4aISNWp2WDoo6YkEZFyNRsMakgSEamsdoOh734MakwSESlTs8HQR01JIiLlajYY1JQkIlJZpMFgZmeY2RozW2tmV1SYnzGzH4XzHzSzeVHWJyIiEQaDmSWBa4EzgYXAYjNbOGCxi4Ed7n4Y8A3gy6NVT1dvYbRWLSIypkV5xHA8sNbd17t7FrgdOHfAMucCt4XDPwFOtVG6YcLty18AIJOq2dY0EZGKUhG+1ixgQ8n4RuCEwZZx97yZ7QKmAltLFzKzS4FLw9EOM1tzgDVNu/DLbL3wAJ88Rk1jwPtZA7TNtUHbPDyHDjYjymAYMe6+BFjyStdjZivcvXUEShoztM21QdtcG0Zrm6NsR9kEzCkZnx1Oq7iMmaWAicC2SKoTEREg2mBYDiwws/lmVgecDywdsMxS4KJw+F3Ave7qIFtEJEqRNSWF5wwuA+4GksAt7r7KzK4CVrj7UuBm4HtmthbYThAeo+kVN0eNQdrm2qBtrg2jss2mL+QiIlJK12qKiEgZBYOIiJSp2WDYX/ccY4WZzTGz+8xstZmtMrPLw+lTzOzXZvZM+O/kcLqZ2bfC7V5pZseWrOuicPlnzOyiwV6zWphZ0sweMbM7w/H5YVcqa8OuVerC6YN2tWJmV4bT15jZO+LZkqExs0lm9hMze8rMnjSzk8b7fjazj4V/10+Y2Q/NrH687Wczu8XMtpjZEyXTRmy/mtlxZvZ4+JxvDelHw+5ecw+Ck9/rgFcBdcBjwMK46zrAbZkJHBsONwNPE3Q58hXginD6FcCXw+GzgF8R3LzuRODBcPoUYH347+RweHLc27efbf848APgznD8x8D54fANwAfD4Q8BN4TD5wM/CocXhvs+A8wP/yaScW/XPrb3NuCScLgOmDSe9zPBD16fBRpK9u/7xtt+Bt4KHAs8UTJtxPYr8JdwWQufe+Z+a4r7TYlpR5wE3F0yfiVwZdx1jdC2/QJ4O7AGmBlOmwmsCYdvBBaXLL8mnL8YuLFketly1fYg+B3Mb4G3AXeGf/RbgdTAfUxwJdxJ4XAqXM4G7vfS5artQfCbnmcJLxgZuP/G435mT08IU8L9difwjvG4n4F5A4JhRPZrOO+pkullyw32qNWmpErdc8yKqZYREx46HwM8CBzk7i+Fs14GDgqHB9v2sfae/BvwaaAYjk8Fdrp7Phwvrb+sqxWgr6uVsbTN84E24Lth89lNZtbEON7P7r4JuAZ4AXiJYL89xPjez31Gar/OCocHTt+nWg2GccfMJgB3AB919/bSeR58VRg31yWb2TuBLe7+UNy1RChF0NxwvbsfA3QSNDH0G4f7eTJBx5rzgUOAJuCMWIuKQRz7tVaDYSjdc4wZZpYmCIXvu/tPw8mbzWxmOH8msCWcPti2j6X35M3AOWb2HEEvvW8DvglMsqArFSivf7CuVsbSNm8ENrr7g+H4TwiCYjzv59OAZ929zd1zwE8J9v143s99Rmq/bgqHB07fp1oNhqF0zzEmhFcY3Aw86e5fL5lV2r3IRQTnHvqmXxhe3XAisCs8ZL0bON3MJoff1E4Pp1Udd7/S3We7+zyCfXevu78HuI+gKxXYe5srdbWyFDg/vJplPrCA4ERd1XH3l4ENZvbacNKpwGrG8X4maEI60cwaw7/zvm0et/u5xIjs13Beu5mdGL6HF5asa3Bxn3SJ8WTPWQRX8KwDPhN3Pa9gO04mOMxcCTwaPs4iaFv9LfAM8BtgSri8EdwwaR3wONBasq4PAGvDx/vj3rYhbv8i9lyV9CqC//Brgf8EMuH0+nB8bTj/VSXP/0z4XqxhCFdrxLytbwBWhPv65wRXn4zr/Qz8X+Ap4AngewRXFo2r/Qz8kOAcSo7gyPDikdyvQGv4/q0Dvs2ACxgqPdQlhoiIlKnVpiQRERmEgkFERMooGEREpIyCQUREyigYRESkjIJBZD/MbJ6ZuZlV1Y3mq7UuGfsUDCL7t4GgM7JHAcxsUfiBPC2qAszsd2b27X3VJTJSIrvns8hY5e4Fgo7MRpyZpT3o7mHYRrMuqW06YpCqEn4zvs7MvmhmW8MbmFxjZolw/nNm9skKz/l2yfhzZvZZM7vVzHab2QYze7cFN7q53cw6wpuZnD7EmvqbbMIebO8LZ7WF028NlzMz+7SZrTOz7vDmKO+tsJ7FZnavmXUD/9vMplpwE5qN4fNWmdn7S553K3AK8OHw+R6ua6+mJDN7qwU3qekxs81m9o2w25chvb/hMn9twU1gus1su5n93sz6eveUGqBgkGr0HiAPvAm4DPgo8O5hruOjBN0iHEtwY5fbCG7qs4yga4k/AP9hZvXDXO8G4Lxw+EiCppzLw/EvEHRn8GGCm8NcDdxoZmcPWMfVwHXhMj8n6MrhYeCd4Tq/GT7v1HD5y4H7ge+GrzeT8i6WATCzWQQ3YnmEoPv1iwn63796wKKDvr9mdjBBx4S3AUcQ3ETme0N4X2Q8ibufED30KH0AvwPuHzDt18BN4fBzwCcrPOfbJePPAT8sGZ9A0J/Ut0qmzQuntQ6hprJlCfpncmBayTJNQDfwlgHP/Tdg2YD1fGIIr3l73zZX2sZB6vpXgr51EiXLvA/oBRqH+P4eG67z0Lj/FvSI76FzDFKNVg4YfxGYcaDrcPcOM+si6HSsz+bw3+GudzALCb75/38zK+2ALE0QVKVWlI6YWZLg3grvJriJSobg1p2/G2YNRwAPuHuxZNofw3Udxp73ZF/v72MEnbY9YWb3hMM/cfe2YdYiY5iCQarRwJOxzp5mzyJBD5Ol0kNcR27AOIxcc2rfev6KoLvofdXSOWD8k8AnCJqMHgc6gC8ycqEF5Td6GfT9dfdCeO7lRIKumy8GrjazU9z9sRGsR0gmmWcAAAGSSURBVKqYgkHGmjaCNnYAwnMEhxO0q0clG/6bLJm2mqDJ5lB3v3eY6zsZ+KW7fw/677HxGmDngNdMVnhuqSeBvzWzRMlRw8nhc9cNtRh3d4JzGveb2VXAKoKjGQVDjVAwyFhzL/ABM1tKEBKfIfq/4+cJvmWfbWa/BLrdfbeZXQNcE36w/4Hg3MaJQNHdl+xjfU8D7zazkwluYP8RgttZlobdc8Dx4VVRHcD2Cuu5juBE8nVm9k2C+xZ8ieDcRNdQNiy8+ctpBDd+2UxwEnsOQfBJjdBVSTLWXE0QDr8A7iFoQ4/yaAEPblL/OYKTvZsJbn4C8M/A5wmahlYRnNQ9D3h2P6v8AsEVVL8iCJRO4PsDlrmG4Jv/aoJAnDtIXWcSfJg/CtxCcBOYfxzG5u0iuH3mnQQnsr8G/Iu7/8cw1iFjnG7UIyIiZXTEICIiZRQMUvPM7Ibw19CVHjfEXZ9I1NSUJDXPzGYALYPMbnf3LVHWIxI3BYOIiJRRU5KIiJRRMIiISBkFg4iIlFEwiIhImf8GZLK3cC6X/60AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# use first reward function initially and then switch to third reward function\n",
    "\n",
    "# Instantiate learner in VW\n",
    "vw = pyvw.vw(\"--cb_explore_adf -q UA --quiet --epsilon 0.2\")\n",
    "\n",
    "num_iterations_per_cost_func = 5000\n",
    "cost_functions = [get_cost, get_cost_new2]\n",
    "total_iterations = num_iterations_per_cost_func * len(cost_functions)\n",
    "\n",
    "ctr = run_simulation_multiple_cost_functions(vw, num_iterations_per_cost_func, users, times_of_day, actions, cost_functions)\n",
    "\n",
    "plot_ctr(total_iterations, ctr)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Without Learning"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAEOCAYAAACNY7BQAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAeoklEQVR4nO3deZhcdZ3v8fe3tt67s4dskAhBjIxIaBlUFEZcEuTCPMM8Sq5el0HRueLoqHceuG4jooiDGyOKuY6jgwqi49UMgxOUZbjOiEMDggkYDSGQDiTprL13bd/7xzmd1OlUd7pD96nqrs/reerJWX516nvqdOpTZ6nfMXdHRERkWKLSBYiISHVRMIiISISCQUREIhQMIiISoWAQEZEIBYOIiETEFgxm9i0z22Nmm0aZb2Z2o5ltNbPHzGx1XLWJiMgRce4xfBtYM8b8tcDK8HEF8PUYahIRkRFiCwZ3vx/YP0aTS4B/8sADwCwzWxRPdSIiMixV6QJKLAF2lIx3htOeG9nQzK4g2KugqanprNNOOy2WAkVEZoqHHnpor7vPLzevmoJh3Nx9PbAeoL293Ts6OipckYjI9GJmT482r5quStoJLCsZXxpOExGRGFVTMGwA3hZenXQOcMjdjzqMJCIiUyu2Q0lmditwPjDPzDqBTwJpAHe/GbgTuBDYCvQD74yrNhEROSK2YHD3dceY78D7YipHRERGUU2HkkREpAooGEREJELBICIiEQoGERGJUDCIiEiEgkFERCIUDCIiEqFgEBGRCAWDiIhEKBhERCRCwSAiIhEKBhERiVAwiIhIhIJBREQiFAwiIhKhYBARkQgFg4iIRCgYREQkQsEgIiIRCgYREYlQMIiISISCQUREIhQMIiISoWAQEZEIBYOIiEQoGEREJELBICIiEQoGERGJUDCIiEiEgkFERCIUDCIiEqFgEBGRCAWDiIhEKBhERCRCwSAiIhGxBoOZrTGzLWa21cyuKjP/RDO718weMbPHzOzCOOsTEZEYg8HMksBNwFpgFbDOzFaNaPYx4HZ3PxO4DPhaXPWJiEggzj2Gs4Gt7r7N3bPAbcAlI9o40BoOtwHPxlifiIgQbzAsAXaUjHeG00r9LfBWM+sE7gTeX25BZnaFmXWYWUdXV9dU1CoiUrOq7eTzOuDb7r4UuBC4xcyOqtHd17t7u7u3z58/P/YiRURmsjiDYSewrGR8aTit1OXA7QDu/iugHpgXS3UiIgLEGwwPAivNbIWZZQhOLm8Y0eYZ4AIAM3sRQTDoWJGISIxiCwZ3zwNXAhuBJwiuPtpsZteY2cVhsw8D7zazR4FbgXe4u8dVo4iIQCrOF3P3OwlOKpdO+0TJ8OPAK+OsSUREoqrt5LOIiFSYgkFERCIUDCIiEqFgEBGRCAWDiIhEKBhERCRCwSAiIhEKBhERiVAwiIhIhIJBREQiFAwiIhKhYBARkQgFg4iIRCgYREQkQsEgIiIRCgYREYlQMIiISISCQUREIhQMIiISoWAQEZEIBYOIiEQoGEREJELBICIiEQoGERGJUDCIiEiEgkFERCIUDCIiEqFgEBGRCAWDiIhEKBhERCRCwSAiIhEKBhERiVAwiIhIhIJBREQiYg0GM1tjZlvMbKuZXTVKmzeZ2eNmttnMvh9nfSIiAqm4XsjMksBNwOuATuBBM9vg7o+XtFkJXA280t0PmNmCuOoTEZFAnHsMZwNb3X2bu2eB24BLRrR5N3CTux8AcPc9MdYnIiLEGwxLgB0l453htFKnAqea2X+Y2QNmtqbcgszsCjPrMLOOrq6uKSpXRKQ2VdvJ5xSwEjgfWAf8HzObNbKRu69393Z3b58/f37MJYqIzGxxBsNOYFnJ+NJwWqlOYIO759z9KeD3BEEhIiIxiTMYHgRWmtkKM8sAlwEbRrT5CcHeAmY2j+DQ0rYYaxQRqXmxBYO754ErgY3AE8Dt7r7ZzK4xs4vDZhuBfWb2OHAv8L/cfV9cNYqICJi7V7qG56W9vd07OjoqXYaIyLRiZg+5e3u5edV28llERCpMwSAiIhEKBhERiVAwiIhIhIJBREQiFAwiIhKhYBARkQgFg4iIRCgYREQkYlzBYGZpM9tlZi+e6oJERKSyxhUM7p4DcsD07j9DRESOaSKHkv4euNrMYrsdqIiIxG8iH/KvAs4DdprZJqCvdKa7X1z2WSIiMq1MJBj2Av88VYWIiEh1GHcwuPs7p7IQERGpDuM+x2Bm95S7/7KZtZrZPZNb1tR7dMdB3vWdDvKFYqVLERGpKhM5lHQ+kCkzvZ7g/MO08oHbHmH7vn52HBhgxbymSpcjIlI1jhkMZra6ZPQlZra/ZDwJvAHYOdmFiYhIZYxnj6GD4PcLDtxVZv4A8P7JLCpO0/3WpiIik208wbACMGAb8DKCq5OGZYE97l6YgtqmlJlVugQRkap0zGBw96cBzOxzwMvc/ebS+Wb2XjNb4u4fn6IaRUQkRhP55fNbgUfKTH8YeNvklBM/HUgSEYmaSDAsALrKTN8LLJyccuKjA0kiIuVNJBieofxlqa8GOienHBERqbSJ/I7hG8CXzCwDDP+g7QLgOuD6yS4sLrooSUQkaiJdYnzBzOYBN3Lkh25Z4Cvu/vmpKG5K6ViSiEhZE+pC292vNrNrgVXhpCfcvXfyyxIRkUqZ8L0V3L0PeHAKaqkQHUsSESlVs/d81pEkEZHyajYYtJ8gIlJezQaDiIiUV7PBoENJIiLl1WwwDNPvGEREomo2GNS7qohIeTUbDCIiUl6swWBma8xsi5ltNbOrxmh3qZm5mbVPdU06kiQiEhVbMJhZErgJWEvwy+l1ZraqTLsW4APAr6e0nqlcuIjINBbnHsPZwFZ33+buWeA24JIy7T5N0CnfYIy1iYhIKM5gWALsKBnvDKcdZmargWXu/q9jLcjMrjCzDjPr6Ooqd4uI8dNVSSIiUVVz8tnMEsAXgQ8fq627r3f3dndvnz9//nG+3nE9TURkxoszGHYCy0rGl4bThrUApwP3mdl24BxgQxwnoEVE5Ig4g+FBYKWZrQhv9nMZsGF4prsfcvd57r7c3ZcDDwAXu3vHVBblui5JRCQitmBw9zxwJbAReAK43d03m9k1ZnZxXHUMM12XJCJS1oTvx/B8uPudwJ0jpn1ilLbnT2kt2lMQESmrak4+i4hIdajZYNChJBGR8mo2GEREpLyaDwb9wE1EJKpmg0E/cBMRKa9mg2GY9hhERKJqPhhERCRKwSAiIhE1Hwz6oZuISFTNBoPu+SwiUl7NBoOIiJRX88Ggq5JERKJqPhhERCSq5oNBpxpERKJqPhh0KElEJKpmg0E7CiIi5dVsMIiISHkKBhERiajZYNBJZxGR8mo2GEREpLyaD4ZcoVjpEkREqkrNBsPmZ7sBuOGuLRWuRESkutRsMAz7w+7eSpcgIlJVaj4YivqBm4hIRM0HA7ofg4hIRM0Hg/YYRESiaj4YXJ0liYhEKBgqXYCISJWp+WAo6liSiEhEzQeDjiSJiETVfDAUlQwiIhE1HwwFBYOISETNB0NRXSWJiETUfDBoj0FEJCrWYDCzNWa2xcy2mtlVZeZ/yMweN7PHzOxuMztpqmsq6KokEZGI2ILBzJLATcBaYBWwzsxWjWj2CNDu7i8BfgR8furqmaoli4hMb3HuMZwNbHX3be6eBW4DLilt4O73unt/OPoAsHSqikkqGUREyoozGJYAO0rGO8Npo7kc+Fm5GWZ2hZl1mFlHV1fXcRWTTIwvGPb0DPLU3r7jeg0RkemoKk8+m9lbgXbg78rNd/f17t7u7u3z588/ztc4etqenkHe8KX72bG/f/h1OPszd/MnN9zH73Z1H9friIhMN3EGw05gWcn40nBahJm9FvgocLG7D01VMcbRyXD2Z+5my+4eXvX5ewHYuHn34Xlrvvz/2LTzEACPdR7kv57aP1WliYhUVCrG13oQWGlmKwgC4TLgv5c2MLMzgW8Aa9x9z1QWc6xTDIf6c7z3uw9Fpl3097+MjN/94fM4eX7zZJcmIlJRse0xuHseuBLYCDwB3O7um83sGjO7OGz2d0Az8EMz+42ZbZiqeo51huGMa+46PPyh151ats0FX/h3/u8jnZNYlYhI5cW5x4C73wncOWLaJ0qGXxtXLTZil+FjP/ntqG3/6oKVXPhHJ/DaL95/1Ly//sGjvPGPFpNJVeXpGhGRCYs1GKpJaSx0Hujnuw88U7bdts9eCMApC1rY/rk3Rua955YONm7ezQ13beHqtacdFTYiItNR7X7NLfkMP/f6e3nXuSsOj7/nvBcAkDBIjHFZ681vPQuA9fdvY8XVd/LQ0wemplYRkRhpjyH0zV8+BcDVa0/jPeedzAsXtnDqwpaxl2HG19+ymr/83sMAXPr1/wTglAXNdPUMcclLF3PBixby6pXztDchItOGTfd7Hre3t3tHR8eEn3fGp+7i0EDuqOkjDxeNx97eIf7q1kf4zyf3jav9T9/3Ss5YNmvCryMiMlnM7CF3by83r3b3GCbxC/y85jq+/+5z6B3Kc/onNx6z/SU3/cfh4YvPWEx/Ns/1l76Euc11k1cUMJgr8PDTB3jpibOoTyXZ3TNIKpEgVyiSMOOEtvpJfb3Jls0X2bqnl/ktddSlE2zZ1UM2X+TQQI6n9vZxaCBH54F+Tl3YwvK5TSxoCd6/VDLBrMY0C1vqaa5PjftX7jPRYK7AQLaAGRzoz9E7mKe5PsWCljoa0kkSCcPd2d+XpejQUp8inUxgjH0YdbJk80We2d+HO+QKTl82T10qQWMm2G49gzkaMylmN6ZJpxK4Q0tdKpbaalntBsMULLO5LsXvr13Lvr4hFrU1HJ7eM5jjFZ+7h57BPIva6nnu0ODheRsefRaAs679xeFpN647kzUvPmHUK522dfWydHYje3uHyOaL3PLA06w5/QQSZiyf28jHf7qJPd1DdBzjnMcLF7Zw5omzeM1pCzh1YQvN9Sl+91wPBXdOnt/E0tmN5ApFftjRycPPHOAlS9toyqTY35cFYPm8JnqHcjTXpXnZ8tl0HhigY/t+eofyPLGrh2LR2dU9yP6+LJlkgmTCWDq7gcZMihPnNNLakAprbmL7vj4yqaDNjv0DbNy8a1xdkcxpyvCzTbvGvEVrYybJgpY6Vp84m9OXtNHWkKbojjtkUgn6snma61I0pJPMb6mjpT5NOml09QwxqzFDYybJ3t4hGtJJkglj8awG6tNJIPh1fKHopJLHPl1XKDoH+7Ps68uyrzdLrlAkk0pQn06SLxQxg3zB2XlwgP5sgc4DA+zuHqRQdPb1DXGwP4cZ9A8V2NeXpbUhRbEYfPgP5YvUp5PMaUrTN1QgmTCy+SK7ugfHrKmlLkV/rnBUL8NmUJ9KUnCnLpmgtSHNvOYMTlBjfzZP92Ce3sE8ZsGdEOtSSYwgXOa31tPWkKalPkVTJslALrjxyf6+ITLJBPv7c+zY33/4b2miEka4vsH2mdOUYW5zEHZFdwwjmQjCZk9P8B4Mr89AtkBLfZpCschQvkh/tkAmmcAJntdUl6Qxk6Ihk6RQdMygtT5NrhCsQyphpJIJ6lIJDg7kKBScRMIYzBUwYDBf4GB/jmTCGMoXyRWKdA/kKHrw3IQZfdk8g7kiCYN0MkEiAalEgqFcgfrw7yxhxnOHBnCCz6vGTIqW+hT5omNAtlDkI69/IX965lg9Cx2fmj2UtPrTPz/qj/Jd567gYxeN7PB1amzaeYjt+/r4+n1PsvnZsbvb+MRFq/jX3z53XCe3X7SolaF8gX29WfKFImctn0PngX5a69NsfvYQucLUbP/GTJLGTJLmuhSL2hoYyhcAODiQo3sgx97esT8QTmit5+KXLmZRWz2FotOfLXDS3Eaawv+wKxc209aQpi6VpHswx65Dg+ztHQIPvhnni0V2dw9yaCDHQLZI54F+Hn7mwDFfd7wWttaRNKNnKE/PYJ76dIKiAw6NdUlOCD8YWxvSHOzPMpArsK2rj/5sYdyvkUwY85vrMIMFrfW01qcwM+pTCeY21zGQzWNmNGSS1KUS9A8V2N+fpbkuRaHoZFIJFrfV09aYwd1pyCSZ11xH31Ce3d1DDOYKHBoIPsAWtdVTl0rQPZgnVyhSKAYfoMmkMZQr0j2Yo6tn6HAANGaStDakaa5L4e4kEwmy+SKO0z2QZ0/PIN0DOXoG8xwayNGQCT7s5jRlGMwVmd2Y5qS5jSxoqWfZnEbSSSOTTNBUlyKbL9KXzZMvOK0NafqzeQ70ZcmH4dUzmCdfLDKYK7K/L0vfUJ4D/UHYDuYKh8/nFYpOMmEsaK0jYXZ4z6mpLkXPYJ5kAjLJBM31abL5I9ulb6hAfzYftg+W1T2YIxOGf77o5AtFBvNF2hqCLxGFItSng/mZVIJZDWmK4RePdNJoa0iTMKPoTr7gNGaSNGSSFIuQKxYpFp18MXhvh/JBUOcKRRbPajj8vIFsgZ7BfLjn5KQSxptetoxXnDzvuP6GdSipjHJ7oq869fj6XToepy9p4/QlbVz0ksWHp+UKRTY/282flhxqArjmjsfHXNasxjRLZzdwyvxm/m3zLr785jM566TZDOYKLJvTOOrz3J1Hdhzk7id284fdvWzb28dpJ7TQ2pAmnTDueOw5BnMFVp80m/X/o52uniH29g2RLzgLW+voPDBAJhUcmnrkmYMkE8Yfr5jDvOY6ls5uGPOE+8H+LH3ZAumksfnZbha3NdCQTpItFJjTVMfsxvS4T9i31qdprU8f82IBd+fZQ4MM5QoUncPf0JvqkgzmivQN5Xnu0AADuQL5gjO3OcPB/hwH+3PMbc4AwTfQZw8OsH1fH0O5IrMa08xrrqM/mwcgmUjQM5g7HEpP7e2jpT7F7MYMl7y0jVMWtNBan2LxrAZSCTsceulUgly+SCoZ7JHMakjTkEnSUp8e13sgMplqdo+h/dpfBN8wS2z+1BtoqquerHR3Hus8xM827eLHD3fy52ct5f2vWUlDJlnp0kRkmtMeQxnlvoxWUyhAcDnsGctmccayWVy19rRKlyMiNaJmf+BWmgsnzW08rstURURmopoNhlKXrp6yG8WJiEw7NRsMiZJjSe897+QKViIiUl1qNhhOW3TkChb1jCoickTNfiK+4uS5lS5BRKQq1WwwDDtjaVulSxARqSo1Hwzff/c5lS5BRKSq1HwwiIhIlIJBREQiFAwiIhKhYBARkQgFg4iIRCgYREQkQsEgIiIRNRsM0/w2FCIiU6Zmg2HYOG8SJiJSM2o+GEREJErBICIiEQoGERGJUDCIiEiEgkFERCIUDCIiEqFgEBGRCAWDiIhEKBhERCQi1mAwszVmtsXMtprZVWXm15nZD8L5vzaz5XHWJyIiMQaDmSWBm4C1wCpgnZmtGtHscuCAu58CfAm4Pq76REQkEOcew9nAVnff5u5Z4DbgkhFtLgG+Ew7/CLjATL0ZiYjEKRXjay0BdpSMdwJ/PFobd8+b2SFgLrC3tJGZXQFcEY72mtmW46xpXtP10WXXgHmgda4BWufa8HzW+aTRZsQZDJPG3dcD65/vcsysw93bJ6GkaUPrXBu0zrVhqtY5zkNJO4FlJeNLw2ll25hZCmgD9sVSnYiIAPEGw4PASjNbYWYZ4DJgw4g2G4C3h8N/DtzjrlvqiIjEKbZDSeE5gyuBjUAS+Ja7bzaza4AOd98A/ANwi5ltBfYThMdUet6Ho6YhrXNt0DrXhilZZ9MXchERKaVfPouISISCQUREImo2GI7VPcd0YWbLzOxeM3vczDab2QfC6XPM7Odm9ofw39nhdDOzG8P1fszMVpcs6+1h+z+Y2dtHe81qYWZJM3vEzO4Ix1eEXalsDbtWyYTTR+1qxcyuDqdvMbM3VGZNxsfMZpnZj8zsd2b2hJm9fKZvZzP76/DvepOZ3Wpm9TNtO5vZt8xsj5ltKpk2advVzM4ys9+Gz7lxXD8adveaexCc/H4SeAGQAR4FVlW6ruNcl0XA6nC4Bfg9QZcjnweuCqdfBVwfDl8I/Aww4Bzg1+H0OcC28N/Z4fDsSq/fMdb9Q8D3gTvC8duBy8Lhm4G/DIf/J3BzOHwZ8INweFW47euAFeHfRLLS6zXG+n4HeFc4nAFmzeTtTPCD16eAhpLt+46Ztp2BVwOrgU0l0yZtuwL/Fba18Llrj1lTpd+UCm2IlwMbS8avBq6udF2TtG4/BV4HbAEWhdMWAVvC4W8A60rabwnnrwO+UTI90q7aHgS/g7kbeA1wR/hHvxdIjdzGBFfCvTwcToXtbOR2L21XbQ+C3/Q8RXjByMjtNxO3M0d6QpgTbrc7gDfMxO0MLB8RDJOyXcN5vyuZHmk32qNWDyWV655jSYVqmTThrvOZwK+Bhe7+XDhrF7AwHB5t3afbe/Jl4G+AYjg+Fzjo7vlwvLT+SFcrwHBXK9NpnVcAXcA/hofPvmlmTczg7ezuO4EbgGeA5wi220PM7O08bLK265JweOT0MdVqMMw4ZtYM/DPwQXfvLp3nwVeFGXNdspldBOxx94cqXUuMUgSHG77u7mcCfQSHGA6bgdt5NkHHmiuAxUATsKaiRVVAJbZrrQbDeLrnmDbMLE0QCt9z9x+Hk3eb2aJw/iJgTzh9tHWfTu/JK4GLzWw7QS+9rwG+AsyyoCsViNY/Wlcr02mdO4FOd/91OP4jgqCYydv5tcBT7t7l7jngxwTbfiZv52GTtV13hsMjp4+pVoNhPN1zTAvhFQb/ADzh7l8smVXavcjbCc49DE9/W3h1wznAoXCXdSPwejObHX5Te304req4+9XuvtTdlxNsu3vc/S3AvQRdqcDR61yuq5UNwGXh1SwrgJUEJ+qqjrvvAnaY2QvDSRcAjzODtzPBIaRzzKwx/DsfXucZu51LTMp2Ded1m9k54Xv4tpJlja7SJ10qeLLnQoIreJ4EPlrpep7HepxLsJv5GPCb8HEhwbHVu4E/AL8A5oTtjeCGSU8CvwXaS5b1F8DW8PHOSq/bONf/fI5clfQCgv/wW4EfAnXh9PpwfGs4/wUlz/9o+F5sYRxXa1R4XV8KdITb+icEV5/M6O0MfAr4HbAJuIXgyqIZtZ2BWwnOoeQI9gwvn8ztCrSH79+TwFcZcQFDuYe6xBARkYhaPZQkIiKjUDCIiEiEgkFERCIUDCIiEqFgEBGRCAWDyDGY2XIzczOrqhvNV2tdMv0pGESObQdBZ2S/ATCz88MP5HlxFWBm95nZV8eqS2SyxHbPZ5Hpyt0LBB2ZTTozS3vQ3cOETWVdUtu0xyBVJfxm/DUz+6yZ7Q1vYHKDmSXC+dvN7CNlnvPVkvHtZvYJM/u2mfWY2Q4ze7MFN7q5zcx6w5uZvH6cNR0+ZBP2YHtvOKsrnP7tsJ2Z2d+Y2ZNmNhDeHOWtZZazzszuMbMB4D1mNteCm9B0hs/bbGbvLHnet4HzgPeFz/dwWUcdSjKzV1twk5pBM9ttZl8Ku30Z1/sbtvkzC24CM2Bm+83s381suHdPqQEKBqlGbwHywCuAK4EPAm+e4DI+SNAtwmqCG7t8h+CmPncSdC1xP/BdM6uf4HJ3AJeGwy8mOJTzgXD8WoLuDN5HcHOY64BvmNkbRyzjOuBrYZufEHTl8DBwUbjMr4TPuyBs/wHgV8A/hq+3iGgXywCY2RKCG7E8QtD9+uUE/e9fN6LpqO+vmZ1A0DHhd4AXEdxE5pZxvC8yk1S6nxA99Ch9APcBvxox7efAN8Ph7cBHyjznqyXj24FbS8abCfqTurFk2vJwWvs4aoq0JeifyYF5JW2agAHgVSOe+2XgzhHL+fA4XvO24XUut46j1PUZgr51EiVt3gEMAY3jfH9Xh8s8qdJ/C3pU7qFzDFKNHhsx/iyw4HiX4e69ZtZP0OnYsN3hvxNd7mhWEXzz/zczK+2ALE0QVKU6SkfMLElwb4U3E9xEpY7g1p33TbCGFwEPuHuxZNovw2WdwpH3ZKz391GCTts2mdld4fCP3L1rgrXINKZgkGo08mSsc+SwZ5Ggh8lS6XEuIzdiHCbvcOrwcv4bQXfRY9XSN2L8I8CHCQ4Z/RboBT7L5IUWRG/0Mur76+6F8NzLOQRdN18OXGdm57n7o5NYj1QxBYNMN10Ex9gBCM8RnEZwXD0u2fDfZMm0xwkO2Zzk7vdMcHnnAv/i7rfA4XtsnAocHPGayTLPLfUE8CYzS5TsNZwbPvfJ8Rbj7k5wTuNXZnYNsJlgb0bBUCMUDDLd3AP8hZltIAiJjxL/3/HTBN+y32hm/wIMuHuPmd0A3BB+sN9PcG7jHKDo7uvHWN7vgTeb2bkEN7B/P8HtLEvDbjtwdnhVVC+wv8xyvkZwIvlrZvYVgvsWfI7g3ET/eFYsvPnLawlu/LKb4CT2MoLgkxqhq5JkurmOIBx+CtxFcAw9zr0FPLhJ/ScJTvbuJrj5CcDHgb8lODS0meCk7qXAU8dY5LUEV1D9jCBQ+oDvjWhzA8E3/8cJAvHEUepaS/Bh/hvgWwQ3gfnfE1i9QwS3z7yD4ET2F4BPu/t3J7AMmeZ0ox4REYnQHoOIiEQoGKTmmdnN4a+hyz1urnR9InHToSSpeWa2AGgdZXa3u++Jsx6RSlMwiIhIhA4liYhIhIJBREQiFAwiIhKhYBARkYj/Dwt55LEP6uI8AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Do not learn\n",
    "# use first reward function initially and then switch to third reward function\n",
    "\n",
    "# Instantiate learner in VW\n",
    "vw = pyvw.vw(\"--cb_explore_adf -q UA --quiet --epsilon 0.2\")\n",
    "\n",
    "num_iterations_per_cost_func = 5000\n",
    "cost_functions = [get_cost, get_cost_new2]\n",
    "total_iterations = num_iterations_per_cost_func * len(cost_functions)\n",
    "\n",
    "ctr = run_simulation_multiple_cost_functions(vw, num_iterations_per_cost_func, users, times_of_day, actions, cost_functions, do_learn=False)\n",
    "\n",
    "plot_ctr(total_iterations, ctr)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Summary\n",
    "\n",
    "This tutorial aimed at showcasing a real world scenario where contextual bandit algorithms can be used. We were able to take a context and set of actions and learn what actions worked best for a given context. We saw that the learner was able to respond rapidly to changes in the world.  We showed that allowing the learner to interact with the world resulted in higher rewards than the no learning baseline.\n",
    "\n",
    "This tutorial worked with simplistic features. VW supports high dimensional sparse features, different exploration algorithms and policy evaluation approaches."
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}