openai-cookbook/examples/Clustering_for_transaction_classification.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Clustering for Transaction Classification\n",
    "\n",
    "This notebook covers use cases where your data is unlabelled but has features that can be used to cluster them into meaningful categories. The challenge with clustering is making the features that make those clusters stand out human-readable, and that is where we'll look to use GPT-3 to generate meaningful cluster descriptions for us. We can then use these to apply labels to a previously unlabelled dataset.\n",
    "\n",
    "To feed the model we use embeddings created using the approach displayed in the notebook [Multiclass classification for transactions Notebook](Multiclass_classification_for_transactions.ipynb), applied to the full 359 transactions in the dataset to give us a bigger pool for learning"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Setup"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "%load_ext autoreload\n",
    "%autoreload "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "import openai\n",
    "import pandas as pd\n",
    "import numpy as np\n",
    "from sklearn.cluster import KMeans\n",
    "from sklearn.manifold import TSNE\n",
    "import matplotlib\n",
    "import matplotlib.pyplot as plt\n",
    "import os\n",
    "\n",
    "openai.api_key = os.getenv(\"OPENAI_API_KEY\")\n",
    "COMPLETIONS_MODEL = \"text-davinci-002\"\n",
    "\n",
    "# This path leads to a file with embeddings created in the notebook linked above\n",
    "embedding_path = 'data/transactions_with_embeddings_359.csv'"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Clustering\n",
    "\n",
    "We'll reuse the approach from the [Clustering Notebook](Clustering.ipynb), using K-Means to cluster our dataset using the feature embeddings we created previously. We'll then use the Completions endpoint to generate cluster descriptions for us and judge their effectiveness"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(359, 2048)"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "embedding_df = pd.read_csv(embedding_path)\n",
    "embedding_df[\"babbage_similarity\"] = embedding_df.babbage_similarity.apply(eval).apply(np.array)\n",
    "matrix = np.vstack(embedding_df.babbage_similarity.values)\n",
    "matrix.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "n_clusters = 5\n",
    "\n",
    "kmeans = KMeans(n_clusters=n_clusters, init=\"k-means++\", random_state=42)\n",
    "kmeans.fit(matrix)\n",
    "labels = kmeans.labels_\n",
    "embedding_df[\"Cluster\"] = labels"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Text(0.5, 1.0, 'Clusters identified visualized in language 2d using t-SNE')"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEICAYAAABS0fM3AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAABVRklEQVR4nO29eXycV3X//z6zL1pGm2VJI2/xFtshTmxlIUkT1iwNmNCUL2UNhORLf0ADbb4FSksohbbQlDRQ2hIaylIKBEMIpFlIAgkhcWI5dpzE8W7H1ki2rF0aSbPf3x9nRhrL2ixb+32/Xnpp5nmeeZ77zHI+955z7rlijMFisVgs8xfHdDfAYrFYLNOLFQKLxWKZ51ghsFgslnmOFQKLxWKZ51ghsFgslnmOFQKLxWKZ58xqIRCRL4jIf093O04HEXlYRD44wr4lImJExDVJ146KyLLsY7+I/EpEukTkpyLyXhH59QTPe5OI/P5M2zQZDH1PR3v/z+AaI34PReQKEdk7wfNO+H2dD0zG+yMii7LfSefZPO9MZ8YLgYi8R0S2ZT+cY9kf8uVn8fyTanyHYoy51hjzvcm+jog8KSIfGXLtAmPMoezTG4FKoMwY88fGmB8aY9462e0aypA2TcX1puT9z7ve08aYVVN1vdmMiFwiIo+JSLuItGQ7KFVT2QZjzNHsdzJ9ts8tIt8VkS+NcUxIRL4jIsdFpEdE9onIZ/L2GxF5WUQcedu+JCLfzT7O2bPokL//M9p1Z7QQiMifA/8C/D1qtBYB/wZsmsZmncRUCcgksBjYZ4xJTXdDLJYsJcA9wBL0+9kD/Nd0NmgauAsoAM4FioG3AweGHFMNvHuM84Sygpb7+8moRxtjZuRf9k2IAn88yjFfAP47+/gqIDJk/2vAm7OPLwK2Ad1AM/C17PajgMleKwpcmt3+YWA30AE8CizOO68BPgbsBw4Dkv0AT2TP/zKwboQ2Pwl8JPvYCdwJtAKHsuc0gCvvPbgXOAY0Al8CnNl9NwG/z76+I9uOa7P7vgykgVj2nv41r93Lgb8FEkAyu//m3Pny2rkaeAxoB/YC78rbVwb8MnuvW4G/y3/tkPt9GPj4kG07gXfmtyn7+DrgVdQANAK359/rkHPkv+4PgR3Z9jQAX8g7bsmQ9zT//d+Z97lHs8ddld13CfAs0Jk97qq8cy4Fnsq28zHgX8l+D4e5/6vI+16i38nbgZeALuAngG+E1w79TO7O3l838AJwxZDfwn3A97Pt2gVszNt/YfY96gF+mr3ul870/c3u/wBwBGgD/oaTf3cO4DPAwez++4DScdqAC4GeCX7vTnrfT8MeDPd9+Tvgmex792ugfDz3PuTat6K/twT6XfvVCO1+BXjHKO+JAT6N2p5cG78EfHe49o/3b9oN/ig3fA2QGu2GOD0h2AK8P/u4ALhkpDcOHXEcQFXZBfw18OyQD+MxoBTwA1ejP8wQKgrnAlUjtPlJBg3RR4E9QG32XL8d8iW8H/gWEAQWZL/8/zfvx5sEbkEF5U+BJkCGXmeEH/fAe5d3vt9nHwfRH/yHsvd/ASpWa7L7f4z+oIPAOtRoj/SD/ADwTN7zNahx9Q7TpmNkjRvaO7xwaNtGuJergPNQo/M69If9jlF+2B8Zpp23Zj+LIqAG/WFflz3nW7LPK/K+S18DvMAfoAbidIRgK9qrK0U7Gx8d4bUn3TfwPtQYuoC/AI6TFZHs5xnLttkJ/APwXHafBzVWtwFu4J2oQRqvEIz2/q5BDdvl2evciX4vc7+724DngHD2/foW8KNx2oBP5u5hAt+7k973idqD7PflILAS/a0/CfzjeO59mDZ9N/eej3LP/4mK+IeAFcPsN8AK1N7k7MgZC8FMdg2VAa3m7LkuksByESk3xkSNMc+NcuxHgX8wxuzOXv/vgfUisjjvmH8wxrQbY/qz5y5Ee9GSfd2xcbTpXcC/GGMajDHt6I8XABGpRH/UnzTG9BpjTqCjjvwh4RFjzLeN+jO/B1ShLrQz5XrgNWPMfxljUsaYHcDPgD/OBtH+CPh8tl2vZK89Evdz8nv3XuDnxpj4MMcmgTUiUmSM6TDGbB9PY40xTxpjXjbGZIwxLwE/Aq4c361CNub0JeDtxphu1OA+ZIx5KHvOx9De43UisgioA/7GGBM3xvwO+NV4r5Xl68aYpuxn/itg/XheZIz5b2NMW/Yz+WfUsObHH36fbXMa+AFwfnb7Jah4fN0YkzTG/BwVo3Exxvt7I9q7/b0xJgF8HjVEOT4KfM4YE8l+5l8AbhzLpSoir8ue6/9ln5/u924sTsce/JcxZl/2t34fg5/XWPc+ET4B/BD4OPCqiBwQkWuHHGPQ0cffiIhnhPO0ikhn3t+5o110JgtBG1B+Fn3wN6OqvkdE6kXk+lGOXQzcnXsTUfeIoD3FHA25B8aY36DugW8CJ0TkHhEpGkebqvPPg/ba8tvgBo7lteNb6Mggx/G8NvRlHxaM47pjsRi4OP+LhBrwhUAFalRGavdJGGN6gP9lUMD+BP2iD8cfoeJ3RESeEpFLx9NYEblYRH6bDTB2ocanfJyvrUV/3B80xuzLbl6Mil7+/V+OCm010GGM6c07zYj3PwLH8x73Mc7PTERuF5Hd2UyvTtR1mH+fQ8/ry/5+qoFGk+0yZsn//Ma67mjv70nf4ez3sC3v5YuB+/Pex92o23LEDouILEddircZY57Obj6t7904OB17MNLnNda9j0o2Uy8XzH04e45+Y8zfG2M2oJ3h+4Cfikhp/muNMQ8BEeD/jnD6cmNMKO9v92htmclCsAWIA+8Y5/G9QCD3JNuDqMg9N8bsN8b8CWpIvwJsFpEgwyt4A+qCyX8j/caYZ/OOOel1xpivZz+8NegX7P+No83HULdQjkVD2hDn5A+0yBizdhznPaV9p0kD8NSQ+y8wxvwp0IK67EZq93D8CPiTrGH3oS6wUxtsTL0xZhP6Gf0C/RHAqZ/twiEv/R/Ud1xrjCkG/gMV7lEREX/2Ov9ijHk4b1cD8IMh9x80xvwj+pmVZL87Oca6/zNGRK4A/hIdRZYYY0JojGHM+0TbXCMi+cfmf35n8v4eQ90+udf6UQOWowGNXeW/lz5jTOMI97kYeBz4O2PMD/J2ne73bqL24HQY696HMtRm/NAMBnOH9vrJjk7/HnWFLR3mfJ8D/oq8+5woM1YIjDFd6FDrmyLyDhEJiIhbRK4Vka8O85J9aA/oD0XEjfr1vbmdIvI+EakwxmRQHzVABv2CZYD8XPb/AD4rImuzry0WkT8eqa0iUpftNbnRL2Ase86xuA/4MxEJi0gJGlTL3f8xNDD1zyJSJCIOETlHRMbr8mgeck+nw4PAShF5f/Y9d2fv8dys2+HnwBeyn8ka4INjnO8htGf4ReAn2c/gJETEk+0hFRtjkmgQL3fcTmCtiKwXER/qXsinEGg3xsRE5CLgPeO8z+8Ae4wxQ79P/w28TUSuFhGniPhE5CoRCRtjjqBuor/Ntvly4G3jvN6ZUIgawhbAJSKfR+MZ42EL2gv/uIi4RGQTGizNcSbv72b0vXp91k3xBU4Wp/8AvpxzDYpIRfb6pyAiNcBv0OSG/8jfN4Hv3UTtwekw1r0PZczfpIj8Tfa35sl+Frdl23fKXBRjzJNocHms39+YzFghAMj6Qf8c/RBb0N7Fx9Fe3NBju4D/Dw22NKIGOZJ3yDXALhGJotkX784Ow/rQLJtnssPXS4wx96O9hB+LSDf6Zp+i2HkUAd9Gs3dyGQT/NI5b/DaakbQT2I5+0fP5ABqEejV77s2oe2I83I36YjtE5OvjfA0w4M55K+rOaUKHxl9h8If0cXR4fBwNgP3XGOeLo/f2ZrR3ORLvB17LvucfRd1RZF02X0R7ivvRbKl8/j/giyLSg3Ye7mN8vBu4QU7Ot77CGNOAJgz8FYPfu//H4O/lPcDFqMvwDjRTZ7J5FHgENXBH0M7GuNw7Wf/1O1F3SCcaA3kQHXGe0ftrjNmF+rV/jPaQo2j2XC4GdDc6mvh19vXPoe/dcHwENZRfyP9M8vaP+3s3UXsw0vlGuMZY9z6Ue9EYWKeI/GKk06L31Yr+9t4C/KExJjrC8X+NJh0MpXPI9/rPR7uXXIaJxWKZR4jI88B/GGNGFfEJnLcAFZsVxpjDZ/PcM53ZfO8zekRgsVjODiJypYgszLqGPoimgT5yls79tqy7JoimUL6MpmrOeebKvVshsFjmB6tQF2QnOgfhRjO+FOfxsAl1YzShOe7vNvPH1TAn7n3SXUPZaP02NH3tehFZivrUytBJEe/P+jAtFovFMg1MxYjgNjR3OMdXgLuMMcvRAOjNU9AGi8VisYzApI4IRCSMzv77Mpr98zY0C2OhMSaVzSv/gjHm6tHOU15ebpYsWTJp7bRYLJa5yAsvvNBqjKkY67jJrpz5L+gkmMLs8zKg0wyWjYhw8mzdYVmyZAnbtm2blAZaLBbLXEVExjX7etJcQ6JTtk8YY16Y4OtvFV2HYFtLS8tZbp3FYrFYckxmjOAy4O0i8hoaHH4jOnEjJIP1g8LoZI9TMMbcY4zZaIzZWFEx5sjGYrFYLBNk0oTAGPNZY0zYGLMEncH5G2PMe9E6MzdmD/sg8MBktcFisVgsYzMd8wg+Dfy5iBxAYwb3TkMbLBaLxZJlqtbpfRJd0AGj69NeNNrxFovFYpk6Zut6uxbLHCUC1KNZ1hXoGjjhUV9hsZwpVggslhlDBA2ZhdB1W46g67OcgxaYFLRSshUIy9nFCoFlnjITe971qAg0o0siHEDXJGlHqy8bdInkPlQwNjH9bbbMBawQWOYwOWO/F61m4kB71A50yYh16Ho5UabXsNaja+Q8CPSgy+n60VGBE10e4E3oSqEHgUvzXmeFwHLmWCGwzDFyxn8fajRr0MKQrcCLqBB0oT3tZ4Bz0ZVFFzJoWA2jLzQ11v7TGW3Uo2vCNKAi0J89fybb7hp06epDwBIGF9MqQEcOFsuZY4XAMoeIoKWtWoA9qEHdgxrh3egiVUl0ca9e1PAmUP97G+py+QJqbO9ieGNvgE+ho4orgK3Za4AKyhWobz+E9ujHGm1szrYpii50l1vcKoOuTNmKikFHtt3F2f1R8pbgtVjOCCsEljlEbiXHVPZ/AjXuzaiBTWe35VxEoP73RtRF1J7dfzdqaK9DDXGuV1+DisDdwFWoUX4a8KEunCbgH4CrGVxXPbes8EhunEZUcBLoKKUr+zh/1FGYbX8rGiPoRsVqvMtXWyyjY4XAMofYhvauX0SNaAA1oEfQuYv96BLQKVQM0qjbpQE10mVojcQoOs8xgi6d3Isuk/1idvsfopPl70cNfc5Qx9Cf1IvABXntGs2NU4MaeA8qACF0dOLOntuVbc/b0eV809n7uhIbH7CcLawQWOYQgo4E/Kjxbcs+7kVHBpnsMQ7Am33uQXvza1ED+wv0Z7EeXS/eg66p/oPs89vQuEJh9vw594wne51ydG31fEZz49wIvMqgWLhQEXBl2/cH6Prodaf3Vlgsp4EVAssc4kLUWJejhrkIdfcsQI1xCO21e9HedwlQhYpAD52d3ezf78fvd+DxXMWSJSk8nl8Bv8qe/41o7OAX2fOVoca/MHu+AOomCqDum4LscZ2M7MapAz6PZg1tQV1Di9HMoA9jBcAyFVghsMwhrkWDxd0Mun2WArWosV4M/Ax1EQWzfwXAJbS3G7ZtC+DxFFFUFCKVirNr11u44IJX8s7/MXREUYcGgM8Hfo2KgJNBQbgNdVE1oyOBsdw4dViDb5lOrBBY5hBh4DNowNaFjgwCqHH+EGps38twqZ2vvXYPHk8BgQD09a2gpOQ5Kiu3Djn/Q8AN2etsyp6nn8GsoWrU1WONumV2YYXAMse4Hk3b3Ixm5Aw1zmGG6523tFRQVhYlkykimSwjEDhEWdmLtLWto6zsb1ERuBcdQdyVd54bJvuGLJZJxwqBZQ5y+q6WTKaOTOYBXC5DdfVPWLDgcY4fv4odO77PtdfWoga/AE0dhZHnGVgssw8rBJZ5xnCzfqGurp69e7tZtuy/KC9/gUjkD3n88X/nzW+uzb5OUOMPVgwscw0rBJZ5QD3qKtqHZvFcgebvPwV8FwhSXn4JsIZ4vJ7GxkJefllFIHySFylfDEJYEbDMFawQWOY49cCdaKqnA52x+yCaSbQEzdVvB16lvLwQDSp3UVOzLXvMUHJiYEXAMnewQmCZ42xGRaAMze6JAcfQyWBhNM3UQVNTkN279/Pii+WUlRVy5ZXNLF060jmtCFjmFtOxZrHFMoU0om6cLrRwWxydbBZHy1MnaW938eyzPpzOLioqwJgomzdXUF8/bY22WKYUKwSWOU4NOrP3GINlHtLoDGAB4hw+XERFRQciRfj93ZSXd9LeXsfmzdPTYotlqrFCYJnj3Ii6gU6gJSdyhdwKs/8X88orF+HzpYjFQiSTAXbv3oRImMbG6Wu1xTKV2BiBZY5TB9wO/D06KlgIvA2tM7QTjRnU8otfvBOXazBFqLMTamqmvLGznkgE6uuhpQUqKqCujiGZV5aZyKQKgYjUAt9Hp3oa4B5jzN0iUgr8BE3beA14lzGmYzLbYpnP1AHfYHBh+FwxuBXAJtasCfPQQ1BWBqGQikBbG3zoQ9PV3tlJJAIPPKDvYWUlRKP6fNMmFQMrEjMXMcZM3slFqoAqY8x2ESkEXgDeAdwEtBtj/lFEPgOUGGM+PdJ5Nm7caLZt2zZp7bScXWbuD37kJSTr62HzZmhs1JHAjTdquy3j5ytfgV/+Eo4e1ecLFoDPBy4XrFkDRUXwutdBQQEcOQKvvALLlsGqVYPv9cz83sxeROQFY8zGMY+bTCE45WIiDwD/mv27yhhzLCsWTxpjVo30OisEs4f8XmFBgfYKX3tNjYIxp/7AZ65oWE6H+np43/t0NOX1Qk+P/vl8sGKFjhA6O+Fd74KSEvjNb1QwurqgtFRFYsUKWL9+8HvT2Tk4mrBMjPEKwZTFCERkCbps0/NApTHmWHbXcdR1NPT4W4FbARYtWjR0t2VGcGoPu74+TCikP2yAPXvg5z+HdBqqqiCRgGQSNmyAP/xD7RkO70o4nQXgLdPN5s0Qj2vvP/cZu936uTc1qfH3+eC++8Dvh2PZX7/XC4WFsHcvNDfryMHhGPz+1NdbIZgKpkQIRKQALQT/SWNMt8jghBxjjBGRU4Ylxph7gHtARwRT0U7L6RBh0Oc+uEh7b+8mCgv1l7tvH/zoR2oUYjH98Xs86go4fBg+8xkVAZdL/fOvfz0sXAgvvxwhHD713CMvAG+ZDvJHc1u2qAFPJvV/JgNOpz4X0Z5/IgGplApGTiwKClQg3G7tCOzfD+Xlev6CAhUHy+Qz6UIgIm5UBH5ojPl5dnOziFTluYZOTHY7LGebetRQ5xZn1//nnFPPoUNhiorg8cf1xx0M6o/e5xvsIQYCcPCgGo3KShWJnTtVDDZtGv7cIy8Ab5lqhgaGi4rUwDscKgA5MfB6obhY96VS2vtPpfTP4VB3Iehx8bgKRo5oVN2FlslnUucRiHb97wV2G2O+lrfrl8AHs48/iHb3LLOKFjT7Jp8CVq9uobMTurvV7eN26+Ncb9Dh0Iyc115TgwG6v6VFn7/0EvT1tdDaeuq59ZqWmUB9PQMuQIcD3vAGfd7fryO84mIVfbeb7Gxt7eGvXav7QJ+DCoDXqx0Gt1sFpLtbYwQ2YD81TPaEssuA9wNvFJEXs3/XAf8IvEVE9gNvzj63zCoqUJdNPlFKSirYtEl7/Mmk/nm9g/7ivj51Ezmdui13TCajAtHbCyUlFRw+fOq5R14A3jLVtLQMGnKAlSvhpptg6VI17Mao0V+9WjsBtbXw9rfDH/0RfPCDGiMyRvdlMuoS/IM/0Nc0N+v3xwaKp45JdQ0ZY37PyBW63jSZ17ZMNrl1e2HoIu3hsP6Af/5z9R37/Wr4T2QdgE6n/vX3D+5Lp1UEzjsP3O464vHhz22ZGVRUqOsmF9QFqK6GO+7QXvzQTDBQV1J3t2YJveENg599MAgbN8I111jDP13YmcWWCZK/bu/wi7SvXasphI2N2vsLBrUH2NU1KAZOp/YIUykNKhoDe/eGWbFi9HNPJnZOwdjU1alhh5PTPa+8koGOwFA2bdL3trlZRwjvfKc1/DMFKwSWM2D49X9zrFqlQ/zjx9X4FxerCyAW0zRCl0sDxydOqCCce64KwTPPwKWXjn7uyaK+Hu68U7OYFi1S43bnnXD77VYM8gmHTzbsFRWDIjDaa6zhn5lYIbBMGnV1aujXrj251/ie9+j+3/1Ot5WUaMpgLn/8sss0K2U6DO/mzSoCZWX6PPd/82YrBEOxhn3uYIVgGpnuWbWRCDzyCGzbpm6ZCy+Ea689e20Yrdf4nveoq+iFF3SUEI9rjODii9WH3NwMDz4I996rYlJdDTffDNdff3baNhKNjToSyCcUGiybYLHMRawQTAORCHzzm+pjTae13sqll6rBm6pMiUgEvvc92LFDjXB/v84C3rcPbrttctqwfz/8+tcaN6ipgSuuUAPb3Ky56OvW6cigu1snnH3rWzpaqK3VkcQdd6jr6G1vG/kaxqioTZSaGr1WbiQAthKpZe5jhWCKiUTg7rtVBHy+wQJcbW3a252qKfX19XDoEHR06CSf0lI10Nu2wcMPwy23nPmIJX/SUXe3ZhF5PFp4rKcHfvADeP/7B8tMFBQM5o/v2KEikO+iOXQIPvlJfZ+GM/bGwKc+pef6whcm9r7ceKPGBMBWIp1LRLoi1DfV09LbQkWwgrpq9fPlbwsXhon0RE46Jlw8P3xfVgimmPp62L1bc+iLitSgORzaK3/xRXWTTAUtLeoGKSjQPH9QQUgkYPv2sUsKj4f8SUc//7mKjcejI4DVq/WYp5+GT3wC/ud/4IkndGSydi00NGjwOIcxGlA+cECN/V13nSwGORG4+24d0Ux0ZFBXp4HhzZt1tFJToyJg4wOzi3zDLyLsb91PwiRIpBJ4XB5+feDXdMW78Ll8VBZUsrd1L881Pke4KMzSkqUc6znGw/sf5pySc1hZvnLOi4IVgimmpUV7w0VFmjLpdmv2TDI56EefCioq1OjnRAD0eTCoRjTfiMP4ioANHUHs2zdo8NvadJuIxgJg0Pd+7JjGCs45Z7AX3tmpI4Vly/RYEfXdZzJq7GFQDIaKwFCROF3q6qzhn61EuiI8cuARnjj8BJUFldQU1PDIwUc41HGI5SXLCReG2d+2n53NOykPlHNR9UVsa9zGyy0vU+wtpthTTHtfO6+2vMqykmV0xjvpS/bxwN4H2LRq05wVAysEU0xFhfa8o1EVBBicYVlaOnUGqK5ORx9Hjw66X3IZPBs3qjGvHFITdrQiYMONIA4e1ElDS5fq/4MH1ZD7/ZpOmkppj3u4TJ3XvU5HJsXFJ4vDv/yLjhzyxeBsioBl9hLpivC9F7/HlsgW+pP9tPa18uzRZ+mOdxNwBdjXto9XWl7BgYNEKkFjTyNPvPYEhe5C3OImnorzcsvLHOk+gsfpoaGrgRN9J+iOd+N2unn4wMPcsuGW6b7NScEKwRRTVwfPPqsukYICLbnQ0aHGbrKCtMMRDg+6WJqbdWSyZIn2wK+5Jpvpc6SX1PF2Yl0xfMU+XAtLEX+Q++8/NW4w3Ahi3TpdfATU+Hd3q2uookLdYMXF8PnPqxEfmqlzwQU6cigsVDdRdfVgfCCXOXT33YOCYEXA8siBR9jXvo94Oo7P5WNf+z56Ej3EU3H8Lj8pk8ItbtpibRhjEISkO0lDTwMucZHMJCn1l9IT78EhDg50HGBp8VJWlqxEEH57+Ldcu/zaOTkqsEIwxYTDarTKy+G3v9XRQF3d9Pih6+rgq18dPiDcfaybLT9uJ1QmFIf8dHamOfpyO+VrHAQC/lPiBi0t8Nxzmo7a06MG/JprtAbNnj3amy8v11hINKq9/iuu0OuNlKlz6aW66hUMup3uuUfb+Rd/MSgCYEXAAtuObaPMX0ZrbyuRnggpk6LQXUh3vJtYKqaL5QLpTBqHw0Eqk6K9rx2XuEhLmrRJ0xHrIGMyeBwe3OLG4XCwv2M/tUW1LChYQH1TvRUCy9khHIZPf1r/ppuRJgVJpJE3X2bYd7yYti4npSGhMp4mGG+nqEhzKfPjBs8+Cz/9qWZCFRXp7OGf/EQXn3E6deUpv1+39/aqAKTT+vqxMnWGup16euADHzi5vcMFkC2zlwf3PMi9O+6lqaeJ6sJqbr7gZq5fPfokEjGCwYCAMQa3w00sGcNkDG6nm2Q6SSqTwoiOBpwOp44KTBJjDOX+cvqT/aRMilg6RrgwjFO0RO6RziO8a927aOmdmxVwrRDMIc7mBLXell4WLS5g8dLBKqA/frgIR2LweWurBoQbGnQk4HRq8FlE/6dSuiTh+943mB0VCOhrt2zREUCuh//+96u7bLhMnXy3kzFanuLJJ1VkfvWrwRgBWDGYCzy450HueOoOSnwl1BbX0hnr5I6n7gAYVQwurL6QZxqeIZFOsKhoEZGeCM3RZvxuP1UFVTRGG8kkM7jFjdftJZaKqSAYJ16nl6A7SDwdp9xXTpG3iL5kHwZD0B0k6Anic/kIuANT9TZMKVYI5giRCHz96/Dqq4OumS1b4M/+bGJiEKwIkogmSCfStO/XOEH6aBX9C0OAisDzz6txD4e1p+92q6FOp1UUyso0/nDeeWrMYTBQvGOHliOurBxcyPycc+CSS04VsFzg2hj4z/9U45+LFYio8QcrBnOFe3fcS4mvhLKA+gpz/+/dce+oQnDt8mtpjjbT2tdKf6qfRaFFYKA0UEpjTyPl/nJSnhS9yV6SJkllsBJB6Ix1UuwrptBXiM/lI+gJEvKFONRxiI1VGwl6ghgMnbFOrlw8NyvgWiGYI/zP/8BTT2nm0YIF6od/6in1y//lX57++QrDhTx/1/N0vNaBO+BGRAg0tfPy7mVEm6O0u6tIFxTiCfhYtUqFJ7fOQM5l1NWl230+LR2xf78Gxo8fV1fR0qUqKK++qim0nZ0aPB86X6GiQsXtvvtUBN72Nl0EPRjU/VYM5hZNPU3UFteetC3kC9HQ1TDq68LFYW5afxOVBZX89vBvWVCwgEVFi8iYDDWFNTgdTrY2bcXn1p590BOkta+VFWUrqAhWUFVQxd7WvcTSMcQhXL7ocgKeAM3RZt6w9A1zNlAMVghmNF2RLprqm+ht6SVYEaS6rpricPGwxz7xhKZ+5mftGKPbRxKCkc7fFemiaWsTLr8Ld8BNd6SbjMlQVhpgY0kTrzWUc9S1kIWlbZx3dQnl5QGuuUZjBD09aqCjUTXqn/iEGvhQSMUgGlUxeP3rtQ379+vxPp8eN9x8hVzJY5dLRwHvepeKzFVXDd5LvhiEQlYEZjPVhdV0xjoHRgIAnbFOqgurx3xtuDjMLRtu4drl11LfVM++1n0c7DhIXU0di0OLWRpaymOHHqOmqIZlJcvwu/zsad1Db7KXgCfAB9Z/gN5EL6+0vMKy0DJWla+a85PJwArBuIjUR9i9eTfdjd0U1RRx7o3nEq6b3C9GpD7C9m9vx6QNwYogqViKnqYeVm1aNawY9PerAcwnldL6RTk/fL7LpSvSxd4H9uIL+SioLCARTbD3gb2s2rSKpvomfCEfLp8LX8iHp8BDtCmKSRtqag2l7UcpLS7FXRlC2tqAAO98p44IfvN4iuONhiJ/ilve1sFf3Bqkh+KTCs+94Q1q+EENekmJtj9/CcP8+Qq54nXV1VqyOhhUERjq8sqJgRWB2c3NF9w8EBMI+UJ0xjrpiHXwyYs/Oe5zhIvDA8Y7N8u4OdrM6xa+jutWXHdSKYkPX/BhYLDcRG1xLe88951z3vjnY4VgDCL1EbbcuQV/mZ/iRcXEOmNsuXMLl95+6aSJQVekix3f3oHD5cC/wE+qP0XLqy1UrKmgqb5pWCFYu1b97Ol4gnRPPx2dhqY2H8uWO3B2d/PYT1r45lEvi6tSXHp9KeeUdFAW8tGRCPLy8z7au5wUuvvpffgEhaaXgsoCvMVe2va14SvxkUlnAMgkMvhDfpYEWjiYrqT1eC81Ge3pb1zXz3vW7mHxUieeAk9WXCKs2rSKG24YbHMuCwjUddTerhPN1q3TbcMtWj7eksdWBGY/uTjAvTvupaGrgerCaj558SfHzBoaiXxRyFHHqbna88nwD8UKwRjs3rwbf5mfQJlmC+T+7968+4yEIFIf4cX/epGWV1pwBpwse8MyznvveRSHi2mqbyKdSlNYUYiI4A64AYg2RXH5hv/IPvxh+PxnE7Qd6cXpcdDV7yXgS7HKE+H5H7RzNFWNv8BBd0+KnT8/wM7iYq54cyG7DgYoCKYpK0nR1+vlkd86uOYNxXii/ZSuKKXphSYSPQnEKWSSGZKxJMWLiikrF6oWtrO3sXCgp3/hgghlASfeIq1bkfs/VLzyy1OHQioE69ZpfCNXdO7KuRmTOysMVz68qkqzrubKqmrXr75+wobfcvpYIRiD7sZuihed3AP3hXx0He2a8Dkj9RGe+uJTxDvi+Ep9ZBIZ9ty/h97WXi657RJ6W7RHnoqlBkTA5XcRPRal5uLh6yHX1cFNb4nwyNMB2qJ+xJXm9ef3kn7lODvbK+n1B0jHBPBTvaiXZGs7Dz91DmtWpSkIaG/f64hTtsBDk5QT7NyFL+Rj2ZuXceChA2RSGbyFXtxBN52vdRLviuNr7eOGWy4gnDU42+7pwpO/ojmoW6l56EL0J/fwc2mv413parKZjnUQxkt++fBEQudkbN2qYrpokabn5gLwn//87BYDy9RhhWAMimqKiHXGBkYCALHOGEU1RaO8anR2b96NSRgC5QGcXif4AYG23W001TcNxARaXtXJKy6fi/6OfsQpVNeNHDA7p6SdT9+SQBxdPLalgP64sLvLweGuMlxpD6m0A5cjwysnKjnHfYRjzQ7OX9kPGTfJWJJkb5Jw3QKi6eBgrKDYx3nvOQ8E2g+207StiUBpgJKlJRRUF9C0tYnCqkKKw8UDKae5kQBAIpogWBEc9f2YSStdPfigrnswdB0EmBliMLR8eFmZzuXo7tYYS2WlCkR7O3znOzpSmM7Fjyyzg2kTAhG5BrgbcAL/aYz5x+lqy2jUXlHLM//wjPrry/24A24yiQzrP7R+zKyekfZ3N3ZjHAaHxzFwrLvATd+JPnpbellx3Qp6mnqoWFNB9FiUnmM9OF1OVr59JQceOcCxbccwYqi+sJrl1y4fuGa+IT5vRYzfPF/Awb4qYmkn3rQDJ4YCX5LWTheukhrOfX0hiUySTGcfvmIfC9YtIO0JUpHVvL7WPo7tOIYYoWpjFaXnlFK5rvIkQx/vjg+4fqrrqtn7wF6AgRhBrDPG4isXT8EndXa4995T10HIbZ8JQjBc+fD+fp3D0dMzOJkvFNKifWdaSnymMNx6AvPZp3+2mRYhEBEn8E3gLUAEqBeRXxpjXp2O9oxEV6SL7iPdnHP1ORx/8TjR41E8AQ8X3XYRhVWFI2bd5FIwR9pfVFNEf0s/mURGRwRAMprEW+hFRGiqbyIejRNrj+Er8VFzUQ2F4UIO/vog7fva8Zf5wUDDMw1Em6Osv2n9KYa4stRw2Zo2HnykgqAnSTLtoKgwRcCRoCctdLgqeO+H/GzdWjOwKExuTeELVnbz4ndfpONAB74SHwgcffoo/Z39rH776pOEIN/1UxwuHhhJRJujBCuCLL5y8Ygpr0OZjuysoTQ16Uggn1BIZ0/PBCoqtGR5fvlwh0OD7fkkk1rX6XRLiQ+lPlLP5t2baexupKaohhvPvZGqwqopNcqRrggP7H2AkC9EZUEl0USUX+z5Be9Y/Y5TMoNybdpYtZHaUO0YZz47bctf90AQWvta6Yh1UOorpSxQhsFgjJnRAjZdI4KLgAPGmEMAIvJjYBMwY4SgK9LF1m9spbell8KFhSy6fBH9rf30HO+h4ekGeo/14gv58BZ5advfxqHHDtG6r5Wn/u4pPAUeAEJLQqy4bgVlRWV4i7z0tfWx9RtbQaCvrQ8jhsKqQjKJDLHOGAvOX0D0RBR3wE3F6oqBHnV1XTVN9U30t/QTKA8MxA3EIfS39g/0yIca4traIOuvKMSRTtC6t4PjrU7i4qVieZCFS/3U1Q26DvL98z31jfS39uMv0xFQrDtGrCNG5+FOdt+/mzV/tAaD4dgLx2jb3Ya4BBEZGJ2M1/DnMx3ZWcNRXT18AbzqsVPYp4S6Ou3d55cPLyrS2kx+vwpCb6+6jpYtU4HPZ7RS4kOpj9Rz55Y7KfOXsah4EZ2xTr741BdZU7GG1y183YBRHlqrfzjxqAtPPFjxyIFH2H5sO409jQhCW38bAVeA6sJqwsXhU4SiJ97DB37xAS5YeAFfu+ZrE77uWORf1ylOnm54mmg8itflpdhbTGNXI7F0jAJvAVfUXjGj1zWYLiGoAfL7WBHg4vwDRORW4FaARUNrFE8yud58X2sfBQsL6G3t5cjTRyg/t1yfN/fS1dDF0jctJdocZef3dtJ7opd4NE46liYRTeDyuug0nbzy41dY9+51+Ep8tOxqIZ1Ms3rTalx+F4cfP0y0OYqvxMfqG1ar4fW7h826iZ6Ikkqk8Jf4B9rp8rvob++nt6VXy+qKnGKIL22Dp5/2seyqItb4NNe/tVUrf8Lw/vltD/XqtUr9xLpjtO1tw+V14SvzEWuLse9/92HShr7WPhwuByW1JaeMTk6XycrOOl1uvnkwJpArgNfRoSWwZwL55cOPH9ey3hs2qLvI69V5FoWFKujLl+sorygvnDVcau5IbN69mTJ/2UmlHg51HGJXyy4uX3w5AEVePXmuKudw4nHnlju5/dLbJyQG9ZF6vvXCt+iOdVPkK8Lr8HKo/RBtsTZeOvES39/5fQRh/cL1LCpehDGG+169jyePPKmlIbK/i8mgvqmekC9EkbeIXSd2UeYvo72/nUQ8wbKSZUR6IgCU+cs42HGQS2svHXjdSCOZ6RoxzNhgsTHmHuAegI0bN5qpuGbOp3/w0YO4fC4y6QwnXjlBd2M3khGijVG8hV6ClUFSfSlaXm6hdW8r/W39YMBkDA6XA3EI6UQaYwwmaWh4poGyVWWIU/B4PDQ+30isK0b1hmoq1lWw4SMbANh2z7aB0UQOT4GH7f+5HQyUnFNyUiZRqj+Fy+MiUB7g0U89ii/k44KPXHBSXOL159fQ3FxES4saNI9HS0Nfc83I70OwIojL4yLVnyLaFMXt0+v5inwULSqi+7VuOo90UrK0hNDSEN4iL8m+5Emjk9NlMrKzJkIuDnDvvaeugzBTGK58eDisGUX5QWEYnK+R7/obb2puY3cji4pP7oSJCD3xnpO2FXgKaI7qMGM48chtH0sIhhrFcGGYb2//NslMkqAnSCqT4nj0OIWeQnoSPXTEOthxfAeVgUpa+lpIppI8evhRth3bxoaqDVwWvmzcIpC79t7WvQNunbGWqGzpbaGyQFdv6op3UeIrIZVODezPPfa7/HTEOk55r4ZzeU3XiGG6hKARyHfghbPbpo18nz4OiEfjnHjlBC6fC5MyiEdoO9CGt9TL0jcuhQwcevwQHYc6wAkmYTBpDQBn0hlM0hDviuPwOuht7SXWFaN1dyvJWBJ/iZ+FFyykuLaYw08cZvk1y0fMuon3xMHA/v/dT82lNZQuLyVQHgADsY4YvnIfz33tORqeaaB0dSmte1sJXxomtDikrqWte7jhravZGykad+ZIdV01J3adoONAB/0d/bgDbhLRBIU1hVRvqIYMxHpiLDhvAeLQH1r+6GQiTEZ21kTJX/xmpjLcSG64VNHcfI2JpObWFNWcUurBGEOht/Ck46KJKBVBHWYMJx4hX4ijXUdHvdZwRvHbO75NV7yLUm8pnfFOehO9eJwe2vrb8Lv8FDoLOdp9lEQ6wbLiZXxrx7do7W+lrrqODQs3cLjrMJGuyJhGNXftdCbN4c7DOMVJZ38nfrefpp6mEQ1zRbCCaCJKkbeIYm8x/al+XE7XwLoHLqea1/5UP8XeYlp7W/l9w+850HaAn+3+GbFkjIUFC/G5fJzoO4HBUFNYg8fpmfKV0KZLCOqBFSKyFBWAdwPvmaa2AAyUVfAWefGH/Jx45YSmPTrApA3x7jjeIi++Yh/B8iDx7jhL37SU4zuPk4wmyaQ0WpdO6KIXOMEddJPqSdHr6CXWGkNE8Jf6ySQzNDzdQPriNKVLS0fNuol3xbnu36/juX9+jufvfh5PwIPb7wYHlK0u49XNr3L8heMsvGAhBdUFdEe6adzaiLfQq4KBri1www3jN6jF4WLW37SeAw8foPtn3aR6UyxYs4CFGxYSLNfRgrfQO+zoZKxU0ZE498Zz2XLnFkBHArHOGP1t/az/0PoJnc+inElq7o3n3sidW3ShiFypB7fDzZqKNXTHuynwFBBNRE+qyjmceHTGOqkpGn7+S458NwuoyymVTnGi9wQNXQ1Ek1EdiQi64pjbTzgYptRXyv6O/RzvPQ5AXXUdl9deTl+yj3UV68a1kEzu2rtO7KLAU0DAHaAv2cfx3uOsrVg74jnqqut4YK8Ouc4pOYenG57G4/DgdXlp62ujyF1ELB2jrb+NteVr2fzqZl5qfolkJglAT7yHjMngFCe1RbWEi8Mc6jhEa2/rlBe4c4x9yNnHGJMCPg48CuwG7jPG7JqOtuTobekdcMuUrigl1hXD4XaAgYrzKvCV+Kg8vxKTUVGIdcaoPL+S8pXluINuXAGXikbSkE6kcXqc6rapCCBGKFlegsPjGJgp7PA4aN7ZTMV5FQO96Fyw1x1wE23WoPGqTasI1Ya4+q6rufi2izn8xGHcfjfX/9v1NPy+YUAEll+3HAz4S/0kuhO0728HVFAm0ksvDhez4ZYNbPrPTax+52rCl4YJlAaId8fxl/upWFtBX2sfid4EiWiC/rZ+/OX+Uec5jEa4Lsylt1+Kt9BL19EuvIXeKQ8UW06mLlzH7ZfeTqG3kKNdRyn0FvL5Kz/PJy7+BAG3VuUMuAMn9ZhvPPdG2vrbaOtrI51J09bXRlt/Gzeee+Oo12rpbaHAc3JkO51Jc6DtwIB7J2mSJFIJHOLAgYP2/nYqAycvrL2xaiM+t4+LwxezOLR4XAvJ5K7dFe/C79IYnM/loyvWRYGnYMRzhIvDbFq1iYA7QNqkuaz2Mq5Zfg3nLTiPQm8h51WexzXLr+Gy2st46fhLHOg4QDKTxOvy6r2kE8QzcXBAa38rx6LHMMaQJk19U/2Y7T6bTFuMwBjzEPDQdF1/KPlumUB5gMrzKtWYCgTLg5TdUEa0ITrQC1585WKa6ptY8sYl4ICOgx2k+lKk0diAv9RP6cpSlr5xKVu/uZWKNRWk42kS3Vrj3+V3kU6mcfvcA71qYMSsGxHh6ruuBuD5u5/n+bufB2DBeQs455pzBgQmlUiRSWWIdcWA8U3oGo3hUkLX37Qe4KQ5DbWX1Z40p2EihOvC1vDPMOrCdSf59usj9Xxj6zcGMoKqC6qpb6rnof0PDQQ7b7/0djbv3szRrqPUFNXwofUfGjM+kO9myRHpjlDoLSRcFGbnsZ24xU3SkcTn1DUDMLCzZedJ54n0RLhm+TWICN3x7gGX1XiunXPvBNwBYqkYxb7ik9xewzFcHaPh2N60HY/Tg8/tI5VOYTIGp8NJOpMmk8kgLqE70U1XTMXori130do3dSODGRssnmqGumVKV5bS29JL7WW1A/52j89zUvXP/Q/tJ7Q4xKq3raJ9fzvHth/D4XbgK/Wx9sa1gE64Ci4IEuuMUVhVSHemG5fHRbIvibfAe1oTrnJikBMBgCVvWEI6nsYRcFBQVUDLrhacLq33kxu5nOmErpHEacNHNsBHzujUlhnM0OCtW9z84OUfDGQENXQ1cMdTd7Bp1SbOX3j+ScHOr7zlK6d1rXw3S87l1JPooa6qjr50H06nk7A/jNfppTPeSU1hDXva9tCT6GFD1QbuuOIOvvb81/jVvl9hjOH/rP0/dMW7xrWQTO7aCwsWsqtlF/3JftKZNIuLF5+1xWh6Ej10xbroT/WTSCUwGFwOlwqByWCMoSfeg8vhoixYRrGvmGcanqE52sxN62+adDGwQpBlaM+3uLaY8O1heiI9I06Oyo0iAuUBAuUBSleU8tqTr+Ev8WMyZmAewIY/3cC2b2wDoGBhAT0NPcQ6Y5z3/vNGLCs9HMYYHv3Uoydta6xvHAggO5wO/KV+HG4H3pB3YORyJr10y/xkuODtd3d8l6rCqgH/fywdo8RXws7mnVxQdcEpqaSnQ87NkisXXRGsYEPVBkSERYFFRLojJNIJAGo8NWTI0BnvZEPVBupvqUdEWF+1nj996E95cP+DGAz/ft2/j6sd+dfuS/YNZA3VFteelXTOSFeEDBk8Tg/JdJIYMZLpJGIEl8OF2+XG7/LTn+rH5/QR8oZYGlqK2+mmtb91Qu/n6WKFII9he76jjGiHjiKcHielK0spWFBwingEy4PsuHcHPU09hJaHuODmC1h9/epxty0nAs/f/TwX33YxV9919cDz/ADyOW89Z1wumpkwi9cycxkueNuX6qM3ORhv6kv2EfKFaO1rHdiWnx55ugx1s4QLwwPB6pWlK3m64WnECF6Xl91tu1lRuoJ/u/bfBmIItaFafvUnv+JTj36Ku5+/m39+7p+56+q7xpVCOl4Xz0Sob6rnkppLiCfj7GrZhcHQnekmQ4YKfwWVhZU0dDdgMCwvWc55ledR7CvGGEN7f/u44hxnihWCM2BY//kHh59Qtfr61adl+PMZTgSGxgwWrFswsH009jy4h2e++gytL7fiKfZQfVE18Z44T33xKRZfsRhjzEBpi/JV5aOuimaZu+TnyOdYEFygRinPZb63fS9ucfPrg7/GYOhL9FERrKA+Un/S4i8T6VnngtWbd2+mO97NZeHLQGDn8Z1sqNrAv137b1xUe9FJrxER7rpal6oL+UKTNpnsdGjpbWFxaDHXr7qe6qJq9rfvJ5aI0Z/uZ2XpSgo8BdRV15HIJAj5QgTcmu3Xn+rH4/KMK85xpogxUzJX64zYuHGj2bZt23Q3Y1oYSQTGuz+fPQ/u4ak7nqKvtQ+cIEZI9adYWLeQRFcCh9tB6fJSHE4HJm2oWFuBw+k4LfeVZfYxXEmISE+EvmTfScHb7U3beezQYywrWUZrXyu7TuyiL9XHhZUX0p3opjPeidfpJegO0p3o5uKai1ldvppIT4QT0RNndd3frQ1b+dmen41YxmIyZxSfLvfvvv+U97I73k3AHeCGc28A1H30vRe/x772fZT5tT5RR6yD5SXLzyhGICIvGGM2jnXctKSPWsaPiOAL+UY08rmRwcW3XYwv5Bv1y7/j3h1aRA5weV24/W5cfhfHXziOL+Sjt7kXb4GXQFlAi8kdj+IL+Wiqb5rUe7RMH7mSED3xHhYVL6In3sOdW+7ELW46Y510x3W96u54N+XBct656p3sbN7JC8deQESoClax88RODncdJpPJkEgl6Ip30dHfwW9f+y3f2fEdjnQcoSfRw49e/hFf/N0XqY+cWWpkfaSef37un09pc/55Z4oIgAajh76XnbFO6qoHhStcHOaD6z/IFYuuoC/VRywZ47Lay6YkUAx2RDBrGKuHM54e0D0X30NxbTEnXj5BJp3B5XGRyWSIHosSfn2Y3uNaAlscoi6ijhjnvPUcos1RNt46ZqfCMgv59GOfpifec9IEsLa+Ngq9hXziok8MW/Ihloqx49gO2vvbcTld9CZ6cYmLpEniFCfJTBJjDIl0gqAnSDKdpLqwGoc4CBfpJLCvvvWrEzZwo7X5dLOVporpqik03hGBjRHMEsYy8uPpARVWF2oaa20hHXs7SJEinUzj8DiIdegEudw8iVR/Cm+x94znIVjGZjoD96OVhBgaQL1/9/10xbvo7O+kP92Pz+WjP9lPZ6wTj1MnYybSCZ1Zn8lWx032kzZpmnubWRJagtfp5XD3YR4+8PCEyyhMtIzFdDKZweizgXUNnSUi9REe+/Rj/Ox9P+OxTz9GpD4y3U06hQtuvoBYRwxBKF5eTCaZIdGdoKquisVXLqbqwiri0Th9bX0kogkKFhYMlMG2TA658tvHdx6n8flGXvyvF/nxO37M77/2+ym5fq4kRD4jlYRo6W0hmU7icDjwu/wkMgl6k72ICGmT1pEABgykSePAocKAAwQEQRxCub+c7U3bp6TNlvFhheAskPsxx3viFC8qJt4TZ8udW2acGKy+fjVX/u2VeAo9pPpSLLxwIe/4wTv4wEMf4OJPXExxbTEly0rwFnoJLQ1RXFtsA8WTzO7Nu4l1xWh+sZlMKoOvVGM4W76yhT0P7pn0659OSYiKYAUuh4uMyeB3+XGKk7RJ43F4cIoTpzhxiIOMyeB1evG5fDpZCkPIEyJlUlpoLbgQIxN3SU+0jIVlZKxr6CwwU2rpj4eR0lgnuqCM5czobuym42AHLr8G70HXoYi1x9hx744JpxyPl/wUzbFKQtRV1/HwgYfxOr24xMXRrqP4nD7KgmX0JnoHYgNd8S68Li8LAgsQEQxaTsHv8rOoeBHxdJwrwldMSZst48MKwVlgptTSt8w+imqKONR1aGAkAJBJZvCGvPQ09YzyyrPH0HpCIxEuDnPLBbfw7e3fZqFrIevNejr6O2jrb6MsUKallvtaKfQW4nV6McawMLiQhQU6AqgtqiXgCRCSEO397Xz0wY/SG++lwFPAktCSMev/T6TNM4WZsgDNSFghOE2GC+zNpFr6ltnFuTeey87/3km8R8ucZ5IZ0vE0wQVBCqsLxz7BFFMXrhtYs3hv614OdR6ipqCGE30n2H5sO8W+YlaUrCCeidPe386Gqg2cU3oOgpAxGUSEA20H2NWyC4/DQ2NPI6lMiu5Y95j1/2crM2kBmpGYt0LQFeniwMMHaNrehBihamPVwAIxIzHSurrLrl3GoYcPAbaWvuX0CNeFueQvLmHLV7YQa4/hDXkJLtAsrQtuvmCaWzc8+Rkw+T3dCxZeMOZC7ffvvp9dLbsoD5RztOsoRb4iMNCT6hmz/v9sZbhyHbntM+U+57QQ5Jae7G3RzAZEl5MUh+gSky39OsFK4MCjB9j/4H76O/vpeq2LjGQoX1XORR+/aMBPO1IsoH13O5fefim7N++m62gXRTVFrP/Q+hkXH7DMTC7/88spX1k+UIuqsLrwtGtRTRenmxbZ0ttCIpWgxF9CX7KPoFtFL5qIDtT/n2itopnKcOU6Ztp9zlkhyF96UpxCw9MNGGNY/AeLOfHyCY7vPE7pilI8BR7i3XF6mnroPNJJojuBN+TFIQ7a97XzxGefADTImh8LaOnysOdYER1RF56+Ls79RJi3fMUafsvEOJNaVGdKpCvCIwceYduxbYgRLqy+cNLq4FcEK/C4PMRSMQLuAIlMAowu6zie+v+zkeHWWphp9zlnhSB/6ckTu07gL9OVhzoOdpBOphGn0N3YTawjRsfBDmLdMWIdMdw+N74CH5mUpr2JCD/96ms0/e9qHn3kTfTGXPi8GYoCKTYsbac22E6iMMgDD+j6sBNdFtBimQ6Gq3EzmXXw66rrePXEq+xr30fIG+Jgx0FSmRS1hbUsDC48a/X/ZxLDrbUw0+5zzs4jyF96Mt4Vx+V34fK5iHXF8BX7tKjTgQ4yyQwItPX5eTm2nK2J9ezsXkpnugCTMRzpX8B9L67i2WehlyCpNLT0eGju9PLs7lKONLlZ/cZqQiFdJNximU3UN9XT0t9CeaCcoCdIgaeAMn/ZQB38s01+TR2n00lNYQ1rK9byuoWvo7a4dkYFUM8W+UtaDre850xgzo4I8pee9BZ7SfWnAPAV+yhdUcrhJw8jTs1k6EgGeKWvCqczSjDTQzxTxEs9S1gXPMT26FLE66G3FwpDbhyF0NycJpEU0n6htXQl5SuCZDLQPHNcfhbLuMj32efwu/yTWgc/XBzmIxs+wkfm0fJ2tsTENFFdV02sM0a8O07JOSX0t/XT19pHyTklOD1OghVBqjZUkYlnaEzXECoXQgscmHQGZ7IXbzrGkb4KOhNBQouKSCTA7QaX301BiQ9nMECotoQ+kw12RaFi5rj8LJZxke+zzzGVdfAtM4M5OyLIXzQm0Zug9rJaEMikM7gDblZevxJ3wI23yMuhR4pZ6uylbU8LHQc6dH1iZxxTXsn555YTiRbi6YJkEjwecDrB4YBEAkpLobsbOjvhypnj8rNYxkW+z35oHfz8MsmWuc2cFQIYvWxCLqsIoKQwRU+3i9JlpZz3J+cRLA/S3Q2BgAZ/v/hFCAahpQXicUiltPfv8cCqVXrclVfaQLFl9pHz2ednDV1We9mkZQ1ZZiaTth6BiPwT8DYgARwEPmSM6czu+yxwM5AG/swY8+hI54HJW48gN8/g0N4Ezx1aQO26IioXB4lGtYefywKqr4fvfAe2bNHtJSVwySXw4Q9Dne00WSyWGcp41yOYTCF4K/AbY0xKRL4CYIz5tIisAX4EXARUA48DK40x6ZHONRUL00QiavBbWrS3X1dne/gWi2V2M+0L0xhjfp339DkgVyN2E/BjY0wcOCwiB1BR2DJZbRkP4bA1/BaLZX4yVTGCDwM/yT6uQYUhRyS77SRE5FbgVoBFixYN3W2xWCxjEumK8PCBh9netB0jho1VGzm/8nwiPZEZWwl0Ojij9FEReVxEXhnmb1PeMZ8DUsAPT+fcxph7jDEbjTEbK2xepsViOU0iXRG+++J3eabhGXxuHwFXgEcOPMIXn/oiDV0NVBZU0pfs44G9DxDpmlmLSE01ZzQiMMa8ebT9InITcD3wJjMYjGgEavMOC2e3WSwWy1mjvqme1v5WyvxlBNxaIDKRSZAwCY73HmdpydIZWQl0Opg015CIXAP8JXClMaYvb9cvgf8Rka+hweIVwNbJaofFYpmf5GZNl/pLB7al0ikcxkFXbHDRqNOuBFpfD5s3Q2Mj1NTAjTfO+vTByYwR/CvgBR4TEYDnjDEfNcbsEpH7gFdRl9HHRssYslgslomQmzXdn+ofGBG4nC7i6TjFvsH5RaNWAh2aTuh2ww9+AGVlsGiR5pN/8Ytw+eWaVz5LUw4nM2to+Sj7vgx8ebKubbFYLHXVdew6sYsDHQd0dTQEj8ODyegSmhmTGb0SaCQCDzwAoRBUVmodme9+FxYuVGF44gno6tISA01N8PnP6zGzsBTxnJ5ZbLFY5i/h4jA3rb/ppKyha5ZfM5A11BxtpiJYwZWLrxw+PlBfryJQlF1HoKgI+vpg/344cgS8Xq01IwJ798K2bXDRRYOvtUJgsVgs00+4OMwtG26BDSdvr2McPv2WFh0J5FNZCU8+qaIQCGjNGbdbC5A99pgKQUHBrCtFPGerj1osFssZUVGhrp581q9XV5DDAcboXzKp7qL2dj1mFpYitkJgsVgsw1FXp8Hg7m7IZPR/WRmcd55Wnuzr02qUFRWQTmuwOFeKeJZlEVnXkMVimTOc1Zph4bAGfevr1dVTUaFlhhcvhjvuUMMfCsGxY/p35ZWzthSxFQKLxTInGC7JZyCBhwg88ogGdEXgwgvh2mvHNtjDFSHLPb/3XmhogOpq+Oxn4frrJ+W+poJJqz56NpmK6qMWi2UWE4lw/zca6Gvpo2hhAFasgPJyug+3Edj9Ajc0fVPdNsuXazC3owMWLNDjjJm1+f9jMd7qozZGYLFYZjfZoUBLq4OChUHN5Hn+edi/j4Jdz9Gyr119/H6/9uAzGV1VautW2LVLhw99fTp8iMzPmkPWNWSxWGY32Xz/ikohGndTFMj2b595lmjlCircEQ3mFhQMTv4yBlwuzfhxOAbnCpxJ/v8sXtTECoHFYpndZPP961Z08cDzmvdf4PURPd5HZ3kBVy7vgKMuFQGPB3p79bHfD8V5S9kOzf83RuMJQ8kZ/BMn1L2UyxAaLkBx0UV6fL44wIwTDCsEFotldpPN9w+XO9h0cTP1+4tpbhYqFgS5cu1BwoWV0FWsReK8Xh0JpFLg82mMIEd+/v8XvqBpoHfdpa/LGW4ROHAAYjH47W8Hs4S6u9UltXChppi2tcHBg/CLX2gdov5+dVf99KdaqG79+mEi2tMnBlYILBbL7OTBBzVz5/Bh7b2/5S2Er7iC8NoTUNOpvfGtW8ETgquugmefhX37YMkSuOKKwRFCJsPAQuVXXqnnikT03Nu2wbnnqpHu74enn9ZRg9Op51q3TrORYjFYtgyOHlXBWLAAWlv1vF1davirqvT1L72k5wiHVTQaGuBLX4Krr5620YEVAovFMvt48MHBXP7ly+H4cfj5z9Ugv+lNg7n8VVXam+/t1fTOfEObc/HkzxEIh3X7+vVqmB99VHv3paWwdq0GlTs7dduSJSoiuRnFqZSOFtrb1fj7/RqDaGrS1y1erPMN3G593toKv/sdrF6tgpQLWE/D6MAKgcVimX3ce6+KQFmZPq+pUVdPYyPccMPgcaMtRj7SvkceUYO+bJka7z179PzptG6PRtUlVFiocYJEQvft3q0CAGrkPR4VB7dbBcHphJ4edU/t3auvB9i5U+8jFFLX0jQUrLNCYLFYZh9NTVBbe/K2UEjdLGdCJKLlpXNzDQoL1eg3NuofqEHPZHQUkkzqnzE6EshlIxmjIpBMalwhHteKpTmKi/V5NKoCtnatHrNrl7qgphgrBBaLZfZRXa0umtyIAPR5dfWZnbe+Xg3/wYOaUrp4sQpAX94iiyIqBLnevcOhz0UGBcDpHBSIVEqPS2fX38oFq3OjBL9f3VNFRfo852qawnRUO6HMYrHMPm6+WXvsbW1qYNva9PnNN5/ZeVtatIeeSxv1+9XI5+PzDfb6/X4dNeRGBum0CgOoyyiVUrEKh3V7MqlxjHAYystVdAIBdSO99tpg8bpcvYy+vimZ8GZHBBaLZfaRq+uTX+/nk58883o/FRXqmlmzRgO727ap3z8QUMMdi2lcIBjUXn4iodeOxdRYi6g45Hr9Tqcek8moOyj3uKFBjwuF9Lq50cPateryGm5RHJi0+IEVAovFMju5/vqzX+itrg4efliN+NGj2ktfvVrFIZehtH07HDqkQgA6Ili3ToPFuZhAbkGb4mJ19XR363Hl5bq/rEwFob9f01ODQRUCp1Pb8NBDpy6KM4kL3ljXkMViseQIh+EjH4Hf/1573xs2wEc/Cuefry6hzk5NLV23ToO76bQa+kRC4wl+vwrGypXwB3+gLh+PRw3/ZZep0V+zRkcP6bSmigYCauDXrh1MHR1uUZxJXPDGjggsFosln4sugre8RQ37296mk8P++q91TkFrqxr94mJ1EeXmF3R1Dc4sXrBAhSAWUxE4dkzPW14+UBWVw4c1CF1UpPMXhgaC6+o0JgA6Esif8DYJTHoZahH5C+BOoMIY0yoiAtwNXAf0ATcZY7aPdg5bhtpisUw5Q2sN5WfxHDyoPfljx9St4/NpsPrYMR0V+Hzq2qmqUreQiE5AyzfqY00cOwtZQ+MtQz2pIwIRqQXeChzN23wtsCL7dzHw79n/FovFMnMYWnAufwLaPfeooV+0CPbvV8NeXKyjiPe+91QDDsPPYh6N0SbDnWUm2zV0F/CXwAN52zYB3zc6FHlOREIiUmWMOTbJbbFYLJazQ86HX16uf6A9/0BgZAM+g0tST1qwWEQ2AY3GmJ1DdtUA+dP/ItltQ19/q4hsE5FtLS0tk9VMi8ViOX2GW9h+Fi5an+OMRgQi8jiwcJhdnwP+CnULTQhjzD3APaAxgomex2KxWM46Iy1sP4N7/aNxRkJgjHnzcNtF5DxgKbBTY8OEge0ichHQCOQXCQlnt1ksFsvsYQp9+JPNpLiGjDEvG2MWGGOWGGOWoO6fC40xx4FfAh8Q5RKgy8YHLBaLZfqYjnkED6GpowfQ9NEPTUMbLBaLxZJlSoQgOyrIPTbAx6biuhaLxWIZG1tiwmKxWOY5VggsFotlnmOFwGKxWOY5VggsFotlnmOFwGKxWOY5VggsFotlnmOFwGKxWOY5VggsFotlnmOFwGKxWOY5VggsFotlnmOFwGKxWOY5VggsFotlnmOFwGKxWOY5VggsFotlnmOFwGKxWOY5VggsFotlnmOFwGKxWOY5VggsFotlnmOFwGKxWOY5VggsFotlnmOFwGKxWOY5VggsFotlnjOpQiAinxCRPSKyS0S+mrf9syJyQET2isjVk9kGi8VisYyOa7JOLCJvADYB5xtj4iKyILt9DfBuYC1QDTwuIiuNMenJaovFYrFYRmYyRwR/CvyjMSYOYIw5kd2+CfixMSZujDkMHAAumsR2WCwWi2UUJlMIVgJXiMjzIvKUiNRlt9cADXnHRbLbTkJEbhWRbSKyraWlZRKbabFYLPObM3INicjjwMJhdn0ue+5S4BKgDrhPRJaN99zGmHuAewA2btxozqSdFovFYhmZMxICY8ybR9onIn8K/NwYY4CtIpIByoFGoDbv0HB2m8VisVimgcl0Df0CeAOAiKwEPEAr8Evg3SLiFZGlwApg6yS2w2KxWCyjMGlZQ8B3gO+IyCtAAvhgdnSwS0TuA14FUsDHbMaQxWKxTB+TJgTGmATwvhH2fRn48mRd22KxWCzjx84stlgslnmOFQKLxWKZ51ghsFgslnmOFQKLxWKZ51ghsFgslnmOFQKLxWKZ51ghsFgslnmOFQKLxWKZ51ghsFgslnmOFQKLxWKZ51ghsFgslnmOFQKLxWKZ51ghsFgslnmOFQKLxWKZ51ghsFgslnmOFQKLxWKZ51ghsFgslnmOFQKLxWKZ51ghsFgslnmOFQKLxWKZ51ghsFgslnnOpAmBiKwXkedE5EUR2SYiF2W3i4h8XUQOiMhLInLhZLXBYrFYLGMzmSOCrwJ/a4xZD3w++xzgWmBF9u9W4N8nsQ0Wi8ViGYPJFAIDFGUfFwNN2cebgO8b5TkgJCJVk9gOi8VisYyCaxLP/UngURG5ExWc12e31wANecdFstuO5b9YRG5FRwwsWrRoEptpsVgs85szEgIReRxYOMyuzwFvAj5ljPmZiLwLuBd483jPbYy5B7gHYOPGjeZM2mmxWCyWkTkjITDGjGjYReT7wG3Zpz8F/jP7uBGozTs0nN1msVgslmlgMmMETcCV2cdvBPZnH/8S+EA2e+gSoMsYc2y4E1gsFotl8pnMGMEtwN0i4gJiZP39wEPAdcABoA/40CS2wWKxWCxjMGlCYIz5PbBhmO0G+NhkXddisVgsp4edWWyxWCzzHCsEFovFMs+xQmCxWCzzHCsEFovFMs+xQmCxWCzzHCsEFovFMs+xQmCxWCzzHCsEFovFMs+xQmCxWCzzHCsEFovFMs+xQmCxWCzzHCsEFovFMs+xQmCxWCzzHCsEFovFMs+xQmCxWCzznMlcmGZuE4lAfT20tEBFBdTVQTg83a2yWCyW08aOCCZCJAIPPAB9fVBZqf8feEC3WywWyyzDCsFEqK+HUAiKisDh0P+hkG63WCyWWYYVgonQ0gIFBSdvKyjQ7RaLxTLLsEIwESoqIBo9eVs0qtstFotllmGFYCLU1UFnJ3R3Qyaj/zs7dbvFYrHMMs5ICETkj0Vkl4hkRGTjkH2fFZEDIrJXRK7O235NdtsBEfnMmVx/2giHYdMmCASguVn/b9pks4YsFsus5EzTR18B3gl8K3+jiKwB3g2sBaqBx0VkZXb3N4G3ABGgXkR+aYx59QzbMfWEw9bwWyyWOcEZCYExZjeAiAzdtQn4sTEmDhwWkQPARdl9B4wxh7Kv+3H22NknBBaLxTJHmKwYQQ3QkPc8kt020naLxWKxTBNjjghE5HFg4TC7PmeMeeDsN2ngurcCtwIsWrRosi5jsVgs854xhcAY8+YJnLcRqM17Hs5uY5TtQ697D3APwMaNG80E2mCxWCyWcTBZrqFfAu8WEa+ILAVWAFuBemCFiCwVEQ8aUP7lJLXBYrFYLOPgjILFInID8A2gAvhfEXnRGHO1MWaXiNyHBoFTwMeMMensaz4OPAo4ge8YY3aNdZ0XXnihVUSOnElbJ0A50DrF15xM7P3MfObaPdn7mX4Wj+cgMcZ6XYZDRLYZYzaOfeTswN7PzGeu3ZO9n9mDnVlssVgs8xwrBBaLxTLPsUIwMvdMdwPOMvZ+Zj5z7Z7s/cwSbIzAYrFY5jl2RGCxWCzzHCsEFovFMs+xQjAEEfk7EXlJRF4UkV+LSHV2u4jI17Pls18SkQunu63jQUT+SUT2ZNt8v4iE8vYNWyp8JjOR0ucznblQml1EviMiJ0TklbxtpSLymIjsz/4vmc42jhcRqRWR34rIq9nv2m3Z7bPyfsaDFYJT+SdjzOuMMeuBB4HPZ7dfi86QXoHWQPr36WneafMYsM4Y8zpgH/BZOKVU+DXAv4mIc9paOX5ypc9/l79xtt5Pto3fRL9fa4A/yd7LbOO76Puez2eAJ4wxK4Anss9nAyngL4wxa4BLgI9lP5PZej9jYoVgCMaY7rynQSAXTd8EfN8ozwEhEama8gaeJsaYXxtjUtmnz6H1nSCvVLgx5jCQXyp8xmKM2W2M2TvMrll5P2gbDxhjDhljEkCuNPuswhjzO6B9yOZNwPeyj78HvGMq2zRRjDHHjDHbs497gN1oleRZeT/jwQrBMIjIl0WkAXgvgyOCuVBC+8PAw9nHc+F+8pmt9zNb2z0eKo0xx7KPjwOV09mYiSAiS4ALgOeZA/czEme6QtmsZKzS2saYzwGfE5HPAh8H7pjSBp4m4ykVLiKfQ4e8P5zKtk2E6Sp9bpk8jDFGRGZVrrqIFAA/Az5pjOnOX4BrNt7PaMxLITiN0to/BB5ChWC00trTylj3IyI3AdcDbzKDE0dm7f2MwIy9nzGYre0eD80iUmWMOZZ1o56Y7gaNFxFxoyLwQ2PMz7ObZ+39jIV1DQ1BRFbkPd0E7Mk+/iXwgWz20CVAV94wccYiItcAfwm83RjTl7drpFLhs5XZej9zuTT7L4EPZh9/EJgVoznRrv+9wG5jzNfyds3K+xkPdmbxEETkZ8AqIAMcAT5qjGnMfjn+Fc2M6AM+ZIzZNn0tHR/Z9aK9QFt203PGmI9m930OjRuk0OHvw8OfZeYwpPR5J/CiMebq7L5Zdz8AInId8C8Mlmb/8vS26PQRkR8BV6GlmpvRUfQvgPuARehv6V3GmKEB5RmHiFwOPA28jNoBgL9C4wSz7n7GgxUCi8VimedY15DFYrHMc6wQWCwWyzzHCoHFYrHMc6wQWCwWyzzHCoHFYrHMc6wQWCwWyzzHCoHFYrHMc/5/SKLpCDe8EKIAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "tsne = TSNE(\n",
    "    n_components=2, perplexity=15, random_state=42, init=\"random\", learning_rate=200\n",
    ")\n",
    "vis_dims2 = tsne.fit_transform(matrix)\n",
    "\n",
    "x = [x for x, y in vis_dims2]\n",
    "y = [y for x, y in vis_dims2]\n",
    "\n",
    "for category, color in enumerate([\"purple\", \"green\", \"red\", \"blue\",\"yellow\"]):\n",
    "    xs = np.array(x)[embedding_df.Cluster == category]\n",
    "    ys = np.array(y)[embedding_df.Cluster == category]\n",
    "    plt.scatter(xs, ys, color=color, alpha=0.3)\n",
    "\n",
    "    avg_x = xs.mean()\n",
    "    avg_y = ys.mean()\n",
    "\n",
    "    plt.scatter(avg_x, avg_y, marker=\"x\", color=color, s=100)\n",
    "plt.title(\"Clusters identified visualized in language 2d using t-SNE\")\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Cluster 0 Theme:  Literary & Archival ItemsOne possible way to group these transactions into a meaningful cluster is by their theme or purpose. A common theme that appears in most of these transactions is literary and archival items, which suggests that these are payments for acquiring, preserving, or accessing various collections of books, manuscripts, papers, and other materials related to literature and history. This theme could be relevant for a library, a museum, a research institution, or a collector.Some possible indicators of this theme are:-\n",
      "JISC SERVICES LTD SUBSCRIPTION ACCOUNT, Annual Subscription\n",
      "   Private Sale, Literary & Archival Items\n",
      "   ALDL, Oct19-Dec19 charge from Agency for Legal Deposit Libraries\n",
      "   ALDL, Legal Deposit Services\n",
      "   Private Sale, Literary & Archival Items\n",
      "   SONYA LEONARD, Literary and personal papers of Tom Leonard 1961 to 2018\n",
      "   FROST AND SULLIVAN LTD, Literary & Archival Items\n",
      "   Cengage Learning (Emea )Ltd, Literary & Archival Items\n",
      "   Agency for the Legal Deposit Libraries, Agency services\n",
      "   Robert Harland, Correspondance and Literary papers - Thomas Carlyle\n",
      "   ----------------------------------------------------------------------------------------------------\n",
      "Cluster 1 Theme: There are different ways to approach the task of grouping transactions into meaningful clusters, but one possible theme that these transactions have in common is:- They are all payments to the same supplier, ECG Facilities Service, for various services related to facility management and maintenance.- They are all relatively large amounts, ranging from £27,013.16 to £125,000.00, with a mean of £54,798.67 and a standard deviation of £19,894.67.-\n",
      "ECG FACILITIES SERVICE, Maintenance contract - all properties\n",
      "   ECG, This payment covers 16 invoices including upgrade to boiler control panel & remedial works following 5 year test\n",
      "   ECG FACILITIES SERVICE, This payment covers multiple invoices for facility management fees\n",
      "   ECG FACILITIES SERVICE, CB Bolier Replacement (1),USP Batteries,Gutter Works & Cleaning of pigeon fouling\n",
      "   ECG FACILITIES SERVICE, Facilities Management Charge\n",
      "   ECG FACILITIES SERVICE, Facilities Management Charge\n",
      "   ECG Facilities Service, Facilities Management Charge\n",
      "   ECG FACILITIES SERVICE, Maintenance contract - August 21 period\n",
      "   ECG Facilities Service, Facilities Management Charge\n",
      "   ECG Facilities Service, Facilities Management Charge\n",
      "   ----------------------------------------------------------------------------------------------------\n",
      "Cluster 2 Theme: There are different ways to approach the task of grouping transactions into meaningful clusters, but one possible theme that these transactions have in common is:- They are all related to building or refurbishment projects at different locations in Edinburgh- They are all paid to either John Graham Construction Ltd or Arthur McKay Building Services, two contractors that provide similar services- They are all relatively large and consistent in value, ranging from around £27,000 to £125,000 per transaction\n",
      "John Graham Construction Ltd, Causewayside Refurbishment\n",
      "   John Graham Construction Ltd, Causewayside Refurbishment\n",
      "   SJS Property Services, George IV Bridge Work\n",
      "   John Graham Construction Ltd, Causewayside Refurbishment\n",
      "   John Graham Construction Ltd, Causewayside Refurbishment\n",
      "   John Graham Construction Ltd, Causewayside Refurbishment\n",
      "   John Graham Construction Ltd, Causewayside Refurbishment\n",
      "   ARTHUR MCKAY BUILDING SERVICES, Causewayside Work\n",
      "   John Graham Construction Ltd, Causewayside Refurbishment\n",
      "   John Graham Construction Ltd, Causewayside Refurbishment\n",
      "   ----------------------------------------------------------------------------------------------------\n",
      "Cluster 3 Theme: There are different ways to approach the task of grouping transactions into meaningful clusters, and different criteria or methods may yield different results. However, one possible theme that these transactions have in common is that they are all related to operational costs or capital expenditures for a large organization, such as a university, a hospital, or a government agency. Some possible sub-themes or categories within this theme are:- Electricity and utility services (EDF, EDF Energy)- IT work, equipment, and software\n",
      "EDF, Electricity\n",
      "   m-hance, IT Work\n",
      "   EDF Energy, This payment covers 5 invoices for utility services\n",
      "   JISC Services Ltd, Managed router service charge annual subscription 01/04/22 to 31/03/23\n",
      "   ALDL, ALDL Charges\n",
      "   COMPUTACENTER UK, IT equipment\n",
      "\n",
      "   BSI, Subscription\n",
      "   PHOENIX SOFTWARE LTD, IT Hardware plus 5 year licence\n",
      "   AM Phillip, Vehicle Purchase\n",
      "   XMA SCOTLAND LTD, Purchase of IT equipment and renewal of maintenance agreement.  This payment covers 2 invoices\n",
      "   ----------------------------------------------------------------------------------------------------\n",
      "Cluster 4 Theme: There are different ways to approach the task of clustering transactions, and different criteria to define meaningful clusters, but one possible theme that these transactions have in common is:- They are all related to property or facility costs, such as rent, rates, service charges, or maintenance.- They are all paid to either Glasgow City Council, Glasgow Life, or City of Edinburgh Council, which are public sector organisations that manage or provide various services in Scotland.- They are all relatively large amounts, ranging from\n",
      "GLASGOW CITY COUNCIL, Kelvin Hall\n",
      "   XMA Scotland Ltd, Kelvin Hall\n",
      "   Creative Video Productions Ltd, Kelvin Hall\n",
      "   GLASGOW LIFE, Oct 20 to Dec 20 service charge - Kelvin Hall\n",
      "   CITY OF EDINBURGH COUNCIL, Rates for 33 Salisbury Place\n",
      "   Glasgow Life, Service Charges\n",
      "   City Of Edinburgh Council, Non Domestic Rates \n",
      "   Glasgow Life, Service charges\n",
      "   Glasgow Life, Kelvin Hall Service Charges\n",
      "   Glasgow City Council, Kelvin Hall\n",
      "   ----------------------------------------------------------------------------------------------------\n"
     ]
    }
   ],
   "source": [
    "# We'll read 10 transactions per cluster as we're expecting some variation\n",
    "transactions_per_cluster = 10\n",
    "\n",
    "for i in range(n_clusters):\n",
    "    print(f\"Cluster {i} Theme:\", end=\" \")\n",
    "\n",
    "    transactions = \"\\n\".join(\n",
    "        embedding_df[embedding_df.Cluster == i]\n",
    "        .combined.str.replace(\"Supplier: \", \"\")\n",
    "        .str.replace(\"Description: \", \":  \")\n",
    "        .str.replace(\"Value: \", \":  \")\n",
    "        .sample(transactions_per_cluster, random_state=42)\n",
    "        .values\n",
    "    )\n",
    "    response = openai.Completion.create(\n",
    "        engine=COMPLETIONS_MODEL,\n",
    "        # We'll include a prompt to instruct the model what sort of description we're looking for\n",
    "        prompt=f'''We want to group these transactions into meaningful clusters so we can target the areas we are spending the most money. \n",
    "                What do the following transactions have in common?\\n\\nTransactions:\\n\"\"\"\\n{transactions}\\n\"\"\"\\n\\nTheme:''',\n",
    "        temperature=0,\n",
    "        max_tokens=100,\n",
    "        top_p=1,\n",
    "        frequency_penalty=0,\n",
    "        presence_penalty=0,\n",
    "    )\n",
    "    print(response[\"choices\"][0][\"text\"].replace(\"\\n\", \"\"))\n",
    "\n",
    "    sample_cluster_rows = embedding_df[embedding_df.Cluster == i].sample(transactions_per_cluster, random_state=42)\n",
    "    for j in range(transactions_per_cluster):\n",
    "        print(sample_cluster_rows.Supplier.values[j], end=\", \")\n",
    "        print(sample_cluster_rows.Description.values[j], end=\"\\n   \")\n",
    "        #print(str(sample_cluster_rows.Value).values[j])\n",
    "\n",
    "    print(\"-\" * 100)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Conclusion\n",
    "\n",
    "We now have five new clusters that we can use to describe our data. Looking at the visualisation some of our clusters have some overlap and we'll need some tuning to get to the right place, but already we can see that GPT-3 has made some effective inferences. In particular, it picked up that items including legal deposits were related to literature archival, which is true but the model was given no clues on. Very cool, and with some tuning we can create a base set of clusters that we can then use with a multiclass classifier to generalise to other transactional datasets we might use."
   ]
  }
 ],
 "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.7.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}