[33bf40]: / myeloid / HSC_plotting_publish.ipynb

Download this file

265 lines (264 with data), 71.8 kB

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Loading required package: limma\n",
      "\n",
      "\n",
      "Attaching package: ‘dplyr’\n",
      "\n",
      "\n",
      "The following objects are masked from ‘package:stats’:\n",
      "\n",
      "    filter, lag\n",
      "\n",
      "\n",
      "The following objects are masked from ‘package:base’:\n",
      "\n",
      "    intersect, setdiff, setequal, union\n",
      "\n",
      "\n",
      "\n",
      "Attaching package: ‘reshape’\n",
      "\n",
      "\n",
      "The following object is masked from ‘package:dplyr’:\n",
      "\n",
      "    rename\n",
      "\n",
      "\n",
      "The following object is masked from ‘package:Matrix’:\n",
      "\n",
      "    expand\n",
      "\n",
      "\n",
      "Loading required package: Rcpp\n",
      "\n"
     ]
    }
   ],
   "source": [
    "library(Matrix)\n",
    "library(edgeR)\n",
    "library(dplyr)\n",
    "library(Seurat)\n",
    "library(biomaRt)\n",
    "library(reshape)\n",
    "library(ggplot2)\n",
    "library(fgsea)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "\n",
    "\n",
    "***HSC priming - plotting***\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "#Import summary scores - \n",
    "#DEGs are calculated between the three committed HSC lineages - myeloid, erythroid and MK and then those with adjusted \n",
    "#p value <0.05 used to form a signature to probe for enrichment in the uncommitted clusters (CD38neg and CD38 pos) - all in python\n",
    "#these scores are them imported for plotting in to R\n",
    "\n",
    "setwd('/home/jovyan/data/covid/HSC_data/')\n",
    "values<-read.csv('summary_scores.csv', row.names = 1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 56,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "`summarise()` ungrouping output (override with `.groups` argument)\n",
      "\n"
     ]
    }
   ],
   "source": [
    "test<-values %>% filter(temp == \"Uncommitted\") %>% \n",
    "group_by(Status_on_day_collection_summary) %>%\n",
    "summarise(erythroid_mean = mean(erythroid_score), MK_mean = mean(MK_score), myeloid_mean = mean(myeloid_score))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 57,
   "metadata": {},
   "outputs": [],
   "source": [
    "keep <-c(\"Healthy\", \"Asymptomatic\", \"Mild\", \"Moderate\", \"Severe\", \"Critical\")\n",
    "test<-filter(test, Status_on_day_collection_summary %in% keep)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 58,
   "metadata": {},
   "outputs": [],
   "source": [
    "test$Status_on_day_collection_summary <- factor(test$Status_on_day_collection_summary, \n",
    "                                                  levels = c(\"Healthy\", \"Asymptomatic\", \"Mild\", \"Moderate\", \n",
    "                                                             \"Severe\", \"Critical\"))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 59,
   "metadata": {},
   "outputs": [],
   "source": [
    "test$erythroid_standard<-test$erythroid_mean / test$erythroid_mean[[3]]\n",
    "test$myeloid_standard<-test$myeloid_mean / test$myeloid_mean[[3]]\n",
    "test$MK_standard<-test$MK_mean / test$MK_mean[[3]]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 60,
   "metadata": {},
   "outputs": [],
   "source": [
    "setwd('/home/jovyan/data/covid/HSC_data/figures')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 61,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0gAAANICAMAAADKOT/pAAACcFBMVEUAAAABAQECAgIDAwMF\nBQUICAgJCQkMDAwNDQ0PDw8QEBAREREUFBQVFRUXFxcZGRkiIiIpKSkzMzM6Ojo7OztAQEBB\nQUFDQ0NHR0dISEhJSUlLS0tNTU1OTk5PT09QUFBRUVFSUlJTU1NUVFRVVVVWVlZXV1dYWFhZ\nWVlaWlpbW1tcXFxdXV1eXl5fX19gYGBhYWFiYmJjY2NkZGRlZWVmZmZnZ2doaGhpaWlqampr\na2tsbGxtbW1ubm5vb29wcHBxcXFycnJzc3N0dHR1dXV2dnZ3d3d4eHh5eXl6enp7e3t8fHx9\nfX1+fn5/f3+AgICBgYGCgoKDg4OEhISFhYWGhoaHh4eIiIiJiYmKioqLi4uMjIyNjY2Ojo6P\nj4+QkJCRkZGSkpKTk5OUlJSVlZWWlpaXl5eYmJiZmZmampqbm5ucnJydnZ2enp6fn5+goKCh\noaGioqKjo6OkpKSlpaWmpqanp6eoqKipqamqqqqrq6usrKytra2urq6vr6+wsLCxsbGysrKz\ns7O0tLS1tbW2tra3t7e4uLi5ubm6urq7u7u8vLy9vb2+vr6/v7/AwMDBwcHCwsLDw8PExMTF\nxcXGxsbHx8fIyMjJycnKysrLy8vMzMzNzc3Ozs7Pz8/Q0NDR0dHS0tLTfX3T09PU1NTV1dXW\n1tbX19fY2NjZ2dna2trb29vc3Nzd3d3e3t7f39/g4ODh4eHi4uLj4+Pk5OTl5eXm5ubn5+fo\n6Ojp6enq6urr6+vs7Ozt7e3u7u7v7+/w8PDx8fHy8vLz8/P09PT19fX29vb39/f4+Pj5+fn6\n+vr7+/v8/Pz9/f3+/v7////QwEwcAAAACXBIWXMAABJ0AAASdAHeZh94AAAgAElEQVR4nO2d\n+38U1bqn98zsM2fO7Jk5ZzwzZ87MrE4ghEAIlxCCkVuI4RpRDDchCBFEhAgCWy5RdrgJCBsB\niRhFBJGBAyqiGMRswCC3EC5J1780a1V3VVd1V3ev7rzVXZ33+/ywrVr1VtUKn372WrWufzAA\nAAPmD/nOAACDAYgEAAEQCQACIBIABEAkAAiASAAQAJEAIAAiAUAARKJk6z+15TsLID9AJErW\n/mFzvrMA8gNEogQisQUiUQKR2AKRKIFIbIFIlEAktkAkSiASWyASJRCJLRCJEojEFohECURi\nC0SiBCKxBSJRApHYApEogUhsgUiUQCS2QCRKIBJbIBIlEIktEIkSiMQWiEQJRGILRKIEIrEF\nIlECkdgCkSiBSGyBSJRAJLZApCy4f/pwW9uhU3cSLkAktkCkGHcXCkm6qGdHpodEhIk777mv\nZSvS//MgqweBfAGRbE6PFhoifVstHJSfdF2ESGyBSFEer426kTqsvVi4OeS8CpHYApEiXKkR\nOiKdNz0a8+evOzvPbKlQx0XfOi5DJLZAJEX/+0qQssZ0Ij0x63VrHkfOHr2mzuY4rkMktkAk\nyYNZyoiZN/enE8kMWGOf9tWr887YdYjEFogkuSp1CP25z0gr0nR5vaI3dn5O3bA7dg6R2AKR\nDFOkcecNI61ID9T1jY6E/jKZsCp2DpHYApEMJdLC++q/6UT6pbFuXPEXzpQp8oaFsVOIxBaI\nJOn8MPLftFU7STjsPKuVNzTGTiESWyCSAx2R3JSjagdMIJKDjEW6GdclC5HYApEcZCzSezK+\nuDt2DpHYApEcZCpS9wgZv8KRAJHYApEcZChS+FUZPqzLkQKR2AKRHGQo0gYVvs+ZApHYApEc\nZCbSDhW9zD7tu3Tp0oI/bs7qxRCp4IFIDjIRKdyiguc+sRO6n5P84+asXgyRCh6I5CADkXoX\nmx49iqX0rF69uvI/bc7qxRCp4IFIDvRFujVNhS5+EpeMbyS2QCQH2iJdMuf0reqPT4dIbIFI\nDnRFah+q5l20Jl6ASGyBSA40RdqjVhEqafe4ApHYApEc6In0gYqq+NbrEkRiC0RyoCWS6VH1\nr57XIBJbIJIDHZGOqpj6u94XIRJbIJIDDZHM9bhm9ia5CpHYApEcpBfp/hgZMbUn2WWIxBaI\n5CC9SE2qneF20ssQiS0QyYGXSNePS25FT86ogC/ib4sBkdgCkRx4ibRXpVlL5asFIV9J8QCI\nxBaIZBi/zo0y0RyIGuXryEWnSF+p4+fnxnPZfhJEYgtEiiy0msjxyEWnSK2egeJL+0kQiS0Q\nCSIBAiASRAIEQCRKIBJbIBIlEIktEIkSiMQWiEQJRGILRKIEIrEFIlECkdgCkSiBSGyBSJRA\nJLZAJEogElsgEiUQiS0QiRKIxBaIRAlEYgtEogQisQUiUQKR2AKRKIFIbIFIlEAktkAkSiAS\nWyASJRCJLRCJEojEFohECURiC0SiBCKxBSJRApHYApEogUhsgUiUQCS2QCRKIBJbIBIlEIkt\nEImSQhMJApMBkSiBSGyBSJRAJLZAJEogElsgEiUQiS0QiRKIxBaIRAlEYgtEogQisQUiUQKR\n2AKRKIFIbIFIlEAktkAkSiASWyASJRCJLRCJEojEFohECURiC0SiBCKxBSJRApHYApEogUhs\ngUiUQCS2QCRKIBJbIBIlQRFJ93kQiQyIRAlEYgtEogQisQUiUQKR2AKRKIFIbIFIlEAktkAk\nSiASWyASJRCJLRCJEojEFohESbxI+fpBQ6ScA5EogUhsgUiUQCS2QCRKIBJbIBIlEIktEIkS\niMQWiETJoBUJxqUDIlECkdgCkSiBSGyBSJRAJLZAJEogElsgEiUQiS0QiRKIxBaIRAlEYgtE\nogQisQUiUQKR2AKRKIFIbIFIlEAktkAkSiASWyASJRCJLRCJEojEFohECURiC0SiBCKxBSJR\nApHYApEogUhsgUiUQCS2QCRKIBJbIBIlEIktEIkSiMQWiEQJRGILRKIEIrEFIlECkdgCkSiB\nSGyBSJRAJLZAJEogElsgEiUQiS0QiRKIxBaIRAlEYgtEogQisQUiUQKR2AKRKIFIbIFIlEAk\ntkAkSiASWyASJRCJLRCJEojEFohECURiC0SiBCKxBSJRApHYApEogUhsgUiUQCS2QCRKIBJb\nIBIlEIktEIkSiMQWiEQJRGILRLLpPXe0bdfRi8/SR948eXDn7iNf3U+4AJHYApEihE80DBEm\npYvPp4y8815lJFCEpu3vdV+DSGyBSCa/zRIOmp8kDQzvK3FGVpx2XYVIbIFIituVwkXD02SR\nqyJlUVXt5FGR0H3OqxCJLRBJ8qxWKTF+26WbXRc2lavjVUki31cXKw+ZX0c3NqnKYMhZJkEk\ntkAkSavSo+lR5KTbrOWd8wy8pep1U+5ap5fVaU1f7DpEYgtEMoyHZdKHObYP98bJ03rPyK2q\nCLoWO9+lnDsTO4dIbIFIhrFP2lDcFTv/TOnxnVdknbwwy3F+PyQTWmLnEIktEMkwXpQ2vOY4\nD4+XCeu9IkfIC+84EybIhIWxU4jEFohk3FHFyqfOlBaZUOUVqiJbnQnTVaUwdgqR2AKRjJOq\nJnfPmXJKpdzyCB0p07c4E6pQIgETiGRskjJMcKV0K5FOeoTOcXtjhFUD+NbYOURiC0QyFksZ\n5ruThsuk9z1C22R6ycPY+YW4ZgmIxBaIZEyTMqxxJ01K0id7r8x1oW92XCseRGILRDLGSBve\ncye9JJMWeMUeUUXQmmjX7Z0FqoD60XHZb5HyFQeR0gKRjKFShzZ3kqrtzfQMblUNd6NWffRZ\n+4dNpWqsuGvUKkRiC0QKq0LmgDtthUya6h1+ZopjcGto0S+uixCJLRDpiTLikDutWSa9kCS+\ne63tUdnOO7H0/l9++WXZf9jsig26IBCJDIj02EMkNVmixjP67hvFzvkWxc32CNbu5yT/uNkV\nHXRBIBIZEKlPCXHQnbZSDfL2Cv6hQl4Z8ubp7qeP/ta+VH0vjbHGsD6YPXv2sP+82RUedEEg\nEhkQyRz3s9edtFQm1XmE3laTlSb/ap3+oAY2jHUs3YBvJLZAJEPJscOdNF8mvewR2qjEuRs7\n7ypz90FBJLZAJKMmcay3mjH7emLkVVULPOJMeU8mDI0VSRCJLRDJmCdlWOJOGpXYR6vYrj6Q\nHjtTfhSukeMQiS0QyVib0LJwX+lxLDFS1exqXSl9Idd4cIjEFohkHFLN2K4FuM4pkToTI9UM\nwHnuJDW64U37DCKxBSIZnUqbC84UtTTDqHBiZL1Mn+FKCasSaa19CpHYApEis/M2OBOmyoQ3\nPALVINUxrpTbSsJt9ilEYgtEiszsq3DU7S7HrQ1ks1ld+NWZctI9BRAisQUiGUZXyDXPtV/N\nMproUbOLfDttdKao5ocQmr8BRFKowd6xFVM3Kl2OW2fXj0ui6zf0P6/aJRwzYs2Fu5pi5xCJ\nLRBJ0q0WNSlq7VHHN19TcjTYBdJeZ+XNrMmVHo1e7NunBrAOvR57EERiC0RSfFmkDClpWNlU\np6p5orLbvuQSKVJaieq3D508eXBNZOX9447nQCS2QCST9qHOyRETbsSuuEUKbw0JF6WufluI\nxBaIFOGn2AZJxWt6HBfcIhnG5Vec05GWY4YsMIFIFj9sn11VUlaz8KM7ruR4kQyjc8+SF0YP\nKRk7ZcVHt+OeAZHYApEogUhsgUiUQCS2QCRKIBJbIBIlEIktEIkSiMQWiEQJRGILRKIEIrEF\nIlECkdgCkSiBSGyBSJRAJLZAJEogElsgEiUQiS0QiRKIxBaIRAlEYgtEogQisQUiUQKR2AKR\nKIFIbIFIlEAktkAkSiASWyASJRCJLRCJEojEFohECURiC0SiBCKxBSJRApHYApEogUhsgUiU\nQCS2QCRKIBJbIBIlEIktEIkSiMQWiEQJRGILRKIEIrEFIlECkdgCkSiBSGyBSJRAJLZAJEog\nElsgEiUQiS0QiRKIxBaIRAlEYgtEogQisQUiUQKR2AKRKIFIbIFIlEAktkAkSiASWyASJRCJ\nLRCJEojEFohECURiC0SiBCKxBSJRApHYApEogUhsgUiUQCS2QCRKIBJbIBIlEIktEIkSiMQW\niEQJRGILRKIEIrEFIlECkdgCkSiBSGyBSJRAJLZAJEogElsgEiUQiS0QiRKIxBaIRAlEYgtE\nogQisQUiUQKR2AKRKIFIbIFIlEAktkAkSiASWyASJRCJLRCJEojEFohECURiC0SiBCKxBSJR\nApHYApEogUhsgUiUQCS2QCRKIBJbIBIlEIktEIkSiMQWiEQJRGILRKIEIrEFIlECkdgCkSiB\nSGyBSJRAJLZAJEogElsgEiUQiS0QiRKIxBaIRAlEYgtEogQisQUiUQKR2AKRKIFIbIFIlEAk\ntkAkSiASWyASJRCJLRCJEojEFohECURiC0SiBCKxBSJRApHYApEogUhsgUiUQCS2QCRKIBJb\nIBIlEIktEIkSiMQWiEQJRGILRKIEIrEFIlECkdgCkSiBSGyBSJRAJLZAJEogElsgEiUQiS0Q\niRKIxBaIZNN77mjbrqMXnw3gERCJLRApQvhEwxBhUrr4vN4tu1X0e64kiMQWiGTy2yzhoPmJ\nxi2/lEAkYAORFLcrhYuGp2lvCc8UEAnYQCTJs1rlxPhtl252XdhUro5Xpb1nt4BIIAZEkrQq\nJZoeRU66zVreuTS3RCp2EAlEgUiG8bBMGjGnzzq9N06e1qe+RVXsimdCJGABkQxjnxSiuCt2\n/pkqar5LecsuGbFuBUQCFhDJMF6UQrzmOA+PlwnrU92hKnZldyESsIFIxp2QFOJTZ0qLTKhK\ncUf/DBlwyIBIwAYiGSdVTe6eM+WUSrmV/I42eXmeAZFADIhkbJI+THCldCuRTia94bqs2JXe\nhEjAAUQyFksf5ruThsuk95PF99fLqwcNiAQcQCRjmvRhjTtpUqo+2b+ooQ9hAyIBBxDJGJPg\ng/GSTFqQJLxzqBDDzMZyiARsIJIhxRBt7iRV25vpHW1W7PaZhxAJ2ECksGpZOOBOU4pM9Q7f\nqUZBhO0oiARMINITYXYKuWiWSS94Rv8sy6+SG5FjiARsINJjD5FWyaQar+D+Onlld/TELVL3\nc5J/3OwKD7ogEIkMiNQnoq3ZDlbKpClewR+oj6dw9MQt0r2Kiop//S+bXeFBFwQikQGRDDVC\naK87aalMqvMI/XmIrNhdt85QtQM2EMlQM/l2uJPmy6SXEyP76lwNfBAJ2EAkoyZxrLeaMft6\nYuT7aqJSv30KkYANRDLmSR+WuJNGJfbRSq7Jit3wX2PnEAnYQCRjbULLwn3V/nAsIXCOGhvU\nGkONLZqrDjrtR0EkrkAk45CaIOtagOucEqkzIbBaJKHdioBIbIFIRqdy4YIzZatMGBVOCIRI\nECkpEMkwqqQLG5wJU2XCG4lxEAkiJQUiRWb2VTjqdpeVHGfS34fGBmADkQyjS3XJbrVP+2fL\n04mJNbsEIBKwgUhGxIjQaetsoyqQjltn149LvNdvgEjABiJJukdKJYpae9TxzdeURw12gbRX\nnXqv3wCRgA1EUnxZpHQpaVjZVKeqeaKy274EkTIKZAtEMmkf6myFm3AjdgUiZRTIFogU4afY\nBknFa3ocFyBSRoFsgUgWP2yfXVVSVrPwozuuZIiUUSBbIBIlEIktEIkSiMQWiEQJRGILRKIE\nIrEFIlECkdgCkSiBSGyBSJRAJLZAJEogElsgEiUQiS0QiRKIxBaIRAlEYgtEogQisQUiUQKR\n2AKRKIFIbIFIlEAktkAkSiASWyASJRCJLRCJEojEFohECURiC0SiBCKxBSJRApHYApEogUhs\ngUiUQCS2QCRKIBJbIBIlEIktEIkSiMQWiEQJRGILRKIEIrEFIlECkdgCkSiBSGyBSJRAJLZA\nJEogElsgEiUQiS0QiRKIxBaIRAlEYgtEogQisQUiUQKR2AKRKIFIbIFIlEAktkAkSiASWyAS\nJRCJLRCJEojEFohECURiC0SiBCKxBSJRApHYApEogUhsgUiUQCS2QCRKIBJbIBIlEIktEIkS\niMQWiEQJRGILRKIEIrEFIlECkdgCkSiBSGyBSJRAJLZAJEogElsgEiUQiS0QiRKIxBaIRAlE\nYgtEogQisQUiUQKR2AKRKIFIbIFIlEAktkAkSiASWyASJRCJLRCJEojEFohECURiC0SiBCKx\nBSJRApHYApEogUhsgUiUQCS2QCRKIBJbIBIlEIktEIkSiMQWiEQJRGILRKIEIrEFIlECkdgC\nkSiBSGyBSJRAJLZAJEogElsgEiUQiS0QiRKIxBaIRAlEYgtEogQisQUiUQKR2AKRKIFIbIFI\nlEAktkAkSiASWyASJRCJLRCJEojEFohECURiC0SiBCKxBSJRApHYApEogUhsgUiUQCS2QCRK\nIBJbIBIlEIktEIkSiMQWiEQJRGILRKIEIrEFIlECkdgCkSiBSGyBSJRAJLZAJEogElsgEiUQ\niS0QiRKIxBaIRAlEYgtEogQisQUiUQKR2AKRbHrPHW3bdfTiswFEQiS2QKQI4RMNQ4RJ6eLz\nWUdCJLZAJJPfZgkHzU+yjIRIbIFIituVwkXD0+wiIRJbIJLkWa1yYvy2Sze7LmwqV8ersouE\nSGyBSJJWZUTTo8hJt1l3O5dVJERiC0QyjIdl0oc5fdbpvXHytD6rSIjEFohkGPukDsVdsfPP\nVEHzXTaREIktEMkwXpQ2vOY4D4+XCeuziYRIbIFIxp2QtOFTZ0qLTKjKJhIisQUiGSdV/eye\nM+WUSrmVRSREYgtEMjZJFya4UrqVHieziIRIbIFIxmLpwnx30nCZ9H4WkRCJLRDJmCZdWONO\nmuTdJ5s2EiKxBSIZY6QL77mTXpJJC7KIhEhsgUjGUOlCmztJ1eFmZhEJkdgCkcKqveCAO22F\nTJqaRSREYgtEeqL0OOROa5ZJL2QRCZHYApEee+ixSibVZBj5+7/8y7/883/b7LoadEEgEhkQ\nqU/pcdCdtlImTckw8m4oFPqf/3WzT7kEAQciGWrcz1530lKZVJdFZHzVDrABIhlqft4Od9J8\nmfRyFpEQiS0QyahJHOut5sG+nkUkRGILRDLmSReWuJNGJfa8akVCJLZAJGNtQsvCfdWqcCyL\nSIjEFohkHFLTXl3Lap1TenRmEQmR2AKRjE4lwwVnylaZMCqcRSREYgtEMowqacMGZ8JUmfBG\nNpEQiS0QKTJfr8JRY7usCp4z2URCJLZAJMPoUh2tW+3T/tnydKJHzS59JERiC0QyIkO4Q6et\ns42qmDlunV0/LrmlE2lAJMZAJEn3SClEUWuPOr75mrKjwS5m9rpWZUgVaUAkxkAkxZdFyomS\nhpVNdaryJiq77UtukVJFGhCJMRDJpF1NfrWZcCN2JU6kFJEGRGIMRIrwU2zbo+I1PY4L8SIl\njzQgEmMgksUP22dXlZTVLPzojis5QaSkkQZEYgxEogQisQUiUQKR2AKRKIFIbIFIlEAktkAk\nStb+YeYRnngOTeQERKLk8H//gyd//NPfeV+I49/96e+14v7wD/+gF/f3f/r3WnF/96c/asX9\n8U//0fvC/8r3P32+gUikHJrtSdVzpd4X4qh77l+14mb/j3/Wi/s/z03Xihv+XJVW3IRkf8g7\n+f6XzzcQKRecfK5FK677OY9FwLz4v/9bL27mcx77pXmw8bkTWnGfP7dO78XsgEi5ACINeiBS\nLoBIgx6IlAseXLqtFdd36breA698rxd3/dIzrbjfLt3Xint4SU9MfkAkAAiASAAQAJEAIAAi\n+cn8k3qfKKDggUh+IkR5y4/a0dc+iR2f33GZPjvAPyCSn5jzaKfvf6AT+/UMURk7OyzEpC8H\nnoHeK0d3blmX3kndOJAMiOQn0SnpQ5vOei6T5yC8XsYNi50fVvet6Y+e7U3Nnt07d2xKfOaj\ng/PMpVrEx+bp2Q33vN+dPG6WBrr/FoMciOQnY+3lHcZv7UoZucaMii3iaookVkXPhAbxTwwf\nGWNdigiyVoz0KuNSxWXzXqZAJD8JX1hbYf3gQg3HniQNbI/ExNb2+rzUTIiuPpnFD7q/OXYp\nIsgEmYeTRjwp4yCSNhDJZ/rPrhpp/ebK1nznHdRXKa+O+MBZaPV1TJRpVX3mSRY/6NejySFL\nkBvqrCShXEwZB5G0gUj+8+zLFcOtn92kXb97RJyUVybGD755OFVY6xeVOhhZGn1Ukb3C3tia\n2rlLE58oQo0fXXtiCdJ/TG2BG7fhYJq4Zg0G/M8zOIBIOeFJx9Jh0V998aJTffGXV8n0HxJu\nuiJiX0kxjpSL0Mzd57plPbG/59qnLbIwq0yosT1VX2f15sA9SxDD6JRFY+h2VnEgHRApV/Se\nWGQVIaPfjdsOcJIQsz1uqZclWHza2/KHf9WZED4uZdgYF6UKmrl2tTAqiHFCHu/OKg6kAyLl\nkIeH54SiLs086mx5kJWpdz3iW4QYFZe0W4j5T+PSblULsdedJIu4IdEyxSFIWH52Lc0qDqQD\nIuWWWzunRFUq3xabuVAkxE6P4FZZEXSn3C0RpQnruxrnhRjm7iOqFWJZ9NAhiPGOEDVZxSWh\nZ9q0+LKQKxApxzw48oLVRlC67XE0sVSILR6xskQqdafsEsLr635afFVMFnEHoodOQQ7HF3G6\ncUkIh8RknTgGQKRc0nOssdjZdFwZHZQj62ZNHuHzhKh2p7wixCGPwA2ywudKkG/5InroFKRD\n1uSyikvCL0KM1IljAETKGX2fL7EbrOu2TY804R01Ly0SYkT8l49h3JQfVAvdSeOEOOXxaFlQ\nVboS5IusfQWdghyPL+J047z/oH+rS6h7sgUi5YhrG+wxDmVrVbPbxXlmZ5BZIOwXHnW7vjki\noRFhqLdIu4UocSWMjZVcTkHahKjKKk7x8NCyuqrRMUaaDSflSf9iXkCkXPDgQJ1dnZv5sfVl\n1F4mTytUk8O9IUKEdvS77rmldnoeEtd9K0X60OP5smo31JUgX7cieugU5CUh6rOKk3xoD9Bw\nsSD5X80KiOQ74TNNdpWuvOVnx5Wf1BiCNnW0QV2cuOeaNUj8wanmISppa9yzxgvxiscraoQY\n70pYL4uoqIMOQa4VCbEhqzjDWOepkRjybZq/ngsQyWd+fa/S/tXN/STuQ+gLmThDHTyaFIko\nrq5veHlObXn0hhfih7mq2t5XCe84IlPd0xk6VFkRKeFigvTUJtytG6e+mryou6D77zDYgUh+\n0ntkjv2bq9h0IzGgXqabB10TPH6mNd3x8btkaml7XOLHqsBrdSWZo2AXmxMKbUE6lR/jn2UX\nZ061mLSyZaP8z9sbW5arEbU1HVibywYi+Yk1vlSEXj6ZMMBOsV1eihzdXZTg0dK7CfG3S9SF\n6fs6re+p8PUDteYL4hbEU80Xonx7pyVI75lms+H9iJFV3GcyadxFdSQPzL7f628Vick/GyAK\nRPKTqBFj37uZJOCgiM1DOPdqyGHR0IXnvW74QFh1wBfnvjy3vnpI9Hx9XFx4YSS9RHX/1tRV\nRR/d2J9d3FtS1R+sPyk6iOLbSlHqmUeWQCQ/Ub/JosZT8b/eGPuEGB07u9++aVH9pOops5Zs\n6+jxvqHPmj/kZk5vfGDPLI+w+oQwzbg6IV62/yRrNFLXOFHqUV3lCUTyEyGqdqScj/DT8eOn\nU11PJLy9JOF3H1rpMfe2b30oPmxjYqevZpz8RNpu/Um2SMbXQjRklvvBC0TykyVn0i16kgW3\nYtPXI7WyJUkWAv/59SJH2LBViTOetONK7I5bGRHr21ogBJq/I0CkAqT/4p4359ZOHDthyszl\nraceJw+8f2JDY+3z1bUNS98/m6SuqBdXbDc/yOIr1q5xUgj2O4xFgUhAg5H28PJSIS7Zybei\n3WAAIgEdJggRXThvjLWOhKIvpDfdggEQyXc623f++Z2WeNSVeZJO6yAJ2b3zzp07Xv1W+zdv\nvpjV8+YLMTNyVO+szfVh9LcFRPKXvn1eQxZEpPdI/fdb6yAJ2b1WeLcCLIufF7hw8uTJOouT\nb5EP/MY8WiFEhd1G+J3WdAsWQCRfedyQSpCci7RRiCmuBNWUnjh1PZHLMq7siBLokDxaYyUv\n8ppuwROI5CvLUwqiLVLGa3AnEWm91MGVoETSWuLfHFRbfN8wbqtep2XmSI27apHW+IXyuAKR\n/OT7qA3lldXxqMtNkl+tgySYz0khmnfR5S1S3xQhilwpNTJQa8Dcl+oVRapXrEkdhaYtWDbT\nHJ4UP3iPKxDJT9RYadF0Pvma35pkINIVE3n+8ZU4LrfPF/HfNCqDXusXJbJZRo5RB7+WOd86\nwXMsLkMgkp+ohRk+IHhOBiKlC6t1PfhupRAj9cbL7R8uppkHXzrWbym9muYuNkAkPxkjRHXy\nEav6ZLAGdzqR4sqfrjpZ0JzRykTP4V2Rg+8nWg+r60x9CyMgkp/Ir4h1OnGE3T6T44egupkR\nX80Mt88Nifo9V1MMNEqg/6vVL02aNHfDN5lmbhADkfxkaHRJhnQIzW4fHR6d+2Dh6CQaVbfG\njetu3bx+9fIZ5rWRlYmtIUAbiOQn8gtkn05cEpESun20kQ88djWe6/e94pKR3YvZApH85FXv\npYgTSCJSQrePNkke6BUHkUiASH5ywB6ilhrv331it482/oh07ZPY8fkd2APdCUTyk95R1lIH\nycis20eb8+fPp5p/ZLNs5Zp1Gzdt9iIh9usZznWRDwsxyWtvZ65AJF/pCInaR6kCUhQJJpFu\nH6ph4gPYhiUs65liWOzc3HZ9DUXb/uAAIvnL/pCo/1uK6+lE2mlHUYxuHcA2LGvMl8Qaz02R\nPHbm5ApE8pkz48SwtReTDhLS6/ahEin7bVjaIy+JrVj5eWTNvuPZPW7wAZF8wO6LqZlc12BO\npCga5TVoVaHV7UMj0gC2YTFXZB3xQZczqUMNcKjCWLsIEMA1WY0AABzjSURBVMkHUhYyXr97\nkabbR3uYuBPKbVjUps0T4xcofjhVOOed8wYi+UA2IpEva5XRNiy9V47u3LIueYP2KnlnYvOj\nam3EV1IEiOQDQRApg21YHh2cF1nYLrKI/tkN9xJCJgkx2+Ml9UJMIs53oQKRfMCzXyZVN41u\nt48++tuwhI+Msa5GRForRiZ0EJUL8a7HW1o0N21mAEQKMNl3++hvw9LfHNMsItIEIULxXz5F\n3vP/WrGKkAVECjDZd/vob8NiLcofskS6oc5KutxRpd6DBluwipAFRAow2Xf7aG/DctKUqPGj\na08skfqPqf0C49Y0qRYisVnQMOYJgekWESCSnwxowt4Aun20t2F5OlZer78eDYxufdk5Urrl\n3kRjkRAjEreyuCmLsYXZZXDQAZH8JElrnOeEvZTdPhvT43qY7jYsqkCa22cFWpsxnxD2Ut9R\n1MZ+CXW7PrWv595Ufz8jIJKfJBHJa8Je6m4f7zY4F66n6W7DskqIIbftQEuk8EQhlrri7g2R\nhdQO9xDVW7NVc/rvBlBAJD9JIpLHhL003T6ZiqS7DUutEMtiebVEMt4Rosadvw3qFRP3XLP2\ne3pwqtlc1m6rxr8CCyCSn3iL5DFhL123T6Yi6W7DUi7EgdgrbJEOJ/QPPZoUeUtxdX3Dy3Nq\ny6MvfWHAS/YNFiCSP2Q2YS9tt4+6XL7w3XWrm1dojbXT3YZFllxfxF5hi9Qhy8K4P6jLazOA\nmm4DRIBI/pCu/HCv05i22ydy0+jVpzRLAN1tWIYKYe1h6xTpuEf/0N1FCX/E0rt6ueEARPKH\ndCK5xwmk7fZpLo3eN3zFaZ1ZqbrbsIy1GyVcIrV57jJx7lXn5KmhCxM6pTgDkfwhs3Ua03f7\n9B6ZYz1xzMZraV+vuw2LfPEKw36zLdJLQtR7PfZ++6ZF9ZOqp8xasq2DeHBgoQORfCKjdRq1\nun26tlVZD5j+YeIAbTea27CsF6Lkd/vNlkjXioTYkOkfzByI5CdCc51GzW6f8Hm7ile85IuU\nc1M1t2HpkCkL+q03R0XqqZXHX2X8x/IGIvmJ0JxnpNvt46riVbSk2gpCbxsWcwr54gfRvEZE\n6lQejX+mk29gA5H8RFck3W4fk66tdhVv6u7k7WZ627CosT+ifHunJVLvmebixIILpAUi+Ynu\nhD3dbp8o4W9W2lW8BZ8lKzu0tmEJL4xcKXlB/k9NXVW0uGvEgnUZApGCgG63TwxHFa983fdp\nHp9qG5aeWR6tIfW9kYtUC1MyACIFAd1uHxeOKt6kNp2tyb3pWx/fUh/aaDUqCuHXtuuDDogU\nBHS7feIIf9NstSQUzc/+7T+/XuQQY9iq2HJBEEkbiOQDszRw36HZ7ZPIszPrqgf+g75/YkNj\n7fPVtQ1L3z/r/KiDSNpAJB9I8bNL8vvT7Pbx4MfdjSWJD1w4efLkHwf+h2S1MCVPIJIPZC6S\nZrdPPJfXVVqRo10XlFvZfzaBzIFIPpCFSJrdPk5+b5tkxRUvOeU2Ton0QCOnr59BOzcNEMkH\nmjVIuEmr28em74uFtnDT9iaMvKuRyR5rbyWghsASVAEBRAogqbp9InS+a4+HLV/vVWqpCYJe\nKzrGYz5i6q60tcABLYfEAohUcPQcmmFZVNR40ntcw1357TTyhuclF9ZzXjn+OF2c9nJIPIFI\nhUX4/MphlkapumG76mSt7Uzax9nPEqXN58LJ45KI5LUcElMgUh7Iek3v263V1g+/7K3Uw2HD\n7XNDon7P1dQlTW/74qG2S5Wbk36TJRHJYzkkrkCkPJDdmt5PP51vDeYJzfsk9eINrZvXr14e\nqQGOrPTeKtCi52hjrJ1w+j7v4eTeInksh8QWiJQHslnT++p6awksUd16O110Bu3ukvuH59mD\nhIob212OZrYcEl8gUs5JtqZ3+Ma5z48nELm2v9au0q3SaSfLTCTJ3QP2aHJRtvpiWOdJJrXe\nz2MHRPIbza1cP19YkuJ3Hz0pbz6VuJS9FxmLJPltT70dVLX1l/RPMtFpZOcARPIZva1ce+en\n/t2ro4rlR366+u3FC+e9cT1v2co16zZuSr9VYBw3/1Ib9+IMl0PiC0TyF82tXJN65C6RUkGU\n4Rs7Z4Zcz8toOSS+QCRf0dzK9bR9pSieSEDuRJLc2zUi/nlCczkkvkAkP9HdynWFCnvlWOej\nZH2iuROp/9+21IYSnifot10fZEAkP9HdylX59WGqB41LD0V+n3S8PtJTTIiUDojkJ7pbuQ6X\ntb3c5syDx581WWsTidCsA65russh8QUi+YnuVq7F3puGD5zeK0d3bll3WSPw5DLbIjHpg5u+\n5GYwA5H8RHcr13Ih9tG//dHB6ICFyAqqZzckWzC8t31pbPTquHdTzSYESYBIfqK7letsHzo2\nw0fGWG5ERForRn7pEdfbviTWFTxi9fkUY8BBciCSn+iu6d0qxHLiV/c3xxr0IiJNkJ8+J+Oi\nHn2yOGbR0CUdcd1CWCBSG4jkJ7preneXiOHEu9+9bjUbWCLdUGclXe6o2ByKUMNfE5sThMBy\nXJpAJD/RXtO7TYgFqRcMypCTphyNH117YonUf0wNH49bJ8+yoXbXb15PgUjaQCQ/0V/Te50Q\nDb8YZDwdK3/i9WZl0q7aGZ2ygAy5Z2CYKlT9OdmEPoikDUTykwzW9G4fLULz91243n3HTXYv\nVgXS3EgRFxPJOCHsqmYUIUatvZS8eQELRGoDkfwkgzW977xH+X/4q4QYEi17HCKFJwqx1BXX\n9CVphZIxEMlXdNf0Du/0no2UrUi18mXRQ4dIxjtC1GT3QJAGiOQrumt6byL+BCkXwhrj4xTp\nsPDeuQwMGIjkL3prev+YfPZcdq8tFuKL6KFTpA5Z4/OM1x5KBJIAkXxGa01vc/pf+Rv7Os6m\nnviqzVAhTkcPnSId91ysRGMoEdHuFoMYiOQ3Omt6T5WJTY8IXzrWHpvkEqnNo5VDaygRdrdI\nB0TKHcnX9JbVvvHJ9lTOijohVkQPnSK9JER9fKa0hhLp7m7BF4gUBGTNaj3pA9cLURIdJOsQ\n6Zp8z4a4SL2hRLq7W/AFIgUBWSLtIX1gh/zhL4hsfRQTqUctEfSVO1BzKJHu7hZ8gUhBoDZ+\nxMFA6VM7+S02K2O2SJ3Ko7gapO5QIt3dLfgCkYLABiHeon3ifrMdcHunJUjvmeZij/4r3aFE\nurtb8AUi5YJ03TTXi8Ro2pUWwwsjnz4lL8j/qamrinZUNcbtdKk7lEh3dwu+QCTf0ZnxvUWI\nFtq39szy6N2t742L0h1KlMHuFkyBSD6jN+M7LIuGFtoyqW99/HCJ0MaEdVF1hxJ5SGlBmunC\nBSL5i143jaR9rBi/4yppZ9LPrxc5fvDDVv2QGKI7lAgipQMi+YteN822jWubG9WF4tHP17gZ\n0Nvvn9jQWPt8dW3D0vfPeq5LpzuUCCKlAyL5SoYzvnP+O9UdSpTl7haMgEh+ktGM73yIpD+U\nCKQGIvmJ/ozvPImkP5QIpAYi+YluN83rq97e8K5nxWmzvxnUHkoE0gCR/CTwM751hxKBdEAk\nP8n1jO+9qdmze+eOTe47NIcSJfAUi6a4gUh+kumM74GS4lsryUeX1lCiX5uvG3H8peaUL39B\nwQKR/CSjGd8EZC6SxlCivq1F1p4aNr0VQrwRP96INRDJTzKY8U1CFiKlHUpkmjYp7qa/qLip\nvxvAAiL5SYbdNH1XP9781vJFK9Ztbe/yup6WUgcjrY3DiuyV8sfW1M6NH9VtpBlK9KjeTI2b\nHvuimVhDvPB/IQOR/CSjbppLq4c7fs9jt91ODMmEI+UiNHP3ue4nhtHfc+3TlkohKj3G+Jmk\nGEq0wsxN6Wl36tNWs0miIW5WBmMgkp9k0E1zrym+Dla8cyB7fr0tSz3Xol/h47KiuTHTx3xu\nZuX1xJkf35mD2j8YQA4HFxDJT/S7aW5XenzOvJww7UGb3ULMj7/7VrUQezN7THiqykeb16Uu\nleNSrNEVBSL5im43Tc8kS54h5RUj7aUkV3g+VIO7JV6/8fPyCyjZPrLenFG5SDJa6Kradzbj\nIm6wApF8RXPGt7HVrMu9cvCHyBfK/W/3RL7xT2T53l1CNHskT0u3yMqzuH7W1TIPzycb5NAm\nL1agYzYCRPIXvRnf91Uzw3x3S92VOuVell/zr8Ta3Z1skG9xnPa7iqcn+18aLkKjXz3imKmr\nPoQ+SfaSvnHyKhZEiQCRfEZrxvcnMrkpvmnhSYNM/Tq7t8qfuNfIA1lQVdon97eMca4TcWGc\nlcFq+6W/ybORyT/UVDm6NbsMDjogku9ozPiWNagRiSsC/zZUiM3ZvXOot0i7ZSXTOj5eJkRD\n7NKZIQ7VP40mfq3aSpK/5VvVIJJdBgcdECkHpJ3xXSvEao/kJnsL2kyRIn3okSyrdkOjhzuV\nMbGhs7+NcBaaxdEd2P8qUrZwPysWYmJ2GRx0QKQgUCnEXz2SDwgxIbsHjhfiFY/kGiHGR45O\nRyqcdsteZI2WiU1rXjFLpimRNoS9Hg2MTiqEKM8ug4MOiBQEhgmRuESXYZwSoiy7B84RXlPz\njsjUWeZRX7XZlrjaGuNzW7W4Fx9Rn2l3zZ7hw2ZyqzxKNhpCMTFxe3auQKQgkOST5mzW0y12\nSQNK2+MSP1aD7lrNQ9W6IeZ125feFY6yR5VOkSrlzjQl0ljHNxdzIFIQGGsVAW6OR7fNzJzb\n5t7O0/d1Ws3n4esHas12hMjMoiWqbHJ0AY2zyyojMkdCmG3xH4qUm1D0FWNPWguIFASmCrHK\nI3mtEFOzfOIHVrNB9YtzX55bX221yUW2YeofLkTRL7HoH4SrDrfdKojUWME1yV/yo7xcm2UG\nBxsQyQ+epMcV/4YQIxObv3vKvf3Soc9amdLNnEhXcJc8bHRE71A1wdgAhquWcT8Ju3XCizaR\nuNg+VyCSH3j+iN244mUdTqyI75ANrxAphhWkI7y9JOGdoZVRf0/Lk4OOYNU24egvCpdK5dRB\nv5rTFLdxtCNMjXui3depcIFIfpCpSOYQoSb3NLk76jtmxMPsM3FrbYXrjSVLvrcuHZanjk1m\nninnnP1O04V43jx4NVWZc0g9FXudR4BIfpCpSJFBqyPWnre6a3vPrjInuLYOKBv9F/e8Obd2\n4tgJU2Yubz3l2Nlon3x0rMnOuCziCp5GIUabB6rBXFz0fvq9USJxDjpbIJIfZCxSzwvR5HG1\ns195qa4yOj6vzqdtvVRThGOIhVqBocxZs1xmNWs/VguVj7vp9YweNajWVUFkDUTyg2Y3qv8m\nLil+ksPNMR6yVfu1uogSyTECXW2E4RoIsdTuaG018+GxgsQdc1j7+OynHg4yIFIOKE8ogRLp\nnp/g0QrPcXkUqKpdbN5feER8JVLmZUTk6Ik5BKJ0f3xLyIlyM4tY3M4CIuUAHZEM4+RLzgkX\nxa9d8C9D6tPnmn2mGrmF622y1jYuevhdZLru87sdc5d6DkYn9K73L4uFBkTKAXoiyfrd8Xde\nrqsZV1P/6qYOkpWukq3vdV5mKDaASA0EGurs2FLN33ZHa3tU79CkVbuOdJz6ePfb0yzhX8P0\nWBuIlAN0RSIm+fpe3fJ8nX2mxg45J84a12RCk33WUSq8WQOPYkCkHJAXkVKu7zVGiAqrneCi\nuuYa6qeGCO2InXbWe2lUke16EoMTiJQD0op0584dr/933795c5I+nPSkXt/rzVjzQp8qkIqd\nNcknY+K+mcKfTI5/Uvn2AXQVD0YgUg5IK5K8/q1H8jLvtYB0SLO+1xn10fOROnpqFlyuIX2b\nlChxq65c3jgl1hRStaoD2yfFAZFyQLYibRRiSpavTLe+l9mbWrtl9zvmkidFvzpuVW3jji8o\nm8dXOg7t23vw08v3s8zToAYi5YBsRVqf9QzZtOt7XXQuyOLU5qq5Et/QAS48zg+IlAOyFKlv\niiwrsntj+vW99jk8mhFr+94USRnYGD+OQKQckEKkKyby+sdX4rjcrsY6ZLkhmcb6XkftVu3F\njhEUB8yUVweyfD9PIFIOSCGSR9OaiywnoOqs73XzvYnFQgxf4For9Tv10iaMoMsYiJQDBiBS\nihUTUqG5vlf49wdxZc+TYjH2CMqjzIFIOSCFSJPjFzR2M+NJkvvSkP36Xn/+LMtXMgci5YBU\njQ2Pzn2wcHQSjapbs61jka/vBdIAkXKATqvdsavxXB9Afw35+l4gDRApB2Tbj5Q9U8nX9wKp\ngUg5IPci0a/vBVIDkfxglhs11C0uaZYr/vz587SzYX1Y3wukBCL5QbpG7XQl1IDxZX0vkAKI\n5Ad5F8mn9b1AUiCSH2Qp0u1Te1reWLqoafWmj755NLAc5Hh9LwCR/GBcehLuub/NNXuueO6h\nAc35ye36XgAiBYS9wxN+9pUDWuwqp+t7AYgUCMJvelX/QlmOtIuSw/W9AEQKBB/av/chZRVl\nsQ3GOwb2WPL1vUAyIFIQuDtMWVO6+tjPkcrXw2vH15SppHJUxgoEiBQE9itpVrul6d2gKmb7\n8pQjkCEQKQjMk8psTEhtE3Fr24PgApGCQLUQFR4L200TYmzuMwOyASIFgVLvPY9lkTQ053kB\nWQGRgkCxEH/xSP4M8/AKBogUBMq9WxXOZDwPb2N6SDIM4oFIQaBeiPc8ko8J0ZDZg7y6deMg\nyTCIByIFgU3ezXNvZzxaGyLlC4gUBDpDojhxx+PeMSL0t8weBJHyBUQKBLLsmZcw2HudEGsz\nfI45GmLhu+tWN69oSgJRjoEbiBQInjUKMee6K6lnjRAzMp0+FCl0Rq8+hcXpcgxECgb9rSWi\nuPETaxOIh+dbRonQmxlPw2u2FvQevuJ0f/pwQAZECgJ/Xte8ZIopQOnEupl1L1SYxxNenTd3\ndpLlUpLRe2SONXdizMZr6eMBERApCGi0Eei3EnRtq7Jumf7hPT+zDWJApCBAKpJhhM/bVbzi\nJV9g7/FcAJGCALFIhquKV9Fy1adsgxgQKQi8tHDpipXNacj0oV1b7Sre1N2YH+szEGkQE/5m\npV3FW/AZNiL3E4g0uHFU8crXfZ/v3AxiINKgx1HFm9R2J9+5GaxAJAaEv2kui6pUND/fmRmk\nQCQePDuzrjqLxj+gC0Riwo+7G0sgkn9ApMAQvnHu8+MJkDz68rpK6zNpNMkDQTwQKSB8vrBk\n4N2wnvzeNsl6WPGSUxjn4A8QKRD0Jq54TyNS3xcLi61HTduLkXe+AZECQVKPBiZS57ujreeU\nr8dAIT+BSEHgtO1NUTzZP7Tn0Az7qY0nMa7BXyBSEFihfu2vHOt8FL9/ctaEz68chm7YHAKR\ngsBE+XP/kPB5t1urLYvK3vqW8MEgGRApCAwXoo7sYU8/nW8NrwvN+wSLN+QGiBQEioXYQvSo\nq+vLrcKouvV2+nhAA0QKAkmWLM6c/bV2lW7VRZInAj0gUhCYLcTAtou1sNq6m089JXke0AUi\nBYFWIZaTPEhZVLH8yE9Xv7144bw3JO8B8UCkINBdIoaTTAZP2q9LNlQCeAORAkGbEAsoRsFB\npHwBkYLBOiEafhn4YyBSvoBIAaF9tAjN33fhevcdN5k9ZVx6fMk9gEgB4c57KEAKGYgUCMI7\nvWcjQaRCASIFgk34pClwIFIQ+DEEkQociBQE1illyt/Y13EW/acFCkQKAlOlR02P8p0LMAAg\nUhAoE2I8prAWNBApCBQJsT7feQADAiIFAVki7cl3HsCAgEhBoFaI3fnOAxgQECkIbBDirXzn\nAQwIiBQErheJ0VhcoaCBSIFgixAt+c4DGAgQKRCEV0mTUCYVMBApILSPFeN3XEVnUqECkYLA\nto1rmxvVMKHi0c/XuMl31oAeECkIYEJrwQORggBEKnggUhCASAUPRAoCr696e8O7mz3Jd9aA\nHhAJAAIgEgAEQCQACIBIQWA+dqYsdCBSEBCivOXHfGcCDASIFATMhu7p+x/kOx8gayBSEIj2\nGQ1tOku2GzPILRApCIy1+1/Hb+3Kd2ZANkCkIBC+sLbCUinUcAzzKQoPiBQQ+s+uGmm5VLbm\nu3xnB2QIRAoOz75cMdxyadKu3/OdHZAJEClQPOlYOiyqUvGiUxR7+IHcAJGCRu+JRUOjLo1+\ntzPfuQGaQKQA8vDwHGt7iplH0fJQEECkYHJr55SoSuXb7uc7MyA9ECmgPDjygtXwULrtcb5z\nA9IBkYJIz7HGYuc02crL+c4RSANEChx9ny+xWhtE3bbpkSa8o/nOFUgNRAoY1zbYYxzK1l6V\nCRfnqeOiL/KdMZASiBQkHhyos6tzMz+2vozay+RpBZocAg1ECgzhM012la685WfHlZ/KZVJb\n3jIGNIBIAeHX9yrtwmjuJ0/dF7+QiTPyky+gB0QKAr1H5tgWVWy6kRhQL9Nzny2gD0QKAqX2\nHIqXT3oOsNsuL+U6UyATIFIQiGo09r2bSQIOYs3VgAORgoDZwN14qj9pwD4hRucwPyBjIFIQ\nEKJqx+1UAT8dP346V5kB2QCRgsCSM1j0pMCBSAAQAJECSf/N6w/znQeQCRApKNxu67AOb61W\ng4Jq/oJ1jAsHiBQM+rcMFeuix9+XRZvDp9zJa55ABkCkQBB+TXrTGDm+PcYe5TAN88wLBYgU\nCPYpbSZHjpc5ZvS9n99sAW0gUhB4MlpV5NrN45/Vuif1X9++slZNSUKRVCBApCDwuXSmITri\nu0U5ZU5F2iaPOlLdBoIDRAoCsvAJ/S1yGB4n9TljHj4pF3YDBAg4ECkIzBRiZvTwivSoKjrO\nYbkQc/KWJ5ARECkIjBFia/SwTYr0TvRY1u3G5SlHIEMgUhAYJsTB6OF8KdJX0eMPhSjLV5ZA\nZkCkIBASItJkZ/QPF6K4N5p8XIiivOUJZARECgKlQkQXrvtWOL6L/irEsHxlCWQGRAoC8htp\nZ+RINXnvtJJ3YqWGggEiBYFZQiw1D8I1UqRrVvIKIerylieQERApCKyVdbi76qBDLfRtpfaN\nFqIpf5kCmQCRgoBat27efcO4qpYr3mSlyk8kcSCf2QL6QKQg8HSUdGbY7Do1zK64K5p4qSw2\n3gEEHYgUCA7EBny3RFJuri6SJ4vzmy2gDUQKBOEmy6P66HjvS+qktCv1bSAwQKRg0LfdXEA/\n1PQomtArq3klWIKrYIBIQaGnvXXL/uux8+dF/ZX85QZkCEQKKscvYa27AgIiBZt7+c4A0AMi\nBZjwN01D8p0HoAdECix3dlZjD4qCASIFk/BXS4rN5vB8ZwToAZGCyO3t9j6Y+c4K0AMiBY6+\nzxtD9g5+8/OdG6AHRAoYXVtiC62O355y0yQQICBSkHjWPs8ujIqxaVIhAZGCQ+fGUXZhVNN2\nN9/ZAZkAkQLCk6OzHWt+X8x3dkCGQKRA8OO6MluiErTVFSAQKf/0HnrRtqji3etNEKkAgUj5\n5rs3Sy2Lil7t6DMMiFSIQKS88nD/VLswmtD6m5kGkQoRiJRPVpbYH0bLz1mN3RCpEIFI+cTS\naPr+B7FEiFSIQKR8Ylo0aXunKxEiFSIQKZ8ojxrO9bkTIVIhApHySaRiN2qta1Y5RCpEIFI+\nGWJ9JFVtvWknQqRCBCLlk7ttE+0JE/Pao7sxQ6RCBCLll/C5pcWWS+Ut5j4UEKkQgUh553dz\nbYYI9Yd6IFJBApECQPjr1+xiqfSNGRCpAIFIweDOBxMcsyggUsEBkYJC+MwSu1gSC0/35zs/\nICMgUoC4836VrVLltlv5zg7IAIgUKMJfLbKLpZA5qQIUBhApaHS3jreLpTFbsEFSgQCRgkf4\n9MIiu1jKd2aAHhApkPxmL7Wa75wAPSBSQOn/ckERRCocIFJwub2tEiIVChApyPSfasx3FoAe\nEAkAAiASAARAJAAIgEgAEACRACAAIgFAAEQCgACIBAABEAkAAv4/U9+p/uVsBs0AAAAASUVO\nRK5CYII=",
      "text/plain": [
       "plot without title"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 420,
       "width": 420
      }
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "ggplot(data = test, aes(x = Status_on_day_collection_summary, y = erythroid_standard)) +\n",
    "geom_bar(stat=\"identity\", fill = \"#d37d7d\", width = 0.5, position=position_dodge())+\n",
    "  scale_fill_brewer(palette=\"Paired\")+ theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),\n",
    "panel.background = element_blank(), axis.line = element_line(colour = \"black\"))+\n",
    " theme(axis.title.x=element_blank(), axis.title.y=element_blank(),\n",
    "      axis.text.x = element_text(angle = 90, size = 30, vjust = 0.5, hjust=1),\n",
    "      axis.text.y = element_text(size = 30))+\n",
    "theme(aspect.ratio = 2/1)\n",
    "\n",
    "ggsave(\"erythroid_enrichment_CD38pos.pdf\", width = 10, height = 10)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 62,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0gAAANICAMAAADKOT/pAAACc1BMVEUAAAABAQECAgIDAwMF\nBQUHBwcICAgJCQkMDAwNDQ0PDw8QEBARERESEhIVFRUXFxcZGRkiIiIpKSkzMzM6Ojo7OztA\nQEBBQUFDQ0NHR0dISEhJSUlLS0tNTU1OTk5PT09QUFBRUVFSUlJTU1NUVFRVVVVWVlZXV1dY\nWFhZWVlaWlpbW1tcXFxdXV1eXl5fX19gYGBhYWFiYmJjY2NkZGRlZWVmZmZnZ2doaGhpaWlq\nampra2tsbGxtbW1ubm5vb29wcHBxcXFycnJzc3N0dHR1dXV2dnZ3d3d4eHh5eXl6enp7e3t8\nfHx9fX1+fn5/f3+AgICBgYGCgoKDg4OEhISFhYWGhoaHh4eIiIiJiYmKioqLi4uMjIyNjY2O\njYuOjo6Pj4+QkJCRkZGSkpKTk5OUlJSVlZWWlpaXl5eYmJiZmZmampqbm5ucnJydnZ2enp6f\nn5+goKChoaGioqKjo6OkpKSlpaWmpqanp6eoqKipqamqqqqrq6usrKytra2urq6vr6+wsLCx\nsbGysrKzs7O0tLS1tbW2tra3t7e4uLi5ubm6urq7u7u8vLy9vb2+vr6/v7/AwMDBwcHCwsLD\nw8PExMTFxcXGxsbHx8fIyMjJycnKysrLy8vMzMzNzc3Ozs7Pz8/Q0NDR0dHS0tLT09PU1NTV\n1dXW1tbX19fY2NjZ2dna2trb29vc3Nzd3d3e3t7f39/g4ODh4eHi4uLj4+Pk5OTl5eXm5ubn\n5+fo6Ojp6enq6urr6+vs7Ozt7e3u7u7v7+/w8PDx8fHy8vLz8/P09PT19fX29vb39/f4+Pj5\n+fn6+vr7+/v8/Pz9/f3+/v7///+COmTaAAAACXBIWXMAABJ0AAASdAHeZh94AAAgAElEQVR4\nnO3d/X8U1Z7g8bu7s7Ozs3d3ZtfZnX06nUAID5GnAMGIYMyNyIMoRhAhyhrhIgKCyBUwykUQ\nQfEiohFRRBCREZSLYlAZQBBBiEjS+ZO2qrur+lvdVd2nT59Uh+Tz/mFud1Wlu/CVz1SlHk79\nrg9A2X5X6RUABgNCAiwgJMACQgIsICTAAkICLCAkwAJCAiwgJEMb/+vWSq8CBhBCMrTydy9U\nehUwgBCSIUKCREiGCAkSIRkiJEiEZIiQIBGSIUKCREiGCAkSIRkiJEiEZIiQIBGSIUKCREiG\nCAkSIRkiJEiEZIiQIBGSIUKCREiGCAkSIRkiJEiEZIiQIBGSIUKCREiGCAkSIRkiJEiEZIiQ\nIBGSIUKCREiGTEP6fyEsrxoqgJAMERIkQjJESJAIyRAhQSIkQ4QEiZAMERIkQjJESJAIqbDz\nH+x8ZdueT67mzSAkSIRUwI8v1qu0xIyd3cF5hASJkKJtr1XC+E8DMwkJEiFF6V2mghJ75WxC\ngkRIUdan6nnwzS+//+bjdWPc18NOidmEBImQIpxIOOnUHcm8uzbfLWmemE9IkAgpwgInnOrs\nJujXaW5J32fnExIkQgr3k7tBek5MeN8NaVf2PSFBIqRwH7ndfCMmdFcFyyIkSIQU7sS69tbZ\nSTllnBPSU9m3hASJkHS5Ia3MviUkSISk6ab7R9O27HtCgkRImva6fzSdzr4nJEiEpOeKe9Vd\ni5hASJAISUuXexppuNggERICCKmY5E9fv9fmHvse8YmcTEiQCKmwRv+a1Ud+CMwgJEiEVFgm\npMkdXWLi1bvvvlv9Z0JCFiEV5m2RRj+051d/4qU7HP9ISMgipMIW3DNtfHW6pYkfyBns2kEi\nJA1XPnzcPdqgXhXTCAkSIen5aqp7j+zh7ARCgkRImi6419pNz17GSkiQCEnXdnfn7qT/lpAg\nEZKui25IW/23hASJkHT1ugfv1vhvCQkSIYXqeWPzykW7A5N6E4SESIQUrs6p5qHAlPPurt0r\n/ltCgkRI4R5xBxG6LKfsdkPa778lJEiEFG6Pm80yMeGWe7FQ1S/+e0KCREjhuscFD9L1POG+\nX5RdgJAgEVKE1CZJtaZPHPV+3OS+q/4uO5+QIBFSlGfTl6qOnbd40eyR6dd/EbMJCRIhRUlu\nTASfRlHzhpxNSJAIKdrJ+SKjqoVnAzMJCRIhFXJ2x6Jp44fXjJ++dOfFnFmEBImQDBESJEIy\nREiQCMkQIUEiJEOEBImQDBESJEIyREiQCMkQIUEiJEOEBImQDBESJEIyREiQCMkQIUEiJEOE\nBImQDBESJEIyREiQCMkQIUEiJEOEBImQDBESJEIyREiQCMkQIUEiJEOEBImQDBESJEIyREiQ\nCMkQIUEiJEOEBImQDBESJEIyREiQCMkQIUEiJEOEBImQDBESJEIyREiQCMkQIUEiJEOEBImQ\nDBESJEIyREiQCMkQIUEiJEOEBImQDBESJEIyREiQCMkQIUEiJEOEBImQDBESJEIyREiQCMkQ\nIUEiJEOEBImQDBESJEIydLuFRMD9i5AMERIkQjJESJAIyRAhQSIkQ4QEiZAMERIkQjJESJAI\nyRAhQSIkQ4QEiZAMERIkQjJESJAIyRAhQSIkQ4QEiZAMERIkQjJESJAIyRAhQSIkQ4QEiZAM\nERIkQjJESJAIyRAhQSIkQ4QEiZAMDdqQKM4IIRkiJEiEZIiQIBGSIUKCREiGCAkSIRkiJEiE\nZIiQIBGSIUKCREiGCAkSIRkiJEiEZIiQIBGSIUKCREiGCAkSIRkiJEiEZIiQIBGSIUKCREiG\nCAkSIRkiJEiEZIiQIBGSoYESku7nEVL/IiRDhASJkAwREiRCMkRIkAjJECFBIiRDhASJkAwR\nEiRCMkRIkAjJECFBIiRDhASJkAwREiRCMkRIkAjJECFBIiRDhASJkAwREiRCMkRIkAjJECFB\nIiRDhASJkAwREiRCMkRIkAjJECFBIqTCrh56a+vW3Qcv580gJEiEVMCtPfclVNrULT8H5xES\nJEKKdrJBCXX7AzMJCRIhReqsVkG75VxCgkRIUY6lOhr/p0+7ug5vGOu+rjopZhMSJEKKcDO1\nX7fi1/S7G4+772aL+YQEiZAi7Ex15L/taXHfd2XnExIkQopwn9PN2O7s+6NuSNuy7wkJEiGF\nu+Z2s05M6B3lTFiWfU9IkAgp3HetzROrP5JTpjshLci+JSRIhBQtmZTvmpyQWrNvCQkSIemq\nY9cO0QhJ0/mcU7KEBImQNL3odFR9KfuekCARkp5Lo52QlooJhASJkLQkH3E6GnFOTCEkSISk\nZa37F9IOOYWQIBGSjpfcjpb4b3u/++67Jf+OkJBFSMUl17gdzbnpT7h0h+MfCQlZhFRU92Op\njm5kp1ybNWvWiP9ESMgipGIu3Ot29NjNnMn8jQSJkIr4InVP37Le3OmEBImQCusc7mSU6Mif\nQUiQCKmg7e4oQjWdIXMICRIhFfKKu1s39mTYLEKCREgFpDpq+CF0HiFBIqRoe92OWq6EzyQk\nSIQUKTUe18zuiLmEBImQolwd73Q043rUbEKCREhR2tzjDBcjZxMSJEKKcNj9A+mj6PmEBImQ\nIrgDQj5cYD4hQSKkcJ+4G6S75uQ64S9ASJAIKVyHCvWxvwAhQSKkcIRESCUhpHCEREglISRD\nhASJkAwREiRCMkRIkAjJECFBIiRDhASJkAwREiRCMkRIkAjJECFBIiRDhASJkAwREiRCMkRI\nkAjJECFBIiRDhASJkAwREiRCMkRIkAjJECFBIiRDhASJkAwREiRCMkRIkAjJECFBIiRDhASJ\nkAwREiRCMkRIkAjJECFBIiRDhASJkAz1d0iVWo6QzBCSIUKCREiGCAkSIRkiJEiEZIiQIBGS\nIUKCREiGCAkSIRkiJEiEZIiQIBGSIUKCREiGCAkSIRkiJEiEZIiQIBGSIUKCREiGCAkSIRki\nJEiEZIiQIBGSIUKCREiGCAkSIRkiJEiEZIiQIBGSIUKCREiGCAkSIRkiJEiEZIiQIBGSIUKC\nREiGCAkSIRkiJEiEZIiQIBGSIUKCREiGCAkSIRkiJEiEZIiQIBGSIUKCREiGCAkSIRkiJEiE\nZIiQIBGSIUKCREiGCAkSIRkiJEiEZIiQIBGSIUKCREiGCAkSIRkiJEiEZIiQIBGSIUKCREiG\nCAkSIRkiJEiEZIiQIBGSIUKCREiGCAkSIRkiJEiEZIiQIBGSIUKCREiGCAkSIRkiJEiEZIiQ\nIBGSIUKCREiGCAkSIRkiJEiEZIiQIBGSIUKCREiGCAkSIRkiJEiEZIiQIBGSIUKCREiGCAkS\nIRkiJEiEZIiQIBGSIUKCREiGCAkSIRkiJEiEZIiQIBGSIUKCREiGCAkSIRkiJEiEZIiQIBGS\nIUKCREiGCAkSIRkiJEiEZIiQIBGSIUKCREiGCAkSIRVxZYFy5E8nJEiEVNihcYqQCKk4Qirk\n15VKERIhaSCkAk41KkIiJC2EFKn35WonoVGthERIxRFSlGsPuAXNPL+TkAipOEKKctrpJ/Gn\nnj5CIiQNhBTFCWniMed/CYmQNBBSlNNqwVX3fwmJkDQQUpSu19P/S0iEpIGQiiEkQtJASMUQ\nEiFpIKRiCImQNBBSMYRESBoIqRhCIiQNhFQMIRGSBkIqJiSk68uXL6//j4SELEIqJiSkS3c4\n/jEY0kAPhJD6FyEVExLSbwcPHpz37wkJWYRUjObfSAM9EELqX4RUDCERkgZCKoaQCEkDIRVD\nSISkgZCKISRC0kBIxRASIWkgpGIIiZA0EFIxhERIGgipGEIiJA2EVAwhEZIGQiqGkAhJAyEV\nQ0iEpIGQIvwwJ2OqG5L35lN/PiFBIqQIp1WYff58QoJESBEIiZBKQUgRCImQSkFIhggJEiEZ\nIiRIhGSIkCARkiFCgkRIhggJEiEZIiRIhGSIkCARkiFCgkRIhggJEiEZIiRIhGSIkCARkiFC\ngkRIhggJEiEZIiRIhGSIkCARkiFCgkRIhggJEiEZIiRIhGSIkCARkiFCgkRIhggJEiEZIiRI\nhGSIkCARkiFCgkRIhggJEiEZIiRIhGSIkCARkiFCgkRIhggJEiEZIiRIhGSIkCARkiFCgkRI\nhggJEiEZIiRIhGSIkCARkiFCgkRIhggJEiEZIiRIhGSIkCARkiFCgkRIhggJEiEZIiRIhGSI\nkCARkiFCgkRIhggJEiEZIiRIhGSIkCARkiFCgkRIhggJEiEZIiRIhGSIkCARkiFCgkRIhggJ\nEiEZIiRIhGSIkCARkiFCgkRIhggJEiEZIiRIhGSIkCARkiFCgkRIhggJEiEZIiRIhGSIkCAR\nkiFCgkRIhggJEiEZIiRIhGSIkCARkiFCgkRIhggJEiEZIiRIhGSIkCARkiFCgkRIhggJEiEZ\nIiRIhGSIkCARkiFCgkRIhggJEiEZIiRIhGSIkCARkiFCgkRIhggJEiEZIiRIhGSIkCARkiFC\ngkRIhggJEiEZIiRIhGSIkCARkiFCgkRIhggJEiEZIiRIhGSIkCARkiFCgkRIhggJEiEZIiRI\nhGSIkCARkiFCgkRIhggJEiEZIiRIhGSIkCARkiFCgkRIhggJEiEZIiRIhGSIkCARkiFCgkRI\nhggJEiEZIiRIhFRY99G9W1/de/xW3gxCgkRIBSTfmztMpdQ+dixnHiFBIqRoPz6ghPabgZmE\nBImQIl2sVwFzf5NzCQkSIUW51eTWM2nTF+fPfb6+zn29TM4mJEiEFKXDbaftRvrNpdRe3lEx\nm5AgEVKEX0Y55czu8d7+PNF52yLmExIkQoqwwwmn+lz2/QfuJunL7HtCgkRIEf7gdPO4eJ+c\n5ExYnX1PSJAIKdzlhNPN+3LKGmfC5OxbQoJESOH2u3tyP8spB90pF/y3hASJkMKtd6qZEphy\nyQ1pv/+WkCARUrjHnGrmByeNdCa97L8jJEiEFO5ep5oVwUnTAudkCQkSIYUb71TzYnDSg86k\nR/13hASJkMINd6rZGpzk7u3N9N8REiRCCpV0jyzsCk5b6kya4b8jJEiEFOqmG9Lu4LR2Z9Ld\n/jtCgkRIoX4NCWmZM6kx9eqnf/7nf/6nfyAkZBFSqB43pDeC055yJk1PvbqSSCT+x395IeTn\nMFQRUjj3CqHXgpMWO5Oa/Xe5u3YY2ggpnHsn30vBSfOdSQ/57wgJEiGFawxe6+1y75h90n9H\nSJAIKdw8p5pFwUl3Bs7REhIkQgq30j+y4LnqHn94J7sAIUEgpHC73RtkAwNwHXVD6vLfEhIk\nQgrX5WbzuZyy0ZlwZ9J/S0iQCCnCZKebtXLCDGfC09m3hASJkCK4d/aNFft2J9xN1OHse0KC\nREgRzrmnZDf6b3tnOW+nZvfsCAkBhBTFvdg7cch7t87dIO0TswkJEiFFuTTGSaeq47r7+vzj\nqcG/xQaJkBBASJE+rnLrqZn7VFuzu5un6i/JuYQEiZCidQ6XD6OY8n1gJiFBIqQCvsk+IKl6\nxfXgPEKCREgF/XXzrMk1oxoXvHk5dw4hQSIkQ4QEiZAMERIkQjJESJAIyRAhQSIkQyt/N3PP\n0HS4+H+cIYiQDL31334X6m9+/7fhM3L8m9//ndZyv/v7v9db7u9+/2+1lvvb3/+N1nJ/8/v/\nED7jf1b6P/2AREimds8KNfmO2vAZOZrv+F9ay8367/+kt9z/ueM+reVG3jFZa7kpUf+Q5yr9\nX35AIiTL9t+xRmu5S3c0F1/I9X//t95yM++4UHwhx7o73tNa7sM7Vul9MVyEZBkhDU2EZBkh\nDU2EZNm1Ly5qLdfzxVm9Dzz1ld5yZ7+4pbXcj19c1Vruly/0wkQKIQEWEBJgASEBFhCSJfP3\n6/2JgsGJkCxRqm7N19pLn3k3+/rYSyf6YX0QL0KyJHUf7X07r+ks++n9qj777i2lpn1c/gp0\nn9q7ZcOq4k3qLoeSEJIlmVvSh7cdSRZZMrnaWW5E9v1b7s+t6M28e62w7du2vLQ+/zNvvDEv\nNVSLejv19sjan8O/O3q5BzRo/qcYkgjJkgn+8A6TNp4ruOSK1FLZQVxTIallmXdKQ+4nJveM\n92alA1mpxoRt4wotZ/K9yCIkS5Kfrxzr/cIl5r5zM3LBzvQy2bG9PqxNTdiXfmfwC93bnp2V\nDmSKsw7787654HKEVB5Csqf3yLIx3u/cqBVfhi/UU+/MHf2K3Gj1HJjqTJvck3pj8Av9ZGZy\nwgvke/ddTd52seByhFQeQrLq1sdLR3q/dtNe/Slkif3OnKm5F9/8MsOZmt401ApjajMfVeWP\nsDehsWnO4vxPVInWN8/c9ALpfcd9BG7OAweLLNeuoez/PIMYIdl288DiEZnf+uqFB3tyZy9z\npv8174dOqexfSVl76lRi5rajl5z9xN7rZ95f42zM6vP22H5z/zprSV245wXS19flbBoTF42W\ngxFC6gfd7y30NiHjnu8Kzpum1KyQH2lxtmC50551fvFPywnJfU4M63KWcjc0c/zdwkwgfe85\nr7cZLQcjhNQ/fnlrdiLT0sy98siDszP1fMjya5S6M2fSNqXm/5Yz7UKDUq8FJzmbuGGZbYoI\nJOn82bXYaDkYIaR+c2HL9ExKdZuydy5UKbUlZOEOZ0cwOOVKjarNG9+175hSI4LniJqUWpJ5\nKQLpe06pRqPlIly/997cbSEEQuo/1/bc7R0jqN30a2ZirVIbQpZ1tki1wSmvKhX21/29ubti\nziZuV+alDOSt3E2c7nIRkgl1j85yQxUh9ZPr77RWy0PH9ZmLcpx9s7aQxecp1RCc8rBSu0MW\nXOvs8AUmON/yUealDOSAsydntFyE75Qao7PcUEVI/aHnw0X+AevmTfelD+HtTc1aqNTo3L98\n+vrOO39QLQhOmqjUwZCPdjZU9YEJzhd5zxWUgezL3cTpLhf+D/qX5rx9T0iEZN+Ztf41DqNW\nuofdjs9LnQxKbRB2qpB9u57ZKu8gwvDwkLYpVROYMCG75ZKBbFVqstFyrl92L2mePC5rTOrA\nSV3kvxiEZNu1Xc3+7tzMt72/jDpHOW/Huoccfh6mVOKl3sDPXHCf9Dws5/StE9LrIZ/v7NoN\nD0xwvm5p5qUM5EGlWoyWc7zuX6AR8Gj0vxqEZFPycJu/S1e35lsx5xv3GoKt7qu17syp2894\nF4lfO9g+zJ20MeezJin1cMhXNCo1KTBhtbOJyjQoAjlTpdRao+X6+laFZqSGnSzyrx/SCMme\nH16s93/r5ryb84fQR87E+90XN6all6huaJn70OymuswP3J17mau7t/dJ3nfscaYGb2c44G4r\n0lu4bCDXm/J+Wnc596+mMM2f6/53GJIIyZLuPbP937mx67/PX6DFmZ56cW5KyK9p46Xc5V91\nptZ25kx8293gdQQmpa6CfSx1Q6EfSJfbx6RbZsulbrWY9tSadc7/PLtuzRPuFbWNBxibqzBC\nssS7vlQlHtqfd4Gda7MzK/3qysK8jhZfyVv+Yo07474dXd7fU8mzu5pSX5AzIJ57+ELVbe7y\nAuk+3J468L6nz2i5D5xJE4+7r5wXqXO/Z5+pUvd824dCCMmSTBETXjwfscAbKnsfwtFHEqKi\n4QuOhf3AK8rbB/zDnIfmtDQMy7xfnbNcckF6eo17+rexeXLmo1t7zZZ7xkn1r94/KXMRxcl6\nVRu6jvAQkiXu72RV68Hc396sHUqNy7672rl+Ycu0hukPLNp04Hr4D/R49w8Fze7OXfD6AyGL\nteQtprlcs1IP+f8k72qkcxNVbcjuKnyEZIlSk18qeD/CN/v2HSo0P19yc03e733iqZB7b3tW\nJ3IXW5d/0ldzOedPpM3eP8kPqe9TpeaWtvZDDCFZsuhwsUFPDFzI3r6e3itbFDEQ+LdPVonF\nRizLv+NJe7ka/8Sts0T23NajSnH4uwBCGth6j2//45ymqROmTJ/5RMfBX6MXvPre2tamuxqa\n5i5++UjEvqLectX+4Qdn85U9rrFfKZ4wVgAhIWiMf3l5rVJf+JMvZE6DIRwhIWiKUuvTr8Z7\n40i4ehJ6t1sMVYRkU1fnlj89tyaXO2eeo8t7EcHsOy9fvhx23mrnCy8cN/q8+UrNTL9qkXtz\nPVz9XRAhWdOzI+ySBZU+e+T+70nvRQSzr1XhRwGW5N4XuOCee+7RGZx8g/OBn6VeLVVqrH+M\n8Eut2y2GLkKy5de5hQKJPaR1Sk0PTHAPpeffup7vhLPcqD1uQLudVyu8yQvDbreAj5BseaJg\nINohlTwGd0RIq50cAhPckLSG+E9dVFt9ta/vonvWaUnqSo0r7iCtuQPlQSAkS77K1FBX35DL\nnd3m+MF7ESH1OQVCC990hYfUM12pqsCURmdBrQvmPna/oso9K9bmvkrc++iSmanLk3Iv3oNA\nSJa410qrtmPRY35rKiGkUynO+7dP5TjROV/l/k3jrmDY+EX5XnCWHO+++GGU/NYpodfiIo2Q\nLHEHZnjFwueUEFKxxZoCH3ylXqkxetfL7Ryp7k29+FiM31J7ushPDW2EZMl4pRqir1jVV8IY\n3MVCytn+nGt2NjSHtVbi+luvpl98NdX7sOauwj8y1BGSJc5fEat0lrN42uee3EtQg+7P3c1M\nds5JqJbtpwtcaJSn95PlD06bNmftZ6Wu3FBDSJYMzwzJUIzSPO2j48bRVxaMi8iooSPnuu6O\nF1Yvf+L+1Lwx9flHQ1AeQrLE+Qtkh85yESHlnfbR5nzgO6dznb0atlwUsy+GREiWPBI+FHGe\niJDyTvtoi/jAsOUIqf8QkiW7/EvUCgv/vc8/7aOtf0I682729bGXeAZ6UYRkSfed3lAHUUo7\n7aPt2LFjhe4/8i15asWqdetfCJO37Kf3y3GR31JqWtiznSEQki0HEqrpRqEFCmwSUtKnfWxd\nJl7GY1iSzn6mGpF9n3rs+gobx/YHMUKyZmdCtfxrgfnFQtriL2Xj6tYyHsOyIvUl2YPnqZBC\nnswJgZDsOTxRjVh5PPIiIb3TPrZCMn8MS2f6S7IjVn6YHrNvn9nHDRGEVB7/XEzjPc1zUzdS\nVN0ZdtGqS+u0j52QyngMS2pE1tGvnJOTDrgXOEzmWrsCCKk8BTcyYb/3qshpH+3LxCWbj2Fx\nH9o8NXeA4l9mKHnfOfIQUnlMQrI+rFVJj2HpPrV3y4ZV0Qe0lzk/mX/40T3ayF9JBRBSeQZC\nSCU8huXGG/PSA9ulB9E/svbnvEWmKTUr5EtalJpmeb0HFUIqT+h5mUKnaXRP++jTfwxLcs94\nb246pJVqTN4Jojqlng/5ljWaD20eqghpYDA/7aP/GJbe9mxm6ZCmKJXI/cunKvz+vw5GESqI\nkAYG89M++o9h8QblT3ghfe++qzkXXKo2/KLBNYwiVBAhDQzmp320H8OyPxVR65tnbnoh9b7j\nPi8wZ0yTBqXyDwv29c1TitstCiAkS8q6Ya+M0z7aj2H5bYIzv+VsZsHMoy+7xjhtBR+isVCp\n0fmPsjjvbMYWmK3g0EBIlkQcjQu9Ya/gaZ91xQU+TPcxLO4GaU6Pt6D3MOb3lD/Ud4b7YL+8\nfbse97merxX69w91hGRJREhhN+wVPu0TfgwuIPBpuo9hWabUsIv+gl5IyalKLQ4s9/MwZyP1\nUvAS1Quz3MPpP/UhEiFZEhFSyA17RU77lBqS7mNYmpRakl1XL6S+55RqDK7fWvcrpm4/4z3v\n6drB9tSwdhs1/isMXYRkSXhIITfsFTvtU2pIuo9hqVNqV/Yr/JDeyjs/dGNa+luqG1rmPjS7\nqS7zpXeXPWTfoEZIZSvthr2ip33c2XULnl+1vH2p1rV2uo9hcbZcH2W/wg/pgLMtzPkHnQt7\nGEDjpT4UQEhlK7b9CI7TWPS0T/qHxi0/qLkF0H0My3ClvGfYypD2hZwfurIw7x+x+Ire2gxZ\nhFS2YiEFrxMoetqnvTbzcyOXHtK5K1X3MSwT/IMSgZC2hj5l4ugj8uap4QvyTkohByGVrbRx\nGouf9uneM9v7xPHrzhT9et3HsDhfvLTP/2Y/pAeVagn72Kud6xe2TGuY/sCiTQcsXxw4KBFS\n+Uoap1HrtM+5TZO9D7jv9fwLtIM0H8OyWqman/xv9kI6U6XU2lL/wchHSJYozXEaNU/7JI/5\nu3jViz4qeG+q5mNYDjhTHu31vjkT0vUm5/UnJf9jkYeQLFGa9xnpnvYJ7OKNXVPoURB6j2FJ\n3UL+2LXMuqZD6nI7mnRLZ71RGCFZohuS7mmflHMb/V28Gduij5vpPYbFvfZH1W3u8kLqPtxe\nnb/hghlCskT3hj3d0z4Zyc+e8nfxHv0gatuh9RiW5IL0nJq7nf/T2Dw5s7lrZcA6GwgpZrqn\nfbLELl7dqq+KfHyhx7BcfyDkaEhLd3qmrYEphypCipnuaZ8AsYs3bavOo8nD9azOPVKfWOcd\nVFSqvx67PjQQUsx0T/vkSH7W7h1JqJpv/u3fPlklwhixLDtcECGVh5DK84CG4E9onvbJd+vw\nqobyf6Gvvre2temuhqa5i18+Iv+oI6TyEFJ5CvzaRfz+aZ72CfH1ttaa/A9ccM8993xd/j/E\naGBK+AipPKWHpHnaJ9eJVfXekuMCM9y2zP9sgiWEVB6DkDRP+0g/bZ3mLVe96GCwODekaxpr\n+uRhjnP3I0IqT7uGvB/SOu3j6/logR/cva/lXXnX6EwOGXsrj3sJrIVdQIQjpIoqdNonret5\n/3rYutVhWy33BsGwER1zpT5ixqtF9wLLGg5p6CKkgez67vu9iqpa94df13DF+dtpzPehswK8\nz3l436/FltMeDgk+QhqwkseeGuFlVOg07LlmZ6/tcNGP8z9L1bYfTUYvFxFS2HBIyCKk/mU8\npvfFjgbvF3/UM4Uvh012zkmolu2nC29pujsfG+63VP9C5N9kESGFDIcEgZD6l9mY3r+9P9+7\nmCcx793Cgzd0vLB6+RPpPcAx9eGPCvRc39uaPU54347wy8nDQwkR55AAABeXSURBVAoZDgkS\nIfUvkzG9T6/2hsBSDR0Xiy1dwnF3x9W35vkXCVW3dgYaLW04JAQQUn+KGtM7+f3RD/flSc/b\n2eTv0i3TOU5WWkiOK7v8q8nVqOXHkzqflNIU/nlwEZJFmo9y/XBBTYHf+8ybuvaD+UPZhyk5\nJMeP21v8hSZv/K74J6XoHGQfsgjJHr1HuXbPL/x7774a+8Seb06fPP75sXCBz1vy1IpV69YX\nf1RgjvN/bsr54hKHQ0IAIVmj+SjXyI6CW6RCLK3w91tmJgKfV9JwSAggJFs0H+V6yJ9TlSu9\nQHwhOX5+dXTu5ynN4ZAQQEiW6D7Kdam72MPvdN2IOicaX0i9/7KhKZH3ecr+Y9eHAkKyRPdR\nrm5frxf6oInF2VjfmweeHBMaJiEZISRLdB/lOtLZ24t3zUL8+kGbNzaRSjywKzBPdzgkBBCS\nJbqPcq0Of2h4+bpP7d2yYdUJjQX3L/ErUtNeOd8vazPkEJIluo9yrVNqh/1vv/FG5oKF9Aiq\nR9ZGDRje3bk4e/XqxOcL3U2IUhCSJbqPcp3VDyc2k3vGe22kQ1qpxnwcslx356LsqeDRy48V\nuAYcJSIkS3TH9O5Q6gnLX93bnj2glw5pivOnz/6cpW68+1i2ouGLDuScFmKAyPIQkiW6Y3pf\nqlEjLT/97knvsIEX0vfuu5pzwaWy91Ak5v4l/3CCUgzHVQ5CskR7TO+tSj1aeMCgEu1PxdH6\n5pmbXki977iXj+eMk+fV0PTqj2GfQkjlISRL9Mf0XqXU3O/6rPltgvMr3pLamfR37fq6nA1k\nIngHRiqFyX+KuqGPkMpDSJaUMKZ35ziVmL/j87OXLgeZfbG7QZqT3sRlQ+p7T/m7mhlK3bny\ni+jDCwwQWR5CsqSEMb0vv2jz/+EvU2pYZtsjQkpOVWpxYLm2j63uUCKIkGzRHdM7uSX8biTT\nkJqcL8u8FCH1PadUo9kHwgQh2aI7pvd6y3+C1CnlXeMjQ3pLhT+5DP2DkKzRG9P76+i758y+\ntlqpjzIvZUgHnD2+0OW1LyVCKQjJHq0xvVO3/9U9vePAkcI3vmobrtShzEsZ0r7QwUo0LiWy\n9HSLoYaQLNIZ03uGM7HthsUvneBfmxQIaWvIUQ6tS4l4uoURQuoX0WN6O7t9k6KeqWykWaml\nmZcypAeVasldKa1LiXSfboEAQoqZs2e12uoHrlaqJnORrAjpjPM9a3OW1LuUSPfpFgggpJg5\nW6TtVj/wgPOL/2j60UfZkK67QwR9ElxQ81Ii3adbIICQYtaUe8VBuXrcJ/k9ltoZ80PqcjvK\n2YPUvZRI9+kWCCCkmK1V6hm7n7gzdRxwc5cXSPfh9uqQ81e6lxLpPt0CAYRkWbHTNGer1Di7\nIy0mF6T/9Km52/k/jc2TMyeqWnOedKl7KZHu0y0QQEg26dzxvUGpNXa/9foDIWd3W7pzltK9\nlKiEp1sgi5Ds0bvjO+lsGtbY3Sb1rM69XCKxLm9cVN1LiUKi9Fhd6UGGkKzRO03j6JygJr10\n2urJpG+frBK/8COW/TV/Ed1LiQjJCCFZo3eaZtO6le2t7ozqcXc1BpX17VffW9vadFdD09zF\nLx8JHZdO91IiQjJCSLaUeMd37L+nupcSGT7dYqgjJEtKuuO7EiHpX0oEA4Rkif4d3xUKSf9S\nIhggJEt0T9M8uezZtc+H7jj1856T9qVEMEFIlgz4O751LyWCEUKyJO47vl8rbPu2LS+tD/6E\n5qVEeX5j0BQNhGRJqXd8l6vA31oRf3RpXUr0Q/vZvhx/bjzYL/+CwYWQLCnpjm8LSg9J41Ki\nno1V3jM1fN1jlXo693oj5CIkS0q449sKg5CKXkqUKm1azg/92V1uxk99KIiQLCnxNE3P6bdf\neOaJhUtXbew8Fza/qFphjPfgsCp/pPwJjU1zcq/q7ityKdGNltTUnNtj/5Ca2Gh54P9Bh5As\nKek0zRfLR4rf5wmbLuYvUoo9dSoxc9vRSzf7+nqvn3l/Tb1S9SHX+KUUuJRoaWptag8Fp/7W\nkTokMTfnrgwEEZIlJZym+bktdx+seks5z/x61tnqBQb9Su5zdjTXlfoxH6ZW5cn8Oz++TF3U\n/koZazgEEJIl+qdpLtaH/DnzUN5tD9q2KTU/96cvNCj1Wmkfk5zhrsfWsFnn3DWuZYyuQgjJ\nFt3TNNenefEMqxs7xh9Kcmnoh2q4UhP2O37M+Qso6jmy4Q67axFxtdBp97mzJW/ihhRCskXz\nju++jal9uYff+Gv6L5SrJ7en/8Z/z/B7X1WqPWTyvcUGWbmVc551ubMOd0Vd5LDVmTmWE7MF\nEJI1end8X3UPM8wPHqk71ey2Z/jX/MPZ4+7SWudbxNvewObp5s4HR6rEuEf2iDt13T+E3o36\nkp6JzlwGRCmAkOzRuuP7XWdyW+6hhZtznamfmn2r8yseduWBs6Gq999c3TBejhPx+URvBRv8\nL/3ReTcm+g81dzu60WwFhwZCsknjjm9nD2p0/ojAPw5XyvDq7+HhIW1zdjK91/tGKTU3O+vw\nMJH6+5mJn7rHSqK/5aR7QMRsBYcGQrKr6B3fTUotD5nc5j+CtlROSK+HTHZ27YZnXm5xi8le\nOvvjaLnRrM48gf0vquAR7lvVSk01W8GhgZBiVq/UX0Im71JqitkHTlLq4ZDJjUpNSr86lN7h\n9I/spcdomdq24uHUlml6+hjCayEHGKWxStWZreDQQEgxG6FU/hBdfX0HlRpl9oGzVditeXuc\nqQ+kXvU0pI4lLveu8bnoHnGv3uP+mXYldWb4rdTkDudV1NUQrqn5j2eHQEgxi/iT5ojx7Rav\nOgXUduZMfNu96K4j9dI9uqHmXfJnPa/EtsfdOqV3KbcU2SJNEH9zIR8hxWyCtwkI2pd5bGbp\nLqae7Xzfji7v8Hny7K6m1HGE9J1Fi9xtkzgFNNHfVvWl75FQqWPxr6uCD6HoqeaZtAURUsxm\nKLUsZPJKpWYYfuIr3mGDhj/MeWhOS4N3TC79GKbekUpVfZdd+q8qsA+32dsQudcKroj+kq+d\n2U2GKzgkEFKZbhYXWP5ppcbkH/6+Xhfel44eb2TKoNnpU8HnnJetYumX3D3B7AUMp73ivlH+\n0Ykw7qUNIbdlwENIZQr9JQ4KLO/sw6mluSdkk0tVgcsKiklursn7zsRTmX4POW/eEAu7xybE\n+aJkrZOc+6LXvacp58HRYjH3uie7z3UaZAipTKWGlLpEqC14m9xl9++Y0b+Yr8SFlWMD31iz\n6Ctv1lvOW/GQmVtuc/K8031K3ZV68Uihbc5u91N51nkBhFSmUkNKX7Q6euUx73Rt95FlqRtc\nO8pajd7j2/84p2nqhCnTZz7RcVA82WiH89HZQ3Z9J1TOhqdVqXGpF+4Bc3U8/NN/vlPl34MO\niZDKVHJI1+/OTJ7YNOvhB5vrM9fnNffTY73cQxHiEos/O29HyT3LJd5h7V/dgconng/7jOvu\nRbWBHUTkIqQytQe5529yJuXe5HB+fH5rqqG/RhdxQxJXoLsPwghcCLHYP9HakVqPkBEkLqcu\na59kfuvhUEBIdtXlbYHyXZqf19HS0OvybHB37bL3/SVH5+5EOusyOv3qZuoSiNqduUdC3qtL\nrSKD2xVESHbphNTXt/9BecNF9eOf998KuX/6nPHfuQe5VeDbnL22iZmXX6Zv171rm7h36fob\nmRt6V/ffKg4KhGSXXkjO/t2+5x5qbpzY2PLI+gNWRrqKGt/rmLNC2QuI3AuBhssTW+7hb/9E\na2cm78S0Za/uOXDw7W3P3usF/zi3xxZGSHbphmRZ9Phel5z3q/x37rVD8sbZvjPOhDb/3YFa\nFW4FHRVBSHZVJKSC43uNV2qsd5zguDsvcKmfe4nQS9m3XS1hGY01HU9iCCEku4qGdPny5bD/\n777zhRcizuEUV3h8rz9mDy/0uBukarkneXN8zt9MyXfvyf2kus1lnCoeMgjJrqIhOfNPhkxe\nEj4WkI4i43u5w2wl3nRf/ZbacAUu6VvvhpIz6sqJddOzh0ImLzvA45N0EJJdpiGtU2q64VcW\nG98rdTa1acO251JDnlT9IH7UPTYu/oLy/XrqwO4dr73x/omrhus09BCSXaYhrTa+Q7bo+F7H\n5YAsMpvTqZH4hpc58DhSCMkuw5B6pjvbCrNvLD6+1w7R0f3ZY9/r01PKu8YPGYRkV4GQTqU4\n898+leNEp3utg+EDyTTG99rrH9V+TFxBsSs15ZFyhu+Hj5DsKhBSyKG1AMMbUHXG9zr/4tRq\npUY+Ghgr9Uv3S9u4gs4OQrKrjJAKjJhQiOb4XsmfruVse25Wqwl72B5ZQkh2FQjpntwBjYPu\nvxnxc0WYj+/1pw8MvxL5CMmuQgcbbhx9ZcG4iIwaOkz3sayP7wUThGSXzlG7d07nOlvG+Rrr\n43vBBCHZZXoeydwM6+N7wQAh2RV/SPbH94IBQirTA0HupW45kx4ILH/s2DG7d8P2w/heKB0h\nlanYQe1iW6iy9cv4XigVIZWp4iH10/heKA0hlckwpIsHt695evHCtuXr3/zsRnlrEPP4XghF\nSGWaWFzez1zdFLh7rnrO7rLu+Yl3fC+EIqT4vTYy79e+vqzBrmId3wuhCCluyT+G7f4lDK+0\ny4hxfC+EIqS4ve7/vg8bNXZU9gHjB8r7WOvje6EkhBSzKyPcamqXv/NteufrlzP7VoxyJ9Wx\nM3Y7I6SY7XSjWR6Mpnutu2O2o0JrBBsIKWbznGTW5U11H4j3cMjSuF0QUswalBobMrDdvUpN\niH9lYA0hxaw2/JnHziZpeOzrAnsIKWbVSv05ZPIH3Id3eyOkmNWFH1U4XPJ9eOuKs7LC0EJI\nMWtR6sWQye8oNbe0Dwo7rZvDygpDCyHFbH344blnS75am5AGFEKKWVdCVec/8bh7vEr8a2kf\nREgDCiHFzdn2zMu72HuVUitL/JzU1RALnl+1vH1pWwRLawwNhBS3W61KzT4bmHR9hVL3l3r7\nUHqjM275QQanGwgIKXa9HTWquvVd7yEQvxxbc6dK/LHk2/DavQG9Ry491Ft8cfQvQorZn1a1\nL5qeCqB2avPM5rvHpl5PeWTenFkRw6VE6d4z27t3Yvy6M8WXR38ipJhpHCPQP0pwbtNk70fu\ne/3n/lxtFEFIMbMaUl9f8pi/i1e96COePV4xhBQzyyH1BXbxxq453U+rjSIIKWYPLli89Kn2\nIkr90HMb/V28Gdu4P7YSCGlwSH72lL+L9+gHPIg8doQ0aIhdvLpVX1V6bYYaQhpMxC7etK2X\nK702QwohDS7Jz9pHZVKqml/plRlKCGnQuXV4VYPBwT+UhZAGn6+3tdYQUswIqRKS3x/9cF8e\nKx99YlW992fSOCsfCC2EFL8PF9SUfxo21E9bp3kfVr3oINc5xIiQ4tadP+K9nZB6PlpQ7X3U\nva9x5V28CClukR2VF1LX8+O8z6lbzYVCsSOkmB3yu6nKZf6h13ff739q636ua6gAQorZUve3\n/eF3um7kPj/ZWPLYUyM4DVtphBSzqc6v++sWP+9iR4NX0ahnTlr8YJSEkGI2Uqlmax/22/vz\nvcvrEvPeZfCGCiKkmFUrtcHSR51eXedtjBo6LhZfHv2IkGIWMWRx6XY2+bt0y45b+USUgZBi\nNkup8h4X6/GOdbcf/M3K56EshBSzDqWesPJBbkVjn9jzzemTxz8/Fs7K90ALIcXsUo0aaeVm\n8MjzutYulUAJCCluW5V61MZVcIQ0oBBS7FYpNfe78j+GkAYUQopf5ziVmL/j87OXLgeV9ikT\ni+uXtUcoQorf5RfZgAw6hBS35Jbwu5EI6bZGSHFbz580gxEhxezrBCENRoQUs1VuMnVP7zhw\nhPOngwkhxWyG01HbjUqvBWwjpJiNUmoSt7AOPoQUsyqlVld6HWAfIcXM2SJtr/Q6wD5CilmT\nUtsqvQ6wj5BitlapZyq9DrCPkGJ2tkqNY3CFwYeQ4rZBqTWVXgdYR0hxSy5zSmKbNNgQUvw6\nJ6hJL53mZNKgQkgx27RuZXure5lQ9bi7GoMqvWooAyHFjBtaBydCihkhDU6EFDNCGpwIKWZP\nLnt27fMvhKr0qqEMhARYQEiABYQEWEBIMZvPkykHJUKKmVJ1a76u9ErAOkKKWepA9307r1V6\nPWAXIcUsc85oeNsRa09jxgBASDGb4J9/nbTxXKVXBtYQUsySn68c66WUmPsO91MMEoQUv94j\ny8Z4LY1a8WWlVwc2EFJF3Pp46UivpWmv/lTp1UHZCKlSbh5YPCKTUvXCgzae4YcKIqQK6n5v\n4fBMS+Oe76r02qAchFRZv7w123s8xcy9HHm4fRFSxV3YMj2TUt2mq5VeGRgipMq7tudu78BD\n7aZfK702MEJIFXb9ndZqeZts/YlKrxFMEFIl9Xy4yDvaoJo33Zc+hLe30msFA4RUOWfW+tc4\njFp52plwfJ77uuqjSq8YSkdIFXJtV7O/Ozfzbe8vo85RztuxHHK4/RBSJSQPt/m7dHVrvhVz\nvqlzJm2t2IrBFCHF74cX6/2N0Zx3fwvO/MiZeH9l1gtlIKSYde+Z7Vc0dv33+Qu0ONPjXy2U\niZBiVuvfQ/HQ/tAL7DY7s+JeKZSNkGKWyWjCi+cjFniDMVdvR4QUs9QB7taDvZEL7FBqXIzr\nAzsIKWZKTX7pYqEFvtm371BcKwNrCClmiw4z6MlgREiABYRUab3nz/5S6XVA2QipAi5uPeC9\nvLDcvSio8c+MY3ybI6TY9W4YrlZlXn81KnM4fPrliq4TykVIcUs+7nTTmn59cbx/lcO93Gd+\nWyOkuO1ws7kn/XqJuKPv5cquFspDSDG7Oc7dketMvf7WHfek5dOLp1a6tySxSbqdEVLMPnSa\nmZu54nuN21TqVqRNzqsDhX4MAxwhxczZ+CT+Nf0yOdHJ53Dq5c065R+AwO2IkGI2U6mZmZen\nnI4mZ65zeEKp2RVbJ5SPkGI2XqmNmZdbnZCey7x29u0mVmiNYAMhxWyEUm9kXs53Qvok8/p1\npUZVapVgASHFLKFU+pBdX+9Ipaq7M5P3KVVVsXVC+QgpZrVKZQauO6nE30V/UWpEpVYJFhBS\nzJy/kbakX7mHvLd4k7cwUsPtjZBi9oBSi1Mvko1OSGe8yUuVaq7YOqF8hBSzlc4+3BX3xQF3\noG9vas84pdoqt1IoGyHFzB23bt7Vvr7T7nDF672pzp9IalcF1wrlIqSY/Xan08yIWc3uZXbV\n5zITvxiVvd4BtyVCituu7AXfa9JTzi+vct48VtnVQnkIKW7JNq+jlsz13l+4b2rPFf4xDGyE\nFLuezakB9BNtNzITup3dvBqG4Lq9EVIFXO/s2LDzbPb9XarlVOXWBjYQ0gCw7wvGurvdEdKA\n8XOlVwBlIKSBIflZ27BKrwPKQEgDweUtDTyD4vZGSBWX/GRRdepweKVXBGUgpAq7uNl/Dmal\nVwVlIKRK6vmwNeE/wW9+pdcGZSCkyjm3ITvQ6qTNBR+ahIGOkCrkVuc8f2NUzUOTbnuEVBFd\n6+70N0aNW69UenVQNkKK3829s8SY38crvTqwgZDi9vWqUX5ENRyrGywIKVbdu//gVzT2+bNt\nhDRYEFKMvvxjrVdR1SMHevr6CGnQIKS4/LJzhr8xmtLxY2oaIQ0ahBSTp2r8P4yeOOod7Cak\nQYOQYuJldN/Oa9mJhDRoEFJMUhVN29wVmEhIgwYhxcTtaO7RnuBEQho0CCkm6R27O1cG7ion\npEGDkGIyzPsjafLG8/5EQho0CCkmV7ZO9W+YmNeZeRozIQ0ahBSb5NHF1V5LdWtSz6EgpEGD\nkOL0U2pshrSW3dcJafAgpHglP33c3yzVPn0/IQ0WhBS7y69MEXdRENLgQEgVkDy8yN8sqQWH\neiu9PigfIVXG5Zcn+ynVb7pQ6dVBuQipUpKfLPQ3S4nUTRW4jRFSBV3qmORvlsZv4AFJtzNC\nqqjkoQVV/map0iuDMhBSpf3oD7Va6TVBGQip8no/frSKkG5zhDQgXNxUT0i3NUIaIHoPtlZ6\nFVAGQgIsICTAAkICLCAkwAJCAiwgJMACQgIsICTAAkICLPj/q60YuD0Fp2EAAAAASUVORK5C\nYII=",
      "text/plain": [
       "plot without title"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 420,
       "width": 420
      }
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "ggplot(data = test, aes(x = Status_on_day_collection_summary, y = MK_standard)) +\n",
    "geom_bar(stat=\"identity\", fill = \"#8e8d8b\", width = 0.5, position=position_dodge())+\n",
    "  scale_fill_brewer(palette=\"Paired\")+ theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),\n",
    "panel.background = element_blank(), axis.line = element_line(colour = \"black\"))+\n",
    " theme(axis.title.x=element_blank(), axis.title.y=element_blank(),\n",
    "      axis.text.x = element_text(angle = 90, size = 30, vjust = 0.5, hjust=1),\n",
    "      axis.text.y = element_text(size = 30))+\n",
    "theme(aspect.ratio = 2/1)\n",
    "\n",
    "ggsave(\"MK_enrichment_CD38pos.pdf\", width = 10, height = 10)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 63,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0gAAANICAMAAADKOT/pAAACbVBMVEUAAAABAQECAgIDAwMF\nBQUICAgJCQkLCwsMDAwPDw8QEBAREREVFRUXFxcZGRkiIiIpKSkzMzM6Ojo7OztAQEBBQUFD\nQ0NHR0dISEhJSUlLS0tNTU1OTk5PT09QUFBRUVFSUlJTU1NUVFRVVVVWVlZXV1dYWFhZWVla\nWlpbW1tcXFxdXV1eXl5fX19gYGBhYWFiYmJjY2NkZGRlZWVmZmZnZ2doaGhpaWlqampra2ts\nbGxtbW1ubm5vb29wcHBxcXFycnJzc3N0dHR1dXV2dnZ3d3d4eHh5eXl6enp7e3t8fHx9fX1+\nfn5/f3+AgICBgYGCgoKDg4OEhISFhYWGhoaHh4eIiIiJiYmKioqLi4uMjIyNjY2Ojo6Pj4+Q\nkJCRkZGSkpKTk5OUlJSVlZWWlpaXl5eYmJiZmZmampqbm5ucnJydnZ2enp6fn5+goKChoaGi\noqKjo6OkpKSlpaWmpqanp6eoqKipqamqqqqrq6usrKytra2urq6vr6+wsLCxsbGysrKzs7O0\ntLS1tbW2tra3t7e4uLi5ubm6urq7u7u8vLy9vb2+vr6/v7/AwMDBwcHCwsLDw8PExMTFxcXG\nxsbHx8fIyMjJycnKysrLy8vMzMzNzc3Ozs7Pz8/Q0NDR0dHS0tLT09PU1NTV1dXW1tbX19fY\n2NjZ2dna2trb29vc3NzdpUrd3d3e3t7f39/g4ODh4eHi4uLj4+Pk5OTl5eXm5ubn5+fo6Ojp\n6enq6urr6+vs7Ozt7e3u7u7v7+/w8PDx8fHy8vLz8/P09PT19fX29vb39/f4+Pj5+fn6+vr7\n+/v8/Pz9/f3+/v7///+Zzp6kAAAACXBIWXMAABJ0AAASdAHeZh94AAAgAElEQVR4nO2d+38U\nxd7nz+4+z+7z7Nl9nvXZffZWk0AIkXANIRi5hRi5RRTDTYhAhIMIkUs4colywk1AOBGQiFFE\nEAMPqLgoBjUHMMgtRCCZ/pu2qme6p3umZ6Zm8u2Znnw/7x+O3dXfqa7wmvepmrr+wQAADJk/\n5LsAAAwHIBIABEAkAAiASAAQAJEAIAAiAUAARAKAAIgEAAEQiZId/9yW7yKA/ACRKNnwh5Z8\nFwHkB4hECURiC0SiBCKxBSJRApHYApEogUhsgUiUQCS2QCRKIBJbIBIlEIktEIkSiMQWiEQJ\nRGILRKIEIrEFIlECkdgCkSiBSGyBSJRAJLZAJEogElsgEiUQiS0QiRKIxBaIRAlEYgtEogQi\nsQUiUQKR2AKRKIFIbIFIlEAktkAkSgpNpP/nQb7LVKBAJEogElsgEiUQiS0QiRKIxBaIRAlE\nYgtEogQisQUiUQKR2AKRKIFIbIFIlEAktkAkSiASWyASJRCJLRCJEojEFohECURiC0SiBCKx\nBSJRApHYApEogUhsgUiUQCS2QCRKIBJbIBIlEIktEIkSiMQWiEQJRGILRKIEIrEFIlECkdgC\nkSiBSGyBSJRAJLZAJEogElsgEiUQiS0QiRKIxBaIRAlEYgtEogQisQUiUQKR2AKRKIFIbIFI\nlEAktkAkSiASWyASJRCJLRCJEojEFohECURiC0SiBCKxBSJRApHYApEogUhsgUiUQCS2QCRK\nIBJbIBIlEIktEIkSiMQWiEQJRGILRKIEIrEFIlECkdgCkSiBSGyBSJRAJLZAJEogElsgEiVB\nEUlXEIhEBkSiBCKxBSJRApHYApEoGbYiwbh0QCRKIBJbIBIlEIktEIkSiMQWiEQJRGILRKIE\nIrEFIlECkdgCkSiBSGyBSJRkKxL19xQi5RyIRAlEYgtEogQisQUiUQKR2AKRKIFIbIFIlEAk\ntkAkSiASWyASJRCJLRCJEojEFohECURiC0SiBCKxBSJRApHYApEogUhsgUiUQCS2QCRKIBJb\nIBIlEIktEIkSiMQWiEQJRGILRKIEIrEFIlECkdgCkSiBSGyBSJRAJLZAJEogElsgEiUQiS0Q\niRKIxBaIRAlEYgtEogQisQUiUQKR2AKRKIFIbIFIlEAktkAkSiASWyASJRCJLRCJEojEFohE\nCURiC0SiBCKxBSJRApHYApEogUhsgUiUQCS2QCRKIBJbIBIlEIktEIkSiMQWiEQJRGILRKIE\nIrEFIlECkdgCkSiBSGyBSJRAJLZAJEogElsgEiUQiS0QiRKIxBaIRAlEYgtEogQisQUiUQKR\n2AKRsqD/wvG2vccvPU14AJHYApFi3F0iJOmiwqfqRwiT0mVdcc8gElsgks3Z8UJDpF/nCgdN\nj10PIRJbIFKU3zdE3UgddrtCuKh/4nwKkdgCkSJcrRY6Ij2tUSGTd16+2XNxW7m6Xut8DJHY\nApEUg+8WSyfKGtKK1KoiGh9FbnrNVt4Fx2OIxBaIJHlgGjHn5qF0Ij0skwHzB6zbe5PkbZ3j\nOURiC0SSXJM+hP48YKQV6aB8XtwTu/9EfeCb2D1EYgtEMkyRJqme7LQivSCfv+a4D0+WCZti\n9xCJLRDJUCItua/+m06kOyH5/GNnSrNMqIzdQiS2QCRJ9/uR/6YT6bR6fs+Zckal3LJvIRJb\nIJKDdCJtk4+nuFJ61SdO27cQiS0QyUE6kZbJx4vcSaNl0rv2HURiC0RykE6kWfLxenfSNNeY\nLERiC0RykE6kCfLxO+6kl2TSYvsOIrEFIjlIJ9JI+bjNnaRae3PsO79FylccREoLRHKQRqSw\nenzYnbZaJs207yASWyCSgzQiPVaPj7rTmmTS8+ZVf0tLS/V/asnqxRCp4IFIDtKI9LuHSGtl\nUrV51fuM5J9aXE+DLghEIgMiOUgj0oB6fMSdtkYmzTCvHre3t8/5+xbX06ALApHIgEgO0nU2\nqBlCB9xJK2RSrX0X/xsp6IJAJDIgkoN0IqmVfLvdSYtk0sv2HURiC0RykE6kavdcb4VaMbvK\nvoNIbIFIDtKJtFA+Xu5OGusao4VIbIFIDtKJtMHuWbC4rz5xIhYAkbgCkRykE+moWiDr2oDr\ngvpEt30LkdgCkRykE6lbPb/oTNkhE8aG7VuIxBaI5CDtUvNK+XyzM2GmTHgjdguR2AKRHKQV\nSa3sG+do211RHzgXu4dIbIFIDtKK1KOGZHfYt4Pz5O3UWMsOIvEFIjnwEunGSYm9K4Oa7B06\na91tUfEnHcEQiS0QyYGXSAdcuzL0jpF3Ra196vrma+pRvaNCgkh8gUiG8cuCKFOVGdbNl5GH\nbpGMz4vUfUn9msZa1cwTFb3OnCASWyBSZKPVRKJNtjiRjI6RzqApP7tygkhsgUgZimT8EDsg\nqXh9nzsniMQWiJSpSIbx3a55lSVl1Us+uBOfE0RiC0SiBCKxBSJRApHYApEogUhsgUiUQCS2\nQCRKIBJbIBIlEIktEIkSiMQWiEQJRGILRKIEIrEFIlECkdgCkSiBSGyBSJRAJLZAJEogElsg\nEiUQiS0QiRKIxBaIRAlEYgtEogQisQUiUQKR2AKRKIFIbIFIlEAktkAkSiASWyASJRCJLRCJ\nEojEFohECURiC0SiBCKxBSJRApHYApEogUhsgUiUQCS2QCRKIBJbIBIlEIktEIkSiMQWiEQJ\nRGILRKIEIrEFIlECkdgCkSiBSGyBSJRAJLZAJEogElsgEiUQiS0QiRKIxBaIRAlEYgtEogQi\nsQUiUQKR2AKRKIFIbIFIlEAktkAkSiASWyASJRCJLRCJEojEFohECURiC0SiBCKxBSJRApHY\nApEogUhsgUiUQCS2QCRKIBJbIBIlEIktEIkSiMQWiEQJRGILRKIEIrEFIlECkdgCkSiBSGyB\nSJRAJLZAJEogElsgEiUQiS0QiRKIxBaIRAlEYgtEogQisQUiUQKR2AKRKIFIbIFIlEAktkAk\nSiASWyASJRCJLRCJEojEFohECURiC0SiBCKxBSJRApHYApEogUhsgUiUQCS2QCRKIBJbIBIl\nEIktEIkSiMQWiEQJRGILRKIEIrEFIlECkdgCkSiBSGyBSJRAJLZAJEogElsgEiUQiS0QiRKI\nxBaIRAlEYgtEogQisQUiUQKR2AKRKIFIbIFIlEAktkAkSiASWyASJRCJLRCJEojEFohECURi\nC0SiBCKxBSJRApHYApEogUhsgUiUQCS2QCRKIBJbIBIlEIktEIkSiMQWiEQJRGILRKIEIrEF\nIlECkdgCkSiBSGyBSJRAJLZAJEogElsgEiUQiS0QiRKIxBaIRAlEYgtEogQisQUiUQKR2AKR\nKIFIbIFIlEAktkAkSiASWyASJRCJLRCJEojEFohECURiC0SiBCKxBSJRApHYApEogUhsgUiU\nQCS2QCRKIBJbIJJN/4XjbXuPX3o6hCwgElsgUoTwqfoRwqR0WVfysMkigQbHY4jEFohk8utc\npxtNj5PFlUEkiOQJRFLcrnDLUf/EOy4cgkgQyROIJHlao4yYvPPyzZ6L28rV9VrvwAfqWd0C\nF9sdzyESWyCSpFX50fgoctNrtvIueAb+TT26nTwjiMQWiGQYD9Uvn/kD1u29Sara8Yy8pkTq\nT54TRGILRDKMg9KO4p7Y/SdKl2+8Irvkg6IUOUEktkAkw3hB6vGa4z6sOrk3eUV2ygflKXKC\nSGyBSMYd1RX3sTOlWSZUeoW2ywdTUmQFkdgCkYzTqiV3z5lyRqXc8gjdL9NrUmQFkdgCkYxt\nCdVMrxLptEfoTjXGlCIriMQWiGQsk3YscieNlknveoSqNt+yFFlBJLZAJGOWtGO9O2lakjHZ\nJjV/yDAGruzfvqFl37mEmUQQiS0QyZgg7XjHnfSSTFrsEbpEpjff2fpsdHZQSeN193OIxBaI\nZIyUSrS5k1Rrb45H6HyZ/oJz4mpoy6DzOURiC0QKKyEOu9NWy6SZHrEzE+esLhpwPIdIbIFI\nj5UOR91p6rfQ8x6xkeVI5c0Xbj55cK1tinn3VvTZg3nz5o36zy2u+KALApHIgEi/e4i0ViZV\ne8SarbqND6N3AztMky5F7nqfkfxTiys+6IJAJDIg0oCS4Yg7bY1MmuER+15ra+tZx705bTw6\nsDT4008/rfwPLa74oAsCkciASIaaIXTAnbRCJtVqfDQ8W5n0q32P30hsgUiGWsm32520SCa9\nrPPZk0qkdvsWIrEFIhnViXO91YrZVTqfva+qs832LURiC0QyFkoXlruTxiaO0SZhjIxcbd9B\nJLZAJGNDQs/CfdVgO6H14ckuCyESWyCScVQtkHVNm7ugROrW+vBo10Q9iMQWiGR0K20uOlPU\n+NDYsM5n7wpXIxAisQUiGUalq8PAiMwEesMj8NyRP68+6ErpUCKdsW8hElsgUmRl3zhH2+6K\nsuOcR2Ctqqlce4PXq91Q+uxbiMQWiGQYPaoPe4d9OzhP3k71atntVoY5q6TPhLvHDyKxBSIZ\nkcneIXvqzxZlx0nr7sZJSXT/ht/UiouSL+3PXTYXJn0XywgisQUiSXrVaFBRq9lEu/maOX/O\nrpAOOPdvMOfWFW37zby5u908v6LZkRFEYgtEUnxeZC54rV/TWGtuk1/Raz9yiTT4ijnfu6hm\n5boVsyMb6tdjPRIwIFKUjpHOxXpTfo49cYlkPFkRt65vrevcCojEFogU4YfYAUnF6/scD9wi\nGcbH1Q6Nas+4c4FIbIFIFt/tmldZUla95IM7ruR4kYxwV0v95LKi8ucWtl6JzwMisQUiUQKR\n2AKRKIFIbIFIlEAktkAkSiASWyASJRCJLRCJEojEFohECURiC0SiBCKxBSJRApHYApEogUhs\ngUiUQCS2QCRKIBJbIBIlEIktEIkSiMQWiEQJRGILRKIEIrEFIlECkdgCkSiBSGyBSJRAJLZA\nJEogElsgEiUQiS0QiRKIxBaIRAlEYgtEogQisQUiUQKR2AKRKIFIbIFIlEAktkAkSiASWyAS\nJRCJLRCJEojEFohECURiC0SiBCKxBSJRApHYApEogUhsgUiUQCS2QCRKIBJbIBIlEIktEIkS\niMQWiEQJRGILRKIEIrEFIlECkdgCkSiBSGyBSJRAJLZAJEogElsgEiUQiS0QiRKIxBaIRAlE\nYgtEogQisQUiUQKR2AKRKIFIbIFIlEAktkAkSiASWyASJRCJLRCJEojEFohECURiC0SiBCKx\nBSJRApHYApEogUhsgUiUQCS2QCRKIBJbIBIlEIktEIkSiMQWiEQJRGILRKIEIrEFIlECkdgC\nkSiBSGyBSJRAJLZAJEogElsgEiUQiS0QiRKIxBaIRAlEYgtEogQisQUiUQKR2AKRKIFIbIFI\nlEAktkAkSiASWyASJRCJLRCJEojEFohECURiC0SiBCKxBSJRApHYApEogUhsgUiUQCS2QCRK\nIBJbIBIlEIktEIkSiMQWiEQJRGILRKIEIrEFIlECkdgCkSiBSGyBSJRAJLZAJEogElsgEiUQ\niS0QiRKIxBaIRAlEYgtEogQisQUiUQKR2AKRKIFIbIFIlEAktkAkSiASWyASJRCJLRCJEojE\nFohECURiC0SiBCKxBSJRApHYApEogUhsgUiUQCS2QCRKIBJbIBIlEIktEIkSiMQWiEQJRGIL\nRKIEIrEFIlECkdgCkSiBSGyBSJRAJLZAJEogElsgEiUQiS0QiRKIxBaIRAlEYgtEogQisQUi\nUQKR2AKRKIFIbIFIlEAktkAkSiASWyCSTf+F4217j196OoRIiMQWiBQhfKp+hDApXdaVdSRE\nYgtEMvl1rnDQ9DjLSIjEFoikuF0hXNQ/yS4SIrEFIkme1ignJu+8fLPn4rZydb02u0iIxBaI\nJGlVRjQ+itz0mm23C1lFQiS2QCTDeFgmfZg/YN3emyRv67KKhEhsgUiGcVDqUNwTu/9EVTTf\nZBMJkdgCkQzjBWnDa4778GSZsCmbSIjEFohk3AlJGz52pjTLhMpsIiESWyCScVq1z+45U86o\nlFtZREIktkAkY5t0YYorpVfpcTqLSIjEFohkLJMuLHInjZZJ72YRCZHYApGMWdKF9e6kad5j\nsmkjIRJbIJIxQbrwjjvpJZm0OItIiMQWiGSMlC60uZNUG25OFpEQiS0QKaz6Cw6701bLpJlZ\nREIktkCkx0qPo+60Jpn0fBaREIktEOl3Dz3WyqTqDCN/+9d//dd/+a8trqdBFwQikQGRBpQe\nR9xpa2TSjAwj74ZCof/xX1p8KiUIOBDJUPN+DriTVsik2iwi45t2gA0QyVDr83a7kxbJpJez\niIRIbIFIRnXiXG+1DnZVFpEQiS0QyVgoXVjuThqbOPKqFQmR2AKRjA0JPQv3Va/CiSwiIRJb\nIJJxVC17dW2rdUHp0Z1FJERiC0QyupUMF50pO2TC2HAWkRCJLRDJMCqlDZudCTNlwhvZREIk\ntkCkyHq9cY4W2xVV8ZzLJhIisQUiGUaPGmjdYd8OzpO3Uz1adukjIRJbIJIRmcIdOmvdbVHV\nzEnr7sZJyS2dSAMiMQYiSXrHSCGKWvvU9c3XlB31djVzwLUrQ6pIAyIxBiIpPi9STpTUr2ms\nVY03UdFrP3KLlCrSgEiMgUgmHWrxq82Un2NP4kRKEWlAJMZApAg/xI49Kl7f53gQL1LySAMi\nMQYiWXy3a15lSVn1kg/uuJITREoaaUAkxkAkSiASWyASJRCJLRCJEojEFohECURiC0SiZMMf\n5rTzxHNqIicgEiXH/tsfPPm7P/6994M4/t0f/0Er7g//+I96cf/wx3+vFff3f/w7rbi/++N/\n9H7wP/P9T59vIBIpR+d5UvlMqfeDOGqf+V9acfP++7/oxf2fZ2ZrxY1+plIrbkqyP+TtfP/L\n5xuIlAtOP9OsFdf7jMcmYF783/+tFzfnGY/z0jzY8swprbhPn9mo92J2QKRcAJGGPRApF0Ck\nYQ9EygUPLt/Wihu4fEMvw6vf6sXduPxUK+7Xy/e14h5e1hOTHxAJAAIgEgAEQCQACIBIfrLo\ntN5PFFDwQCQ/EaK8+Xvt6Osfxa67dl+hLw7wD4jkJ+Y62tmHHujEfvmiqIjdHRNi2udDL0D/\n1eN7tm9M76RuHEgGRPKT6JL0kY3nPbfJcxDeJONGxe6Pqc+tH4zeHUjN/n17dm9LzPPRkYXm\nVi3iQ/P2/OZ73u9OHjdXA91/i2EORPKTifb2DpN39KSMXG9GxTZxNUUSa6N3QoP4HMPtE6xH\nEUE2iDFedVyquGzeyxSI5CfhixvGWV+4UP2Jx0kDOyIxsb29Pi01E6K7T2bxhR5sij2KCDJF\nluG0EU/KOIikDUTymcHza8dY37my9d94Bw1UyKfPvuestAY6p8q0ygHzJosv9KpocsgS5Gd1\nV5JQL6aMg0jaQCT/efr56tHW127a3t88Ik7LJ1PjJ988nCms/YtKHYwpjWZVZO+wN7G6ZsGK\nxBxFqOGD648tQQZPqCNw4w4cTBPXpMGQ/3mGBxApJzzuXDEq+q0vXnpmIP7xWpn+XcKHrorY\nr6QY7eUiNGffhV7ZThzsu/5xs6zMKhJabE/Ur7M6c+KeJYhhdMuqMXQ7qziQDoiUK/pPLbWq\nkPFb444DnCbEPI+P1MkaLD7tLfnFv+ZMCJ+UMmyJi1IVzQK7WRgVxDglr/dlFQfSAZFyyMNj\n80NRl+Ycd/Y8yMbUVo/4ZiHGxiXtE2LRk7i0W1VCHHAnySpuRLROcQgSlj+7VmQVB9IBkXLL\nrT0zoiqV74ytXCgSYo9HcKtsCLpT7paI0oT9XY0uIUa5x4hqhFgZvXQIYrwtRHVWcUnomzUr\nvi7kCkTKMQ/an7f6CEp3/h5NLBViu0esrJFK3Sl7hfD6dT8rvikmq7jD0UunIMfiqzjduCSE\nQ2K6ThwDIFIu6TvRUOzsOq6ITsqRbbNGj/CFQlS5U14R4qhH4GbZ4HMlyLd8Fr10CtIpW3JZ\nxSXhJyHG6MQxACLljIFPl9sd1rU7Z0e68I6bj5YK8Wz8Lx/DuCl/UC1xJ00S4oxH1rKiqnAl\nyBdZ5wo6BTkZX8Xpxnn/Qf9Wm9D2ZAtEyhHXN9tzHMo2qG63SwvNwSCzQjgkPNp2A/NFQifC\nSG+R9glR4kqYGKu5nIK0CVGZVZzi4dGVtZXjY4wxO07Kk/7FvIBIueDB4Vq7OTfnQ+uXUUeZ\nvB2nuhzujRAitHvQ9Zlb6qTnEXHDt1Kk9z3yl027ka4E+brV0UunIC8JUZdVnOR9e4KGi8XJ\n/2pWQCTfCZ9rtJt05c0/Op78oOYQtKmrzerh1P3XrUniD840jVBJO+LymizEKx6vqBZisith\nk6yiog46BLleJMTmrOIMY6OnRmLE12n+ei5AJJ/55Z0K+1u34KO4H0KfycQX1cWjaZGI4qq6\n+pfn15RHP/B8/DRX1dr7IuEd7TLVvZyhU9UVkRouJkhfTcKndePUryYvai/q/jsMdyCSn/S3\nz7e/c+O2/ZwYUCfTzYueKR5f0+re+Pi9MrW0Iy7xQ1XhtbqSzFmwy8wFhbYg3cqPyU+zizOX\nWkxb07xF/uetLc2vqxm11Z3Ym8sGIvmJNb9UhF4+nTDBTrFLPopc3V2a4NGKuwnxt0vUg9kH\nu63fU+Ebh2vMF8RtiKe6L0T5rm5LkP5zTWbHe7uRVdwnMmnSJXUlL8yx3xtvFonpPxogCkTy\nk6gRE9+5mSTgiIitQ7jwashh0cglXV4feE9YbcAXFry8oK5qRPR+U1xceEkkvUQN/1bXVkaz\nbhjMLu5Nqep31p8UnUTxdYUo9SwjSyCSn6jvZFHDmfhvb4yDQoyP3d3v2La0blrVjLnLd3b2\neX9gwFo/5GZ+f3xg31yPsLqEMM24WiFetv8kazZSzyRR6tFc5QlE8hMhKnenXI/ww8mTZ1M9\nTyS8qyThex9a47H2dmBTKD5sS+Kgr2ac/Im0y/qTbJGML4Woz6z0wxeI5CfLz6Xb9CQLbsWW\nr0daZcuTbAT+46oiR9iotYkrnrTjSuyBWxkRG9taLAS6vyNApAJk8NL+Py2omTpxyow5r7ee\n+T154P1TmxtqnquqqV/x7vkkbUW9uGK7+0FWX7F+jdNCsD9hLApEAhqMsaeXlwpx2U6+FR0G\nAxAJ6DBFiOjGeROsfSQUAyG95RYMgEi+092x589vN8ejniyUdFsXScjunXfu3PEatzrU0nIp\nq/wWCTEnclXnbM0NYPa3BUTyl4GDXlMWRGT0SP33a+siCdm9Vnj3AqyMXxe4ZPr06Tqbk2+X\nGX5lXq0WYpzdR/iN1nILFkAkX/m9PpUgORdpixAzXAmqKz1x6XoiV2RcWbsS6Ki8Wm8lL/Va\nbsETiOQrr6cURFukjPfgTiLSJqmDK0GJpLXFvzmptvi+YdxWo04rzZkad9UmrfEb5XEFIvnJ\nt1Ebyiuq4lGPGyW/WBdJMPNJIZp31eUt0sAMIYpcKdUyUGvC3OfqFUVqVKxRXYVmLV45x5ye\nFD95jysQyU/UXGnR2JV8z29NMhDpqom8//BqHFc6Fon43zSqgF77FyXSIiMnqItfypxvneI5\nF5chEMlP1MYM7xHkk4FI6cJqXBnfrRBijN58uUOjxSzz4nPH/i2l19J8ig0QyU8mCFGVfMaq\nPhnswZ1OpLj6p6dWVjTntArRd2xv5OLbqVZmtd2pP8IIiOQn8lfERp04wmGf6fFTUN28GN/M\nDHcsCIm6/ddSTDRKYPCLdS9Nm7Zg81eZFm4YA5H8ZGR0S4Z0CM1hHx0eXXhvyfgkGlW1xs3r\nbm3ZtO71F81nYyoSe0OANhDJT+QvkIM6cUlEShj20UZmeOJaPDfue8UlI7sXswUi+cmr3lsR\nJ5BEpIRhH22SZOgVB5FIgEh+ctieopYa7+994rCPNv6IdP2j2HXXbpyB7gQi+Un/WGurg2Rk\nNuyjTVdXV6r1RzYr16zfuGVbixcJsV++6NwX+ZgQ07zOduYKRPKVzpCoeZQqIEWVYBIZ9qGa\nJj6EY1jCsp0pRsXuzWPX11P07Q8PIJK/HAqJur+leJ5OpD12FMXs1iEcw7LefEms89wUyeNk\nTq5AJJ85N0mM2nAp6SQhvWEfKpGyP4alI/KS2I6Vn0b27DuZXXbDD4jkA/ZYTPX02npzIUXR\nWK9JqwqtYR8akYZwDIu5I+uz7/U4kzrVBIdKzLWLAJF8IGUl4/W9F2mGfbSniTuhPIZFHdo8\nNX6D4oczhXPdOW8gkg9kIxL5tlYZHcPSf/X4nu0bk3dor5WfTOx+VL2N+JUUASL5QBBEyuAY\nlkdHFkY2totson9+872EkGlCzPN4SZ0Q04jLXahAJB/wHJdJNUyjO+yjj/4xLOH2CdbTiEgb\nxJiEAaJyIbZ6vKVZ89BmBkCkAJP9sI/+MSyDTTHNIiJNESIU/8unyHv9Xyt2EbKASAEm+2Ef\n/WNYrE35Q5ZIP6u7kh53VKn3pMFm7CJkAZECTPbDPtrHsJw2JWr44PpjS6TBE+q8wLg9TaqE\nSOwWNIyFQmC5RQSI5CdDWrA3hGEf7WNYnkyUz+tuRAOjR192j5FuuQ/RWCrEs4lHWdyU1diS\n7Ao47IBIfpKkN85zwV7KYZ8t6XFlpnsMi6qQFgxYgdZhzKeEvdV3FHWwX0LbbkCd63kg1d/P\nCIjkJ0lE8lqwl3rYx7sPzoUrN91jWNYKMeK2HWiJFJ4qxApX3L0RspLa7Z6iemue6k7/zQAK\niOQnSUTyWLCXZtgnU5F0j2GpEWJlrKyWSMbbQlS7y7dZvWLq/uvWeU8PzjSZ29rt0PhXYAFE\n8hNvkTwW7KUb9slUJN1jWMqFOBx7hS3SsYTxoUfTIm8prqqrf3l+TXn0pc8Pecu+4QJE8ofM\nFuylHfZRj8uXbN24rmm11lw73WNYZM31WewVtkidsi6M+4N6vA4DqO41QASI5A/p6g/3Po1p\nh30iHxq/7oxmDaB7DMtIIawzbJ0infQYH7q7NOGPWJuwmMgAABqRSURBVHFXrzQcgEj+kE4k\n9zyBtMM+TaXRz41efVZnVaruMSwT7U4Jl0htnqdMXHjVuXhq5JKEQSnOQCR/yGyfxvTDPv3t\n860cJ2y5nvb1usewyBevNuw32yK9JESdV7b3O7YtrZtWNWPu8p2dxJMDCx2I5BMZ7dOoNezT\ns7PSymD2+4kTtN1oHsOySYiS3+w3WyJdLxJic6Z/MHMgkp8IzX0aNYd9wl12E694+Wcp16Zq\nHsPSKVMWD1pvjorUVyOvv8j4j+UNRPITobnOSHfYx9XEG9ec6igIvWNYzCXkyx5EyxoRqVt5\nNPmpTrmBDUTyE12RdId9THp22E28mfuS95vpHcOi5v6I8l3dlkj955qKEysukBaI5Ce6C/Z0\nh32ihL9aYzfxFn+SrO7QOoYlvCTypOR5+T/VtZXR6q4BG9ZlCEQKArrDPjEcTbzyjd+myT7V\nMSx9cz16Q+r6Iw+pNqZkAEQKArrDPi4cTbxpbTpHk3szsCm+pz60xepUFMKvY9eHHRApCOgO\n+8QR/qrJ6kkoWpT9239cVeQQY9Ta2HZBEEkbiOQDczVwf0Jz2CeRp+c2Vg39C33/1OaGmueq\naupXvHve+aMOImkDkXwgxdcuyfdPc9jHg+/3NZQkZrhk+vTp3w/9D8lqY0qeQCQfyFwkzWGf\neK5srLAix7seKLey/9kEMgci+UAWImkO+zj5rW2aFVe8/IzbOCXSA42SrjqHfm4aIJIPNGmQ\n8CGtYR+bgc+W2MLNOpAw865aJnvsvZWAmgJL0AQEECmApBr2idC91Z4PW77Jq9ZSCwS9dnSM\nx8xi5t60rcAhbYfEAohUcPQdfdGyqKjhtPe8hrvyt9OYnz0fubDyeeXk7+nitLdD4glEKizC\nXWtGWRqlGobtqZWttnNps7PzEqVNF8LJ45KI5LUdElMgUh7Iek/v261V1he/7M3U02HDHQtC\nom7/tdQ1TX/HspG2SxUtSX+TJRHJYzskrkCkPJDdnt5PPl5kTeYJLfwo9eYNrS2b1r0eaQGO\nqfA+KtCi73hDrJ9w9kHv6eTeInlsh8QWiJQHstnT+9omawssUdV6O110Bv3ukvvHFtqThIob\nOlyOZrYdEl8gUs5Jtqd3+OcLn55MIPLsUI3dpFur00+WmUiSu4ft2eSibN2lsE5OJjXe+bED\nIvmN5lGuny4pSfG9j96UN51J3Mrei4xFkvy6v84OqtzxU/qcTHQ62TkAkXxG7yjX/kWpv/fq\natzr7T9c+/rSxS5vXPmtXLN+45Zt6Y8KjOPmX2riXpzhdkh8gUj+onmUa1KP3DVSKogK/POe\nOSFXfhlth8QXiOQrmke5nrWfFMUTCcidSJJ7e5+Nz09obofEF4jkJ7pHua5WYa+c6H6UbEw0\ndyIN/tv2mlBCfoL+2PVhBkTyE92jXJVf76fKaFJ6KMr7uHPVGE8xIVI6IJKf6B7lOlq29nJb\nMg9+/6TR2ptIhOYedj3T3Q6JLxDJT3SPci32PjR86PRfPb5n+8YrGoGnV9oWiWnv3fSlNMMZ\niOQnuke5lgtxkP7tj45EJyxEdlA9vznZhuH9HStis1cnbU21mhAkASL5ie5RrvN8GNgMt0+w\n3IiItEGM+dwjrr9jeWwo+Nl1XSnmgIPkQCQ/0d3Tu1WI14lfPdgU69CLiDRF/vQ5HRf16KNl\nMYtGLu+MGxbCBpHaQCQ/0d3Tu7dEjCY+/W6V1W1gifSzuivpcUfF1lCE6v+a2J0gBLbj0gQi\n+Yn2nt5tQixOvWFQhpw25Wj44PpjS6TBE2r6eNw+eZYNNXt/9coFImkDkfxEf0/vjULU/2SQ\n8WSi/IrXmY1Ju2lndMsKMuRegWGqUPnnZAv6IJI2EMlPMtjTu2O8CC06ePFG7x032b1YVUgL\nIlVcTCTjlLCbmlGEGLvhcvLuBWwQqQ1E8pMM9vS+8w7l/+GvFWJEtO5xiBSeKsQKV1zj56QN\nSsZAJF/R3dM7vMd7NVK2ItXIl0UvHSIZbwtRnV2GIA0QyVd09/TeRvwTpFwIa46PU6Rjwvvk\nMjBkIJK/6O3p/X3y1XPZvbZYiM+il06ROmWLzzNeeyoRSAJE8hmtPb3N5X/lbxzsPJ964as2\nI4U4G710inTSc7MSjalERKdbDGMgkt/o7Ok9UyY2PiJ86UR7bpJLpDaPXg6tqUQ43SIdECl3\nJN/TWzb7Jic7UzkraoVYHb10ivSSEHXxhdKaSqR7ugVfIFIQkC2rTaQZbhKiJDpJ1iHSdfme\nzXGRelOJdE+34AtECgKyRtpPmmGn/OIvjhx9FBOpT20R9IU7UHMqke7pFnyBSEGgJn7GwVAZ\nUCf5LTMbY7ZI3cqjuBak7lQi3dMt+AKRgsBmId6kzfGQ2Q+4q9sSpP9cU7HH+JXuVCLd0y34\nApFyQbphmhtFYjztTovhJZGfPiXPy/+prq2MDlQ1xJ10qTuVSPd0C75AJN/RWfG9XYhm2rf2\nzfUY3a3rj4vSnUqUwekWTIFIPqO34jssq4Zm2jppYFP8dInQloR9UXWnEnlIaUFa6MIFIvmL\n3jCNpGOimLz7Gulg0o+rihxf+FFrv0sM0Z1KBJHSAZH8RW+YZueWDU0N6kHx+Oeq3Qzp7fdP\nbW6oea6qpn7Fu+c996XTnUoEkdIBkXwlwxXfOf+e6k4lyvJ0C0ZAJD/JaMV3PkTSn0oEUgOR\n/ER/xXeeRNKfSgRSA5H8RHeYZtXatzZv9Ww4tfhbQO2pRCANEMlPAr/iW3cqEUgHRPKTXK/4\nPpCa/fv27N7m/oTmVKIEnmDTFDcQyU8yXfE9VFL81kryo0trKtEvTTeMOP5SfcaXv6BggUh+\nktGKbwIyF0ljKtHAjiLrTA2b/nFCvBE/34g1EMlPMljxTUIWIqWdSmSaNi3uQ39RcTN/M4AF\nRPKTDIdpBq592PLm60tXb9zR0eP1PC2lDsZYB4cV2TvlT6yuWRA/q9tIM5XoUZ2ZGrc89gUz\nsZp44/9CBiL5SUbDNJfXjXZ8nyfuvJ0Ykgnt5SI0Z9+F3seGMdh3/ePmCiEqPOb4maSYSrTa\nLE3pWXfqk1azS6I+blUGYyCSn2QwTHOvMb4NVrxnKGd+vSVrPdemX+GTsqG5JdNsPjWLsipx\n5cc35qT294ZQwuEFRPIT/WGa2xUeP2deTlj2oM0+IRbFf/pWlRAHMssmPFOVo83rUY8qcSn2\n6IoCkXxFd5imb5olz4jycWPsrSRXe2aqwd0Sr+94l/wFlOwcWW/OqVIkmS10TZ07m3EVN1yB\nSL6iueLb2GG25V458l3kF8r9r/dHfuOfyvK9e4Vo8kielW6Tladx46zrZBmeSzbJoU0+HIeB\n2QgQyV/0VnzfV90Mi9w9dVdrlXtZ/pp/Jdbv7mSzfIvjdtBVPT0+9NJoERr/artjpa76IfRR\nspcMTJJPsSFKBIjkM1orvj+SyY3xXQuP62Xql9m9VX7FvWYeyIqqwr65v32Cc5+Ii5OsAlbZ\nL/1V3o1J/kNN1aM7sivgsAMi+Y7Gim/Zgno2cUfgX0cK0ZLdO0d6i7RPNjKt65NlQtTHHp0b\n4VD942jil6qvJPlbvlYdItkVcNgBkXJA2hXfNUKs80hutI+gzRQp0vseybJpNzJ6uUcZE5s6\n++uzzkqzOHoC+19Fyh7up8VCTM2ugMMOiBQEKoT4q0fyYSGmZJfhZCFe8UiuFmJy5OpspMFp\n9+xF9miZ2rj+FbNmmhHpQzjg0cHoZJwQ5dkVcNgBkYLAKCESt+gyjDNClGWX4XzhtTSvXabO\nNa8Gqsy+xHXWHJ/bqse9uF39TLtrjgwfM5Nb5VWy2RCKqYnHs3MFIgWBJD9pzme93GKvNKC0\nIy7xQzXprtW8VL0bYmGv/WircNQ9qnaKNCn3pKmRJjp+czEHIgWBiVYV4OZk9NjMzLltnu08\n+2C31X0evnG4xuxHiKwsWq7qJscQ0CS7rjIiaySE2Rf/vkh5CMVAMc6ktYBIQWCmEGs9kjcI\nMTPLHN+zug2qXljw8oK6KqtPLnIM0+BoIYp+ikV/J1xtuF1WRaTmCq5P/pLv5eOaLAs43IBI\nfvA4Pa74N4QYk9j93Vfu7ZcOA9bOlG7mR4aCe+RlgyN6t2oJxiYwXLOM+0HYvRNetInEzfa5\nApH8wPNL7MYVL9twYnX8gGx4tUgxrSAd4V0lCe8MrYn6e1beHHEEq74Jx3hRuFQqpy4G1Zqm\nuIOjHWFq3hPtuU6FC0Tyg0xFMqcINbqXyd1Rv2OefZh9IW5tGOd6Y8nyb61Hx+St45CZp8o5\n57jTbCGeMy9eTVXnHFW54qzzCBDJDzIVKTJp9dkNXdZwbf/5teYC19YhFWPw0v4/LaiZOnHK\njDmvt55xnGx0UGYd67Izroi4iqdBiPHmheowF5e8c783ViSuQWcLRPKDjEXqez6aPKlm3isv\n1VZE5+fV+nSsl+qKcEyxUDswlDlbliutbu3f1Ublk2565dGnJtW6GoisgUh+0ORGjd/EJcUv\ncrg5wUO2Kr92F1EiOWagq4MwXBMhVtgDra1mOTx2kLhjTmufnP3Sw2EGRMoB5Qk1UCK9ixI8\nWu05L48C1bSLrfsLPxvfiJRleTZy9dicAlF6KL4n5FS5WURsbmcBkXKAjkiGcfol54KL4tcu\n+lcg9dPnun2nOrmF622y1TYpevlNZLnuc/sca5f6jkQX9G7yr4iFBkTKAXoiyfbdybdfrq2e\nVF336rZOkp2uku3v1SULFJtApCYCjXQObKnub3ugtSOqd2ja2r3tnWc+3PfWLEv417A81gYi\n5QBdkYhJvr9Xr7zfaN+puUPOhbPGdZnQaN91lgpv1sOjGBApB+RFpJT7e00QYpzVT3BJPXNN\n9VNThHbHbrvrvDQal+1+EsMTiJQD0op0584dr/93P9TSkmQMJz2p9/f6U6x7YUBVSMXOluTj\nCXG/mcIfTY/PqXzXEIaKhyMQKQekFUk+/9ojeaX3XkA6pNnf65z60fOBunpiVlyuKX3blChx\nu65c2TIj1hVSubYTxyfFAZFyQLYibRFiRpavTLe/lzmaWrN939vmlidFvzg+qvrGHb+gbH6/\n2nn04IEjH1+5n2WZhjUQKQdkK9KmrFfIpt3f65JzQxanNtfMnfhGDnHjcX5ApByQpUgDM2Rd\nkd0b0+/vddDh0Yuxvu9tkZShzfHjCETKASlEumoin394NY4rHWquQ5YHkmns73Xc7tVe5phB\ncdhMeXUo2/fzBCLlgBQieXStuchyAarO/l4335laLMToxa69Ur9RL23EDLqMgUg5YAgipdgx\nIRWa+3uFf3sQV/c8LhYT21EfZQ5EygEpRJoev6GxmxcfJ/lcGrLf3+vPn2T5SuZApByQqrPh\n0YX3loxPolFVa7ZtLPL9vUAaIFIO0Om1O3EtnhtDGK8h398LpAEi5YBsx5GyZyb5/l4gNRAp\nB+ReJPr9vUBqIJIfzHWjprrFJc11xXd1ddGuhvVhfy+QEojkB+k6tdPVUEPGl/29QAogkh/k\nXSSf9vcCSYFIfpClSLfP7G9+Y8XSxnXbPvjq0dBKkOP9vQBE8oNJ6Un4zP2drtVzxQuODmnN\nT2739wIQKSAcGJ3wta8Y0mZXOd3fC0CkQBD+k1fzL5TlTLsoOdzfC0CkQPC+/X0fUTauLHbA\neOfQsiXf3wskAyIFgbujlDWl6078GGl8Pbx+cn2ZSipHY6xAgEhB4JCSZp1bmv7NqmF2ME8l\nAhkCkYLAQqnMloTUNhG3tz0ILhApCFQJMc5jY7tZQkzMfWFANkCkIFDqfeaxrJJG5rwsICsg\nUhAoFuIvHsmfYB1ewQCRgkC5d6/CuYzX4W1JD0mBQTwQKQjUCfGOR/IJIeozy8hrWDcOkgKD\neCBSENjm3T33VsaztSFSvoBIQaA7JIoTTzzunyBCf8ssI4iULyBSIJB1z8KEyd4bhdiQYT7m\nbIglWzeua1rdmASiEgM3ECkQPG0QYv4NV1LfeiFezHT5UKTSGb/uDDanyzEQKRgMtpaI4oaP\nrEMgHnY1jxWhP2W8DK/J2tB79Oqzg+nDARkQKQj8eWPT8hmmAKVTa+fUPj/OvJ7y6sIF85Js\nl5KM/vb51tqJCVuup48HRECkIKDRR6DfS9Czs9L6yOz37/lZbBADIgUBUpEMI9xlN/GKl3+G\ns8dzAUQKAsQiGa4m3rjmaz4VG8SASEHgpSUrVq9pSkOmmfbssJt4M/dhfazPQKRhTPirNXYT\nb/EnOIjcTyDS8MbRxCvf+G2+SzOMgUjDHkcTb1rbnXyXZrgCkRgQ/qqpLKpS0aJ8F2aYApF4\n8PTcxqosOv+ALhCJCd/vayiBSP4BkQJD+OcLn55MgCTrKxsrrJ9J40kyBPFApIDw6ZKSoQ/D\nevJb2zQrs+LlZzDPwR8gUiDoT9zxnkakgc+WFFtZzTqAmXe+AZECQVKPhiZS99bxVj7lmzBR\nyE8gUhA4a3tTFE/2mfYdfdHOteE05jX4C0QKAqvVt/2VE92P4s9Pzppw15pRGIbNIRApCEyV\nX/f3CfO73VplWVT25teEGYNkQKQgMFqIWrLMnny8yJpeF1r4ETZvyA0QKQgUC7GdKKtrm8qt\nyqiq9Xb6eEADRAoCSbYszpxDNXaTbu0lkhyBHhApCMwTYmjHxVpYfd1NZ56Q5Ad0gUhBoFWI\n10kyUhaNe739h2tfX7rY5Q3Je0A8ECkI9JaI0SSLwZOO65JNlQDeQKRA0CbEYopZcBApX0Ck\nYLBRiPqfhp4NRMoXECkgdIwXoUUHL97oveMms1wmpceX0gOIFBDuvIMKpJCBSIEgvMd7NRJE\nKhQgUiDYhp80BQ5ECgLfhyBSgQORgsBGpUz5Gwc7z2P8tECBSEFgpvSo8VG+SwGGAEQKAmVC\nTMYS1oIGIgWBIiE25bsMYEhApCAga6T9+S4DGBIQKQjUCLEv32UAQwIiBYHNQryZ7zKAIQGR\ngsCNIjEemysUNBApEGwXojnfZQBDASIFgvBaaRLqpAIGIgWEjoli8u5rGEwqVCBSENi5ZUNT\ng5omVDz+uWo3+S4a0AMiBQEsaC14IFIQgEgFD0QKAhCp4IFIQWDV2rc2b23xJN9FA3pAJAAI\ngEgAEACRACAAIgWBRTiZstCBSEFAiPLm7/NdCDAUIFIQMDu6Zx96kO9ygKyBSEEgOmY0svE8\n2WnMILdApCAw0R5/nbyjJ9+FAdkAkYJA+OKGcZZKofoTWE9ReECkgDB4fu0Yy6Wy9d/kuzgg\nQyBScHj6+erRlkvT9v6W7+KATIBIgeJx54pRUZWKl56hOMMP5AaIFDT6Ty0dGXVp/NbufJcG\naAKRAsjDY/Ot4ynmHEfPQ0EAkYLJrT0zoiqV77yf78KA9ECkgPKg/Xmr46F05+/5Lg1IB0QK\nIn0nGoqdy2QrruS7RCANEClwDHy63OptELU7Z0e68I7nu1QgNRApYFzfbM9xKNtwTSZcWqiu\niz7Ld8FASiBSkHhwuNZuzs350Ppl1FEmb8ehyyHQQKTAED7XaDfpypt/dDz5oVwmteWtYEAD\niBQQfnmnwq6MFnz0xP3wM5n4Yn7KBfSASEGgv32+bdG4bT8nBtTJ9NwXC+gDkYJAqb2G4uXT\nnhPsdslHuS4UyASIFASiGk1852aSgCPYczXgQKQgYHZwN5wZTBpwUIjxOSwPyBiIFASEqNx9\nO1XADydPns1VYUA2QKQgsPwcNj0pcCASAARApEAyePPGw3yXAWQCRAoKt9s6rctb69SkoOq/\nYB/jwgEiBYPB7SPFxuj1t2XR7vAZd/JaJpABECkQhF+T3jRErm9PsGc5zMI680IBIgWCg0qb\n6ZHrlY4Vfe/mt1hAG4gUBB6PVw25DvP6R7XvSd2Xt69uUEuSUCUVCBApCHwqnamPzvhuVk6Z\nS5F2yqvOVB8DwQEiBQFZ+YT+FrkMT5L6nDMvH5cLuwMCBByIFATmCDEnenlVelQZnefwuhDz\n81YmkBEQKQhMEGJH9LJNivR29Fq27SblqUQgQyBSEBglxJHo5SIp0hfR6/eFKMtXkUBmQKQg\nEBIi0mVnDI4Worg/mnxSiKK8lQlkBEQKAqVCRDeu+1o4fhf9VYhR+SoSyAyIFATkb6Q9kSvV\n5b3HSt6DnRoKBogUBOYKscK8CFdLka5byauFqM1bmUBGQKQgsEG24e6qi0610beVOjBeiMb8\nFQpkAkQKAmrfuoX3DeOa2q54m5UqfyKJw/ksFtAHIgWBJ2OlM6Pm1appdsU90cTLZbH5DiDo\nQKRAcDg24bs5knJzXZG8WZbfYgFtIFIgCDdaHtVF53tfVjelPak/BgIDRAoGA7vMDfRDjY+i\nCf2ymVeCLbgKBogUFPo6WrcfuhG7f07UXc1faUCGQKSgcvIy9rorICBSsLmX7wIAPSBSgAl/\n1Tgi32UAekCkwHJnTxXOoCgYIFIwCX+xvNjsDs93QYAeECmI3N5ln4OZ76IAPSBS4Bj4tCFk\nn+C3KN+lAXpApIDRsz220erkXSkPTQIBAiIFiacdC+3KqBiHJhUSECk4dG8Za1dG1W13810c\nkAkQKSA8Pj7Psef3pXwXB2QIRAoE328ssyUqQV9dAQKR8k//0Rdsi8ZtvdEIkQoQiJRvvvlT\nqWVR0audA4YBkQoRiJRXHh6aaVdGU1p/NdMgUiECkfLJmhL7h9HrF6zObohUiECkfGJpNPvQ\ng1giRCpEIFI+MS2atqvblQiRChGIlE+UR/UXBtyJEKkQgUj5JNKwG7vBtaocIhUiECmfjLB+\nJFXuuGknQqRCBCLlk7ttU+0FEws7oqcxQ6RCBCLll/CFFcWWS+XN5jkUEKkQgUh55zdzb4YI\ndUf7IFJBApECQPjL1+xqqfSNFyFSAQKRgsGd96Y4VlFApIIDIgWF8LnldrUklpwdzHd5QEZA\npABx591KW6WKnbfyXRyQARApUIS/WGpXSyFzUQUoDCBS0OhtnWxXSxO244CkAgEiBY/w2SVF\ndrWU78IAPSBSIPnV3mo13yUBekCkgDL4+eIiiFQ4QKTgcntnBUQqFCBSkBk805DvIgA9IBIA\nBEAkAAiASAAQAJEAIAAiAUAARAKAAIgEAAEQCQACIBIABPx/YyHs7wuh0oYAAAAASUVORK5C\nYII=",
      "text/plain": [
       "plot without title"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 420,
       "width": 420
      }
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "ggplot(data = test, aes(x = Status_on_day_collection_summary, y = myeloid_standard)) +\n",
    "geom_bar(stat=\"identity\", fill = \"#dda54a\", width = 0.5, position=position_dodge())+\n",
    "  scale_fill_brewer(palette=\"Paired\")+ theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),\n",
    "panel.background = element_blank(), axis.line = element_line(colour = \"black\"))+\n",
    " theme(axis.title.x=element_blank(), axis.title.y=element_blank(),\n",
    "      axis.text.x = element_text(angle = 90, size = 30, vjust = 0.5, hjust=1),\n",
    "      axis.text.y = element_text(size = 30))+\n",
    "theme(aspect.ratio = 2/1)\n",
    "\n",
    "ggsave(\"myeloid_enrichment_CD38pos.pdf\", width = 10, height = 10)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "R_4",
   "language": "R",
   "name": "r_4"
  },
  "language_info": {
   "codemirror_mode": "r",
   "file_extension": ".r",
   "mimetype": "text/x-r-source",
   "name": "R",
   "pygments_lexer": "r",
   "version": "3.6.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}