section1-fig2.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import matplotlib as mpl\n",
    "mpl.rcParams['agg.path.chunksize'] = 10000"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "\n",
    "plt.rcParams['font.size']=20\n",
    "plt.rcParams['axes.linewidth']=2\n",
    "plt.rcParams['xtick.major.size']=10\n",
    "plt.rcParams['xtick.minor.size']=5\n",
    "plt.rcParams['xtick.major.width']=2\n",
    "plt.rcParams['xtick.minor.width']=1\n",
    "plt.rcParams['xtick.direction']='in'\n",
    "plt.rcParams['ytick.major.size']=10\n",
    "plt.rcParams['ytick.minor.size']=5\n",
    "plt.rcParams['ytick.major.width']=2\n",
    "plt.rcParams['ytick.minor.width']=1\n",
    "plt.rcParams['ytick.direction']='in'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "kindle1 = pd.read_csv('ratings_Kindle_Store.csv')\n",
    "kindle = kindle1.rename(columns={'A2GZ9GFZV1LWB0':'User', '1603420304':'Item', '4.0':'Rating'})"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "A13QTZ8CIMHHG4    2628\n",
       "A320TMDV6KCFU     1980\n",
       "A328S9RN3U5M68    1202\n",
       "A2WZJDFX12QXKD    1148\n",
       "A2G5IFYYHFIQNB    1133\n",
       "A5JZNVV4TEAWU      979\n",
       "A37LY77Q2YPJVL     976\n",
       "A3PTWPKPXOG8Y5     951\n",
       "A26J2TRDPZH7RZ     899\n",
       "A3KI93AT8FAXI1     891\n",
       "A3A7FF87LEVCQ1     818\n",
       "A1JLU5H1CCENWX     806\n",
       "A1E481YU3ZPK6S     756\n",
       "A3LXRIY0HE71K      727\n",
       "A2VXSQHJWZAQGY     724\n",
       "A3EPHBMU07LZ50     717\n",
       "A2JZCZYHNQHSCP     692\n",
       "A1UFU1Q1KSSWIY     684\n",
       "A3GSSX9XHONVCD     658\n",
       "AR9RKLQQC2L6K      649\n",
       "ANSX922QNYA67      643\n",
       "A2YJ8VP1SSHJ7      622\n",
       "A28REJAM6899K8     620\n",
       "AFVQZQ8PW0L        613\n",
       "ADDT2MU773IYL      612\n",
       "A1E1LEVQ9VQNK      611\n",
       "A20R37WRPLUM1D     609\n",
       "A2ZJHEB8K8JZR3     601\n",
       "A1DA6E4FNRSAWN     601\n",
       "ANOSVLTGRKABQ      600\n",
       "                  ... \n",
       "A1GQ2CYSRRTH6M       1\n",
       "A2W8728Y9PO4C1       1\n",
       "A26CLD0P9ML3J5       1\n",
       "A203MV2H9X6IH9       1\n",
       "A162C24RBF31M5       1\n",
       "A2UJQF5WYF8M9B       1\n",
       "A3T26TDX0ZEINP       1\n",
       "AWBUMDF1XO50W        1\n",
       "A3MYJ8I0WBLBD        1\n",
       "A3FS9RGMLW77MD       1\n",
       "A29XOYCEEOXYFK       1\n",
       "A1KELIVWEB8LWB       1\n",
       "A1LY1ZJTMUR6ED       1\n",
       "AIJTTVQLQQQM         1\n",
       "A3KOE4F83L6HWL       1\n",
       "AKGWNIO3APJ6F        1\n",
       "A2FYOTASUYZ34F       1\n",
       "A3KO98AQ3EIA5D       1\n",
       "A282DB3I18I7VM       1\n",
       "A3J08XKK8OIRU3       1\n",
       "A19OZRCBC7PT5J       1\n",
       "A2SBQT7MQFVP50       1\n",
       "A3HZLQG28MEDIP       1\n",
       "A142G3IVHI8IGG       1\n",
       "A2SWYN6LTCERHM       1\n",
       "A2MDC2RKAAXA4O       1\n",
       "A34HN8NK0KQEDN       1\n",
       "A2DNM5SL97RBP9       1\n",
       "A1BY1RNCSAUJX        1\n",
       "ARLXSUMR9LQ8I        1\n",
       "Name: User, Length: 1406889, dtype: int64"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "\n",
    "rating_list = kindle['User'].value_counts()\n",
    "max_rating = np.amax(rating_list)\n",
    "# rating_list.remove(max_rating)\n",
    "sorted_list = rating_list.sort_values(ascending=False)\n",
    "sorted_list"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 186,
   "metadata": {},
   "outputs": [],
   "source": [
    "#Rating Chart for the Most Active User\n",
    "most_rated = kindle[ kindle['User']== 'A2G5IFYYHFIQNB'] \n",
    "user_sort = most_rated.groupby('Item')\n",
    "rating = user_sort.Rating.mean()\n",
    "height1, bin_edge1 = np.histogram(rating, bins=bins)\n",
    "\n",
    "#Rating Chart for Another Active User\n",
    "sec_most_rated = kindle[ kindle['User']== 'A5JZNVV4TEAWU'] \n",
    "user_sort      = sec_most_rated.groupby('Item')\n",
    "rating         = user_sort.Rating.mean()\n",
    "height2, bin_edge2 = np.histogram(rating, bins=bins)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 183,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA+gAAAIHCAYAAAAMxnYHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nOzdedwdZXn4/89FCCBgAJXFL1skCiiVTWQRlACKiFUwYqtWtq/+LIIFFVsX8EvA0qJSFygUl2JArVhbRKmIIJsigggooigUCIuAS5CAISCQ6/fHzJCTw9nmyXmeZ/Kcz/v1mtdkZu65z3UmydznmuW+IzORJEmSJEmTa6XJDkCSJEmSJJmgS5IkSZLUCCbokiRJkiQ1gAm6JEmSJEkNYIIuSZIkSVIDmKBLkiRJktQAJuiSGici5kdERsTsyY5FkjR8ETGzPM873u+Q2YZKKzYTdI2kiJhX/TCIiMcjYr0+5fdvKZ8RccgEhdotnkMiYm5EbDvEOrdp+X63DKvets/Ytoz7kPGofzJExOXlMZs3zLKS1K6tLbposuPpJiJml+f6/Sc7loliGzo2tqHS05mgS7Ay8NY+ZQ6aiEBqOAQ4Dhhagl7WWXlBROwyxLor21LEfUifcrcBvwYeGYcYJGlFdXDLn/eKiI0mLZLeZlOc60cmQcc2VNKQmKBr1N1Vzrsm4BHxLOC1wJ+AByYiqIkWEa0XKb5Qzg/uUnzcZeZembllZv54smKQpCaJiGdTtEWPAP9B8RvubZMalADbUEnDZYKuUfcjiivN20XEVl3KvBlYBfhvYPFEBTbBXgOsR3E8/hFI4K8jYrVJjUqSVHkrMB34JvDZct2kJYFahm2opKExQZfgS+W82130av3Z/SqKiFUj4n0RcU1ELIyIxRHx64j4ZERs0GO/bSLi7LJjl8ci4uGIuD0iLoyI90TE6mW5Q8oOdXYvd/1i27vx8wf7yk9T/cj7SmbeCfwQWBt4/QDfeeOI+JeIuKmM++GI+GVE/HtE7NFSLoEvlou7t8W9TGc2nTq4iYgvlOv+q088HyrLXd9l++si4psRcX9E/DkifhcR50fEq/t91/EUEetFxCfK47goIh6NiLsj4qqIOCEiNu2y37oR8c8R8fOI+FO5700RcWL59EenfZ46vhGxYUScXv57eywifjq+31TSGD11ngZ+QPEE2JYRsWO3Hcr3lZ96ZzciDi7bp4cj4qGIuCwiXtXrQyNi/fIc/6uIeKRs234cEUdHxKptZWeW5/rjqpg7nOtndvmcv4iIc8pz86Pl530kIlbpE99u5X73lOewBRHxvYh4S0REh/KzW9vLiHhNRHynbAuWRMR7en1eF7ahtqG2oRqezHRyGrkJmEdxhfscYLPyz/cAK7WV27zcdhfFBa17yuVDOtS5LnB9uT2BR4GHWpYfAHbusN++wJ/b9lvYspzAlmXZvwbubym/sFyupmvHcCyeBTwGPA6sW647rKz/2332fSPF45ZVnIuBh1uW57eUvb/le/25Le77gZe1lJ1flpvdsm7Pls+Y0SOmG8ty729bPx34cttxbT/OHx/D8bu83HfeWMsCmwL3tsTxRPnvZUnLusM61LcbsKClzGNtfx93AVt02K86vu8Efl/+eRHFaxw/nez/n05OTstOwFbl/9M/ANPLdSeV607rsd/c6pxD8eh1dX5pPfc9Cbyxy/47tp1jHirPwdXyT4H1WspvXJ7P/9Ryvm4/129clp3ZUs/eLeeuB8uYqm3n9fh+H2s7hz/Utu9XeXq7PrvcNh84uvzzEuCP5bF5T82/G9tQ21DbUKehTpMegJPTZEy0JOjl8pXl8ivbylWPqv1zudwrQf8OSxPxNwHTyvU7tDR49wPPadvvtnLb+cDmLetnAC8HPgfMbNvn8m5xjOFYHF7WdUHLumeXPwCeADbost8uFD9IErgUeCkQ5bZ1KToHOrNtn0PK8pf3ialq/Ga3rFsJ+E25/uAu+1U/YpcAG7Vt+1S57Q7gLcCa5fo1y0a2+qHxlprHr/q7mDfWssCZ5fpby7/zlcr1qwJ/AXwU2L9tn00pflAm8Hlgi/IYRXkcqn+Pv6j+LXY4vg+X/zZbf9g9f7L/fzo5OS07AR8v/8+e3rLuxeW6BcAqXfabW5b5I0VidhiwerntecAV5fZ7gZXb9l2HpUnPjcBLy/XTgAMo2roELu7xufN6fKeZLE2E/gh8jbKtA9YAPsjSBGvfDvsfVW77HfAuYO1y/WoUbXAV+4fa9pvN0kT1CeA0YP2WfTfqFnOX72EbahtqG+o01GnSA3BymoyJpyfo7yyXz24pEy0n4ReW6zom6GWDUP3Q2KfD563P0h8zJ7SsX69lv/VrxH95pzjGeCyuKet6W9v6b5Xrj+6z3xWUd3QG+Kwx/7go13+yXP/dLvv9Y6f6gRdQ3FX5I7BZl33/qtz3pprHbxg/Ln5Zrv/rGp9b3cn4TJftq1Dc3UrggC7H9491/t05OTlN/ESREFfJ5m5t26qLv93ugM9taWP+psP251LcNUzgFW3bPtJynnhakklx17uqe88unzuvx/ea2bL/RZTJaVuZ88vt7Ynq2hTJ0ePAjl3q35ki0XyAlgsYLE3QE/iPIfz92IambaiT0zAn30GXCv9J8Wj5nIhYo1y3O8UV1p9k5s199j+gnP8kMy9s35iZvwXOKBf/qmXTwxQ/IKD4oTShIuKFFI8wPgKc17b5K+X84A77bVnuB/APmfn4uAW5rP8o53tF57Hr39xWrnIQxZXx8zLz9i51n0vxQ3WriJjov4uHyvlAnxsRz6C4QwTFD66nycw/A9W7ht3eMT27/Lcpqbn2pjg3VO82t+p6nm5zF08/L5KZ9wFVT99/0ba5ate+kJn3d9j3IopO0WDZdm0sTsrM7LC+apfaY3sjxZ3bK7NLT+WZeTVwO8WTAC/p8rmfGEOsT7ENXYZtqDQkJugSkJkPUlypX4Oi4YcancMB25fzy3qUubScb15dBMjMxRRXzwG+GxHHRsS2ETFt4OCXT/XD4VuZ+ae2bd+iuIDw4ohoH29953L+QGZeM54BtsrMnwC3UNxRWuYHYUTsBMyiuKPS3gnOy8r5AWXHNk+bKJ6OmF6W23jcvkRnF5Tzj0XEaRGxR/kDopsdKK7uA1zT4zv9fVmm2/f5UZf1kpqjOk9/tUMS+1WKO3mviYh1e9Txky4JMBSPPUORyAJQdsxWJcWDtGvb9ygziGu7rH9abKXqnL5Tt/NfeQ7cpCzX6Ry4GPjZ8oVtG2obKg2fCbq0VJWIH1ie2N9I0VB9dYB9qx9Gv+lR5p5yHsBzWta/A7iZ4nH3jwI3AA9GxLcj4m1RjK86dBHROoZupzsri4FvlIvtdwDWL+d3MfGqv4+3tK2vli/MzPbx6qur6mtSxN5tqs6Jqw8z4AF8jOLH3CoU7zNeCjxU9j779xGxdlv51rsEvb7PjLJMt+/z++GEL2k8RMRawH7lYqfz9F0UPbq3jsPdycM9tj1azqe3rHsWS8+Hg7RrvS4O9JWZ3eLrFBssPQc+g97nwGq/TufABZm5pMP6gdiG2oZiG6pxYoIuLXUhRWczewLvpjgxfycz/1CjjlX7F1lW+bjY1sAbKDqEu5miEdyXYgi4ayJizbr1DuCVwIbln7/VYciWZOlTBG9tu1DwtKFrJlD12OAuUQ7XU/5Qqu4GPO2HEkvPdUdlZgwwXV4jnuoHZK+r9ZWqkV/cujIzH8vM/Sg6Dfo4cDXFXbFq+ZaI2KbD9/njgN9ndpd4nhwgZkmT568pOi4DuLHLefoV5fZ+j7mPVe12bQJU58BPDXgOnNehjuU9/9mG2obahmpcmKBLpcx8gmLYtZWAE8vVX+q+xzKqq6ib9iizUfVRFEPlLPPZmXleZv5tZr6I4uru31M0XNuzdEzZYarzY2494DUty9X7iJt0KDuuMvNW4DqKHzjV+3KzKY7ZIoqr6O2qd8ReNA4hLSjng7z7VpVZ0GljZl6dmR/IzF0oHul8C8UdlnUphkiqVN9nnYjYoH7IklYQdc7T20XEi4f0udUQVTBYuzbRdxLH85w+KNvQ4bANldqYoEvLqh5zn07RO+f5A+53fTnfPSK6XRnfs5zfkpmLelWWmfdn5snAp6t624pUP5zGdBU+ImZQ3LEHeB1FQ9ZtOq0s1/pj5Opy/qyI2JnBLVfcLaor/NUjndWjeedl5iMdylfvib0uItoflVxeN5Tz7SJitW6FImJjlv6YvaFbuUpmLsrMcyhGGAB4SUsHhj+hGL4HYE79kCU1XUQ8n6Xv/m5L7/N01VYN5S562UHWTeXiHj2KVu3a9W3rh3Wu76Y6p+8eEc8ep8/oyjZ0qGxDpTYm6FKLzLyOYniYfwHek5mPDbhr1aHKVix9X/ApEbE+xfizUPQYX62f3iOhh6WPcbU/Ylj1WNr+XtWg3kTxONkCivfNHuw2UYxNC0XD/CyAzPwVS3v+/XiNBnt5466cQ/FD5cURsR1LO/br9GgewFll+f8DfKhXxRHR3hlRP98o614T+Nse5Y4u5w8DF7d95ipPL/6U6t9AUHZqU76v+d/l+mPLf18dRcTK4/SKhKTxVSV0P8vMn/U5T3+9LPs3Q+xktGrXDunUK3dE7E3xCDG0tGulYZ3ru/k6xd3e1ejTE/sYzumDsA3twjZUWn4m6FKbzDw+M9+fmYP03l7t8wOKd9gBzoyIA6ofSRHxEooxXteheKzqMy27bgXcFBHviYjNq2S9TNzfCLyvLPfdto/8RTmfU3YiVFdrz7NP9CxZDOvzW4qG7c0t699HcQX65cCFEbFDtSEinhMRb46Ir7CsKu4XlT3Gjklm3svS3u+/QHFs/0BxnDuVv5mlTyMcX/byullLvGtGxKsi4kss/aE7aCy3AWeWiydHxIejpTfliNg4Ij4JHFWuOikzH2qr5qaI+KeIeGn1QyMKOwKnlmWuzcw/tuzzQYrHUJ8LXBURb4iIpy7kRMTzI+I9FH0a7ICkFUbZFhxYLp47wC7nU3RqugHw6iGF8a/AfRSJ6FPn+IiYVrZP55TlvpeZl7btW53rd4uIFwwpnqdk5gKWJoqHRsR/RsRTQ7FFxGoRsVtEnMbTh6YbBttQ21Bp/GQDBmN3cproCZhH8S74OTX3u6fc75AO29aleOwqy2kxxdXuavkBYJe2fbZt2Z4U75wvoOh4pFp3LTCjbb8tKcYbTYofZb8B5lOMCdvvO2xGcbU6gdcO+L3PKMtf07b+zWXMVayPUFzdrpbnd6jripbtC8q45wM7t5SZX26f3SOmd7Qdu9P7fIdpwOlt+zxE8SrDkpZ1l43h39MzKMbAba37wba//wT+DVipw/4PtpR5ojwuf25Z93tg6w77vbT8u6/KPU7xI+vRts/dvW2/vsfXyclp8iaKx8qr/79bDbjPhWX5r7Wsm1uum9djv3llmbkdtu1I0Xa1njMXtyz/DFivw37Tgf8tyyyh6IC1OtdvVJaZWdXTI7bZ3dqScvuxbefvRWW8rW3oHXXqHOA424bahvY9vk5OyzN5B10aksz8PcXjfkdTvN/0OMUV81sprjxvlZntY2beDBxA0XjfQNHIzKBolK4E/g7YNduuFmfxeNyrKH6QLaS4a7IpS9/P6uUgike9Hga+N+DXqx4F2zEitmyJ4xzghRR3Wm4pVy8pv9cXWNqDbas5FI38HRSPtG1aTl3fPevivyguUlS6PZpXxfpkZh4O7AZ8GbiT4u/nGRSdyHyD4q7I/jXjIIvhdN5A8XrDf1NcyFmNYuijOyh6zd09M9+VnYf12Q/4Z4o7LfdSHJc/AzcCJ1H827mxw+deS3Gx5gPAVRR/p2tT/ID+CcXQMy/NzCva95XUaAeX81sy8xc9Sy5Vnaf3i6cPKzUmmfljio7BPkVxjp9OkQD9hKIj050y83cd9nsc2Iuio9XfUNyhrc71Qxs6NDP/EdiGYgSUWynatjUo7vx/B3gXMOY7zV3YhtqGSuMqMnOyY5AkSZIkaeR5B12SJEmSpAYwQZckSZIkqQFM0CVJkiRJaoChddSxoooIX8KXJEmSpCkuM2OyY+jHO+iSJEmSJDXAyN9Br9ibvSRpVEQUNxBs+yRJo6Bq91YE3kGXJEmSJKkBTNAlSZIkSWoAE3RJkiRJkhrABF2SJEmSpAYwQZckSZIkqQFM0CVJkiRJagATdEmSJEmSGsAEXZIkSZKkBjBBlyRJkiSpAUzQJUmSJElqABN0SZIkSZIawARdkiRJkqQGMEGXJEmSJKkBTNAlSZIkSWoAE3RJkiRJkhrABF2SJEmSpAYwQZckSZIkqQFM0CVJkiRJagATdEmSJEmSGsAEXZIkSZKkBjBBlyRJkiSpAUzQJUmSJElqABN0SZIkSZIawARdkiRJkqQGMEGXJEmSJKkBTNAlSZIkSWqAlSc7AEmSJtK5rzt3skMYkznnz5nsECRJ0jjzDrokSZIkSQ1ggi5JkiRJUgOYoEuSJEmS1AAm6JIkSZIkNYAJuiRJkiRJDdCYBD0iNoqIMyPi3oh4LCLmR8SnI2KdGnXMj4jsMt0/nvFLkiRJklY8TcpFGzHMWkTMAq4C1gO+CfwK2BE4CtgnInbNzAUDVrcQ+HSH9X8aRqySJEmSpKmhabloIxJ04HSKA3JkZp5arYyITwLvBU4EDhuwrgczc+7QI5QkSZIkTTWNykUn/RH3iNgM2BuYD5zWtvk4YBFwYESsMcGhSZIkSZKmqCbmok24g75nOb8oM5e0bsjMhyPihxQHbWfgkgHqWzUi3gZsQnFAbwS+n5lPDjFmSZIkSdKKrXG5aBMS9C3K+S1dtt9KcVA2Z7CDsgHwpbZ1d0TEoZl5RbedIqJvxZk5wMdLkjR5BmnP6pS17ZMkNVmddq+DRuSirSb9EXdgrXK+sMv2av3aA9T1RWAvigOzBvBi4LPATOA7EbHN2MOUJEmSJE0hjctFm3AHvZ/qkkjfS/iZeXzbqpuAwyLiT8DRwFzgDV32XY4QJUlqhkHas+pug22fJGlFV6fdG4MJyUVbNeEOenVVYq0u22e0lRuLM8r5K5ajDkmSJEnS1NG4XLQJCfqvy/nmXba/oJx3ey9gEL8r5/YEL0mSJEmCBuaiTUjQLyvne0fEMvFExDOBXYHFwNXL8Rm7lPPbl6MOSZIkSdLU0bhcdNIT9My8DbiI4uX5I9o2H09xpeHszFwEEBHTI2LLiJjVWjAitoqIZ7XXHxGbAv9aLn55yOFLkiRJklZATcxFm9JJ3OHAVcApEbEXcDOwE7AHxeMEx7SU3bDcfifFgay8CfhgRFwG3AE8DMwCXgusBlwAnDyu30KSJEmStCJpVC7aiAQ9M2+LiB2AE4B9gH2B+4BTgOMz84EBqrmMYhy77SgeI1gDeBC4kmIsui+l3dVKkiRJkkpNy0Vj1HPWiEhwqBlJGhXnvu7cyQ5hTOacP2dodTnMmqSmWBHPycM8H2titLR7Yx5vbaJM+jvokiRJkiTJBF2SJEmSpEYwQZckSZIkqQFM0CVJkiRJagATdEmSJEmSGsAEXZIkSZKkBjBBlyRJkiSpAUzQJUmSJElqABN0SZIkSZIawARdkiRJkqQGMEGXJEmSJKkBTNAlSZIkSWoAE3RJkiRJkhrABF2SJEmSpAYwQZckSZIkqQFM0CVJkiRJagATdEmSJEmSGsAEXZIkSZKkBjBBlyRJkiSpAUzQJUmSJElqABN0SZIkSZIawARdkiRJkqQGMEGXJEmSJKkBTNAlSZIkSWoAE3RJkiRJkhpg5ckOQJJUOPd15052CLXNOX/OZIcgSZI0ZXgHXZIkSZKkBjBBlyRJkiSpAUzQJUmSJElqABN0SZIkSZIawARdkiRJkqQGMEGXJEmSJKkBTNAlSZIkSWoAE3RJkiRJkhrABF2SJEmSpAYwQZckSZIkqQFM0CVJkiRJagATdEmSJEmSGsAEXZIkSZKkBjBBlyRJkiSpAUzQJUmSJElqABN0SZIkSZIawARdkiRJkqQGMEGXJEmSJKkBTNAlSZIkSWoAE3RJkiRJkhrABF2SJEmSpAYwQZckSZIkqQFM0CVJkiRJagATdEmSJEmSGsAEXZIkSZKkBlh50IIRMQ1YNTMfaVu/J7Af8Ajwucy8Y7ghSpIkSZI09dW5g34y8EBErFWtiIg3AxcDfwd8APhxRGw83BAlSZIkSZr66iTorwAuy8yFLeuOAx4EDgL+AVgbeN/wwpMkSZIkaTTUSdA3Bv63WoiIzYAtgFMz88uZeTLwHWCf4YYoSZIkSdLUVydBnwE81LK8K5DAhS3rfgFsNIS4JEmSJEkaKXUS9PuA57UsvxJYDFzXsm5N4IkhxCVJkiRJ0kgZuBd34Grg9RHxl8CjwAHAJZn5eEuZzYDfDDE+SZIkSZJGQp076P9Ulv8m8F1gFeDEamNEzABmA9cMMT5JkiRJkkbCwHfQM/PnEbETcHC56muZeW1Lka2Bi4CvDjE+SZIkSZJGQp1H3MnMnwPv77LtSuDKYQQlSZIkSdKoqfOIuyRJkiRJGicD30GPiIMGKLaEYii2mzPz1jFHJUmSJEnSiKnziPs8inHPBxIRvwCOyMwf1A1KkiRJkqRRU+cR90OBbwEBXAIcDxxezi8t138T+AfgHOCFwHcjYptBKo+IjSLizIi4NyIei4j5EfHpiFinRoztdR4YEVlO7xhrPZIkSZKkqalJuWidO+i/B14DvCYzv9shgH2A84DPZ+a/RMSZFL26fwB4a5/gZwFXAetRJPm/AnYEjgL2iYhdM3NBjViJiI2BU4E/AWvW2VeSJEmSNPU1LRetcwf9GODcTsk5QGZeCJwLHFsuXwJcDOw+QN2nUxyQIzNz/8z8YGbuCXwK2IKW8dYHEREBfBFYAJxRZ19JkiRJ0shoVC5aJ0HfBri9T5nbKcZDr/wCeE6vHSJiM2BvYD5wWtvm44BFwIERsUaNWI8E9qR4LH9Rjf0kSZIkSSOgiblonQT9z8CL+5TZGni8ZXll4JE+++xZzi/KzCWtGzLzYeCHwOrAzoMEGREvBE4CPpOZ3x9kH0mSJEnSyGlcLlonQb8M+MuIeGeXYA4D/pKiw7jKlsA9ferdopzf0mV7NVzb5v0CjIiVgS8BdwEf7le+bd++kyRJTVenPbPtkySt6JazLWtELtqqTidxHwRmA/8WEUcDPwJ+C6wP7AI8H3iwCiYi1gf2AP6tT71rlfOFXbZX69ceIMb/B2wH7JaZiwcoL0mSJEkaTY3LRQdO0DPz1oh4GcVL9HsAL2grcjnFuOfV1YffAc9k2Ufex6K65NFzDPaI2JHi4sC/ZOaP6n5I5sBDvEuS1FiDtGfV3QTbPknSiq5OuzcGE5KLtqpzB53M/DWwV0RsBGxLccXhIeCGzLynrWwCjw1QbXVVYq0u22e0lXualscJbgE+MsBnSpIkSZJGW+Ny0VoJeqVMxvu9Wz6oX5fzbs/1V3fqu70XAMXYctX+j3a5QvL5iPg8xQv776kdpSRJkiRpKmlcLjqmBH3ILivne0fESq2950XEM4FdgcXA1T3qeAz49y7btqd4F+BKir+A5XrkQJIkSZI0JTQuF62VoEfEs4D/C+wIrANM61AsM3OvQevMzNsi4iKK8eeOAE5t2Xw8sAbw2cxcVMYwHZgFPJ6Zt5V1LAbe0SXmuRQH5azM/MKgcUmSJEmSpq4m5qIDJ+gRsSVFR3DrsvRl+U7G0uPM4cBVwCkRsRdwM7ATRWd0twDHtJTdsNx+JzBzDJ8lSZIkSRI0LBetMw76ycB6wMeAzYDpmblSh6nTXfWeyqsPOwDzKA7G0RRXJk4BdsnMBXXrlCRJkiSpl6blonUecX858O3MHPOg671k5t3AoQOUm0/vO/jt5ecCc8calyRJkiRp6mpSLlrnDnoAv6xTuSRJkiRJGkydBP06YIvxCkSSJEmSpFFWJ0E/Adg3ImaPUyySJEmSJI2sOu+gbwx8E7goIr5KcUf9wU4FM/PsIcQmSZIkSdLIqJOgz6MYQi2AA8upfUi1KNeZoEuSJEmSVEOdBL1vr3aSJEmSJGlsBk7QM/Os8QxEkiRJkqRRVqeTOEmSJEmSNE5M0CVJkiRJaoCuj7hHxO0UHb69MjPvKJcHkZk5ayjRSZIkSZI0Inq9g74Sy/bS3r7cTSxXRJIkSZIkjaCuCXpmzuy1LEmSJEmShsd30CVJkiRJaoCBE/SIuDQiDupT5m0RcenyhyVJkiRJ0mipcwd9NjCzT5lNgd3HGowkSZIkSaNq2I+4PwN4Ysh1SpIkSZI05fXqxb2Tjr24R0QAmwD7Ancvb1CSJEmSJI2annfQI2JJRDwZEU+Wq+ZWy60TxV3z24FtgXPGOWZJkiRJkqacfnfQv8/Su+avAO4C5nco9ySwALgE+MKwgpMkSZIkaVT0TNAzc3b154hYAnwxM08Y76AkSZIkSRo1dd5Bfx7w4HgFIkmSJEnSKBs4Qc/MO8czEEmSJEmSRlndXtyJiOcCewEbAqt2KJKZ+dHlDUySJEmSpFFSK0GPiOOBD7btFyztSK76swm6JEmSJEk19BxmrVVE/A3wEeAHwAEUyfhZwFuBzwNLKIZY23P4YUqSJEmSNLXVuYP+LuAeYJ/MfCIiAOZn5jnAORHxDeDbwFeHH6YkSZIkSVPbwHfQgRcDF2TmEy3rplV/yMzvAt8F/n5IsUmSJEmSNDLqJOjTgQUty4uBtdrK3ARss7xBSZIkSZI0auok6PcBz21ZvgvYuq3MhsATSJIkSZKkWuok6DdQPOZeuRR4eUQcGBFrRMRrgTeW5SRJkiRJUg11EvT/AbaKiOeVyycBC4F5wEPAtyh6dj92mAFKkiRJkjQKBu7FPTPnUSTj1fLdEfFS4GhgFjAfOD0zfz7cECVJkiRJmvrqDLP2NJl5B/DuIcUiSZIkSdLIqvOIe18RsWVEfH2YdUqSJEmSNAqW6w56JSI2BeYCb2PISb8kSZIkSaOgbzIdEbtExCUR8VBELIiI8yLi+eW21SLiZODXwMHA74CjxjdkSZIkSZKmnp530CNia+ASYLWW1a8Hto+IlwHnA9sAvwc+RtFJ3KPjFKskSZIkSVNWvzvo/5ib5+8AACAASURBVECRnH8W2BHYCfh3YCPgB8DWwMnArMz8pMm5JEmSJElj0+8d9N2AazLzXS3rro2IbYCXAMdk5knjFp0kSZIkSSOi3x30DYAfdlj/g3L+78MNR5IkSZKk0dQvQV8FeKjD+ocAMvP3Q49IkiRJkqQR5JBokiRJkiQ1wCDjoO8fETPb1m0LEBFndiifmfn25YxLkiRJkqSRMkiCvm05dXJIh3UJmKBLkiRJklRDvwT90AmJQpIkSZKkEdczQc/MsyYqEEmSJEmSRpmdxEmSJEmS1AAm6JIkSZIkNYAJuiRJkiRJDWCCLkmSJElSA5igS5IkSZLUACbokiRJkiQ1gAm6JEmSJEkNYIIuSZIkSVIDrNxtQ0QcNNZKM/Psse4rSZIkSdIo6pqgA/OArFlflPuYoEuSJEmSVEOvBP3QCYtCkiRJkqQR1zVBz8yzJjIQSZIkSZJGmZ3ESZIkSZLUACbokiRJkiQ1QK930J8mItYADgdeDWwIrNqhWGbmrCHEJkmSJEnSyBg4QY+ItYErgRcBDwEzgIXAKsAzymL3Ao8POUZJkiRJkqa8Oo+4H0uRnL8dWKdc9ylgTeBlwPXAbcALhxmgJEmSJEmjoE6C/nrg+5n5xcx8anz0LFwN7AtsCRwz5BglSZIkSZry6iToG1PcJa8soeUd9Mz8HfAd4M1jCSQiNoqIMyPi3oh4LCLmR8SnI2Kd/ns/VcfHIuKSiLg7IhZHxAMRcUNEHBcRzx5LXJIkSZKkqatJuWidBP0R4MmW5YXABm1lfkvReVwtETELuA44FPgxxaPztwNHAT+q8YXeC6wBXAx8BvgK8AQwF7gxIjauG5skSZIkaWpqWi5apxf3uynuold+CbwiIqZlZpW47wbcX6POyunAesCRmXlqtTIiPknxRU8EDhugnhmZ+Wj7yog4Efgw8CGKXuglSZIkSWpULlrnDvoVwO4REeXy14BZwLcj4oiI+DqwM3BBjTqJiM2AvYH5wGltm48DFgEHlkO89dTpgJT+s5y/oE5skiRJkqSpqYm5aJ0E/SzgPGCjcvmMcnlv4FTgjcBVFL2917FnOb8oM5e0bsjMh4EfAqtTJP9j9bpyfuNy1CFJkiRJmjoal4sO/Ih7Zl4PvKtl+QlgTkS8BHg+xVWHa9u/2AC2KOe3dNl+K8VFgM2BSwapMCLeTzH821rADhSP3t8InNRjn771tnReL0lSIw3SntUpa9snSWqyOu1eB43IRVvVeQe9o8y8juKl+rFaq5wv7LK9Wr92jTrfD6zfsnwhcEhm/r5mbJIkSZKkqalxuejAj7hHxO0RcWSfMkdExO2D1jnoR5fzgS/hZ+YGmRkUvczPATYDboiI7Xvs03eSJKnp6rRntn2SpBXdOLdlE5KLtqrzDvpM+l85WBvYtEadsPSqxFpdts9oKzewzPxtZn6D4rGEZwNn161DkiRJkjQlNS4XrZOgD2JN4M819/l1Od+8y/aqt7tu7wX0lZl3UgwLt1VEPGes9UiSJEmSpozG5aI930GPiE3aVq3dYR3ANGAT4ACKQd3ruKyc7x0RK7V2MhcRzwR2BRYDV9est93/KedP9iwlSZIkSRoFjctF+91Bnw/cUU4AR7Ust07/C1xKMS765+tEmpm3ARdRPEJ/RNvm44E1gLMzcxFAREyPiC0jYlZrwXLdBu31R8RK5eDw6wFXZeYf68QnSZIkSZp6mpiL9uvF/WyKF+IDOIiie/ifdij3JLAAuCQzL+r3oR0cTjGG+ikRsRdwM7ATsAfF4wTHtJTdsNx+J8WBrOwDfCIivg/cVsazPrA7xYv59wP/3xhikyRJkiRNTY3KRXsm6Jl5SPXniDgI+EZmnjBIxXVk5m0RsQNwAsWX2xe4DzgFOD4zHxigmu8Bn6N4DGEbig7rFlEc1C8BpwxYjyRJkiRpBDQtFx14HPTMHHaHcu313w0cOkC5+Szt7r51/U08/bEESZIkSZK6alIuOnCC3ioiNgK2o7gysBC4PjPvGUZAkiRJkiSNoloJetmD++eAV3XYdjFwWHlVQZIkSVphnfu6cyc7hNrmnD9nskOQtJwGTtDLXul+SPFi/Hzg+xTP5j8X2I1iAPYrI2KHzLx/+KFKkiRJkjR11bmD/hGK5PwDwCcz86kx3CJiGvBe4OPAscC7hxmkJEmSJElTXZ2O314LXJSZn2hNzgEy88nMPJliDLm/HGaAkiRJkiSNgjoJ+gbAdX3KXFeWkyRJkiRJNdRJ0BcCm/Yps0lZTpIkSZIk1VAnQb8SOCAiXtZpY0TsBLypLCdJkiRJkmro2UlcRJwJnJeZ3wL+ieI99Csi4hzgMope3DcAZgNvAZaU5SRJkiRJUg39enE/hGJItW9l5nURcQAwD/gb4K0t5QJ4APi/mdnvPXVJkiRJktSmzjBrZOb/RMSmwH7A9sBaFO+c30Bxp33R8EOUJEmSJGnqq5WgA5RJ+H+UkyRJkiRJGoI6ncRJkiRJkqRxMsgd9G0j4qA6lWbm2WOMR5IkSZKkkTRIgr5fOdVhgi5JkiRJUg2DJOg/KydJkiRJkjROBknQz8vME8Y9EkmSJEmSRpidxEmSJEmS1AAm6JIkSZIkNYAJuiRJkiRJDdAvQb8CmD8BcUiSJEmSNNJ6dhKXmXtMVCCSJEmSJI0yH3GXJEmSJKkBTNAlSZIkSWoAE3RJkiRJkhrABF2SJEmSpAbomqBHxNYRsd5EBiNJkiRJ0qjqdQf9BuCwaiEiLo2Ig8Y/JEmSJEmSRk+vBH0JMK1leTYwczyDkSRJkiRpVPVK0O8Btp2oQCRJkiRJGmUr99h2PvDuiLgZuK9cd0hEzO5TZ2bmXsMITpIkSZKkUdErQT8GWAV4LbA7kBSPuM/sU2cOIzBJkiRJkkZJ10fcM/PhzDwsMzfOzGlAAHMzc6U+07RudUqSJEmSpM7qjIN+BTB/nOKQJEmSJGmk9XrEfRmZucd4BiJJkiRJ0igbOEGvRMTqwBxgO2BtYCFwPfCNzFw03PAkSZIkSRoNtRL0iNgXOAt4FsU76ZUEPhURh2bm/wwxPkmSJEmSRsLACXpEbA+cC0wDvgJcSjH82nOBPYG3AP8VEbtm5nXjEKskSZIkSVNWnTvox1DcKX95Zl7dtm1eRJwGXA58GHjjcMKTJEmSJGk01OnF/eXA1zsk5wBk5jXAf5XlJEmSJElSDXUS9LWAu/uUuQuYMfZwJEmSJEkaTXUS9HuBHfuU2YHivXRJkiRJklRDnQT9AmDPiPhgRExr3RARK0XE0cAry3KSJEmSJKmGOp3EfRTYHzgR+NuI+AHF3fINgN2AmcD9wD8OOUZJkiRJkqa8gRP0zLw/InYFPgu8Cti0rcjFwGGZ6SPukiRJkiTVVOcOOpk5H3h1RGwIbEfRcdxC4IbM/M3ww5MkSZIkaTTUStArZTJuQi5JkiRJ0pDU6SROkiRJkiSNExN0SZIkSZIawARdkiRJkqQGMEGXJEmSJKkBTNAlSZIkSWoAE3RJkiRJkhrABF2SJEmSpAaoPQ56RKwH7ACsA0zrVCYzz17OuCRJkiRJQ3Du686d7BBqm3P+nMkOYVIMnKBHxHTgDOAgut95DyABE3RJkiRJkmqocwf9o8ChwG3AV4C7gSfGIyhJkiRJkkZNnQT9rcAtwHaZuXic4pEkSZIkaSTV6SRuPeACk3NJkiRJkoavToJ+FzBjvAKRJEmSJGmU1UnQ5wGviYi1xikWSZIkSZJGVp0E/STgSuB7EbFHRHg3XZIkSZKkIanTSdzj5TyA7wFERKdymZm1x1eXJEmSJGmU1Umkf0AxxrkkSZIkSRqygRP0zJw9jnEQERsBJwD7AM8G7gPOA47PzD8OsP+zgTcArwVeDGwI/Bn4OfBF4IuZuWR8opckSZIkrYialIs24lH0iJgFXEUxlNs3gV8BOwJHAftExK6ZuaBPNW8C/o3iYF5G0ev8+sAc4AsUHdy9KTN9CkCSJEmS1LhcdEwJekRMB7YE1gYWAjdn5uO99+rpdIoDcmRmntryOZ8E3gucCBzWp45bgNcD3269OhERHwZ+DLyR4gD993LEKUmSJEmaOhqVi9bpxZ2ImBERZwAPAj8FLgduAB6MiDMiYu069ZV1bgbsDcwHTmvbfBywCDgwItboVU9mXpqZ57c/OpCZ9wNnlIuz68YnSZIkSZp6mpiLDpygl8Oq/RB4J/AERadx/1nOHy/XXzmG4df2LOcXdfhCD5efuTqwc816W1V3959YjjokSZIkSVNH43LROnfQPwRsRfFs/aaZOTsz31J2HrcpxRWHF5Xl6tiinN/SZfut5XzzmvUCEBErAweVixf2KNd3kiSp6eq0Z7Z9kqQV3XK2ZY3IRVvVSdDnAFdn5hGZ+WDrhsxcmJl/B/yI4vn6OtYq5wu7bK/W1358vnQS8BfABZn53THWIUmSJEmaWhqXi9ZJ0DeheOe8lyuAjWvUOYjqkkft3tcj4kjgaIqe+A7sVTYz+06SJDVdnfbMtk+StKIb57ZsQnLRVnUS9EcoerfrZd2yXB3VVYm1umyf0VZuIBFxBPAZ4JfAHpn5QM24JEmSJElTV+Ny0ToJ+rXAmyLiBV2CmAX8VVmujl+X827P9Vef1+29gE6xvAf4V+AmigNyf82YJEmSJElTW+Ny0TrjoH8CuAi4NiJOpRiA/T5gA4ou4/8OWBM4uU4AZT0Ae0fESm3jxj0T2BVYDFw9SGUR8QGKZ/1/CrwqM/9QMx5JkiRJ0tTXuFx04DvomXkJcDiwGvBh4GKKqwLfAz4CrAG8OzO/VyeAzLyNIvGfCRzRtvn4st6zM3MRQERMj4gtyzv2y4iIj1AckOuAvUzOJUmSJEmdNDEXrXMHncz8bER8h+Il9+0ontVfCNwAfDkz7xxLEBSJ/1XAKRGxF3AzsBOwB8XjBMe0lN2w3H4nxYEEICIOBk4AnqQYm/3IDl3qz8/MeWOMUZIkSZI0tTQqF62VoANk5l3AiXX361PnbRGxA8WX2gfYl+Lx+VOA4wd8qf555Xwa8J4uZa4A5i1ftJIkSZKkqaBpuWjtBH28ZObdwKEDlJvP0u7uW9fPBeYOOy5JkiRJ0tTVpFy0a4IeEa8o//jjzHy0ZbmvzPz+ckcmSZIkSdII6XUH/XKKAdlfSPHsfbU8iGnLFZUkSZIkSSOmV4J+AkVC/oe2ZUmSJEmSNGRdE/TyOfquy5IkSZIkaXgGHgc9IjaJiBl9yjwzIjZZ/rAkSZIkSRotAyfowB3AUX3KHFmWkyRJkiRJNdRJ0IMOXcpLkiRJkqTlVydBH8T6wKIh1ylJkiRJ0pTXqxd3IuKgtlXbdlgHxbBqmwAHAj8fUmySJEmSJI2Mngk6MI+lQ6slsF85tasefX8EOH4okUmSJEmSNEL6JeiHlvMAzgTOA77ZodyTwALgR5n54PDCkyRJkiRpNPRM0DPzrOrPEXEwcF5mnj3uUUmSJEmSNGL63UF/SmbuMZ6BSJIkSZI0ygbuxT0ito+I/xcR63fZvkG5fdvhhSdJkiRJ0mioM8za+4F3AL/rsv23wNuB9y1vUJIkSZIkjZo6CfouwGWZmZ02lusvBXYdRmCSJEmSJI2SOgn6BsA9fcrcCzx37OFIkiRJkjSa6iTojwDr9imzLvDY2MORJEmSJGk01UnQfwrsFxFrdtoYETOA/cpykiRJkiSphjoJ+uco7pBfHBFbt26IiG2Ai4DnlOUkSZIkSVINdcZB/1pEvAY4CLghIn4L/AbYEFgfCOCszPzquEQqSZIkSdIUVucOOpl5CHAY8EuKTuNeUs5/AbwzMw8ddoCSJEmSJI2Cge+gVzLzc8DnImJ1YG3gwcx8ZOiRSZIkSZI0Qmon6JUyKTcxlyRJkiRpCGo94i5JkiRJksZHrTvoEbEGcDjwaorO4VbtUCwzc9YQYpMkSZIkaWQMnKBHxNrAlcCLgIeAGcBCYBXgGWWxe4HHhxyjJEmSJElTXp1H3I+lSM7fDqxTrvsUsCbwMuB64DbghcMMUJIkSZKkUVAnQX898P3M/GJmZrUyC1cD+wJbAscMOUZJkiRJkqa8Ogn6xhR3yStLaHkHPTN/B3wHePNwQpMkSZIkaXTUSdAfAZ5sWV4IbNBW5rcUncdJkiRJkqQa6iTod1PcRa/8EnhFRExrWbcbcP8wApMkSZIkaZTUSdCvAHaPiCiXvwbMAr4dEUdExNeBnYELhhyjJEmSJElTXp1x0M+iGFJtI4q76WcAewL7A3uXZX5I0du7JEmSJEmqYeAEPTOvB97VsvwEMCciXgI8H5gPXJuZS4YdpCRJkiRJU93ACXpEvAJ4KDN/2ro+M68Drht2YJIkSZIkjZI676BfBrxzvAKRJEmSJGmU1UnQ/wAsHq9AJEmSJEkaZXUS9MuBl41THJIkSZIkjbQ6CfqxwBYR8dGImD5eAUmSJEmSNIrqDLP2IeAm4MPA2yPiZ8D9QLaVy8x8+5DikyRJkiRpJNRJ0A9p+fMG5dRJAibokiRJkiTVUCdBf964RSFJkiRJ0ojrmaBHxEHATzPzxsy8c4JikiRJkiRp5PTrJG4esH/riog4OCIuHbeIJEmSJEkaQXV6ca/MBHYfchySJEmSJI20sSTokiRJkiRpyEzQJUmSJElqABN0SZIkSZIaYJAEPcc9CkmSJEmSRtwg46DPjYi57Ssj4sku5TMz64yvLkmSJEnSyBskkY6addYtL0mSJEnSyOuZoGem76hLkiRJkjQBTMAlSZIkSWoAE3RJkiRJkhrABF2SJEmSpAYwQZckSZIkqQFM0CVJkiRJagATdEmSJEmSGsAEXZIkSZKkBjBBlyRJkiSpAUzQJUmSJElqABN0SZIkSZIawARdkiRJkqQGMEGXJEmSJKkBGpOgR8RGEXFmRNwbEY9FxPyI+HRErFOjjgMi4tSI+EFEPBQRGRFfHs+4JUmSJEkrribloivX3WE8RMQs4CpgPeCbwK+AHYGjgH0iYtfMXDBAVccC2wB/Au4BthyfiCVJkiRJK7qm5aJNuYN+OsUBOTIz98/MD2bmnsCngC2AEwes573A5sAM4F3jEqkkSZIkaapoVC466Ql6RGwG7A3MB05r23wcsAg4MCLW6FdXZl6WmbdmZg49UEmSJEnSlNHEXHTSE3Rgz3J+UWYuad2QmQ8DPwRWB3YezyAiou8kSVLT1WnPbPskSSu65WzLGpGLtmpCgr5FOb+ly/Zby/nmExCLJEmSJGk0NC4XbUIncWuV84Vdtlfr1x7PIHwqXpI0FQzSnlV3E2z7JEkrujrtXgeNyEVbNeEOej/V0fRXhCRJkiRpokx4LtqEBL26KrFWl+0z2spJkiRJkrS8GpeLNiFB/3U57/Zc/wvKebf3AiRJkiRJqqtxuWgTEvTLyvneEbFMPBHxTGBXYDFw9UQHJkmSJEmashqXi056gp6ZtwEXATOBI9o2Hw+sAZydmYsAImJ6RGwZEbMmNFBJkiRJ0pTRxFy0Cb24AxwOXAWcEhF7ATcDOwF7UDxOcExL2Q3L7XdSHMinRMT+wP7l4gblfJeImFf++Q+Z+f5xiF+SJEmStOJpVC7aiAQ9M2+LiB2AE4B9gH2B+4BTgOMz84EBq9oWOLht3WblBMWBNEGXJEmSJDUuF21Egg6QmXcDhw5Qbj5Lu7tv3zYXmDvMuCRJkiRJU1eTctFJfwddkiRJkiSZoEuSJEmS1Agm6JIkSZIkNYAJuiRJkiRJDWCCLkmSJElSA5igS5IkSZLUACbokiRJkiQ1gAm6JEmSJEkNYIIuSZIkSVIDmKBLkiRJktQAJuiSJEmSJDWACbokSZIkSQ1ggi5JkiRJUgOYoEuSJEmS1AAm6JIkSZIkNYAJuiRJkiRJDWCCLkmS/v/27j7Wkrq+4/j7I1AwaAGxQKxaBBUoTahoVJTIQyuFGqoFTPoHFAhQDdiV2k19SBUQiJj6RNUCirCATWpTQKy1ainPrg1tBSlSFZVFW0AFBBQEFb79Y+aut9d73QfOOfO7w/uVTGbvzOyc757sPd/5nJn5jSRJaoABXZIkSZKkBhjQJUmSJElqgAFdkiRJkqQGGNAlSZIkSWqAAV2SJEmSpAYY0CVJkiRJaoABXZIkSZKkBhjQJUmSJElqgAFdkiRJkqQGGNAlSZIkSWqAAV2SJEmSpAYY0CVJkiRJaoABXZIkSZKkBhjQJUmSJElqgAFdkiRJkqQGGNAlSZIkSWqAAV2SJEmSpAYY0CVJkiRJaoABXZIkSZKkBhjQJUmSJElqgAFdkiRJkqQGGNAlSZIkSWqAAV2SJEmSpAYY0CVJkiRJaoABXZIkSZKkBhjQJUmSJElqgAFdkiRJkqQGGNAlSZIkSWqAAV2SJEmSpAYY0CVJkiRJaoABXZIkSZKkBmw6dAGS2nfJwZcMXcIGO+QfDxm6BEmSJGmDeAZdkiRJkqQGGNAlSZIkSWqAAV2SJEmSpAYY0CVJkiRJaoABXZIkSZKkBhjQJUmSJElqgAFdkiRJkqQGGNAlSZIkSWqAAV2SJEmSpAYY0CVJkiRJaoABfRlJQpKhyxg93+fZOfTTh3Lopw8duozR832eDd/n6fAzeTZ8n2fDz4nZ8H2eDd/n6TCgS5IkSZLUAAO6JEmSJEkNMKBLkiRJktQAA7okSZIkSQ1oJqAneWaS85LckeSRJGuSfCDJNkPsR5IkSZI0fi1l0U03vPzJS7IzsBrYDrgM+CrwYuCNwIFJXl5V98xqP5IkSZKk8Wsti7ZyBv1v6P4hK6rqNVX1lqraH3g/sAtw+oz3I0mSJEkav6ay6OABPclOwAHAGuDDC1afBDwIHJFky1nsR5IkSZI0fk1m0aoadAKOBQo4Z4n1n+vX/8409tMvc3JycnJycnJycnJychrx1FoWXWwa/Aw63el+gK8vsf7Wfv78Ge1HkiRJkjR+zWXRFgaJ26qf37/E+rnlW09jP1WVdexXkiRJkjQ+g2bRxbRwBn1d5gJ0NbIfSZIkSdL4zTyLthDQ575N2GqJ9b+6YLtp70eSJEmSNH7NZdEWAvrX+vlS1+M/r58vdT3/pPcjSZIkSRq/5rJo+lHlBtM/0P0bdEPS71xVj81b91TgTrovEn6tqh6c9n4kSZIkSePXYhYd/Ax6VX0T+DywI3DCgtWnAFsCF879Q5JslmTX/k3Y6P0sF0kOS/LBJNcmeSBJJfn40HWNTZJtkxyb5NIk30jy4yT3J7kuyTFJBv9dGYsk707yr0m+07/P9ya5IclJSbYdur6xSnJE//lRSY4dup6xSLJm3vu6cLpr6PqWI/vebNj3Zse+Nxx733SMqfe1mEUHP4MOa79xWA1sB1wG/DfwEmA/ussAXlZV9/Tb7gjcBtxeVTtu7H6WiyQ3AnsAPwL+B9gV+NuqOnzQwkYmyeuBs+i+3boS+DawPXAI3b0kFwOvrRZ+YZa5JD8BvgTcAnyP7gPrpcCLgDuAl1bVd4arcHySPAv4L2AT4CnAcVV17rBVjUOSNXQjsn5gkdU/qqr3zLai5c++Nxv2vdmx7w3D3jc9Y+t9rWXRJgI6rP0leidwILAtXcP4JHBKVd07b7sdWeJN2ZD9LBdJ9qM7QPkGsA9dE/VAZcKS7E/XMP9pwSUpOwDXA88CDquqiwcqcTSSbFFVDy+y/HTgbcBZVXX87CsbpyQB/gV4DnAJsBIPUiamP0hhsX6kjWPfmw373uzY92bP3jddY+x9LWXRFp6DDkD/zeHR67HdGn4+TP1G72e5qKor5/7cfdZoGqrqiiWW35XkbOB0YF+6Mwp6HBY7SOn9Pd2ByvOWWK+NswLYn+7/7/7DliKtm31vNux7s2PfG4S9TxukpSzaTECXGvbTfv6zQasYv4P7+U2DVjEiSXYDzgDOrKpr+jNmmrzNkxwOPBt4kO7/8DVV9eiwZUkbzb43G/a9KbD3zYy9b0oM6NIvkWRT4I/7Hz87ZC1jk2Ql3T1hW9Hdh7c33Yf7GUPWNRb9/92L6O4rfdvA5YzdDnTv9Xy3JTm6qq4eoiBpY9n3pse+N332vpmy902JAV365c4Afgv4TFV9buhiRmYl3YBEcz4LHFVV3x+onrF5B/ACYO+q+vHQxYzY+cC1wFeAHwI7AW8A/gT45yR7VdWXB6xP2lD2vemx702fvW827H1T5CM0pCUkWQH8OfBV4IiByxmdqtqhqkL3DewhdB/uNyTZc9jKlr8kL6Y7c/Deqvri0PWMWVWdUlVXVNV3q+qhqrq5ql4PvA94MnDysBVK68++N132vemy982OvW+6DOjSIpKcAJxJ90iU/ZbjEwCWi/7D/VLgALrRLi8cuKRlbd7lfV8H3j5wOU9kZ/fzVwxahbSe7HuzY9+bPHtfM+x9E2BAlxZIciLwIeBmuoOUuwYu6Qmhqm6nOzDcPcnTh65nGXsK8HxgN+DhJDU3ASf123y0X7bY80s1Gd/r51sOWoW0Hux7w7DvTZS9rw32vgnwHnRpniRvprv/7kbglVV198AlPdE8o587AujGewT42BLr9qS7N+864GuAlwBOz179/FuDViGtg31vcPa9ybD3tcHeNwEGdKmX5O3AO4H/BA7w8r7JS7IrcN/CszNJngScCmwHrK6qHwxR3xj0g+Icu9i6JCfTHaRcUFXnzrKuMUqyO3Dnws+KJL9BdzYS4OMzL0xaT/a96bPvzYa9b3bsfdNnQG9cktcAr+l/3KGf75VkVf/nu6tq5cwLG5kkR9IdpDxKNyrliiQLN1tTVatmXNrYHAj8VZJrgG8C99CNaLsP3WA5dwHHDVeetEFeC7wlyZXAbXQj2e4MvArYAvgM8J7hylue7HuzYd+bGfuexsbeN2UG9Pb9NnDkgmU79RPA7XSP7dDj85x+vglw4hLbXA2smkk143U58BHg5cAewNbAg3SDulwE/LVncLSMXAnsQndmZi+6e+7uo7uM8iLgBlLTVwAABPxJREFUoqqq4cpbtux7s2Hfmw37nsbG3jdl8f2TJEmSJGl4juIuSZIkSVIDDOiSJEmSJDXAgC5JkiRJUgMM6JIkSZIkNcCALkmSJElSAwzokiRJkiQ1wIAuSZIkSVIDDOiSmpDk5CSVZN+ha5EkaRbsfZIWMqBLoj84mD89muTeJFclOSpJJvAaR/X7PmoCJUuS9LjY+yS1aNOhC5DUlFP6+WbAc4E/BPYBXgS8Ycqv/SHg74BvT/l1JEmaz94nqRkGdElrVdXJ839O8nLgGuD4JO+tqtum+Np3A3dPa/+SJC3G3iepJV7iLmlJVfUF4KtAgBfOX5fkhUnOTPLl/pLAh5PcmuS9SbZZsO1VwPn9j+cvuKRwx36bRe/D65ddleTpST6S5M4kjyT5SpKjF6s7yeb9/r7Vb3tbktP65dXXI0nSL7D3SRqSZ9AlrcvcPXg/XbD8OLrLAK8GLgc2AfYE3gQclOQlVfXDfttVwH3Aq4HLgBvn7ee+9ahha+ALwE+AfwC2AA4DzkvyWFVdsLbY7p7Bi4FXAbfSXT64GXAUsPt6vJYkSfY+SYMwoEtaUpJXALvQHRxcv2D1u4ATqurRBX/nGOBc4Hjg3QBVtaofa+fVwCeratUGlrIH8DHgdXOvl+T9wE3Am4EL5m17ON0ByrXA71bVT/rt3wH82wa+riTpCcbeJ2lIXuIuaa3+0riTk5ye5BN0ZwcCrKyqO+dvW1W3LzxA6Z0HPAD83gRLewh40/zXq6pb6M4s7JbkqfO2PbKf/+XcAUq//X3AqROsSZI0AvY+SS3xDLqk+U5a8HMBx1TV+Qs3TLIZ8Drgj4DfBLbi/3/p9+sTrOvWqnpgkeXf6edbA3OXFL4AeAxYvcj2102wJknSONj7JDXDgC5praoKQJItgb3oLq07O8ntVXXFgs0/QXcf3rfo7q27C3ikX3cisPkES1vqXr2f9fNN5i3bCri3qn62yPbfnWBNkqQRsPdJaokBXdIvqKoHgcuTHAx8CbggyS5V9RBAkhfRHaBcDvx+Va0dRCfJk4C/GKDsOQ8AT0uy6SIHKtsPUZAkqX32Pkkt8B50SUuqqpuAjwLPBP5s3qrn9vNPzT9A6b0YePIiu5u7h26TRdZN0g10n20vW2Td3lN+bUnSMmfvkzQkA7qkdTkNeBhYOe8Zr2v6+b7zN0yyHfDhJfZzTz9/9oTrW+jCfn5akl+ZW5hkK+DtU35tSdI42PskDcJL3CX9UlX1v0nOAd5Id/neW4F/pxtF9pAkq+kGoNkeOAj4GnDHIrv6It2ItCcmeRo/vyfug1V1/wRLvpBu8J4DgZuTfIruWbCHAv9B9+icxyb4epKkkbH3SRqKZ9AlrY930R1grEiyff/Ilz8AzgKeAaygu4TuXLpHzCy89I+q+gHdgcItwNF0j305Fdhm4baPR1UV3T2Cp9IdnPwp3TNoLwBO6DdbbFRcSZLms/dJmrl0v8+SNH5JXgl8Hjijqt46dD2SJE2bvU9aXjyDLml0kjxjkWXbAmf0P14624okSZoue580Dt6DLmmM3pdkD2A18H26kXgPAp4GnFNV1w9ZnCRJU2Dvk0bAgC5pjC6hG7jnYGBrupF4vwKcR3evoCRJY2Pvk0bAe9AlSZIkSWqA96BLkiRJktQAA7okSZIkSQ0woEuSJEmS1AADuiRJkiRJDTCgS5IkSZLUgP8D8D4qTftJc5IAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 1152x576 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "#Plotting Histograms for the two users\n",
    "f, ax = plt.subplots(1, 2, gridspec_kw = {'wspace':0.00, 'hspace':0.00},figsize=(16,8))\n",
    "\n",
    "ax[0].bar(bin_edge1[:-1], height1/np.sum(height1), width = 0.5, color='darkmagenta',alpha=0.7)\n",
    "ax[1].bar(bin_edge2[:-1], height2/np.sum(height2), width = 0.5, color='darkmagenta',alpha=0.7)\n",
    "ax[0].set_ylabel('Fraction of Total Ratings')\n",
    "ax[0].set_title('Most Active User')\n",
    "ax[1].set_title('Another Active User')\n",
    "plt.setp(ax, 'ylim', (0.,0.55))\n",
    "ax[0].tick_params('y', which='both', right=True, left=True, labelright=False, labelleft=True)\n",
    "ax[1].tick_params('y', which='both', left=True, labelleft=False, labelright=True, right=True)\n",
    "plt.setp(ax,'xlabel', 'Rating')\n",
    "plt.savefig('user_histograms.pdf', bbox_inches='tight', dpi=160)"
   ]
  }
 ],
 "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.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}