statistics_datasets.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Calculate stats for the full TA-ORF dataset (after QC)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Calculate wells/sites/cells per pathway"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "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>Metadata_Plate</th>\n",
       "      <th>Metadata_Well</th>\n",
       "      <th>Metadata_Site</th>\n",
       "      <th>DNA</th>\n",
       "      <th>RNA</th>\n",
       "      <th>ER</th>\n",
       "      <th>AGP</th>\n",
       "      <th>Mito</th>\n",
       "      <th>gene_name</th>\n",
       "      <th>pert_name</th>\n",
       "      <th>broad_sample</th>\n",
       "      <th>pert_name_replicate</th>\n",
       "      <th>val</th>\n",
       "      <th>Cells</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>41744</td>\n",
       "      <td>a01</td>\n",
       "      <td>1</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s1_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s1_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s1_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s1_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s1_...</td>\n",
       "      <td>EMPTY</td>\n",
       "      <td>EMPTY_</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>38</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>41744</td>\n",
       "      <td>a01</td>\n",
       "      <td>2</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s2_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s2_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s2_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s2_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s2_...</td>\n",
       "      <td>EMPTY</td>\n",
       "      <td>EMPTY_</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>34</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>41744</td>\n",
       "      <td>a01</td>\n",
       "      <td>3</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s3_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s3_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s3_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s3_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s3_...</td>\n",
       "      <td>EMPTY</td>\n",
       "      <td>EMPTY_</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>50</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>41744</td>\n",
       "      <td>a01</td>\n",
       "      <td>4</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s4_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s4_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s4_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s4_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s4_...</td>\n",
       "      <td>EMPTY</td>\n",
       "      <td>EMPTY_</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>68</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>41744</td>\n",
       "      <td>a01</td>\n",
       "      <td>5</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s5_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s5_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s5_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s5_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a01_s5_...</td>\n",
       "      <td>EMPTY</td>\n",
       "      <td>EMPTY_</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>27</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",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17216</th>\n",
       "      <td>41757</td>\n",
       "      <td>p24</td>\n",
       "      <td>5</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s5_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s5_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s5_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s5_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s5_...</td>\n",
       "      <td>EMPTY</td>\n",
       "      <td>EMPTY_</td>\n",
       "      <td>NaN</td>\n",
       "      <td>175</td>\n",
       "      <td>0</td>\n",
       "      <td>37</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17217</th>\n",
       "      <td>41757</td>\n",
       "      <td>p24</td>\n",
       "      <td>6</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s6_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s6_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s6_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s6_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s6_...</td>\n",
       "      <td>EMPTY</td>\n",
       "      <td>EMPTY_</td>\n",
       "      <td>NaN</td>\n",
       "      <td>175</td>\n",
       "      <td>0</td>\n",
       "      <td>44</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17218</th>\n",
       "      <td>41757</td>\n",
       "      <td>p24</td>\n",
       "      <td>7</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s7_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s7_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s7_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s7_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s7_...</td>\n",
       "      <td>EMPTY</td>\n",
       "      <td>EMPTY_</td>\n",
       "      <td>NaN</td>\n",
       "      <td>175</td>\n",
       "      <td>0</td>\n",
       "      <td>62</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17219</th>\n",
       "      <td>41757</td>\n",
       "      <td>p24</td>\n",
       "      <td>8</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s8_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s8_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s8_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s8_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s8_...</td>\n",
       "      <td>EMPTY</td>\n",
       "      <td>EMPTY_</td>\n",
       "      <td>NaN</td>\n",
       "      <td>175</td>\n",
       "      <td>0</td>\n",
       "      <td>53</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17220</th>\n",
       "      <td>41757</td>\n",
       "      <td>p24</td>\n",
       "      <td>9</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s9_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s9_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s9_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s9_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p24_s9_...</td>\n",
       "      <td>EMPTY</td>\n",
       "      <td>EMPTY_</td>\n",
       "      <td>NaN</td>\n",
       "      <td>175</td>\n",
       "      <td>0</td>\n",
       "      <td>34</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>17221 rows × 14 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "       Metadata_Plate Metadata_Well  Metadata_Site  \\\n",
       "0               41744           a01              1   \n",
       "1               41744           a01              2   \n",
       "2               41744           a01              3   \n",
       "3               41744           a01              4   \n",
       "4               41744           a01              5   \n",
       "...               ...           ...            ...   \n",
       "17216           41757           p24              5   \n",
       "17217           41757           p24              6   \n",
       "17218           41757           p24              7   \n",
       "17219           41757           p24              8   \n",
       "17220           41757           p24              9   \n",
       "\n",
       "                                                     DNA  \\\n",
       "0      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s1_...   \n",
       "1      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s2_...   \n",
       "2      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s3_...   \n",
       "3      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s4_...   \n",
       "4      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s5_...   \n",
       "...                                                  ...   \n",
       "17216  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s5_...   \n",
       "17217  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s6_...   \n",
       "17218  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s7_...   \n",
       "17219  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s8_...   \n",
       "17220  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s9_...   \n",
       "\n",
       "                                                     RNA  \\\n",
       "0      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s1_...   \n",
       "1      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s2_...   \n",
       "2      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s3_...   \n",
       "3      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s4_...   \n",
       "4      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s5_...   \n",
       "...                                                  ...   \n",
       "17216  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s5_...   \n",
       "17217  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s6_...   \n",
       "17218  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s7_...   \n",
       "17219  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s8_...   \n",
       "17220  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s9_...   \n",
       "\n",
       "                                                      ER  \\\n",
       "0      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s1_...   \n",
       "1      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s2_...   \n",
       "2      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s3_...   \n",
       "3      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s4_...   \n",
       "4      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s5_...   \n",
       "...                                                  ...   \n",
       "17216  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s5_...   \n",
       "17217  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s6_...   \n",
       "17218  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s7_...   \n",
       "17219  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s8_...   \n",
       "17220  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s9_...   \n",
       "\n",
       "                                                     AGP  \\\n",
       "0      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s1_...   \n",
       "1      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s2_...   \n",
       "2      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s3_...   \n",
       "3      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s4_...   \n",
       "4      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s5_...   \n",
       "...                                                  ...   \n",
       "17216  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s5_...   \n",
       "17217  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s6_...   \n",
       "17218  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s7_...   \n",
       "17219  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s8_...   \n",
       "17220  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s9_...   \n",
       "\n",
       "                                                    Mito gene_name pert_name  \\\n",
       "0      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s1_...     EMPTY    EMPTY_   \n",
       "1      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s2_...     EMPTY    EMPTY_   \n",
       "2      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s3_...     EMPTY    EMPTY_   \n",
       "3      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s4_...     EMPTY    EMPTY_   \n",
       "4      41744/taoe005-u2os-72h-cp-a-au00044859_a01_s5_...     EMPTY    EMPTY_   \n",
       "...                                                  ...       ...       ...   \n",
       "17216  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s5_...     EMPTY    EMPTY_   \n",
       "17217  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s6_...     EMPTY    EMPTY_   \n",
       "17218  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s7_...     EMPTY    EMPTY_   \n",
       "17219  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s8_...     EMPTY    EMPTY_   \n",
       "17220  41757/taoe005-u2os-72h-cp-a-au00044858_p24_s9_...     EMPTY    EMPTY_   \n",
       "\n",
       "      broad_sample  pert_name_replicate  val  Cells  \n",
       "0              NaN                    1    0     38  \n",
       "1              NaN                    1    0     34  \n",
       "2              NaN                    1    0     50  \n",
       "3              NaN                    1    0     68  \n",
       "4              NaN                    1    0     27  \n",
       "...            ...                  ...  ...    ...  \n",
       "17216          NaN                  175    0     37  \n",
       "17217          NaN                  175    0     44  \n",
       "17218          NaN                  175    0     62  \n",
       "17219          NaN                  175    0     53  \n",
       "17220          NaN                  175    0     34  \n",
       "\n",
       "[17221 rows x 14 columns]"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "taorf_meta = pd.read_csv(\"/raid/data/cellpainting/TAORF/inputs/metadata/index_new_taorf_minus6wells.csv\")\n",
    "cells_taorf = []\n",
    "for i in taorf_meta.index:\n",
    "    cells_taorf.append([\n",
    "        taorf_meta.loc[i, \"Metadata_Plate\"],\n",
    "        taorf_meta.loc[i, \"Metadata_Well\"],\n",
    "        taorf_meta.loc[i, \"Metadata_Site\"],\n",
    "        sum(1 for line in open('/raid/data/cellpainting/TAORF/inputs/locations/{}/{}-{}-Nuclei.csv'.format(taorf_meta.loc[i, \"Metadata_Plate\"], taorf_meta.loc[i, \"Metadata_Well\"], taorf_meta.loc[i, \"Metadata_Site\"]))) - 1\n",
    "    ])\n",
    "counts = pd.DataFrame(data=cells_taorf, columns=['Metadata_Plate', 'Metadata_Well', 'Metadata_Site', 'Cells'])\n",
    "taorf_meta = pd.merge(taorf_meta, counts, right_on=['Metadata_Plate', 'Metadata_Well', 'Metadata_Site'], left_on =['Metadata_Plate', 'Metadata_Well', 'Metadata_Site'])\n",
    "\n",
    "taorf_meta"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "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>Metadata_Plate</th>\n",
       "      <th>Metadata_Well</th>\n",
       "      <th>Metadata_Site</th>\n",
       "      <th>DNA</th>\n",
       "      <th>RNA</th>\n",
       "      <th>ER</th>\n",
       "      <th>AGP</th>\n",
       "      <th>Mito</th>\n",
       "      <th>gene_name</th>\n",
       "      <th>pert_name</th>\n",
       "      <th>broad_sample</th>\n",
       "      <th>pert_name_replicate</th>\n",
       "      <th>val</th>\n",
       "      <th>Cells</th>\n",
       "      <th>Var1</th>\n",
       "      <th>Metadata_moa.x</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>41744</td>\n",
       "      <td>a03</td>\n",
       "      <td>1</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s1_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s1_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s1_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s1_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s1_...</td>\n",
       "      <td>ATF4</td>\n",
       "      <td>ATF4_WT.1</td>\n",
       "      <td>ccsbBroad304_00117</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>40</td>\n",
       "      <td>ccsbBroad304_00117</td>\n",
       "      <td>ER Stress/UPR</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>41744</td>\n",
       "      <td>a03</td>\n",
       "      <td>2</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s2_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s2_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s2_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s2_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s2_...</td>\n",
       "      <td>ATF4</td>\n",
       "      <td>ATF4_WT.1</td>\n",
       "      <td>ccsbBroad304_00117</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>16</td>\n",
       "      <td>ccsbBroad304_00117</td>\n",
       "      <td>ER Stress/UPR</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>41744</td>\n",
       "      <td>a03</td>\n",
       "      <td>3</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s3_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s3_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s3_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s3_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s3_...</td>\n",
       "      <td>ATF4</td>\n",
       "      <td>ATF4_WT.1</td>\n",
       "      <td>ccsbBroad304_00117</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>31</td>\n",
       "      <td>ccsbBroad304_00117</td>\n",
       "      <td>ER Stress/UPR</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>41744</td>\n",
       "      <td>a03</td>\n",
       "      <td>4</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s4_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s4_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s4_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s4_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s4_...</td>\n",
       "      <td>ATF4</td>\n",
       "      <td>ATF4_WT.1</td>\n",
       "      <td>ccsbBroad304_00117</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>16</td>\n",
       "      <td>ccsbBroad304_00117</td>\n",
       "      <td>ER Stress/UPR</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>41744</td>\n",
       "      <td>a03</td>\n",
       "      <td>5</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s5_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s5_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s5_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s5_...</td>\n",
       "      <td>41744/taoe005-u2os-72h-cp-a-au00044859_a03_s5_...</td>\n",
       "      <td>ATF4</td>\n",
       "      <td>ATF4_WT.1</td>\n",
       "      <td>ccsbBroad304_00117</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>10</td>\n",
       "      <td>ccsbBroad304_00117</td>\n",
       "      <td>ER Stress/UPR</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",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9246</th>\n",
       "      <td>41757</td>\n",
       "      <td>p16</td>\n",
       "      <td>5</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s5_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s5_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s5_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s5_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s5_...</td>\n",
       "      <td>ARNTL</td>\n",
       "      <td>ARNTL_WT</td>\n",
       "      <td>ccsbBroad304_00105</td>\n",
       "      <td>5</td>\n",
       "      <td>0</td>\n",
       "      <td>27</td>\n",
       "      <td>ccsbBroad304_00105</td>\n",
       "      <td>Circadian Rhythm</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9247</th>\n",
       "      <td>41757</td>\n",
       "      <td>p16</td>\n",
       "      <td>6</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s6_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s6_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s6_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s6_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s6_...</td>\n",
       "      <td>ARNTL</td>\n",
       "      <td>ARNTL_WT</td>\n",
       "      <td>ccsbBroad304_00105</td>\n",
       "      <td>5</td>\n",
       "      <td>0</td>\n",
       "      <td>72</td>\n",
       "      <td>ccsbBroad304_00105</td>\n",
       "      <td>Circadian Rhythm</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9248</th>\n",
       "      <td>41757</td>\n",
       "      <td>p16</td>\n",
       "      <td>7</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s7_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s7_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s7_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s7_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s7_...</td>\n",
       "      <td>ARNTL</td>\n",
       "      <td>ARNTL_WT</td>\n",
       "      <td>ccsbBroad304_00105</td>\n",
       "      <td>5</td>\n",
       "      <td>0</td>\n",
       "      <td>48</td>\n",
       "      <td>ccsbBroad304_00105</td>\n",
       "      <td>Circadian Rhythm</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9249</th>\n",
       "      <td>41757</td>\n",
       "      <td>p16</td>\n",
       "      <td>8</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s8_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s8_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s8_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s8_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s8_...</td>\n",
       "      <td>ARNTL</td>\n",
       "      <td>ARNTL_WT</td>\n",
       "      <td>ccsbBroad304_00105</td>\n",
       "      <td>5</td>\n",
       "      <td>0</td>\n",
       "      <td>31</td>\n",
       "      <td>ccsbBroad304_00105</td>\n",
       "      <td>Circadian Rhythm</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9250</th>\n",
       "      <td>41757</td>\n",
       "      <td>p16</td>\n",
       "      <td>9</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s9_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s9_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s9_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s9_...</td>\n",
       "      <td>41757/taoe005-u2os-72h-cp-a-au00044858_p16_s9_...</td>\n",
       "      <td>ARNTL</td>\n",
       "      <td>ARNTL_WT</td>\n",
       "      <td>ccsbBroad304_00105</td>\n",
       "      <td>5</td>\n",
       "      <td>0</td>\n",
       "      <td>78</td>\n",
       "      <td>ccsbBroad304_00105</td>\n",
       "      <td>Circadian Rhythm</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>9251 rows × 16 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      Metadata_Plate Metadata_Well  Metadata_Site  \\\n",
       "0              41744           a03              1   \n",
       "1              41744           a03              2   \n",
       "2              41744           a03              3   \n",
       "3              41744           a03              4   \n",
       "4              41744           a03              5   \n",
       "...              ...           ...            ...   \n",
       "9246           41757           p16              5   \n",
       "9247           41757           p16              6   \n",
       "9248           41757           p16              7   \n",
       "9249           41757           p16              8   \n",
       "9250           41757           p16              9   \n",
       "\n",
       "                                                    DNA  \\\n",
       "0     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s1_...   \n",
       "1     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s2_...   \n",
       "2     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s3_...   \n",
       "3     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s4_...   \n",
       "4     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s5_...   \n",
       "...                                                 ...   \n",
       "9246  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s5_...   \n",
       "9247  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s6_...   \n",
       "9248  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s7_...   \n",
       "9249  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s8_...   \n",
       "9250  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s9_...   \n",
       "\n",
       "                                                    RNA  \\\n",
       "0     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s1_...   \n",
       "1     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s2_...   \n",
       "2     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s3_...   \n",
       "3     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s4_...   \n",
       "4     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s5_...   \n",
       "...                                                 ...   \n",
       "9246  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s5_...   \n",
       "9247  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s6_...   \n",
       "9248  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s7_...   \n",
       "9249  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s8_...   \n",
       "9250  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s9_...   \n",
       "\n",
       "                                                     ER  \\\n",
       "0     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s1_...   \n",
       "1     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s2_...   \n",
       "2     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s3_...   \n",
       "3     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s4_...   \n",
       "4     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s5_...   \n",
       "...                                                 ...   \n",
       "9246  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s5_...   \n",
       "9247  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s6_...   \n",
       "9248  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s7_...   \n",
       "9249  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s8_...   \n",
       "9250  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s9_...   \n",
       "\n",
       "                                                    AGP  \\\n",
       "0     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s1_...   \n",
       "1     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s2_...   \n",
       "2     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s3_...   \n",
       "3     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s4_...   \n",
       "4     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s5_...   \n",
       "...                                                 ...   \n",
       "9246  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s5_...   \n",
       "9247  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s6_...   \n",
       "9248  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s7_...   \n",
       "9249  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s8_...   \n",
       "9250  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s9_...   \n",
       "\n",
       "                                                   Mito gene_name  pert_name  \\\n",
       "0     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s1_...      ATF4  ATF4_WT.1   \n",
       "1     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s2_...      ATF4  ATF4_WT.1   \n",
       "2     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s3_...      ATF4  ATF4_WT.1   \n",
       "3     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s4_...      ATF4  ATF4_WT.1   \n",
       "4     41744/taoe005-u2os-72h-cp-a-au00044859_a03_s5_...      ATF4  ATF4_WT.1   \n",
       "...                                                 ...       ...        ...   \n",
       "9246  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s5_...     ARNTL   ARNTL_WT   \n",
       "9247  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s6_...     ARNTL   ARNTL_WT   \n",
       "9248  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s7_...     ARNTL   ARNTL_WT   \n",
       "9249  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s8_...     ARNTL   ARNTL_WT   \n",
       "9250  41757/taoe005-u2os-72h-cp-a-au00044858_p16_s9_...     ARNTL   ARNTL_WT   \n",
       "\n",
       "            broad_sample  pert_name_replicate  val  Cells                Var1  \\\n",
       "0     ccsbBroad304_00117                    1    0     40  ccsbBroad304_00117   \n",
       "1     ccsbBroad304_00117                    1    0     16  ccsbBroad304_00117   \n",
       "2     ccsbBroad304_00117                    1    0     31  ccsbBroad304_00117   \n",
       "3     ccsbBroad304_00117                    1    0     16  ccsbBroad304_00117   \n",
       "4     ccsbBroad304_00117                    1    0     10  ccsbBroad304_00117   \n",
       "...                  ...                  ...  ...    ...                 ...   \n",
       "9246  ccsbBroad304_00105                    5    0     27  ccsbBroad304_00105   \n",
       "9247  ccsbBroad304_00105                    5    0     72  ccsbBroad304_00105   \n",
       "9248  ccsbBroad304_00105                    5    0     48  ccsbBroad304_00105   \n",
       "9249  ccsbBroad304_00105                    5    0     31  ccsbBroad304_00105   \n",
       "9250  ccsbBroad304_00105                    5    0     78  ccsbBroad304_00105   \n",
       "\n",
       "        Metadata_moa.x  \n",
       "0        ER Stress/UPR  \n",
       "1        ER Stress/UPR  \n",
       "2        ER Stress/UPR  \n",
       "3        ER Stress/UPR  \n",
       "4        ER Stress/UPR  \n",
       "...                ...  \n",
       "9246  Circadian Rhythm  \n",
       "9247  Circadian Rhythm  \n",
       "9248  Circadian Rhythm  \n",
       "9249  Circadian Rhythm  \n",
       "9250  Circadian Rhythm  \n",
       "\n",
       "[9251 rows x 16 columns]"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "taorf_moa_matches = pd.read_csv(\"/raid/data/cellpainting/TAORF/nikita_experiments/TAORF_MOA_MATCHES.csv\")\n",
    "taorf_meta_moa = pd.merge(taorf_meta, taorf_moa_matches[['Var1', 'Metadata_moa.x']], right_on=['Var1'], left_on =['broad_sample'])\n",
    "taorf_meta_moa"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Calculate perturbations per pathway"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(array([2., 2., 0., 6., 0., 4., 0., 3., 0., 1., 0., 2., 1., 0., 1., 0., 1.,\n",
       "        0., 0., 0., 1., 0., 2., 0., 0., 0., 0., 0., 0., 0., 0., 1., 0., 1.,\n",
       "        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 1.]),\n",
       " array([ 1.  ,  1.54,  2.08,  2.62,  3.16,  3.7 ,  4.24,  4.78,  5.32,\n",
       "         5.86,  6.4 ,  6.94,  7.48,  8.02,  8.56,  9.1 ,  9.64, 10.18,\n",
       "        10.72, 11.26, 11.8 , 12.34, 12.88, 13.42, 13.96, 14.5 , 15.04,\n",
       "        15.58, 16.12, 16.66, 17.2 , 17.74, 18.28, 18.82, 19.36, 19.9 ,\n",
       "        20.44, 20.98, 21.52, 22.06, 22.6 , 23.14, 23.68, 24.22, 24.76,\n",
       "        25.3 , 25.84, 26.38, 26.92, 27.46, 28.  ]),\n",
       " <a list of 50 Patch objects>)"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWoAAAD4CAYAAADFAawfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAALXklEQVR4nO3cb4hld33H8c/XbKRtDNWSIYQk07F/CIigCUNKMQSbokS31BZEDFS0tEwfmBKh0G77pPZBYSmt2AcibGvalEZF1LRiqFUwYoU27W6aav5oK2HFhJhNCMXEBw3Gbx/M3WQy3MncZOfM/Obu6wXL3nvuyb3fM4d95+y552x1dwAY1ysOegAAXpxQAwxOqAEGJ9QAgxNqgMEdmeJNL7nkkl5bW5virQGW0qlTp57o7pV5r00S6rW1tZw8eXKKtwZYSlX1nZ1ec+oDYHBCDTA4oQYYnFADDE6oAQYn1ACDWyjUVfXqqvp0VX2zqh6sql+cejAANi16HfVfJvlCd7+zql6Z5CcmnAmALXYNdVX9ZJLrk7wvSbr7mSTPTDsWAGctckT92iSPJ/mbqnpDklNJbunuH2xdqao2kmwkyerq6l7P+ZKtHbtz7vLTx4/u8yQA52aRc9RHklyT5KPdfXWSHyQ5tn2l7j7R3evdvb6yMvd2dQBehkVC/XCSh7v77tnzT2cz3ADsg11D3d3fS/LdqrpqtuiXkzww6VQAPGfRqz5+N8ntsys+Hkrym9ONBMBWC4W6u+9Nsj7xLADM4c5EgMEJNcDghBpgcEINMDihBhicUAMMTqgBBifUAIMTaoDBCTXA4IQaYHBCDTA4oQYYnFADDE6oAQYn1ACDE2qAwQk1wOCEGmBwQg0wOKEGGJxQAwxOqAEGJ9QAgxNqgMEdWWSlqjqd5Kkkzyb5YXevTzkUAM9bKNQzv9TdT0w2CQBzOfUBMLhFQ91JvlhVp6pqY94KVbVRVSer6uTjjz++dxMCnOcWDfV13X1NkrcleX9VXb99he4+0d3r3b2+srKyp0MCnM8WCnV3PzL7/UySO5JcO+VQADxv11BX1UVVdfHZx0nemuS+qQcDYNMiV31cmuSOqjq7/se7+wuTTgXAc3YNdXc/lOQN+zALAHO4PA9gcEINMDihBhicUAMMTqgBBifUAIMTaoDBCTXA4IQaYHBCDTA4oQYYnFADDE6oAQYn1ACDE2qAwQk1wOCEGmBwQg0wOKEGGJxQAwxOqAEGJ9QAgxNqgMEJNcDghBpgcEINMLiFQ11VF1TVf1bV56ccCIAXeilH1LckeXCqQQCYb6FQV9UVSY4m+etpxwFguyMLrvfhJL+f5OKdVqiqjSQbSbK6unruk+2ztWN3zl1++vjRfZ4E4IV2PaKuql9Jcqa7T73Yet19orvXu3t9ZWVlzwYEON8tcurjTUl+tapOJ/lkkhuq6u8nnQqA5+wa6u7+w+6+orvXkrw7yZe7+zcmnwyAJK6jBhjeol8mJkm6+ytJvjLJJADM5YgaYHBCDTA4oQYYnFADDE6oAQYn1ACDE2qAwQk1wOCEGmBwQg0wOKEGGJxQAwxOqAEGJ9QAgxNqgMEJNcDghBpgcEINMDihBhicUAMMTqgBBifUAIMTaoDBCTXA4IQaYHC7hrqqfqyq/r2q/quq7q+qP9mPwQDYdGSBdf4vyQ3d/XRVXZjka1X1T939bxPPBkAWCHV3d5KnZ08vnP3qKYcC4HmLHFGnqi5IcirJzyX5SHffPWedjSQbSbK6urqXMw5p7didc5efPn50nycBlt1CXyZ297Pd/cYkVyS5tqpeP2edE9293t3rKysrez0nwHnrJV310d3/m+SuJDdOMw4A2y1y1cdKVb169vjHk7wlyTenHgyATYuco74syW2z89SvSPKp7v78tGMBcNYiV318PcnV+zALAHO4MxFgcEINMDihBhicUAMMTqgBBifUAIMTaoDBCTXA4IQaYHBCDTA4oQYYnFADDE6oAQYn1ACDE2qAwQk1wOCEGmBwQg0wOKEGGJxQAwxOqAEGJ9QAgxNqgMEJNcDghBpgcEINMLhdQ11VV1bVXVX1QFXdX1W37MdgAGw6ssA6P0zye919T1VdnORUVX2pux+YeDYAssARdXc/2t33zB4/leTBJJdPPRgAmxY5on5OVa0luTrJ3XNe20iykSSrq6sve6C1Y3e+pPVPHz/6sj9rdDv9LJZ5m3fiZ8H5bOEvE6vqVUk+k+QD3f397a9394nuXu/u9ZWVlb2cEeC8tlCoq+rCbEb69u7+7LQjAbDVIld9VJKPJXmwuz80/UgAbLXIEfWbkrwnyQ1Vde/s19snnguAmV2/TOzuryWpfZgFgDncmQgwOKEGGJxQAwxOqAEGJ9QAgxNqgMEJNcDghBpgcEINMDihBhicUAMMTqgBBifUAIMTaoDBCTXA4IQaYHBCDTA4oQYYnFADDE6oAQYn1ACDE2qAwQk1wOCEGmBwQg0wuF1DXVW3VtWZqrpvPwYC4IUWOaL+2yQ3TjwHADvYNdTd/dUkT+7DLADMcWSv3qiqNpJsJMnq6upeve1SWTt259zlp48fnfR9pl7/5diPzzgMdvo5JAe7f3hx+70P9uzLxO4+0d3r3b2+srKyV28LcN5z1QfA4IQaYHCLXJ73iST/muSqqnq4qn5r+rEAOGvXLxO7+6b9GASA+Zz6ABicUAMMTqgBBifUAIMTaoDBCTXA4IQaYHBCDTA4oQYYnFADDE6oAQYn1ACDE2qAwQk1wOCEGmBwQg0wOKEGGJxQAwxOqAEGJ9QAgxNqgMEJNcDghBpgcEINMDihBhicUAMMbqFQV9WNVfWtqvp2VR2beigAnrdrqKvqgiQfSfK2JK9LclNVvW7qwQDYtMgR9bVJvt3dD3X3M0k+meQd044FwFnV3S++QtU7k9zY3b89e/6eJL/Q3TdvW28jycbs6VVJvrXl5UuSPLFXQw/I9h1+y76Ntm98P93dK/NeOLJXn9DdJ5KcmPdaVZ3s7vW9+qzR2L7Db9m30fYdbouc+ngkyZVbnl8xWwbAPlgk1P+R5Oer6rVV9cok707yuWnHAuCsXU99dPcPq+rmJP+c5IIkt3b3/S/xc+aeElkitu/wW/ZttH2H2K5fJgJwsNyZCDA4oQYY3KShPh9uPa+q01X1jaq6t6pOHvQ856qqbq2qM1V135ZlP1VVX6qq/5n9/pqDnPFc7LB9H6yqR2b78N6qevtBznguqurKqrqrqh6oqvur6pbZ8mXahztt49Lsx+0mO0c9u/X8v5O8JcnD2bx65KbufmCSDzwgVXU6yXp3H/aL7ZMkVXV9kqeT/F13v3627M+SPNndx2f/w31Nd//BQc75cu2wfR9M8nR3//lBzrYXquqyJJd19z1VdXGSU0l+Lcn7sjz7cKdtfFeWZD9uN+URtVvPD6Hu/mqSJ7ctfkeS22aPb8vmH4pDaYftWxrd/Wh33zN7/FSSB5NcnuXahztt49KaMtSXJ/nulucPZzl/mJ3ki1V1anYb/TK6tLsfnT3+XpJLD3KYidxcVV+fnRo5tKcFtqqqtSRXJ7k7S7oPt21jsoT7MfFl4l64rruvyea/Lvj+2V+tl1Zvnitbtms6P5rkZ5O8McmjSf7iYMc5d1X1qiSfSfKB7v7+1teWZR/O2cal249nTRnq8+LW8+5+ZPb7mSR3ZPOUz7J5bHZe8Oz5wTMHPM+e6u7HuvvZ7v5Rkr/KId+HVXVhNgN2e3d/drZ4qfbhvG1ctv241ZShXvpbz6vqotmXGamqi5K8Ncl9L/5fHUqfS/Le2eP3JvnHA5xlz50N2Myv5xDvw6qqJB9L8mB3f2jLS0uzD3faxmXaj9tNemfi7PKYD+f5W8//dLIPOwBV9TPZPIpONm/H//hh38aq+kSSN2fzn418LMkfJ/mHJJ9KsprkO0ne1d2H8gu5Hbbvzdn863InOZ3kd7aczz1Uquq6JP+S5BtJfjRb/EfZPIe7LPtwp228KUuyH7dzCznA4HyZCDA4oQYYnFADDE6oAQYn1ACDE2qAwQk1wOD+H4CF8fS7yPKWAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "pert_taorf_counts = taorf_meta_moa[['Metadata_moa.x','pert_name']].drop_duplicates().groupby(['Metadata_moa.x']).size()\n",
    "plt.hist(pert_taorf_counts, bins = 50)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Calculate single cells per well/site/perturbation"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(array([2., 2., 2., 3., 1., 4., 2., 1., 0., 1., 2., 0., 0., 2., 0., 0., 1.,\n",
       "        0., 0., 0., 1., 0., 2., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n",
       "        1., 1., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 1.]),\n",
       " array([ 3799.  ,  5588.06,  7377.12,  9166.18, 10955.24, 12744.3 ,\n",
       "        14533.36, 16322.42, 18111.48, 19900.54, 21689.6 , 23478.66,\n",
       "        25267.72, 27056.78, 28845.84, 30634.9 , 32423.96, 34213.02,\n",
       "        36002.08, 37791.14, 39580.2 , 41369.26, 43158.32, 44947.38,\n",
       "        46736.44, 48525.5 , 50314.56, 52103.62, 53892.68, 55681.74,\n",
       "        57470.8 , 59259.86, 61048.92, 62837.98, 64627.04, 66416.1 ,\n",
       "        68205.16, 69994.22, 71783.28, 73572.34, 75361.4 , 77150.46,\n",
       "        78939.52, 80728.58, 82517.64, 84306.7 , 86095.76, 87884.82,\n",
       "        89673.88, 91462.94, 93252.  ]),\n",
       " <a list of 50 Patch objects>)"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXsAAAD4CAYAAAANbUbJAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAARf0lEQVR4nO3df4xlZX3H8fenrOCv1gWZknUXu2tLbKhpC5kghsYYsYpAWJoYAzHtqjSbVmz9leAiSU2TmkA1/kpa7BbQtaH8KGrZoJZSxJgmdXXwB/JDZESE3SzsUBVbTarUb/+4Z+l13d2ZuefO3pl53q/kZs55znPveZ6cmc8897nnnJuqQpK0uv3SpBsgSVp6hr0kNcCwl6QGGPaS1ADDXpIasGbSDQA4/vjja+PGjZNuhiStKHfeeefjVTW1kLrLIuw3btzIzMzMpJshSStKku8utK7TOJLUAMNekhpg2EtSA+YN+yTXJNmX5O6DbHtHkkpyfLeeJB9OMpvkriSnLkWjJUmLs5CR/ceAsw4sTHIi8Erg4aHiVwMndY+twJX9myhJ6mvesK+qLwDfO8imDwCXAMN3UtsMfLwGvgisTbJuLC2VJI1spDn7JJuBPVX19QM2rQceGVrf3ZUd7DW2JplJMjM3NzdKMyRJC7TosE/yTOBdwF/02XFVba+q6aqanppa0DUBkqQRjXJR1a8Dm4CvJwHYAHwlyWnAHuDEobobujJJ0gQtOuyr6hvAr+5fT/IQMF1VjyfZCbw5yfXAi4EnqmrvuBortWTjtk8ftPyhy885wi3RarCQUy+vA/4DeGGS3UkuOkz1zwAPArPA3wNvGksrJUm9zDuyr6oL59m+cWi5gIv7N0uSNE5eQStJDTDsJakBhr0kNcCwl6QGGPaS1ADDXpIaYNhLUgOWxXfQSquJV75qOXJkL0kNMOwlqQFO40gjOtR0jbQcObKXpAYY9pLUAMNekhpg2EtSA/yAVjpCPP9ek+TIXpIaYNhLUgMMe0lqgGEvSQ0w7CWpAfOGfZJrkuxLcvdQ2XuTfDPJXUk+lWTt0LZLk8wmuT/Jq5aq4ZKkhVvIyP5jwFkHlN0GvKiqfhv4FnApQJKTgQuA3+qe87dJjhpbayVJI5k37KvqC8D3Dij716p6slv9IrChW94MXF9V/1NV3wFmgdPG2F5J0gjGMWf/RuCz3fJ64JGhbbu7sl+QZGuSmSQzc3NzY2iGJOlQeoV9ksuAJ4FrF/vcqtpeVdNVNT01NdWnGZKkeYx8u4QkrwfOBc6squqK9wAnDlXb0JVJkiZopJF9krOAS4DzqurHQ5t2AhckOSbJJuAk4Ev9mylJ6mPekX2S64CXAccn2Q28m8HZN8cAtyUB+GJV/UlV3ZPkRuBeBtM7F1fV/y5V4yVJCzNv2FfVhQcpvvow9d8DvKdPoyRJ4+UVtJLUAMNekhpg2EtSAwx7SWqAYS9JDTDsJakBhr0kNWDk2yVIrdi47dOTboLUmyN7SWqAYS9JDTDsJakBhr0kNcCwl6QGGPaS1ADDXpIa4Hn2asrhzpl/6PJzjmBLpCPLkb0kNcCwl6QGGPaS1ADDXpIaYNhLUgPmDfsk1yTZl+TuobLjktyW5IHu57FdeZJ8OMlskruSnLqUjZckLcxCRvYfA846oGwbcHtVnQTc3q0DvBo4qXtsBa4cTzMlSX3MG/ZV9QXgewcUbwZ2dMs7gPOHyj9eA18E1iZZN67GSpJGM+pFVSdU1d5u+VHghG55PfDIUL3dXdleDpBkK4PRP89//vNHbIa08vnlKDoSen9AW1UF1AjP215V01U1PTU11bcZkqTDGDXsH9s/PdP93NeV7wFOHKq3oSuTJE3QqNM4O4EtwOXdz5uHyt+c5HrgxcATQ9M9ksbgUNM+3ttHhzNv2Ce5DngZcHyS3cC7GYT8jUkuAr4LvLar/hngbGAW+DHwhiVos7QknDvXajZv2FfVhYfYdOZB6hZwcd9GSZLGy1scS6uE0zs6HMNeq5JTMtLP8944ktQAw16SGmDYS1IDDHtJaoBhL0kNMOwlqQGGvSQ1wLCXpAYY9pLUAMNekhpg2EtSAwx7SWqAYS9JDTDsJakBhr0kNcCwl6QGGPaS1ADDXpIaYNhLUgN6hX2StyW5J8ndSa5L8vQkm5LsSjKb5IYkR4+rsZKk0Ywc9knWA38OTFfVi4CjgAuAK4APVNVvAN8HLhpHQyVJo+s7jbMGeEaSNcAzgb3Ay4Gbuu07gPN77kOS1NPIYV9Ve4D3AQ8zCPkngDuBH1TVk1213cD6gz0/ydYkM0lm5ubmRm2GJGkB+kzjHAtsBjYBzwOeBZy10OdX1faqmq6q6ampqVGbIUlagD7TOK8AvlNVc1X1U+CTwBnA2m5aB2ADsKdnGyVJPfUJ+4eB05M8M0mAM4F7gTuA13R1tgA392uiJKmvPnP2uxh8EPsV4Bvda20H3gm8Pcks8Fzg6jG0U5LUw5r5qxxaVb0bePcBxQ8Cp/V5XUnSeHkFrSQ1wLCXpAYY9pLUAMNekhpg2EtSAwx7SWqAYS9JDeh1nr2k5W/jtk8ftPyhy885wi3RJDmyl6QGGPaS1ADDXpIaYNhLUgMMe0lqgGEvSQ0w7CWpAZ5nrxXBc8WlfhzZS1IDDHtJaoBhL0kNMOwlqQGGvSQ1oNfZOEnWAlcBLwIKeCNwP3ADsBF4CHhtVX2/VyulQzjUWTqSfl7fkf2HgH+pqt8Efge4D9gG3F5VJwG3d+uSpAkaOeyTPAd4KXA1QFX9pKp+AGwGdnTVdgDn922kJKmfPiP7TcAc8NEkX01yVZJnASdU1d6uzqPACQd7cpKtSWaSzMzNzfVohiRpPn3Cfg1wKnBlVZ0C/IgDpmyqqhjM5f+CqtpeVdNVNT01NdWjGZKk+fQJ+93A7qra1a3fxCD8H0uyDqD7ua9fEyVJfY0c9lX1KPBIkhd2RWcC9wI7gS1d2Rbg5l4tlCT11vdGaH8GXJvkaOBB4A0M/oHcmOQi4LvAa3vuQ5LUU6+wr6qvAdMH2XRmn9eVJI2XV9BKUgMMe0lqgGEvSQ0w7CWpAX4toZYVb2wmLQ1H9pLUAEf2UqNGeRflF7yvXI7sJakBhr0kNcCwl6QGGPaS1ADDXpIaYNhLUgMMe0lqgGEvSQ0w7CWpAYa9JDXAsJekBhj2ktQAw16SGmDYS1IDeod9kqOSfDXJLd36piS7kswmuSHJ0f2bKUnqYxz3s38LcB/wK936FcAHqur6JB8BLgKuHMN+tAId6p7p3hddOrJ6jeyTbADOAa7q1gO8HLipq7IDOL/PPiRJ/fWdxvkgcAnws279ucAPqurJbn03sP5gT0yyNclMkpm5ubmezZAkHc7IYZ/kXGBfVd05yvOrantVTVfV9NTU1KjNkCQtQJ85+zOA85KcDTydwZz9h4C1SdZ0o/sNwJ7+zZQk9THyyL6qLq2qDVW1EbgA+FxVvQ64A3hNV20LcHPvVkqSehnH2TgHeidwfZK/Ar4KXL0E+9AKd6izdLS8eXbVyjWWsK+qzwOf75YfBE4bx+tKksbDK2glqQGGvSQ1wLCXpAYY9pLUAMNekhpg2EtSAwx7SWqAYS9JDTDsJakBhr0kNWAp7o2jBnmvG2l5c2QvSQ0w7CWpAU7jSOrNWx8vf47sJakBhr0kNcCwl6QGGPaS1ADDXpIaYNhLUgMMe0lqwMhhn+TEJHckuTfJPUne0pUfl+S2JA90P48dX3MlSaPoM7J/EnhHVZ0MnA5cnORkYBtwe1WdBNzerUuSJmjksK+qvVX1lW75v4D7gPXAZmBHV20HcH7fRkqS+hnLnH2SjcApwC7ghKra2216FDjhEM/ZmmQmyczc3Nw4miFJOoTe98ZJ8mzgE8Bbq+qHSZ7aVlWVpA72vKraDmwHmJ6ePmgdTY73OpFWl14j+yRPYxD011bVJ7vix5Ks67avA/b1a6Ikqa+RR/YZDOGvBu6rqvcPbdoJbAEu737e3KuFWlb8khKNg+8cj7w+0zhnAH8IfCPJ17qydzEI+RuTXAR8F3htvyZKWqkcHCwfI4d9Vf07kENsPnPU15UkjZ9X0EpSAwx7SWqAYS9JDTDsJakBhr0kNcCwl6QGGPaS1IDe98bR8ne4C1u8YlFqgyN7SWqAYS9JDXAap3Heu0RqgyN7SWqAYS9JDXAaR9Ky4X3ul44je0lqwIof2XsOuSTNb8WH/Wq22Le0nlmj1crpnf6cxpGkBjiyX4ClHlUsdkTuCF46vMX+jbTwDqHJsJ9UeLfwCyVpeXIaR5IakKpamhdOzgI+BBwFXFVVlx+q7vT0dM3MzIy0n0lOafhBqbS6Lfd340nurKrphdRdkmmcJEcBfwP8PrAb+HKSnVV171Lsb1IMdWl1G+VvfLGDwCP1D2WppnFOA2ar6sGq+glwPbB5ifYlSZrHUn1Aux54ZGh9N/Di4QpJtgJbu9X/TnL/ErXlSDseeHzSjVgiq7lvsLr7Z9+OkFwx1vrz9e3XFrqfiZ2NU1Xbge2T2v9SSTKz0Dm0lWY19w1Wd//s28o0zr4t1TTOHuDEofUNXZkkaQKWKuy/DJyUZFOSo4ELgJ1LtC9J0jyWZBqnqp5M8mbgVganXl5TVfcsxb6WoVU3NTVkNfcNVnf/7NvKNLa+Ldl59pKk5cMraCWpAYa9JDXAsF+kJCcmuSPJvUnuSfKWrvy4JLcleaD7eWxXniQfTjKb5K4kp062B4eX5KgkX01yS7e+Kcmurv03dB+4k+SYbn22275xku1eiCRrk9yU5JtJ7kvyklV03N7W/T7eneS6JE9fqccuyTVJ9iW5e6hs0ccpyZau/gNJtkyiLwc6RN/e2/1O3pXkU0nWDm27tOvb/UleNVR+Vlc2m2TbgnZeVT4W8QDWAad2y78MfAs4GfhrYFtXvg24ols+G/gsEOB0YNek+zBP/94O/CNwS7d+I3BBt/wR4E+75TcBH+mWLwBumHTbF9C3HcAfd8tHA2tXw3FjcBHjd4BnDB2z16/UYwe8FDgVuHuobFHHCTgOeLD7eWy3fOwy7dsrgTXd8hVDfTsZ+DpwDLAJ+DaDE16O6pZf0P0efx04ed59T7rzK/0B3MzgHkD3A+u6snXA/d3y3wEXDtV/qt5yezC4HuJ24OXALd0f0ONDv4gvAW7tlm8FXtItr+nqZdJ9OEzfntMFYg4oXw3Hbf8V68d1x+IW4FUr+dgBGw8IxEUdJ+BC4O+Gyn+u3nLq2wHb/gC4tlu+FLh0aNut3XF86lgerN6hHk7j9NC9/T0F2AWcUFV7u02PAid0ywe7dcT6I9TExfogcAnws279ucAPqurJbn247U/1q9v+RFd/udoEzAEf7aaprkryLFbBcauqPcD7gIeBvQyOxZ2snmMHiz9OK+b4HeCNDN6pwJj7ZtiPKMmzgU8Ab62qHw5vq8G/2xV1TmuSc4F9VXXnpNuyRNYwePt8ZVWdAvyIwXTAU1bicQPo5q83M/iH9jzgWcBZE23UElqpx2k+SS4DngSuXYrXN+xHkORpDIL+2qr6ZFf8WJJ13fZ1wL6ufKXcOuIM4LwkDzG4S+nLGXwfwdok+y++G277U/3qtj8H+M8j2eBF2g3srqpd3fpNDMJ/pR83gFcA36mquar6KfBJBsdztRw7WPxxWknHjySvB84FXtf9M4Mx982wX6QkAa4G7quq9w9t2gns/8R/C4O5/P3lf9SdNXA68MTQ29Flo6ouraoNVbWRwYd2n6uq1wF3AK/pqh3Yr/39fU1Xf9mOtqrqUeCRJC/sis4E7mWFH7fOw8DpSZ7Z/X7u79uqOHadxR6nW4FXJjm2e+fzyq5s2cngi54uAc6rqh8PbdoJXNCdPbUJOAn4EqPejmbSH1astAfwewzeQt4FfK17nM1gzvN24AHg34Djuvph8EUu3wa+AUxPug8L6OPL+P+zcV7Q/YLNAv8EHNOVP71bn+22v2DS7V5Av34XmOmO3T8zOEtjVRw34C+BbwJ3A//A4AyOFXnsgOsYfPbwUwbvyC4a5TgxmP+e7R5vmHS/DtO3WQZz8Pvz5CND9S/r+nY/8Oqh8rMZnAn4beCyhezb2yVIUgOcxpGkBhj2ktQAw16SGmDYS1IDDHtJaoBhL0kNMOwlqQH/B5vLBGfhlzmWAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXEAAAD4CAYAAAAaT9YAAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAOw0lEQVR4nO3dXYxc9XnH8e9Tm5cUKLZha1kYZ42CiHzRAFoREChqoVBeIuACIaMotVpHltpEArVSahqpUqRcmF7kpVLVxAq0vkh4KUlqC9QS1xBVrSqTNa8Gx7WhRrFls04KIe1FW5OnF/M3GdxZdvZlZueB70dazTn/c2bnt3D48d8z58xGZiJJqulXFjuAJGnuLHFJKswSl6TCLHFJKswSl6TClg7zxc4///wcHx8f5ktKUnl79uz5SWaO9do21BIfHx9ncnJymC8pSeVFxGvTbfN0iiQV1tdMPCIOAT8H3gZOZOZERKwAHgbGgUPAnZn5xmBiSpJ6mc1M/Lcy89LMnGjrm4FdmXkxsKutS5KGaD6nU24DtrXlbcDt848jSZqNfks8ge9HxJ6I2NTGVmbm0bZ8DFjZ64kRsSkiJiNi8vjx4/OMK0nq1u/VKddk5pGI+HVgZ0T8qHtjZmZE9PwkrczcCmwFmJiY8NO2JGkB9TUTz8wj7XEK+B5wBfB6RKwCaI9TgwopSeptxhKPiLMi4pyTy8ANwF5gB7Ch7bYB2D6okJKk3vo5nbIS+F5EnNz/25n5DxHxQ+CRiNgIvAbcObiYkqReZizxzHwV+FiP8Z8C1w0ilDRs45sf7zl+aMstQ04izY53bEpSYZa4JBVmiUtSYZa4JBVmiUtSYZa4JBVmiUtSYZa4JBVmiUtSYZa4JBVmiUtSYZa4JBVmiUtSYZa4JBVmiUtSYZa4JBVmiUtSYZa4JBVmiUtSYZa4JBVmiUtSYZa4JBVmiUtSYZa4JBVmiUtSYZa4JBVmiUtSYZa4JBVmiUtSYZa4JBVmiUtSYZa4JBVmiUtSYZa4JBXWd4lHxJKIeDYiHmvrayNid0QcjIiHI+L0wcWUJPUym5n43cC+rvX7gK9k5keAN4CNCxlMkjSzvko8IlYDtwDfbOsBXAs82nbZBtw+iICSpOn1OxP/KvB54Bdt/Tzgzcw80dYPAxf0emJEbIqIyYiYPH78+LzCSpLebcYSj4hPAlOZuWcuL5CZWzNzIjMnxsbG5vItJEnTWNrHPlcDt0bEzcCZwK8BXwOWRcTSNhtfDRwZXExJUi8zzsQz897MXJ2Z48B64MnM/BTwFHBH220DsH1gKSVJPfUzE5/OnwAPRcSXgGeB+xcmkjR/45sf7zl+aMstQ04iDdasSjwzfwD8oC2/Clyx8JEkSf3yjk1JKswSl6TCLHFJKswSl6TCLHFJKswSl6TCLHFJKswSl6TCLHFJKswSl6TCLHFJKswSl6TCLHFJKswSl6TCLHFJKmw+fxRCWnTT/fGHhdpfGnXOxCWpMEtckgqzxCWpMEtckgqzxCWpMEtckgqzxCWpMEtckgqzxCWpMEtckgqzxCWpMEtckgqzxCWpMEtckgqzxCWpMEtckgqzxCWpsBlLPCLOjIinI+L5iHgpIr7YxtdGxO6IOBgRD0fE6YOPK0nq1s9M/L+BazPzY8ClwI0RcSVwH/CVzPwI8AawcXAxJUm9zFji2fGfbfW09pXAtcCjbXwbcPtAEkqSptXXOfGIWBIRzwFTwE7gFeDNzDzRdjkMXDDNczdFxGRETB4/fnwhMkuSmr5KPDPfzsxLgdXAFcBH+32BzNyamROZOTE2NjbHmJKkXmZ1dUpmvgk8BVwFLIuIpW3TauDIAmeTJM2gn6tTxiJiWVv+EHA9sI9Omd/RdtsAbB9USElSb0tn3oVVwLaIWEKn9B/JzMci4mXgoYj4EvAscP8Ac0qSepixxDPzBeCyHuOv0jk/LklaJN6xKUmFWeKSVJglLkmFWeKSVJglLkmFWeKSVJglLkmFWeKSVJglLkmFWeKSVJglLkmFWeKSVJglLkmFWeKSVJglLkmF9fNHIaShGd/8eM/xQ1tuGXISqQZn4pJUmCUuSYVZ4pJUmOfEpffgOXqNOmfiklSYJS5JhVniklSYJS5JhVniklSYJS5JhVniklSYJS5JhXmzj0qY7qYb6YPOmbgkFWaJS1JhlrgkFWaJS1JhM5Z4RFwYEU9FxMsR8VJE3N3GV0TEzog40B6XDz6uJKlbPzPxE8AfZ+Y64ErgsxGxDtgM7MrMi4FdbV2SNEQzlnhmHs3MZ9ryz4F9wAXAbcC2tts24PZBhZQk9Tar68QjYhy4DNgNrMzMo23TMWDlNM/ZBGwCWLNmzVxz6n3G676lhdH3G5sRcTbwHeCezHyre1tmJpC9npeZWzNzIjMnxsbG5hVWkvRufZV4RJxGp8C/lZnfbcOvR8Sqtn0VMDWYiJKk6fRzdUoA9wP7MvPLXZt2ABva8gZg+8LHkyS9l37OiV8NfBp4MSKea2N/CmwBHomIjcBrwJ2DiShJms6MJZ6Z/wzENJuvW9g4kqTZ8I5NSSrMEpekwixxSSrMEpekwixxSSrMEpekwixxSSrMEpekwixxSSrMEpekwixxSSpsVn8UQlLHdH/U4tCWW4acRB90zsQlqTBLXJIKs8QlqTBLXJIKs8QlqTBLXJIKs8QlqTBLXJIKs8QlqTBLXJIKs8QlqTBLXJIKs8QlqTBLXJIKs8QlqTBLXJIKs8QlqTBLXJIKs8QlqTBLXJIKs8QlqbAZSzwiHoiIqYjY2zW2IiJ2RsSB9rh8sDElSb30MxP/G+DGU8Y2A7sy82JgV1uXJA3ZjCWemf8E/Mcpw7cB29ryNuD2Bc4lSerDXM+Jr8zMo235GLBygfJIkmZh6Xy/QWZmROR02yNiE7AJYM2aNfN9OWmkjW9+fNpth7bcMsQk+qCY60z89YhYBdAep6bbMTO3ZuZEZk6MjY3N8eUkSb3MtcR3ABva8gZg+8LEkSTNRj+XGD4I/CtwSUQcjoiNwBbg+og4APx2W5ckDdmM58Qz865pNl23wFkkSbPkHZuSVJglLkmFWeKSVNi8rxOXYPrro702WhosZ+KSVJglLkmFWeKSVJjnxNW39/pcEEmLw5m4JBVmiUtSYZa4JBXmOXFpSLyWXoPgTFySCrPEJakwS1ySCrPEJakw39h8HxnFN868QUgaLGfiklSYJS5JhVniklSYJS5JhVniklSYJS5JhVniklSY14lLi2wUr+9XHc7EJakwS1ySCrPEJakwz4mPgEGfE53t55d4Llaqw5m4JBVmiUtSYZa4JBVmiUtSYWXe2Jztm3+D3n+uz5mNxboJxD/k8P7izUTDsxj/rOc1E4+IGyNif0QcjIjNCxVKktSfOZd4RCwB/hK4CVgH3BUR6xYqmCRpZvOZiV8BHMzMVzPzf4CHgNsWJpYkqR+RmXN7YsQdwI2Z+Zm2/mng45n5uVP22wRsaquXAPu7Np8P/GROARZXxdwVM0PN3BUzQ83cFTPD7HN/ODPHem0Y+BubmbkV2NprW0RMZubEoDMstIq5K2aGmrkrZoaauStmhoXNPZ/TKUeAC7vWV7cxSdKQzKfEfwhcHBFrI+J0YD2wY2FiSZL6MefTKZl5IiI+BzwBLAEeyMyXZvltep5mKaBi7oqZoWbuipmhZu6KmWEBc8/5jU1J0uLztntJKswSl6TCFq3ER+mW/Yh4ICKmImJv19iKiNgZEQfa4/I2HhHxFy33CxFxeddzNrT9D0TEhgFnvjAinoqIlyPipYi4u0juMyPi6Yh4vuX+YhtfGxG7W76H25vlRMQZbf1g2z7e9b3ubeP7I+J3Bpm7vd6SiHg2Ih4rlPlQRLwYEc9FxGQbG/VjZFlEPBoRP4qIfRFxVYHMl7R/xie/3oqIe4aSOzOH/kXnjdBXgIuA04HngXWLkaXl+QRwObC3a+zPgc1teTNwX1u+Gfh7IIArgd1tfAXwantc3paXDzDzKuDytnwO8G90Pv5g1HMHcHZbPg3Y3fI8Aqxv418H/qAt/yHw9ba8Hni4La9rx80ZwNp2PC0Z8HHyR8C3gcfaeoXMh4DzTxkb9WNkG/CZtnw6sGzUM5+SfwlwDPjwMHIP/Aea5oe8Cniia/1e4N7FyNKVYZx3l/h+YFVbXgXsb8vfAO46dT/gLuAbXePv2m8I+bcD11fKDfwq8AzwcTp3ry099figc/XTVW15adsvTj1muvcbUNbVwC7gWuCxlmGkM7fXOMT/L/GRPUaAc4F/p110USFzj5/hBuBfhpV7sU6nXAD8uGv9cBsbJSsz82hbPgasbMvTZV+0n6n9un4ZnVntyOdupyWeA6aAnXRmpG9m5okeGd7J17b/DDhvEXJ/Ffg88Iu2fl6BzAAJfD8i9kTnIzBgtI+RtcBx4K/bqatvRsRZI575VOuBB9vywHP7xmYfsvO/xJG8FjMizga+A9yTmW91bxvV3Jn5dmZeSmd2ewXw0UWO9J4i4pPAVGbuWewsc3BNZl5O59NGPxsRn+jeOILHyFI6pzb/KjMvA/6LzmmId4xg5ne090VuBf721G2Dyr1YJV7hlv3XI2IVQHucauPTZR/6zxQRp9Ep8G9l5ner5D4pM98EnqJzKmJZRJy8+aw7wzv52vZzgZ8y3NxXA7dGxCE6n9Z5LfC1Ec8MQGYeaY9TwPfo/E9zlI+Rw8DhzNzd1h+lU+qjnLnbTcAzmfl6Wx947sUq8Qq37O8ATr4zvIHOOeeT47/b3l2+EvhZ+3XpCeCGiFje3oG+oY0NREQEcD+wLzO/XCj3WEQsa8sfonMefx+dMr9jmtwnf547gCfbjGYHsL5dCbIWuBh4ehCZM/PezFydmeN0jtUnM/NTo5wZICLOiohzTi7T+Xe7lxE+RjLzGPDjiLikDV0HvDzKmU9xF788lXIy32BzD+NE/zQn/2+mc0XFK8AXFitHy/IgcBT4XzozgY10zmHuAg4A/wisaPsGnT+G8QrwIjDR9X1+HzjYvn5vwJmvofOr2QvAc+3r5gK5fwN4tuXeC/xZG7+ITqEdpPOr6Blt/My2frBtv6jre32h/Tz7gZuGdKz8Jr+8OmWkM7d8z7evl07+d1bgGLkUmGzHyN/RuUpjpDO31zuLzm9c53aNDTy3t91LUmG+sSlJhVniklSYJS5JhVniklSYJS5JhVniklSYJS5Jhf0ftEf4xGuU2ZUAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAD4CAYAAAD8Zh1EAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAARTUlEQVR4nO3dfYxldX3H8ffHZXlQqaBMlO6Dg4GYoFHACWJoGgJVVzBsEzFZ0wpYySZWIrYmBjTByF/QNGoRI9kAFagVLBK78hBCC0b9g9VZ5HlBR6WyBMvytEhV7Npv/7hHvDvO7L0ze2dn9zfvV3Iz5/zO7577ncNvP5w59zykqpAk7ftettgFSJJGw0CXpEYY6JLUCANdkhphoEtSI/ZbrA8+7LDDanx8fLE+XpL2SZs3b36qqsZmWrZogT4+Ps7k5ORifbwk7ZOS/NdsyzzkIkmNMNAlqREGuiQ1wkCXpEYY6JLUCANdkhoxdKAnWZbkh0lummHZAUmuTzKVZFOS8VEWKUkabC576OcBW2ZZ9mHg2ao6Evg8cMnuFiZJmpuhAj3JSuA04IpZuqwFru6mbwBOSZLdL0+SNKxhrxT9AvBJ4OBZlq8AHgOoqh1JtgOvAZ7q75RkPbAeYPXq1fOpd9GNn3/zjO2PXnzaHq5EknY2cA89yXuBJ6tq8+5+WFVtqKqJqpoYG5vxVgSSpHka5pDLicDpSR4FrgNOTvIv0/o8DqwCSLIf8Crg6RHWKUkaYGCgV9UFVbWyqsaBdcAdVfXX07ptBM7qps/o+viwUknag+Z9t8UkFwGTVbURuBK4NskU8Ay94Jck7UFzCvSq+jbw7W76wr723wDvH2VhkqS58UpRSWqEgS5JjTDQJakRBrokNcJAl6RGGOiS1AgDXZIaYaBLUiMMdElqhIEuSY0w0CWpEQa6JDXCQJekRhjoktQIA12SGmGgS1IjhnlI9IFJvp/k3iQPJvnsDH3OTrItyT3d65yFKVeSNJthnlj0InByVb2QZDnwvSS3VtVd0/pdX1Xnjr5ESdIwBgZ697DnF7rZ5d3LB0BL0l5mqGPoSZYluQd4Eri9qjbN0O19Se5LckOSVSOtUpI00FCBXlW/q6pjgJXA8UnePK3Lt4DxqnoLcDtw9UzrSbI+yWSSyW3btu1O3ZKkaeZ0lktVPQfcCayZ1v50Vb3YzV4BvG2W92+oqomqmhgbG5tPvZKkWQxzlstYkkO66YOAdwIPT+tzeN/s6cCWURYpSRpsmLNcDgeuTrKM3v8Avl5VNyW5CJisqo3Ax5KcDuwAngHOXqiCJUkzG+Ysl/uAY2dov7Bv+gLggtGWJkmaC68UlaRGGOiS1AgDXZIaYaBLUiMMdElqhIEuSY0w0CWpEQa6JDXCQJekRhjoktQIA12SGmGgS1IjDHRJaoSBLkmNMNAlqREGuiQ1wkCXpEYM80zRA5N8P8m9SR5M8tkZ+hyQ5PokU0k2JRlfiGIlSbMbZg/9ReDkqnorcAywJskJ0/p8GHi2qo4EPg9cMtoyJUmDDAz06nmhm13evWpat7XA1d30DcApSTKyKiVJAw18SDRAkmXAZuBI4EtVtWlalxXAYwBVtSPJduA1wFPT1rMeWA+wevXq3at8RMbPv3nG9kcvPm0PVyJJu2eoL0Wr6ndVdQywEjg+yZvn82FVtaGqJqpqYmxsbD6rkCTNYk5nuVTVc8CdwJppix4HVgEk2Q94FfD0KAqUJA1nmLNcxpIc0k0fBLwTeHhat43AWd30GcAdVTX9OLskaQENcwz9cODq7jj6y4CvV9VNSS4CJqtqI3AlcG2SKeAZYN2CVSxJmtHAQK+q+4BjZ2i/sG/6N8D7R1uaJGkuvFJUkhphoEtSIwx0SWqEgS5JjTDQJakRBrokNcJAl6RGGOiS1AgDXZIaYaBLUiMMdElqhIEuSY0w0CWpEQa6JDXCQJekRhjoktQIA12SGjHMM0VXJbkzyUNJHkxy3gx9TkqyPck93evCmdYlSVo4wzxTdAfwiaq6O8nBwOYkt1fVQ9P6fbeq3jv6EiVJwxi4h15VT1TV3d30L4EtwIqFLkySNDdzOoaeZJzeA6M3zbD4HUnuTXJrkjfN8v71SSaTTG7btm3OxUqSZjd0oCd5JfAN4ONV9fy0xXcDr6+qtwJfBL450zqqakNVTVTVxNjY2HxrliTNYKhAT7KcXph/tapunL68qp6vqhe66VuA5UkOG2mlkqRdGuYslwBXAluq6nOz9Hld148kx3frfXqUhUqSdm2Ys1xOBD4I3J/knq7tU8BqgKq6HDgD+EiSHcCvgXVVVQtQryRpFgMDvaq+B2RAn8uAy0ZVlCRp7rxSVJIaYaBLUiMMdElqhIEuSY0w0CWpEQa6JDXCQJekRhjoktQIA12SGmGgS1IjDHRJaoSBLkmNMNAlqREGuiQ1wkCXpEYY6JLUCANdkhoxzDNFVyW5M8lDSR5Mct4MfZLk0iRTSe5LctzClCtJms0wzxTdAXyiqu5OcjCwOcntVfVQX5/3AEd1r7cDX+5+SpL2kIF76FX1RFXd3U3/EtgCrJjWbS1wTfXcBRyS5PCRVytJmtUwe+gvSTIOHAtsmrZoBfBY3/zWru2Jae9fD6wHWL169dwq7TN+/s3zfu9Cma2mRy8+ba9b/0LXulha/b2kYQ39pWiSVwLfAD5eVc/P58OqakNVTVTVxNjY2HxWIUmaxVCBnmQ5vTD/alXdOEOXx4FVffMruzZJ0h4yzFkuAa4EtlTV52bpthE4szvb5QRge1U9MUtfSdICGOYY+onAB4H7k9zTtX0KWA1QVZcDtwCnAlPAr4APjb5USdKuDAz0qvoekAF9CvjoqIqSJM2dV4pKUiMMdElqhIEuSY0w0CWpEQa6JDXCQJekRhjoktQIA12SGmGgS1IjDHRJaoSBLkmNMNAlqREGuiQ1wkCXpEYY6JLUCANdkhphoEtSI4Z5puhVSZ5M8sAsy09Ksj3JPd3rwtGXKUkaZJhnin4FuAy4Zhd9vltV7x1JRZKkeRm4h15V3wGe2QO1SJJ2w6iOob8jyb1Jbk3yptk6JVmfZDLJ5LZt20b00ZIkGE2g3w28vqreCnwR+OZsHatqQ1VNVNXE2NjYCD5akvR7ux3oVfV8Vb3QTd8CLE9y2G5XJkmak90O9CSvS5Ju+vhunU/v7nolSXMz8CyXJF8DTgIOS7IV+AywHKCqLgfOAD6SZAfwa2BdVdWCVSxJmtHAQK+qDwxYfhm90xolSYvIK0UlqREGuiQ1wkCXpEYY6JLUCANdkhphoEtSIwx0SWqEgS5JjTDQJakRBrokNcJAl6RGGOiS1AgDXZIaYaBLUiMMdElqhIEuSY0w0CWpEQMDPclVSZ5M8sAsy5Pk0iRTSe5Lctzoy5QkDTLMHvpXgDW7WP4e4KjutR748u6XJUmaq4GBXlXfAZ7ZRZe1wDXVcxdwSJLDR1WgJGk4Ax8SPYQVwGN981u7tiemd0yynt5ePKtXrx7BRy+c8fNvXtD1PHrxaQv6uXvCXH+3vW39S9F8xpfbe/729Bjeo1+KVtWGqpqoqomxsbE9+dGS1LxRBPrjwKq++ZVdmyRpDxpFoG8EzuzOdjkB2F5Vf3S4RZK0sAYeQ0/yNeAk4LAkW4HPAMsBqupy4BbgVGAK+BXwoYUqVpI0u4GBXlUfGLC8gI+OrCJJ0rx4pagkNcJAl6RGGOiS1AgDXZIaYaBLUiMMdElqhIEuSY0w0CWpEQa6JDXCQJekRhjoktQIA12SGmGgS1IjDHRJaoSBLkmNMNAlqREGuiQ1YqhAT7ImySNJppKcP8Pys5NsS3JP9zpn9KVKknZlmGeKLgO+BLwT2Ar8IMnGqnpoWtfrq+rcBahRkjSEYfbQjwemquqnVfVb4Dpg7cKWJUmaq2ECfQXwWN/81q5tuvcluS/JDUlWzbSiJOuTTCaZ3LZt2zzKlSTNZlRfin4LGK+qtwC3A1fP1KmqNlTVRFVNjI2NjeijJUkwXKA/DvTvca/s2l5SVU9X1Yvd7BXA20ZTniRpWMME+g+Ao5IckWR/YB2wsb9DksP7Zk8HtoyuREnSMAae5VJVO5KcC9wGLAOuqqoHk1wETFbVRuBjSU4HdgDPAGcvYM2SpBkMDHSAqroFuGVa24V90xcAF4y2NEnSXHilqCQ1wkCXpEYY6JLUCANdkhphoEtSIwx0SWqEgS5JjTDQJakRBrokNcJAl6RGGOiS1AgDXZIaYaBLUiMMdElqhIEuSY0w0CWpEQa6JDViqEBPsibJI0mmkpw/w/IDklzfLd+UZHzUhUqSdm1goCdZBnwJeA9wNPCBJEdP6/Zh4NmqOhL4PHDJqAuVJO3aMHvoxwNTVfXTqvotcB2wdlqftcDV3fQNwClJMroyJUmDpKp23SE5A1hTVed08x8E3l5V5/b1eaDrs7Wb/0nX56lp61oPrO9m3wg8AhwG7NRPbpNp3B47c3vsbKltj9dX1dhMC/bbk1VU1QZgQ39bksmqmtiTdezt3CY7c3vszO2xM7fHHwxzyOVxYFXf/MqubcY+SfYDXgU8PYoCJUnDGSbQfwAcleSIJPsD64CN0/psBM7qps8A7qhBx3IkSSM18JBLVe1Ici5wG7AMuKqqHkxyETBZVRuBK4Frk0wBz9AL/WFtGNxlyXGb7MztsTO3x87cHp2BX4pKkvYNXikqSY0w0CWpEYsa6INuKbAvS7IqyZ1JHkryYJLzuvZXJ7k9yY+7n4d27Ulyabct7ktyXN+6zur6/zjJWX3tb0tyf/eeS/f2i7mSLEvywyQ3dfNHdLeKmOpuHbF/1z7rrSSSXNC1P5Lk3X3t+9xYSnJIkhuSPJxkS5J3LPHx8Xfdv5UHknwtyYFLfYzMWVUtyoveF6w/Ad4A7A/cCxy9WPUswO93OHBcN30w8CN6t074B+D8rv184JJu+lTgViDACcCmrv3VwE+7n4d204d2y77f9U333vcs9u89YJv8PfCvwE3d/NeBdd305cBHuum/BS7vptcB13fTR3fj5ADgiG78LNtXxxK9q6vP6ab3Bw5ZquMDWAH8DDiob2ycvdTHyFxfi7mHPswtBfZZVfVEVd3dTf8S2EJv0PbfJuFq4C+76bXANdVzF3BIksOBdwO3V9UzVfUscDuwplv2J1V1V/VG8jV969rrJFkJnAZc0c0HOJnerSLgj7fFTLeSWAtcV1UvVtXPgCl642ifG0tJXgX8Ob0zxKiq31bVcyzR8dHZDziou5bl5cATLOExMh+LGegrgMf65rd2bc3p/hw8FtgEvLaqnugW/QJ4bTc92/bYVfvWGdr3Vl8APgn8Xzf/GuC5qtrRzffX/9Lv3C3f3vWf6zbamx0BbAP+uTsMdUWSV7BEx0dVPQ78I/BzekG+HdjM0h4jc+aXogssySuBbwAfr6rn+5d1e07Nnzea5L3Ak1W1ebFr2YvsBxwHfLmqjgX+h94hlpcslfEB0H1XsJbe/+j+FHgFsGZRi9oHLWagD3NLgX1akuX0wvyrVXVj1/zf3Z/DdD+f7Npn2x67al85Q/ve6ETg9CSP0vtT92Tgn+gdNvj9xW399c92K4m5bqO92VZga1Vt6uZvoBfwS3F8APwF8LOq2lZV/wvcSG/cLOUxMmeLGejD3FJgn9Udz7sS2FJVn+tb1H+bhLOAf+9rP7M7m+EEYHv3p/dtwLuSHNrtxbwLuK1b9nySE7rPOrNvXXuVqrqgqlZW1Ti9/853VNVfAXfSu1UE/PG2mOlWEhuBdd0ZDkcAR9H74m+fG0tV9QvgsSRv7JpOAR5iCY6Pzs+BE5K8vKv399tjyY6ReVnMb2TpfXP/I3rfPn96MWtZgN/tz+j9uXwfcE/3OpXecb7/BH4M/Afw6q5/6D1I5CfA/cBE37r+ht6XO1PAh/raJ4AHuvdcRnfl7978Ak7iD2e5vIHeP7Yp4N+AA7r2A7v5qW75G/re/+nu932EvrM29sWxBBwDTHZj5Jv0zlJZsuMD+CzwcFfztfTOVFnSY2SuLy/9l6RG+KWoJDXCQJekRhjoktQIA12SGmGgS1IjDHRJaoSBLkmN+H9dF+JUtTqtNAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure()\n",
    "plt.hist(taorf_meta.groupby(['Metadata_Plate','Metadata_Well']).agg('sum')['Cells'], bins = 50)\n",
    "plt.figure()\n",
    "#Cut EMPTY, \n",
    "plt.hist(taorf_meta[~taorf_meta['pert_name'].isin(['EMPTY_', 'Luciferase_CTRL', 'LacZ_CTRL', 'eGFP_CTRL'])].groupby(['pert_name']).agg('sum')['Cells'], bins = 50)\n",
    "plt.figure()\n",
    "plt.hist(taorf_meta_moa.groupby(['Metadata_moa.x']).agg('sum')['Cells'], bins = 50)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Calculate stats for the full BBBC022 dataset (after QC)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "bbbc022_meta = pd.read_csv(\"/dgx1nas1/cellpainting-datasets/BBBC022/inputs/metadata/index_after_qc.csv\")\n",
    "#the number of cells is inferred from locations\n",
    "cells_bbbc022 = []\n",
    "for i in bbbc022_meta.index:\n",
    "    cells_bbbc022.append([\n",
    "        bbbc022_meta.loc[i, \"Metadata_Plate\"],\n",
    "        bbbc022_meta.loc[i, \"Metadata_Well\"],\n",
    "        bbbc022_meta.loc[i, \"Metadata_Site\"],\n",
    "        sum(1 for line in open('/dgx1nas1/cellpainting-datasets/BBBC022/inputs/locations/{}/{}-{}-Nuclei.csv'.format(bbbc022_meta.loc[i, \"Metadata_Plate\"], bbbc022_meta.loc[i, \"Metadata_Well\"], bbbc022_meta.loc[i, \"Metadata_Site\"]))) - 1\n",
    "    ])\n",
    "counts = pd.DataFrame(data=cells_bbbc022, columns=['Metadata_Plate', 'Metadata_Well', 'Metadata_Site', 'Cells'])\n",
    "bbbc022_meta = pd.merge(bbbc022_meta, counts, right_on=['Metadata_Plate', 'Metadata_Well', 'Metadata_Site'], left_on =['Metadata_Plate', 'Metadata_Well', 'Metadata_Site'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "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>Metadata_Plate</th>\n",
       "      <th>Metadata_Well</th>\n",
       "      <th>Metadata_Site</th>\n",
       "      <th>Plate_Map_Name</th>\n",
       "      <th>DNA</th>\n",
       "      <th>ER</th>\n",
       "      <th>RNA</th>\n",
       "      <th>AGP</th>\n",
       "      <th>Mito</th>\n",
       "      <th>broad_sample_Replicate</th>\n",
       "      <th>Treatment</th>\n",
       "      <th>Compound</th>\n",
       "      <th>Concentration</th>\n",
       "      <th>Replicate_Use</th>\n",
       "      <th>Cells</th>\n",
       "      <th>Var1</th>\n",
       "      <th>Metadata_moa.x</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>20585</td>\n",
       "      <td>A01</td>\n",
       "      <td>1</td>\n",
       "      <td>H-BIOA-002-1</td>\n",
       "      <td>20585/IXMtest_A01_s1_w164FBEEF7-F77C-4892-86F5...</td>\n",
       "      <td>20585/IXMtest_A01_s1_w3FD54227A-6BE4-460B-8E2C...</td>\n",
       "      <td>20585/IXMtest_A01_s1_w22CA7BDF3-1425-4AE6-8E57...</td>\n",
       "      <td>20585/IXMtest_A01_s1_w47457436C-5BE7-4DA2-B966...</td>\n",
       "      <td>20585/IXMtest_A01_s1_w58DA6AF62-1799-4B85-BE17...</td>\n",
       "      <td>1</td>\n",
       "      <td>BRD-K98763141-001-06-8@8.858401198</td>\n",
       "      <td>BRD-K98763141-001-06-8</td>\n",
       "      <td>8.858401</td>\n",
       "      <td>Training</td>\n",
       "      <td>54</td>\n",
       "      <td>BRD-K98763141-001-06-8</td>\n",
       "      <td>cyclooxygenase inhibitor</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>20585</td>\n",
       "      <td>A01</td>\n",
       "      <td>2</td>\n",
       "      <td>H-BIOA-002-1</td>\n",
       "      <td>20585/IXMtest_A01_s2_w1AA6B1894-F561-42EE-9D1D...</td>\n",
       "      <td>20585/IXMtest_A01_s2_w3A597237B-C3D7-43AE-8399...</td>\n",
       "      <td>20585/IXMtest_A01_s2_w246FFAEE1-BEB6-4C81-913B...</td>\n",
       "      <td>20585/IXMtest_A01_s2_w46657239A-5AFE-4C29-BB9B...</td>\n",
       "      <td>20585/IXMtest_A01_s2_w50F1562CD-EBCF-408E-9D8D...</td>\n",
       "      <td>1</td>\n",
       "      <td>BRD-K98763141-001-06-8@8.858401198</td>\n",
       "      <td>BRD-K98763141-001-06-8</td>\n",
       "      <td>8.858401</td>\n",
       "      <td>Training</td>\n",
       "      <td>87</td>\n",
       "      <td>BRD-K98763141-001-06-8</td>\n",
       "      <td>cyclooxygenase inhibitor</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>20585</td>\n",
       "      <td>A01</td>\n",
       "      <td>3</td>\n",
       "      <td>H-BIOA-002-1</td>\n",
       "      <td>20585/IXMtest_A01_s3_w108F3A6FF-A865-40A1-8488...</td>\n",
       "      <td>20585/IXMtest_A01_s3_w30B946780-C923-4A40-8EC5...</td>\n",
       "      <td>20585/IXMtest_A01_s3_w22AD5BD08-58A3-466D-A364...</td>\n",
       "      <td>20585/IXMtest_A01_s3_w47B5F2FF6-A5A2-4101-AC2A...</td>\n",
       "      <td>20585/IXMtest_A01_s3_w5489716D4-2D06-45BE-8D68...</td>\n",
       "      <td>1</td>\n",
       "      <td>BRD-K98763141-001-06-8@8.858401198</td>\n",
       "      <td>BRD-K98763141-001-06-8</td>\n",
       "      <td>8.858401</td>\n",
       "      <td>Training</td>\n",
       "      <td>68</td>\n",
       "      <td>BRD-K98763141-001-06-8</td>\n",
       "      <td>cyclooxygenase inhibitor</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>20585</td>\n",
       "      <td>A01</td>\n",
       "      <td>4</td>\n",
       "      <td>H-BIOA-002-1</td>\n",
       "      <td>20585/IXMtest_A01_s4_w189967984-D8CB-4449-9CE5...</td>\n",
       "      <td>20585/IXMtest_A01_s4_w3F582BADD-2249-48D6-BC19...</td>\n",
       "      <td>20585/IXMtest_A01_s4_w231D1EC99-22A8-47A8-9691...</td>\n",
       "      <td>20585/IXMtest_A01_s4_w47E34AB5F-23B4-42E1-8B5B...</td>\n",
       "      <td>20585/IXMtest_A01_s4_w5D0566314-B9B2-4299-98AC...</td>\n",
       "      <td>1</td>\n",
       "      <td>BRD-K98763141-001-06-8@8.858401198</td>\n",
       "      <td>BRD-K98763141-001-06-8</td>\n",
       "      <td>8.858401</td>\n",
       "      <td>Training</td>\n",
       "      <td>59</td>\n",
       "      <td>BRD-K98763141-001-06-8</td>\n",
       "      <td>cyclooxygenase inhibitor</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>20585</td>\n",
       "      <td>A01</td>\n",
       "      <td>5</td>\n",
       "      <td>H-BIOA-002-1</td>\n",
       "      <td>20585/IXMtest_A01_s5_w1DCC72C3D-6ABF-4D49-897E...</td>\n",
       "      <td>20585/IXMtest_A01_s5_w3EB7ABD45-EE05-41A0-979E...</td>\n",
       "      <td>20585/IXMtest_A01_s5_w2B3F07F02-876D-455B-94E4...</td>\n",
       "      <td>20585/IXMtest_A01_s5_w4AFA01300-7373-42C4-872F...</td>\n",
       "      <td>20585/IXMtest_A01_s5_w513201EB0-3500-4A0C-8FAC...</td>\n",
       "      <td>1</td>\n",
       "      <td>BRD-K98763141-001-06-8@8.858401198</td>\n",
       "      <td>BRD-K98763141-001-06-8</td>\n",
       "      <td>8.858401</td>\n",
       "      <td>Training</td>\n",
       "      <td>65</td>\n",
       "      <td>BRD-K98763141-001-06-8</td>\n",
       "      <td>cyclooxygenase inhibitor</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",
       "      <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>43758</th>\n",
       "      <td>20646</td>\n",
       "      <td>M02</td>\n",
       "      <td>5</td>\n",
       "      <td>H-BIOA-006-1</td>\n",
       "      <td>20646/IXMtest_M02_s5_w101C2C7A9-A498-464F-B018...</td>\n",
       "      <td>20646/IXMtest_M02_s5_w3D2C72899-1373-4194-B8D4...</td>\n",
       "      <td>20646/IXMtest_M02_s5_w29D1FDC1A-C092-4862-A792...</td>\n",
       "      <td>20646/IXMtest_M02_s5_w4CB980F20-9EF6-4E98-B1BB...</td>\n",
       "      <td>20646/IXMtest_M02_s5_w5B6C1B344-D5CB-4B43-972D...</td>\n",
       "      <td>4</td>\n",
       "      <td>BRD-A28970875-300-02-3@5</td>\n",
       "      <td>BRD-A28970875-300-02-3</td>\n",
       "      <td>5.000000</td>\n",
       "      <td>Training</td>\n",
       "      <td>10</td>\n",
       "      <td>BRD-A28970875-300-02-3</td>\n",
       "      <td>protein synthesis inhibitor</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>43759</th>\n",
       "      <td>20646</td>\n",
       "      <td>M02</td>\n",
       "      <td>6</td>\n",
       "      <td>H-BIOA-006-1</td>\n",
       "      <td>20646/IXMtest_M02_s6_w105CAF1F0-CA39-4087-9B02...</td>\n",
       "      <td>20646/IXMtest_M02_s6_w3C278E8AD-280F-4C4B-9250...</td>\n",
       "      <td>20646/IXMtest_M02_s6_w20EAA76C6-97A5-4E8F-9172...</td>\n",
       "      <td>20646/IXMtest_M02_s6_w48350A6C8-F059-4449-9324...</td>\n",
       "      <td>20646/IXMtest_M02_s6_w5B8591209-F255-4CA9-A1C9...</td>\n",
       "      <td>4</td>\n",
       "      <td>BRD-A28970875-300-02-3@5</td>\n",
       "      <td>BRD-A28970875-300-02-3</td>\n",
       "      <td>5.000000</td>\n",
       "      <td>Training</td>\n",
       "      <td>15</td>\n",
       "      <td>BRD-A28970875-300-02-3</td>\n",
       "      <td>protein synthesis inhibitor</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>43760</th>\n",
       "      <td>20646</td>\n",
       "      <td>M02</td>\n",
       "      <td>7</td>\n",
       "      <td>H-BIOA-006-1</td>\n",
       "      <td>20646/IXMtest_M02_s7_w1BC4B5401-C5BA-4A34-8779...</td>\n",
       "      <td>20646/IXMtest_M02_s7_w38E3FCCB1-150E-4298-A258...</td>\n",
       "      <td>20646/IXMtest_M02_s7_w2856A96A2-AED3-4692-A050...</td>\n",
       "      <td>20646/IXMtest_M02_s7_w4F7E71E60-F806-43D6-A338...</td>\n",
       "      <td>20646/IXMtest_M02_s7_w56505AF21-F372-4C13-AFE0...</td>\n",
       "      <td>4</td>\n",
       "      <td>BRD-A28970875-300-02-3@5</td>\n",
       "      <td>BRD-A28970875-300-02-3</td>\n",
       "      <td>5.000000</td>\n",
       "      <td>Training</td>\n",
       "      <td>4</td>\n",
       "      <td>BRD-A28970875-300-02-3</td>\n",
       "      <td>protein synthesis inhibitor</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>43761</th>\n",
       "      <td>20646</td>\n",
       "      <td>M02</td>\n",
       "      <td>8</td>\n",
       "      <td>H-BIOA-006-1</td>\n",
       "      <td>20646/IXMtest_M02_s8_w1B0FC866D-211E-43EC-B1F9...</td>\n",
       "      <td>20646/IXMtest_M02_s8_w39B4556A4-AB25-4047-94BB...</td>\n",
       "      <td>20646/IXMtest_M02_s8_w2514D5EAF-8D8B-4E88-8BEE...</td>\n",
       "      <td>20646/IXMtest_M02_s8_w4BF8EDC77-D1C2-4160-9B2D...</td>\n",
       "      <td>20646/IXMtest_M02_s8_w549BDA3B3-959A-4881-8E24...</td>\n",
       "      <td>4</td>\n",
       "      <td>BRD-A28970875-300-02-3@5</td>\n",
       "      <td>BRD-A28970875-300-02-3</td>\n",
       "      <td>5.000000</td>\n",
       "      <td>Training</td>\n",
       "      <td>4</td>\n",
       "      <td>BRD-A28970875-300-02-3</td>\n",
       "      <td>protein synthesis inhibitor</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>43762</th>\n",
       "      <td>20646</td>\n",
       "      <td>M02</td>\n",
       "      <td>9</td>\n",
       "      <td>H-BIOA-006-1</td>\n",
       "      <td>20646/IXMtest_M02_s9_w18A5C2635-593B-45AE-9000...</td>\n",
       "      <td>20646/IXMtest_M02_s9_w37C024823-09AD-4A7F-8881...</td>\n",
       "      <td>20646/IXMtest_M02_s9_w293E08447-1672-41B2-949E...</td>\n",
       "      <td>20646/IXMtest_M02_s9_w488F09CF1-3664-4621-B749...</td>\n",
       "      <td>20646/IXMtest_M02_s9_w5EA7793F2-2FAC-478E-B02B...</td>\n",
       "      <td>4</td>\n",
       "      <td>BRD-A28970875-300-02-3@5</td>\n",
       "      <td>BRD-A28970875-300-02-3</td>\n",
       "      <td>5.000000</td>\n",
       "      <td>Training</td>\n",
       "      <td>11</td>\n",
       "      <td>BRD-A28970875-300-02-3</td>\n",
       "      <td>protein synthesis inhibitor</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>43763 rows × 17 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "       Metadata_Plate Metadata_Well  Metadata_Site Plate_Map_Name  \\\n",
       "0               20585           A01              1   H-BIOA-002-1   \n",
       "1               20585           A01              2   H-BIOA-002-1   \n",
       "2               20585           A01              3   H-BIOA-002-1   \n",
       "3               20585           A01              4   H-BIOA-002-1   \n",
       "4               20585           A01              5   H-BIOA-002-1   \n",
       "...               ...           ...            ...            ...   \n",
       "43758           20646           M02              5   H-BIOA-006-1   \n",
       "43759           20646           M02              6   H-BIOA-006-1   \n",
       "43760           20646           M02              7   H-BIOA-006-1   \n",
       "43761           20646           M02              8   H-BIOA-006-1   \n",
       "43762           20646           M02              9   H-BIOA-006-1   \n",
       "\n",
       "                                                     DNA  \\\n",
       "0      20585/IXMtest_A01_s1_w164FBEEF7-F77C-4892-86F5...   \n",
       "1      20585/IXMtest_A01_s2_w1AA6B1894-F561-42EE-9D1D...   \n",
       "2      20585/IXMtest_A01_s3_w108F3A6FF-A865-40A1-8488...   \n",
       "3      20585/IXMtest_A01_s4_w189967984-D8CB-4449-9CE5...   \n",
       "4      20585/IXMtest_A01_s5_w1DCC72C3D-6ABF-4D49-897E...   \n",
       "...                                                  ...   \n",
       "43758  20646/IXMtest_M02_s5_w101C2C7A9-A498-464F-B018...   \n",
       "43759  20646/IXMtest_M02_s6_w105CAF1F0-CA39-4087-9B02...   \n",
       "43760  20646/IXMtest_M02_s7_w1BC4B5401-C5BA-4A34-8779...   \n",
       "43761  20646/IXMtest_M02_s8_w1B0FC866D-211E-43EC-B1F9...   \n",
       "43762  20646/IXMtest_M02_s9_w18A5C2635-593B-45AE-9000...   \n",
       "\n",
       "                                                      ER  \\\n",
       "0      20585/IXMtest_A01_s1_w3FD54227A-6BE4-460B-8E2C...   \n",
       "1      20585/IXMtest_A01_s2_w3A597237B-C3D7-43AE-8399...   \n",
       "2      20585/IXMtest_A01_s3_w30B946780-C923-4A40-8EC5...   \n",
       "3      20585/IXMtest_A01_s4_w3F582BADD-2249-48D6-BC19...   \n",
       "4      20585/IXMtest_A01_s5_w3EB7ABD45-EE05-41A0-979E...   \n",
       "...                                                  ...   \n",
       "43758  20646/IXMtest_M02_s5_w3D2C72899-1373-4194-B8D4...   \n",
       "43759  20646/IXMtest_M02_s6_w3C278E8AD-280F-4C4B-9250...   \n",
       "43760  20646/IXMtest_M02_s7_w38E3FCCB1-150E-4298-A258...   \n",
       "43761  20646/IXMtest_M02_s8_w39B4556A4-AB25-4047-94BB...   \n",
       "43762  20646/IXMtest_M02_s9_w37C024823-09AD-4A7F-8881...   \n",
       "\n",
       "                                                     RNA  \\\n",
       "0      20585/IXMtest_A01_s1_w22CA7BDF3-1425-4AE6-8E57...   \n",
       "1      20585/IXMtest_A01_s2_w246FFAEE1-BEB6-4C81-913B...   \n",
       "2      20585/IXMtest_A01_s3_w22AD5BD08-58A3-466D-A364...   \n",
       "3      20585/IXMtest_A01_s4_w231D1EC99-22A8-47A8-9691...   \n",
       "4      20585/IXMtest_A01_s5_w2B3F07F02-876D-455B-94E4...   \n",
       "...                                                  ...   \n",
       "43758  20646/IXMtest_M02_s5_w29D1FDC1A-C092-4862-A792...   \n",
       "43759  20646/IXMtest_M02_s6_w20EAA76C6-97A5-4E8F-9172...   \n",
       "43760  20646/IXMtest_M02_s7_w2856A96A2-AED3-4692-A050...   \n",
       "43761  20646/IXMtest_M02_s8_w2514D5EAF-8D8B-4E88-8BEE...   \n",
       "43762  20646/IXMtest_M02_s9_w293E08447-1672-41B2-949E...   \n",
       "\n",
       "                                                     AGP  \\\n",
       "0      20585/IXMtest_A01_s1_w47457436C-5BE7-4DA2-B966...   \n",
       "1      20585/IXMtest_A01_s2_w46657239A-5AFE-4C29-BB9B...   \n",
       "2      20585/IXMtest_A01_s3_w47B5F2FF6-A5A2-4101-AC2A...   \n",
       "3      20585/IXMtest_A01_s4_w47E34AB5F-23B4-42E1-8B5B...   \n",
       "4      20585/IXMtest_A01_s5_w4AFA01300-7373-42C4-872F...   \n",
       "...                                                  ...   \n",
       "43758  20646/IXMtest_M02_s5_w4CB980F20-9EF6-4E98-B1BB...   \n",
       "43759  20646/IXMtest_M02_s6_w48350A6C8-F059-4449-9324...   \n",
       "43760  20646/IXMtest_M02_s7_w4F7E71E60-F806-43D6-A338...   \n",
       "43761  20646/IXMtest_M02_s8_w4BF8EDC77-D1C2-4160-9B2D...   \n",
       "43762  20646/IXMtest_M02_s9_w488F09CF1-3664-4621-B749...   \n",
       "\n",
       "                                                    Mito  \\\n",
       "0      20585/IXMtest_A01_s1_w58DA6AF62-1799-4B85-BE17...   \n",
       "1      20585/IXMtest_A01_s2_w50F1562CD-EBCF-408E-9D8D...   \n",
       "2      20585/IXMtest_A01_s3_w5489716D4-2D06-45BE-8D68...   \n",
       "3      20585/IXMtest_A01_s4_w5D0566314-B9B2-4299-98AC...   \n",
       "4      20585/IXMtest_A01_s5_w513201EB0-3500-4A0C-8FAC...   \n",
       "...                                                  ...   \n",
       "43758  20646/IXMtest_M02_s5_w5B6C1B344-D5CB-4B43-972D...   \n",
       "43759  20646/IXMtest_M02_s6_w5B8591209-F255-4CA9-A1C9...   \n",
       "43760  20646/IXMtest_M02_s7_w56505AF21-F372-4C13-AFE0...   \n",
       "43761  20646/IXMtest_M02_s8_w549BDA3B3-959A-4881-8E24...   \n",
       "43762  20646/IXMtest_M02_s9_w5EA7793F2-2FAC-478E-B02B...   \n",
       "\n",
       "       broad_sample_Replicate                           Treatment  \\\n",
       "0                           1  BRD-K98763141-001-06-8@8.858401198   \n",
       "1                           1  BRD-K98763141-001-06-8@8.858401198   \n",
       "2                           1  BRD-K98763141-001-06-8@8.858401198   \n",
       "3                           1  BRD-K98763141-001-06-8@8.858401198   \n",
       "4                           1  BRD-K98763141-001-06-8@8.858401198   \n",
       "...                       ...                                 ...   \n",
       "43758                       4            BRD-A28970875-300-02-3@5   \n",
       "43759                       4            BRD-A28970875-300-02-3@5   \n",
       "43760                       4            BRD-A28970875-300-02-3@5   \n",
       "43761                       4            BRD-A28970875-300-02-3@5   \n",
       "43762                       4            BRD-A28970875-300-02-3@5   \n",
       "\n",
       "                     Compound  Concentration Replicate_Use  Cells  \\\n",
       "0      BRD-K98763141-001-06-8       8.858401      Training     54   \n",
       "1      BRD-K98763141-001-06-8       8.858401      Training     87   \n",
       "2      BRD-K98763141-001-06-8       8.858401      Training     68   \n",
       "3      BRD-K98763141-001-06-8       8.858401      Training     59   \n",
       "4      BRD-K98763141-001-06-8       8.858401      Training     65   \n",
       "...                       ...            ...           ...    ...   \n",
       "43758  BRD-A28970875-300-02-3       5.000000      Training     10   \n",
       "43759  BRD-A28970875-300-02-3       5.000000      Training     15   \n",
       "43760  BRD-A28970875-300-02-3       5.000000      Training      4   \n",
       "43761  BRD-A28970875-300-02-3       5.000000      Training      4   \n",
       "43762  BRD-A28970875-300-02-3       5.000000      Training     11   \n",
       "\n",
       "                         Var1               Metadata_moa.x  \n",
       "0      BRD-K98763141-001-06-8     cyclooxygenase inhibitor  \n",
       "1      BRD-K98763141-001-06-8     cyclooxygenase inhibitor  \n",
       "2      BRD-K98763141-001-06-8     cyclooxygenase inhibitor  \n",
       "3      BRD-K98763141-001-06-8     cyclooxygenase inhibitor  \n",
       "4      BRD-K98763141-001-06-8     cyclooxygenase inhibitor  \n",
       "...                       ...                          ...  \n",
       "43758  BRD-A28970875-300-02-3  protein synthesis inhibitor  \n",
       "43759  BRD-A28970875-300-02-3  protein synthesis inhibitor  \n",
       "43760  BRD-A28970875-300-02-3  protein synthesis inhibitor  \n",
       "43761  BRD-A28970875-300-02-3  protein synthesis inhibitor  \n",
       "43762  BRD-A28970875-300-02-3  protein synthesis inhibitor  \n",
       "\n",
       "[43763 rows x 17 columns]"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "bbbc022_moa_matches = pd.read_csv(\"./bbbc022/data/BBBC022_MOA_MATCHES_fixed.csv\")\n",
    "bbbc022_moa_matches = bbbc022_moa_matches.assign(var1=bbbc022_moa_matches['Metadata_moa.x'].str.split('|')).explode('var1')\n",
    "bbbc022_moa_matches.reset_index(inplace=True, drop = True)\n",
    "bbbc022_moa_matches.drop(columns=['Unnamed: 0', 'Unnamed: 0.1', 'Metadata_moa.x'], inplace= True)\n",
    "bbbc022_moa_matches.rename(columns={\"var1\": \"Metadata_moa.x\"}, inplace=True)\n",
    "\n",
    "bbbc022_meta_moa = pd.merge(bbbc022_meta, bbbc022_moa_matches[['Var1', 'Metadata_moa.x']], right_on=['Var1'], left_on =['Compound'])\n",
    "bbbc022_meta_moa"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Calculate wells/sites/cells per MoA"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(array([221.,  71.,  32.,  16.,  15.,   2.,   5.,   2.,   6.,   3.,   1.,\n",
       "          0.,   0.,   2.,   0.,   0.,   1.,   0.,   1.,   1.,   0.,   0.,\n",
       "          1.,   0.,   2.,   1.,   0.,   1.,   1.,   0.,   0.,   0.,   3.,\n",
       "          0.,   0.,   0.,   0.,   0.,   1.,   0.,   0.,   0.,   0.,   0.,\n",
       "          0.,   0.,   0.,   0.,   0.,   1.]),\n",
       " array([  1.  ,   5.34,   9.68,  14.02,  18.36,  22.7 ,  27.04,  31.38,\n",
       "         35.72,  40.06,  44.4 ,  48.74,  53.08,  57.42,  61.76,  66.1 ,\n",
       "         70.44,  74.78,  79.12,  83.46,  87.8 ,  92.14,  96.48, 100.82,\n",
       "        105.16, 109.5 , 113.84, 118.18, 122.52, 126.86, 131.2 , 135.54,\n",
       "        139.88, 144.22, 148.56, 152.9 , 157.24, 161.58, 165.92, 170.26,\n",
       "        174.6 , 178.94, 183.28, 187.62, 191.96, 196.3 , 200.64, 204.98,\n",
       "        209.32, 213.66, 218.  ]),\n",
       " <a list of 50 Patch objects>)"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD4CAYAAADvsV2wAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAOzUlEQVR4nO3dXYycZ3nG8f/VGDjgo3Fq17KclA3IRUoPGtIVjVRAVGlDEloc2ioKqohLU7mVEgnUVpUpUuEEKbSCqqhtkCkRTsVHgiCKJWghtVBRDwJsUhPyQYgTHMWWYxuCCBIVbcLdg3mWDGbt3Z3ZmdnN8/9Jo3nnnndm7nlm5tp3nnlnNlWFJOn57+dm3YAkaToMfEnqhIEvSZ0w8CWpEwa+JHVi06wbANiyZUvNzc3Nug1J2lDuueee71TV1pWuvy4Cf25ujoWFhVm3IUkbSpLHV7O+UzqS1AkDX5I6YeBLUicMfEnqhIEvSZ0w8CWpEwa+JHXCwJekThj4ktSJdfFN23HM7f3cGc87ctObptiJJK1vbuFLUicMfEnqhIEvSZ0w8CWpEwa+JHXCwJekThj4ktQJA1+SOmHgS1InDHxJ6oSBL0mdMPAlqRMGviR1wsCXpE4Y+JLUiWUDP8kFSb6U5MEkDyR5R6ufl+SuJI+0482tniQfSnI4yX1JLpn0nZAkLW8lW/jPAH9RVRcBlwI3JLkI2AscrKqdwMF2GuBKYGc77AFuXvOuJUmrtmzgV9Xxqrq3Lf8AeAjYAewC9rfV9gNXt+VdwK01cDdwbpLta965JGlVVjWHn2QOeDXwFWBbVR1vZz0JbGvLO4Anhi52tNVOv649SRaSLJw6dWqVbUuSVmvFgZ/kJcBngHdW1dPD51VVAbWaG66qfVU1X1XzW7duXc1FJUkjWFHgJ3kBg7D/eFV9tpVPLE7VtOOTrX4MuGDo4ue3miRphlayl06AjwIPVdUHh846AOxuy7uBO4fq17W9dS4Fvj809SNJmpFNK1jnN4C3Ad9IcqjV/hq4Cbg9yfXA48A17bzPA1cBh4EfAm9f044lSSNZNvCr6r+AnOHsy5ZYv4AbxuxLkrTG/KatJHXCwJekThj4ktQJA1+SOmHgS1InDHxJ6oSBL0mdMPAlqRMGviR1wsCXpE4Y+JLUCQNfkjph4EtSJwx8SeqEgS9JnTDwJakTBr4kdcLAl6ROGPiS1AkDX5I6YeBLUicMfEnqhIEvSZ0w8CWpEwa+JHXCwJekThj4ktQJA1+SOmHgS1InDHxJ6oSBL0mdMPAlqRMGviR1wsCXpE4Y+JLUCQNfkjph4EtSJwx8SerEsoGf5JYkJ5PcP1R7b5JjSQ61w1VD570ryeEkDyd546QalyStzkq28D8GXLFE/e+r6uJ2+DxAkouAa4FfaZf55yTnrFWzkqTRLRv4VfVl4KkVXt8u4FNV9aOq+jZwGHjNGP1JktbIOHP4Nya5r035bG61HcATQ+scbbWfkWRPkoUkC6dOnRqjDUnSSowa+DcDrwQuBo4DH1jtFVTVvqqar6r5rVu3jtiGJGmlRgr8qjpRVc9W1Y+Bj/DctM0x4IKhVc9vNUnSjI0U+Em2D518C7C4B88B4NokL0pyIbAT+Op4LUqS1sKm5VZI8kngDcCWJEeB9wBvSHIxUMAR4E8BquqBJLcDDwLPADdU1bOTaV2StBrLBn5VvXWJ8kfPsv77gPeN05Qkae35TVtJ6oSBL0mdMPAlqRMGviR1wsCXpE4Y+JLUCQNfkjph4EtSJwx8SeqEgS9JnTDwJakTBr4kdcLAl6ROGPiS1AkDX5I6YeBLUicMfEnqhIEvSZ0w8CWpEwa+JHXCwJekThj4ktQJA1+SOmHgS1InDHxJ6oSBL0mdMPAlqRMGviR1wsCXpE4Y+JLUCQNfkjph4EtSJwx8SeqEgS9JnTDwJakTBr4kdcLAl6ROGPiS1IllAz/JLUlOJrl/qHZekruSPNKON7d6knwoyeEk9yW5ZJLNS5JWbiVb+B8Drjitthc4WFU7gYPtNMCVwM522APcvDZtSpLGtWzgV9WXgadOK+8C9rfl/cDVQ/Vba+Bu4Nwk29eqWUnS6Eadw99WVcfb8pPAtra8A3hiaL2jrfYzkuxJspBk4dSpUyO2IUlaqbE/tK2qAmqEy+2rqvmqmt+6deu4bUiSljFq4J9YnKppxydb/RhwwdB657eaJGnGRg38A8DutrwbuHOofl3bW+dS4PtDUz+SpBnatNwKST4JvAHYkuQo8B7gJuD2JNcDjwPXtNU/D1wFHAZ+CLx9Aj1LkkawbOBX1VvPcNZlS6xbwA3jNiVJWnt+01aSOmHgS1InDHxJ6oSBL0mdMPAlqRMGviR1wsCXpE4Y+JLUCQNfkjph4EtSJwx8SeqEgS9JnTDwJakTBr4kdcLAl6ROGPiS1AkDX5I6YeBLUicMfEnqhIEvSZ0w8CWpE5tm3cAkze393JL1Ize9acqdSNLsuYUvSZ0w8CWpEwa+JHXCwJekThj4ktQJA1+SOmHgS1InDHxJ6oSBL0mdMPAlqRMGviR1wsCXpE4Y+JLUCQNfkjph4EtSJwx8SerEWP8AJckR4AfAs8AzVTWf5DzgNmAOOAJcU1XfG69NSdK41mIL/zer6uKqmm+n9wIHq2oncLCdliTN2CSmdHYB+9vyfuDqCdyGJGmVxg38Ar6Y5J4ke1ptW1Udb8tPAtuWumCSPUkWkiycOnVqzDYkScsZ95+Yv7aqjiX5ReCuJN8cPrOqKkktdcGq2gfsA5ifn19yHUnS2hlrC7+qjrXjk8AdwGuAE0m2A7Tjk+M2KUka38iBn+TFSV66uAxcDtwPHAB2t9V2A3eO26QkaXzjTOlsA+5Isng9n6iqf0/yNeD2JNcDjwPXjN+mJGlcIwd+VT0G/OoS9e8Cl43TlCRp7flNW0nqhIEvSZ0w8CWpEwa+JHXCwJekThj4ktQJA1+SOmHgS1InDHxJ6oSBL0mdMPAlqRPj/h7+hjS393NL1o/c9KYpdyJJ0+MWviR1wsCXpE4Y+JLUCQNfkjph4EtSJwx8SeqEgS9JnTDwJakTBr4kdcLAl6ROGPiS1AkDX5I6YeBLUicMfEnqRJc/j3wm/myypOczA38F/EMg6fnAKR1J6oSBL0mdMPAlqRMGviR1wsCXpE4Y+JLUCQNfkjrhfvjrgPv5S5oGA38CDHBJ65FTOpLUiYlt4Se5AvgH4BzgX6rqpknd1kZxpi3/tboe30H0zeeFljORwE9yDvBPwG8DR4GvJTlQVQ9O4vY0OkPiOY6Fnu8mtYX/GuBwVT0GkORTwC7geRX4a7XFvlZmGVirve3Vjt1q78NajsV6u28biX9En7MexiJVtfZXmvwBcEVV/Uk7/Tbg16vqxqF19gB72slXAQ+PeHNbgO+M0e4k2dto7G009jaajdzby6tq60qvbGZ76VTVPmDfuNeTZKGq5tegpTVnb6Oxt9HY22h66m1Se+kcAy4YOn1+q0mSZmRSgf81YGeSC5O8ELgWODCh25IkrcBEpnSq6pkkNwJfYLBb5i1V9cAkbos1mBaaIHsbjb2Nxt5G001vE/nQVpK0/vhNW0nqhIEvSZ3YsIGf5IokDyc5nGTvDG7/giRfSvJgkgeSvKPV35vkWJJD7XDV0GXe1fp9OMkbJ9zfkSTfaD0stNp5Se5K8kg73tzqSfKh1tt9SS6ZYF+vGhqbQ0meTvLOWY5bkluSnExy/1Bt1WOVZHdb/5EkuyfY298l+Wa7/TuSnNvqc0n+Z2gMPzx0mV9rz4fDrf9MqLdVP46TeC2fobfbhvo6kuRQq09t3M6SG9N5vlXVhjsw+CD4UeAVwAuBrwMXTbmH7cAlbfmlwLeAi4D3An+5xPoXtT5fBFzY+j9ngv0dAbacVvtbYG9b3gu8vy1fBfwbEOBS4CtTfByfBF4+y3EDXg9cAtw/6lgB5wGPtePNbXnzhHq7HNjUlt8/1Nvc8HqnXc9XW79p/V85od5W9ThO6rW8VG+nnf8B4G+mPW5nyY2pPN826hb+T366oar+F1j86YapqarjVXVvW/4B8BCw4ywX2QV8qqp+VFXfBg4zuB/TtAvY35b3A1cP1W+tgbuBc5Nsn0I/lwGPVtXjZ1ln4uNWVV8GnlridlczVm8E7qqqp6rqe8BdwBWT6K2qvlhVz7STdzP4nssZtf5eVlV31yAtbh26P2va21mc6XGcyGv5bL21rfRrgE+e7TomMW5nyY2pPN82auDvAJ4YOn2Us4ftRCWZA14NfKWVbmxvv25ZfGvG9Hsu4ItJ7sngZywAtlXV8bb8JLBtRr0tupafftGth3FbtNqxmlWff8xgC3DRhUn+O8l/Jnldq+1o/Uyrt9U8jrMYt9cBJ6rqkaHa1MfttNyYyvNtowb+upHkJcBngHdW1dPAzcArgYuB4wzeOs7Ca6vqEuBK4IYkrx8+s22xzGyf3Ay+kPdm4NOttF7G7WfMeqzOJMm7gWeAj7fSceCXqurVwJ8Dn0jysim3tW4fxyFv5ac3NKY+bkvkxk9M8vm2UQN/Xfx0Q5IXMHjQPl5VnwWoqhNV9WxV/Rj4CM9NP0y156o61o5PAne0Pk4sTtW045Oz6K25Eri3qk60PtfFuA1Z7VhNtc8kfwT8DvCHLSBo0yXfbcv3MJgb/+XWx/C0z8R6G+FxnPa4bQJ+D7htqOepjttSucGUnm8bNfBn/tMNbR7wo8BDVfXBofrw3PdbgMW9BA4A1yZ5UZILgZ0MPhCaRG8vTvLSxWUGH/Ld33pY/DR/N3DnUG/XtT0CLgW+P/T2clJ+aitrPYzbaVY7Vl8ALk+yuU1jXN5qay6Dfy70V8Cbq+qHQ/WtGfwvCpK8gsFYPdb6ezrJpe15e93Q/Vnr3lb7OE77tfxbwDer6idTNdMctzPlBtN6vo3zifMsDww+vf4Wg7/G757B7b+Wwduu+4BD7XAV8K/AN1r9ALB96DLvbv0+zBrsJXGW3l7BYG+HrwMPLI4P8AvAQeAR4D+A81o9DP5hzaOt9/kJj92Lge8CPz9Um9m4MfjDcxz4PwZzodePMlYM5tMPt8PbJ9jbYQbzt4vPuw+3dX+/Pd6HgHuB3x26nnkG4fso8I+0b9lPoLdVP46TeC0v1Vurfwz4s9PWndq4cebcmMrzzZ9WkKRObNQpHUnSKhn4ktQJA1+SOmHgS1InDHxJ6oSBL0mdMPAlqRP/D3/91Y7W/8PVAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAMxElEQVR4nO3dT4ic933H8fendppDErBdqULIonKDKCiHKmZxDTXFwTSR5YMcKMY+xCK4qAcZEshFaQ7JJaBLEmpoDEptLJfUriExFti0cUXA9ODE6+DK/+paSWQsIUubujiGQFo73x72ER7Lu17tzs6O9qv3C4aZ+T3P7P72x7NvzT47O0pVIUnq5femPQFJ0uoz7pLUkHGXpIaMuyQ1ZNwlqaHLpz0BgA0bNtS2bdumPQ1JWleeffbZX1XVxoW2XRRx37ZtG7Ozs9OehiStK0leW2ybp2UkqSHjLkkNGXdJasi4S1JDxl2SGjLuktSQcZekhoy7JDVk3CWpoYviL1THse3A44tuO3HwljWciSRdPHzmLkkNGXdJasi4S1JDxl2SGjLuktSQcZekhoy7JDVk3CWpIeMuSQ0Zd0lqyLhLUkPGXZIaMu6S1JBxl6SGjLskNbRk3JNsTfLjJC8leTHJl4bxq5I8meTV4frKYTxJ7klyPMmxJNdO+ouQJL3fhTxzfwf4SlXtAK4H9ifZARwAjlbVduDocB/gZmD7cNkH3Lvqs5Ykfagl415Vp6vqZ8Ptt4GXgS3AHuDwsNth4Nbh9h7gwZr3NHBFks2rPnNJ0qKWdc49yTbg08BPgE1VdXrY9Aawabi9BXh95GEnh7HzP9a+JLNJZufm5pY5bUnSh7nguCf5OPAD4MtV9evRbVVVQC3nE1fVoaqaqaqZjRs3LuehkqQlXFDck3yE+bB/v6p+OAyfOXe6Zbg+O4yfAraOPPzqYUyStEYu5NUyAe4DXq6qb49sOgLsHW7vBR4bGb9zeNXM9cBbI6dvJElr4PIL2OfPgS8Azyd5bhj7W+Ag8EiSu4DXgNuGbU8Au4HjwG+AL67qjCVJS1oy7lX170AW2XzTAvsXsH/MeUmSxuBfqEpSQ8Zdkhoy7pLUkHGXpIaMuyQ1ZNwlqSHjLkkNGXdJasi4S1JDxl2SGjLuktSQcZekhoy7JDVk3CWpIeMuSQ0Zd0lqyLhLUkPGXZIaMu6S1JBxl6SGjLskNWTcJakh4y5JDRl3SWrIuEtSQ8Zdkhoy7pLUkHGXpIaMuyQ1ZNwlqSHjLkkNGXdJasi4S1JDxl2SGjLuktSQcZekhoy7JDVk3CWpoSXjnuT+JGeTvDAy9o0kp5I8N1x2j2z7apLjSV5J8rlJTVyStLgLeeb+ALBrgfHvVNXO4fIEQJIdwO3Ap4bHfDfJZas1WUnShVky7lX1FPDmBX68PcDDVfXbqvolcBy4boz5SZJWYJxz7ncnOTactrlyGNsCvD6yz8lh7AOS7Esym2R2bm5ujGlIks630rjfC3wS2AmcBr613A9QVYeqaqaqZjZu3LjCaUiSFrKiuFfVmap6t6p+B3yP9069nAK2jux69TAmSVpDK4p7ks0jdz8PnHslzRHg9iQfTXINsB346XhTlCQt1+VL7ZDkIeBGYEOSk8DXgRuT7AQKOAH8DUBVvZjkEeAl4B1gf1W9O5mpS5IWs2Tcq+qOBYbv+5D9vwl8c5xJSZLG41+oSlJDxl2SGjLuktSQcZekhoy7JDVk3CWpIeMuSQ0Zd0lqyLhLUkPGXZIaMu6S1JBxl6SGjLskNWTcJakh4y5JDRl3SWrIuEtSQ8Zdkhoy7pLUkHGXpIaMuyQ1ZNwlqSHjLkkNGXdJasi4S1JDxl2SGjLuktSQcZekhoy7JDVk3CWpIeMuSQ0Zd0lqyLhLUkPGXZIaMu6S1JBxl6SGjLskNWTcJamhJeOe5P4kZ5O8MDJ2VZInk7w6XF85jCfJPUmOJzmW5NpJTl6StLALeeb+ALDrvLEDwNGq2g4cHe4D3AxsHy77gHtXZ5qSpOVYMu5V9RTw5nnDe4DDw+3DwK0j4w/WvKeBK5JsXq3JSpIuzErPuW+qqtPD7TeATcPtLcDrI/udHMY+IMm+JLNJZufm5lY4DUnSQsb+hWpVFVAreNyhqpqpqpmNGzeOOw1J0oiVxv3MudMtw/XZYfwUsHVkv6uHMUnSGlpp3I8Ae4fbe4HHRsbvHF41cz3w1sjpG0nSGrl8qR2SPATcCGxIchL4OnAQeCTJXcBrwG3D7k8Au4HjwG+AL05gzpKkJSwZ96q6Y5FNNy2wbwH7x52UJGk8/oWqJDVk3CWpIeMuSQ0Zd0lqyLhLUkPGXZIaMu6S1JBxl6SGjLskNWTcJakh4y5JDRl3SWrIuEtSQ8Zdkhoy7pLUkHGXpIaMuyQ1ZNwlqSHjLkkNGXdJasi4S1JDl097ApO07cDjC46fOHjLGs9EktaWz9wlqSHjLkkNGXdJasi4S1JDxl2SGjLuktSQcZekhoy7JDVk3CWpIeMuSQ0Zd0lqyLhLUkPGXZIaMu6S1JBxl6SGjLskNTTWf9aR5ATwNvAu8E5VzSS5CvhnYBtwAritqv5nvGlKkpZjNZ65f6aqdlbVzHD/AHC0qrYDR4f7kqQ1NInTMnuAw8Ptw8CtE/gckqQPMW7cC/hRkmeT7BvGNlXV6eH2G8CmhR6YZF+S2SSzc3NzY05DkjRq3P8g+4aqOpXkD4Enk/zn6MaqqiS10AOr6hBwCGBmZmbBfSRJKzPWM/eqOjVcnwUeBa4DziTZDDBcnx13kpKk5Vlx3JN8LMknzt0GPgu8ABwB9g677QUeG3eSkqTlGee0zCbg0STnPs4/VdW/JHkGeCTJXcBrwG3jT1OStBwrjntV/QL40wXG/xu4aZxJSZLG41+oSlJDxl2SGjLuktSQcZekhoy7JDVk3CWpIeMuSQ0Zd0lqyLhLUkPGXZIaMu6S1NC47+e+Lm078PiC4ycO3rLGM5GkyfCZuyQ1ZNwlqSHjLkkNGXdJasi4S1JDxl2SGjLuktSQcZekhoy7JDVk3CWpIeMuSQ0Zd0lqyLhLUkPGXZIauiTf8ncxvhWwpC6M+wUw+pLWG0/LSFJDxl2SGjLuktSQcZekhoy7JDVk3CWpIeMuSQ35OveLgK+jl7TajPsEGGtJ0+ZpGUlqaGLP3JPsAv4OuAz4h6o6OKnP1ZU/AWg5PF40aiJxT3IZ8PfAXwIngWeSHKmqlybx+aZlsW+m1dp/LRiE97gW6mRSz9yvA45X1S8AkjwM7AFaxX1aLsYILXdOy/2Hbrlf22qu0cX2ta0nF+OxOi1rvRapqtX/oMlfAbuq6q+H+18A/qyq7h7ZZx+wb7j7J8ArK/hUG4BfjTndjlyXhbkuC3NdFnexr80fVdXGhTZM7dUyVXUIODTOx0gyW1UzqzSlNlyXhbkuC3NdFree12ZSr5Y5BWwduX/1MCZJWgOTivszwPYk1yT5feB24MiEPpck6TwTOS1TVe8kuRv4V+ZfCnl/Vb04gU811mmdxlyXhbkuC3NdFrdu12Yiv1CVJE2Xf6EqSQ0Zd0lqaN3GPcmuJK8kOZ7kwLTnM01JTiR5PslzSWaHsauSPJnk1eH6ymnPc9KS3J/kbJIXRsYWXIfMu2c4fo4luXZ6M5+sRdblG0lODcfMc0l2j2z76rAuryT53HRmPXlJtib5cZKXkryY5EvDeItjZl3GfeTtDW4GdgB3JNkx3VlN3WeqaufIa3IPAEerajtwdLjf3QPArvPGFluHm4Htw2UfcO8azXEaHuCD6wLwneGY2VlVTwAM30e3A58aHvPd4futo3eAr1TVDuB6YP/w9bc4ZtZl3Bl5e4Oq+l/g3Nsb6D17gMPD7cPArVOcy5qoqqeAN88bXmwd9gAP1ryngSuSbF6bma6tRdZlMXuAh6vqt1X1S+A4899v7VTV6ar62XD7beBlYAtNjpn1GvctwOsj908OY5eqAn6U5NnhbR0ANlXV6eH2G8Cm6Uxt6hZbB48huHs4vXD/yGm7S3JdkmwDPg38hCbHzHqNu97vhqq6lvkfG/cn+YvRjTX/etdL/jWvrsP73At8EtgJnAa+Nd3pTE+SjwM/AL5cVb8e3baej5n1Gnff3mBEVZ0ars8CjzL/Y/SZcz8yDtdnpzfDqVpsHS7pY6iqzlTVu1X1O+B7vHfq5ZJalyQfYT7s36+qHw7DLY6Z9Rp3395gkORjST5x7jbwWeAF5tdj77DbXuCx6cxw6hZbhyPAncMrIK4H3hr5Uby9884Vf575Ywbm1+X2JB9Ncg3zvzz86VrPby0kCXAf8HJVfXtkU49jpqrW5QXYDfwX8HPga9OezxTX4Y+B/xguL55bC+APmP9N/6vAvwFXTXuua7AWDzF/iuH/mD8fetdi6wCE+Vdc/Rx4HpiZ9vzXeF3+cfi6jzEfrc0j+39tWJdXgJunPf8JrssNzJ9yOQY8N1x2dzlmfPsBSWpovZ6WkSR9COMuSQ0Zd0lqyLhLUkPGXZIaMu6S1JBxl6SG/h8MMAsiqI1wAQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "site_moa_counts = bbbc022_meta_moa[['Metadata_moa.x','Metadata_Plate', 'Metadata_Well', 'Metadata_Site']].drop_duplicates().groupby(['Metadata_moa.x']).size()\n",
    "plt.figure()\n",
    "plt.hist(site_moa_counts, bins = 50)\n",
    "well_moa_counts = bbbc022_meta_moa[['Metadata_moa.x','Metadata_Plate', 'Metadata_Well']].drop_duplicates().groupby(['Metadata_moa.x']).size()\n",
    "plt.figure()\n",
    "plt.hist(well_moa_counts, bins = 50)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Calculate treatments per MoA"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(array([292.,  32.,  16.,  13.,   3.,   4.,   4.,   6.,   0.,   3.,   1.,\n",
       "          0.,   1.,   1.,   0.,   1.,   0.,   1.,   0.,   1.,   0.,   0.,\n",
       "          1.,   0.,   2.,   1.,   1.,   1.,   0.,   0.,   0.,   2.,   1.,\n",
       "          0.,   0.,   0.,   0.,   1.,   0.,   0.,   0.,   0.,   0.,   0.,\n",
       "          0.,   0.,   0.,   0.,   0.,   1.]),\n",
       " array([ 1.  ,  2.08,  3.16,  4.24,  5.32,  6.4 ,  7.48,  8.56,  9.64,\n",
       "        10.72, 11.8 , 12.88, 13.96, 15.04, 16.12, 17.2 , 18.28, 19.36,\n",
       "        20.44, 21.52, 22.6 , 23.68, 24.76, 25.84, 26.92, 28.  , 29.08,\n",
       "        30.16, 31.24, 32.32, 33.4 , 34.48, 35.56, 36.64, 37.72, 38.8 ,\n",
       "        39.88, 40.96, 42.04, 43.12, 44.2 , 45.28, 46.36, 47.44, 48.52,\n",
       "        49.6 , 50.68, 51.76, 52.84, 53.92, 55.  ]),\n",
       " <a list of 50 Patch objects>)"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAPa0lEQVR4nO3df8idZ33H8fdnSf2Bytrax5Al2dJpNokwU3noKvpHbVFjlaWCKy2bBumIf7RQwTGq/6hjhQrTbsJWiLYYh1qDWhu0bGax4PzD1ic1tk1iMWpKE9Lk0bZakXUkfvfHcwWP8UmeH+f5kXP5fsHh3Pd1X/c53yu58zk317nPnVQVkqS+/MFyFyBJWniGuyR1yHCXpA4Z7pLUIcNdkjq0crkLALjkkktq/fr1y12GJI2UvXv3/rSqxqbbdl6E+/r165mYmFjuMiRppCR54mzbZpyWSfKiJA8l+X6S/Uk+2tovTfJgkkNJvpjkBa39hW39UNu+fqEGIkmandnMuT8PXFVVrwU2AZuTXAF8DLijql4FPAPc2PrfCDzT2u9o/SRJS2jGcK8pv2yrF7RHAVcBX2rtO4Br2/KWtk7bfnWSLFjFkqQZzepqmSQrkuwDTgC7gR8Bz1bVydblCLCmLa8BngRo238OvHya19yWZCLJxOTk5HCjkCT9llmFe1WdqqpNwFrgcuDVw75xVW2vqvGqGh8bm/bLXknSPM3pOveqehZ4AHg9cGGS01fbrAWOtuWjwDqAtv0PgZ8tSLWSpFmZzdUyY0kubMsvBt4MHGQq5N/Vum0F7mvLu9o6bfs3y1tPStKSms117quBHUlWMPVhsLOqvpbkAHBPkn8Cvgfc1frfBfxHkkPA08D1i1C3JOkcZgz3qnoEuGya9h8zNf9+Zvv/An+9INVJkublvPiF6jDW3/r1s247fPvbl7ASSTp/eOMwSeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SerQjOGeZF2SB5IcSLI/yS2t/SNJjibZ1x7XDOzzwSSHkjye5K2LOQBJ0u9aOYs+J4EPVNXDSV4G7E2yu227o6r+ebBzko3A9cBrgD8C/jvJn1XVqYUsXJJ0djOeuVfVsap6uC0/BxwE1pxjly3APVX1fFX9BDgEXL4QxUqSZmdOc+5J1gOXAQ+2ppuTPJLk7iQXtbY1wJMDux3h3B8GkqQFNutwT/JS4MvA+6vqF8CdwCuBTcAx4ONzeeMk25JMJJmYnJycy66SpBnMKtyTXMBUsH+uqr4CUFXHq+pUVf0a+BS/mXo5Cqwb2H1ta/stVbW9qsaranxsbGyYMUiSzjCbq2UC3AUcrKpPDLSvHuj2TuCxtrwLuD7JC5NcCmwAHlq4kiVJM5nN1TJvAN4NPJpkX2v7EHBDkk1AAYeB9wFU1f4kO4EDTF1pc5NXykjS0pox3Kvq20Cm2XT/Ofa5DbhtiLokSUPwF6qS1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUodmDPck65I8kORAkv1JbmntFyfZneSH7fmi1p4kn0xyKMkjSV632IOQJP222Zy5nwQ+UFUbgSuAm5JsBG4F9lTVBmBPWwd4G7ChPbYBdy541ZKkc5ox3KvqWFU93JafAw4Ca4AtwI7WbQdwbVveAny2pnwHuDDJ6gWvXJJ0VnOac0+yHrgMeBBYVVXH2qangFVteQ3w5MBuR1rbma+1LclEkonJyck5li1JOpdZh3uSlwJfBt5fVb8Y3FZVBdRc3riqtlfVeFWNj42NzWVXSdIMZhXuSS5gKtg/V1Vfac3HT0+3tOcTrf0osG5g97WtTZK0RGZztUyAu4CDVfWJgU27gK1teStw30D7e9pVM1cAPx+YvpEkLYGVs+jzBuDdwKNJ9rW2DwG3AzuT3Ag8AVzXtt0PXAMcAn4FvHdBK5YkzWjGcK+qbwM5y+arp+lfwE1D1iVJGoK/UJWkDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktShGcM9yd1JTiR5bKDtI0mOJtnXHtcMbPtgkkNJHk/y1sUqXJJ0drM5c/8MsHma9juqalN73A+QZCNwPfCats+/J1mxUMVKkmZnxnCvqm8BT8/y9bYA91TV81X1E+AQcPkQ9UmS5mGYOfebkzzSpm0uam1rgCcH+hxpbb8jybYkE0kmJicnhyhDknSm+Yb7ncArgU3AMeDjc32BqtpeVeNVNT42NjbPMiRJ05lXuFfV8ao6VVW/Bj7Fb6ZejgLrBrqubW2SpCU0r3BPsnpg9Z3A6StpdgHXJ3lhkkuBDcBDw5UoSZqrlTN1SPIF4ErgkiRHgA8DVybZBBRwGHgfQFXtT7ITOACcBG6qqlOLU7ok6WxmDPequmGa5rvO0f824LZhipIkDcdfqEpShwx3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHZox3JPcneREkscG2i5OsjvJD9vzRa09ST6Z5FCSR5K8bjGLlyRNbzZn7p8BNp/Rdiuwp6o2AHvaOsDbgA3tsQ24c2HKlCTNxYzhXlXfAp4+o3kLsKMt7wCuHWj/bE35DnBhktULVawkaXbmO+e+qqqOteWngFVteQ3w5EC/I63tdyTZlmQiycTk5OQ8y5AkTWfoL1SrqoCax37bq2q8qsbHxsaGLUOSNGC+4X789HRLez7R2o8C6wb6rW1tkqQlNN9w3wVsbctbgfsG2t/Trpq5Avj5wPSNJGmJrJypQ5IvAFcClyQ5AnwYuB3YmeRG4Angutb9fuAa4BDwK+C9i1CzJGkGM4Z7Vd1wlk1XT9O3gJuGLUqSNBx/oSpJHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOmS4S1KHDHdJ6pDhLkkdMtwlqUOGuyR1yHCXpA4Z7pLUIcNdkjpkuEtShwx3SeqQ4S5JHTLcJalDK4fZOclh4DngFHCyqsaTXAx8EVgPHAauq6pnhitTkjQXC3Hm/qaq2lRV4239VmBPVW0A9rR1SdISWoxpmS3Ajra8A7h2Ed5DknQOw4Z7Ad9IsjfJtta2qqqOteWngFXT7ZhkW5KJJBOTk5NDliFJGjTUnDvwxqo6muQVwO4kPxjcWFWVpKbbsaq2A9sBxsfHp+0jSZqfoc7cq+poez4B3AtcDhxPshqgPZ8YtkhJ0tzMO9yTvCTJy04vA28BHgN2AVtbt63AfcMWKUmam2GmZVYB9yY5/Tqfr6r/TPJdYGeSG4EngOuGL1OSNBfzDveq+jHw2mnafwZcPUxRkqTh+AtVSeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6NOy9Zc5r62/9+rTth29/+xJXIklLyzN3SeqQ4S5JHTLcJalDhrskdchwl6QOGe6S1CHDXZI6ZLhLUocMd0nqkOEuSR3q+vYDZ+NtCST1zjN3SeqQ4S5JHfq9nJY5G6drJPXCM3dJ6pDhLkkdMtwlqUOGuyR1yC9Ul5Bf2EpaKosW7kk2A/8KrAA+XVW3L9Z7LZezhbUkLbdFCfckK4B/A94MHAG+m2RXVR1YjPdbbMsV4p7pL61R+fMelTq1vBbrzP1y4FBV/RggyT3AFmAkw/18M9d/3PP5cJrray1X/7NZyABcqA/3xQ7f5Qx9P3BmttR/RqmqhX/R5F3A5qr6u7b+buAvq+rmgT7bgG1t9c+Bx2fx0pcAP13gcs8nvY8P+h+j4xt9ozTGP6mqsek2LNsXqlW1Hdg+l32STFTV+CKVtOx6Hx/0P0bHN/p6GeNiXQp5FFg3sL62tUmSlsBihft3gQ1JLk3yAuB6YNcivZck6QyLMi1TVSeT3Az8F1OXQt5dVfsX4KXnNI0zgnofH/Q/Rsc3+roY46J8oSpJWl7efkCSOmS4S1KHRiLck2xO8niSQ0luXe56FkKSu5OcSPLYQNvFSXYn+WF7vmg5axxGknVJHkhyIMn+JLe09i7GmORFSR5K8v02vo+29kuTPNiO1S+2CwpGVpIVSb6X5GttvbfxHU7yaJJ9SSZaWxfH6Hkf7gO3MngbsBG4IcnG5a1qQXwG2HxG263AnqraAOxp66PqJPCBqtoIXAHc1P7eehnj88BVVfVaYBOwOckVwMeAO6rqVcAzwI3LWONCuAU4OLDe2/gA3lRVmwaube/iGD3vw52BWxlU1f8Bp29lMNKq6lvA02c0bwF2tOUdwLVLWtQCqqpjVfVwW36OqYBYQydjrCm/bKsXtEcBVwFfau0jOz6AJGuBtwOfbuuho/GdQxfH6CiE+xrgyYH1I62tR6uq6lhbfgpYtZzFLJQk64HLgAfpaIxtymIfcALYDfwIeLaqTrYuo36s/gvwD8Cv2/rL6Wt8MPWB/I0ke9stUaCTY9T7uZ+nqqqSjPx1qkleCnwZeH9V/WLq5G/KqI+xqk4Bm5JcCNwLvHqZS1owSd4BnKiqvUmuXO56FtEbq+poklcAu5P8YHDjKB+jo3Dm/vt0K4PjSVYDtOcTy1zPUJJcwFSwf66qvtKauxojQFU9CzwAvB64MMnpk6ZRPlbfAPxVksNMTYVexdT/z9DL+ACoqqPt+QRTH9CX08kxOgrh/vt0K4NdwNa2vBW4bxlrGUqbn70LOFhVnxjY1MUYk4y1M3aSvJip/7vgIFMh/67WbWTHV1UfrKq1VbWeqX9z36yqv6GT8QEkeUmSl51eBt4CPEYvx+go/EI1yTVMzf+dvpXBbctc0tCSfAG4kqnbix4HPgx8FdgJ/DHwBHBdVZ35petISPJG4H+AR/nNnO2HmJp3H/kxJvkLpr5sW8HUSdLOqvrHJH/K1JnuxcD3gL+tqueXr9LhtWmZv6+qd/Q0vjaWe9vqSuDzVXVbkpfTwzE6CuEuSZqbUZiWkSTNkeEuSR0y3CWpQ4a7JHXIcJekDhnuktQhw12SOvT/+r5jZQ6jGogAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "compound_bbbc022_counts = bbbc022_meta_moa[['Metadata_moa.x','Treatment']].drop_duplicates().groupby(['Metadata_moa.x']).size()\n",
    "#compound_bbbc022_counts[compound_bbbc022_counts>30]\n",
    "plt.hist(compound_bbbc022_counts, bins = 50)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Calculate single cells per well/treatment/MoA"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(array([100., 162.,  43.,  29.,  13.,  10.,   3.,   5.,   1.,   3.,   4.,\n",
       "          1.,   1.,   0.,   1.,   0.,   0.,   0.,   0.,   2.,   0.,   1.,\n",
       "          0.,   1.,   0.,   0.,   0.,   0.,   3.,   1.,   0.,   0.,   0.,\n",
       "          0.,   1.,   0.,   1.,   1.,   0.,   1.,   0.,   1.,   0.,   0.,\n",
       "          0.,   0.,   0.,   0.,   0.,   1.]),\n",
       " array([4.0000000e+00, 3.1173800e+03, 6.2307600e+03, 9.3441400e+03,\n",
       "        1.2457520e+04, 1.5570900e+04, 1.8684280e+04, 2.1797660e+04,\n",
       "        2.4911040e+04, 2.8024420e+04, 3.1137800e+04, 3.4251180e+04,\n",
       "        3.7364560e+04, 4.0477940e+04, 4.3591320e+04, 4.6704700e+04,\n",
       "        4.9818080e+04, 5.2931460e+04, 5.6044840e+04, 5.9158220e+04,\n",
       "        6.2271600e+04, 6.5384980e+04, 6.8498360e+04, 7.1611740e+04,\n",
       "        7.4725120e+04, 7.7838500e+04, 8.0951880e+04, 8.4065260e+04,\n",
       "        8.7178640e+04, 9.0292020e+04, 9.3405400e+04, 9.6518780e+04,\n",
       "        9.9632160e+04, 1.0274554e+05, 1.0585892e+05, 1.0897230e+05,\n",
       "        1.1208568e+05, 1.1519906e+05, 1.1831244e+05, 1.2142582e+05,\n",
       "        1.2453920e+05, 1.2765258e+05, 1.3076596e+05, 1.3387934e+05,\n",
       "        1.3699272e+05, 1.4010610e+05, 1.4321948e+05, 1.4633286e+05,\n",
       "        1.4944624e+05, 1.5255962e+05, 1.5567300e+05]),\n",
       " <a list of 50 Patch objects>)"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAD4CAYAAAAAczaOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAQ6klEQVR4nO3df4ylVX3H8fenIOCvsvxYt9vdtYOVtDFNFLJBiKaxUJVfcWmCBGLKQmk2aW2i1UQXTdqY9I+lbURMGnQjtqtBhaKUDdpaCpimf4guiIgiZcRFdgPsagFriY3Ub/+4Z/GyzrIzs3fuvbPn/Uom93nOc+7c75zMfObcc5/73FQVkqQ+/MqkC5AkjY+hL0kdMfQlqSOGviR1xNCXpI4cOekCAE488cSamZmZdBmStKzcfffdP6yqlQu5z1SE/szMDDt27Jh0GZK0rCR5ZKH3cXlHkjpi6EtSRwx9SeqIoS9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I6MhXvyJU0fjObvzhn+84t5425Eo2TM31J6oihL0kdMfQlqSOGviR1xNCXpI549o6k5/GsnsObM31J6si8Qj/JziTfSnJvkh2t7fgktyV5qN0e19qT5KNJZpPcl+TUpfwBJEnzt5CZ/u9V1euqan3b3wzcXlUnA7e3fYBzgJPb1ybg2lEVK0k6NIeyvLMB2Na2twEXDLV/qga+CqxIsvoQHkeSNCLzDf0C/jXJ3Uk2tbZVVfVY234cWNW21wCPDt13V2t7niSbkuxIsmPv3r2LKF2StFDzPXvnjVW1O8krgNuSfHf4YFVVklrIA1fVVmArwPr16xd0X0nS4sxrpl9Vu9vtHuBm4DTgiX3LNu12T+u+G1g3dPe1rU2SNGEHDf0kL03y8n3bwFuA+4HtwMbWbSNwS9veDlzazuI5HXh6aBlIkjRB81neWQXcnGRf/89U1b8k+TpwY5IrgEeAi1r/LwHnArPAM8DlI69akrQoBw39qnoYeO0c7T8CzpqjvYB3jqQ6SdJI+Y5cSeqIoS9JHTH0Jakjhr4kdcTQl6SOGPqS1BE/REU6TBzow0/G8f39gJXlw5m+JHXE0Jekjhj6ktQRQ1+SOmLoS1JHDH1J6oihL0kdMfQlqSOGviR1xNCXpI4Y+pLUEUNfkjpi6EtSRwx9SeqIoS9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I6YuhLUkcMfUnqiKEvSR2Zd+gnOSLJN5Lc2vZPSnJXktkkNyQ5qrUf3fZn2/GZpSldkrRQC5npvwt4YGj/KuDqqno18CRwRWu/AniytV/d+kmSpsC8Qj/JWuA84BNtP8CZwE2tyzbggra9oe3Tjp/V+kuSJmy+M/2PAO8Dft72TwCeqqpn2/4uYE3bXgM8CtCOP936S5Im7KChn+R8YE9V3T3KB06yKcmOJDv27t07ym8tSTqAI+fR5w3A25KcCxwD/CpwDbAiyZFtNr8W2N367wbWAbuSHAkcC/xo/29aVVuBrQDr16+vQ/1BpMPNzOYvztm+c8t5Y65Eh5ODzvSr6sqqWltVM8DFwB1V9Q7gTuDC1m0jcEvb3t72acfvqCpDXZKmwKGcp/9+4D1JZhms2V/X2q8DTmjt7wE2H1qJkqRRmc/yznOq6ivAV9r2w8Bpc/T5KfD2EdQmaQ4HWvaR5sN35EpSRwx9SeqIoS9JHVnQmr4kzcXTS5cPZ/qS1BFDX5I6YuhLUkcMfUnqiKEvSR0x9CWpI4a+JHXE0Jekjhj6ktQRQ1+SOmLoS1JHDH1J6oihL0kdMfQlqSOGviR1xNCXpI4Y+pLUET85S5qwA33qlLQUnOlLUkcMfUnqiKEvSR0x9CWpI4a+JHXE0Jekjhj6ktQRQ1+SOnLQ0E9yTJKvJflmkm8n+VBrPynJXUlmk9yQ5KjWfnTbn23HZ5b2R5Akzdd8Zvr/C5xZVa8FXgecneR04Crg6qp6NfAkcEXrfwXwZGu/uvWTJE2Bg4Z+Dfyk7b6ofRVwJnBTa98GXNC2N7R92vGzkmRkFUuSFm1ea/pJjkhyL7AHuA34HvBUVT3buuwC1rTtNcCjAO3408AJoyxakrQ487rgWlX9H/C6JCuAm4HfPtQHTrIJ2ATwyle+8lC/naQpdKCLye3cct6YK9E+Czp7p6qeAu4EzgBWJNn3T2MtsLtt7wbWAbTjxwI/muN7ba2q9VW1fuXKlYssX5K0EPM5e2dlm+GT5MXAm4EHGIT/ha3bRuCWtr297dOO31FVNcqiJUmLM5/lndXAtiRHMPgncWNV3ZrkO8DnkvwV8A3gutb/OuDTSWaB/wIuXoK6pWXH6+ZrGhw09KvqPuCUOdofBk6bo/2nwNtHUp0kaaR8R64kdcSPS5RGzGUcTTNn+pLUEUNfkjpi6EtSRwx9SeqIoS9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I6YuhLUkcMfUnqiKEvSR0x9CWpI4a+JHXE0Jekjhj6ktQRQ1+SOmLoS1JHDH1J6oihL0kdMfQlqSOGviR15MhJFyAtVzObvzjpEpatA43dzi3njbmS/jjTl6SOONOXXoCzeR1unOlLUkcMfUnqyEFDP8m6JHcm+U6Sbyd5V2s/PsltSR5qt8e19iT5aJLZJPclOXWpfwhJ0vzMZ03/WeC9VXVPkpcDdye5DbgMuL2qtiTZDGwG3g+cA5zcvl4PXNtupanl2r16cdCZflU9VlX3tO3/Bh4A1gAbgG2t2zbggra9AfhUDXwVWJFk9cgrlyQt2ILW9JPMAKcAdwGrquqxduhxYFXbXgM8OnS3Xa1NkjRh8w79JC8DPg+8u6p+PHysqgqohTxwkk1JdiTZsXfv3oXcVZK0SPMK/SQvYhD411fVF1rzE/uWbdrtnta+G1g3dPe1re15qmprVa2vqvUrV65cbP2SpAWYz9k7Aa4DHqiqDw8d2g5sbNsbgVuG2i9tZ/GcDjw9tAwkSZqg+Zy98wbgD4FvJbm3tX0A2ALcmOQK4BHgonbsS8C5wCzwDHD5SCuWJC3aQUO/qv4DyAEOnzVH/wLeeYh1SZKWgO/IlaSOeME1HZa8dK80N2f6ktQRQ1+SOuLyjrriNXbUO2f6ktQRQ1+SOuLyjpYFz8aRRsOZviR1xNCXpI4Y+pLUEUNfkjpi6EtSRzx7R8uab7aSFsaZviR1xNCXpI4Y+pLUEUNfkjpi6EtSRwx9SeqIoS9JHTH0Jakjhr4kdcTQl6SOeBkGTRUvqyAtLWf6ktQRQ1+SOmLoS1JHDH1J6oihL0kdMfQlqSMHPWUzySeB84E9VfU7re144AZgBtgJXFRVTyYJcA1wLvAMcFlV3bM0pWs589RMaTLmM9P/B+Ds/do2A7dX1cnA7W0f4Bzg5Pa1Cbh2NGVKkkbhoDP9qvr3JDP7NW8A3tS2twFfAd7f2j9VVQV8NcmKJKur6rFRFSzp8HWgZ4A7t5w35koOX4td0181FOSPA6va9hrg0aF+u1rbL0myKcmOJDv27t27yDIkSQtxyC/ktll9LeJ+W6tqfVWtX7ly5aGWIUmah8WG/hNJVgO02z2tfTewbqjf2tYmSZoCiw397cDGtr0RuGWo/dIMnA487Xq+JE2P+Zyy+VkGL9qemGQX8JfAFuDGJFcAjwAXte5fYnC65iyDUzYvX4KaJUmLNJ+zdy45wKGz5uhbwDsPtShJ0tLwevoaCd9sJS0PXoZBkjpi6EtSR1zekTT1fKfu6DjTl6SOGPqS1BFDX5I6YuhLUkcMfUnqiGfvaEF8E5a0vDnTl6SOGPqS1JHDdnnHN3NI0i9zpi9JHTH0Jakjhr4kdeSwXdNfjIWejujrA5KWG0N/jJbTi8uejy8dngz9jhnsWu6W00RqWhj6HTDcJe3jC7mS1JHuZvrjmPUu9DFGVZNPaSUdTHehP0rTtmwybfVImj7LPvQNOkmaP9f0Jakjhr4kdWTZL+9I0v48f//AnOlLUkcMfUnqiMs7krrhss8ShX6Ss4FrgCOAT1TVlqV4HEkahZ6usDvy5Z0kRwB/B5wDvAa4JMlrRv04kqSFW4qZ/mnAbFU9DJDkc8AG4DtL8FiSNHajfFPouJ81LEXorwEeHdrfBbx+/05JNgGb2u5Pkjy4iMc6EfjhIu43SdY8HtY8HtZ8iHLVvLodqObfWOjjTeyF3KraCmw9lO+RZEdVrR9RSWNhzeNhzeNhzeMxypqX4pTN3cC6of21rU2SNGFLEfpfB05OclKSo4CLge1L8DiSpAUa+fJOVT2b5M+ALzM4ZfOTVfXtUT9Oc0jLQxNizeNhzeNhzeMxsppTVaP6XpKkKedlGCSpI4a+JHVk2YZ+krOTPJhkNsnmSdezT5J1Se5M8p0k307yrtZ+fJLbkjzUbo9r7Uny0fZz3Jfk1AnVfUSSbyS5te2flOSuVtcN7UV5khzd9mfb8ZkJ1bsiyU1JvpvkgSRnLIMx/vP2O3F/ks8mOWYaxznJJ5PsSXL/UNuCxzbJxtb/oSQbx1zv37TfjfuS3JxkxdCxK1u9DyZ561D72DJlrpqHjr03SSU5se2Pdoyratl9MXiB+HvAq4CjgG8Cr5l0Xa221cCpbfvlwH8yuBzFXwObW/tm4Kq2fS7wz0CA04G7JlT3e4DPALe2/RuBi9v2x4A/adt/CnysbV8M3DChercBf9y2jwJWTPMYM3jT4veBFw+N72XTOM7A7wKnAvcPtS1obIHjgYfb7XFt+7gx1vsW4Mi2fdVQva9peXE0cFLLkSPGnSlz1dza1zE4CeYR4MSlGOOx/uKPcMDOAL48tH8lcOWk6zpArbcAbwYeBFa3ttXAg23748AlQ/2f6zfGGtcCtwNnAre2X64fDv3RPDfe7RfyjLZ9ZOuXMdd7bAvQ7Nc+zWO8753qx7dxuxV467SOMzCzX4guaGyBS4CPD7U/r99S17vfsT8Arm/bz8uKfeM8iUyZq2bgJuC1wE5+EfojHePlurwz16Ue1kyolgNqT8lPAe4CVlXVY+3Q48Cqtj0NP8tHgPcBP2/7JwBPVdWzc9T0XL3t+NOt/zidBOwF/r4tSX0iyUuZ4jGuqt3A3wI/AB5jMG53M93jPGyhYzvxMR/yRwxmyjDF9SbZAOyuqm/ud2ikNS/X0J96SV4GfB54d1X9ePhYDf4tT8W5sknOB/ZU1d2TrmUBjmTw1PjaqjoF+B8GSw7PmaYxBmhr4BsY/MP6deClwNkTLWqRpm1sX0iSDwLPAtdPupYXkuQlwAeAv1jqx1quoT/Vl3pI8iIGgX99VX2hNT+RZHU7vhrY09on/bO8AXhbkp3A5xgs8VwDrEiy7817wzU9V287fizwozHWC4MZza6quqvt38Tgn8C0jjHA7wPfr6q9VfUz4AsMxn6ax3nYQsd24mOe5DLgfOAd7R8VL1DXpOv9TQYTgm+2v8W1wD1Jfu0FaltUzcs19Kf2Ug9JAlwHPFBVHx46tB3Y9+r6RgZr/fvaL22v0J8OPD30NHrJVdWVVbW2qmYYjOMdVfUO4E7gwgPUu+/nuLD1H+usr6oeBx5N8lut6SwGl+6eyjFufgCcnuQl7XdkX81TO877WejYfhl4S5Lj2rOct7S2scjgg5zeB7ytqp4ZOrQduLidHXUScDLwNSacKVX1rap6RVXNtL/FXQxOCHmcUY/xUr5QscQvgpzL4MyY7wEfnHQ9Q3W9kcFT3/uAe9vXuQzWY28HHgL+DTi+9Q+DD535HvAtYP0Ea38Tvzh751UM/hhmgX8Ejm7tx7T92Xb8VROq9XXAjjbO/8Tg7IWpHmPgQ8B3gfuBTzM4g2Tqxhn4LIPXHX7WwueKxYwtg7X02fZ1+ZjrnWWw3r3vb/BjQ/0/2Op9EDhnqH1smTJXzfsd38kvXsgd6Rh7GQZJ6shyXd6RJC2CoS9JHTH0Jakjhr4kdcTQl6SOGPqS1BFDX5I68v9YvWfQ9W9A0gAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAOsUlEQVR4nO3db4xU13nH8e9T8J8kbgOYFaJgulhBjVDV1mjlErmKItMmjokCLyyLKEpWLhVSm7ZOUynFjVSrUl/gqkriSFUiZNwSyY3tEksgu21EMFHUFyEB4/gP1GFNcAwCs6mNk7RSHZqnL+aQTle7wMyd3dk58/1Iq7n33Dsz5x60P545e++dyEwkSXX5hX53QJLUe4a7JFXIcJekChnuklQhw12SKrSw3x0AWLp0aY6Ojva7G5I0UI4cOfLDzByZbtu8CPfR0VEOHz7c725I0kCJiFdm2ua0jCRVyHCXpAoZ7pJUIcNdkipkuEtShQx3SaqQ4S5JFTLcJalChrskVWheXKEqSQCj25+atv3Ujo1z3JPBZ+UuSRWycpc071nRd87KXZIqZLhLUoWuGO4R8XBEnI+IF9ralkTE/og4UR4Xl/aIiC9ExEREPBcR62az85Kk6V1N5f4PwB1T2rYDBzJzDXCgrAN8EFhTfrYBX+xNNyVJnbhiuGfmN4HXpzRvAnaX5d3A5rb2L2fLt4BFEbG8V52VJF2dbs+WWZaZZ8vyOWBZWV4BvNq23+nSdpYpImIbreqeVatWddkNScPMs2hm1vgPqpmZQHbxvJ2ZOZaZYyMj034FoCSpS92G+2uXplvK4/nSfga4qW2/laVNkjSHug33fcB4WR4H9ra1f7ycNbMeeLNt+kaSNEeuOOceEV8B3gcsjYjTwP3ADuDxiNgKvALcXXb/Z+BOYAL4L+CeWeizJOkKrhjumfmRGTZtmGbfBD7RtFOSpGa8QlWSKmS4S1KFDHdJqpDhLkkVMtwlqUKGuyRVyHCXpAoZ7pJUIcNdkipkuEtShQx3SaqQ4S5JFTLcJalChrskVchwl6QKGe6SVCHDXZIqZLhLUoUMd0mqkOEuSRW64hdkS9KgGd3+1LTtp3ZsnOOe9I+VuyRVyMpd0pybqbJW71i5S1KFDHdJqpDhLkkVMtwlqUKGuyRVyHCXpAoZ7pJUoUbnuUfEnwK/DyTwPHAPsBx4FLgROAJ8LDPfathPSfOYV4TOP11X7hGxAvgTYCwzfw1YAGwBHgA+l5nvAt4Atvaio5Kkq9d0WmYh8LaIWAi8HTgL3A7sKdt3A5sbvockqUNdh3tmngH+FvgBrVB/k9Y0zIXMvFh2Ow2saNpJSVJnmkzLLAY2AauBXwbeAdzRwfO3RcThiDg8OTnZbTckSdNoMi3zO8D3M3MyM38KPAHcBiwq0zQAK4Ez0z05M3dm5lhmjo2MjDTohiRpqibh/gNgfUS8PSIC2AAcAw4Cd5V9xoG9zbooSepU16dCZuahiNgDPANcBI4CO4GngEcj4q9L265edFRSf3mb3sHS6Dz3zLwfuH9K80ng1iavK0lqxitUJalChrskVchwl6QKGe6SVCHDXZIqZLhLUoUMd0mqkOEuSRUy3CWpQoa7JFXIcJekChnuklQhw12SKmS4S1KFDHdJqpDhLkkVMtwlqUKGuyRVyHCXpAoZ7pJUoUZfkC1JlzO6/al+d2FoWblLUoUMd0mqkOEuSRVyzl3S/+M8eR2s3CWpQoa7JFXIcJekChnuklQhw12SKuTZMtKQ8qyYujWq3CNiUUTsiYh/j4jjEfGeiFgSEfsj4kR5XNyrzkqSrk7TaZkHgX/NzHcDvwEcB7YDBzJzDXCgrEuS5lDX4R4R7wTeC+wCyMy3MvMCsAnYXXbbDWxu2klJUmeaVO6rgUng7yPiaEQ8FBHvAJZl5tmyzzlg2XRPjohtEXE4Ig5PTk426IYkaaom4b4QWAd8MTNvAf6TKVMwmZlATvfkzNyZmWOZOTYyMtKgG5KkqZqE+2ngdGYeKut7aIX9axGxHKA8nm/WRUlSp7oO98w8B7waEb9amjYAx4B9wHhpGwf2NuqhJKljTc9z/2PgkYi4FjgJ3EPrP4zHI2Ir8Apwd8P3kCR1qFG4Z+azwNg0mzY0eV1JUjPefkCSKmS4S1KFvLeMNE/NdO+XUzs2znFPNIis3CWpQlbu0oCxotfVsHKXpAoZ7pJUIcNdkipkuEtShQx3SaqQ4S5JFTLcJalChrskVchwl6QKGe6SVCHDXZIq5L1lpMrNdC8a1c3KXZIqZLhLUoUMd0mqkOEuSRUy3CWpQp4tI/WZZ7NoNli5S1KFrNylSvgJQO2s3CWpQlbuUo/NVEGf2rFxjnuiqYbp38bKXZIqZLhLUoUMd0mqkOEuSRVqHO4RsSAijkbEk2V9dUQcioiJiHgsIq5t3k1JUid6UbnfCxxvW38A+Fxmvgt4A9jag/eQJHWgUbhHxEpgI/BQWQ/gdmBP2WU3sLnJe0iSOte0cv888GngZ2X9RuBCZl4s66eBFdM9MSK2RcThiDg8OTnZsBuSpHZdh3tEfAg4n5lHunl+Zu7MzLHMHBsZGem2G5KkaTS5QvU24MMRcSdwPfBLwIPAoohYWKr3lcCZ5t2UJHWi68o9M+/LzJWZOQpsAZ7OzI8CB4G7ym7jwN7GvZQkdWQ2znP/c+BTETFBaw5+1yy8hyTpMnpy47DM/AbwjbJ8Eri1F68rSeqOV6hKUoUMd0mqkOEuSRUy3CWpQn4Tk3QZl/te0hq/vWdY1fgNTVbuklQhK3epS5er6qV+s3KXpAoZ7pJUIcNdkirknLuE8+eqj5W7JFXIyl1DxQpdw8LKXZIqZLhLUoUMd0mqkOEuSRUy3CWpQp4toyp5VoyGnZW7JFXIcJekChnuklQhw12SKmS4S1KFDHdJqpDhLkkVMtwlqUKGuyRVyCtUpTniVbOaS1buklQhw12SKtR1uEfETRFxMCKORcSLEXFvaV8SEfsj4kR5XNy77kqSrkaTyv0i8GeZuRZYD3wiItYC24EDmbkGOFDWJUlzqOtwz8yzmflMWf4xcBxYAWwCdpfddgObm3ZSktSZnpwtExGjwC3AIWBZZp4tm84By2Z4zjZgG8CqVat60Q0NIc9AkabX+A+qEXED8FXgk5n5o/ZtmZlATve8zNyZmWOZOTYyMtK0G5KkNo3CPSKuoRXsj2TmE6X5tYhYXrYvB84366IkqVNNzpYJYBdwPDM/27ZpHzBelseBvd13T5LUjSZz7rcBHwOej4hnS9tfADuAxyNiK/AKcHezLkqSOtV1uGfmvwExw+YN3b6uJKk5r1CVpAoZ7pJUIe8KKUkzmOk6ilM7Ns5xTzpn5S5JFTLcJalChrskVcg593lgkOf15or3kJE6Y+UuSRWqtnKfi2rYilvSfGXlLkkVMtwlqUKGuyRVyHCXpAoZ7pJUoWrPltH84BlFUn9YuUtShazcr4JXR0oaNFbuklShoavcL1eFOw/cf35KknrDyl2SKjR0lfvlzHbV2K/Xn4+fSKzQpdll5S5JFRr4yt0KsHvd/P1hkD4dSLNlEH4PrNwlqUIDX7nXrFfVQTefbvxEJA02K3dJqpCV+wCaj1X1fOyTNNfm03U0Vu6SVCHDXZIqZLhLUoUMd0mq0KyEe0TcEREvRcRERGyfjfeQJM2s52fLRMQC4O+A3wVOA9+JiH2ZeazX7zVfeeaIpKnm+qrW2ajcbwUmMvNkZr4FPApsmoX3kSTNYDbOc18BvNq2fhr4rak7RcQ2YFtZ/UlEvNTl+y0Fftjlc2viODgGlzgOLQMxDvFAo6f/ykwb+nYRU2buBHY2fZ2IOJyZYz3o0kBzHByDSxyHlmEfh9mYljkD3NS2vrK0SZLmyGyE+3eANRGxOiKuBbYA+2bhfSRJM+j5tExmXoyIPwK+BiwAHs7MF3v9Pm0aT+1UwnFwDC5xHFqGehwiM/vdB0lSj3mFqiRVyHCXpAoNdLjXfJuDiHg4Is5HxAttbUsiYn9EnCiPi0t7RMQXyjg8FxHr2p4zXvY/ERHj/TiWJiLipog4GBHHIuLFiLi3tA/NWETE9RHx7Yj4bhmDvyrtqyPiUDnWx8oJDETEdWV9omwfbXut+0r7SxHxgf4cUTMRsSAijkbEk2V9KMfhijJzIH9o/bH2ZeBm4Frgu8Dafverh8f3XmAd8EJb298A28vyduCBsnwn8C9AAOuBQ6V9CXCyPC4uy4v7fWwdjsNyYF1Z/kXge8DaYRqLciw3lOVrgEPl2B4HtpT2LwF/UJb/EPhSWd4CPFaW15bfk+uA1eX3Z0G/j6+L8fgU8I/Ak2V9KMfhSj+DXLlXfZuDzPwm8PqU5k3A7rK8G9jc1v7lbPkWsCgilgMfAPZn5uuZ+QawH7hj9nvfO5l5NjOfKcs/Bo7Tugp6aMaiHMtPyuo15SeB24E9pX3qGFwamz3AhoiI0v5oZv53Zn4fmKD1ezQwImIlsBF4qKwHQzgOV2OQw3262xys6FNf5sqyzDxbls8By8ryTGNR1RiVj9W30Kpch2osylTEs8B5Wv8xvQxcyMyLZZf24/n5sZbtbwI3MuBjUHwe+DTws7J+I8M5Dlc0yOE+1LL1+XJozmONiBuArwKfzMwftW8bhrHIzP/JzN+kdcX3rcC7+9ylORcRHwLOZ+aRfvdlEAxyuA/jbQ5eK1MMlMfzpX2msahijCLiGlrB/khmPlGah3IsMvMCcBB4D60pp0sXIrYfz8+PtWx/J/AfDP4Y3AZ8OCJO0ZqGvR14kOEbh6syyOE+jLc52AdcOstjHNjb1v7xcqbIeuDNMmXxNeD9EbG4nE3y/tI2MMoc6S7geGZ+tm3T0IxFRIxExKKy/DZa35VwnFbI31V2mzoGl8bmLuDp8ulmH7ClnEWyGlgDfHtujqK5zLwvM1dm5iit3/enM/OjDNk4XLV+/0W3yQ+tMyO+R2v+8TP97k+Pj+0rwFngp7TmBLfSmi88AJwAvg4sKfsGrS9IeRl4Hhhre53fo/UHowngnn4fVxfj8Nu0plyeA54tP3cO01gAvw4cLWPwAvCXpf1mWqE0AfwTcF1pv76sT5TtN7e91mfK2LwEfLDfx9ZgTN7H/50tM7TjcLkfbz8gSRUa5GkZSdIMDHdJqpDhLkkVMtwlqUKGuyRVyHCXpAoZ7pJUof8FqS3i9JNFHfoAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYMAAAD4CAYAAAAO9oqkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAASqElEQVR4nO3de4yld13H8ffHji0WlW3ZoS67i7NoiylGpI61BDWFIiyFsDUhZBsvC9RslIIgRGwhsWpCUsDIJSiwQmHRWlhrtQ14KxVFE9k6LbT0Sode6Gy27CCCRhNk5esfz2/p2WF2Z+Y858zMwvuVnMzz/J7L77u/Ped85rmcOakqJEnf2b5rrQuQJK09w0CSZBhIkgwDSRKGgSQJmFjrAgA2btxYU1NTa12GJJ1Qbrnlli9V1eQo9rUuwmBqaoqZmZm1LkOSTihJHhrVvjxNJEkyDCRJhoEkCcNAkoRhIEnCMJAkYRhIkjAMJEksIwySXJXkUJI7FrS/Ksk9Se5M8paB9suTzCa5N8nzxlG0JGm0lvMJ5A8C7wI+dKQhybOAHcDTquprSZ7Q2s8GdgJPBZ4IfDzJWVX1f6MufDmmLvvYou0PXvmCVa5Ekta3JY8MquqTwJcXNP8acGVVfa2tc6i17wA+XFVfq6oHgFng3BHWK0kag2GvGZwF/EyS/Un+KclPtvbNwMMD6821tm+RZHeSmSQz8/PzQ5YhSRqFYcNgAjgdOA/4TWBfkqxkB1W1p6qmq2p6cnIkf3RPkjSkYcNgDriuOjcD3wA2AgeArQPrbWltkqR1bNgw+CvgWQBJzgJOBr4E3ADsTHJKkm3AmcDNoyhUkjQ+S95NlOQa4HxgY5I54ArgKuCqdrvp/wK7qqqAO5PsA+4CDgOXrtWdRJKk5VsyDKrq4mMs+sVjrP8m4E19ipIkrS4/gSxJMgwkSYaBJAnDQJKEYSBJwjCQJGEYSJIwDCRJGAaSJAwDSRKGgSQJw0CShGEgScIwkCRhGEiSMAwkSRgGkiSWEQZJrkpyqH3F5cJlr0tSSTa2+SR5Z5LZJLcnOWccRUuSRms5RwYfBLYvbEyyFXgu8IWB5ucDZ7bHbuDd/UuUJI3bkmFQVZ8EvrzIorcBrwdqoG0H8KHqfArYkGTTSCqVJI3NUNcMkuwADlTVbQsWbQYeHpifa22L7WN3kpkkM/Pz88OUIUkakRWHQZJTgTcAv92n46raU1XTVTU9OTnZZ1eSpJ4mhtjmh4BtwG1JALYAtyY5FzgAbB1Yd0trkyStYys+Mqiqz1bVE6pqqqqm6E4FnVNVjwA3AL/c7io6D/hqVR0cbcmSpFFbzq2l1wD/CjwlyVySS46z+l8D9wOzwB8DrxhJlZKksVryNFFVXbzE8qmB6QIu7V+WJGk1+QlkSZJhIEkyDCRJGAaSJAwDSRKGgSQJw0CShGEgScIwkCRhGEiSMAwkSRgGkiQMA0kShoEkCcNAkoRhIElied90dlWSQ0nuGGh7a5J7ktye5C+TbBhYdnmS2ST3JnneuAqXJI3Oco4MPghsX9B2I/CjVfVjwOeAywGSnA3sBJ7atvmjJCeNrFpJ0lgsGQZV9Ungywva/r6qDrfZTwFb2vQO4MNV9bWqeoDuu5DPHWG9kqQxGMU1g5cDf9OmNwMPDyyba23fIsnuJDNJZubn50dQhiRpWL3CIMkbgcPA1Svdtqr2VNV0VU1PTk72KUOS1NPEsBsmeSnwQuCCqqrWfADYOrDaltYmSVrHhjoySLIdeD3woqr6n4FFNwA7k5ySZBtwJnBz/zIlSeO05JFBkmuA84GNSeaAK+juHjoFuDEJwKeq6ler6s4k+4C76E4fXVpV/zeu4iVJo7FkGFTVxYs0v/84678JeFOfoiRJq8tPIEuSDANJkmEgScIwkCRhGEiSMAwkSRgGkiQMA0kShoEkCcNAkoRhIEnCMJAkYRhIkjAMJEkYBpIkDANJEoaBJIllhEGSq5IcSnLHQNvpSW5Mcl/7eVprT5J3JplNcnuSc8ZZvCRpNJZzZPBBYPuCtsuAm6rqTOCmNg/wfODM9tgNvHs0ZUqSxmnJMKiqTwJfXtC8A9jbpvcCFw20f6g6nwI2JNk0qmIlSeMxMeR2Z1TVwTb9CHBGm94MPDyw3lxrO8gCSXbTHT3wpCc9acgyYOqyjw29rSSp0/sCclUVUENst6eqpqtqenJysm8ZkqQehg2DLx45/dN+HmrtB4CtA+ttaW2SpHVs2DC4AdjVpncB1w+0/3K7q+g84KsDp5MkSevUktcMklwDnA9sTDIHXAFcCexLcgnwEPCStvpfAxcCs8D/AC8bQ82SpBFbMgyq6uJjLLpgkXULuLRvUZKk1eUnkCVJhoEkyTCQJGEYSJIwDCRJGAaSJAwDSRKGgSQJw0CShGEgScIwkCRhGEiSMAwkSRgGkiQMA0kShoEkiZ5hkOQ3ktyZ5I4k1yR5TJJtSfYnmU3ykSQnj6pYSdJ4DB0GSTYDvw5MV9WPAicBO4E3A2+rqh8G/gO4ZBSFSpLGp+9pognge5JMAKcCB4FnA9e25XuBi3r2IUkas6HDoKoOAL8PfIEuBL4K3AJ8paoOt9XmgM19i5QkjVef00SnATuAbcATgccC21ew/e4kM0lm5ufnhy1DkjQCfU4TPQd4oKrmq+rrwHXAM4EN7bQRwBbgwGIbV9WeqpququnJyckeZUiS+uoTBl8AzktyapIAFwB3AZ8AXtzW2QVc369ESdK49blmsJ/uQvGtwGfbvvYAvwW8Nsks8Hjg/SOoU5I0RhNLr3JsVXUFcMWC5vuBc/vsV5K0uvwEsiTJMJAkGQaSJAwDSRKGgSQJw0CShGEgScIwkCRhGEiSMAwkSRgGkiQMA0kShoEkCcNAkoRhIEnCMJAkYRhIkugZBkk2JLk2yT1J7k7yjCSnJ7kxyX3t52mjKlaSNB59jwzeAfxtVf0I8DTgbuAy4KaqOhO4qc1LktaxocMgyeOAn6V94X1V/W9VfQXYAextq+0FLupbpCRpvPocGWwD5oEPJPl0kvcleSxwRlUdbOs8Apyx2MZJdieZSTIzPz/fowxJUl99wmACOAd4d1U9HfhvFpwSqqoCarGNq2pPVU1X1fTk5GSPMiRJffUJgzlgrqr2t/lr6cLhi0k2AbSfh/qVKEkat6HDoKoeAR5O8pTWdAFwF3ADsKu17QKu71WhJGnsJnpu/yrg6iQnA/cDL6MLmH1JLgEeAl7Ssw9J0pj1CoOq+gwwvciiC/rsV5K0uvwEsiTJMJAkGQaSJAwDSRKGgSQJw0CShGEgScIwkCRhGEiSMAwkSRgGkiQMA0kShoEkCcNAkoRhIEnCMJAk0f+bzk5IU5d9bNH2B698wSpXIknrQ+8jgyQnJfl0ko+2+W1J9ieZTfKR9pWYkqR1bBSniV4N3D0w/2bgbVX1w8B/AJeMoA9J0hj1CoMkW4AXAO9r8wGeDVzbVtkLXNSnD0nS+PU9Mng78HrgG23+8cBXqupwm58DNi+2YZLdSWaSzMzPz/csQ5LUx9BhkOSFwKGqumWY7atqT1VNV9X05OTksGVIkkagz91EzwRelORC4DHA9wPvADYkmWhHB1uAA/3LlCSN09BHBlV1eVVtqaopYCfwD1X1C8AngBe31XYB1/euUpI0VuP40NlvAa9NMkt3DeH9Y+hDkjRCI/nQWVX9I/CPbfp+4NxR7FeStDr8cxSSJMNAkmQYSJIwDCRJGAaSJAwDSRKGgSQJw0CSxHfoN50di9+AJuk7lUcGkiTDQJJkGEiSMAwkSRgGkiQMA0kShoEkiR5hkGRrkk8kuSvJnUle3dpPT3Jjkvvaz9NGV64kaRz6HBkcBl5XVWcD5wGXJjkbuAy4qarOBG5q85KkdWzoMKiqg1V1a5v+L+BuYDOwA9jbVtsLXNS3SEnSeI3kmkGSKeDpwH7gjKo62BY9ApxxjG12J5lJMjM/Pz+KMiRJQ+odBkm+F/gL4DVV9Z+Dy6qqgFpsu6raU1XTVTU9OTnZtwxJUg+9wiDJd9MFwdVVdV1r/mKSTW35JuBQvxIlSePW526iAO8H7q6qPxhYdAOwq03vAq4fvjxJ0mro8yesnwn8EvDZJJ9pbW8ArgT2JbkEeAh4Sb8SJUnjNnQYVNW/ADnG4guG3a8kafX5CWRJkmEgSTIMJEn4HcjL4ncjS/p255GBJMkwkCQZBpIkvGbQi9cSJH278MhAkmQYSJIMA0kShoEkCcNAkoR3E60q7z6StF55ZCBJ8sjgROQRhqRRMwzG4Fhv1mu1H60PhrjWs7GFQZLtwDuAk4D3VdWV4+pL3z58w5TWxliuGSQ5CfhD4PnA2cDFSc4eR1+SpP7GdWRwLjBbVfcDJPkwsAO4a0z9ndDW8rTSsX7jPlF+Q19pnccbo7X6t63VWA/T70qfY6N6fp0oz8fjWe//hlTV6HeavBjYXlW/0uZ/CfipqnrlwDq7gd1t9inAvUN2txH4Uo9yx20912dtw7G24VjbcI5X2w9W1eQoOlmzC8hVtQfY03c/SWaqanoEJY3Feq7P2oZjbcOxtuGsVm3j+pzBAWDrwPyW1iZJWofGFQb/BpyZZFuSk4GdwA1j6kuS1NNYThNV1eEkrwT+ju7W0quq6s5x9MUITjWN2Xquz9qGY23DsbbhrEptY7mALEk6sfi3iSRJhoEkCaiqE/YBbKf7fMIscNkY+9kKfILuQ3N3Aq9u7acDNwL3tZ+ntfYA72x13Q6cM7CvXW39+4BdA+0/AXy2bfNO2im8ZdZ3EvBp4KNtfhuwv+3rI8DJrf2UNj/blk8N7OPy1n4v8LxRjTGwAbgWuAe4G3jGOhq332j/n3cA1wCPWauxA64CDgF3DLSNfZyO1ccyantr+z+9HfhLYMOw4zHMmB+vtoFlrwMK2Lhexq21v6qN3Z3AW9Zi3BZ9Dq70xb1eHnRvgJ8HngycDNwGnD2mvjYdeeIA3wd8ju7PbLzlyH8OcBnw5jZ9IfA37cl3HrB/4Al0f/t5Wps+8gK/ua2btu3zV1Dfa4E/49Ew2AfsbNPvAX6tTb8CeE+b3gl8pE2f3cbvlPYE+3wb395jDOwFfqVNn0wXDms+bsBm4AHgewbG7KVrNXbAzwLncPQb7tjH6Vh9LKO25wITbfrNA7WteDxWOuZL1dbat9LdwPIQj4bBehi3ZwEfB05p809Yi3Fb9Dk4yjfN1XzQ/Yb5dwPzlwOXr1Lf1wM/R5fWm1rbJuDeNv1e4OKB9e9tyy8G3jvQ/t7Wtgm4Z6D9qPWWqGULcBPwbOCj7Un7JR59oX5znNqL4xlteqKtl4Vjd2S9vmMMPI7uDTcL2tfDuG0GHqZ7A5hoY/e8tRw7YIqj3zjGPk7H6mOp2hYs+3ng6sX+nUuNxzDP1+XURnc0+jTgQR4NgzUfN7o38Ocsst6qj9vCx4l8zeDIi/mIudY2VkmmgKfTHX6dUVUH26JHgDOWqO147XOLtC/H24HXA99o848HvlJVhxfZ1zf7b8u/2tZfab3LtQ2YBz6Q5NNJ3pfksayDcauqA8DvA18ADtKNxS2sn7GD1RmnY/WxEi+n+615mNqGeb4eV5IdwIGqum3BovUwbmcBP5Nkf5J/SvKTQ9Y28nE7kcNg1SX5XuAvgNdU1X8OLqsuhmuV63khcKiqblnNfldggu4w+d1V9XTgv+kOqb9pLcYNIMlpdH88cRvwROCxdOdm16XVGKdh+kjyRuAwcPVYilqhJKcCbwB+e7X6XOG4TdAdjZ4H/CawL0nGVdtKnMhhsKp/8iLJd9MFwdVVdV1r/mKSTW35JrqLRcer7XjtWxZpX8ozgRcleRD4MN2poncAG5Ic+UDh4L6+2X9b/jjg34eod7nmgLmq2t/mr6ULh7UeN4DnAA9U1XxVfR24jm4818vYweqM07H6WFKSlwIvBH6hvSEOU9u/s/IxP54fogv429rrYgtwa5IfGKK2cYzbHHBddW6mO6LfOERtox63E/qawQTdhZ5tPHph5alj6ivAh4C3L2h/K0dfRHpLm34BR1+ourm1n053Dv209ngAOL0tW3ih6sIV1ng+j15A/nOOvrD0ijZ9KUdfWNrXpp/K0Rev7qe7cNV7jIF/Bp7Spn+njdmajxvwU3R3c5zatt1Ld5fHmo0d33p+eezjdKw+llHbdrq76yYXrLfi8VjpmC9V24JlD/LoNYP1MG6/Cvxemz6L7nRO1mLcvqXWlbyw19uD7u6Az9FdbX/jGPv5abrDwNuBz7THhXTn4W6iu73s4wNPoNB9uc/n6W5Lmx7Y18vpbvmaBV420D5Nd4vj54F3sYJbJNv25/NoGDy5PYln2xPmyJ0Lj2nzs235kwe2f2Pr+14G7sjpO8bAjwMzbez+qr3Y1sW4Ab9Ld4vfHcCftBfimowd3a2tB4Gv0/32eMlqjNOx+lhGbbN0b2RHXg/vGXY8hhnz49W2YPmDHH1r6VqP28nAn7Z93go8ey3GbbGHf45CknRCXzOQJI2IYSBJMgwkSYaBJAnDQJKEYSBJwjCQJAH/DwK8dBGigTxqAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure()\n",
    "plt.hist(bbbc022_meta.groupby(['Metadata_Plate','Metadata_Well']).agg('sum')['Cells'], bins = 50)\n",
    "plt.figure()\n",
    "#Cut DMSO as it has 1144258 cells \n",
    "plt.hist(bbbc022_meta[~bbbc022_meta['Treatment'].isin(['DMSO@0'])].groupby(['Treatment']).agg('sum')['Cells'], bins = 50)\n",
    "plt.figure()\n",
    "plt.hist(bbbc022_meta_moa.groupby(['Metadata_moa.x']).agg('sum')['Cells'], bins = 50)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Calculate stats for the full CDRP dataset (after QC)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "cdrp_meta = pd.read_csv(\"/dgx1nas1/cellpainting-datasets/CDRP/inputs/metadata/index.csv\")\n",
    "#the number of cells is inferred from locations\n",
    "cells_cdrp = []\n",
    "for i in cdrp_meta.index:\n",
    "    cells_cdrp.append([\n",
    "        cdrp_meta.loc[i, \"Metadata_Plate\"],\n",
    "        cdrp_meta.loc[i, \"Metadata_Well\"],\n",
    "        cdrp_meta.loc[i, \"Metadata_Site\"],\n",
    "        sum(1 for line in open('/dgx1nas1/cellpainting-datasets/CDRP/inputs/locations/{}/{}-{}-Nuclei.csv'.format(cdrp_meta.loc[i, \"Metadata_Plate\"], cdrp_meta.loc[i, \"Metadata_Well\"], cdrp_meta.loc[i, \"Metadata_Site\"]))) - 1\n",
    "    ])\n",
    "counts = pd.DataFrame(data=cells_cdrp, columns=['Metadata_Plate', 'Metadata_Well', 'Metadata_Site', 'Cells'])\n",
    "cdrp_meta = pd.merge(cdrp_meta, counts, right_on=['Metadata_Plate', 'Metadata_Well', 'Metadata_Site'], left_on =['Metadata_Plate', 'Metadata_Well', 'Metadata_Site'])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Calculate wells/sites/cells per MoA"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "cdrp_moa_matches = pd.read_csv(\"./cdrp/data/CDRP_MOA_MATCHES_fixed.csv\")\n",
    "cdrp_moa_matches = cdrp_moa_matches.assign(var1=cdrp_moa_matches['Metadata_moa.x'].str.split('|')).explode('var1')\n",
    "cdrp_moa_matches.reset_index(inplace=True, drop = True)\n",
    "cdrp_moa_matches.drop(columns=['Unnamed: 0', 'Unnamed: 0.1', 'Metadata_moa.x'], inplace= True)\n",
    "cdrp_moa_matches.rename(columns={\"var1\": \"Metadata_moa.x\"}, inplace=True)\n",
    "\n",
    "cdrp_meta_moa = pd.merge(cdrp_meta, cdrp_moa_matches[['Var1', 'Metadata_moa.x']], right_on=['Var1'], left_on =['Compound'])\n",
    "cdrp_meta_moa"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Calculate wells/sites/cells per MoA"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "site_moa_counts_cdrp = cdrp_meta_moa[['Metadata_moa.x','Metadata_Plate', 'Metadata_Well', 'Metadata_Site']].drop_duplicates().groupby(['Metadata_moa.x']).size()\n",
    "plt.figure()\n",
    "plt.hist(site_moa_counts_cdrp, bins = 50)\n",
    "well_moa_counts_cdrp = cdrp_meta_moa[['Metadata_moa.x','Metadata_Plate', 'Metadata_Well']].drop_duplicates().groupby(['Metadata_moa.x']).size()\n",
    "plt.figure()\n",
    "plt.hist(well_moa_counts_cdrp, bins = 50)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Calculate treatments per MoA"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "compound_cdrp_counts = cdrp_meta_moa[['Metadata_moa.x','Treatment']].drop_duplicates().groupby(['Metadata_moa.x']).size()\n",
    "plt.hist(compound_cdrp_counts, bins = 50)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Calculate single cells per well/treatment/MoA"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "plt.figure()\n",
    "plt.hist(cdrp_meta.groupby(['Metadata_Plate','Metadata_Well']).agg('sum')['Cells'], bins = 50)\n",
    "plt.figure()\n",
    "#Cut NA (DMSO)\n",
    "plt.hist(cdrp_meta[~cdrp_meta['Treatment'].isin(['NA@NA'])].groupby(['Treatment']).agg('sum')['Cells'], bins = 50)\n",
    "plt.figure()\n",
    "plt.hist(cdrp_meta_moa.groupby(['Metadata_moa.x']).agg('sum')['Cells'], bins = 50)"
   ]
  },
  {
   "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.7.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}