B_cells_versus_T_cells.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## B cells versus T cells\n",
    "The data used in this notebook is lymphocyte data for one patient's B cells and T cells. We use this data to show the proteome variation between the cell types. Here, we note presence/absense changes and quantify fold changes using a T test.\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "First, we import our loader module. This brings the functions defined there (in our repository at ~/load_data.py) into scope so we can use them in this script. Then we can load our data and store it as <code>data</code>.\n",
    "\n",
    "Calling <code>head</code> shows the first several lines of the dataframe, which provides an idea of the type of data present and the structure of the dataframe."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Intensity B_C10</th>\n",
       "      <th>Intensity B_C11</th>\n",
       "      <th>Intensity B_C12</th>\n",
       "      <th>Intensity B_C13</th>\n",
       "      <th>Intensity B_C9</th>\n",
       "      <th>Intensity T_D10</th>\n",
       "      <th>Intensity T_D11</th>\n",
       "      <th>Intensity T_D12</th>\n",
       "      <th>Intensity T_D13</th>\n",
       "      <th>Intensity T_D9</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Protein IDs</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>P0DPI2;A0A0B4J2D5</th>\n",
       "      <td>265300</td>\n",
       "      <td>248290</td>\n",
       "      <td>1878900</td>\n",
       "      <td>289210</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>654710</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>A0AVT1</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>440250</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>A6NGU5;P19440</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>2590400</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>A6NHQ2</th>\n",
       "      <td>126410</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>89903</td>\n",
       "      <td>0</td>\n",
       "      <td>263160</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>110120</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>A6NI72;P14598;A8MVU1</th>\n",
       "      <td>1146300</td>\n",
       "      <td>1354400</td>\n",
       "      <td>2254700</td>\n",
       "      <td>2912100</td>\n",
       "      <td>1045700</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9Y5J7</th>\n",
       "      <td>313760</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>259020</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>237540</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9Y5S9</th>\n",
       "      <td>766910</td>\n",
       "      <td>415980</td>\n",
       "      <td>458400</td>\n",
       "      <td>452850</td>\n",
       "      <td>862000</td>\n",
       "      <td>343680</td>\n",
       "      <td>495150</td>\n",
       "      <td>551340</td>\n",
       "      <td>766580</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9Y5X3</th>\n",
       "      <td>291730</td>\n",
       "      <td>748590</td>\n",
       "      <td>1022100</td>\n",
       "      <td>133040</td>\n",
       "      <td>510570</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9Y6C9</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>311590</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>416530</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9Y6N5</th>\n",
       "      <td>0</td>\n",
       "      <td>121830</td>\n",
       "      <td>440970</td>\n",
       "      <td>0</td>\n",
       "      <td>524620</td>\n",
       "      <td>0</td>\n",
       "      <td>1352300</td>\n",
       "      <td>660590</td>\n",
       "      <td>193530</td>\n",
       "      <td>1067300</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>1211 rows × 10 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                      Intensity B_C10  Intensity B_C11  Intensity B_C12  \\\n",
       "Protein IDs                                                               \n",
       "P0DPI2;A0A0B4J2D5              265300           248290          1878900   \n",
       "A0AVT1                              0                0           440250   \n",
       "A6NGU5;P19440                       0                0                0   \n",
       "A6NHQ2                         126410                0                0   \n",
       "A6NI72;P14598;A8MVU1          1146300          1354400          2254700   \n",
       "...                               ...              ...              ...   \n",
       "Q9Y5J7                         313760                0                0   \n",
       "Q9Y5S9                         766910           415980           458400   \n",
       "Q9Y5X3                         291730           748590          1022100   \n",
       "Q9Y6C9                              0                0           311590   \n",
       "Q9Y6N5                              0           121830           440970   \n",
       "\n",
       "                      Intensity B_C13  Intensity B_C9  Intensity T_D10  \\\n",
       "Protein IDs                                                              \n",
       "P0DPI2;A0A0B4J2D5              289210               0                0   \n",
       "A0AVT1                              0               0                0   \n",
       "A6NGU5;P19440                       0               0                0   \n",
       "A6NHQ2                          89903               0           263160   \n",
       "A6NI72;P14598;A8MVU1          2912100         1045700                0   \n",
       "...                               ...             ...              ...   \n",
       "Q9Y5J7                              0               0                0   \n",
       "Q9Y5S9                         452850          862000           343680   \n",
       "Q9Y5X3                         133040          510570                0   \n",
       "Q9Y6C9                              0               0                0   \n",
       "Q9Y6N5                              0          524620                0   \n",
       "\n",
       "                      Intensity T_D11  Intensity T_D12  Intensity T_D13  \\\n",
       "Protein IDs                                                               \n",
       "P0DPI2;A0A0B4J2D5                   0                0                0   \n",
       "A0AVT1                              0                0                0   \n",
       "A6NGU5;P19440                       0                0          2590400   \n",
       "A6NHQ2                              0                0           110120   \n",
       "A6NI72;P14598;A8MVU1                0                0                0   \n",
       "...                               ...              ...              ...   \n",
       "Q9Y5J7                         259020                0                0   \n",
       "Q9Y5S9                         495150           551340           766580   \n",
       "Q9Y5X3                              0                0                0   \n",
       "Q9Y6C9                         416530                0                0   \n",
       "Q9Y6N5                        1352300           660590           193530   \n",
       "\n",
       "                      Intensity T_D9  \n",
       "Protein IDs                           \n",
       "P0DPI2;A0A0B4J2D5             654710  \n",
       "A0AVT1                             0  \n",
       "A6NGU5;P19440                      0  \n",
       "A6NHQ2                             0  \n",
       "A6NI72;P14598;A8MVU1               0  \n",
       "...                              ...  \n",
       "Q9Y5J7                        237540  \n",
       "Q9Y5S9                             0  \n",
       "Q9Y5X3                             0  \n",
       "Q9Y6C9                             0  \n",
       "Q9Y6N5                       1067300  \n",
       "\n",
       "[1211 rows x 10 columns]"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import load_data\n",
    "#data_raw = load_data.load_FragPipe()\n",
    "data_raw=load_data.load_max_quant()\n",
    "\n",
    "data_raw"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "These are some methods that assist with the data analysis, such as checking if a value was measured or not and if the protein was measured in each type."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "import data_utils"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Before any analysis, we log normalize and median normalize across runs."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "data_log2_normalized = data_utils.normalize(data_raw)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Our next steps will require separating the runs by cell type, to compare B cells to T cells. As you see above, the pattern for naming the columns is \"Intensity \", B or T indicating cell type, and the replicate number. Here we define the list of cell types that the functions should look for based on the naming pattern."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "cell_types=[' B_', ' T_']"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Here we record which proteins are only identified in one cell type. We'll set these aside for now, in the <code>data_in_one_type</code> dataframe. Note that <code>min_reps</code> is the minimum number of replicates required to list a protein as present in one type. This can be changed, i.e., <code>min_reps=5</code> to determine the proteins identified in all 5 replicates of a single type; do not set it to a number greater than the replicates run or no proteins will be returned."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Proteins identified only in:\n",
      " T_    49\n",
      " B_    38\n",
      "dtype: int64\n"
     ]
    }
   ],
   "source": [
    "present_only_in = data_log2_normalized.apply(data_utils.check_presence_absence, axis=1, \n",
    "                                             min_reps=3, cell_types=cell_types).dropna()\n",
    "print(\"Proteins identified only in:\")\n",
    "print(present_only_in.value_counts())\n",
    "data_in_one_type =data_log2_normalized.loc[present_only_in.index]\n",
    "data_in_one_type['Type'] = present_only_in"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Next, we select the proteins that are measured in at least three samples from each group, allowing the calculations to proceed without imputed zero-handling."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "indecies = data_log2_normalized.apply(data_utils.check_three_of_each_type, axis=1, cell_types=cell_types)\n",
    "data = data_log2_normalized[indecies]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Total groups identified:\n",
      "1211\n",
      "Proteins identified in at least 3 of each cell type:\t\n",
      "589\n"
     ]
    }
   ],
   "source": [
    "print (\"Total groups identified:\")\n",
    "print(data_raw.shape[0])\n",
    "\n",
    "print (\"Proteins identified in at least 3 of each cell type:\\t\")\n",
    "print (data.shape[0])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now we calulate p values using a ttest. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "D:\\Anaconda3\\lib\\site-packages\\scipy\\stats\\_distn_infrastructure.py:1932: RuntimeWarning: invalid value encountered in less_equal\n",
      "  cond2 = cond0 & (x <= _a)\n"
     ]
    }
   ],
   "source": [
    "#Call ttest function\n",
    "tstats = data.apply(data_utils.ttest_wrapper, axis=1, cell_types=cell_types)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Here we calculate fold change."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "tstats['log2(B)-log2(T)'] = data.apply(data_utils.get_fold_changes, axis=1, cell_types=cell_types)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Fold change and p values allow us to create a volcano plot."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXgAAAEWCAYAAABsY4yMAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3de3yU1bXw8d+amcxkcoEESJAKilVEUxvEUEXoab20VivVWlCrgvVWoNTXtqfeeizHcw56XhU9nnoqBaxXvFSKWq3Wausr9RRvBakUg0gRLShIwARyn2RmvX/MhZlkJplc5pr1/Xz4lMw88zx7qFnPftbee21RVYwxxuQfR6YbYIwxJjUswBtjTJ6yAG+MMXnKArwxxuQpC/DGGJOnLMAbY0yesgBv8oqIfCAiXwn9/d9E5OFMtylVROQSEflz1M8qIkdksk0mu1iAN2kRCrytItIkIvUi8pyIjMt0u1JBRJaGvmeTiPhEpCPq5+cz3T4zdFiAN+n0DVUtAcYAnwD/k+H2pISqzlfVktB3/U/g8fDPqnpGpttnhg4L8CbtVLUNWAVUJTpGREaIyP0i8nGox/+bqPdmiMhfRaRBRF4VkererikihSLysIjsDX3uLyIyOs5x14vIqi6v/UxE7gr9/RIReV9EGkVkm4hc1Jfv3kP7zg59p/0islVETg+9PlxE7hWRnSLykYjcJCLOJM73dRGpDbXzIxG5ejDaaXKLK9MNMEOPiBQB5wOv93DYCqAJ+Fzof6eFPnsccB/wDWAtMBt4RkQmqmp7D+f7DjAcGAe0A8cCrXGOewz4VxEZpqr7Q8H0POAcESkG7gK+oKqbRWQMMCLJr52QiBwPPATMAl4i+IRTGnr7QYJPO0cAxcCzwHZgWS+nvRc4T1X/V0TKgcMG2k6Te6wHb9LpNyLSAOwHvgosjndQKHCeAcxX1XpV7VDVP4Xe/i6wTFXfUFW/qj5IMGBP7eXaHcBI4IjQ59ap6v6uB6nqh8BbwDdDL50CtKhq+GYUAI4REa+q7lTVd5L98j24HLhPVf+gqgFV/UhV3w09YZwB/FBVm1V1N3An8O0kztkBVIVuVPWq+tYgtNPkGAvwJp2+qaplgAe4EviTiBwU57hxwKeqWh/nvUOBH4fSLA2hG8Y44DO9XHsF8ALwq1Da5zYRKUhw7KPABaG/Xxj6GVVtJvjkMR/YGRooPqqX6yZjHLA1zuuHAgWha4W/6zKgMolzzgS+DnwoIn8SkRMHoZ0mx1iAN2kX6kE/CfiBL8Y5ZDswQkTKErx3s6qWRf0pUtXHerlmh6r+u6pWEUz3zAAuTnD4r4GTRGQscA6hAB86zwuq+lWCaZR3gXt6/rZJ2Q4cnuD1dmBU1Hcdpqqf6+2EqvoXVT2b4M3gN8DKQWinyTEW4E3aSdDZQDmwqev7qroTeB5YIiLlIlIgIl8KvX0PMF9ETgidp1hEzhSR0q7n6XLNk0Xk86Gc+n6CKQx/vGNVtQ5YDdwPbFPVTaFzjBaRs0K5+HaCYwNxz9FH9wKXisipIuIQkYNF5KjQv8OLwB0iMiz03uEi8uVevqtbRC4SkeGq2hH6voPRTpNjLMCbdPqtiDQRDDg3A9/pIYc9h2AQfhfYDfwQQFXXEszD/xyoB/4OXJLEtQ8iOHNnP8Gbyp+AnhZBPQp8hajeO8Hflx8DHwOfAl8GFiRx7R6p6pvApQTz6/tCbTs09PbFgBuoJfh9VxF8eujNHOADEdlPMKU0e6DtNLlHbMMPY4zJT9aDN8aYPGUB3hhj8pQFeGOMyVMW4I0xJk9lVamCUaNG6fjx4zPdDGOMyRnr1q3bo6oV8d7LqgA/fvx41q5dm+lmGGNMzhCRDxO9ZykaY4zJUxbgjTEmT1mAN8aYPGUB3hhj8pQFeGOMyVNZNYvGDB2BgLK32Yev04/b5WRksRuHQzLdLGPyigV4k3aBgLL5k0a++9BadtS3Mrbcyz0XT2Hi6FIL8sYMIkvRmLTb2+yLBHeAHfWtfPehtext9mW4ZcbkFwvwJu18nf5IcA/bUd+Kr9P2pDBmMFmAN2nndjkZW+6NeW1suRe3y5mhFhmTnyzAm7QbWezmnounRIJ8OAc/stid4ZYZk19skNX0yWDMfnE4hImjS3lqwXSbRWNMClmAN0kbzNkvDodQUepJUUuNMWApGtMHNvvFmNxiAd4kzWa/GJNbLMCbpNnsF2NyiwV4kzSb/WJMbrFBVpM0m/1iTG6xAG/6xGa/GJM7LMAbY/LeUK1eagHeGJPXhnL1UhtkNcbktaG8fsMCvDEmrw3l9RsW4I0xeW0or9+wAG+MyWtDef2GDbIaY/LaUF6/YQHeGJP3hur6jZQGeBH5AGgE/ECnqk5J5fWMMcYckI4e/MmquicN1zHGGBPFBlmNMSZPpTrAK/CiiKwTkbnxDhCRuSKyVkTW1tXVpbg5xhgzdKQ6wE9X1eOAM4Dvi8iXuh6gqstVdYqqTqmoqEhxc4wxZuhIaYBX1Y9D/7sbeAo4PpXXM8YYc0DKAryIFItIafjvwGnAxlRdzxhjTKxUzqIZDTwlIuHrPKqqv0/h9YwxxkRJWYBX1feBSak6vzHGmJ7ZNEljjMlTFuCNMSZPWYA3xpg8ZQHeGGPylFWTNGkzVDc+NiZTLMCbtBjKGx8bkymWojFpMZQ3PjYmUyzAm7QYyhsfG5MpFuBNWgzljY+NiRYIKHWN7XxU30JdYzuBgKbsWhbgTVoM5Y2PjQkLj0Wds2QN0299mXOWrGHzJ40pC/Kimrq7R19NmTJF165dm+lmmBSxWTRmqKtrbOecJWti0pVjy708tWB6v/eMFZF1ibZDtVk0Jm2G6sbHxoSleyzKUjTGGJMm6R6LsgBvjDFpku6xKEvRGGNMmjgcwsTRpTy1YDq+Tj8FLgcuh7BzX2tKxqUswBtjzAD0dfJAeCwqHau7LUVjjDH9NJBpj+lY3W0B3hhj+mkgQTodM2oswBtjTD8NJEinY0aNBXhjjOmngQTpdMyosZWsxhjTTwMdKB2M1d22ktUYY1Kg67THvgbpVK/utgBvjDEDkM0lOCwHb4wxecoCvDHG5CkL8MYYk6dSHuBFxCki60Xk2VRfyxhjzAHp6MH/ANiUhusYY4yJktIALyJjgTOBX6byOsYYY7pLdQ/+v4FrgUCiA0RkroisFZG1dXV1KW6OMcYMHSkL8CIyA9itqut6Ok5Vl6vqFFWdUlFRkarmGGPMkJPKhU7TgbNE5OtAITBMRB5W1dmDeRHbyNkYY+JLWYBX1Z8APwEQkZOAq1MR3FNdMN8YY3JVTs+DT0fBfGOMyVVpqUWjqquB1YN93nQUzDfGmFyV0z34dBTMN8aYXJXTAT4dBfNN/wUCSl1jOx/Vt1DX2J7UPpXGmMGT0+WCB1qL2aSODYAbk3k53YOHA7WYDy4voqLUY8EjS9gAuDGZl9M9eJO9+jIAbmsZjEkNC/AmJcID4NFBPt4AuKVyjEmdnE/RmOyU7AD4QFM5NpBrTGLWgzcpkewA+EDWMljv35ieWQ/epEwyA+ADWctgA7nG9MwCvMmogaxlsJXMxvTMUjQmowayliHRQK6IEAiopWnMkGc9eJNx/V3LEK/3f+vMav7tmY1s/qQx4YCrDcyaocJ68CZnhXv/K+edyMcNrext9nH7C5tZv72B2p2NPLVgOhWlnpjP2MCsGUoswJuc5nAIqsqspa/FvJ4oF59oYDbezcCYXGcpGpPz+jITxwZmzVBiAd7kvL7MxLES02YoEdXsGWCaMmWKrl27NtPNMDko2Xo2A8nBW80ck41EZJ2qTon3nuXgTV4Iz8RJ5rj+TMu0wVmTiyxFY4ac/kzLtFWzJhdZD96YBKJTMgAVJZ6YAVobnDXZzgK8GXKSyaXHS8ksnlXNbb8PzrMHG5w12c9SNCbnDGQlajhwn7NkDdNvfZlzlqyJu+o1XkrmmlUbuOrUCYDt/2tyg/XgTU4Z6GBnsgudEs2XP7yyhDXXnWyzaExOsB68ySkDHexMdqFTovny3gKn7f9rckavAV5ECkXkn0XkSRF5QkR+JCKF6WicMV0NdCVqsgudBlLG2JhskUyK5iGgEfif0M8XACuAc1PVKJPfBrJgKNm9XhMJB+6uKZ6ugXsgZYyNyRa9rmQVkbdVdVJvr8X5XCHwCuAheCNZpao39vQZW8ma/waaQx+MBUe93WACAWVPczttHX6cInjdTsq8FtxNdhroStb1IjJVVV8PnewEYE0Sn2sHTlHVJhEpAP4sIs+Hz2OGpoFWc3Q4hAkVJaycdyKd/gAup4PKkr7lw3ta9ZpoeuToYYWMH1lsQd7klGQGWU8AXhWRD0TkA+A14Msi8jcR2ZDoQxrUFPqxIPQnewrfmIwYaA49EFC21DVx3rLX+NLi1Zy37DW21DUN2qYdiaZHfri3xVatmpyTTA/+9P6eXEScwDrgCOBuVX0jzjFzgbkAhxxySH8vZXLEQHPo4QBcUeJh4YwqyrwF7NrXxuhhHkYUD7yee6IbUJHbaatWTc7ptQevqh/29KeXz/pV9VhgLHC8iBwT55jlqjpFVadUVFT0/5uYnDDQ2Sm+Tj8VJR6u/tpEFj1by/nLX2fh0xvZ2dA2KL34RLNsWnx+W7Vqck7aygWLyI1As6renugYG2QdGgYyi6ausZ2NH+1j4dMbuz0FDMauTJaDN7kmI+WCRaQC6FDVBhHxAl8Bbk3V9UzuSLa0bzwji90cVlGcsl2ZwtMjn1wwjbaOAE7BZtGYnJXKUgVjgAdDeXgHsFJVn03h9cwQ4YAB5fF7Pb9DqCy1tXwm9/Ua4EWkke6zX/YBa4Efq+r78T6nqhuAyQNuoTFR9jb7uOm5Wm6dWc11T2yIpFGWzamxVaZZyHbByqxkevD/BXwMPAoI8G3gIGAzcB9wUqoaZ0xXvk4/L9bupq7RF5lF09DawSgLHFnHdsHKvGTmwZ+uqstUtVFV96vqcuDrqvo4UJ7i9hkTIzzLZf32BuatWMf5y19n0bO1OBxWNy/b2C5YmZfMb0VARM4TEUfoz3lR79nCJZNy0fXfnQ6sCFiOGOiiNjNwyaRoLgJ+BiwJ/fwaMDs0M+bKVDXMGIj/mP/QZcfz5IJpdHQGLK+bxQa6qM0MXDILnd5X1W+o6qjQn2+o6t9VtVVV/5yORpqhK95j/sX3vYkgWVGXfSC7S+WDnr6/lVzOvGRm0YwlWCp4OsGUzJ+BH6jqjhS3zZisfswf6oOIvX1/K7mcecnk4O8HngE+AxwM/Db0mjEpl+wGHZkw1AcRk/n+4UVt2fC0NRQlE+ArVPV+Ve0M/XkAsKIxJi2y+TE/m58u0iEbv/9QT5l1lcwg6x4RmQ08Fvr5AmBv6ppkzAHZ/Jg/1AcRs+37D/WUWTzJ9OAvA84DdgE7gVmh14xJi2x9zM/mp4t0yLbvP9RTZvH02oNX1X8AZ6WhLcakTCqWzGfz00U6ZNv3z8aUUaYlDPAi8j/0sJBJVa9KSYuMGWSpfHQfSGXMfND1+4dz4JkI+NmWMsoGPfXgrTC7yQsD3QfWJCedOfB4T2ThlFHX6w+VlFk8CQO8qj7Y9TUROUhVd6W2ScYMLnt0T4903Uh7upFkU8ooG/S1QtPvUtIKkxI2ZSwom+fS55N03Uh7GkzN1gH5TOlrgB/a/1o5JNzLOWfJGqbf+jLnLFnD5k8ah2SQz7bZHvkqXTdSeyJLXl8D/D0paYUZdDZl7IDo2R5rrjuZpxZMH9Jzowdb+EkxEAiwbE5Nym+k9kSWvD5t2aeqS3o/ymSDfOzlDGSqY6LZLrbj0MB0zYefVlXJo1ecgNMhKfv3tMHU5KVyT1aTQfk2ZSwVMzRs5WPPkrn5dX1SfLF2N7U7G1M6Qynb5t9nM9sGJ0/lW945FSknS2MlluwYTqaeFG0wNTnWg89T+dbLSUUgycc01mBJdspjvj0p5hvrweexfOrlpGJgzQbrEkv25jeYT4o2rXfwWQ/e5IS+Dqwlkz+2wbrEku2ZD9aToo2HpIaoZs9dcsqUKbp2rVVIMPElO+OlL8HCZtHE19kZ4N1PGpn/8LrIv+HS2TUcNboUl2vwH/zrGts5Z8mabjcUKyfROxFZp6pT4r1nPXiTM5It7NWXJfNDvVhYIvWtHdz10nssnFFFmbeAhtDPN59Tzchi96DfFG08JDUswJu8EwgEIoEpoIo/oBQ4Hfg6/QQCOqBplQ2tPlp9fvyqFBY4GVWc22Mbifg6/bxYu5sXa3fHvH7jN/wpSaXYYG1qpGyQVUTGicjLIrJJRN4RkR+k6lrGhAUCyp5mH4uereWW598F4Pon/8b5y1/n/OWv97tcQyCgfLC3mc27Gjl/+et86bbVfGvJq3lb/iHRALSIpGRqaa5O6832geGU5eBFZAwwRlXfEpFSYB3wTVWtTfQZy8EPLanIf0fncpfNqWHRs7WDkteta2xn40f7WPj0xiGRJ040jjGiqIAT/u//63b8mutO5uDyogFfM5fGQ7JlYDgjOXhV3Ulwiz9UtVFENgEHAwkDvBk6UvXLEZ3LLfMW9DmvmyjI+Dr9FLmdQyZPnGh2zN5mX8pSKbk2HpIL+wykJQcvIuOBycAb6bieyX6p+uWIzuU2tHZwWlUlM2vGRQYKn1i3PWEw6umm43Y5afH5GVvupaLEw/yTDqfMW0CLz4/XnZ954ngB16aWHpALA8MpD/AiUgI8AfxQVffHeX8uMBfgkEMOSXVzTJZI9MsRCAQGtOVbdAB6qfYTrjxlAgseeStmql+5tyDugCkQCVyTx5Ux/6TDaW7vZOe+VipLPBw6soi7L5yMrzPAj1a+fSDAzZlCmTe70wmDJd9WSA9ELgwMp3QevIgUAM8CL6jqf/V2vOXgh454855Pq6rkB185knkr1vU7bRMduH3+AHPufbPbL+CTC6bR1NbJJ/vbuGbVhsi1Hr78BE66fTWTx5Vx9dcmct0TG2JuDBMrS9jT4uPcpa/FPWdlaeHg/QOZiGzNzedCDj6Vs2gEuBfYlExwN0NLvFkTPz2zKhLcoe8zNMK/cGf9PHjj2L2/Pe5TQltHgA/3tkSCe/j1bXuaGVvuZf5Jh0eCe/i9+Q+vY0+zjw5/IOE5zeDL5o1rcmGfgVTWopkOzAFOEZG/hv58PYXXMzkk+pfjjZ+cwsp5J6LAwhlVTB5XFjmuLznN6Lx+Q2tHJGcebWy5F6cQd8D0rpe2sGx2DSOL3XGDeIc/gFMk4TnN4Mv2ip/ZXu8pZQFeVf+sqqKq1ap6bOiP7elqIhwOYWSxm09bOjhv2Wt8efFqFj1by9VfmxgJ8n3JaUbn9Zeu3kp5cQGLZ1XHzq2eMwWv2xk3+Nc1tTOmrJAxwwvjBnGHQ3CIcPeFk2POuXhWdd4OtGZaLgxkZjNbyWoyKl4P7bonNrBwRhWLnq3t0wwNt8sZM2umqa2T8iI3j313Ku2dfvY0+Rg93EOZ182hI4tYPKs6Jgd/z8XBwVKlvdt7i2dV838eXU9dUztLZ9dw53nHElClxedn9LBCyrxDbxZJX/Q3j55oIFNEBrQqeaiwYmMmLRL9gn9U38L0W1/udvwr15yE1+3q04BavAJZSy46jodf+5CV63YEz3vtyTgFClwO3E6hud2PX6GwwBEpOxBetfrh3hbGjfCyp8mHquIQiUy1vOmbn6fDH8j7kgWDYSCDkfE+e+vMah58dRs/+urErMt5Z4IVGzMZ1dv88ng9NK/b1et8+K43DUUjwR2CTwMLHnmLhTOqWLluB2PLvWzd3cSlD/yFseVels2pYVSxG7fDEXMjcTiE8SOLKS0soMPvR4Cro3rzv7joOOqa2gc022coGciah/BYzcp5J/JxQyt7m33c/sJm1m9vSPnWgPnANvwwKdfTQFl/a5DEm13R0h4/X1vmLYikWe56aQuTx5WxcEYVrT4/+9s6ueGpDd1mZoQHzwThx79+O6bte5p8A5rtM9QMNI/ucAiqyqylrzFvxTrWb2/o8zmGKuvBm5Tr6Re8vwtn4t00wtMcuz4NjC338vjcqVz56HqAbnPcb51ZzZ1/2MzN51R36w12BLpPixxKJQsGw2AsCMqFRUXZyHrwpl/6UkWvt63x+jPVLN5NIzzNsevTwJjhwWvVNbXHneN+3RMbmFkzLm6AjjctMtH0y4IUbISRaYNRLXEwKkXmarXJTLNBVtNnfR00S8WKv0Q7AD1z5XT8Abo9DYTb0Nzeyaylr3U736r5J3LoyOJum1k4HMp7u5piZtTcfeFkXE5HTA5+8axqRg8rZPzIYoCsXHnZV4P5/9tgrEbN1hWtmdbTIKsFeNNn/dlebbB/OfsTfAIBZdf+Ns5b1r3UwKNXnMBnhnvZUtcUc86HLjsegA/3tlAUmj9/6MgiRhQXsP4f+yhyO2lo7WDp6q3UNbXzzJXT+WR/e8aXrw8G20YvN9gsGjOo+jNoFl2ZcDCCfTK5+3jXqSzxsHR2Tbe9RscMK6S+taNbXv/i+97kmSunU1o4HF+nHxHBKdDc7ueul7ZEBvzCWn3+rC8hmyxbZJT7LMCbPhvIgNdgPvb3VD880XVGD/PE3Wv0pnM+T1tH/IDWGioJ3NDSwbyHY9Myt/0+OGVv8rgyrjp1Ap0BzZugmMqBTUu3pEf+jQqZlBvIgNdAa4skO+gX7zp3/mEzrT4/l3/xswDc8vy7zFuxjhdrd9PS7mfr7ua4g6f+gPL29n2R4B4+3zWrNnDVqROYPK6Ma0+fyMKnN/LursYeB5RzSaoGNrO5gFi+sR686bOB1AQfyGN/X3r/Xa8zeVwZ35l2GOcvfz1meuTtL2ymrqmdbXuaueulLdw6szpmCuWyOTXc9Fwtl3/xs3HbfXhlCT+/cHLkvEtXb+12jlyd7ZGq2u/p2AnJnhCCLMCbfunv9moDeexvaPWxa18bd5w7KTKw2TUwhH+x/arcf8kXInnyRNMjF519DJWlHn76m42s397A7S9sjqRvxpV7cTqEy7/4WUYUuzmtqpIXa3fHtNtb4Iy5mUSf46iDSilwOqgsyZ4yBn0NfKnYRi/Vuf1sqdOeDSxFY9JqICtXdza0sfDpjZy//PVI1cmKEk8kMEQ/+n/pttUsfHoj154erEyZqATw4ZUljCkrpK6pHQgG6Hkr1nHvn9+nvrWDmUtf4/zlr3PpA3/hylMmcFpVZbd2d53nv357A4uereXdXY2ct+w1ttQ1ZUX6IVtSI72tixiobC8xnE42TdKkXX8enxNN2Vt09jEcc/BwKko9CY95fO5URCTu9MinFkxnZLG7W4/v0StO4MJfvtHt+AcuPR6XQyjyOGOKk8UriBWumZItUwuzZdpjqnvYiQrYrbnuZA4uLxrw+bONTZM0WaU/j/2JHusPG1Uc6f0nOgbgoGGFcTeLdjiUHQ0tDCt08eT3pgWrRjocCc/ldjkYW+aNCUTRuepWXyebdjVGgnv4c9kwiyZd0x57u4Gnel9XK2twgAV4kxMS/dIWeZyRwNDTL3bXoFLgctDU1snZP381zmrUQvY2++JXuSxwxg1E4ZtWXSMserY2K4NLOgJfsr3zVOT2w6I3Xs/1ge6Bshy8yQmJcvejij29HhP+xQ7vIOV2OenoDHDxfW/G5GmvWbWBD/e2sLfZR7m3gKVd6tosnV1DubeAQED5tDk4VfMfnzazu7EtksfO5pop6WhbNuS/c2Gv1HTJrhz8kaW6dklNppthMkw1WMVRFUSgwOFAJPHrXT/r8wdQVUQEt/PAMarQ0uHnvV2NHF5RTO3O/d2uXTVmGJ4CJyKwbU8zFSUeXA6hM6DUNbVz2Khi/AHF1xlg1/42Kko8FDgdFDgFj8uZdDszJdVta+8MsP4f9d1en3xIOZ48LMaWDeSrf7IcvMl+4eATCCitHQE+qm+hw68ceVApRaGg63YmDhLRAby904/H5eSoMaU4HUK4H7OjvoX2Tj+dAcXjctIelX/2uJz4Q4FPFeqbfdR36XkeOqKYto4AuxvbGDOskK11zZFrHTm6hCK3K9LO8Pfx+QORYAqZDf69/RsOxvnj/btmyw1uqMmuAD9sInxldaZbYTKgp5kodX9rT2qmx54us0Qmjyvj2iMnxlSCDJ+THXD9GUdxdWgzj3AOfm+pB5/HhV/h6re7z7p5/Nyp7KhvpcMf4PIn/xb7fp2XJ783LbRfaIA9zb5uuz55XI5Iaiif5mdHD6yODCg3PVfLi7W7I9/RNboUcvw7Zq/E/67ZlaKxaZJDVqIpfAtnVDFvxbqkprh1nR63bE5N3AHPFZcdz3u7m3ip9hO+VTOWzwwvxOkQHAI3PvMOL9bu5rSqSq469ciYomThWjY76lspL3LzcUNrZMFVeMbMy1efxJx734hsGh5vWuelD/wl5rVsmEI5EPFuzuHtEB1dtkM0g8+mSZqsl2gKX3i7vWRmenSdJVLmLYh7zt2N7Sx6tpbFs6oZVeLG43LgV7jxmY2Rlap1jT5afX5+NXcqgahNuf3+AA4RLrgnfsmDD/Y0R9od79pFbme317JhCuVAxBtYnbdiXc7fuPKBjXqYrJBodWOLz5/0TI+us0QS7bzU4Q+wcEYVBU4HDhHmP/wW5y17je9MO4zJ48qYPK6Mq782kXv+dyvbP23BHwjQ6vNT39LO7qZ22jr8LJxRxeRxZZGSB1edOoFfXHQcd720BYAOfyDh9+n6WianUA7Gjk1WVjh7WQ/eZIV4c5eXza5hTFkhZd7kHvG7znX3up3dznnHuZPwFDi4PpQ/j+6BX/fEBhbOqALgwVe3seDkI2j1+Zlz75tUlHi49vT4+fz12xs4dGQRHpeDG848mrYOP6NK3SyeVR1z/PI5NcGFUqGnjExPoextznqyK45tYVH2shy8yRqpqAAYPmdrh5+tu5soLHBEgm7Y2HIv91/yBT5t9lE5zEOnPwAIO+pbWfj0RnbUtybM54dz7Y9ccQIXhUobhAdsn3rrI06tGk2Zt4AWn58Sj4v//N0mrjp1AoeNKo4pd5AJPZUuiFe+IdGAsBX3yizLwZuckIrVjeFzdnYGaGzroNUXP52wr7WD85e/zmlVlfyfU4+kvaOTIrez13z+yGI3S2fXcPNztTE56GtWbYgMECziTNUAABvtSURBVIc9Pncq67c3cOkDf4nUyNnb7MvYIGRPqZW+lPQdjNID0Td3r9tJZ0Dp6AwM6VK/gyFlOXgRuU9EdovIxlRdwwwdA8kVBwLKlromfvbH9xgeGrSNNrbcG1lpObNmHN97eB17m30xOfyG1o64nxvuLUCEmDLCcGCAOPrYhtaOmPd31LdmdLOLnqo69jWvHr6RHlxeREVp355KoqtcXvnoejbvauRbS161zUAGQSoHWR8ATk/h+c0QMdAyt+He6Iu1u7l21QZunVkds1z/1pnVLF29FTjQU1+6eiujSoJ59LHlXpau3hr5e/Tnro3KsUeLHlANp2zC1wi/1tDakdFStj2VLkh1Sd9o0U8L8086PCaFNpRL/Q6GlKVoVPUVERmfqvOboWNPc3vcdMGTC6ZRWVoY9zPRj/wAFSXB+evRG3JMqCzB6RB27mvj+jOOoqG1g4BqJLD5VXE5HKy47HgCCkUeJyvnTuXjfW3sbfZFBlifWLe920bed5w7iYAqj8+dSoc/gNftjNScj1lwReZmnPSUWklnwa7op4VEqbB4/z62a1PvMp6DF5G5wFyAQw45JMOtMdko0WbYbR2BuMfHG/SL3iA7vCHH4lnVOERiVrPeed4k7rtkCh/Vt3Hlo+u7DUAunlUNHKgYObbcy6XTD+OhVz9g0dnHMH5UMYUuR8yceoDTqir51dypqMLfdzfFlBMezJ7xYO3YlOqSvtGiZ+GEU2G9zcixgd3kZHwevKouV9UpqjqloqIi080xWcgpEjdd4Iz6PY7O0e/a39atxx/eIDv82WVzahARfhwK7uHjfrTybfa3djJuhDfuTcUhwm2/38yis4/hlWtOYtHZx3Db7zezct0O7nppCx/sacbnD3DDmVUxuz99Z9ph/Mdv38FT4OCg4YUxvfnB6hkPJJUVb4xjIHn1vohOFcVLhcX798mGqpW5IOM9eGN643U7u80pXzyrGm9oVWjX3tyq+SfGDc6HV5aw5rqTcbuclHsL2NHQGvc4f0DZ/mlrt57kaVWVjCh286/fqKLE48LllEjZgfDiqOjNtpfOrmHhjCraOgLc88r7vFi7mxu/EUhZz7i/m1lnujccb/3Ckwum9TiLxhZXJccCvMl6ZV43o4cVsujsYyhyO2nx+Rk9LLgACroHtp4264gOdEWe+At0Wnx+xo0Ipmt+tDLYwz+tqpIrT5nApQ/8JSaAn1ZVSV2jj9tmVUfeg2Cwmf/wusg8+VtnVtPQ6otsPpKKJfzJBL14KZxEYxx9LTUwkJx4X/9NbHFVclIW4EXkMeAkYJSI7ABuVNV7U3U9k78cDmH8yGJKCwviBo+ugW3p6q3cOrM6pjcd7zF/VLEn7urZilIPP/3N36hr9LFwRhWVpR5GlXgi9WfgQAB/+PIT2LmvFV+o/EGZtyCmANmEyhIWzqjiwVe38dMzq1K6arW3oBevp/7QZcfjEBlwbzjRU8CEihLqWzsG/WnFdm1KTipn0VyQqnOboSdeDy/cYwRiAtv67Q08+Oq2yOrUseVexgz3dgss8bbxczmExrZOZtaMY+nqrZGFSk8tmBY/naPK/Wu28f2TJ8QMvN46s5oHX92GP6CMLHbzkzOOpsgdf7u/wdJb0IuXwvlwb0u3f7/wzwV92KAj3rnv/MNmfvCVI7uVTB6M1E86B4FzmaVoTE6K7jFWlHi65ei/M+0wrl21gbqmYC35RL/44RtHT/Xo129voKEl/uwOB8HFUd9/9K2Y4HbdExt45IoTuDmqLvqy2TWMSGFpgnDQe3LBNNo6AjglOH4RvhG2+DpZOKMqprxxkdvJLc+/yx3nTooMOIfHOJraOhlVrEm1NxDo/gQzs2ZcJLiH/136k/rp6ftatcqeWYA3OSm6x7ijvjUys+WzFcW8X9ccKd8b7sGGg1wgEMCvoKoxvb54PdAHX93Gzy6YTKc/gMfl4KFLj+fi+9+MuQHs2t/GyGJ33N59Q0tHZKrkjvpW5j2cnhK6e5sOfJd4de2jb1wtPj91Te0EVCNjHA2tHdz2+82Rm2P4Bpgovx4IKHuafd2eYEo8LhsIzTAL8CYndc27h2u8rLnuZI45eDg/v3ByJBABbP6kkTv/sJnvTDusW25+4ujSbuebPK6M70w7jAuj6r4vnV3D01dOo9UXwB9Qbn6ulrpGH3ecNylu737X/raYNqcjuHW9Uc2sGRcJ7uE2XPfEBu6/5AssfuFdDh1ZxD0XT6G5vZML7nmj2/l8nf6EG3pMrCzF5XKwN2rnquhr/GruVBsIzbCMz4M3pj96Wkrfde52OOjNrBkXCe4QO3dausy1n3/S4d2Onf/wOvx+cIjQ0NLBNV87iv84+3MUuR0suei4mLnbS2fX8MS67XHbB4NThz2erjeqRCtD97V28IOvHMkh5UVMHF3KZ8q8Cf89E23o8fG+VgIBTTh7p8AhCUshmPSwAG9yUrw6Kstm1+B00C1YhgNQZaknYcrAKcTUqImXdqko8bC7sZ3zlr3G2XevYfEL71LicdHsC66offjyE1g1/0QevOx41m7bw/dPnhC3fZ2dgQHV1ulJ1xtfoo1H2jr8zFuxjvrWDhwO4aBhhQmDcaIAvruxnb3NvoQ3W4fDERkIXXPdyTy1YLqtNE0zS9GYnBQ9oNjS7mfbnmZ++puNVJS6+emZVTgdEknRhANQiccVN2UgIjgcDh58dVtkoDBcdTL62KtOnRBJd4RTOHOiNtAOl0OA4Ibed7+8hYUzglMjR5a4+cXLW3n1/b0sm13Dz156j4oST+R6u/a1MXqYhxHFA8vPd51J43HFXyQGsSmjnmaleN1O7r/kC5H8/NLVW6lrCgb3McMLGTPcm3D2jg2EZpZt+GFyWvSmFfFWk4bnYm+pa8IhsKfJF/P+rTOrOaKimIrSQjbt2h/JJYcXNi145K3IsSsuP56Tb/8TkHhD73AN+Mnjyph/0uFMPKiUbXXN3PXSlpjaM4tnVRNQYtqybHYNR48Z1q8ebvQgaHi6Z6vPj1+VHzz2V+afdHjMDJfrzziKH//67cjmHvE+G17xu6WuqVtdH6/byZKX/87N51T3OghrUss2/DB5Kzp9EC9vHp6WN6GihI/2tcb00htaO3jw1W3cfE41Docwqtgd897Dr33IwhlVHHVQKQVOB4QqTe6oT7ypdrgGfLig2WPfnRopZxBWUeJhbHlRt4VT/Z1lk2gBU0mhC3+Hn6tOndDtBuMQ4dErTqC908+O+hZuiprOGX4SqWtqZ9mcGn72x/di2nnNqg3cfu4kfvTViZF8uvXUs5MFeJPToldv9lRqtr4Vbn6uNmYWzWlVldxwZhW+Tj91je0UuBzdeuWvvr+XlfNORFUpcDoi6Y5EVQ816u9LLjqOXfvaYo6bPK6Ma0+fyCf72wZtCmHXQdBpnx1JU3snF8dJH9U1tfPApV+gvTPAhVFbDN46s5q6Rh/rtzfE7EY1b0Ww3EJ0Zcwd9a3BQdmy7ovHTHaxAG9yWnTOuadSs75OPy/W7o6UH/jM8EICSsw+qvdcPIWHLjs+JjAuueg4/i1U+ve0qkp+OqOKBy49Ho9LWDa7hnlR88sXz6qmvKiA3145ncICJ4tfeJe6Rl9M2YSrTp0QCaCDNYUw+ikmnBq65flNMU8j96/Zxn+dP4n3Pmmiqa2TKx9bH9MrD284Hk5RhZ9EdtS3dpv1Eq7rY8E9+1mANzktenAwEAiwbE5Nt6Xx4Rzz2HIv67c3MG/Fum459HA659fzTows+BlR7GbxC+9Geq91jT4+2dfGj1a+TUWJh7suODbu4qBwYA1/7vYXNrN4VjUHDSvE4RAWzqhiX3N7t01ComcB9SV4Rj/FXHv6RFxO6Tbf/9aZwZWp81as44///KUe00vh3abCf68s9UTOb1Mdc4sFeJPzovO/FaWFSe1QlGj1aYc/EMmZPz53aiRITx5Xxn9/+1jer2vmjnMnMaLYzf62zm75dQjOPQ/v5xq+RkCJzLiZ90/jmXHsWO566b3ILJuKUg/PrP+Ix9ft6HO9lvB3u/MPmynxuPB1arexiPDiJjhQX7/r00P4CSiczgkvaPK4HKycdyJOAYfDYQOoOcQCvMkrye5QBPELbDkdB4JfOOBVlHi48awq9rd2sPDpjZGe7N0XHsdpVZUx+emx5V4qh3nYta8tUt+l6+DvrCmHREoLhz87ttzL/Zd8gTv+uKVf9VpGlri54cwq3q9rZlRp/OqQTe2dAOza39Zt6uTS2TWMLC7gye9Nw+UUfn7hZPwBjRl8tR2Tco8tdDJDRvQORS6HxN182xW1+jK8u9BVp06gvrmD7z0SW1Ds+4++xU/OODp2MdOcGgoLHBQ4HQRUuf3cSUw8qDQm4Dod8QOwMxQ4+zLYGp5B860lr/LlxatZ+PRGit3xFx6VFhYweVwZ96/ZRkWph9vPncRLP/4yi84+hoW/2ci5y15nb7OPYZ4CRIQLf/lGTC0d2zEp91gP3gxJXRc2RU+ZjM7piwitHX527Ys/6wUhZiMSVGlq66TDH+D+Ndu49vSjKHAIq+afyN5mH0tXb8Uf0LhPDwVOB8vmBEscJDvY2tDqizwthOe43/L8Ju6+8LhIhcvwzeu232/ipm8ew5iyQoZ5Cihyuzhv2Wsx7QiX+G312Y5J+cACvBmSRha7+dFXJyZcfVnuLeDjfcHl+MPj5NQhGJQ/2NMSk4cfW+5l0dnHALDg5COob/Zxycq3YwLtmi27+cXsGr4XNcD6i4uO49PmdkYWB1Mtwz29B/hAQNnZ0BaTNgpXinQIPHDp8TS3d1LkdtLW4WdmzTh+9tJ73HxONS6XA6eQsMTvYM7yMZljK1nNkJVo9WUgoN1WtV79tYnsbfLF5K1/cdFx/OvT70QWEIU9PncqAC0+fyT4ho0t9/KruVPxdQRwuRz4A4rLITT7OrjiwdgZNb2tao1exRt9/tvPncQh5V46AopDhF372vjP322irqk94crd6BK/Z9+9Ju6q4GVzajj6oP6ttDWpYytZjYkj0YBs1/K34Tz0DWdW8au5UwkoOAS2f9pCXVN7zGfHlgf3dPX5AwmLm+3a18aspa8FbxKzawgEAnz/0dh56cmsak1UBGz0sELOW36gzPHiWdXceFYV//5MLdc9sYGV805MWOL3se9OjUwnvf2FzZFZPsO9BQwrdFlwzzE2yGpMF/ECZ3CRVDsel5NDRhQhwP1rtnHHuZNiBlkXz6qmvLiAl2o/iRQsiza23BsZqNxR38r3Hl7HiOL4N4K2js4eK0y6XU5Oq6pk2ZwaHp87NbgG4J/G4w8od5w7iWVzaqgo8XDNqg3UN3cw/6TD2VHfimriEr/N7R0snlUdCfKLnq2lvTPA4hfexeGwcJFrrAdvTBeJNq+uLPVEFvh43U4unX4Y9/75fW751uf5TFlwkHTXvjb+/Zlarjp1Arc8v6nb5t/hHHlYeO59vOu1dwZTRRMrSyN17X2dfkQEp4CnwBGzW9NpVZVcc/pRbP+0hSK3E7fTwY1nVfHo6/9g/MgiFLj/ki/gdTvxB+JPE925r51DRxbx2Hen8sn+NvY2+3jw1W0xdWdM7rAcvDFd9LaDUfiYD/Y28+HeYDBVYEJlMW0dATpDve4vL17N5HFl3Darmn2tHQz3FsSsjIVgUH30uyfQ6uvE6XDikOCiKLdLUFUu+uWb/HreVBpaO2O24bv+jKNxOSRSTwbgse+eANCtNPDYci8X3BNbkiFcYTPmO86uYUxZIWXeA5t0W3XI7Gc5eGP6oKfa6NHHjB9ZTGlhQeSY8qIDx+ze3xZJc1y7agNXf20ii194l8u/+FlqdzbGBGBVxelwRHreLT4/Y0d4GVEcLJ7W1hmIbC5+94WTcTocXHzfm9xx7qSYmvKfKfN2q1B5zaoNrLj8+LgVNnv7jlYdMvdZgDcmjmTK3yY6JhBQmto7ufO8Sfxo5dus397Ag69u4/ozjsbjcnSrX/MvXz+azkAgZrrj4lnVlHpcnFZViT+gVJR4+L8zj6GixEOLL8DDV5xAgUO4ZeYxPLluB7OmHEJANW5e3d8ljx+ez24lfvOfBXhjBtneZh8X3/dmpHd9RGUJ/9jbwo9Xvs31ZxzVrX7NqBI3c6KOD9eyAbjxG58joMpN3zyGkcUePt7XHjt/fnYNF049lIt++Sa3fOvzcfPqe5piV5/afPahw4bFjRlk4Rkq4cqVV698myK3k7qm9kh9m2h+DW4CcvXXJrLo2VrOX/46C5/eyJ4mH//+23fwdQZAoLGtMxLc4cAsnIAG/37Hi+91m9Xzs28fy8jigpjXrBrk0GE9eGMGWddZOOu3N3Dvn9/ngUuPp8jt6FZH3uNycNWpE7pVgJz/8Doeuux49rd2sGT13/mXM6vil0sgWO1y/fYGbnn+3cguVIGAsqfJh69TI5uW2IDp0JLSHryInC4im0Xk7yJyfSqvZUy2CJfvje41/+DUIxlRXMBBw7xMHF3KI1ecwKr5J7JwRhWKMn5UUcICZKNKPfzL16siZX6jhcslzD/pcODAVoG+zgAn3/Enzl32GsUeFwcNK+Tg8iIqSj0W3IeQlPXgRcQJ3A18FdgB/EVEnlHV2lRd05hs0NssHIdDOHh4cF/UEcVuHCLgiD8v3R/QyK5T8/5pfLdNQsLz6m848+jIZ35x0XHc88r7kZ+LPLb70lCVyhTN8cDfVfV9ABH5FXA2YAHe5L3eZqi4XA7GDCtk5/42vG4H7R2BbjXal1x0HLc8vykS9Jf97wecOengmFk4t78Q3EVqzPBC/nTNSTgdwopXt7Fy3Y5Ivn1Usc2UGapSGeAPBrZH/bwDOKHrQSIyF5gLcMghh6SwOcZkF5fLwcFlXhpafYgE552Hg3flMA9NbZ0xi6IAbnz6Hf7j7M9FatOPLfeyfE4NLocE5+J7C7jiS0dw8bTDLN9uUhrg4/1X1W3ZrKouB5ZDcCVrCttjTNZxOIQRxR4CAaXe2U6R20l7p+IPDZB2TdvUNbVTUerh8blTCajiKXAyqthjC5RMXKkcZN0BjIv6eSzwcQqvZ0zOcjiEkSWFVJYUUux2MtzrYtyIA6UG4MCOUS6nMGa4l3EjiqksLbQeukkolT34vwATROQw4CPg28CFKbyeMTnP5XJQOawQgNICF6UeV7BEcUBxuxyMLHJTUGCLlExyUhbgVbVTRK4EXgCcwH2q+k6qrmdMviksdFFYaEtVTP+l9L8eVf0d8LtUXsMYY0x8VqrAGGPylAV4Y4zJUxbgjTEmT1mAN8aYPJVVW/aJSB3wYT8+OgrYM8jNSQdrd3pZu9PL2p0eh6pqRbw3sirA95eIrE20J2E2s3anl7U7vazdmWcpGmOMyVMW4I0xJk/lS4BfnukG9JO1O72s3ell7c6wvMjBG2OM6S5fevDGGGO6sABvjDF5Ku8CvIhcLSIqIqMy3ZZkiMgiEdkgIn8VkRdF5DOZblMyRGSxiLwbavtTIlKW6TYlQ0TOFZF3RCQgIlk/FS4XN64XkftEZLeIbMx0W/pCRMaJyMsisin038gPMt2mgcqrAC8i4whu8v2PTLelDxararWqHgs8C/xrphuUpD8Ax6hqNfAe8JMMtydZG4FvAa9kuiG9idq4/gygCrhARKoy26qkPACcnulG9EMn8GNVPRqYCnw/R/69E8qrAA/cCVxLnK0Bs5Wq7o/6sZgcabuqvqiqnaEfXye4Y1fWU9VNqro50+1IUmTjelX1AeGN67Oaqr4CfJrpdvSVqu5U1bdCf28ENhHcWzpn5c1uAiJyFvCRqr4tkltbmInIzcDFwD7g5Aw3pz8uAx7PdCPyUFIb15vBJyLjgcnAG5ltycDkVIAXkT8CB8V56wbgX4DT0tui5PTUblV9WlVvAG4QkZ8AVwI3prWBCfTW7tAxNxB8tH0knW3rSTLtzhFJbVxvBpeIlABPAD/s8oSdc3IqwKvqV+K9LiKfBw4Dwr33scBbInK8qu5KYxPjStTuOB4FniNLAnxv7RaR7wAzgFM1ixZU9OHfO9vZxvVpJiIFBIP7I6r6ZKbbM1A5FeATUdW/AZXhn0XkA2CKqmZ9RTgRmaCqW0I/ngW8m8n2JEtETgeuA76sqi2Zbk+eso3r00iCvcN7gU2q+l+Zbs9gyLdB1lx0i4hsFJENBFNMuTI16+dAKfCH0BTPpZluUDJE5BwR2QGcCDwnIi9kuk2JhAaxwxvXbwJW5sLG9SLyGPAaMFFEdojI5ZluU5KmA3OAU0L/Tf9VRL6e6UYNhJUqMMaYPGU9eGOMyVMW4I0xJk9ZgDfGmDxlAd4YY/KUBXhjjMlTFuBN1hGRpgF89pFQ9cWNoaqGBVHvfVNE/jX0938TkY9CU+HeFZFfiIgj9N7tInJKgvOP72+VRBEpEpHnQtd7R0Ru6fL+D0XkYhG5O9SuWhFpjZqyN6unthnTlQV4k28eAY4CPg94gSui3rsWWBL1852hKp5VoeO/HHr9f4BUlea9XVWPIljnZLqInAEgIi6CNX0eVdXvh9r1dWCrqh4b+rMqxW0zecYCvMlaErQ41Bv/m4icH3rdISJLQr3gZ0XkdyIyC0BVf6chwJuEqlyKyJFAe4LVzW6gEKgPneNDYKSIxKtnE92+QhG5P9S29SJycuj1IhFZGaqV/7iIvCEiU1S1RVVfDl3DB7zFgSqcpwBvRVXojCvZthkDFuBNdvsWcCwwCfgKsFhExoReH0+w130FwVWpMUKpmTnA70MvTScYUKP9SET+CuwE3lPVv0a991boMz35PoCqfh64AHhQRAqBBUB9qFb+IqAmTvvKgG8AL0W1b10v1+tL24yxAG+y2heBx1TVr6qfAH8CvhB6/deqGggVk3s5zmeXAK+o6v+Gfh4D1HU5JpyiqQSKReTbUe/tBnrbXeuLwAoAVX0X+BA4MvT6r0KvbwQ2RH8olI55DLhLVd/voX2JJNM2YyzAm6yWqLB/jwX/ReRGoAL456iXWwmmYbpR1Q6CPf0vRb1cCLSKyAlRg5xnDUb7gOXAFlX972TaF0dh6HhjemQB3mSzV4DzRcQpIhUEA/CbwJ+BmaFc/GjgpPAHROQK4GvABaoaiDrXJuCIeBcJVRGcBmyNevlIYKOqvhE1yPlMnPZdFDrHkcAhwOZQ+84LvR4ewA1f6yZgOPDDLudK2L44jiS49aAxPbIAb7LZUwTTG28D/w+4NpSSeYJgrfSNwDKCu+7sC31mKTAaeC3U6w7vcfsKMFlit/sK5+A3EiydvQQi+fsjgLW9tG8J4BSRvxHc0eoSVW0PvV4RqhB6Xeg77BORsQQ3p6kiuF/BX0M3JIDniX2CiKsPbTPGqkma3CQiJaraJCIjCfbqp/e2uYuI/Az4rar+sZfjzgGOU9WF/WybEyhQ1TYROZzgQOqRoZkzPX3uKYI3sS09HDOgtpmhJS82/DBD0rOhmShuYFGSO3f9J8ntaeoC7hhA24qAl0O9bQG+11twD7me4GBrwgA/CG0zQ4j14I0xJk9ZDt4YY/KUBXhjjMlTFuCNMSZPWYA3xpg8ZQHeGGPy1P8HFlbRqLRs73YAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "#graph volcano\n",
    "import matplotlib.pyplot as plt\n",
    "import seaborn as sns\n",
    "from numpy import log10\n",
    "tstats['-log p'] = -log10(tstats['pvalue'])\n",
    "\n",
    "figure = sns.scatterplot(x='log2(B)-log2(T)',y='-log p',data=tstats)\n",
    "figure.set_title(\"B cells vs T cells\")\n",
    "#plt.axvline(1, color ='orange')\n",
    "#plt.axvline(-1, color ='orange')\n",
    "plt.axhline(-log10(.05), color ='orange')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Next we will identify the altered proteins."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "altered = tstats.apply(data_utils.is_altered, axis=1, change_factor=False).dropna()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "These are significantly altered (p<0.01) without restricting by the magnitude of change:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "T cells    30\n",
       "B cells    15\n",
       "dtype: int64"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "altered.value_counts()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Here are counts for other p and change cutoffs:\n",
      "P<0.05, >150% of the other type\n",
      "T cells    65\n",
      "B cells    22\n",
      "dtype: int64\n",
      "\n",
      "P<0.05, >200% of the other type\n",
      "T cells    35\n",
      "B cells    14\n",
      "dtype: int64\n",
      "\n",
      "P<0.05, >500% of the other type\n",
      "B cells    2\n",
      "T cells    2\n",
      "dtype: int64\n",
      "\n",
      "P<0.01, >150% of the other type\n",
      "T cells    28\n",
      "B cells    15\n",
      "dtype: int64\n",
      "\n",
      "P<0.01, >200% of the other type\n",
      "T cells    19\n",
      "B cells     9\n",
      "dtype: int64\n",
      "\n",
      "P<0.01, >500% of the other type\n",
      "B cells    2\n",
      "T cells    2\n",
      "dtype: int64\n"
     ]
    }
   ],
   "source": [
    "print(\"Here are counts for other p and change cutoffs:\")\n",
    "print(\"P<0.05, >150% of the other type\")\n",
    "print(tstats.apply(data_utils.is_altered, axis=1, pvalue=.05,change_factor=1.5).value_counts()) #50% increase from either\n",
    "print(\"\\nP<0.05, >200% of the other type\")\n",
    "print(tstats.apply(data_utils.is_altered, axis=1, pvalue=.05,change_factor=2).value_counts()) #Double\n",
    "print(\"\\nP<0.05, >500% of the other type\")\n",
    "print(tstats.apply(data_utils.is_altered, axis=1, pvalue=.05,change_factor=5).value_counts()) #5 times larger\n",
    "\n",
    "print(\"\\nP<0.01, >150% of the other type\")\n",
    "print(tstats.apply(data_utils.is_altered, axis=1, pvalue=.01,change_factor=1.5).value_counts()) #50% increase from either\n",
    "print(\"\\nP<0.01, >200% of the other type\")\n",
    "print(tstats.apply(data_utils.is_altered, axis=1, pvalue=.01,change_factor=2).value_counts()) #Double\n",
    "print(\"\\nP<0.01, >500% of the other type\")\n",
    "print(tstats.apply(data_utils.is_altered, axis=1, pvalue=.01,change_factor=5).value_counts()) #5 times larger"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Differentially expressed proteins"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(array([-1.,  0.,  1.,  2.,  3.,  4.,  5.,  6.]),\n",
       " <a list of 8 Text major ticklabel objects>)"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXsAAAEKCAYAAADzQPVvAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3de3yU5Znw8d/1zCGZJEACJIAGKlpEU4xiqHJoFXXVWqmsBbVFsOLKoei6dVsP3Zbqu7R9tWhdqSKHd+sB0YK4rud6oPp2q3ggUBFR5LWKYAUCJCHnyczc7x9zYCYzk0wOc8pc388nH81kMs+dANfcz3Vf93WLMQallFL9m5XuASillEo+DfZKKZUDNNgrpVQO0GCvlFI5QIO9UkrlAHu6BxDL0KFDzXHHHZfuYSilVFaprq4+aIwpjfW1jAz2xx13HJs3b073MJRSKquIyO54X9M0jlJK5QAN9koplQM02CulVA7QYK+UUjlAg71SSuWAjKzGUbnF5/XibaxBvG0YWx62olIsmy3dw1KqX9Fgr9LK5/Xi278Dx/pZUPc5FI/Cc/ljMKxCA75SfUjTOCqtvI012IOBHqDuc+zrZ+FtrEnvwJTqZzTYq7QSb9vRQB9U9znidadnQEr1UxrsVVoZWx4Uj4p8sHgUxuZMz4CU6qc02Ku0shWV+nP0wYAfyNnbimK291BK9ZAu0Koe6asKGstmg2EVtM99BfG6MTanVuMolQQa7FW39XUFjWWzYQ0anoSRKqWCNI2juk0raJTKPhrsVbdpBY1S2UeDveo2raBRKvtosFfdphU0SmUfXaBV3aYVNEplHw32qke0gkap7KLBXimVU3K1y6oGe6VUzsjlLqu6QKuUyhm5vEdEg71SKmfk8h4RDfZKqZyRy3tENNgrpXJGLu8R0QVapVTOyOU9Ikmb2YvI7SJiOnzsS9b1lFIqEZbNhmPQcOyDR+EYNDwnAj0kf2a/E5ga9rk3yddTSikVQ7KDvccYo7N5pZRKs2Qv0B4vIl+IyKci8gcROT7J11NKKRVDMoP928DVwEXAPGA48KaIDIn1ZBGZLyKbRWRzTU3/3+CglFKpJMaY1FxIpAj4G3CHMea3nT13woQJZvPmzSkZl1JK9RciUm2MmRDraymrszfGNAIfAGNSdU2llFJ+KQv2IpIPnAR8maprKqWU8ktmnf1dInK2iIwWkTOBDUAh8HCyrqmUUiq2ZJZelgOPA0OBGuAtYKIxZncSr6mUUiqGpAV7Y8z3kvXaSimlukcboSmlVA7QYK+UUjlAg71SSuUAbXGsUi5XD3xWKp002KuUyuUDn5VKJ03jqJTK5QOflUonDfYqpXL5wGel0knTOCqlQgc+hwf8HDnwWamOUrl+pTN7lVK5fOCzUuFC61cPno99WSWOB8/Ht38HPm9yDvTTmb1KqVw+8FmpcN7GmqOFChBav2qf+wrWoOF9fj0N9irlLJstKX+ZlcomqV6/0jSOUkqlQWj9KlwS16802CulVBqkev1K0zhKKZUGsdavxLLhq9+LNwmVORrslVKqj3S3lDK4fhWszLEncWe5pnGUUqoP9KaUMhU7yzXYK6VUH+hNwE5FZY4Ge6WU6gO9CdipqMzRYK+UUn2gNwE7FZU5ukCrlFJ9IBiwOy6yJhKwU7GzXIO9Ukr1gd4G7GTvLNdgr5RSfSSTW4Fozl4ppXKABnullMoBGuyVUioHpCzYi8i/iYgRkftSdU2llFJ+KQn2IjIRmAdsS8X1lFJKRUp6sBeRQcBa4J+A2mRfTymlVLRUzOxXARuMMX/q7EkiMl9ENovI5pqavmv+o5RSKsnBXkTmAV8FFnf1XGPMKmPMBGPMhNJSPXxaKaX6UtI2VYnIWODXwDeNMck5VDFMd/tIK6VULknmDtpJwFBgu4gEH7MBZ4nIQqDQGNPWFxcK9ZFOYuN/pZTKZslM4/w3cApwWtjHZuAPgf/vs9l+Khr/K6VUNkvazN4YUwfUhT8mIk3AYWPM9r68Vioa/yulVDbrFztoU9H4XymlsllKg70xZqox5vq+ft1UNP5XPefzemmv34fn8G7a6/cldCanUqpv9YsWx6lo/K96RhfPlcoM/SKNA/6A7xg0HPvgUTgGDddAkiF08VypzNAvZvYqc3V38Vz3SyiVHBrsVVKFFs/DA36cxXNN+SiVPP0mjaMyU3cWz3uT8tFFYKU6pzN7lVTdWTzv6X4JvSNQqms6s1dJl+jieU/3S+gisFJd02CvMkZP90voDmqluqZpHJUxerpfIt4isM2C9vp9WtGjFDqzVxmmJ/slYt0RMH05smEujgfPx7d/R8wFW13UVblEZ/Yq64XfEdi9rcihj2Hj7bB3MwD29bNon/sK1qDhoe/RRV2Va3Rmr/qF4B2BF4G1l4UCPRAzf6+LuirXaLBX/UqiFT26qKtyjQZ71a8kWtGjbbFVrtGcvepXEq3oCb4p2Dvk7BNpi639e1Q20mCv+h3LZotYjI33nJ6UeerCrspWmsZROasnZZ66sKuylc7slepCeNrGbgkUlUUu7urCrsoCGuxVzkok9x4rbcP05RF1/Lqwq7KBpnFUVuvpLthQEH/wfOzLKuPutI2VtuHpRXD2Lf7P9bxjlSV0Zq+yVm8WS72NNUe/D0K59447bePV45shJ+K94X0971hlDZ3Zq6zVm8XSRDdVxavH99jy9bxjlVW6DPYiki8i/yoi/yUiT4rIjSKSn4rBKdWZ3uyCTXRTVU/bLiuVaRJJ4zwCNAC/C3z+fWANcFmyBqVyR282KHXnfNuOEt1U1dN6fKUyTSLBfqwx5tSwz18Tkfe6+iYRuQ5YABwXeOgD4JfGmOe7PUrVL/V2g1JvdsEmEsSDb0SWtw2bZcPYHJhe/cRKpU8iwX6riEw0xrwFICJnAm8k8H17gVuAXfjTRT8A/ltEqowx23o6YNV/JLpI2qkBwzFXPw8+Hx5bfrdm3Z3ttI1ZcnnJfdjeXonn7J/qjlmVdRJZoD0TeFNEPhORz4BNwNki8r6IxA3axpinjTEvGmP+nzHmY2PMz/Cngyb1ychV1utNzj0YjO3/eS7yH6cgj3wHaeq7XawxSy6fuR5O+77umFVZKZGZ/bd6exERseHP8RcBb8Z5znxgPsCoUaNiPUX1M73JuUfcFZRPgCk/wu5pwtfwJb4BI3o96473RoSrRHfMqqzUZbA3xuzu6YuLyCn47wTygUbgUmPM+3GuswpYBTBhwgRNjeaA3uTcQ8G4fAKc+wv/rLvuc6w+akwW742IllrdMauyUrLr7HcCpwETgQeAh0VkXJKvqbKEZbNhBRZJPTe878/VJxikQ8F4yo9CgR7os8ZkMc+1veQ++OvjWnqpspIYk7pJtIi8Cuw2xvxTZ8+bMGGC2bx5c2dPUTkulLM37bD6nKive254H/vg3qUDj1bjuBHLwogNH6KllypjiUi1MWZCrK+lul2CBeSl+JqqH7JsNnxlJ2GO7EV6mPdP5Bodq3U0xKtslcgO2gYROdLhY4+IPCUix3fyfXeIyDdF5DgROUVE/jcwFVjbh+NXOcw0HUJe+pk/vRKWbvFdsVbTLBmmpw3rVN9JZGb/W+DvwGOAAN8DhuPPx/8efwCPZTjwaOC/9cA24CJjzEu9G7JSfuJtg53PQ9N+uPDX/kqZllpMQSk2TbNkDD3dKzN0mbMXkbeNMWd2eOwtY8xEEXmvw+7aPqE5e5WI9vp9OB48PyqF0z73FRyJbspSSad/TqnTWc4+kWocn4hcLiJW4OPysK9piaRKmY6pACkcok3KskBvNs+pvpNIGudK4F5geeDzTcBsEXEB1ydrYEqFi5sKKDtJm5RluN5snlN9J6Wll4nSNI7qKNNTAb3p3tkfdPbzh8pkO7xRJ7qnQiWuV6WXIlKOv73xFPxpm78A/2KM2duno1SqE5mcCsj1Bciufn5tE50ZEsnZPwg8AxwDHAs8G3hMqZRJ9LCRdOjNiVn9QSI/v2Wz4Rg0XE/3SqNEgn2pMeZBY4wn8PEQoCtgKqUy+cSoTL7rSIVM/fm1tj9SIgu0B0VkNvB44PPvA4eSNySlomVyKiDXFyAz8efP9dRaLInM7K8BLgf2AV8CMwOPKZVSmZoKyOS7jlTIxJ8/11NrsSTS4vhz4JIUjEWppEtG1Uwm33WkQib+/JmaWkqnuMFeRH5HJ5umjDE3JGVESiVJMm/tOzviMBd0/PmD+fJ0laJmYmop3Tqb2Wuhu+pX+uTMW9WlVOfLY92t9eZgnP4qbrA3xjzc8TERGW6M2ZfcISmVHHprnxqpfFON98ZiDasIHYyTKamldOvuSVUvJGUUKqm0BM0vk2v1+5NUvql2thCbqQv66dLdYC9JGYVKmtDM58HzsS+rxPHg+fj278jJgJ+JVSP9USrfVPVuLXHdDfarkzIKlTRagnZUb868VV0L3kFaGHxXrE3Jm6rerSWuW8cSGmOWd/0slUn648ynN+WTnVXN5Hozs96Iyp2PvRhz1TP4sOFLYr5cF2ITl+ozaFWK9bcStGRVeuiOy/gSeROMWpTd+Tyy/318Se5Kmok1/pmqu2kclWX6W546WWkpTXfFluiaTzrvIHUhNjE6s+/n+tvMJ1lBpT+mu/pComWU/e0Osj/SmX0O6E8zn2QtyOlCX2yJvgn25R2klgonh87sVVbp7oJcoouuutAXW6Iz9r66g9S1k+TRYwlV1jkawDsPKt09Di/R180lPq8Xs38HtrDfIdOX4y0cjgw5vs9/P5l+/GSm69WxhEplmkSbjnV3236uNzOLxbLZ8AwYDhffDY4CaKmFjbdjazxA+9xXoKi0T8tVde0keTTYq37LwsCFvwZXCRgf+Lxgc2A3bnxeb+9LNZsPY/O2gM+Hx5bff+8E2pth7WVRD1ted5+nXHShN3mStkArIj8VkXdF5IiI1IjIsyIyLlnXUyqcz+tFmmvgpX+DV28DBJ69AR66GHno4l61jPB5vZhDf8N+8CPkoYuRZafiePB8TD9tQxFv8Vosq8/LVbO1VDgbFpWTWY0zFVgOTAbOBTzAqyIyOInXVFnI5zPUNLTxRW0zNQ1t+Hy9X0fyNtZgrbvSH4im/AieXtRnQcnbWIOt/tOo17T107r8eAHYiK3PUy7Z2NIiW/pPJS2NY4y5MPxzEZkD1ANTgGeTdV2VXXw+w879Dcx7ZDN7a1soL3Gx+qoJjB02AMvqed+9iNyvq6TbQcnnMxxqcuP2eHHabQwpdIbGI942f/46R3LL8SptvI012JKQcsm2tZNsOSchlXX2AwLXq03hNVWGO9TkDgV6gL21Lcx7ZDOHmnoXNCNSDy213aqhD74BXbr8Dabc+RqXLn+DnfsbQnccxpbnz2PnUF1+rL0a2Zpy6WvZsqicymB/L/BXYFOsL4rIfBHZLCKba2r6362wis3t8YYCfdDe2hZ8Pl+vUjsRgeiN/4Dpy6OCkhQOBcDj8fH3uhZ2H2ri73UtHGxsi3gDKi3KY199K3tqmznQ0IoUDsVbPDrqNb1hr5kLsjHlkgzZsiEvJdU4IvJb4BvAN4wxMRNZxphVwCrw19mnYlwq/Zx2G+UlroiAf0FFGQeb3CxYU93z1I5YNBWPIf/ql7F8bnYfAc+0pyi0eTnQbFj16hF+eamXEpfFR/sbWPjo0Wut+aczQuMZP7KYn1w4llue3HZ0LHMmcMygkdQ0ORk5+zkwXj6t93HPq4dZcqmHsgG5U+SW6pRLZ+m1dMmWDXlJ/1spIvcA3wPOMcb8LdnXU9llSKGT1VdNiMjZ//ziCmb9n7ejUjtPLZpC6YC8Ll8zfB3g7stOBZxcseqtqOf94jteDjT6QoE+eK3PDjaH3oAWTj0hFOhDY1mzmXXzJ/IPK7ZHvebPpvl68dtQnUnW+k5vZUv/qaSmcUTkXmAWcK4x5qNkXktlJ8sSxg4bwFOLpvD2T89l/YJJGGDxtArGjywOPW9vbQtuT2LVDeHrAHUt7TS7vZSXuCKeU17iwmm30e71RaWRlm3cxYrZVZSXuCh2OWKmmbw+E/M1bXqWW9Ika32nL2RD/6lk1tnfD8wFvg/UisjwwEdRsq6pspNlCUMKnRxubufylZs4e+nrLHluBz+5cGwo4AeDcyLC1wFWvP4JJYUOls6sDAXn8hIXK2ZXUeJy4LBZUUG7prGNsgF5rF8wiWOKXbGDuiXcP2t8xGsunVmJy5l5/8j7i3jrO4lOAnJdMmf2i/BX4GwEvgz7+EkSr6myVKxZ2y1PbmPh1BNCt+tDChNb8AquA4wfWczCqSdgDBxb4uIP8yfy6r+exR3fPYVlGz+mtqWdsqK80Cwejr4RDC5wMnxgPnaLqDeKO2dU8r+e/QCbZXHP5aexbv5Elkwfx7CB+RS7MmtRLtP0Zk9F8M81XHmJC4ddm/cmQhuhqZSKt8D2RW0zU+58Ler5f75pKi6nvVsLcT6f4bNDTew/0spNG44urN592anc8eJHbN1T53/tm8/BJuCw+YNFa7sXu82irCgPeyCAeDw+PjvcBAhNbR4KnDZa2738vb6VJ6v3cPsl47AJeA0YYzJm0TAT9TbnHuv7l86sZNjAfI4bUqi/c7QRmsoQnf1jj1WVU17iwuW0d7koG+sNpCjfzlW/j1xY/fET77F4WgUL1lRTXuLikwONzH3oXcpLXKycXcWIYv/MPDxo2O0Wxw8tYn9DC4cavVz32JaINw+nTahpdGfcomEmipdzT3Th3bKEYQPzWDJ9HAVOG3Ut7fzmjzupaWxL+DVymQZ7lTKd/WOPVZWTSOom3hvIwHx7zPxuscsRmhE+teULVs6potjl4EBDG81uL0X5bVGB2rIEjPDjJ96LevNYN39irwJYLumLnHuL28vch96N+dqqc5rsUinT2T/28KqcN245h6cWTUlodhzvDUREYuZ3y0tcrJs/kae2fMH08cey5LkdXLHqLRY/vZ18h8U9r+yMWd3R7ouu2tlb24LHZ3TRMEHxcu6JLrz31WvkKg32qtcSXXTr6h+qZQmlA/I4tqSA0gF5CaVB4r2B2ARWXzUhYmF19VUTGDHIf72LThkRVT//w7VbmFE1MmagtsV587BbsR/vb8GnL5rVBe/eOv6ZJLrw3levkas0jaN6pTuLbj1N1XQmXq7fsqzQnULHxeAhhU5GDy2M+SYxpNAZCtThawEOm8XSmZURC75LZ1ZSkGexYnZVxA7cYFlnx9fI1sXbvtrMFH731tPfR1+8Rq7SahzVKzUNbVy6/I2oYBsvZ93Xwa+ngehAQyvfXf5m1Lgfu/ZMyksKACJe94KKMm7+1knsOdxCgdNGs9vLV4YUUFzg4NYntzGjaiTFLgd1Le08Wb2HX11ayZBCZ0bu+Oyu7v4Zq/TRahyVNN1ddAumaoKC6YFkz/Q6vskMdkXfZSy/8nQk8G0d1wJe3nEAgNsvGYcxBofdwm4JDa0eZlSNZMXrn4RKOgFu+46319UnmUI3M/UPGuxVr8RLoySSs+7L9EBnwTPedUYMiizju+3pD6hpbGP9gkkx2yi8vOMAP59myLMJNQ1tEY3a7pxRyV0v7WTrnjouqChDRGh2e/pFkOzNn3Ei+kOqKxvoAq3qld4smPW210mii4YdrxNsWVzf4sHt9XHHix+xYE01W/fUsbe2hb/XtfDRvoaYC6+fHGjkr3vrQ4E+OO7gbt8LKsq44bwTuXzlprivkW2Lt8lcFO3q7ADVd3Rmr3qlNwtmvUkPdOeuIPw6sVoWh8/Ky0tcHGpys+L1T7hzRmXE85bOrOQ3f9zJrRedFHPcJw8fwO2XjOPylZvYW9sS8zWysXIkmYuiqUh16Z2DnwZ71WtdpVHi6Wl6wOcz7DvSSlObh8XTKkL58o5BIviP3GsMD179dZZt3BWzZfEtT25j8bQKljy3IyLw3/XSThZPq+Dk4QOwW8Luw83cetFJDC50ckFFWSiPHxy3y2mPeGMJf42Thg/AEWjFkAmBprsBsKd/xl1J9npAprZFTgdN46i06Ul6IPiP9/KVm5i5YlNEd8zwIBGeHjjrN6+z+Ont3PytsZQNyIs7K1+/YBIPv/lpaKF16546ljy3g6J8O7Ut7dy0YRtXrHqLuQ+9y/XnjuGCirKocXfcSxB8jY/2+ce8q6Yx7SmKTEqdJHuTVCa3RU41DfYqbXqya7ar7pjBIBHreTdt2EbZgLyYwcXltDN8YD43nj826s3H4zNROfpFa7fwi+98LWrcsd7A7pxRyYrXP8mYQJNJATDZm6S0kugoTeOotOpueiDeP95g0AgGiXjPs1sSc2OXZRn21jUzMN/Of/1wMsYYLMtiSKGTL+tbYr6WAMcGavLDf57gG1iL28OH+xpCaaHg96U70KQqACaSKkr2JqlkVxJlEw32KqvE+8d7TLGL4QPzQ0Ei0Z21DrtFY6uH6fe9GaNtrv/1uhswgm9gNQ2w5LkdGRdoUhEAu5MrT9Z6ACRn13a20jSOyirxbvvDA31nzwtvmeC022j3+Ljq9+9EpXt2H2oOpTVKXI6YB5yUuBz4fIbDTf7yz88PN3GgoTWU+87UPi6pGFempIp62mCvP9KZvcposVIBid72DxuYx7r5E/EayHdYDC30V8J0PJA8VkqjwGkLpTVqW9pZtvFjFk+rCLVEWLbxY3556Sk0tnpoaG3ncFP70TYKgws4bmhhxvZxScW4MilXnsw7h2yiwV5lnPAA7/UZfvn8Dl7ecSAiFdDdHbOPXHMGAO0eHyLCPa/sZG+t/0DyWCmNZrc3lNZwe7y8vONARKklwM+n+TjU6Mbj87H46e0RbRfqW9soKcgPBZrgz/RlfUsouAJpq/9OdgDUXHnm0WCvMkqsQH3njEpqGtwxa+ljibVjdv+R1tDJVeGvueL1T7j7slNDB5MEc/b+AO2jpqENh92KGbhsAkOLnMzpkAZatHYLa689kwFOH7Ut7fh8Pg42uSPaK6y+agJ5diuUQuov9d/hb9SPXXtm1Bt1ulNYuUy7XqqMEq/DYvA4QYA3bjknqgomXMfzbFfOqYq5UPrg1V/ncJMbnzG4HDYGFzqxWYIlcNszH4SC1CPXnEGbxxezt05Dq5e/1/nvEMKboT193RTsNmHBmurQhq2O118yfVzEqUvZ3kky1hv1yjlVDC10hiqbsvmNLBto10uVNeLleosD/eETSQV0TCEUuxwxX7O+pZ0rVr1FeYmLB2ZX+fP0XsNtz2wPpWxKi/LYfaiZk0YUsX7BpNCh4iUuB7sONDJvTeQdyF0v+c9ELXDamPvQu6Gxx1sX6PhYussyeyPWouyCNdVZ/QbWn2iwVxklXq43mFtPJBXQsdyu2e2N+ZqHmtyMH1nMwqkn0Nbupba5naUvfcQPJo+mpsFfNRLsozP5+CHMP/sEHDZBvD7qW93sO9LK3ZedGprV3/LkNpZMH4fTbuEOdM0cP7KYwYXOuOsC4dKd0+5tD5lMWpRV0TSNozJKX6UCwgOXy2lj/5G2qHWAp7f6z6Ht2BTt4Tc/ZUbVSMBfJz/5+CHMnvQVFq3dEjrI5IbzTow4nSo4q196WSVPvPs5syeN5nCTmwH5dta9s5uzxg6LuM7KOVXk2S2ufvDdjMjZd1UXn8gbgR5ykn6dpXE02KuMk4wuheGvKSLc/sx2ZlSNjJlLX3PNGSDg8fr428FmKkYM5Pur3wo9L94awJLp4zhxWBF1ze0s6PBG8PTWLzivYhhDCp0MH5TPvz/7ATUNbm44bwyjhxZSkGcLlYamQ2eBOtETt7TpWPp1FuyTuqlKRM4SkWdE5AsRMSJydTKvp/qHnhw8nuhrjhjkIt9h8S/nnciQQmfMtMOBhjb+9wsf0uoxLHluR9RBJvFy8F8ZUoDPmFCgDz5+y5PbOK9iGAvWVDNzxSb21bfy8o4DbN1Tx9yH3mX2f76Nx5veSVdnKZhEN0j1xQamjmcUeDy+Xh90rvySnbMvArYDjwQ+lOqV3sz6w2eepUV5LL3s1LjrAzOqRvLDQND2+kzE8+LV5v+9roV8h63LBeaOQXJvrf/AlPqW9rTNgjuri+9OLr439fsd7wxipcv0TqHnkjqzN8a8YIz5N2PMBsCXzGup/i9Wa94P9x3B40nsr1b4DHXrnjpueuI9ls6sjNmhMnz2vvrPf+OBK08PPe/J6j080KF9wp0zKrn75Y851OSO2VUz+AZx/6zTebJ6T9TXg2NLV0fMzlooJLsNcVDHO4gZVSNDgR5yuz1xX9BqHJU14pX2PXbtmZSXFMSc7YXfCYC/lDL4/Vv31PGbP+5kzT+dgc0S3B5DQ2s7C6eegM8cnc3vOtCIwy48NPcMLAGfgcI8i8fnTWT/kVYONblDnS2frN7DitlVEbPR+2edjiWweFoFa9/azQ8mj2bHlw1Ri7vprFzprIVCqpqJdbyDiJcu6/g70pOoEpMxwV5E5gPzAUaNGpXm0ahMFC+dcKChDZfTHpU+iLVgGDxaMLj5qaaxjX31rVgiEbto77n8VH5/9QSueWgzC6eewLUPV0elOO65/DTynbbQYm15iYu5U0bzyJufsWT6OEYPLSTPbkXU7QPUtbj5w/yJHGr0l2+GH4nY150n++I0qlT1+OmYSoqXLgv/HemicOIypuulMWaVMWaCMWZCaWlpuoejMlC8dEL4zB2OLvJ9Wd8S8wCTG84bE/relXOqkLBAH3zejevf40iLJ3SkYKw3GZ8x/G7jx6ybP5E/3zSVJdPH8Zs/7mR99V6WbdzFpwebaPf6+NnFFRGnWv1g8mj+/dkPKClwsOS5HaFA35ez5d6cRhXrIPdkLJp31DGVFLxL6qw7Z6Z018wGGTOzV6orQwqdrJxTFdFjJlgXf/qoSiByphevo+UJZUW8ccs5oZ2we+tiH07i9fmrcdYvmBQxwxw/spgbzhvD0AF53HThSYiAgVDrg1iHmq+YXcXiaRW0tvtY/ee/8fKOA9z2na8lbbbc04O80zlTjnUHUeJydPo70o1cidNgr7KGZQljywbw2LVncqChjUNNbh5+81NuPH9sRBfJYKCKlwZwOWwRAa8gL3YlSrPby+o5E/D6fCydWclNG7ZRWpTHzd8ay00bjgby+2edzpAiB+UlLkqL8nE7c1oAABqASURBVPjNzMpQqwTwB5+Fj1ZHHGpe1+LGabclbbNRV0EwXornYFNbj94kwvUmhx4rldTZdbW7ZuKSGuxFpAj4auBTCxglIqcBh40xnyfz2qp/ststyksKcDntjBiUz+mjKiOCSXiQW/H6J9w5ozJq52rHVMnQwryoBciVs6sYUZwPwCX3vUFpUR53fPcUvjKkMGKD1d7aFq57bAt3fPcU7p81nma3F7fXF9H7PtggbUxZEYunVfDwm5/y84srktoBsrMgGG/2Pqa0iOa23s2UO7szgL5v6awnUSUuqTtoRWQq8FqMLz1sjLk63vfpDlrVHR13x16+clNUymXkYBef1DRxWvkghg1ydfoaDruF3RJaAr1rrn9sa2hB96lFk7l0+ZtR379u/kTavT4e2fQZ150zhuse2xKVarrpwpOob2lncKGTgfl2hg7IT+rvJF7QPdTkjrlbdv2CSezc1xDqzR/+tURn9vF24j5z/ZSolhV9lR7Sapyj0tb10hjzOpCbv3WVErE24oSXPtY0tuG0W9z0xDZqGtt4atGUmK8TfshIZxU8Bxra4m7EKnY5mFE1MhTo4egO2rXXnsmvwnq7r5xdxeAktkewLGFMqb9Tp8frw26zKCvyB2u3xxvRwG3rnjr21rbQ7vXx4vtfsvzK00N9gIIdQUsCm8K64vZ4KS3Ki7qzaXF7e50e6uxn1d47XdOcvcpqHRcigyWO6xdMwu3x8enBplDb4dVXTaDE5QhUmPjwGkIti4OzwVgLmzdt2MYf5p9Ju8dgICpwB+vkF049IW4Lhrrm9tDY9ta2sODR5Lb+9fkMu2oao07r6tiXP7wts8NmcdEpI7jvT7sigvXvNn7Mry6tjDhxK94s2uW0Ra1pLJ1Zic8YXUhNMw32KqvFWoj0V7oYRg0uoDDPzn2zxh/tQV/TyD2v7OQHk0dH5PKDKYVYr1dalEddsydio9TKOVX8+/RxtHl8/Or5HaENVT+7uCLmzH/fkdaI10x2oIv1prX7UHNEiiZ413HXZacyyOWgrCgP91BfzCMYb/uON3ZH0sDaRrHLH/Q9PhMK9MFr3LRhW1RFU/D3ogupqZMxdfZK9URnW/k71obXtrQz75HNzKgaGQr0EFmbHev1bjhvTNS2/QVrqnHaLZw2i5u/dTKv/fhsfjGtgme2fsHysNYKwbLLWC0SgoulyWj0FetNq7gg9o7UYQPzGVNahN1uhSqTYo015g7mR6t5b099qIa/3eOLeQ2bELcdg0oNDfYqq8Xq6bJyTlXo/Njw4BkMgGUD8uKmFIYUOlnZYSPPqCEFMWf7X9a1ctnKTdz0xHt8dqiZdp/hmyeW4rQLS6aP46lFk3l83kQ2f3qQ684ZEznG2VUU59t7vPGpKx3ftMaPLGaQyxEzkH92sInalnbgaGVSrKAcr5yzwGnr9M2yvMSFZVm97oipekfTOCqrddyI4/WZ0CHXF1SU8fOLK7BZgtNuw+X0B6KiPHvMlIKIYFnC8OI8lkwfR4HTRl1LO1/WtUQ9/4bzxrDg0WpKi/KiNlAtnVnJso272Lqnjsurypk96SuhPPiQQieDC52se2c3l5xWzr0bP45Y0NxX38qwgXkMLuxdLr9jSeIN543hjhc/jCpFvX/W6dz+zAfcN2t8zN9neF7eabdxQUUZM6pGhvL5T1bvoa6lPfRmOWKQK24ppC6kppcGe5X1gkEkvOxv/MhifjB5NLP+z9sRQeeRa86gpd0bFfTunFGJLTDJFIR8hxVqoXBBRVlUhcpxQ/2z/cXTKqJSQjdt2BY6IH199V52HWhk2ffHh5qm/er5D9m6p47nt+8PLF4SuRdgdlUoB95d4QuoQ4qcPHP9FFrcXrzG8PKOA9Q0uCMWXy0hULFki/p+EcHn83Goyc2QQiclLkdUy+HlV57Oo5t2R6TOUtFHR3WfBnvVb4SnGRZOPSFmXv6/Fk3G5bBx76sfRwS9h9/8lF9eegoALW4vv/njzoivP7ppN3+YPxGfgXyHhcfjo7zEFbczY3FYqWJNYxsAM1dsinreMcUurgy8IQUf62mlTqwF1EeuOYOifDvedi8PXv11lm3cxYI11cDRZm6PXXsmbo+XAw2tNLZ6uOr370TtEbjx/LEMG5gXtXaxaO0Wlkwfx9xvjA7l33UGn5k02Kt+I3zXaLwg3O7x4XLamDtldKhq5IKKMm696GTa2r3UNLThctqoaWwLBUXwB8Yb5UScdqHd48Nhs1g6szLuYeYm7P/vvuxUwEQ974KKsj4tSey4gDr5+CE0tkUG7+CegZrGNh6a+3XaPL6Iu5+lMytDbaCD1TqLp1Uw75HNrJs/MeZYTygrorzYpbP3DKcLtKrfCF+sDfbFCRdMNYTP3J+9fgrXnTOGq37/Dt/8zetcuvwN9h9p45FrzohYpLz7slM50NDKd5e/yZQ7X+Pn//0+x5a4OKGsMGpBd+nMSgYXOnj2+iksnlaBz5hQT5zw5/3s4go+O9jcZweDhN/ZjB9ZzMKpJ7Bo7RZKi/JYOaeKuy87ldZ2Hw/MPp0l08exr7411FQOjqagFk49IfSa4W+cXkPMsbocNg30WUCDveo3wvPFp5UPYuWc2O1xnfajM/e/17dG7Xid98hmivLtrF8wiQ0LJ4UC9nWPbQ09r6bBzf76Vmatfhu318uS6eNYN38ii6dV8Js/7uSahzZT39LOgjXVWCK8vOMAd720k6UzK3ntx2fz0NyvY7OE7Xvroko1V86uwmbR7aqc8EqYm781Fp8htIC85LkdXLHqLRY/vZ2ahjaWbdyFw2Z1mYIKf+PMd1haPpnFNI2j+pXwfHHpgPwuT17qLN0DR/PswRTG+JHF3PytsZSXFLD/SKu/wqYoj0uXR+bjAUYEevAEUz3gP+VqTiCtckFFGT+5cCwtbm+o+qfZ7cVhF259chs3nj+2W+WJwZ/rnld2UpRnZ8/hZm44b0zU2sUP125h8bSKuF1BmwM9gcJz9quvmoDdEgYG3gRtApZl6eJrFtFgr/qtRE5eavPEzrmLCA6bFfpaXUs7F1SUseicr9Li9oY6XwY3TV1QURax69SfirFYN38iPmO4+7JTaXZ7IwLvjKqRfFHbGrPxWDBP3p2F2mA/nNu+8zW+rG8l32ExbGB+3Nn7HS9+xP2zTo9o2vbA7CoGuez8z83nYLMEm8AvLz2FxlYPl9z3Rp83MVOpo2kclZOCbwQOS6Jy6cEyzPA1gBWvf8KtF51MbVN7VDuAhY9Wc+tFJ0e8xvIrT8flsCjKt3Ok1cOfPtzH6NLCiMBb7HJQ4LTFDcbdXagN9sO5YtVbzFyxiZs2bMMSiZlnN/jP4H1h2xc8Pm8if/rx2SyZPo5f/Pd2Zq1+m8Y2D8MH5lM6IB+P14QWeYPj09Ogso/O7FVOsyyLh9/8NKoM81eXVkbcAfh8Pto8vrjB2bKER645g8NNbg41ubnvT7uY983jKS5wMLjAycwJI3HahAev/npos5bPGFrbfTHvLAYXOrmgoqxbC7V1LW721bdGdLT89Qs7ovYILJ1ZydAiJ2//9Fwsy8IYE9GjH2DeI5t5+vrJHDjipqnNo03M+gEN9iqnDSl0cuP5Yzs9/MJgaGn34fGZuKWWdkv4Xli9PMCOLxtYMn0cAE674PaaUMrGX+N+KqNLC0KnYIVvVHIHzq4dlJdYsPf5DF/WtUa8frCjZUmBgyXTx1Fc4KAoz05ru5c9h1s4deQgBhfm8UVtc+y2xG3+DpmLp8Vu7qZNzLJLUg8v6Sk9vESlUry2vbF65f/kwrEcanRHBOeVs6tA4OJlf4l67XXzJwL+RdpYufm1156JiH/Xrs8YjIFfvxDZ9/7kEQO7zI3HOzTkvu+Pp3RAHm6Pv6f9o5s+ZeX/fBbx2nUtbnbua4hqSzxsYD7n3v1/Y56pu3JOFScP73pcKrXSdniJUtkg3kJuvF75P7u4gjXXnIHN8i/iHm52h4JgrMqWfIfFyMGumKmQmoY2Zq7YFJrR3/enXT3qex+vSdmgAidXrHor4q6httnD+uq9odc2xG5L/HBgr8HWPXXc9dLOUG+fQS4HA/PtGuizjC7QKhVHvF75++pbmfP7d3A57aEdsE9W7+Huy06N2lw1ojgPS4Q9h1tiLpQGFzmDrQdmVI2MeM7e2hZa2z1d1twHm5StnFPFuvkTWTmnigXfPA6vz18JtHJOFaVFeSxau4V5Zx0fem23xxu3LXFruze0YWzrnjqWPLeDNo+PpS99hGVp6Mg2OrNXKo54h3Y3u72hvP6hJjdPVu/hB5NH859/+Rt3fPcUjil24bBZ7Ktv5WCDm1v/631Ki/J44MrT+eHaLVE59aC9tS1RG5TKS1x4fIYP9x1hbNkA7HYrqlmZTSDPYfHP553IDwO9axZ88zhmTTyOmoa2UHfKn1w4lqe3fkG+w8a6+RNpdntxOW14fcT8OQfk2RkxMJ/Hrj2TAw1tHGpyh/rk6Eaq7KPBXqk4OrYJjnUyU3CB955XdjKjaiQD8h047RYFTotjivPx+Eyoz0xjmye0eWpwoZOlL30UOsgc/AG2bGAeTyyYxNAiJ14DDpsgAve++jE/v7iCYwa5QscNlhblccN5Yxg1pICWdl8o0I8fWcy3K48NNVgL3xx160UnR+wRWH3VBMaUFkX/nHOqOGaQC7vdorykAJfTzohB+Zw+qlI3UmUpXaBVqhNdnbna1XM6tl0OLnSWFuXxb98+iRvXvxexKPqVoQW0tfs41OgO3TX883kn4vP5aPcajil2cfsz25k7ZTRF+Y5QgN+wcBK/ev5DFk49gRPLikK7dIOCG7WGFDojum+Wl7h4atGU0F2KtiXObrpAq1QPJdKut7PnlLgcrJhdxcJHq9m6p46H3/yUR645g/qWdgYXOiMOSXlqyxfMnvSViJr4O2dU8ruNH7N42tfYW9sMGH4weTQlBQ68Pnh83kS8PoPdEn57RSUHG9qxLImZgw8G9I6Puz1ebUucAzTYK5VEtS3tLNvo751fNiCPoUV5LHnuA17ecYB18ycy96F3Q89dOacq1KUyWPPe7PYyd8povMYwemghXp9h1756xlQew/4Gd0RK5oHZVQx02fhbTVPMHHzpgDxW/fmTiPFpvXzu0CV1pZLI7fHy8o4DLFhTzaXL3+R3G3eFzqPt2Ia52OWI2aWyKM9Onl042OimKN/iGyeW0eYxoRQOBBqcPVqNzbKxbOOuqBYQ9886ndc+3Mf150aehatdK3NH0mf2IrIIuAkYAXwA/MgY8z/Jvq5SmaBjRc/66r3Utbh5aO4ZFDgtVs6uYkEgaDe7vXG7VD4+byLP/nUv//SN46lpaOOY4th1+zaLiLr4YpeDY4pdHGltZ8qYMjZs/pz1CyZhjNHcfI5J6sxeRK4A7gV+DYwH3gReFJFRybyuUpkivJka+GfT/3LeiQwudDB8oIvywfmsvfZMNiycRL7DCp1tG25vbQsNre3MmTyaFo+P1nYfdit2gzN7oP596546Fqyp5sdPvMeOL49w8bK/UN/Szj+ePpLhA/M5tqSA0gF5GuhzSLJn9v8KPGSMWR34/J9F5FvAD4GfJvnaSqVdVwdwD8hzcqTFw5AiJyA4bBIz316YZ+d7YTth182fGPPQdEuO1syXl7j4jytO41fPf0h5iYtjil0MH5ivAT5HJS3Yi4gTqALu6vCll4HJybquUpmms0oXyxKOLS7gYFMbPp/hSGt7qHonvF/+r57fEZHa+Xh/I4+/szuqW+dt3/kaj8/z99CvbXLzv57dQU1jG6uvmqCBPsclc2Y/FLAB+zs8vh/4h45PFpH5wHyAUaM0y6Nyh2UJZQPy8Xh8HGpyc7ixhbXXnklNYNfqIJc94mAUgGUbd/Hv078WsSN31Zwq7JZgWRYlLgcFTjv3zRqvuXkFpKb0suOuLYnxGMaYVcAq8G+qSsG4lMoodrvFycMHcrAoD5sF4L8bMESndmoa2ygdkMcfAidh5dttDC2KzMFr3bwKl8wF2oOAFxje4fEyomf7SikCs/yB+ZQU5DHIZeeYQfnk24UHZkcenv5A4FByp92ivLiAMk3RqC4kbWZvjHGLSDVwPvBE2JfOB55M1nWV6g8sSyguyKO4wN+OwWkX1s2fiCewW7bYZeFyajWNSlyy0zi/BdaIyDvAG8BC4BhgRZKvq1S/4Q/8+RQXpHskKpslNdgbY9aJyBDg5/g3VW0Hvm2M2Z3M6yqllIqU9AVaY8xyYHmyr6OUUio+7Y2jlFI5QIO9UkrlAA32SimVAzLypCoRqQF6uog7FH+Nf7bRcaeWjju1dNyp8RVjTGmsL2RksO8NEdkc71iuTKbjTi0dd2rpuNNP0zhKKZUDNNgrpVQO6I/BflW6B9BDOu7U0nGnlo47zfpdzl4ppVS0/jizV0op1YEGe6WUygEa7JVSKgf0+2Avfn8UESMiM9M9nq6IyGoR+UREWkSkRkSeFpGT0z2uzojIYBH5nYh8FBj3HhF5INDxNKOJyHwReU1E6gJ/R45L95hiEZFFIvKpiLSKSLWIfDPdY+qKiJwlIs+IyBeB3+3V6R5TV0TkpyLyrogcCfz7e1ZExqV7XH2h3wd74Mf4T8zKFpuBq4GTgQvxH+P4qog40jmoLhwDHAvcDJwCzAbOAh5P56ASVAC8DNye5nHEJSJXAPcCvwbGA28CL4pIph/WXIS/rfm/AC1dPDdTTMXfpXcycC7gwf/vb3A6B9UX+nU1johMAJ4CqvAfhXiZMWZDekfVPSJSCbwHnGSM2Znu8SRKRL4NPAcUG2OOpHs8XQn8XXkXGG2M+SzNw4kgIm8D24wx88Ie2wVsMMb8NH0jS5yINALXG2MeSvdYukNEioB64B+NMc+mezy90W9n9iIyAP/McoEx5kC6x9MTIlIIzAU+Bz5L72i6bSDQBjSneyDZTESc+CcrL3f40sv4Z58quQbgj5O16R5Ib/XbYI//6MM/GmNeSPdAuiuQn20EGoGLgPOMMW1pHlbCRKQYWAKsNsZ40j2eLDcUsOG/Mw23Hxie+uHknHuBvwKb0j2Q3sqqYC8ivwws9HT2MVVE5gCnAjele8yQ+LjDvmUt/tzs2cDHwBMikvITSHsw7uDdyLPAF/hz+CnXk3FngY75VonxmOpDIvJb4BvADGNMNq37xZRVOXsRGYp/ptOZz/EvsFwF+MIetwU+32SM+UZyRhhbouM2xkSlPAK38bXAQmPMmmSML57ujjuQ33wBfyC6yBjTmOQhxtST33em5uwDf/7NwPeNMU+EPX4/MM4Yc3baBtcN2ZazF5F7gO8B5xhjPkr3ePpC0s+g7UvGmIMk0FtaRH4G3NXh4feBnwBPJ2FonUp03HFI4COv70aUmO6MO7BG8iL+sX4rXYEeev37zijGGLeIVAPnA0+Efel84Mn0jKp/E5F78Qf6qf0l0EOWBftEGWO+wJ9GCBERgD3GmL+lZVAJEJGvAjOAV4EaoBy4Ff9C53NpHFqnAoH+ZfyLsv8IFAbSOQCHjTHutA2uCyIyHH/u+8TAQxWBNYfPjTGH0zeyCL8F1ojIO8AbwEL85a4r0jqqLgTu9L4a+NQCRonIafj/TnyevpHFF7hjmoP/73Ft4O8HQGM6JzB9whiTEx/485sz0z2OLsY4Ev/s+ADgBvbgz9+flO6xdTHuqYHfb6yPqekeXxdjvz3OuK9O99g6jHMR/oqsNqAaOCvdY+rF34uH0j22TsYc7+/x7ekeW28/sipnr5RSqmeyqhpHKaVUz2iwV0qpHKDBXimlcoAGe6WUygEa7JVSKgdosFdKqRygwV5ltMA2+55+71oR2Ski20Xk9+FnAojIP4rILwL/f3vggI2/Bg5geUBErMDX7hKRc+O8/nEisr2HYysQkecD1/tARO7o8PUfichVInJ/YFw7xH8wzF8DHzM7G5tSHWmwV/3ZWuAk/AequIBrw752M/4eSkH3GGNOAyoCzw/2nPkd/l3MyXCXMeYk/E3vpojIRQAiYgeuAR4zxlwXGNe3gU+MMacFPjYkeWyqn9Fgr7KC+C0NzNLfF//pTYiIJSLLA7Pj50TkBQkcP2mMecEEAO/gbz+BiJwItBl/D52OnEA+gf7lxpjdwJCwbfPxxpcvIg8GxrZVRM4JPF4gIutFZJuIrBORt0VkgjGm2RjzWuAabmBLcHz4T0jaYrpoD53o2JQCDfYqe3wXOA1/6+p/AJaKyIjA48fhn41fC0zq+I2B9M0c4I+Bh6bgD67hbhSRvwJfAh8bY/4a9rUtge/pzHUAxphTgO8DD4tIPv42B7XGmEr8Pf6rYoyvGPgOsDFsfNVdXK87Y1NKg73KGt8AHjfGeI0x+4H/C3w98PgTxhifMWYf8FqM710O/NkY8z+Bz0fgbzQXLpjGKcPfyO17YV87gL/xWFfjWwNg/J0Sd+NvrvYN4A+Bx7cD28K/KZCyeRxYZo426Ys1vngSGZtSGuxV1pBuPu7/oshtQCnwr2EPt+BP1UQxxrTjvwM4K+zhfKBFRM4MWyC9pC/GB6wCdhlj/iOR8cWQT/Yc5q3SSIO9yhZ/Bq4QEZuIlOIPxu8AfwFmBHL3w/B3WgRARK4FLsR/8Ef4QTYfcrT1bgTx98KeDHwS9vCJwHZjzNthC6TPxBjflYHXOBEYBewMjO/ywOPBxd/gtX4JDAJ+1OG14o4vhhOBHlUEqdyiwV5li6fwp0DeA/4E3BxI2zwJ7MUf8FYCbwP1ge9ZAQwDNgVm478IPP5nYHwgsAcFc/bb8Z/zsBxC+f6vApu7GN9ywCYi7wPr8LdIbgs8Xioi24BbAj9DvYiUAz/DX/2zJTC+YLXQi0TeWcTUjbEppS2OVfYTkSJjTKOIDME/258SeCPo7HvuBZ41xrzaxfMuBU43xizu4dhsgMMY0yoiJ+BfhD3RdHGgi4g8hf8NbVeyxqZyS788qUrlnOcCFS1OYElXgT7g18CZCTzPDtzdi7EVAK8FZuEC/LCrQB9wK/6F2rjBvg/GpnKIzuyVUioHaM5eKaVygAZ7pZTKARrslVIqB2iwV0qpHKDBXimlcsD/B67ENWIMVarxAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "figure = sns.scatterplot(x='log2(B)-log2(T)',y='-log p',data=tstats)\n",
    "figure2 = sns.scatterplot(x='log2(B)-log2(T)',y='-log p',data=tstats.loc[altered.index])\n",
    "plt.rc('axes',labelsize=16)\n",
    "plt.xticks(fontsize=14)\n",
    "plt.yticks(fontsize=14)\n",
    "#plt.title(\"B cells vs T cells\", fontsize=20)\n",
    "\n",
    "#plt.savefig(\"data/BvsT.png\", dpi=300)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This shows the differentially expressed proteins. Recall that this dataset is restricted to those expressed in at least three of each cell type.\n",
    "\n",
    "The next step is to identify those proteins by name. Here we get the protein names from the fasta file, which remain the same regardless of which software the data was quantified using."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "protein_names = load_data.load_fasta()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "protein_names_by_group={}\n",
    "for ID in data_raw.index:\n",
    "    if ';' in ID:\n",
    "        ids = ID.split(';')\n",
    "        names=''\n",
    "        for i in ids:\n",
    "            if i in protein_names.index:\n",
    "                names=names + protein_names[i] + ';'\n",
    "        names = names[:-1]\n",
    "        protein_names_by_group[ID] = names\n",
    "    else:\n",
    "        if ID in protein_names.index:\n",
    "            protein_names_by_group[ID]=protein_names[ID]\n",
    "protein_names_by_group = pd.Series(list(protein_names_by_group.values()), index=protein_names_by_group.keys())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Here we add the names to the <code>tstats</code> dataframe, which contains the fold changes for each protein. Then we can list the proteins that are altered. Optionally, use <code> head()</code> to display the first few or <code>to_csv</code> to write out to a file. We now split up the proteins heightened in B cells from those heightened in T cells."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "tstats['Name']=protein_names_by_group\n",
    "tstats['GeneID'] = tstats.Name.apply(lambda n: str(n).split('_HUMAN')[0])\n",
    "\n",
    "altered_proteins = tstats.loc[altered.index].sort_values('log2(B)-log2(T)',\n",
    "                        ascending=False)[['pvalue','log2(B)-log2(T)','Name','GeneID']]\n",
    "higher_in_B_cells = altered_proteins.loc[altered_proteins['log2(B)-log2(T)'] >0]\n",
    "higher_in_T_cells = altered_proteins.loc[altered_proteins['log2(B)-log2(T)'] <0]\n",
    "\n",
    "#To write out to a file, use:\n",
    "#altered_proteins.to_csv('data/filename.tsv', sep='\\t')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Next we do the same for those only identified in one cell type."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [],
   "source": [
    "data_in_one_type['Name'] = protein_names_by_group\n",
    "data_in_one_type['GeneID']=data_in_one_type.Name.apply(lambda n: str(n).split('_HUMAN')[0])\n",
    "\n",
    "B_cell_only_geneset = data_in_one_type[(data_in_one_type.Type==' B_')]\n",
    "T_cell_only_geneset= data_in_one_type[(data_in_one_type.Type==' T_')]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Proteins Higher in B Cells:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>pvalue</th>\n",
       "      <th>log2(B)-log2(T)</th>\n",
       "      <th>Name</th>\n",
       "      <th>GeneID</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Protein IDs</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>P09960</th>\n",
       "      <td>0.000045</td>\n",
       "      <td>2.626130</td>\n",
       "      <td>LKHA4_HUMAN Leukotriene A-4 hydrolase</td>\n",
       "      <td>LKHA4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q96KP4</th>\n",
       "      <td>0.000029</td>\n",
       "      <td>2.567266</td>\n",
       "      <td>CNDP2_HUMAN Cytosolic non-specific dipeptidase</td>\n",
       "      <td>CNDP2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P29350</th>\n",
       "      <td>0.000048</td>\n",
       "      <td>2.302902</td>\n",
       "      <td>PTN6_HUMAN Tyrosine-protein phosphatase non-re...</td>\n",
       "      <td>PTN6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q16666;Q6K0P9</th>\n",
       "      <td>0.009769</td>\n",
       "      <td>2.267724</td>\n",
       "      <td>IF16_HUMAN Gamma-interferon-inducible protein ...</td>\n",
       "      <td>IF16</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P09874</th>\n",
       "      <td>0.000006</td>\n",
       "      <td>1.805240</td>\n",
       "      <td>PARP1_HUMAN Poly [ADP-ribose] polymerase 1</td>\n",
       "      <td>PARP1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P53396</th>\n",
       "      <td>0.002844</td>\n",
       "      <td>1.671055</td>\n",
       "      <td>ACLY_HUMAN ATP-citrate synthase</td>\n",
       "      <td>ACLY</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P58546</th>\n",
       "      <td>0.003921</td>\n",
       "      <td>1.543587</td>\n",
       "      <td>MTPN_HUMAN Myotrophin</td>\n",
       "      <td>MTPN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P50395</th>\n",
       "      <td>0.006946</td>\n",
       "      <td>1.291659</td>\n",
       "      <td>GDIB_HUMAN Rab GDP dissociation inhibitor beta</td>\n",
       "      <td>GDIB</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P15311</th>\n",
       "      <td>0.001631</td>\n",
       "      <td>1.049916</td>\n",
       "      <td>EZRI_HUMAN Ezrin</td>\n",
       "      <td>EZRI</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P51149</th>\n",
       "      <td>0.003642</td>\n",
       "      <td>0.968890</td>\n",
       "      <td>RAB7A_HUMAN Ras-related protein Rab-7a</td>\n",
       "      <td>RAB7A</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P22695</th>\n",
       "      <td>0.005077</td>\n",
       "      <td>0.901173</td>\n",
       "      <td>QCR2_HUMAN Cytochrome b-c1 complex subunit 2, ...</td>\n",
       "      <td>QCR2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P39687;O43423;O95626</th>\n",
       "      <td>0.002319</td>\n",
       "      <td>0.881030</td>\n",
       "      <td>AN32A_HUMAN Acidic leucine-rich nuclear phosph...</td>\n",
       "      <td>AN32A</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P50552</th>\n",
       "      <td>0.008336</td>\n",
       "      <td>0.821799</td>\n",
       "      <td>VASP_HUMAN Vasodilator-stimulated phosphoprotein</td>\n",
       "      <td>VASP</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P59666;P59665</th>\n",
       "      <td>0.007289</td>\n",
       "      <td>0.697647</td>\n",
       "      <td>DEF3_HUMAN Neutrophil defensin 3 ;DEF1_HUMAN N...</td>\n",
       "      <td>DEF3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9UJU6</th>\n",
       "      <td>0.005791</td>\n",
       "      <td>0.679712</td>\n",
       "      <td>DBNL_HUMAN Drebrin-like protein</td>\n",
       "      <td>DBNL</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                        pvalue  log2(B)-log2(T)  \\\n",
       "Protein IDs                                       \n",
       "P09960                0.000045         2.626130   \n",
       "Q96KP4                0.000029         2.567266   \n",
       "P29350                0.000048         2.302902   \n",
       "Q16666;Q6K0P9         0.009769         2.267724   \n",
       "P09874                0.000006         1.805240   \n",
       "P53396                0.002844         1.671055   \n",
       "P58546                0.003921         1.543587   \n",
       "P50395                0.006946         1.291659   \n",
       "P15311                0.001631         1.049916   \n",
       "P51149                0.003642         0.968890   \n",
       "P22695                0.005077         0.901173   \n",
       "P39687;O43423;O95626  0.002319         0.881030   \n",
       "P50552                0.008336         0.821799   \n",
       "P59666;P59665         0.007289         0.697647   \n",
       "Q9UJU6                0.005791         0.679712   \n",
       "\n",
       "                                                                   Name GeneID  \n",
       "Protein IDs                                                                     \n",
       "P09960                           LKHA4_HUMAN Leukotriene A-4 hydrolase   LKHA4  \n",
       "Q96KP4                  CNDP2_HUMAN Cytosolic non-specific dipeptidase   CNDP2  \n",
       "P29350                PTN6_HUMAN Tyrosine-protein phosphatase non-re...   PTN6  \n",
       "Q16666;Q6K0P9         IF16_HUMAN Gamma-interferon-inducible protein ...   IF16  \n",
       "P09874                      PARP1_HUMAN Poly [ADP-ribose] polymerase 1   PARP1  \n",
       "P53396                                 ACLY_HUMAN ATP-citrate synthase    ACLY  \n",
       "P58546                                           MTPN_HUMAN Myotrophin    MTPN  \n",
       "P50395                  GDIB_HUMAN Rab GDP dissociation inhibitor beta    GDIB  \n",
       "P15311                                                EZRI_HUMAN Ezrin    EZRI  \n",
       "P51149                          RAB7A_HUMAN Ras-related protein Rab-7a   RAB7A  \n",
       "P22695                QCR2_HUMAN Cytochrome b-c1 complex subunit 2, ...   QCR2  \n",
       "P39687;O43423;O95626  AN32A_HUMAN Acidic leucine-rich nuclear phosph...  AN32A  \n",
       "P50552                VASP_HUMAN Vasodilator-stimulated phosphoprotein    VASP  \n",
       "P59666;P59665         DEF3_HUMAN Neutrophil defensin 3 ;DEF1_HUMAN N...   DEF3  \n",
       "Q9UJU6                                 DBNL_HUMAN Drebrin-like protein    DBNL  "
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "higher_in_B_cells"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Proteins Only in B Cells:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Name</th>\n",
       "      <th>GeneID</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Protein IDs</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>A6NI72;P14598;A8MVU1</th>\n",
       "      <td>NCF1B_HUMAN Putative neutrophil cytosol factor...</td>\n",
       "      <td>NCF1B</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O00151</th>\n",
       "      <td>PDLI1_HUMAN PDZ and LIM domain protein 1</td>\n",
       "      <td>PDLI1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O15400</th>\n",
       "      <td>TX7_HUMAN Syntaxin-7</td>\n",
       "      <td>TX7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O60262</th>\n",
       "      <td>GBG7_HUMAN Guanine nucleotide-binding protein ...</td>\n",
       "      <td>GBG7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O60493</th>\n",
       "      <td>NX3_HUMAN Sorting nexin-3</td>\n",
       "      <td>NX3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P00352</th>\n",
       "      <td>AL1A1_HUMAN Retinal dehydrogenase 1</td>\n",
       "      <td>AL1A1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P01871;P0DOX6</th>\n",
       "      <td>IGHM_HUMAN Immunoglobulin heavy constant mu</td>\n",
       "      <td>IGHM</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P01903;P01906;P52735</th>\n",
       "      <td>DRA_HUMAN HLA class II histocompatibility anti...</td>\n",
       "      <td>DRA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P01911</th>\n",
       "      <td>2B1F_HUMAN HLA class II histocompatibility ant...</td>\n",
       "      <td>2B1F</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P04233</th>\n",
       "      <td>HG2A_HUMAN HLA class II histocompatibility ant...</td>\n",
       "      <td>HG2A</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P04440</th>\n",
       "      <td>DPB1_HUMAN HLA class II histocompatibility ant...</td>\n",
       "      <td>DPB1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P04839</th>\n",
       "      <td>CY24B_HUMAN Cytochrome b-245 heavy chain</td>\n",
       "      <td>CY24B</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P09917</th>\n",
       "      <td>LOX5_HUMAN Arachidonate 5-lipoxygenase</td>\n",
       "      <td>LOX5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P11586</th>\n",
       "      <td>C1TC_HUMAN C-1-tetrahydrofolate synthase, cyto...</td>\n",
       "      <td>C1TC</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P11836</th>\n",
       "      <td>CD20_HUMAN B-lymphocyte antigen CD20</td>\n",
       "      <td>CD20</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P13284</th>\n",
       "      <td>GILT_HUMAN Gamma-interferon-inducible lysosoma...</td>\n",
       "      <td>GILT</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P13762</th>\n",
       "      <td>DRB4_HUMAN HLA class II histocompatibility ant...</td>\n",
       "      <td>DRB4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P16885</th>\n",
       "      <td>PLCG2_HUMAN 1-phosphatidylinositol 4,5-bisphos...</td>\n",
       "      <td>PLCG2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P27707</th>\n",
       "      <td>DCK_HUMAN Deoxycytidine kinase</td>\n",
       "      <td>DCK</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P28072</th>\n",
       "      <td>PSB6_HUMAN Proteasome subunit beta type-6</td>\n",
       "      <td>PSB6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P34897</th>\n",
       "      <td>GLYM_HUMAN Serine hydroxymethyltransferase, mi...</td>\n",
       "      <td>GLYM</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P49006</th>\n",
       "      <td>MRP_HUMAN MARCKS-related protein</td>\n",
       "      <td>MRP</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P55327</th>\n",
       "      <td>TPD52_HUMAN Tumor protein D52</td>\n",
       "      <td>TPD52</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P56279</th>\n",
       "      <td>TCL1A_HUMAN T-cell leukemia/lymphoma protein 1A</td>\n",
       "      <td>TCL1A</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P79483</th>\n",
       "      <td>DRB3_HUMAN HLA class II histocompatibility ant...</td>\n",
       "      <td>DRB3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P80723</th>\n",
       "      <td>BASP1_HUMAN Brain acid soluble protein 1</td>\n",
       "      <td>BASP1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q06187</th>\n",
       "      <td>BTK_HUMAN Tyrosine-protein kinase BTK</td>\n",
       "      <td>BTK</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q08209;P48454;P16298</th>\n",
       "      <td>PP2BA_HUMAN Serine/threonine-protein phosphata...</td>\n",
       "      <td>PP2BA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q14011</th>\n",
       "      <td>CIRBP_HUMAN Cold-inducible RNA-binding protein</td>\n",
       "      <td>CIRBP</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q6JBY9</th>\n",
       "      <td>CPZIP_HUMAN CapZ-interacting protein</td>\n",
       "      <td>CPZIP</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q6NUK1</th>\n",
       "      <td>CMC1_HUMAN Calcium-binding mitochondrial carri...</td>\n",
       "      <td>CMC1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q8NBS9</th>\n",
       "      <td>TXND5_HUMAN Thioredoxin domain-containing prot...</td>\n",
       "      <td>TXND5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q96D96</th>\n",
       "      <td>HVCN1_HUMAN Voltage-gated hydrogen channel 1</td>\n",
       "      <td>HVCN1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9H8H3</th>\n",
       "      <td>MET7A_HUMAN Methyltransferase-like protein 7A</td>\n",
       "      <td>MET7A</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9NY33</th>\n",
       "      <td>DPP3_HUMAN Dipeptidyl peptidase 3</td>\n",
       "      <td>DPP3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9UH65</th>\n",
       "      <td>WP70_HUMAN Switch-associated protein 70</td>\n",
       "      <td>WP70</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9ULV4</th>\n",
       "      <td>COR1C_HUMAN Coronin-1C</td>\n",
       "      <td>COR1C</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9Y5X3</th>\n",
       "      <td>NX5_HUMAN Sorting nexin-5</td>\n",
       "      <td>NX5</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                                                   Name GeneID\n",
       "Protein IDs                                                                   \n",
       "A6NI72;P14598;A8MVU1  NCF1B_HUMAN Putative neutrophil cytosol factor...  NCF1B\n",
       "O00151                        PDLI1_HUMAN PDZ and LIM domain protein 1   PDLI1\n",
       "O15400                                            TX7_HUMAN Syntaxin-7     TX7\n",
       "O60262                GBG7_HUMAN Guanine nucleotide-binding protein ...   GBG7\n",
       "O60493                                       NX3_HUMAN Sorting nexin-3     NX3\n",
       "P00352                             AL1A1_HUMAN Retinal dehydrogenase 1   AL1A1\n",
       "P01871;P0DOX6              IGHM_HUMAN Immunoglobulin heavy constant mu    IGHM\n",
       "P01903;P01906;P52735  DRA_HUMAN HLA class II histocompatibility anti...    DRA\n",
       "P01911                2B1F_HUMAN HLA class II histocompatibility ant...   2B1F\n",
       "P04233                HG2A_HUMAN HLA class II histocompatibility ant...   HG2A\n",
       "P04440                DPB1_HUMAN HLA class II histocompatibility ant...   DPB1\n",
       "P04839                        CY24B_HUMAN Cytochrome b-245 heavy chain   CY24B\n",
       "P09917                          LOX5_HUMAN Arachidonate 5-lipoxygenase    LOX5\n",
       "P11586                C1TC_HUMAN C-1-tetrahydrofolate synthase, cyto...   C1TC\n",
       "P11836                            CD20_HUMAN B-lymphocyte antigen CD20    CD20\n",
       "P13284                GILT_HUMAN Gamma-interferon-inducible lysosoma...   GILT\n",
       "P13762                DRB4_HUMAN HLA class II histocompatibility ant...   DRB4\n",
       "P16885                PLCG2_HUMAN 1-phosphatidylinositol 4,5-bisphos...  PLCG2\n",
       "P27707                                  DCK_HUMAN Deoxycytidine kinase     DCK\n",
       "P28072                       PSB6_HUMAN Proteasome subunit beta type-6    PSB6\n",
       "P34897                GLYM_HUMAN Serine hydroxymethyltransferase, mi...   GLYM\n",
       "P49006                                MRP_HUMAN MARCKS-related protein     MRP\n",
       "P55327                                   TPD52_HUMAN Tumor protein D52   TPD52\n",
       "P56279                 TCL1A_HUMAN T-cell leukemia/lymphoma protein 1A   TCL1A\n",
       "P79483                DRB3_HUMAN HLA class II histocompatibility ant...   DRB3\n",
       "P80723                        BASP1_HUMAN Brain acid soluble protein 1   BASP1\n",
       "Q06187                           BTK_HUMAN Tyrosine-protein kinase BTK     BTK\n",
       "Q08209;P48454;P16298  PP2BA_HUMAN Serine/threonine-protein phosphata...  PP2BA\n",
       "Q14011                  CIRBP_HUMAN Cold-inducible RNA-binding protein   CIRBP\n",
       "Q6JBY9                            CPZIP_HUMAN CapZ-interacting protein   CPZIP\n",
       "Q6NUK1                CMC1_HUMAN Calcium-binding mitochondrial carri...   CMC1\n",
       "Q8NBS9                TXND5_HUMAN Thioredoxin domain-containing prot...  TXND5\n",
       "Q96D96                    HVCN1_HUMAN Voltage-gated hydrogen channel 1   HVCN1\n",
       "Q9H8H3                   MET7A_HUMAN Methyltransferase-like protein 7A   MET7A\n",
       "Q9NY33                               DPP3_HUMAN Dipeptidyl peptidase 3    DPP3\n",
       "Q9UH65                         WP70_HUMAN Switch-associated protein 70    WP70\n",
       "Q9ULV4                                          COR1C_HUMAN Coronin-1C   COR1C\n",
       "Q9Y5X3                                       NX5_HUMAN Sorting nexin-5     NX5"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "B_cell_only_geneset[['Name', 'GeneID']]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Proteins higher or only in B cells by GeneID\n",
    "<code>B_cell_set</code> is defined here as both the above lists of names, those higher or only in B cells."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [],
   "source": [
    "B_cell_set = list(higher_in_B_cells.GeneID)+list(B_cell_only_geneset.GeneID)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Proteins Higher in T Cells:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>pvalue</th>\n",
       "      <th>log2(B)-log2(T)</th>\n",
       "      <th>Name</th>\n",
       "      <th>GeneID</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Protein IDs</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>P26447</th>\n",
       "      <td>0.000140</td>\n",
       "      <td>-3.987732</td>\n",
       "      <td>10A4_HUMAN Protein S100-A4</td>\n",
       "      <td>10A4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q09666</th>\n",
       "      <td>0.004426</td>\n",
       "      <td>-2.451097</td>\n",
       "      <td>AHNK_HUMAN Neuroblast differentiation-associat...</td>\n",
       "      <td>AHNK</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9UI08</th>\n",
       "      <td>0.000003</td>\n",
       "      <td>-1.994234</td>\n",
       "      <td>EVL_HUMAN Ena/VASP-like protein</td>\n",
       "      <td>EVL</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P07195;Q6ZMR3;P07864</th>\n",
       "      <td>0.000193</td>\n",
       "      <td>-1.752858</td>\n",
       "      <td>LDHB_HUMAN L-lactate dehydrogenase B chain ;LD...</td>\n",
       "      <td>LDHB</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P08670;P41219;P17661;Q16352;P07197;P07196</th>\n",
       "      <td>0.000023</td>\n",
       "      <td>-1.658728</td>\n",
       "      <td>VIME_HUMAN Vimentin ;PERI_HUMAN Peripherin ;DE...</td>\n",
       "      <td>VIME</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O75347</th>\n",
       "      <td>0.000920</td>\n",
       "      <td>-1.605149</td>\n",
       "      <td>TBCA_HUMAN Tubulin-specific chaperone A</td>\n",
       "      <td>TBCA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O75083</th>\n",
       "      <td>0.000002</td>\n",
       "      <td>-1.590065</td>\n",
       "      <td>WDR1_HUMAN WD repeat-containing protein 1</td>\n",
       "      <td>WDR1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P21333;Q14315</th>\n",
       "      <td>0.000232</td>\n",
       "      <td>-1.520571</td>\n",
       "      <td>FLNA_HUMAN Filamin-A ;FLNC_HUMAN Filamin-C</td>\n",
       "      <td>FLNA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P60660;P14649</th>\n",
       "      <td>0.002014</td>\n",
       "      <td>-1.357658</td>\n",
       "      <td>MYL6_HUMAN Myosin light polypeptide 6 ;MYL6B_H...</td>\n",
       "      <td>MYL6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9UGI8</th>\n",
       "      <td>0.005253</td>\n",
       "      <td>-1.319492</td>\n",
       "      <td>TES_HUMAN Testin</td>\n",
       "      <td>TES</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P50995</th>\n",
       "      <td>0.001022</td>\n",
       "      <td>-1.276039</td>\n",
       "      <td>ANX11_HUMAN Annexin A11</td>\n",
       "      <td>ANX11</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q12906;Q96SI9</th>\n",
       "      <td>0.000723</td>\n",
       "      <td>-1.264212</td>\n",
       "      <td>ILF3_HUMAN Interleukin enhancer-binding factor...</td>\n",
       "      <td>ILF3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q14005</th>\n",
       "      <td>0.001427</td>\n",
       "      <td>-1.263467</td>\n",
       "      <td>IL16_HUMAN Pro-interleukin-16</td>\n",
       "      <td>IL16</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q14847</th>\n",
       "      <td>0.001216</td>\n",
       "      <td>-1.242621</td>\n",
       "      <td>LASP1_HUMAN LIM and SH3 domain protein 1</td>\n",
       "      <td>LASP1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q92522</th>\n",
       "      <td>0.000271</td>\n",
       "      <td>-1.206540</td>\n",
       "      <td>H1X_HUMAN Histone H1x</td>\n",
       "      <td>H1X</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P35579;P35580;REV__Q9UKV3;P12883;P13535;Q9UKX3;Q9Y623;P13533;P12882;P11055;Q9UKX2;A7E2Y1</th>\n",
       "      <td>0.000345</td>\n",
       "      <td>-1.155930</td>\n",
       "      <td>MYH9_HUMAN Myosin-9 ;MYH10_HUMAN Myosin-10 ;MY...</td>\n",
       "      <td>MYH9</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P62917</th>\n",
       "      <td>0.009131</td>\n",
       "      <td>-1.102065</td>\n",
       "      <td>RL8_HUMAN 60S ribosomal protein L8</td>\n",
       "      <td>RL8</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O15144</th>\n",
       "      <td>0.007176</td>\n",
       "      <td>-1.074163</td>\n",
       "      <td>ARPC2_HUMAN Actin-related protein 2/3 complex ...</td>\n",
       "      <td>ARPC2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P63000</th>\n",
       "      <td>0.009289</td>\n",
       "      <td>-1.028629</td>\n",
       "      <td>RAC1_HUMAN Ras-related C3 botulinum toxin subs...</td>\n",
       "      <td>RAC1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9UHD8</th>\n",
       "      <td>0.003477</td>\n",
       "      <td>-0.989198</td>\n",
       "      <td>EPT9_HUMAN Septin-9</td>\n",
       "      <td>EPT9</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P07741</th>\n",
       "      <td>0.000368</td>\n",
       "      <td>-0.921665</td>\n",
       "      <td>APT_HUMAN Adenine phosphoribosyltransferase</td>\n",
       "      <td>APT</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q16778;P33778;P23527;P06899;Q6DRA6;Q6DN03</th>\n",
       "      <td>0.000637</td>\n",
       "      <td>-0.920235</td>\n",
       "      <td>H2B2E_HUMAN Histone H2B type 2-E ;H2B1B_HUMAN ...</td>\n",
       "      <td>H2B2E</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P31946</th>\n",
       "      <td>0.005321</td>\n",
       "      <td>-0.849186</td>\n",
       "      <td>1433B_HUMAN 14-3-3 protein beta/alpha</td>\n",
       "      <td>1433B</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q15027</th>\n",
       "      <td>0.001002</td>\n",
       "      <td>-0.753567</td>\n",
       "      <td>ACAP1_HUMAN Arf-GAP with coiled-coil, ANK repe...</td>\n",
       "      <td>ACAP1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P26038</th>\n",
       "      <td>0.000987</td>\n",
       "      <td>-0.713509</td>\n",
       "      <td>MOES_HUMAN Moesin</td>\n",
       "      <td>MOES</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q01518;P40123</th>\n",
       "      <td>0.000613</td>\n",
       "      <td>-0.708052</td>\n",
       "      <td>CAP1_HUMAN Adenylyl cyclase-associated protein...</td>\n",
       "      <td>CAP1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P67936</th>\n",
       "      <td>0.000497</td>\n",
       "      <td>-0.625443</td>\n",
       "      <td>TPM4_HUMAN Tropomyosin alpha-4 chain</td>\n",
       "      <td>TPM4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P13073</th>\n",
       "      <td>0.009222</td>\n",
       "      <td>-0.605568</td>\n",
       "      <td>COX41_HUMAN Cytochrome c oxidase subunit 4 iso...</td>\n",
       "      <td>COX41</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P31146</th>\n",
       "      <td>0.009155</td>\n",
       "      <td>-0.583292</td>\n",
       "      <td>COR1A_HUMAN Coronin-1A</td>\n",
       "      <td>COR1A</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P31947</th>\n",
       "      <td>0.002352</td>\n",
       "      <td>-0.508935</td>\n",
       "      <td>1433S_HUMAN 14-3-3 protein sigma</td>\n",
       "      <td>1433S</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                                      pvalue  log2(B)-log2(T)  \\\n",
       "Protein IDs                                                                     \n",
       "P26447                                              0.000140        -3.987732   \n",
       "Q09666                                              0.004426        -2.451097   \n",
       "Q9UI08                                              0.000003        -1.994234   \n",
       "P07195;Q6ZMR3;P07864                                0.000193        -1.752858   \n",
       "P08670;P41219;P17661;Q16352;P07197;P07196           0.000023        -1.658728   \n",
       "O75347                                              0.000920        -1.605149   \n",
       "O75083                                              0.000002        -1.590065   \n",
       "P21333;Q14315                                       0.000232        -1.520571   \n",
       "P60660;P14649                                       0.002014        -1.357658   \n",
       "Q9UGI8                                              0.005253        -1.319492   \n",
       "P50995                                              0.001022        -1.276039   \n",
       "Q12906;Q96SI9                                       0.000723        -1.264212   \n",
       "Q14005                                              0.001427        -1.263467   \n",
       "Q14847                                              0.001216        -1.242621   \n",
       "Q92522                                              0.000271        -1.206540   \n",
       "P35579;P35580;REV__Q9UKV3;P12883;P13535;Q9UKX3;...  0.000345        -1.155930   \n",
       "P62917                                              0.009131        -1.102065   \n",
       "O15144                                              0.007176        -1.074163   \n",
       "P63000                                              0.009289        -1.028629   \n",
       "Q9UHD8                                              0.003477        -0.989198   \n",
       "P07741                                              0.000368        -0.921665   \n",
       "Q16778;P33778;P23527;P06899;Q6DRA6;Q6DN03           0.000637        -0.920235   \n",
       "P31946                                              0.005321        -0.849186   \n",
       "Q15027                                              0.001002        -0.753567   \n",
       "P26038                                              0.000987        -0.713509   \n",
       "Q01518;P40123                                       0.000613        -0.708052   \n",
       "P67936                                              0.000497        -0.625443   \n",
       "P13073                                              0.009222        -0.605568   \n",
       "P31146                                              0.009155        -0.583292   \n",
       "P31947                                              0.002352        -0.508935   \n",
       "\n",
       "                                                                                                 Name  \\\n",
       "Protein IDs                                                                                             \n",
       "P26447                                                                    10A4_HUMAN Protein S100-A4    \n",
       "Q09666                                              AHNK_HUMAN Neuroblast differentiation-associat...   \n",
       "Q9UI08                                                               EVL_HUMAN Ena/VASP-like protein    \n",
       "P07195;Q6ZMR3;P07864                                LDHB_HUMAN L-lactate dehydrogenase B chain ;LD...   \n",
       "P08670;P41219;P17661;Q16352;P07197;P07196           VIME_HUMAN Vimentin ;PERI_HUMAN Peripherin ;DE...   \n",
       "O75347                                                       TBCA_HUMAN Tubulin-specific chaperone A    \n",
       "O75083                                                     WDR1_HUMAN WD repeat-containing protein 1    \n",
       "P21333;Q14315                                             FLNA_HUMAN Filamin-A ;FLNC_HUMAN Filamin-C    \n",
       "P60660;P14649                                       MYL6_HUMAN Myosin light polypeptide 6 ;MYL6B_H...   \n",
       "Q9UGI8                                                                              TES_HUMAN Testin    \n",
       "P50995                                                                       ANX11_HUMAN Annexin A11    \n",
       "Q12906;Q96SI9                                       ILF3_HUMAN Interleukin enhancer-binding factor...   \n",
       "Q14005                                                                 IL16_HUMAN Pro-interleukin-16    \n",
       "Q14847                                                      LASP1_HUMAN LIM and SH3 domain protein 1    \n",
       "Q92522                                                                         H1X_HUMAN Histone H1x    \n",
       "P35579;P35580;REV__Q9UKV3;P12883;P13535;Q9UKX3;...  MYH9_HUMAN Myosin-9 ;MYH10_HUMAN Myosin-10 ;MY...   \n",
       "P62917                                                            RL8_HUMAN 60S ribosomal protein L8    \n",
       "O15144                                              ARPC2_HUMAN Actin-related protein 2/3 complex ...   \n",
       "P63000                                              RAC1_HUMAN Ras-related C3 botulinum toxin subs...   \n",
       "Q9UHD8                                                                           EPT9_HUMAN Septin-9    \n",
       "P07741                                                   APT_HUMAN Adenine phosphoribosyltransferase    \n",
       "Q16778;P33778;P23527;P06899;Q6DRA6;Q6DN03           H2B2E_HUMAN Histone H2B type 2-E ;H2B1B_HUMAN ...   \n",
       "P31946                                                         1433B_HUMAN 14-3-3 protein beta/alpha    \n",
       "Q15027                                              ACAP1_HUMAN Arf-GAP with coiled-coil, ANK repe...   \n",
       "P26038                                                                             MOES_HUMAN Moesin    \n",
       "Q01518;P40123                                       CAP1_HUMAN Adenylyl cyclase-associated protein...   \n",
       "P67936                                                          TPM4_HUMAN Tropomyosin alpha-4 chain    \n",
       "P13073                                              COX41_HUMAN Cytochrome c oxidase subunit 4 iso...   \n",
       "P31146                                                                        COR1A_HUMAN Coronin-1A    \n",
       "P31947                                                              1433S_HUMAN 14-3-3 protein sigma    \n",
       "\n",
       "                                                   GeneID  \n",
       "Protein IDs                                                \n",
       "P26447                                               10A4  \n",
       "Q09666                                               AHNK  \n",
       "Q9UI08                                                EVL  \n",
       "P07195;Q6ZMR3;P07864                                 LDHB  \n",
       "P08670;P41219;P17661;Q16352;P07197;P07196            VIME  \n",
       "O75347                                               TBCA  \n",
       "O75083                                               WDR1  \n",
       "P21333;Q14315                                        FLNA  \n",
       "P60660;P14649                                        MYL6  \n",
       "Q9UGI8                                                TES  \n",
       "P50995                                              ANX11  \n",
       "Q12906;Q96SI9                                        ILF3  \n",
       "Q14005                                               IL16  \n",
       "Q14847                                              LASP1  \n",
       "Q92522                                                H1X  \n",
       "P35579;P35580;REV__Q9UKV3;P12883;P13535;Q9UKX3;...   MYH9  \n",
       "P62917                                                RL8  \n",
       "O15144                                              ARPC2  \n",
       "P63000                                               RAC1  \n",
       "Q9UHD8                                               EPT9  \n",
       "P07741                                                APT  \n",
       "Q16778;P33778;P23527;P06899;Q6DRA6;Q6DN03           H2B2E  \n",
       "P31946                                              1433B  \n",
       "Q15027                                              ACAP1  \n",
       "P26038                                               MOES  \n",
       "Q01518;P40123                                        CAP1  \n",
       "P67936                                               TPM4  \n",
       "P13073                                              COX41  \n",
       "P31146                                              COR1A  \n",
       "P31947                                              1433S  "
      ]
     },
     "execution_count": 22,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "higher_in_T_cells.sort_values('log2(B)-log2(T)',ascending=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Proteins Only in T Cells:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Name</th>\n",
       "      <th>GeneID</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Protein IDs</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>P47813;O14602</th>\n",
       "      <td>IF1AX_HUMAN Eukaryotic translation initiation ...</td>\n",
       "      <td>IF1AX</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O15117</th>\n",
       "      <td>FYB1_HUMAN FYN-binding protein 1</td>\n",
       "      <td>FYB1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O60610</th>\n",
       "      <td>DIAP1_HUMAN Protein diaphanous homolog 1</td>\n",
       "      <td>DIAP1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>O94804;Q9H2G2</th>\n",
       "      <td>TK10_HUMAN Serine/threonine-protein kinase 10 ...</td>\n",
       "      <td>TK10</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P02545</th>\n",
       "      <td>LMNA_HUMAN Prelamin-A/C</td>\n",
       "      <td>LMNA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P04083</th>\n",
       "      <td>ANXA1_HUMAN Annexin A1</td>\n",
       "      <td>ANXA1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P05107</th>\n",
       "      <td>ITB2_HUMAN Integrin beta-2</td>\n",
       "      <td>ITB2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P06127</th>\n",
       "      <td>CD5_HUMAN T-cell surface glycoprotein CD5</td>\n",
       "      <td>CD5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P06239</th>\n",
       "      <td>LCK_HUMAN Tyrosine-protein kinase Lck</td>\n",
       "      <td>LCK</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P06729</th>\n",
       "      <td>CD2_HUMAN T-cell surface antigen CD2</td>\n",
       "      <td>CD2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P07305</th>\n",
       "      <td>H10_HUMAN Histone H1.0</td>\n",
       "      <td>H10</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P07766</th>\n",
       "      <td>CD3E_HUMAN T-cell surface glycoprotein CD3 eps...</td>\n",
       "      <td>CD3E</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P09104</th>\n",
       "      <td>ENOG_HUMAN Gamma-enolase</td>\n",
       "      <td>ENOG</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P09326</th>\n",
       "      <td>CD48_HUMAN CD48 antigen</td>\n",
       "      <td>CD48</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P17655</th>\n",
       "      <td>CAN2_HUMAN Calpain-2 catalytic subunit</td>\n",
       "      <td>CAN2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P23743</th>\n",
       "      <td>DGKA_HUMAN Diacylglycerol kinase alpha</td>\n",
       "      <td>DGKA</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P26885</th>\n",
       "      <td>FKBP2_HUMAN Peptidyl-prolyl cis-trans isomeras...</td>\n",
       "      <td>FKBP2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P30203</th>\n",
       "      <td>CD6_HUMAN T-cell differentiation antigen CD6</td>\n",
       "      <td>CD6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P32455</th>\n",
       "      <td>GBP1_HUMAN Guanylate-binding protein 1</td>\n",
       "      <td>GBP1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P43403</th>\n",
       "      <td>ZAP70_HUMAN Tyrosine-protein kinase ZAP-70</td>\n",
       "      <td>ZAP70</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P49588</th>\n",
       "      <td>YAC_HUMAN Alanine--tRNA ligase, cytoplasmic</td>\n",
       "      <td>YAC</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P49863</th>\n",
       "      <td>GRAK_HUMAN Granzyme K</td>\n",
       "      <td>GRAK</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P51812</th>\n",
       "      <td>KS6A3_HUMAN Ribosomal protein S6 kinase alpha-3</td>\n",
       "      <td>KS6A3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P53634</th>\n",
       "      <td>CATC_HUMAN Dipeptidyl peptidase 1</td>\n",
       "      <td>CATC</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P57737</th>\n",
       "      <td>CORO7_HUMAN Coronin-7</td>\n",
       "      <td>CORO7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P59768</th>\n",
       "      <td>GBG2_HUMAN Guanine nucleotide-binding protein ...</td>\n",
       "      <td>GBG2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P60981</th>\n",
       "      <td>DEST_HUMAN Destrin</td>\n",
       "      <td>DEST</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P61221</th>\n",
       "      <td>ABCE1_HUMAN ATP-binding cassette sub-family E ...</td>\n",
       "      <td>ABCE1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>P84095</th>\n",
       "      <td>RHOG_HUMAN Rho-related GTP-binding protein RhoG</td>\n",
       "      <td>RHOG</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q13094</th>\n",
       "      <td>LCP2_HUMAN Lymphocyte cytosolic protein 2</td>\n",
       "      <td>LCP2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q14764</th>\n",
       "      <td>MVP_HUMAN Major vault protein</td>\n",
       "      <td>MVP</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q15942</th>\n",
       "      <td>ZYX_HUMAN Zyxin</td>\n",
       "      <td>ZYX</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q16401</th>\n",
       "      <td>PSMD5_HUMAN 26S proteasome non-ATPase regulato...</td>\n",
       "      <td>PSMD5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q5SXH7</th>\n",
       "      <td>PKHS1_HUMAN Pleckstrin homology domain-contain...</td>\n",
       "      <td>PKHS1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q7Z5R6</th>\n",
       "      <td>AB1IP_HUMAN Amyloid beta A4 precursor protein-...</td>\n",
       "      <td>AB1IP</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q8IUI8</th>\n",
       "      <td>CRLF3_HUMAN Cytokine receptor-like factor 3</td>\n",
       "      <td>CRLF3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q8NHV1</th>\n",
       "      <td>GIMA7_HUMAN GTPase IMAP family member 7</td>\n",
       "      <td>GIMA7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q8WWP7</th>\n",
       "      <td>GIMA1_HUMAN GTPase IMAP family member 1</td>\n",
       "      <td>GIMA1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q96B97</th>\n",
       "      <td>H3K1_HUMAN SH3 domain-containing kinase-bindin...</td>\n",
       "      <td>H3K1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9BUH6</th>\n",
       "      <td>PAXX_HUMAN Protein PAXX</td>\n",
       "      <td>PAXX</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9H1E5</th>\n",
       "      <td>TMX4_HUMAN Thioredoxin-related transmembrane p...</td>\n",
       "      <td>TMX4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9H4E7</th>\n",
       "      <td>DEFI6_HUMAN Differentially expressed in FDCP 6...</td>\n",
       "      <td>DEFI6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9H9B4</th>\n",
       "      <td>FXN1_HUMAN Sideroflexin-1</td>\n",
       "      <td>FXN1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9NUV9;Q6P9H5</th>\n",
       "      <td>GIMA4_HUMAN GTPase IMAP family member 4 ;GIMA6...</td>\n",
       "      <td>GIMA4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9NXG2</th>\n",
       "      <td>THUM1_HUMAN THUMP domain-containing protein 1</td>\n",
       "      <td>THUM1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9P270</th>\n",
       "      <td>LAI2_HUMAN SLAIN motif-containing protein 2</td>\n",
       "      <td>LAI2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9UBW5</th>\n",
       "      <td>BIN2_HUMAN Bridging integrator 2</td>\n",
       "      <td>BIN2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9UNZ5</th>\n",
       "      <td>L10K_HUMAN Leydig cell tumor 10 kDa protein ho...</td>\n",
       "      <td>L10K</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Q9Y3Z3</th>\n",
       "      <td>AMH1_HUMAN Deoxynucleoside triphosphate tripho...</td>\n",
       "      <td>AMH1</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                                            Name GeneID\n",
       "Protein IDs                                                            \n",
       "P47813;O14602  IF1AX_HUMAN Eukaryotic translation initiation ...  IF1AX\n",
       "O15117                         FYB1_HUMAN FYN-binding protein 1    FYB1\n",
       "O60610                 DIAP1_HUMAN Protein diaphanous homolog 1   DIAP1\n",
       "O94804;Q9H2G2  TK10_HUMAN Serine/threonine-protein kinase 10 ...   TK10\n",
       "P02545                                  LMNA_HUMAN Prelamin-A/C    LMNA\n",
       "P04083                                   ANXA1_HUMAN Annexin A1   ANXA1\n",
       "P05107                               ITB2_HUMAN Integrin beta-2    ITB2\n",
       "P06127                CD5_HUMAN T-cell surface glycoprotein CD5     CD5\n",
       "P06239                    LCK_HUMAN Tyrosine-protein kinase Lck     LCK\n",
       "P06729                     CD2_HUMAN T-cell surface antigen CD2     CD2\n",
       "P07305                                   H10_HUMAN Histone H1.0     H10\n",
       "P07766         CD3E_HUMAN T-cell surface glycoprotein CD3 eps...   CD3E\n",
       "P09104                                 ENOG_HUMAN Gamma-enolase    ENOG\n",
       "P09326                                  CD48_HUMAN CD48 antigen    CD48\n",
       "P17655                   CAN2_HUMAN Calpain-2 catalytic subunit    CAN2\n",
       "P23743                   DGKA_HUMAN Diacylglycerol kinase alpha    DGKA\n",
       "P26885         FKBP2_HUMAN Peptidyl-prolyl cis-trans isomeras...  FKBP2\n",
       "P30203             CD6_HUMAN T-cell differentiation antigen CD6     CD6\n",
       "P32455                   GBP1_HUMAN Guanylate-binding protein 1    GBP1\n",
       "P43403               ZAP70_HUMAN Tyrosine-protein kinase ZAP-70   ZAP70\n",
       "P49588              YAC_HUMAN Alanine--tRNA ligase, cytoplasmic     YAC\n",
       "P49863                                    GRAK_HUMAN Granzyme K    GRAK\n",
       "P51812          KS6A3_HUMAN Ribosomal protein S6 kinase alpha-3   KS6A3\n",
       "P53634                        CATC_HUMAN Dipeptidyl peptidase 1    CATC\n",
       "P57737                                    CORO7_HUMAN Coronin-7   CORO7\n",
       "P59768         GBG2_HUMAN Guanine nucleotide-binding protein ...   GBG2\n",
       "P60981                                       DEST_HUMAN Destrin    DEST\n",
       "P61221         ABCE1_HUMAN ATP-binding cassette sub-family E ...  ABCE1\n",
       "P84095          RHOG_HUMAN Rho-related GTP-binding protein RhoG    RHOG\n",
       "Q13094                LCP2_HUMAN Lymphocyte cytosolic protein 2    LCP2\n",
       "Q14764                            MVP_HUMAN Major vault protein     MVP\n",
       "Q15942                                          ZYX_HUMAN Zyxin     ZYX\n",
       "Q16401         PSMD5_HUMAN 26S proteasome non-ATPase regulato...  PSMD5\n",
       "Q5SXH7         PKHS1_HUMAN Pleckstrin homology domain-contain...  PKHS1\n",
       "Q7Z5R6         AB1IP_HUMAN Amyloid beta A4 precursor protein-...  AB1IP\n",
       "Q8IUI8              CRLF3_HUMAN Cytokine receptor-like factor 3   CRLF3\n",
       "Q8NHV1                  GIMA7_HUMAN GTPase IMAP family member 7   GIMA7\n",
       "Q8WWP7                  GIMA1_HUMAN GTPase IMAP family member 1   GIMA1\n",
       "Q96B97         H3K1_HUMAN SH3 domain-containing kinase-bindin...   H3K1\n",
       "Q9BUH6                                  PAXX_HUMAN Protein PAXX    PAXX\n",
       "Q9H1E5         TMX4_HUMAN Thioredoxin-related transmembrane p...   TMX4\n",
       "Q9H4E7         DEFI6_HUMAN Differentially expressed in FDCP 6...  DEFI6\n",
       "Q9H9B4                                FXN1_HUMAN Sideroflexin-1    FXN1\n",
       "Q9NUV9;Q6P9H5  GIMA4_HUMAN GTPase IMAP family member 4 ;GIMA6...  GIMA4\n",
       "Q9NXG2            THUM1_HUMAN THUMP domain-containing protein 1   THUM1\n",
       "Q9P270              LAI2_HUMAN SLAIN motif-containing protein 2    LAI2\n",
       "Q9UBW5                         BIN2_HUMAN Bridging integrator 2    BIN2\n",
       "Q9UNZ5         L10K_HUMAN Leydig cell tumor 10 kDa protein ho...   L10K\n",
       "Q9Y3Z3         AMH1_HUMAN Deoxynucleoside triphosphate tripho...   AMH1"
      ]
     },
     "execution_count": 23,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "T_cell_only_geneset[['Name', 'GeneID']]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Proteins higher or only in T cells by GeneID\n",
    "<code>T_cell_set</code> is defined here as both the above lists of names, those higher or only in T cells."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [],
   "source": [
    "T_cell_set = list(higher_in_T_cells.GeneID)+list(T_cell_only_geneset.GeneID)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}