Skip to content

Instantly share code, notes, and snippets.

@levitsky
Last active May 31, 2016 12:52
Show Gist options
  • Save levitsky/519401cfc2edb7f25300 to your computer and use it in GitHub Desktop.
Save levitsky/519401cfc2edb7f25300 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Populating the interactive namespace from numpy and matplotlib\n"
]
}
],
"source": [
"import pandas as pd\n",
"%pylab --no-import-all inline"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"from pyteomics import pylab_aux as pa"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"/home/lab006/Data/External/Glio/Glioblastoma_20151002/results_DeMix\n"
]
}
],
"source": [
"%cd ~/Data/External/Glio/Glioblastoma_20151002/results_DeMix"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"r1 = pd.read_table('Glioblastoma_sc_S01_r01.mzML.demix_proteins.csv')\n",
"u = pd.read_table('Glioblastoma_sample1_union_proteins.csv')"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"i = pd.merge(r1, u, on='dbname') # пересечение по белкам"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.text.Text at 0x7f872c8e4550>"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAEPCAYAAACKplkeAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X2cjXX++PHXezDNFLlJhCShG5GoTS1l+pYZqUT92mq3\nDCm1bTebaRerMt1IKbZ0r9BIuVmiqY1DMqJVktvcZIhQiKTcDIbz/v3xuWacGWfGjDkzZ86Z9/Px\nOA/X+Zzr5nOda2bePveiqhhjjDGhEBPuDBhjjIkeFlSMMcaEjAUVY4wxIWNBxRhjTMhYUDHGGBMy\nFlSMMcaETKkHFRGJE5GvRGSJiHwrIqleemMvPVNExotIFS/9BBGZ4KV/KSKNAs7V30tfLSKJpZ13\nY4wxxVPqQUVV9wNXquqFwIVAJxFpCzwHDFXVZsCvQC/vkF7AL176v739EJHmwC1Ac6AT8JqIWEnL\nGGPKkTL5o6yq+7zNWKAKoMCVwCQvPQ3o6m138d4DTAau8rZvAMaparaqbgDWApeUbs6NMcYUR5kE\nFRGJEZElwDZgBrAO2KWqfm+XH4EG3nYDYBOAqh4CfhORU4D6wOaA024OOMYYY0w5UFYlFb9X/XU6\n0BY4N1SnDtF5jDHGhEDlsryYqv4mIrOBy4AaIhLjlVZO50gp5EfgDOAnEakMVFfVX0TkR6BhwOlO\n9/bNQ0Qs0BhjzHFQVSnpOcqi91dtEanhbccDHYFVwGzgZm+3ZOBDbzvdew/w/4BZAem3ikisiDQG\nmgELgl1TVaP2NXDgwLDnwe7N7s/uL/peoVIWJZV6QJqIVMIFsQmq+l8RWQmMF5GngUXASG//kcC7\nIpIJ/ALcCqCqK0VkIrASOATcp6H8JowxxpRYqQcVVV0OtAmSvh7XvpI//QDwpwLO9QzwTKjzaIwx\nJjRsnEeESUhICHcWSk003xvY/UW6aL+/UJFoq0ESEasVM8aYYhIRNAQN9WXa+yucREr8XRljisH+\nc1cxVZigAvZDbkxZsf/EVVzWpmKMMSZkLKgYY4wJGQsqxhgT4TZtgl9/DXcuHAsqxhgToVRhzBi4\n6CKYOzfcuXEqVEO9McZEi59/hnvugXXrYOZMaNUq3DlyrKRijDHllM/nIzHxJhITb8Ln8+Wmf/CB\nCyLnnQdff11+AgpUoMGP3sCeMOTIGCczM5Nvv/2WZcuWcf3119OmzVGzF0UN+30rOZ/PR7duyWRl\nPQdAfHxfxox5j6lTr2LBAkhLg8suC931QjX40UoqxgQxZ84csrKy2L9/P3NDVFn98ccf06BBA/r0\n6cMLL7wQknOWpdL4TkzBhg4d4QWUZCCZrKx3uf32VtSqBUuWhDaghJK1qVQw77//Plu2bGHBggV0\n69aNW2+9Neh+o0aN4qeffqJKlSqcc845dO3qVnv+6KOP2Lx5M/v376dRo0bceOONeY5btGgRM2bM\noF+/fqV+LwBTp05l5cqVxMTE0KBBA+64445i77dkyRLGjh2b5w99cnIyGzdu5NRTT+Wtt94KSV4f\nfvhhAFauXEnjxo1Dcs5givqdFPSMC0ovje/EFMdlnH/+CwwfnhrujBQu3HP4l8KaABpMQekVSWZm\npg4fPlxVVbdv3641atTQ77///qj9li1bpu3bt899f/XVV2tWVpZu3LhRn3/++dz0Xr166e7du3Pf\nHz58WDt37qxPPPFEKd7FEbt27dI2bdrkvr/00kt1+/btxdpv6NCh2q1bN+3Ro0eeY0aMGKEbN27U\nQ4cOFStP33zzjV5//fV6xRVX6Ntvv62vvvqq9u7dWzMyMlRV1e/366BBg3Tv3r3FOm9RFfU7KegZ\nF5SuWrzvxH7fSm769OkaG9tF4XeFNRoX11SnT59eatfznlmJ/wZb9VcFsmLFCoYMGQJA7dq1adq0\nKd98881R+02fPj3P/6Tr1KnDF198wY4dO/j00085ePAgACeddBKxsbG5+02ePJkrr7yyzOrSP//8\nc5o3b577vlWrVsyePbtY+/Xp04cbbrjhqGNiY2Np2LAhlSpVKlae2rRpQ7Vq1ejTpw+9evXivvvu\n45prruHBBx8EXEnvgQce4Mcfj1q0NCSK+p0U9IwLSofj/05M8e3fDzNnJlG16gQuvPAVOnbsx9Sp\nr5CUlHTUvgU15oeLVX9FoO+//77Q6odLL7006B/Kzp07M23aNMCVULds2ULTpk2P2q9atWpkZ2fn\nvs/KymL16tX87W9/w+/384c//IHevXuTmJiYG1R27NhBTEwMp556Knv37i3pLRbJ5s2bqVGjRu77\nGjVqkJmZWez9ggXBBQsW4Pf72blzJ82aNaNLly5Fztf8+fNzn8/BgwcZO3Ysffr0YcqUKTzzzDO8\n/PLLJCQkMGDAgALPcbzPuKjfSf5nvH//flatWlXgs7/qqqtK9J2Yolu4ELp3hxYt4Lvv4qhdu3+B\n++ZvzJ83L5kpU9KCBp+yYkGlHNi8eTMLFizgP//5D+PGjePgwYNce+21zJw5M+j+Z511FoMHDy72\ndapUqUKLFi0A+O9//8vFF1/MhRdeeNR+N954I6NGjUJV2bNnD2vWrOGSSy4BoF+/fgwePJhHHnmE\nF198MfeYDz74gLvvvpsxY8bkOVd6ejqVKlVi7ty5tGzZkunTpzNgwADOPffcYuc/v127dhEXF5f7\nPjY2lj179hR7v2CTH/bq1Su3d9aFF17IFVdckeePdUFWrVpF9erVmTt3LuvXr+frr79m2LBhnHHG\nGQB069atSPd2vM+4qN9J/mf83Xff8Yc//IE777yzwGd/vN+JKZrsbHj6aXjjDXjpJbjlFjjWvJx5\nG/MhK8ulWVCp4FatWsUll1yS+0d6/vz5NGrUqNSu99tvv/HOO+8wduzYoJ/XqVOH0aNH89Zbb1Gv\nXj1atmxJnTp1WLNmDRkZGcycOZNPP/2Unj170rJlSypXrkzbtm2P6ka6ceNGmjdvTtOmTXn88cfp\n168f1atXz/0Dm9+QIUPIysoK+llycjJnnnlmnrRq1arxyy+/5L7Pysqibt26Rx17rP2ClVQCg23N\nmjXJyMjIbbAuzGeffcYNN9yQ+0udnp7Oli1bCrznUCvqd1LQMy4oHY7/OzHHtmKFK53UrQuLF0P9\n+uHO0fGzoFIOdOzYkaeffprbb78dgFmzZhX6P43jrRoB9wd0yJAhvP3221StWpUffvghaABr3rw5\n559/PgBPPvkkTz75JOnp6dx8880AXH311aSlpTFv3jzi4+PZt28f06dP54svviArK4v09PTc6pFt\n27ZRrVo1atSowXXXXVdgvv/5z38W+FkwTZo0YeHChbnvd+zYEXTsx7H2y19SGTt2LJ988gnvv/8+\nAHv27KFy5aL9qmRkZHDXXXflvt+5cyfr16+nbdujVs4u1PE+46J+J3D0M37qqacKTC/Jd2IKdvgw\nDBsGQ4bA4MHQq9exSyeBUlJ6M29eMjn/F4uP70tKSlrpZLaoQtHaX55eRGjvr+uuu06//fZbVVVt\n37697ty5s1SuM3z4cF24cKFu2bJFv/rqq9xeSWvXrlW/36+qquvXr9cLLrhAVVVXrlypN910k6qq\nTp48Wd9///3cc33yySc6Z86cPOcfOHCgpqam5h67ZMkSHTVqlD722GOqqvrRRx+F7F727NmjLVq0\nyH1/wQUX6LZt2466n8L2U1UdPXp0nt5fc+fOzf1e9uzZo40aNcrtrZWcnHxUT7Ecfr9f69Spo7/+\n+mtu2oknnqg+n09/+umnkt5ukRT1OynoGReUXth3Ekx5/30rD9auVW3fXrVDB9X164//PNOnT9eO\nHW/Ujh1vLFHvMELU+8tG1JcTo0aNYvv27cTFxTF+/Hjmz58f8mvMmzePDh065H4PIsLGjRtp0KAB\nbdq0YeTIkbRu3Zrs7Gyefvpp6tatS2ZmJo8//jg1a9YE4KWXXmLv3r2cdNJJ1KhRg+Tk5NzzT5w4\nkWeffRYRoV+/fmzZsoXdu3dTr149Vq9ezWWXXUaDBg1y6+hD4d133+WHH37A7/fTpEkT/vKXvwDk\nuZ/C9nvllVeYOHEimzZtokePHjz88MOcfPLJvPfee2zfvp3169fz5z//ObekcfXVV3PbbbfRq1ev\nPPlYtmwZ48aN49VXX+XFF1/kzjvvBOCBBx7g1FNPpVGjRnm+q9JUlO+koGdc2LMv6DsJprz/voWT\nqms3efxxGDAAHnwQYspBP9xQjai3oHIMP//8M/fc04elS1fQosW5vPnmMOrVqxeKrOaaNWsWn376\nKYMHD+aJJ57g7LPP5rbbbgvpNUzJHTx4kNatW7Ns2TLrVnsMFlSC27wZ7rwTfvvNTbMSgv4qIWPT\ntIRIRkYGl19+LRdffBWvvz4izy9CdnY2l1/eiY8/rsv69W8xbdqZtGvXkQMHDuQ5x4EDB/j6669Z\nunQpfr+/2HmoXbs2zZo1Y+zYsZxzzjkWUMqp2NhYVqxYYQHFFJsqvPsutGkDHTrAF1+Ur4ASUqGo\nQyvsBTQEZgMrgG+BB730VGAzsNh7XRNwTH8gE1gNJAakd/LSMoG+BVyvsPrCPL766is98cRTFcYo\n/FdPPPE8femlV3I/X758uVat2lTBr+7Hwq/VqrXQhQsX5u6zZcsWbdy4hVar1lJPOuksbd8+KXcE\nsjEVVUG/hxXRtm2q3bqptmihunhxuHNTMCJoRH028LCqng9cCvxNRM4DFBimqq291zQAEWkO3AI0\n94LIa+JUAl7x0poDt3nnOW6jRr3Hvn19gDuAzuzb9yYvvzw69/O4uDj8/n3AQS/lEIcP7+aEE07I\n3eevf32ETZuuZffupezd+x3ffBPHCy/8uyTZMsZEiSlT3LT055zjBjUGGRYWdUq9T6CqbgW2ett7\nRGQV0MD7OFj93Q3AOFXNBjaIyFqgrbfvWlXdACAi4719Vx1v3ipXrgTsD0jZ76U5TZo0oUOHy8jI\n6EJW1o3Ex6fTtm3LPNNgrFjxHYcO3e9lrzJZWdexdOm8482SMSYK7NoFDzwAX34JkyfDH/8Y7hyV\nnTJtUxGRM4HWwJde0v0islRERopIztDc+rhqsRybcUGoPrApSPpxu+++uzjppFcReQ4YyYkn9uKx\nx/4emF8+/HAcTz7Zidtu+5rU1CuZPn0yMQFdNS688HxiY9/HFbwOEB8/iYsvblGSbBljIlDOHFwX\nXfQkZ5+dRY0abor6ihRQoAwHP4pIVWAS8JBXYnkdeNL7+ClgKNCroOOLIzU1NXc7ISGBhISEoPs1\nb96c+fM/49lnh7N373J69XqV66+/Ps8+VapU4ZFHHi7wWq+99gIrVnTmhx+a4Pfvp0OHy+jT56FQ\n3IYxphzy+XwMHToCcIMPk5KS8Pl8dO36V/bvnwY0IDa2O9dddxcnnRS+6VKOJSMjg4yMjJCft0y6\nFItIFeBjYJqqvhjk8zOBj1S1pYj0A1DVZ73PpgMDcfVLqarayUvvD/hV9bl859Jg91SaXRwPHTpE\nZmYmsbGxnHXWWUHnkjKmIonWLsXBVmOcMiWNRx+dy8KFfYFq3p5pdOyYzowZk8OW1+IKVZfiUi+p\niPsLOxJYGRhQRKSeqm7x3nYDlnvb6cD7IjIMV73VDFgAVAKaeQHoJ1xjfrnoe1u5cmXOO69EfQbK\nXLCFqXIUttCW3++nZs2aeaoAExMTmTBhQrEXfTIm0hw9gWMM99yzny1b+gBfA/8XzuyVC2VR/dUO\nuB1YJiKLvbR/4XpvXYhrjFgP3AOgqitFZCKwEjgE3OcVPQ6JyP2ADxdgRqrqcTfSV2TDhg1j3rx5\nVK9e/ajP/H4/jz32WIGjpX/44Qdef/11/vjHPyIiTJ06lcTERJYvX87o0aNzl5nt2LEjnTp1IjMz\nM2h64Ey6xkSu68nKWsm77y6hR48/5ynBhH0OrjApi95f8wjeIWBaIcc8AzwTJH1aYceZounTpw+n\nnHJK0PrUnIW2CloT5YQTTqBr166ceOKJ/Prrr1SpUoXzzjuP559/PujiTosWLQqaftVVV4X8vowp\nbSkpvZk7txf797cCzqFKlX6kpXWjU6ckqldPC2hrCe+aJuFk04yWE3v27KF///6sXr2a7du307hx\nY6655hp69+4ddP+SzFQMwad7L8pCW/UD5uR+8803c9ddL+qiTznpFlRMJGrYMImGDVewc+c6mjd/\ngAEDbs4NHklJSRU2kASyoFIO+P1+7rzzToYOHUrDhg25+uqrGTduXKFVRMe7iFOOYJ0JClpoK5id\nO3eyY8eO3IGgRV30KSfdmEhy+DD8+9/w3HPwzDPVueuuNoi8He5slUsWVMqB8ePHk5CQQMOGDQG3\nBOzBgwdLtd0hf0llwYIFQRfaKsiECRPydE44nkWfjIkE69ZBz55unZMFCyCgNtcEYUGlHPjiiy/o\n2bMn4Cax3LNnDyeffHKhx5S0+it/SeWrr74qdKGt/GbPnk337t3zpBVn0SdjyjtVePNNeOwx+Ne/\n4KGHyscU9eWdBZVy4IILLuDw4cMAvPHGGzz22GPHPKak1V/5SyMPPPBA7nZqaioikhtQ1q1bd9T4\nm8zMTOLj43Pfb9iwgRtuuIGlS5eyatUqGjVqRNOmTQtMN6Y8+/FHtwrjzp3w+ecQYSMGwsribjlw\nzz338P333zNp0iQ6dOhAu3btSvV6r7zyCqNGjSIjI4MnnniC33//PfeziRMnkp6eTnp6OpMmTQLg\n5ptvZsmSJXnOUatWLRo0ODJLToMGDejatSuvvfYaI0aMyC1FFZRuTHmkCmPHQuvW0K4d/O9/FlCK\nyxbpMsaEXCT+vm3fDvfeC2vWwJgxLrBUJLZIlzHGhMjUqXDBBdCsmZuivqIFlFCyNhVjTIW1a5db\nI37+/Io3RX1psZKKMaZCmjnTlU5OPrliTlFfWqykYoypUPbsgX/+Ez7+GEaNgquvDneOoouVVIwx\nFca8eW5J3337YNkyCyilwUoqxpiot38/PP646y78+utQyLhgU0IVKqjY4lnGVDyLFkH37nDuubB0\nKZx6arhzFN0qTFCJtD7zxpiSyc6GwYPh1VfdZJC33ebm7zKlq8IEFWNMxbFyJSQnQ+3arqQSMPmD\nKWXWUG+MiRqHD8OwYdChA9x9N3zyiQWUsmYlFWNMVPj+e+jRw21/9RWcdVZYs1NhWUnFGBPRcqao\nb9sWunaFjAwLKOFkJRVjTMTKmaL+l19sivrywkoqxpiIowrvvWdT1JdHVlIxxkSU7dvhr3+F1ath\n+nRo0ybcOTKBrKRijIkYH34IrVpBkyZuinoLKOVPqQcVEWkoIrNFZIWIfCsiD3rptURkpoisEZEZ\nIlIj4JjhIpIpIktFpHVAerK3/xoR6R7sesaY6LNrlxt3kpIC//kPPPccxMWFO1cmmLIoqWQDD6vq\n+cClwN9E5DygHzBTVc8GZnnvEZHOQFNVbQb0Bl730msBjwOXeK+BgYHIGBOdPv3UTVFftaqbor6U\nV9s2JVTqbSqquhXY6m3vEZFVQAOgC9DB2y0NyMAFlhu896jqVyJSQ0ROA64EZqjqLgARmQl0AsaX\n9j0YY8re3r1uivqPPoK334bExHDnyBRFmbapiMiZQGvgK6Cuqm7zPtoG1PW26wObAg7bjAtC9b3t\n/OnGmCjzv/+5Ker37HFT1FtAiRxl1vtLRKoCk4GHVHV34IzBqqoiEjjjY4mmfUtNTc3dTkhIICEh\noSSnM8aUkf37YeBAGDPGTVHftWu4cxS9MjIyyMjICPl5pSxm7xWRKsDHwDRVfdFLWw0kqOpWEakH\nzFbVc0XkDSBDVccH7NcBV/2VoKr3eulvAp+p6oR811KbkdiYyLN4MdxxB5xzDrzxhk1RX9ZEBFUt\n8TzOZdH7S4CRwMqcgOJJB5K97WRgakB6d+/YS4FdXjXZDCDRa2OpCXQEfKWdf2NM6crOhqeegqQk\n6N8fJk2ygBLJyqL6qx1wO7BMRBZ7af2BZ4GJItIL2AD8CUBVPxGRziKyFtgL9PTSd4rIU8DX3jme\nyGm0N8ZEplWr3AJap5zipqg//fRw58iUVJlUf5Ulq/4yJvx8Ph9Dh44AICWlN0lJSXk+9/vhxRfd\nIlpPPw29e9sCWuEWquovm6bFGBNSPp+Pbt2Sycp6DoB585KZMiUtN7CsX++mqPf74csv3eh4Ez2s\npGKMCanExJuYObMLR5pM06hV6ynee+9VNm5MYsAA6NcP/v53qFQpnDk1gaykYoyJGDt3nse111am\nSZPfmDOnOs2bhztHprTYhJLGmJBKSelNfHxf3MQYacCXwCT8/lM444y7LaBEOQsqxpiQSkpKYsqU\nNGrUeBVoCDwPnAAsJSbmcHgzZ0qdVX8ZY0LuwIEkYmI6ULnySA4d+hHwEx/fl5SUtHBnzZQyCyrG\nmJD57TfXAP/55/Dhh3Hs3ds0oGtx2lFdi030sd5fxpiQmDUL7rwTrr0WhgxxU9WbyGG9v4wx5cLe\nvdC3r1uVceRIm1G4orOGemPMccuZov7332H5cgsoxkoqxpjjcOCAm6I+LQ1eew26dQt3jkx5YUHF\nGFMsixe7SSCbNYOlS6FOnXDnyJQnVv1ljDmKz+cjMfEmEhNvwudzK0wcOuQmf0xKcm0okydbQDFH\ns5KKMSaPYBNCDh8+iREj2lOzpk1RbwpnXYqNMXkcPSHkAqpUOYfhw6tzzz02RX20ipiVH40xka4R\nl1zSn3vvtYBijq1I1V8i0g44M2B/VdUxpZUpY0z49OnTm4yMGWRn7we+JS6uC489Njrc2TIR4pjV\nXyIyFjgLWALkzganqg+UbtaOj1V/GXP8fvoJ7r4b1qz5jVNPHUjVqpuCrtxook+oqr+KElRWAc0j\n5S+1BRVjik8Vxo9383b99a8wYABUqRLuXJmyVJbTtHwL1AN+KunFjDHlz44dcN99sGIFfPIJXHRR\nuHNkIllRGupPBVaKyAwR+ch7pZd2xowxpe+jj+CCC6BRI/jmGwsopuSKUlJJ9f7NqVOSgG1jTDni\n8/kCppovuC0kcIr6CRPg8svLMpcmmh2zpKKqGcBq4GSgGrBSVeeUcr6MMcWUM2hx5swuzJzZhW7d\nknNHwweaNcuVTuLi3DQrFlBMKB0zqIjIn4CvgJuBPwELROTmol5AREaJyDYRWR6Qlioim0Vksfe6\nJuCz/iKSKSKrRSQxIL2Tl5YpIn2Len1jKoqhQ0d4o+CTATciPqfUArBvHzzwAPToASNGwOuv25on\nJvSKUv31KPAHVf0ZQEROBWYB/yniNUYDLwOB41oUGKaqwwJ3FJHmwC1Ac6AB8KmINMMFv1eAq4Ef\nga9FJF1VVxUxD8ZUaPPnQ3IytG0Ly5ZBzZrhzpGJVkUJKgJsD3j/i5dWJKo6V0TOLOC8+d0AjFPV\nbGCDiKwF2nr7rlXVDQAiMt7b14KKMZ6UlN7Mm5dMVpZ7Hx/flwceGEP//jB6tJui/sYbw5tHE/2K\nElSmAz4ReR/3x/0WYFoIrn2/iHQHFgIpqroLqA98GbDPZlyJBWBTvvS2IciDMVEjKSmJKVPScqu8\nunb9gAED/kiTJq50YjMKm7JQlKDyT+BGoD2u2upNVZ1Swuu+DjzpbT8FDAV6lfCcuVJTU3O3ExIS\nSEhICNWpjSnXkpKSuOqqJJ57DlJTYehQuP12m7PLHC0jI4OMjIyQn7dMZin2qr8+UtWWhX0mIv0A\nVPVZ77PpwEBcCSlVVTt56f0Bv6o+F+R8NqLeVFirV7u2k+rVoXv3OYwZMxwovHuxMVAGsxSLyBfe\nv3tEZHe+1+8luaiI1At42w3I6RmWDtwqIrEi0hhoBizAVZE1E5EzRSQWVwVnAzCN8fj98NJLrntw\njx7Qp4+P3r1vOWb3YmNCrcDqL1Vt5/1bok6HIjIO6ADUFpFNuJJHgohciKtOWw/c411rpYhMBFYC\nh4D7vGLHIRG5H/ABlYCR1vPLGGfDBujZE7KzXS+vpk0hMTGwezFkZbkux1ZaMaWtqFPfVwLqBu6v\nqhuLcqyq3hYkeVQh+z8DPBMkfRqh6SBgTFRQhZEjoX9/+Oc/oU8fqFQp3LkyFd0xg4qIPIArXfxM\nwNT3wFHtI8aYYyvqVCqF2bLFTVG/ZQtkZMD55+f9PFj34pSUtBLm3JhjK8rU9+uAS1T1l7LJUslY\nQ70pz/Kv/x4f35cpU9KKFVjGj4eHHoJ774VHHy14ivpQBC9TcZTleiqzgURvQGK5Z0HFlGdHr/+e\nRseO6cyYMfmYx+7YAX/7GyxfDmPGwMUXl2pWTQVTlmvUrwdme3NypXivPiW9sDHmaD6fj8TEm0hM\nvClPb62PP4ZWreD0090U9RZQTHlVlIb6jd4r1nvZ1PfGHKf8bR2xsX9nx45WJCbeRIcObRg06OXc\nqrF585IZO3Ys//3v1cyeDePGwRVXhDHzxhRBmQx+LEtW/WXC7VhtGTmf79ixjRUr1nDw4PMAxMQ8\njN9/J/BCzp7ExbUmObkOzz8P1aqV4U2YCqes21TyU1X9v5JevDRYUDHhVJyG+GDtK/AGMN97v5fW\nrYexaNFjZZBzU9GV5Rr1/wjYjgNuwg1MNMbkk3dNk+IPOhQ5GdXfgB3ExXVm8ODhpZdZY0rBMYOK\nqi7MlzRPRL4upfwYE3V27NhGYuJNQN7qsLztKzFUrryLuLjJnHXW69St+yUpKcOtG7CJOEWp/qoV\n8DYGuBh4SVXPKc2MHS+r/jLhlL/6Kzb270CV3HaT/NVhPp+Phx6ayNq1j1C16nbeeOMwt956Zbiy\nbyqwsmxT2cCR3l6HgA3AE6o6r6QXLw0WVEy4BTbU79jxC4sX9yTYuJRDh+Duu9fwzju1cOvNfX9c\ngyGNCYUya1NR1TNLehFjKpKkpKTcoJBT7ZXfd9+5KeozM7OA2cDNwOU28aOJeEUZ/GhMhZZ/QGJB\nAxSDSUnpTXx8X1zPrjTi4vrRpMnTtG8P3btDmzZPAfvK4jaMKRNFmqXYmIoqfxvJnDl3ANkcPPgi\n4AYoFlZdFbjEb1bWqeze/S3Llp3C//4HzZpBkyZ388UXNvGjiR42+NGYQhx7LEnBc3fltK2oQosW\nTzB2bAss2geNAAAbCklEQVT+8Q9ISck7Rb1N/GjKgzJrUxGR9sASVd0jIncArXG9v34o6cWNiVZH\nSjjDgcuYNWsXr732Bffe2+6ofQPbYIyJdEVpU3kd2CsirYA+wDpgTKnmyphyIn+bSGzsP4iNXZ37\nPj6+Lx06tDmqjcVVd00A/gQ0RHUJH3wwLGz3YUxZKUqbyiFVVRHpCryqqm+LSK/Szpgx5UFgmwhA\nSsq7APTvP5gffthMzZp1ePLJF/K0saSlvc+yZQ8DLQLO5C/jnBsTHkUZp/I5MB3oCVwObMdVh5XL\nlR+tTcWUtvyN9/AIMBZIAj7lhBNa07nzbqZNS2D//ieA41uMy5iyVJaDH+sBfwYWqOpcETkDuFJV\ny2UXFQsqprQFb7yfAHwC7Oaii4awcOFT1gBvIkqZBZVIY0HFlLajg4oPOA/4kbi425k69TULICbi\nlHrvLxHZjZueJdhFVFVPLunFjYlEKSm9mTXrNvz+SkAboCWQSuXKH/Loow9aQDEVWmG9v2bhJiR6\nCmihqtUCXhZQTIWUU6UVF9ceaAf86P3bjkOHhjBo0MvHHGVvTDQrMKioaldcy+MOYISIzBGRv+Wb\ntfiYRGSUiGwTkeUBabVEZKaIrBGRGSJSI+Cz4SKSKSJLRaR1QHqyt/8aEelerLs0phgKmobF5/PR\ntetdzJw5kH37JgADgQeBVFxVmGu8z2lHMaYiKnSciqruUtVRwDXACOAJoEcxrzEa6JQvrR8wU1XP\nxpWI+gGISGegqao2A3rjxsjkTL//OHCJ9xoYGIiMCZWcnl0zZ3Zh5swudOuWnBtYUlN97N+/FLgA\niAeuonLlneHMrjHlTqHjVESkHXArcAUwD+imqnOLcwGvx9iZ+ZK7AB287TQgAxdYbvDeo6pfiUgN\nETkNuBKYoaq7vHzNxAWq8cXJizHHEmzlxueff5tFi5L45puBwHKgfe7+LVuex+rVfW3uLmM8hTXU\n/wD8iusreTdwGFARaQOgqotKcN26qrrN294G1PW26wObAvbbDDTw0jcHSTcmpHbs+CVfysnMn58K\nwMiRS7jnnlvyrD8/eLALIEe6DttYFFOxFVZSWe/9m+i98gvJ8nTeaP3APsAl7tKWmpqau52QkEBC\nQkJJT2kqjEO4wYwA5wKXU6vWGGbMOJ+YmA7UqZMWNIBYIDGRJiMjg4yMjJCft0zGqXjVXx/ljMIX\nkdVAgqpu9QZXzlbVc0XkDSBDVccH7NcBF8ASVPVeL/1N4DNVnRDkWjZOxRw3NwalDXA9EAtMpWPH\nr4POQmxMNAnVOJUCG+pF5A/eH/yc98kiku71zipWD7Ag0jkyciwZmBqQ3t273qXALq+abAaQ6LWx\n1AQ64kacGRMyqtCyZSpwL5ANLCAmZggdOrQJb8aMiSCF9f4aARwAEJErgGdxjei/e58ViYiMA/4H\nnCMim0Skp3eujiKyBvg/7z2q+gnwvYisBd4E7vPSd+LGy3wNLACeyGm0NyYUtm6FLl3gs89acv/9\nk4mJSQJex++/08aeGFMMBVZ/ichSVW3lbb8KbFfV1PyflTdW/WWKa+JEePBBuPtueOwxuO66o+f2\nKmghLmOiRalXfwGVRKSKt301MDvgM1uG2ESEwtaT/89/PuO00+Zx552b+de/vuSppyA2NkwZNSZK\nFFZSGQBcixtR3xC4SFX9ItIMeEdVj17CrhywkorJkX+K+sDp55966hsef7w+rtf8ImJiHqRVq+YM\nHvwYQIHHGROtymSWYhG5DDgNN/Bwr5d2NlC1hONUSo0FFZOjTZsEFi/uSWA11pVXzqRJk7GMHbuN\n/fsXc2SyB7f2fHz8eqZMyT/2xKatN9GvLGYpjgcuBZoCdURkpKoeUtU1Jb2oMaXN5/OxdOm3+VL3\nMWfOYNasmcnZZw9j2bJb831en6ysexk6dAQzZky2QGLMcSisbSQNOAjMBToDzYGHyiJTxpTU0KEj\n8Pt7AH2BSsAZQBf8/lX8+ONPbN/+JbGx33DwYM4ROevQbw1Lfo2JFoUFlfMCBiuOxHXnNSaCtMQN\ngWqI+/H9EbgNgIMHoXXrt4DRLF36rReAttrcXcaUUGG9vw7lbKjqoUL2M6ZMFdajK0e7dhcDPwOt\ngEygJ67gfUTt2nVZtCiDTz55j44d19OxY7o1yBtTQoX1/joM7AtIige8uVjL78qP1lAf3Qrr0ZVj\n+XJo1+57du/OAoYBu3Dzln4LvFTgccZUZKXeUK+qlUp6cmNCLdjU9EOHjiApKYnDh+Huu9cwdmxd\nYmPfxE2sPRI3o8/twF3AG8TEZDJgwMMWUIwpBTaI0USFzEzo0mUXa9Zsw+9fRnZ2c9yqjABfAC+Q\nE4j8/jTmzElnwIAwZdaYKGZBxUSMQYMGMX/+fNxioU5cXD+aNfuMP/4RTjllEn5/FY6MS4FatZ4C\nYKct0GhMmbCgYiLCoEGDePTRIcBw3OqLf+eMMy6lVq1vWbToFL74Au6/fxrffdcYuMk7qjEXXdSK\nlJTeXjuMS7UeXsaUnkLXqDemvBg2bDQuoCTjqrI+YtOmsdxyyynMmwdnn403Rf1buNWquwBv0aFD\nG5KSkpgyxU0KaT28jCldVlIx5ZbP58udKuXgwX35Pm1OtWp30Lr1g1xzjdvHLQWcE3icnLaTpKQk\nCyTGlAELKqZcyQkkO3b8wvLlSzh0yHUBrlRpNjAT+BNu3Mn/0a3bdXm6F8fEpIQr28YYjwUVU27k\nH4Pi1oo/DUji8OHLqVQJTjrpdipXXkqfPg8zZ86iPN2L/f7lxMQ8jN/vjra2E2PKngUVE1aBVVw7\ndmzLEyScz4AkYD/x8d347bfvcj+ZM+cm8mpJq1bNqV07HYCUFGs7MaasWVAxYZO/ZBIT8zCuZ1eO\ny4AE4HXgcaBKnuNTUnozb17eXl2DB1sgMSacLKiYMhNYKklJ6X3U6HhXbfUwbqHRdrhJIFcCY4CD\n1K1bL8/5cnp1HTmnBRRjwq3QRboikc39VT4Fm7Orfv06rFuXwpHqrrHAb7jG+BHAxbip6B8Ffic2\nNob09PEWOIwpBWWy8mMksqBSPiUm3sTMmV0IXIUR/o4bKjUMOAVoBmwH7gb6efs9BLQAHgO20rFj\nOjNmTC7LrBtTIYQqqNjgR1Mmvv/++yCpJwCKm5L+j5xwwjBgHfAikA68AdQD5uEa640x5Z0FFVOq\nfD4fbdq0Z926tbguwmne6xHgEtyYk4bAyzRvvoL4+L64Kq8uxMauIzZ2a+4xrotw7/DciDGmSMJa\n/SUiG4DfcXOUZ6vqJSJSC5gANAI2AH9S1V3e/sOBa3DrvPRQ1cVBzmnVX+XEkXaUxl5KO2A97v8y\nPYD2QHUgDZGHmTZtHECexvz87609xZjSERVtKiKyHrhIVXcGpA0BdqjqEBHpC9RU1X4i0hm4X1U7\ni0hb4CVVvTTIOS2olAM+n4+bb76b3buzgWxcEBkLvIYLJqtxVV8dcEGlD9OmvW9Bw5gwiaY2lfw3\n0QVX34H3b1dv+4acdFX9CqghInXLJIemSAYNGsQppzTl5JMb0bnzjeze/RTwLC6ovI0LKFfiGuY7\n4wqiaUBfVHvmlkiMMZEr3ONUFJghIgq8qapvAXVVdZv3+TYgJ3DUBzYFHLsZON3bx4RZ3qnpwS2Q\nNRY3X9eJwKlAbVwvruuBQ7iuwpfgAstWXNWYMSaShTuotFPVLSJyKjBTRFYHfqiq6gWcHPlLNUHr\nuVJTU3O3ExISSEhICE1uTR6Bgxm//HIh+WcIdl2G38It5fsSMBRoCsQRGzsGyObgwS7AVpuny5gy\nlpGRQUZGRsjPW27GqYjIQGAPbpBCgqpuFZF6wGxVPVdE3gAyVHW8t/9qoENAqSbnPNamUgZ8Ph/X\nXfcXDh0a6qX0wVVr5QSVCUAdXAnlXmAxrvrLTcfy5JMpXHzxxdYIb0w5EfEN9SJyIlBJVXeLyEnA\nDOAJ3Bwdv6jqcyLSD6gRpKH+UuBFa6gPn6ZNW7NuXRNgiZdyAq5GcjhwBnA+bnqVuoAf14V4LG68\nSZoNYjSmnAlVUAln9VddYIqI5OTjPVWdISILgYki0guvSzGAqn4iIp1FZC2wF+gZnmxXbDlVXm7c\nyffkbUOJwbWfnOW9vwa4I+DoEdggRmOiW9iCiqquBy4Mkr4TV1oJdsz9pZ0vUzCfz0eXLndw8GB3\noAaun4Rb7wQaAM1xI+L7AQtxQSXQTxwZxGjtJ8ZEo3LTphIqVv1VOty4k97s3n0yrpqrF/AFsAX4\nADgb1yDfDdeu4vPevwBAbOw/OP/8s6ldu661nxhTDkVD9ZeJEEdKKM/j1jsZDYwCnsZVb03BjTvZ\njQs0TkzMAVq1Gk3t2qeQkvKuBRJjKgALKqZAOe0n8+bN9wLKaUBfXPfg1kAt3NxdCuwH/gp0xLWd\n/ESrVhewaFFGWPJujAmP8jCi3pRDOfN2zZzZmKysA17qCNzI+NtxvbvmAIO8z87DlVKSgMnAvdSu\nbRMeGFPRWFAxR/H5fPz5z38jK0tx1VwK3A9cC1wVsOcBXON7X9xkkWuwGYWNqdisod7k4fP5uOaa\nm1CNx02lci5wE65aaxsuuAzw9n4QN3vOtcTEjKJx43qcfPKpXhuKNcYbE0msod6EhM/no3//wfzw\nw2ZgDzt35lR1HQRexnUTbgzMAv6Fm9/zUdxQoZbUqrWViy5aT0rKOAsixhgLKhVN4Hxd9etXY8yY\nj1A9GzcifjtuMOMbwD+AG72jJgH/9rZb4tpONhAfv5b330+zYGKMyWVBpQI5smhWe1xgOIAriRzG\njTfJmRDyAK6LcI69uPVPrsJNt3KQqlVPYtKkdyygGGPysKBSgQwdOsILKDNxAWQ5rlQSx5EpVsDN\njJOA6yIMru0kFtfbqzaxsVuZNGm0BRRjzFEsqFQ4n+MCymm4qenjcaPeY3GBZBHwLbDKS/8J10j/\nKU2anMFZZzUjJeVFCyjGmKCs91cF4fP56Nz5Bvz+yrh2kR9wXYVb40bI18Wtd/IMUBPX0+tMXA+w\nbTz9dAoDBgwIdmpjTBSIpuWETSnq0aMHMTGn0KlTZ/z+KrgG+Xa4JX7fxw1UHI9beaA2UAm3OuPp\nQFVgK8nJXSygGGOKxEoqUax+/fps2bIfyMIFixNwEz/WAF4Bmnh7puHaVlYCJwO/ExtbmapVa9Kn\nT08LKMZUADZOxRwlsLvwN99ksHPnIVxV1xLcox6GG6zYHLcWSpOAo7/DBR8/yck38c4775Rhzo0x\n0cJKKlEib3fhz4DfgWrep9WAwcBfcvYG/owLMuB6d2UDwtNP/8tKJsZUQFZSMbmOzNV1GJiKawup\nhlvzBFy1VpeAI7Z6/z6E+xE4QJMmZ/Pqq89bry5jTIlYSSXCDRo0iEcffRrXZrIf1yZyEtAW+BoY\nAlwJPAwkekc94u2bjciJTJtmU6wYU9FZScXQo0cP0tI+xHUNVlwwyZlOZTSQgRtv8grwCa4dBWAf\nrivxCgsoxpiQsqASgVx1V7I3+WM84Mc1vt+Lm2YF3Kj4+3FdhX8H7sZNzfId0BqRFTz11CMWUIwx\nIWXjVCLMoEGD6NTpWnbuzAJeBBrhRsjXx3UVzvERrnSyBzd311vACuAgrVtXZtq08dYgb4wJOWtT\niSAdO3bk008/w7Wb9ALWA0uBbsBduMGL3+GW+H3EOyobF1iUJk3Ot8Z4Y0xQoWpTsaASIVxA+ZQj\nkz6eAJyD69WVANTBjYzPGYuyETf78AFcV+GBVjIxxhSowjbUi0gnXL1PJeBtVX0uzFkqdT6fzwso\nJ3gpsbgG+XOAZsBc3OqMB3A1mlVxc3bto0mTllY6McaUmYgqqYhIJVz9ztXAj7g+s7ep6qqAfaKu\npCIiuGBRw/v3Qdz8XABTcN2GwfX0ygYqUbVqJSZNGm/BxBhTJBW1pHIJsFZVNwCIyHjgBtxf0yhW\nCdddeBhwFm7d+ExcKeV3XJw9gAso2UyfPt2CiTEmLCKt91cDYFPA+81eWpQ7mSOrMl4OzAe64yaC\nfBD4DVDi4mItoBhjwirSSipFqtdKTU3N3U5ISCAhIaGUshMuv+JKJ38Hfqd169YsWrQozHkyxkSS\njIwMMjIyQn7eSGtTuRRIVdVO3vv+gD+wsT5621RySivgSie/A1jJxBgTEhWyS7GIVMb9F/0q3Dq3\nC6gADfWQE1hqAgLspEmTJqxduzbMuTLGRIsK2VCvqodE5H7c3O2VgJGBASWaRWOgNMZEn4gqqRRF\ntJZUjDGmNNka9cYYY8odCyrGGGNCxoKKMcaYkLGgYowxJmQsqBhjjAkZCyrGGGNCxoKKMcaYkLGg\nYowxJmQsqBhjjAkZCyrGGGNCxoKKMcaYkLGgYowxJmQsqBhjjAkZCyrGGGNCxoKKMcaYkLGgYowx\nJmQsqBhjjAkZCyrGGGNCxoKKMcaYkLGgYowxJmQsqBhjjAkZCyrGGGNCJixBRURSRWSziCz2XtcE\nfNZfRDJFZLWIJAakd/LSMkWkbzjybYwxpnDhKqkoMExVW3uvaQAi0hy4BWgOdAJeE6cS8IqX1hy4\nTUTOC1PewyojIyPcWSg10XxvYPcX6aL9/kIlnNVfEiTtBmCcqmar6gZgLdAWuARYq6obVDUbGO/t\nW+FE8w92NN8b2P1Fumi/v1AJZ1C5X0SWishIEanhpdUHNgfssxlo4KVvCpJujDGmHCm1oCIiM0Vk\neZBXF+B14CzgQmALMLS08mGMMabsiKqGNwMiZwIfqWpLEekHoKrPep9NBwbiqspSVbWTl94f8Kvq\nc0HOF94bMsaYCKWqwZoliqVyKDJSXCJST1W3eG+7Acu97XTgfREZhqveagYsACoBzbwA9BOuMf+2\nYOcOxZdijDHm+IQlqADPiciFuF5g64F7AFR1pYhMBFYCh4D71BWlDonI/YAPF2BGquqq8GTdGGNM\nQcJe/WWMMSZ6ROyIehF5XkRWeT3IPhCR6gGfRd0AykjOew4RaSgis0VkhYh8KyIPeum1vI4da0Rk\nRkBvQERkuHfPS0WkdfhyXzQiUskb0PuR976xiHzl3cN4EanipZ8gIhO89C9FpFF4c35sIlJDRCZ5\nv3crRaRtlD27h72fy+Ui8r73jCL2+YnIKBHZJiLLA9KK/bxEJNnbf42IdD/mhVU1Il9ARyDG234W\neNbbbg4sAaoAZ+LGugiu2mytl1bF2+e8cN9HEe81YvOe7z5OAy70tqsC3wHnAUOAf3rpfQOeZWfg\nE2+7LfBluO+hCPfYB3gPSPfeTwT+5G2/Dtzrbd8HvOZt3wKMD3fei3BvacCd3nZloHq0PDtcG+73\nwAne+wlAciQ/P+ByoDWwPCCtWM8LqAWsA2p4r3VAjUKvG+4bD9GX1w0Y6233B/oGfDYduBS4DJge\nkN4P6BfuvBfx/iI278e4r6nA1cBqoK6Xdhqw2tt+E7glYP/c/crjCzgd+BS4EtejEWA7R/7zc2nO\nc/R+Ltt625WB7eHO/zHurTrwfZD0aHl2DYCNQE3veXwEJEb688P9RzQwqBTneZ2G6xD1ekD6G8Ct\nhV0zYqu/8rkT+MTbjsYBlA2I3LwH5fXkaw18hfsh3+Z9tA2o620He2anl1EWj8e/gX8AfgAROQXY\npap+7/MfOfLccp+pqh4CfhORWmWb3WJpDGwXkdEiskhE3hKRk4iSZ6eqP+LGy23E9TDdBXxD9Dy/\nHMV5Xjl/N4P9PS1QuQ4qhQygvD5gnwHAQVV9P4xZLW1R1ZtCRKoCk4GHVHV34Gfq/jsUeL/5u4iX\ny+9CRK4DflbVxRzJczR1b68MtMFV+bQB9uJKzLki9dkBiEhNoAvuf/b1cdWz1xR2TKQrwvM6LuHq\nUlwkqtqxsM9FpAeuLvCqgOQfgYYB70/HRVfJl96QvBG4PMt/T5GU9zy8hs7JwLuqOtVL3iYip6nq\nVhGpB/zspQd7lj+WXW6L5Y9AFxHpDMQBJwMvAjVEJMb7327OzyK4+zgD+ElEKgPVVXVnGPJdVJuB\nzar6tfd+Eq6qeWsUPDtw1bDrVfUXABH5AGhH9Dy/HMX5XdvspScEpDcEPivsAuW6pFIYEemEq2q4\nQVX3B3yUDtwqIrEi0pgjAygX4g2gFJFYXONaelnn+zhFct5ziYgAI4GVqvpiwEfpuEZRvH+nBqR3\n9469FFcVsY1ySFX/paoNVbUxcCvwmareDswGbvZ2SwY+9LYD7/n/AbPKMr/FpapbgU0icraXdDWw\nAtf2ENHPzvMDcKmIxHs/pzn3FxXPL0Bxf9dmAIlez7+auA5SvkKvEO6GpBI0QGXifhAWe6/XAj77\nF6631GogKSD9GlyPo7VA/3DfQzHvN2LzHnAP7XHtDUsCnlsnXA+TT4E13g9xjYBjXvHueSnQJtz3\nUMT77MCR3l+Nce1GmbgeRVW89BNwPYsygS+BM8Od7yLcVyvga+9ZfIBrvI+aZwekAqtwM3yk4Xpa\nRuzzA8bh2ocO4tpLeh7P8/KOy/Reyce6rg1+NMYYEzIRW/1ljDGm/LGgYowxJmQsqBhjjAkZCyrG\nGGNCxoKKMcaYkLGgYowxJmQsqBhzDCJy2JvOfrmITBSReC99gDdV+lLv8z946Rki8kO+c0wVkd3B\nzl/EPNwsbsmAwyLSpmR3ZEzpsaBizLHtU9XWqtoSN5DsXhG5DLgWaK2qrXBTBQVOnfOriLQDtw4J\nUI/jnPtKRGJwA/K6AZ8f/20YU/osqBhTPHOBprhpwXeoajaAqu5U1S3ePoobfX2r9/5G3HxnAiAi\n9UTk84DST/v8FxGRDSLyrIh8A/w/VV2tqmtK99aMKTkLKsYUkTdxYGdgGW6Ki4Yi8p2IvCoiV+Tb\nfRZwhVfKuAUXZHL8GbcuR2vgAty0NfkpLmhdpKoTQ30vxpQWCyrGHFu8iCzGzXu1ARipqnuBi4De\nuIWcJohIcsAxh4F5uEWO4lQ1sI1lAdBTRAYCF6jqngKuO6GAdGPKrXI99b0x5USWV6rIQ9106HOA\nOd464Mm4iQjBlTTGA1OAgfmOmysilwPXAe+IyDBVfTfIdfeG8B6MKRNWUjHmOIjI2SLSLCCpNa4U\nk0tV5wLP4GaLDTz2DNzys28Db3vHFuvyxc6wMWXESirGHFuwXltVgZe9nl2HcNOC9z7qQNVhQc5z\nJfCIiGQDu/HWsSjsmiLSDRgO1Ab+KyKLVTWqVyY0kcmmvjfGGBMyVv1ljDEmZCyoGGOMCRkLKsYY\nY0LGgooxxpiQsaBijDEmZCyoGGOMCRkLKsYYY0LGgooxxpiQ+f/rjeIAlTf2ZgAAAABJRU5ErkJg\ngg==\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f872c8b3828>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"pa.scatter_trend(i['PSMs_x'], i['PSMs_y'])\n",
"pylab.xlabel('PSMs r1')\n",
"pylab.ylabel('PSMs union')"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAEACAYAAABWLgY0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXl4VNX5xz9nCIEBAiHsiEaJIC6IQav0h5W0Gqh7gS4u\ntBEXam1BICimCCIkRVuhVmuLWAUqorUqFrsk4hIEqrUFRKyAsogiakFcggZDmPP747135s5kJpkk\nk23yfp5nHmbuPffcc2+Y9575nncx1loURVGU5MTX1ANQFEVRGg418oqiKEmMGnlFUZQkRo28oihK\nEqNGXlEUJYlRI68oipLE1MvIG2O+Z4z5rzHmiDFmaKIGpSiKoiSG+s7kNwOjgZcSMBZFURQlwaTU\n52Br7VYAY0xiRqMoiqIkFNXkFUVRkpgaZ/LGmFVA7yi7fm6tfSbxQ1IURVESRY1G3lqbW58TGGM0\nOY6iKEodsNbWWwtPpFwTczDW2mb1uu2225p8DC1lXDomHVNrGFdzHFOiqK8L5WhjzHvAMOBvxph/\nJGZYiqIoSiKor3fNCmBFgsaiKIqiJJhW6V2Tk5PT1EOISnMcl44pPnRM8dMcx9Ucx5QoTCK1n6gn\nMMY29DkURVGSDWMMNgELr/WSa+qDBlApSuOik63WSZMZedD/dIrSWOikqvXSKjV5RVGU1oIaeUVR\nlCRGjbyiKEoSo0ZeURQliVEjryiKksSokVcURUlimtSFUlGamrfffps33niD119/nYsvvpihQ7WK\npZJc6ExeadX89a9/5aijjmLq1KncddddTT2cmKxevZry8nIOHTrEmjVrmno4SgtCZ/LNhNdee41l\ny5ZVa2iefvpp3nzzTXw+H0cddRQ//OEP67S9oYn3vLUd90MPPcTevXtp27YtJ5xwAt/5znfqPdYp\nU6YA8Oabb3LcccfVu79YxHtPYl1jXl4e7777Lj169OCBBx4AIBAI0LVrV3y+0FwtNzeXxx9/vMGu\nQ2mBNEJOZBuNWNtbI/Pnz7ejR4+2V111Vcw2n376qR06dGjw87Bhw+z+/ftrtX3fvn0NcwE1jDPa\neWs77tdff92effbZwe3nnXeeLS8vj2tM69evtxdffLE955xz7B/+8Ad733332QkTJtjS0lJrrbWB\nQMAWFRXZL774otbXGw/x3pPqrnHRokX23XfftZWVlcH9O3futI888ojdtWuXfeedd+zdd99t33zz\nzahj0O9by8P5m9XbBqtc0wyYOnUql156abVtXnrpJU466aTg5yFDhvDCCy/UavuLL76Y+MHHMc5o\n563tuIuLi8Nm2j179mTdunVxjWno0KGkpaUxdepUrrnmGm644QbOP/98Jk2aBMAzzzzDxIkTef/9\n92t9vfEQ7z2p7hpTU1M5+uijadOmTXB/u3bt+M53vsOxxx5L586dadu2LSeeeGKDXIPSclG5JkHs\n3Lkz+DM6GsOGDavWkNsa8vjs2bOH9PT04Of09HTefvttMjIyarW9oYk1znjbVXc9hw8fDm4/dOgQ\nW7Zs4dxzz41rXC+//HLw71NRUcGyZcuYOnUqK1as4Be/+AX33nsvOTk5zJgxI2Yfdf0bx3tP0tLS\nwq6xvLycrVu3cu655/Lqq68SCAQ4cOAAAwYM4JJLLqFv377Btvfff39QelIUL2rkY7Bnzx5effVV\n/vznP/Poo49SUVHBhRdeyKpVq6K279+/P/Pmzavz+WpKIPXpp5/Svn374OfU1FQOHjyIMaZW2wFW\nrlxJmzZtWLNmDYMHD6a4uJgZM2YwaNCgOo+/pnHW53q++OILrr32Wh566CGstRw8eJBt27bxta99\nLa4xbdmyhS5durBmzRp27drFv//9bxYsWMAxxxwDwOjRo+Pqp65/43jvyZgxY8Ku8a233uLMM88E\n4Jprrgl6/px22mmcc845wQfHgQMH2L9/P+3atav12JTkR418DLZs2cKZZ57J3XffDchMMDMzs8HO\nV9NMPi0tjY8//jj4uby8nF69etG+fftabX/33Xc58cQTGTBgALNmzeKWW26hS5cuQYMXyS9/+UvK\ny8uj7svLy+PYY4+Na5z1uZ7evXvTs2dPFi9ezAMPPECfPn0YPHgwPXv2jHG3wnnhhRe49NJLGTVq\nFCAPuQ8++CDmNSeaeO9Jddd42mmnBdt17dqV0tLS4KLsn/70J5VplJiokY9Bbm4uhYWFjBs3DoDn\nn38+aCSiUV+5pqaZfFZWFv/5z3+Cnz/++GOGDh1Kenp6XNv379/P0KFDg4bto48+Ii0tjfT0dC66\n6KKY57355purHVdN43TPW9fr8R5/0kkncfLJJwMwZ84c5s6dG9eYSktLufbaa4OfDxw4wK5duzjr\nrLNqdW11/RvHe08g+jUuW7aMv//97yxfvhyAgwcPkpIS+uq++OKL/OhHP6rVtSitiESs3lb3ogV7\n11x00UX2jTfesNZae/bZZ9sDBw402LkWL15cxbtm+/btNhAIWGutPXjwoD3llFOC+0499VT70Ucf\n1Xr7li1b7GuvvWYfeughO3PmTGuttc8880zCriPWeet7Pbt27bKnnnqqtdbaN998044dOzbYJi8v\nL6ZnUiAQsD179rSffPJJcFuHDh1sSUmJ3bt3bwKuuGbivSexrnHNmjVBT6CDBw/azMzMME+g0047\nzT733HPVjqElfN+UcEiQd02Tlf9zSls16Lnry0MPPcS+ffto3749jz32GC+//HKDnOe3v/0tjz/+\nOO+99x5XXXUVU6ZMoXPnzgwdOpQHH3yQ7OxsAB5++GF2795NIBAgKyuLK6+8stbb77nnHsrKyujT\npw9bt27l61//OkcddVRQ+00EscZTn+s5fPgwhYWF9OrVi7fffptZs2bRtWtXAM477zwuv/xyrrnm\nmrBxvP766zz66KPcd9993H333Vx99dUATJw4kR49epCZmUleXl7Crru+96S6a3zkkUfYt28fu3bt\n4oorrgj7FXLuuedy3333Vbum0hK+b0o4iSr/1yKN/MMPP8Kdd/6OQCDAtGk/5uqrr6rnKKvy/PPP\n89xzzzFv3jxuv/12Bg4cyOWXX57w8yj1o6KiguzsbF5//fUw90IlHDXyLY+kNfJffvklN988izVr\nXqV//2O45555HH300cH9TzzxJHl5+Xz55SLAR4cOP+b+++cybtwVYf28++67vPfeewwcOJAePXrU\netybNm1i/fr1pKamkpKSwmWXXVbrPhSluaBGvuWRtEZ+5MjvsGZNWw4duoE2bUrp0WMZ27ZtpHPn\nzgCMGvVdnn12NHClc8QTfOMbS3nppWeCfcyffw8zZ84hNfV4Dh/ezuOPL+XCCy9siMtTlBaBGvmW\nR6KMfLOKeP3kk08oLX2eQ4ceAb7JkSO388UXx/LSSy8F2/j97YFPPEd9it8f8g/etm0bM2cWUV6+\nkc8+e4Uvv/wrP/jBjzh06FCjXYeiKEpzoVkZedFUA0CFs8UCh8K01hkzbqRDh9uBecCddOhQwKxZ\noUi/7du3k5qaDbgSzzCs9fPhhx82xiUoipJElJSUMHLkWEaOHEtJSUlTD6dONDu55oorruEvf9nF\nl19eQ2rqajIzX2XTppfx+/3BNhs3buT3v3+IQMBy/fVXccYZZwT3bd++nVNP/T/Ky/8JHA+8SFra\nD9i37z2NCFRaLSrX1J6SkhJGj86jvPxOAPz+6axYsbTaeJlEkrSafGVlJfPn/4bVq19l4MBMbrut\nIOhGFi+LFj3IjTdOIzX1KAKBj3j66cfiznGiKMmIGvnaM3LkWFatugRw3WyXkpu7kmeffbJRzt8s\njLwx5lfARYi+sgMYb639LKJNk/jJ79+/n/fff5/+/fuTlpbWYOdRlJaAGvmaKSoqYsGCxQBMnTqe\n1as3RBj5P5Kb+5dWZ+RzgeettQFjzB0A1tpbItq02GCoxqQ2BT6qKzCyYcMGnn32WW65Rf4MK1eu\npKysjB07dtC9e3duuOGGBrsGpfmi37fqKSoq4tZbfwnc42yZxKhRt7Bq1TcIBIYAm2jf/gc8/fRD\nrVeuMcaMBsZaa8dFbFcjXwOfffYZ3/rWt1i/fj0AX//613nmmWfo3r17lbYLFixg7dq1dOnShcWL\nF4ftCwQCXHzxxZx11lnMmjWLTz/9lN69e/Ppp5/Srl07unfvzoYNGxo00ZrSPNHvW/V063Y8Bw7M\nJDRr/wCf7zCTJ3/Gpk1z8fmOkJ8/odEMPCTOyCcyQdnVwKMJ7K/VEKuoxPe+970qbadOnUq3bt0o\nLS2tsu/JJ5/km9/8Jl988QUgecvXr18fTHNbWVmpX3RFiYk3M+gOunS5lvnztwItu5xijUbeGLMK\n6B1l18+ttc84bWYAFdba5QkeX5Ny8OBBCgoK2Lp1K/v27eO4447j/PPPZ8KECVHbN3RRCZdohnr/\n/v34fD569OgRNPJAMKPhmjVrGDFiRJXUwIrS2lm9Grp0eY4DByywFtgJTCQ/v3YZWJsrNRp5a21u\ndfuNMVcBFwAx3Vdmz54dfJ+Tk0NOTk6842syAoEAV199NfPnz+foo4/mvPPO49FHHw0r/hBJQxeV\ncImWlvipp57iuuuu449//GOVfY8++ihPPfUU8+fPr/XYFCVZWb0abr8ddu+GWbOO5d135/Gb3zwI\nwNSpN1dbJawhKC0tjfoLvb7US64xxnwbuAkYYa2NGVLqNfIthccee4ycnJxg3pz09HQqKiqqNfJ1\nJd6iEi6RM/lXX32Vs846K6buevnll3PRRReRnZ3Nc889p7N5pVXjNe4zZ8KVV0LbtgAFzJpV0GTj\nipwA33777Qnpt76a/L1AKrDKmV2+bK1NCveNdevWMX78eAAOHz7MwYMHg/lzYtEYRSWg6kz+X//6\nF19++SXFxcWsW7eO8vJy/vKXv5CSksIvfvEL1q1bR1paGr169eKJJ55g2rRp1V6HorRUSkpKmD9/\nEUBwodTdduDASVRWTqasrFuEcU9u6mXkrbUDEjWQ5sapp57KkSNHAFi4cCEzZ86s8Zi6yjXnnHNO\nWAWmDRs2cOedEmW3Y8cO+vfvH2bYI2frEydODL6fPXs2xhguvfRSiouLgzMDay3vvfcep556aq3H\npyjNFa9RHzFiKEVF9wYjVNeuzWPGjInMmfMfKip+B3SibdtbWbHiQi68cGQTjrpxaXYRr82JRx99\nlLZt2zJw4MAGN47xFtqIVWAE4PHHH+eOO+7AGENBQQHf/e53+d3vfseRI0fYvXs3AwYM4Mc//nGD\nXofSPGkJ37faUFJSQkHBXDZtepNA4NfO1huBawE3fqSYlJQuVFaeArgBkY0btVofmp2ffMwTtGAj\nryjJQjJ838Swz2P79l2UlR0AjgWm4U07IAZ+s/O5jI4d5/LFF6cAPwq2aW1GvllloVQURYlGSUkJ\nl1zyQzZuHE9Z2RygPfB+RKtewH3AfmAt7dufSEFBF/z+m5EHwFL8/unk50d3gU5WEhkMpSiK0iAU\nFMyjouJXhGbtALchEk0vYAgwCHgFyCMjow3Llz/IqFGjOOOMMzyLsY2XRbK5oEZeUZRmi6u9b9z4\nZpS9pwA3A/2BJ4GBwFf4fGUsX/5I0JiPGjWq1Rl2L2rkFUVploTyuR+HZE2Z7uzpBRwHLEFm848B\nC5B8M1OYMye/VRv1SHThVVFaAS3x+3b88dns2DEZWIkY9XJgHNAPeBjoR0bGnGBaYKDRk4g1JOpd\noyhK3LSk71tRURFz5txLRUUF8C1gKHA2IsusA24Avonfv7ZRKzU1Nklh5BVFaTyak5F3g5h27tzK\n3r2fcPjwEXr06AQE+OCDA0he916AHzgG2IvUJboZGEdGxtMsX35f0hp4aJ6phmtFc/oPpyhK4+G6\nQ1ZUnAPswS3U8cEHU5Eic2ORgnPdgELEyP8IGA5YYCGnnz4kqQ18ItGFV0VRGgV39r527RoqKnoB\npYiBd90iewHdgROAmUhy278B14f1Y8xb5OfPbqRRt3zUyCuK0uDI7P0yKiq6AEecrZFKxHAkaOkM\nIBvxnHkdmBRsYcxk5s6dprP4WqBGXlGUhBOZOOzOO++nouIQEmQ/EDHoryNeM2XO++8ibpEu28jK\nOpbx469g9eqVAOTnP6YGvpY02cKroijJSci//WzgRSCA1BR6HrgbkWU6IRGqc4ANyHxzk7MfYBKF\nhY1fuKM50eIXXhVFSQ4iZ+133HEf5eUHgVW4i6owBSgipL+vA74DvAn82tleAswG3mr1Bj6RqJFX\nFKXOFBUVMXPmXVh7DbCOVaueBw4D7ZAoVdeojySkxQNsRwx8Gl7NHd5UA59gVK5RFKVOlJSUcP75\nV2DteOBBRH4BeAOJSs1Ciscdh8zcxwC/dNpMolOntpSX+zhy5HN8vk507NiR6dMnqIF3aPHBUIqi\ntDzchGG7d3/I4cOHKSvrC3wAfEGoWMeTQAGSQOwWoCPwAFBOampnOnXqzNSp49WY14Bq8oqiNBol\nJSX89Ke3sGPHW0hZ5xTgKySfzGFkwdSVZsYgQUxfB85ECnscJiMjnY8//rDRx97aUSOvKEpMioqK\nuPPO+ykrOwichBh41wNmGnAOYui9NVNXAP/yfA4A3yYzc38jjFiJRI28oihAyEtm//6P+PzzL9i7\ndw/l5RWEPGQmIUZ+JTABuB2RZI5FZu5nIWkHJgHXIYFNkxH3ydVIJKvS2KiRV5RWSqTrY1HRvZSX\njwNeQvT1hUhKAW81psnIjN4PnI8kDitFDPpjSKBTLrAYaAOcjBj+PLp339XwF6VUQY28orRCQgFL\n4xDXx78D6UjE6WGkrF6biKN6Ib7vPYH3gB8CawklDrsD2Af0Bcbj8y0hELgOwKmturShL0uJghby\nVpRWyPz5ixwDvww4CpFh7gDmOy0OAd9EZu7FiAfNICTd70jEz/1IRK8fIgFOl+DzPcScOVPIzV1J\nbu7KpM773txRF0pFaUW4Es0rr/yHsrJDSCGOHYC3SPZS4FZElrkV8Xn/HZJQ7H5EjlmFLKj+1jnG\n1et74/N9wJw5U9RFsp6oC6WiKNUiPu3z2L17D5mZvRk79nxHd78TCVB6ANHcb404siewHDHue4E1\nwCKgM2LY/4UY+n8gM/1uZGX1o39/CYbKz79LZ+3NCJ3JK0oSEirM8SMk2nQb8CWSwrcvEpX6c2T2\nXoREoj4GnIb4v98B/B9SrANkdj8FOIgxaU7Rn28BF+P3T1c5pgFI1ExeNXlFacGUlJQwcuRYRo4c\nS0lJSXDbFVf81DHwy5DZ+tWIR8w+JGeM12d9BrKAeiLwC0SL/xuymOrFkJU1iEDgY4qLHyU316re\n3gLQmbyitEDc9AKvvfZfrJXgpNTUm5g160bmzLmLigqDeMcc7/y7ldCi6jRgBHAQcXXsBMwC2iKa\ne1sgh/AskpNISTnCX//6pBr0RqLJc9cYY+YClyAi3f+Aq6y1H0Rpp0ZeURJIyP3xOML92JcCbp1U\nEON+LfAIErA0GxiFaOlHIS6RhcBGYBiSZKwC+InT1whEf/+CTp0sTzyhBTsak+Zg5NOstWXO+4nA\nSdban0Rpp0ZeURLIyJFjWbXqOOBpoAdivEFm6B8gXjBdEHfGZcCdzv6ngIcQL5g/Iml+8519vYGd\nSIDT+0ilpnakpLTlyivPZ8mSJQ17UUoVmty7xjXwDp2QGb2iKA1ISUkJL7/8MhJlusDZ+gNksXSh\n83ka8BGy4Hon4YnD7kHcHq9EHgB5yMOiE2Bp06YUY9pw5ZXfU8OeJNTLhdIYU4SEvX2GiHiKotQD\nb6qB/PwJYfJISUkJF100lsrKtoRnfQSRaXojcgxIWoKuVE0c9ltEi18HjEPcKCXPTHb2UDZsWNsQ\nl6U0IdUaeWPMKuR/TiQ/t9Y+Y62dAcwwxtwCTCT0uzGM2bNDm3NycsjJyanjcBUleRG3x8uoqBB/\n8+ee+y79+x9D5849+Pzzz9ixYwuis0emGwDJGbMIMfI9EePdm/DEYdMRw/484ka5DvkBvo7U1Erm\nzZvZoNenVE9paSmlpaUJ7zch3jXGmGOAv1lrB0fZp5q80mqpbmYeydChZ7Nx4zZCxTcmIfr6tYhB\n3oLklTmClNdzPV9c430QuA2Rbu5ENPediFEfhOSYcWfug4FJ9OmTzimnnFHj2JTGpzksvA6w1r7t\nvJ8IfMNa+/0o7dTIK62SkBeMLHzWFDTUrdvxHDgwk5AMczSSQ+Yw8Btn243OtnbIjD4VMdrnISkK\nfok4u61H8suAaPQPAUcw5ius7UibNm0YN+4C1d2bMU2+8ArMM8acgPzeewfx5VIUxUGSgIUWPsvL\nZVssI5+Z2Y8DB9xPRcDnSP4Yr5vkA8BmZBbfEymtdzTwe8TAn4FEpvby9DwYOIrc3IE8++yTibo8\npYVQH++a7yZyIIrSmogm43Tr1hYx0LcCZUh+mDXO5z2IAX8TST9wvtPTOsT34V1CRbJTkchWN7Xv\njaSmGvLzXRlIaU1oxKuiNBCx5BogbLsxN2LtV4SklxJgFyLHeKsynYcUxh6E6O0ghnwh8hAodNpV\nABcCLwJH6NOnC4sXx/4FoTRPmlyTj/sEauSVVoy3pB6k0L17N/bv/5iNG8cj3i/zkFS/nyMz91Ik\n/UBkVaa9iPH+FzLbn+dsn4YkHgsgWSJTgU+Advj9bZgx40ZN+dtCaQ6avKIoUYgsq7d//8ds2rSV\nQOAqXK8WkV5uQvK4gxjr5wn5v69ENHeXXUhGyM+d/XchkalHEAP/NOI+uZTs7MVs2FDagFeotCR0\nJq8oCaSoqIhZs35NIDCAqi6LNyFpCN4B2iMRq+5MfRqSLGwgYuiPR3T595GZ/kRCOdz7AEOQYtof\nIh434n2jaX+TB5VrFKWRqcnnvaSkhAsuuJxA4NfOFtd//WFggLNtE1JL9StkFu9KNpuBOcBoQrKM\n12GtAuiASDGdCLlUTqJt2wDHHHM8/fv3V3/3JEKNvKI0ItUtohYUzGX37g85fPgwZWVzCM8KeRcy\n2z6MlG8oR2bp+5CZegckWdgQp00ZsAFJHFYOVAIXITP6bUji1+7AYjIyujJ16njV3JMU1eQVpREJ\nFb5eCUB5+TgKCuby3//+l4oKHyK9LIxy5HuIl8xvkNn6A4g0A+IK+RtEogEx9n9BDLlF/N/LEa3e\nIpGvoukXFt6sxl2JCzXyihKFyMXTV175N+L5Mh4xtNN4++22VFT0Q4x2byQSdTJizAcjWnkfxJ1x\nJfAK4hIZmRXSNdYW8aKZRKhYdiVwAiLvLKVNG8Ptt6uBV+JH5RpFiSAkzYxDyuDtJTxPzFLEAN+P\nzJMMYqDvRgz8/YgrY6WzzwDXIB4y9yLJwUBm7pMBV8OfjDws3gWOcf49H3jC2b+U3NyVGrXaSlC5\nRlHipDZJwsArzSwDjkMCkLxpfecB/0QCkjogeWRmIgZ6MuI5400y5nrGnIwsrp6OPBSmIdr7XMTr\n5jFcN0iZ+d+Cz5dPICDav98/nfx8N4pVUeJDjbyS1EQumK5dm1eji6EELr2KZHJcHKXFRqQotlu0\n46eIHr8PMfTeh0IusqC6B+gGXIBk5H4HCWLqhxj46wnlgg8xZMgpdO8u6wD5+eoaqdQeNfJKi6M2\nM/N4k4R5I1M3b95KyOWxktBCKYRqqJ6GzLbTkK+R6+54I+Ih47ILKcN3F+I3PwrxtpkGZDv9HXKO\nc5kG5OH3T2fePDXsSv1QI6+0KCJn5qtX/5CTTx5I9+69wgy+a7TXr9+EeKvE2+dCJK1AbyTxVxah\nCFQQeSZAyKhPRjxkImfu65Agpnwk3cBkzxmnOe0fcPZNRYKcbiU19StOPvkEunffpTN3JTFYaxv0\nJadQlMSQmzvGwhIL1nktsTDMwhLr9/eyxcXFtrCw0Pp83ZztYy10dtqF2sTu031faCHdeXUMHg9d\nI84/zMLfPZ/XWjjN6WeMhXynTUenr2EWij1jP8VCr2D/Pl/XKuNTWieO7ay3DdaZvJIE9AXyKC+H\ngoJ5bNr0uhN1uhlYAvSgffupdOjQlcxMKa1XVFTEggWLOXz4MF99VU5otj8B+C7yI/duZ9tkREev\npGq9+qXIwqs7c5+EuDu2d/YXA22R/O6fU1V734s3vUEgUH3OeUWpLWrklWZHdZp7fv4E1q4Vgy5M\nQ7xgSoCFbN68k0DgakRumY5IL5s5dOgBDh2ayYEDOMWwLXAqUjJvPyH/dhCj/GvCPWomI26QAMsR\nWSYNca3ch7g7vgWc4vTjyjmTCAUx3YBIM4Ixk2nfPsVzLYqSeNTIK82KeDT3FSuWOoukH7N5czmV\nlZORxcy7qawEMfyvIN4xvYFZSIWlPcAGKisN4Yulk4FzkepKKcgvg4WIDj/BaTMImbEXAQeQWqp/\nAs5BUg5sQXzj3yY84Amk6MefaN++PSeeeCK7d88lM7Mf8+Y9Bri55aWlukkqCScRmk91L1STV+Kk\nuLjYZmRk1ai5e9unpvYI7g8/pp+jx6d79PTOHo0836Obj7WQYaGn876H55juFr5tYZuFD53+l3p0\n9s4WOjnHh8ZadfzptrCwMOZ15+aOsbm5Y1SPV4KQIE1ejbxSJxJtmIqLi63f3yuGkRwTfJ+bOyZ4\nTGjB1LtwWuz0keYx8JF9DXeM9xLH2IcWZsVYuw+AnzvG/T0LP7Vwk9Oml+e4ThHnKfb07T4kOti8\nvLx63yOldZEoI69yjVJr6hJgVBMhf/beVNXCp0U/KMgE55jNhDI/3oWk5I1kExKEdJdzzFjC5ZU7\ngK3AIkRzn4lo/jlIOb2BSPrgZUie+Med/jY7fQGMAKaSltaR448/gXnzZupCqtJkqJFXak28AUZ1\nww3rn41o6NcgXi6f4PcvIz9/KSUlJRQUzGP79l3A35HUvT2AB522vYEPgKuQxVeXSUBXJLOjy0eE\n9PebkQjXkwnVUM0GngIuRlL9vo1o9kuRdYAHnfYP4K3Hmpc3miVLltT7bihKfVEjrzQLqnrNbEdm\ny6OAwWRkzGX5clmQvOiiH1BZeSJibPcTcklcjBjsNoiBXwacjSx8foF4uFyHGPubkNn3Fqedm8/9\nTmAYUmq7tYxhAAAgAElEQVTPpQuhWqpHEHfLD5EHyDVI6b3wxda9e1eiKM0BNfJKrYk0yInwCHG9\nZgoK5rJ589tUVvYO25+Z2dvxgX+NQKAtIc+YaUi06MeEcslMRgz1RESWcf3db0Jm+fcgM+9diKH3\nZoV8HHgESSAGkm6gLZJ6IBuRexYSPpt/ul7XrigNiRp5pdZ43RghsYmztm7dTmXlL51P44A8UlMf\nZPPmNlRWzkdmz652vwiRat5G/OG9Wv4URJO/m6oZJKcjD4S+iKF3jbxFcta8iTwo+iE+7g8gM/0t\nwI+Rmf/1iIF3c7+H1g3UDVJpTqiRV+rEqFGj6mzYYwU7RWr9AGlps4A0ysrmIhp9BTL7do09hOeF\ncTkBWUD10hPR+o9Dgp0eJDwt8BQk18z1zvu9wENOm91ItOtgZAa/yNnfHngfY76if/+7nTqrmnNG\naT6okVcalZKSEi655DIqKiS9wOrVl7Fy5WOMGjWKnTt3VmlfVlaGBBm5KQquRgzsyYgGX4nMuG/w\nHOUt7HEjYtyHIAFRdyCLqtmIgb4QMegnOH3fi/yCCCASzc8IPQSm4fNNcVImXEJq6k2eQK3ZatiV\nZokaeaVOxJqN15QGuKBgLhUVofwvFRUBCgrmArBzp5v7xa2V+rrzbyoy63b/u7YhPG3AMUg06vVI\nCuCJyIOgO/AkcCKimx+PaPUrAVcSmkt4CoPByC+D85DaqoM9ox/MkCEnBfO7jxhxI6tXb6jxXilK\nk5IIZ/vqXmgwVIsnMvApFLgUntmxuu3u8e3bZ1QJFvL5ujiRrv08wUjDLHSw0MXTtosNZYSMFnjU\nxUJbJ0L1fQufW/iZE8QUPbhKzhkZMHVaMDAqNTUUMeuNuI11rYqSKGguEa+IiBkAMmLsb8DboDQ0\n0YxZdvbwKobRNeKh7YUW+tk2bbrZlBSvoe5gJb3uGMdIL7HQ24bS+ObbUOrdaNGvbrreMRZGROz/\nu4UXHeMeeUxkagM3urVdxIPEfch0tzDIZmePiBrZGy3lsTcaV1HqS6KMfL3kGmPM0Yhrwe769KM0\nX6IFPu3ePbeGo4oQOeQejhwB8TxxXSI7EPJEyUP82H2Ivt4GSfrlni+ar7mPkFQzNWLf2cDtSEZI\nr597P0QC2oos3C4hVGT7Z0iA1FwkRXAlkja4EriQ7t13aeFspUVTX01+ARIm+JcEjEVpIWRm9qO8\nfLrjJ/8M8CKvvNKJMWO+id8/nfLytlTNxLjI+fcuz/bNhEeKTkVS/7oMRTR3l0lIMJN7fBmix5ch\n+v13kJQCywj5uU8mvED2QuQh4eZ9B1lsdT11JgGfAtcEI2yj0RCxAorSINT1JwBwKfBr5/0uVK5p\nsVSXbKw6nT0r6yQbntyrs83Ly7PGZFSRMkT77udIIe72aHJML4/O7s0WOcKz3W3/uYUbrCQcG+ZI\nLz1sSNNPt5JVMroe7/f3dSpIhY8hIyMrrsRrmj1SaUhoDLnGGLOK0O9sLzOAAmCkt3msfmbPnh18\nn5OTQ05OThyPH6UxqCnZWHWBT5988hWRM/bHH5+GtYcITyrmzsAHE5qZD0aKbETSCZmRr0SKcQxG\nFME8JBDJO3O/BPGcOYJIMQaRadY5n08GViMzeAgVGBH8fj+DBg1g48bwEZx++pC4JJr6xAooSiSl\npaWUlpYmvuO6PBmQ8jcfITP4XUjSj3eAnlHaNvDzTqkPsojqLmQWV1lArG62Gj33uzszLnb6HGZh\nUIxZ/VhPezctb6ENLby66XwnWPGW+dDCnyx8z4bqp/axkGVlMbe9s3Da2/OLwU093NvKom/oV0dh\nYaF6ySjNFpqLd41VuaZZURsJobi42Pp8XT2GT/Kku0Y+Ly/PiseLyCY+X1ebnT082G9hYWEVucbv\n7xPF8A+L+NzDwtFWcrF3tiHXSdcI5zufR1l4zsI7Vgpkfy+i73wbKqztes7k2WhukcZk2Ly8PJuR\nkWUzMrLCCnio7KI0R5qbkd+pRr7pqe2sNJoboM/XzRYXF3sMuKtvd3Nm3sPCdPmUlI7OLLm7hfY2\nI6NHFcMvxjtSJ3f3tXO2u5r6KRbGWXjeiua+zjnGHYvXLz5aUZB+NlSxaYlzTV1jVmWKvH9q7JXm\nQrMy8tWeQI18oxHLdzuW8YrWPjt7hLXWlWK8PuuuURUfeXdGXLVNO+fVz3kwZDpGvWtEPyOsyDiu\n3/tYK0FML1rY7Rh3r/F2KzANch4E6TGMfIbTpoP1+/vEbbBVtlGaG4ky8prWIMnZv//jmAur0dwA\n582Tohyff16GLGCGJwwTf3Y4cKCH89lt0xvJ8NgR8aztDVyGLJK+Tyh1QBHi6/4xkuDrHqTghx/J\nL7MVON1p+3+e856CZIH8FMkz8yGSeybSxbI9kE5KygesWLE47oXRhi2EoihNh6+pB6Akjvz8Cfj9\nbnKupfh8U9iyZTPl5cchxvkBysuPC+aKcT1ncnNXkpu7khUrxAtl9Og8Kit/SHTvl68Q3/M2iLHd\nhvjKXw6MR8rjgQQXDQKOctoAlCApgRcgni5nIwnCcoCXgeHIQ+EuJFfNUuc13dnnBi+VIsa4DPG8\nuRXxsR9MWlp7cnP78te//imqgS4pKWHkyLGMHDmWkpKSeG+torRcEvFzoLoXKtc0KsXFxTY7e4Tj\n/z3Whof0dw9KJ9nZI6LKEVVTE3hD/t3C1fmObJNuZQG1i+eYYkcTT4/Q4TvYqn7xn1v4pQ0tlHqL\nZGdZ76JvaLv0YUwXJw9ObM+gaPcmliSjco3S3EA1eSUWIUM9IopmHXJZjPSWCT/WNdidbWhhtIOj\nd+fZ8MXVyOAn1yB7tfrOFs634grptvujjZ4gzM1h4547PM+NMRl1cn+sKd+MLrwqzYlEGXnV5JOW\nzc4rkgG4unMgABs3LmT06Fg6/UJEM9+DyD0ZSNqBlYSCoEqQfO2TnPOtQ/LDlBCu5+c62+8EzkBy\n2k0mumKYgQQt3YxUdgpVYerUqSMDBgzmjDPOAGDQoOPZvXsumZn9mDevfsU6NLhJSUoS8aSo7oXO\n5BNGvDPN4uJia0xkRkfXqyU/YtacZWGYzc4eHjxe3CczbChgKdIl0p19F3tm7T1tuMtlVwv3es61\nzsLZFgY7vzAynL4HRel/kNO3dfrrZyHd+nyhYKbU1B4x0wBXd19UklFaCqhc03oQnX14WOBSTQYq\nLe0YjyEeY8XtsF0Ug5offJ+VdVLwQSLbCm3I9dH7YOjtMcaDbFXt/G+OLPOOlQjVP3oeBt2cMY21\nouV7ZZlhFtIixtjFdurUx2ZlDY4yjvAgq3hS/aoko7QU1Mi3EkKzz6rJvKpLpCXpCrzGs6Nj6N2E\nX66Pu9doZtiUlI7W7+/rHDPchnK9e9t1teC3Ib91d6HUTT/gtv2jM2OPtnjq6vhVdXI4pcq1RdPT\n62LkFaWlkCgjry6UzYxIF7+Q/3bfKm0PHOjBqlXHccEFVzJ0aE7QJVD83D9xWl3vvNogOdIHIyXx\nhhBe2g6gA5WVqZSXH4O4LL6GaOeTCLkzTgK+hfjDn4to5f91+iwAlgA/QDR5t/A1SPm9pc7nbc6+\nWOyha9d2YeUDI91DU1NvIjV1a/CzpPqdUE2fitI6MfLAaMATGGMb+hzJQmRGSL9/OoMGHc/Gjdch\nwUV5hPKeT0MWLkO50P3+6cyYMZGionspL+8KpCMPhwlIUeuHkFxyxzp97AR+57yfjORg/43z+SbE\nJ/4e3JzxskjaD+iP+Kx/gfizpwE/Bx5F/ODd8VUg2ScfIJSF8iYkU+QfgLbONS0lVCx7CvIQeR+/\nf1dYRszI+rFAtfVkFaUlY4zBWhszu2+8qHdNMyJa1OXnn9+Nz5dPIDAAOBufL5+OHTtQVpYHbMDr\nwVJeDgsWzKW8fBxi0G9xev4hYrCvQwyqpAFu02YKR478DIkm7QbMJDy69S5C1ZnSkOCjyYQiVM9w\ntgP8BzHw3uOnIrnrcpGZ/DokWGowEpla7mzr6PTbGakm2Q9YSXn59dVGnao3jKLUjMo1jUB9oix3\n7XqPQGA+cD0+3wvMmTOFP//5Afz+ZUhagGisI5RGIA/4FXASYnDvQn4VrOTIkRPIyOhCauoO4ECU\nfvYiOdsvQdIJPOz0dwESofp0DaPv4By7FrgPkY2OAJNo2/YLOnXqSJ8+/8OYzxBXyTnO+H6K/PoI\nUVRUxAUXXMmqVXtZteo4Ro/O04hVRYmHRAj71b1o5QuvtXHbi2wr3jRjncVJyZ/uLi7G8rgpLCyM\nWu0oFFQULUhprMfLJZYr42obXiA7sn2HKMe3tyG3TXfBdriF/GAitOgLqlJMOzU13WZnj7DZ2SOs\nMZ08fYenRFaUZAQNhmoZ1CbxlbcK0/79H7Nli+HQoeeRWS7ANPbvPyHYdtSoURE6dUi/vvXWGz09\nT0YWOocjMo47y3e5GZFyQPLEHHY+fwtoB3yJaPqXIUnIQGSfckRD74gof9cRkneuQ/T6fk7b7ohm\nL7VWd++u7lfAQDp1Wk5FRVs2bhzvbPsv8gvEvW8LibYYrShKOGrk4yRy0a+htGC339Gj8zh0aCAi\ncXgN8uLgeAoK5rF79x4yM3szduz5zJ+/iPnzF5GfP4GsrKPYsWMqooGfC/wDkXHCrgoxlgHEgF+D\nlN+7B5FkAP6JGNYNzueFSNRrGaLL/9rZ7maDdMvmLUUWZ6d59q9CvHGmkZkpD6v8/Ak8//zlBILO\nNtOBcaSm7uPgwcg1gkW4Rt7ne5v8/Nkx7qKiKEES8XOguhdJINfUJ1KyLseGJIzY+eFTU3tESCXp\nVvzgB1ljMpzkXZk2FKjkJhTrbKMX4Ohu4TwrkanvW7jPQolHZulsoY0NFeSIVoTbmwzNPYd3f3cL\nw2xqanrYPQhJTOJP7/f3stnZI2JITvEXAVGUlgwaDNV41JTYqiZqG2UZOp+bNiAUyp+dPdxTrMM6\nbdzIUTfIyBtN2t2Gol7doCVvbVX3mtwI1UlWgphco9zBwkk2pMG7wVHRtPTONpQOIbPK/rS0o2Pe\ng8h7FPlwdK9dI1WV1oIa+UYkHiOfyHD5cAOXb32+bjYr67SwXC0h4z3CVl1MdfPWuLPh4Y6xznBe\nbpm9v3muaZ2FvjFm5xnODN7NdVP1ASTjGeQ5vth6fymkpvYIuy/x3C9NQaC0ZtTINyI1SS7xSjK1\nLbKdmzvG8S5xZ+/R5It+Nrp0MsyGUhe0teGeL7kWVjkz93U2lH6gU4x+vF44x3r6cgtpn2JDvxa8\nx+eHpSfwXlNtk4spSmtDjXwjU52BjnemX1ttPvwYN/9LyJ0yZJg72ap517s4s/Autqo0U2ZFlrnZ\nhhfp7mDDi4S4vxa8s/olTtsMCz2qnb1XV5Qjsm91h1SUcBJl5NW7Jk7qG11Zlxqicsw4xC1xF/Am\n4vkC4q2Si3jNpCBeOABXINGtbRFvmcGI62Kup+cVwFPAe87npUh0airiibMQKf03npDLIrhBVMJ4\n4BUkh/w0T5tDGDOF0047NSy/e+T1CyFvGUVRGgY18gkgWkHs/Pylteojmovm/v0fIfVMxzv/uoU6\nQIzrn5ECGz9ADPAixGXyUySlwZ+df88HCoGvEUo45jX6OMddj7gwLkXcHR8klGDsBiQVgZsqYRri\n++4D+iB1Vg3wBNauYvfup4PXE/tBtpdQcrHa3S9FUeIkET8HqnuRJHJNTdSkt9elvmhW1mk25Mky\nLKYsIhq516VylIV9jizzO0dzb+9ILN2tLMTGlk5C2noPK8VAetnoeeV7RTkuumSj3jKKUjtIkFyj\nWSgbEe9sfcSIoaxeLQFG+/d/5GSadGfpS8nNXcn69Zs4cGAmIpEcRyhadSHhQVJfj/h8EMkz0wnJ\nMvkmMsvuDJyJzNKvQ4KjtgFXE8oCudTpf4tzjBttO4XwSNmlSCStu3+Sp8/ribyWZ599stECyhQl\nGdAslC0QV9cvKipi1qxfO5klh+PzvUi0eqyZmb05cGAScAwi13RGDPC+iJaDqaq5/xHYisg4fYEP\nEGnlRcQY73K2H4U8PFxZxk1TYIBBhFIJbHb2CcZMpnfvdL76Suqrjh17M6tXb2D9+n0ciJbrDM0a\nqShNgRr5BiZa+oFZs+YTCLjpAKYTCFyNz7eEQGAwsBmfbwn795/Cqacez8aN7sLmZsTAf4IY6WlI\nyt/TgBnAL4HTCRXI7o0UCSl3jvEj+d37IgbdO3P/r3NMPyRfzSpCC7x5hIp99MLny6dXr258+OER\nPvigNzCc8vJlnHFGATNmzPDkxJejVW9XlKZF5ZoGpKSkhEsu+SEVFb9ytkzDmCNYGyl7LCQ7ux2f\nf76PHTvexy3cIXnk5xNeMGQz8AYwCzgaWIYYdwvMBt4BsoBeiMSzDtjuHP8HxPD7CRn5SfTp052R\nI7/Bww//lUDAjyzSho/P59vGkCEnRX1IwThyc3fx7LNPBq9bZRlFqR+Jkmt04bWOxBPYFD2Nbj8b\n6e/u83WLkSLYXWyN7KfMwi8snB2lf7d2qrugOjyijywrfvUZ1pgMm5eXFzHWqmP2+/vYtLSjbUZG\nVsyC2urnriiJhab2kzfGzAauJSQQF1hri+v91GkBRJbpW7s2jxkzJgYXUqufvR5CyuGF/N2/9a0z\nWb16g6PRexkO3Ah8n6qa+wLEhfEmz3Z38dPri744os8DwG+RXwSLef31dyKKb0wIOz4lJZ/y8nLc\nMoGyRhC+fqAZIRWlGVPXpwNwGzA1jnYN9qRrKqLN0EOz8JDbYGFhoTUmw4aiVbtHdUX0+/vYNm16\nOC6J3ojTkTY19SULeyzcaEOJw9JtKAWBm78mw4ZHn4Zm2XJeN8FYvo3MO+MWG/Hmy4GuNivrNNup\nU58ofYYKlWhGSEVpGGjqmbxD/fWiJiZR+rHMwmUGXF4OBQVz2bp1O9YucFrcSFbWUfzvfwcpK/Me\nuZny8kPACcjM/QHgL8BE2rR5kCNHfoUENV2BzKa3IYunJ+LzTXG08a+Qmb2r3btMQsrt+ZHc8ZXO\n9kVE1oZdvXolM2ZMZNYst57s1ezdu4yUlHZVrjUtrRPDhq107tmjqrkrSnOmrk8HZCa/C9iEhEam\nx2jXoE+7+lDXPPHRy/SF506PllAsN3eMLSwstOEJvryJwy60sNXCbguP265d3Zm5O1vv6MzA+9mU\nlI62sLDQ5uaOsWlpx3jOVWxDKYejzezdXxZVxxbtF0pW1kk2sqyfztwVpeGhMWbyxphVyPQwkhnA\n75HKywBzgflIspQqzJ49O/g+JyeHnJycuB9CDUld8slAeJk+gL59L+Hhh5cQCKwDhuP3LyMzc1AV\nf/H9+z9m9eoNZGX14513buLIEUt4qoKxwL3AA/j9X3LsscfzySeincsM3iIzcKisnMTixcvp338Q\n06dPYPbsfCorNyPeNG8hRbTTo4y+Lz7fFqy9EWvd9tvo2/cS9u4tq9K6f/9BjB9/BQsWzAVg6tSb\nmTFjRrX3R1GU2lNaWkppaWniO07EkwI4FtgcY1/DPerqSX2LgVgbfVZfWFgYNYzfm15XNPzLLLxs\n4SMLbziae7rNyjrJ5uXlWb+/R8QsOjL9QL9g3z5fOxuu5/ewklXSm/4gw6akdHQ8ZNrZyEpOcs66\nVcBSFCWx0NSphoE+nvdTgOUx2jXgbagf9S3rl5s7Jqosk509wlorZe0yMrIc18OTItpttOH53KUK\nU3Z2tpOzpqMNVWHyyi1jPO+zqhj88Lan2E6d+ti0tGNsWtrRNivrJKds4LCokk1GRpYW6lCUZkKi\njHx9Fl7vNMachmgIu4Af16OvJiFSdsnPXxrXImK4C+XeKvs3bXqDoqIiioruDbpZfvrpFOBzT6t/\nIBGmPwT+DwjQp88v2Lp1r5NeeA/yAykSydwoi6o3e7ZXVmlpzF6eeGI5IDLU+vWbqKj4EfLnqjpu\n0NQDipJ0JOJJUd2LJp7J12dmGuvYcJmnOEL26GXdikjhM+U9VhZUZeZeNfDJu1gbvcaruC52sn5/\nX5uS0jG4XWbn/ghpJt1mZPR2XCAjC4EUOmMOz2SpC6qK0nygGczkmz3RgpZWrKjLbH0zzz9/JUOG\nnMK8eQURLUcha9MLEdfGpcCHwNNIbhmX3QwZMouePbsA0LfvRTz88BQCAdnr908nM/P4KMm9BgG3\nYsxnTiGOmYwaNSrC9fNhCgrmsXHjGYSKepzLgQMvAO2QmLU8T58LkTXy+2nT5iY6dGjP9Om6oKoo\nSUkinhTVvWjCmXx9FlZDx9YUOLTEpqamO7Np9/OF9vjjt1tj3g3O3L151bOzh3vcLocF0xqE1ghc\n18rQTDslpVu1+dcji39Xt2DrrbuqKErzBJ3JNxbRA4fCtfzHAJg5czU7d36fNm2+z803++nT51nu\nued+p81S/vOf/zgphv14c7gHAkudYKQZrFix1JmV+5z9ct7KSti4cSFwSdgvkpKSEn7605vYvfsj\n2rb14fffRHl5O8JdM0GSl32I3z+d5cvj+zWjKEoSkIgnRXUvmnAmX1/vGTk2euCQl3XrrB050trM\nTGsXLbL2q6+i9yez9+iza3eWnps7xmZnD496Xq9nTXb2CMcLR5KNebV1OCXKsZKmQGfvitIyoKld\nKOM+QQtfeA1JK1UfFP/8Z83G3SV6RsphUeUeOd9YW32Jvq5hfYT32yliUbVLMNukoigtg0QZ+aSX\na6pzCawpb417bHi7pXTuPIpRo2DbNpgxA/LyIDW1LqN7i+zsxcAQNm4cjyuvBAJuLvmrgIUYs402\nbQJUVn5IyH3SzTa5skqvaWldOf74Y9i9W6o2zZv3J5VnFKW1kognRXUvmmkwVF2kHHfmfswx1t5/\nf9WZe3W/GmJFxlobfZafnT0irC+375SUnjY8J03Vot7qCqkoLR9UrqkftfG8ice4VyfreNtFewjE\neuBIvyNsRkaWzc4ebouLix1f+vwIKaeDhXTr9/dWA68oSYIa+XoSj5GvybhbW7sF2uqIfABEy0Wf\nmppu8/LybCgv/DAL6bZPn2N1QVVRkoxEGfmk1+RjkZ8/gbVroxecfvllmD0btm4Vzf2qq8I1d69G\nv3//x07AVFVtvDZ41w5KSkqYNWs+UgsWpI5qHhUV69i7t4zCwptZsGAxAFOnTtMgJkVRYtJqjXy0\nvDXugmos4w5Vo2h9vnxnT3jZPO9DoyZKSkooKJjH7t17yMzsDaQ4xUC8fu4Lg+9mzJihhl1RlPhI\nxM+B6l40U7nGS7yyTCjrpLdASL5Hi5ei3NnZI2qUT9z+srOH25SUUOlA6G6NqVoiEDJsamq6yjKK\n0kpANfn6E49xt7bqwmh1gUzxGOHCwkInQdmwGIFLgyICp7rYrKyT1MArSisiUUY+qeWaWH7wNWnu\nkURWkBJm46YJmDcv/jQBrt4ucgxAfpRW6QwZ0o3u3d06qurnrihK3UhaIx8tA2VR0dMUFw+Latxr\nW9A7I2Mfp5++Mu4c9C7z5y+K0Ns3Azd6WkwjNbWSefMeU8OuKEr9ScTPgepeNJFcU9VFco9t3/5/\nMf3cqwuMqk8OnJrHtcQak2azsk4L84dXFKV1g2ry1VPVmK6z5577/TjbVvVxr20OnHgDn7yRr4qi\nKC5q5GsgZEyX28h87pHGNxEFvaOfO/YvA62jqihKdSTKyBvpq+EwxtiGPkcsInX2UD73AcBw/P5l\nrFghvuxe/d7vnx53Balo55k/fxGrVl1CSHdfSm7uSp599slEXp6iKEmMMQZrral3P8ls5L2UlJRw\nwQWXe7xapgPjyM3dxbPPPlnrhVdvv5EPiEGDBoVllVQjryhKbUmUkU9a75pIqnq1QKgua/UpiWvq\n1+teKWkSHsDvnx41ZYKiKEpj4mvqATQlPt/b5OdPiLt9SUkJI0eOZeTIsZSUlMRs1717L1askNl7\nbu7KWkk/iqIoiaRVyTXhOWemMGdOftw5YKLJMonQ8xVFUaKhmnwdqIvu7h6zfv0mDhyYSTSdva56\nvqIoSizUyDcC4bP3hcD16GKqoiiNgS68NgLhi6q9gXHBfbqYqihKS6BVL7y6xLegOgrIIyNjri6m\nKorSYmj1ck1RUVFYVkjvwmmsxVY17oqiNDTNQpM3xkwEbgCOAH+z1k6P0qbZGnkJkLqSQGA+sbR2\nXVRVFKUpaHJN3hjzTeAS4FRr7WFjTI/6DqaxkQCpAdW2qWuQlKIoSnOgPguvPwHmWWsPA1hr9yVm\nSI3NcCTFgeDzTSE//9GmG46iKEoCqc/C6wDgHGPMK8aYUmPMGYkaVGORnz8Bv38Z4jWzEJ8vnzlz\n8nXmrihK0lCtJm+MWYX4DkYyAygCXrDW3miM+RrwJ2tt/yh9NFtNHlRzVxSledLkC6/GmH8Ad1hr\nVzuftwNnWWs/jmhnb7vttuDnnJwccnJy6jxgRVGUZKS0tJTS0tLg59tvv73JjfyPgb7W2tuMMQOB\n56y1x0Rp16xn8oqiKM2RJveuAR4CHjLGbAYqgB/VdzCKoihKYmn1wVCKoijNkUTN5DWtgaIoShKj\nRl5RFCWJUSOvKIqSxKiRVxRFSWLUyCuKoiQxauQVRVGSGDXyiqIoSYwaeUVRlCRGjbyiKEoSo0Ze\nURQliVEjryiKksSokVcURUli1MgriqIkMWrkFUVRkhg18oqiKEmMGnlFUZQkRo28oihKEqNGXlEU\nJYlRI68oipLEqJFXFEVJYtTIK4qiJDFq5BVFUZIYNfKKoihJjBp5RVGUJEaNvKIoShKjRl5RFCWJ\nUSOvKIqSxKiRVxRFSWJS6nqgMeYx4ATnYzrwqbU2OyGjUhRFURJCnWfy1trLrLXZjmF/0nm1CEpL\nS5t6CFFpjuPSMcWHjil+muO4muOYEkW95RpjjAG+Dzxa/+E0Ds31D9ocx6Vjig8dU/w0x3E1xzEl\nikRo8t8APrLW7khAX4qiKEoCqVaTN8asAnpH2fVza+0zzvvLgeWJHpiiKIpSf4y1tu4HG5MC7AGG\nWjG58ZEAAASzSURBVGv3xmhT9xMoiqK0Yqy1pr591Nm7xuE8YEssAw+JGaSiKIpSN+qryf+AFrTg\nqiiK0tqol1yjKIqiNG8aJOLVGDPRGLPFGPOGMebOGG2+bYzZaox52xgzvSHG4TnXbGPMHmPMRuf1\n7Rjt3jHGvO60ebWZjKnR7lPEefONMQFjTEaM/Uc8Y3+6mYwpzxjzlvP6UQOPZa4xZpNz/SXGmD4x\n2jXafarFmBrzPv3KsQWbjDFPGWO6xGjXaN+9Wo6rMe3U94wx/3X+zwytpl3t7pW1NqEv4JvAKqCt\n87lHlDZtgO3AsUBb4DXgxESPxXO+24CpcbTbBWQ01DhqO6bGvk+e8x4NFFd3P4CyxrhP8Y4JyAB2\nINHX6e77BhxPmuf9ROD3TX2f4hlTE9ynXMDnvL8DuCNGu0b77sU7riawU4OAgcCLiDNLrHa1ulcN\nMZP/CTDPWnsYwFq7L0qbM4Ht1tp3nHaPAZc2wFi8xLsA3JgLxTWdqynuE8AC4OZGOE9tqGlMo4Bn\nrbWfWms/RSYaUX8dJQJrbZnnYycg0FDnipc4x9TY92mVtdYdx7+AftU0b7TvXpzjatTvn7V2q7X2\nrTibx32vGsLIDwDOMca8YowpNcacEaXNUcB7ns97nG0Nyc+cn2YPGmPSY7SxwLPGmP8YY65r4PHE\nM6ZGv0/GmEuBPdba12to2t4Y829jzMvOMU09pr7I/XFpjHtVZIx5F7gCmBWjWaPdpzjH1Oj3ycPV\nwN9j7Gvs756XWONqCjsVD7W6V3VyoawmSGqG02dXa+0wY8zXgMeB/lEGmVBqGNPvgTnO57nAfOCa\nKG2HW2s/MMb0AFYZY7Zaa9c04ZgaZFW8hnEVACO9zWN0c4xzr44DXjDGbLbW7mziMSWUmoIBrbUz\ngBnGmFsQeWR2lLaNdZ9qM6aEUtOYnDYzgAprbazAyYR+9xI0rsa0U94A05qo1b2qk5G31ubG2meM\n+QnwlNPu385CWTdr7ceeZu8jGqvL0YTPLhI6pojx/QGIejOttR84/+4zxqxAfq7V+T9aAsaU8PtU\n3biMMacAxwGbjDEgP2HXG2POtNb+L6IP917tMsaUAtlAnY1XAsb0PpDj+Xw08EJdx1PdmKKwHPgb\nUQxqY92nWoyp0e+TMeYq4ALg3Gr6SOh3L0HjajI7VUMftbtXDbB48GPgduf9QODdKG1SkAWfY4FU\nGn5Bo4/n/RRgeZQ2HXAWroCOwDpgZBOPqVHvU5Tzx1rkTAfaOe+7A28Bg5p4TF0R45nufd+A4xjg\neT8ReLyp71OcY2rs+/Rt4L9A92raNOp3rxbjapLvH7Lwenqi7lVDDLAt8DCwGVgP5Djb+wJ/87Q7\nH9iGrF4XNPBN+yPwOrAJeBroFTkmRFJ6zXm90RzG1Nj3KcoYd7oGFTgdeMB5/3/O2F9z/h3f1GNy\nPo8H3nZeeQ08jiec/+ObgL/gPLSBM5rqPlUzpqa8T28Du4GNzut3zvYm++7FOy7nc2PaqdHIGkA5\n8CHwj0TcKw2GUhRFSWK0/J+iKEoSo0ZeURQliVEjryiKksSokVcURUli1MgriqIkMWrkFUVRkhg1\n8oqiKEmMGnlFUZQk5v8BnvWFYVcFjOIAAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f872c6e90f0>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"pa.scatter_trend(np.log10(i['LFQ(NSAF)_x']), np.log10(i['LFQ(NSAF)_y']))"
]
}
],
"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.5.0"
}
},
"nbformat": 4,
"nbformat_minor": 0
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment