{ "cells": [ { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "# Python_version 2.7.15\n", "\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import pandas as pd\n", "import seaborn as sb\n", "import os\n", "import scipy.io\n", "from scipy import stats\n", "%matplotlib inline\n", "import octopus as oct" ] }, { "cell_type": "code", "execution_count": 23, "metadata": {}, "outputs": [], "source": [ "# Import Data\n", "\n", "R24E09_Ftrace_mean = np.load('R24E09_Ftrace_mean.npy')\n", "R24E09_ArcLight_MOI = np.load('R24E09_MOI.npy')\n", "R24E09_ArcLight_LDir = np.load('R24E09_LDir.npy')\n", "R24E09_R_deg = np.load('R24E09_R_deg.npy')\n", "\n", "R20D01_Ftrace_mean = np.load('R20D01_Ftrace_mean.npy')\n", "R20D01_ArcLight_MOI = np.load('R20D01_MOI.npy')\n", "R20D01_ArcLight_LDir = np.load('R20D01_LDir.npy')\n", "R20D01_R_deg = np.load('R20D01_R_deg.npy')\n", "\n", "R38G02_Ftrace_mean = np.load('R38G02_Ftrace_mean.npy')\n", "R38G02_ArcLight_MOI = np.load('R38G02_MOI.npy')\n", "R38G02_ArcLight_LDir = np.load('R38G02_LDir.npy')\n", "R38G02_R_deg = np.load('R38G02_R_deg.npy')\n", "\n", "SplitT4c_Ftrace_mean = np.load('SplitT4c_Ftrace_mean.npy')\n", "SplitT4c_ArcLight_MOI = np.load('SplitT4c_MOI.npy')\n", "SplitT4c_ArcLight_LDir = np.load('SplitT4c_LDir.npy')\n", "SplitT4c_R_deg = np.load('SplitT4c_R_deg.npy')\n", "\n", "rotations = [0,45,90,135,180,225,270,315,360]" ] }, { "cell_type": "code", "execution_count": 21, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "slope: -0.974543 intercept: 0.000114 r_value : -0.995686\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAMUAAADCCAYAAAD5PQZIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAGctJREFUeJzt3Xl0FHW2wPHvJUDYxo04ssQHPMDloYiIIIPLxAF0GHVYcgRxPKI4pQPIIOpRUEQ44oiAD3BQ0AHjuACuaCAP4zNqnMeILIqKyqIGiAgKIyIZgoHc90dVN52mO+kk3V1V3b/POTl2dddyg7n9q1/9NlFVDMM4qoHbARiG15ikMIwwJikMI4xJCsMIY5LCMMKYpDCMMCYpDCOMSQrDCGOSwjDCNHQ7gHgTkRGqmhfl47Rtvp89ezYA48aNczkSV0ksO6VcUgDt3Q7Aiw4dOuR2CL4hqdb3SUTuV9X7o3ycWr+sUVsxlRSmTmEYYUxSpIlZs2Yxa9Yst8PwBZMUhhHG1CmMdGLqFIZRFyYp0sT06dOZPn2622H4Qiq2U9TZkCFDEBFeeuklt0OJu8zMTLdD8A2TFCF2795NeXk5ZWVlNG/e3O1w4irNW7JrJW2TIjc395j3vv32WyorK7nkkkvIzs6mYcOj/zypWHoYkZk6RYjGjRvToIH9T1JaWsrPP//sckTxM23aNKZNm+Z2GL6QtiVFpG/+QOkxcuRIPvjgAwD++Mc/0qZNm+Bnfi0xjj/+eLdD8I20TYpIQv/gmzVrxjvvvMOTTz7J9ddf72JU8TFmzBi3Q/AN03hXjf79+/Ovf/0LgP3795ORkUGXLl2Cn/u11Ehjadt1PG6OO+44GjRowJ49e6ioqHA7nHqZMmUKAJMnT3Y5Eu8zSVGNQEmwefNmhg8fDsAdd9zBBRdc4GZYdXLKKae4HYJvmKSIwWmnnUarVq3YtWsXb7zxBgcPHiQnJ8ftsGrllltucTsE3zCPZGO0fPly8vPzASguLqagoMDliIxE8XxSiEiJiHwiIh+JyFo3Y2ndujWjR48GYM2aNb6qaE+aNIlJkya5HYYveD4pHDmq2k1Ve7gdSFZWVrDLxMaNG8nLy3M3oBi1a9eOdu3auR2GL3j+kayIlAA9VHVPjPsnZTxFWVkZM2fOBOxECZQgAWeeeSYAn3/+ebwuadRfyoynUKBQRNaJiOV2MAHNmzfn7rvvBmDPnj089NBDLkdkxIsfkqKPqnYHfguMFpGLw3cQEUtE1jp1jvOSFVhmZiYTJ04E7ClkpkyZgldL3gkTJjBhwgS3w/AFz98+hRKR+4EDqjqzun2SPRz1yJEjdOrUKbidmZnJTz/9BNgNgKHcup0K1H1GjBjhyvU9wv8t2iLSHGigqj85r/sDU10O6xgZGRlkZmYGJxzz4sRjaZ4MteLppABOAV4VEbBjfV5VV7obUmRffPEFqsrUqXbO5uXlkZmZaSraPuTpOoWqfqWq5zg/XVTV0wMCRITJkyfTokULwC4xysrKXI7Kduedd3LnnXe6HYYveDop/Or2228Pvp45cyYHDhxwMRpb165d6dq1q9th+IKvKtqx8NK8T88//zxbtmwBYOzYsZx44onJvLxxrJRpp/Ct4cOHc8455wAwd+5cvvvuO5cjMmJhkiLBBg4cSO/evQF4/PHHKS0tdSWO2267jdtuu82Va/uNSYok6N+/P5deeikACxcu5Kuvvkp6DD179qRnz55Jv64fmTpFEq1du5YVK1YAcPXVVwf7RxlJY+oUXtOjRw+GDBkCwAsvvMCHH37ockRGJCYpkuyss87i2muvBeD1119n1apVSbnu2LFjGTt2bFKu5XcmKVzQqVMnbrjhBgDefPNN3nrrrYRf86KLLuKiiy5K+HVSgalTuGjXrl0sWLAAgO7du3PllVe6HFHKM3UKr2vVqhW33norAOvXr2fp0qUuR2SASQrXnXTSScH2gy+++IJFixYl5DqjRo1i1KhRCTl3qjFJ4QHHHXdcsLPejh07mDNnTtyv0a9fP/r16xf386aihNYpnDEQ5ap6JGEXOfaavqlThDt06FBwWGvDhg2ZOHEiTrd5Iz6SX6cQkQYiMlxEVojId8AXwLcislFEZohI53heL9VkZmZyzz33AHD48GGmTp3q2eGtqSzet09vAx2BCUArVT1VVX8JXAS8DzwkIn+I8zVTSsOGDavMzzR16lQqKyvrfV7LsrAsz8z74GlxvX0SkUaqWu1MxLHsU88YfHv7FCp0FB/APffcU2VlpdoKzG6Y5o99Y7p9ijkpROQqVX29VhGIXAA8CGQCM1R1WW2Or4tUSQqwE2PatGkcOWJXySZMmEDjxo1djsrX6p4UIjI4wsnmAaMAVPWViCcTaaWqu0K2XwBudI5fpapnxxR6PaRSUgTMnj2bH3/8EbCHlTZr1qzW5zh8+DBAvUqbFFCvinbgj/kK4Ernv81DXkczX0QmiUgTZ3sfMBwYCuyPJSDjWOPGjSM7OxuAGTNmBKfPqQ3TThG7aEnRG2gKrAFuVNUbgD2qeoOq3hjtZKo6EPgIWC4i1wHjgEqgGTAwrpGnuNzc3CoruI4cOZLTTz8dgEceeSS4wlKsBg4cyMCB5n9BLCImhaquAfoBjYEiEelJjLceqpoPXAacALwCbFLVuar6fXxCTl/Dhg3j3HPPBeDRRx9l9+7dMR87YMAABgwYkKjQUkrUR7KqWqmqc4BrgTtiOZmIXCUi/wCKgE+BYcAgEVksIh3jEXC6u+qqq+jTpw8A8+fPZ8eOHTEdV15eTnl5eSJDSxk1tlOo6k5VvRroE8P5HsAuJYYA01V1n6qOB+4DPD1nk5/07duXvn37ArBo0SK2bt1a4zFmPEXsavMoYgXQvYZ9fsQuHZoCwakrVHWL875Rg0A9YuPGjVW2oepqrH369KFp06bk5+fz3HPPkZubW2Xl1mjnNWpWm6SI5XHWIOAaoAL7qZORQN27d6dJkya8+OKLvPTSS5SXl3PeeZEnXe/fv3+So/Ov2jTejVLVx2rYR7SGE8ayT32kSjtF4Js9liXEvvzyS5599lkAfvOb33DhhRces09glsLAlJ5pqn4dAp3OfQ2c142B90XkpBrO97aI3Coi/xF2rsYicqmIPA1cH0tgRuw6duzIyJEjAXjrrbd48803j9ln/PjxjB8/Ptmh+VLEpBCRgcC3wDci8nvgPWAm8LGIVNd55nLgCLBYRHaKyGci8hWwBfu26r9VNS+ev4Bhy87ODi4LvGrVKl577bUqnw8bNoxhw0y1LhbRunl8iL1yUFNgA3C+qm4SkXbAy7EsyCgijYAs4KCq7otv2NVeNyVun+rqhx9+YO7cuQB07tyZ4cNN1S5E/W6fVHWXqn4NbFfVTc5726o7Juz4ClX9tr4JISKXi8gmEdkqInfX51zp4MQTTwzeJm3ZsoUnn3wSgH379rFvX9K+m3ytupLiPFWtFJGeqvqB834GsEFVz0pKcPb1NmO3rpdidzu5RlU/q+aYtC4pAg4ePMjDDz8M2MNdA4vHPPHEE26G5bZ6lRQWdhcPAgnhOBVI5jKgPYGtzuItPwNLgN8n8fq+1bRp0+DCj/v376dhw4bBSdiM6kXt+6Sqx/QJUNUSVX022smcp0vx1BYI7cdQ6rxXJ5ZlBQfbHD58GMuyKCgoAOxuEJZlUVhYCNiPMC3LoqioCLBvPyzLori4GLCXCbYsKzjD365du7Asi9WrV9uBlpZiWRbr1q0DoKSkBMuy2LBhAwBbt27FsqxgI92mTZuwLItNmzYBduOdZVnB1uoNGzZgWRYlJSUArFu3DsuygrOYr169Gsuy2LXL7rm/atUqxowZE6x8qyoPPvggP/zwAwBFRUVYlhV8VFtYWIhlWcGuIAUFBViWFexynp+fX2Xk3quvvlql1+2LL75YpcV88eLFVWY5f+aZZ6qspJSXl1dltda//e1vVUYczp8/nylTpgS3//rXvzJt2tFOEbNnz2b69OnB7VmzZjFr1qzg9vTp05k9e3ZwO/TYmsR7OGpwqRwRKYzD+SIVd8fcArm1ZLAfZGRkMGnSJCoqKqisrGTGjBlxGd6ayuI9HHW9s+Y1IvKhqp5bz/P1Bu5X1cuc7QkAqvqXao4xdYoILMti27Zt/OpXvwLqP7zVp1yZIbCViIwQkXNjDaAGa4DOItLBaUAcBtRqSKxhGzFiBJMnTw4OZ502bZonlzb2gmhPn2ar6jgRySfCt6uqXhXxZCIW9i3U2UAX7AbAjc7PZ6r6cq0DFBkAzAYygEU1rZBqSoqaPfroo8FBSnUd3upT9RqjfZ6qrhORSyIdpKrvRjyZfbvzfqBvk4hkczRJzlLV62IMvs5MUkQWqIC3atUKgKeeeort27cD9nDX448/3rXYkqh+s3k4t0AdgY2qGtMK6SIyH/sx6mZgJbAydCKDZDBJEVngyVFoO8ULL7wQbL8YM2YMLVu2dCW2JKpXSXEf8AdgHdAL+IuqPhnzlUXOwO4mchlwPPYkaSuB/0v0FJomKSILPCru1atXlfeXL18efGxsWRatW7dOemxJVK+k2Ijd3+nfItIS+xv//DpFIdIUyMFOkt6x9JuqD5MUtVdUVMR7770H2BXydu3auRxRwtTr6VO5qv4bQFX3VrNf1StGaLxT1YOqWqCqtyY6IYzoSktLoy5XfOmllwYHIeXl5bF58+ZkhuY50UqKfUBxYBN7LtjAdnVPn4JtEyJSqKpJH+5lSorIItUpwn300UfBLueDBw/m7LMTPnddssVUUkRrvQnvXzQzxouG/tGdHOMxRhLcfPPNNe7TrVs3mjRpwtKlS3nllVc4ePBgWq69XesWbRFZqqpDo3y2E5iIPQbjKVXtVv8Qa8eUFPX39ddf8/e//x2AnJwcLr74YpcjipuEtWj3ruaz+4Ee2I1t2SLyiYgscabSHFKHaxlxUlJSEuxMWJMOHTpw0003AfD222+zcuXKBEbmPXUpKbar6n/UvKdpvPOSWOoU4b7//nsee8yeq+Lss89m8ODwebd9p+51ChGJNr+TAI1ijUBVS4FSEdkPxJRIRmKMHj261secfPLJ/PnPf2bOnDl88sknlJWVcd11Cf9ec120p09vV3eQqubUeGKRbthzP10N7AbOUNWE9yUwJUX8HThwIDhW4ZRTTgmO0fCh+C7aEtPJRE7D7sk6HDgAvAgsVdUSEflaVTvE7WLRYzBJEUFgsFKnTp3qdHx5eXlwUE+zZs2qDBjyEVeSohK7u/dIVf007LOvVPU/43ax6DGYpIigLnWKcBUVFTz44IPB7fvuu89vq7e6khSDsEuKi4H/xV78ZaWqVpikcFdg2Gt1883G4siRIzzwwAPB7Y8//hgRiWkmQw9I/iAjVX3VacPohN0B8GbsivZTwHHxvJZRO126dKl3QsDR4a0B27ZtS7lljRMyHlFVy4DngOdEpA32smDtE3EtIzaBCRECqyHVR4MGDfj444/Ztm0be/fuZe/evQwePJgGDRr4pcSoVryHo0ayXFUXxPLEykic8Nku6ktEaN++fXB7+/btKTMhQjJGrvuqJpaqbr/99rieL1Ai5Obm8s0333D55ZcDUFZWRvPmzeN6rWRLRkkR8+AkI3FOP/30uNw6RdK2bVs6dLCfts+cOdP303PG9emTF5inT5HF6+lTdV5++WU+/dR+Ej9q1ChOPtlzHaVdmeLG8Kg5c+YwZ86chF5jyJAhwa7mjz32GDt37kzo9RLFlBRpor4t2rXxzjvv8O679oQv119/fZUKucuS33jnBSYpvGH16tXBLudDhw7ljDPOcDkiwNw+GaE2bNgQnNw5GXr16sWgQYMAWLp0aVKvXV8mKdLEvHnzmDdvXlKv2bVr1+CSYsuWLeP9999P6vXrytw+pYnAqDs37u9LSkp4+ml7opeLL76YnBzX2nHN7ZNxVPv27V2r8LZv3z7YS7e4uDi4JohXmaRIE+vWrQvOBOiG1q1bB0f/rVmzxtN9pExSpIkFCxawYMECV2PIyspi3LhxgN2YmJeX52o80Zg6RZoIzA6YnZ3tciR2/6iZM+2pxLKysuo0fryO/F2nEJFfi8iPIvKR83Of2zH5WXZ2ticSAqB58+bcddddgL124EMPJXNt0Zp5Nikc76lqN+dnqtvB+Nnq1auDM497QZMmTZg4cSIAhw4dYsqUKZ4ZrOT1pDDiZOHChSxcuNDtMKpo1KgR9957b3B76tSpnhiT4dk6hYj8GngZe5ngncAdqroxhuNMnSKC8JWMvERVmTr16I3Avffey9Ch9syscX5K5e++TyJyHFCpqgecde/mqGrnKPtaQGCR529V9coop/XmL2sckxgbNmxIxPBW/1W0RWR0oGINtFDVAwCqWgA0EpGsSMep6hOq2sNZ/8K9h/EetmrVKlatWuV2GFGJCJMnT6ZFixZA9OGtubm55ObmJjQWTy2krKrzgHkAItJKRERVVUR6YifwXlcD9LFAm0BgHW0vCvyx79y5MzghwqBBg8jIyADifisVlaeSIkwu8CcROQwcBIapV+/1fCB0EjOva9OmDbt376ayspIdO3bQtm1bGjWKeQrjevNsnaKuTEU7NeTm5rJnzx5ycnLIz8+nTZs2NG7c+JhhtbUsPfxXpzASp7i4mOLi4pp39JCsrCwuuOACwL6lOnToUFKu6+XbJyOOnn32WQDfrUp02WWX0axZM4qKigCorKykadOmCa1fmNunNBGYduaEE05wOZK6Wbt2LStWrCA/P5/9+/fTtWvXuiRGvRaCNFKMX5MhoEePHjRp0oT8/HwqKir46aefEnYtU6dIE0VFRcFbEL8666yzWLx4MS1btmTv3sQ9nTclRZpYsmQJYC8k71eBdowDBw6gqlUa8eJZxzBJkSYeeeQRt0OIGxFJ6GIxpqJt+E6ghEhURdvUKdJEYWEhhYWFbofhC+b2KU0EvlX79+/vciTeZ26f0kR5eTlgj3hLY6adwjgqzZOhVkydIk0UFBR4fhIyrzAlRZpYtmwZAAMGDHA5Eu8zdYo0cfjwYQAaNkzr70FTpzCOSvNkqBVTp0gT+fn55Ofnux2GL5ikSBMmKWKXinWKEaqa53Ychn+lXFIYRn2Z2yfDCGOSwjDCmKQwjDAmKQwjjEkKwwiTVs2cIjICaJ/AS5yHtyd4Tvf4SmJ5XG8eycaRiKx1Zj73JBNfbMztk2GEMUlhGGFMUsTXE24HUAMTXwxMncIwwpiSwjDCmKQwjDAmKeJARH4tIj8GFrEUkfvcjimUiFwuIptEZKuI3O12PJGISImIfOL8+611M5a0arxLsPdU9Qq3gwgnIhnYi2v2w16TfI2IvK6qn7kbWUQ5qrrH7SBMSZH6egJbVfUrVf0ZWAL83uWYPM0kRfz0FpENIvI/ItLF7WBCtAV2hGyXOu95jQKFIrJORCw3AzG3T/GxHminqgdEZACwDOjsckwBkaZ18eJz+D6qulNEfgm8KSJfqKorK1eakqKORGR0oGINtFDVAwCqWgA0EpEsdyMMKgVODdnOBna6FEtUqrrT+e93wKvYt32uMElRR6o6T1W7qWo3oFKcVUREpCf2v2vi1p+qnTVAZxHpICKNgWHA6y7HVIWINBeRXwReA/2BT92Kx9w+xUcu8CcROQwcBIapR7oKqOphERkDvAFkAItUdaPLYYU7BXjV+V5pCDyvqivdCsZ08zCMMOb2yTDCmKQwjDAmKQwjjEkKnxKRd0SkR9h7gT5YH4rI5yIy2Xm/h4jMDdv3fBE5IiK5RCAit4T0RfqHiPxX4n4bbzFPn1LPe6p6hfNo8yMRWa6qa4FgJzunP9R07CdS0TyvqvOd/a8CHgEuT2DcnmFKihSlqmXYM2N0dEqQ5SEf3wq8DHxXzfH7Qzab481W8IQwSZGiRKQlcAGwMez9tsAgYH4M5xgtIl8CDwNjExGnF5mkSD0XiciHQCHwUISGutnAXap6pKYTOa32HYG7gHvjH6o3mTpF6qlpXEcPYInTepwFDHBa4s8HfgfgdF0JtQR4PAGxepJJijSjqh0Cr0UkD1iuqsuwe/beE/JZZ1Xd4mz+DthCmjBJ4W8rRKTCef1P7BF28TJGRPoCFcAPwPVxPLenmb5PhhHGVLQNI4xJCsMIY5LCMMKYpDCMMCYpDCOMSQofEJFBIqIickbY+zeLyLchMxM+E/LZAhHpU90+Ifs2FZF3RSRDRE52esV+KiIDQ/Z5TUTaOK8bi0ixiKTkI32TFP5wDXYv12Fh73cF7g1MoKCq14V81gt4v4Z9Am4EXnG6flwDPA30Bu4EEJErgfUhM278DLwFDI3bb+ghJik8TkRaAJcAI7H/YEOdDXwU4Zgzgc3OH3nEfcJcC7zmvK4AmgKZ2LOUNATGATPCjlnmHJdyTOOdx4nIH4Dfquq1IrIeuElV1zuf7QW+ASqBPara13l/PLBPVRdF2yfk/I2B7araytk+Hngee4aNu4AuwI+q+nTYcRnALlU9OVG/u1tS8p4wxVzD0RV+XnC214vIqdh/lF0jHHMZcEMN+wRkAfsCG6r6I07HQBE5ETsxBovIk8CJwCxV/aeqHhGRn0XkF6r6Uz1/R08xt08eEzrzoDP2oScQmANpKTDUmXitK2FjJZzjmwEnOPf/EfcJcxBoEuWz+4Bp2Im4Drvu8WDI55lAeUy/mI+YpPCYsJkHrwAKVPWQ89nXwC7gQuy6QqQ/+Bzgbed1tH1Cr/cDkCEiVRJDRDoDbVT1XaAZ9u2X4iSQM4jpe1WtIMWYpPC2a4ArnQVNSkSkBDgTGI79Bx9pjYnfcrRkibZPuELsRAs1jaMDixYDI7CfZs103ssBCmL6LXzGVLRTjFMZ71Wbb3ARORcYH+VxbbRjXgEmqOqmOoTpaaakSDGq2r22tzSq+iHwtvNEqUbOE6tlqZgQYEoKwziGKSkMI4xJCsMIY5LCMMKYpDCMMCYpDCOMSQrDCGOSwjDC/D/ivCUA88nTWQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "slope, intercept, r_value, p_value, std_err = \\\n", " scipy.stats.linregress(np.mean(R38G02_Ftrace_mean,1), np.mean(R20D01_Ftrace_mean,1))\n", "print(\"slope: %f intercept: %f r_value : %f\" % (slope, intercept, r_value))\n", "\n", "fig = plt.figure(figsize = (2.3, 2.3))\n", "ax = fig.add_subplot(111)\n", "\n", "plt.errorbar(np.mean(R38G02_Ftrace_mean,1), np.mean(R20D01_Ftrace_mean,1),\n", " xerr = stats.sem(R38G02_Ftrace_mean,1), yerr = stats.sem(R20D01_Ftrace_mean,1), fmt='none', color = '0.1')\n", "plt.plot(np.mean(R38G02_Ftrace_mean,1),\\\n", " intercept + slope*(np.mean(R38G02_Ftrace_mean,1)), color = '0.5', LineWidth = '1.5', zorder=0)\n", "\n", "ax.axhline(y=0, color='0.2', linestyle=':')\n", "ax.axvline(x=0, color='0.2', linestyle=':')\n", "ax.set_ylabel(' LPi 3-4 \\n-$\\Delta F/F$ (%)', size=10)\n", "ax.set_xlabel(' LPi4-3 \\n-$\\Delta F/F$ (%)', size=10)\n", "ax.spines[\"top\"].set_visible(False)\n", "ax.spines[\"right\"].set_visible(False)\n", "ax.spines['left'].set_position(('outward', 7))\n", "ax.spines['bottom'].set_position(('outward', 7))\n", "yticks = np.array([-0.05,0,0.05,0.1])\n", "ax.set_yticks(yticks)\n", "ax.set_yticklabels(('-5','0','5','10'), size=10)\n", "xticks = np.array([-0.05,0,0.05,0.1])\n", "ax.set_xticks(xticks)\n", "ax.set_xticklabels(('-5','0','5','10'), size=10)\n", "for axis in ['top','bottom','left','right']:\n", " ax.spines[axis].set_linewidth(0.5)\n", "ax.tick_params(axis='both', which='major', labelsize=10)\n", "plt.xlim([-0.08,0.08])\n", "plt.ylim([-0.08,0.08])\n", "\n", "#plt.tight_layout()\n", "#plt.savefig('Fig1_R20D01_R38G02_Regression.pdf',bbox_inches='tight', dpi=600, transparent=True)\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 27, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "1.728869090460403e-08" ] }, "execution_count": 27, "metadata": {}, "output_type": "execute_result" } ], "source": [ "p_value" ] }, { "cell_type": "code", "execution_count": 24, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAK0AAAC5CAYAAABa6L/OAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJztnXt8XHWZ/9/PyWSSTCZJm1vvSXrjJqA0KQUpchWLqKwtC+IFFlzZRVh/yur+XF0V7/5kdXfZdRFcUcQVxAuIiogKRUCkV24tFGpJmza9JWnT3G/n+f3xPaedDsnMaeacSWZy3q/XvGbOmTPn+52ZZ77znOf7PJ+vqCohIbmENdEdCAk5VkKjDck5QqMNyTlCow3JOUKjDck5QqMNyTkCM1oRuVNE9onIi2M8LyJyq4hsFZHnRWRJUH0JyS+CHGm/D6xI8fzFwGLndh1wW4B9CckjAjNaVf0j0JHikEuBH6jhz8A0EZkVVH9C8oeJ9GnnAC0J2zudfSEhKYlMYNsyyr5R55RF5DqMCwFQrKonj3HOcE7awbaV+9a1sKa5w3wqAqc3VHJ50zwsa7SPflLgqWMTOdLuBOYlbM8FWkc7UFXvUNUmVW0C+jNt2Fab5s5m1u5ZS3NnM7bamZ5y0tHc3sOa5g7mTCthXmWMOdNKWNPcQXN7z0R3LWPSjrQicjdwo6p2Otv1wJ2qekGGbT8I3Cgi9wLLgE5V3Z3hOdNiq83PX/05G/ZuML9rhSUzlrBy8UosyZ8I4L6uAVCwRFBVOvuG2H9ogOdaDtJQVTqZR9u0eHEPngSeEZGbMD7nJ4B/TPciEbkHOBeoFpGdwOeAQgBV/TbwEPB2YCvQC1wzjv4fMzsO7WDD3g3Mjs/GEgtbbTbs3UDTjCYaKhqy0YWsUFtWBAIjts1Lu7vYeaCPvqFhHnpxDwPD9mR3E1KS1mhV9XYR2QQ8BrQBp6nqHg+vuzLN8wrc4LWjfrG/bz8Ih0dVSywQaOtry47R2jZ0bIPuPRCfCZULwPJ/hG+oKuX0hkr+8PJeXtnbRSxawOLaMo6fEWdNcwenz69kQU3c93azgRf34APAZ4CrgFOBh0TkGlV9LujOBUFNSQ2ocRPckRaF6pLq4Bu3bdj4A9jxNKiAKNSdCadd5avh2rbS3N5DfWWMpvpK+gdHmF8TZ1pJISICCvu7BvLXaIFVwHJV3QfcIyL3A3cBbwq0ZwFRV17HkhlLXufT1pXXBd94xzZjsBV1IBaobbbrl0P1Il+aSI4adA0M0z0wQnlxBBHBVgWBmrIiX9qbCLy4B3+VtL1GRE4PrkvBYonFysUraZrRRFtfG9Ul1dSV12XnIqx7jzPCOm2JZba79/pmtIlRA0uEEdumvWuALXu7KCsqPBz6aqgq9aW9icCLe3AcZop1hqqeLCKnAu8CvhR054LCEuvwyLq/bz9Adgw3PtO4BGofGWlFIT7DtyYSowYABZZFQ3UpZx9XzeyKEmrKiqZE9OA7mIjB7QCq+ryI/IgcNtoJC3tVLjA+bLJPW7nAtybcqIGtenik7RoYAiUvDBa8GW3McQkS9w0H1J+skNWwlxstOLQbRgZh2nwonweRKJTN8j164EYN1jR3oLbS3N4LwEMv7Ob+jbtY2jCd689ZRCSSuzFpL0bbJiILcaZIReQyIPBJgCDJWtjLjRZsfxo6/gI9+yFWC1Xzof7NUH+W7+EuyxIub5rH6fMrea7lIN0v7GbYtmk50Ieqct/6nYBww3mLcnbE9fKJ3YBxDU4QkV3AR4HrA+1VwCSGvYDgwl5utCAah6FeiM+C4V6Ilpn9Hdv8bc/BsoQFNXFmTStBBHZ3DlBeHGFaLEppNMLaHJ/O9RI92AZcKCKlgKWqXcF3K1iyFvZyowVDvebessz/lbvtY9RgNGrLiugfslFVxJnOBSiKWPkZp3WmbUfbD4CqfjOgPgVO1sJebrSgsNTc27b5kRTGYKjb16jBaDRUlbK0YTr3rd+J9A0BMKuimHhxJG/jtGXO/fHAUkyCC8A7gT8G2alsYIlFQ0VDsFO3brRg+9PGULt3G592sMv4tD5GDUbDsoTrz1kECGua27FtZXjEZkF1KXXTY4G2HSSSThZJRB4BVrlugYiUAT9R1VSlNIEhIuucFMXRmFz5tLYN7VuhdSOoQqwKCosDiRqkYnjY5rbH/8La5g6KIhbx4gjL5ldNxqQZT53xEj2oAwYTtgeBhnF0aGqRnGeAQtUiqFsWcLMm72Bf1wC1Tlx2x4FeXmvv4Y3zpmE5U7m5nDTjxWjvBtY4OQcKvBv4QaC9ygcS8wwQ2P0sPPtDM+oWlwWWKDNatUJ9ZeyoWTIrx5NmvEQPviwiDwPLnV3XqOrGYLuVByTmGfR2wKGd5oKsuAIq5vqeKAOvzztwR9RZFcVHzZLletKM1xqxZzETChEAEalT1R2B9SofSMwzGOwGG7DExGwDSJSB1+cduCNqNGLRVD+dx7bso3/Qpjhqcd7xtTmbNJP2v0lE/gHYC/wO+BXwa+c+LSKyQkS2OIIcnxzl+ToReUxENjqCHW8/xv5PXtzIQecO6O+E4R4onwMl0wNJlAETl1WU9p4Bdh7opb1nAEUTRlTB2POkuvg6ZrxED7YCy1S1/ZhOLFIAvAK8FVPEuBa4UlU3JxxzB7BRVW8TkZOAh1S1Ic15cyt60LENunYbd6B9K2AFlvw9PGzzTz97nrXbO8wgL9D+8hoqWp/hzGs/d5TbsOtgHzeet2iy+bS+RQ9agM5xdOB0YKszo4ZTwHgpsDnhGAXKnccVjFGNm7NYlvn7r15k8gw6thmXID4Dpjf4Xnaz40AvinLOcTX0D45QHC3g/i3CUOnMqXUhBmwDVovIr4EBd6eHGbHRxDiS4z03A484LkgpcOFoJ0rSPchCXUwAOAZ87mV/i6A89o1rfC+72dc1gCBUlRaZTxNMfBjy6kLMyye0A+PPRjGzZO4tHV7EOK4Evq+qczGVuXeLvH4uNUn3oM1D25OaOSUDR8Jh0+vNvQ8JNNXxKF0DQ7R09HKgd5AR2wYRoj17OL2hkl0H+2jp6GXXwb6crl5I69OO+8QiZwI3q+rbnO1/BlDVryYcswlYoaotzvY24AynHm2s8/ri09pqs+PQDvb37aempCbz3IObK8b/2jHPmdora/jkr31trvlrl/h6vnGQmU8rIv+uqh8VkV8yijGo6rvSnHstsFhE5gO7gPcA7006ZgdwAfB9ETkRKAb2e+l4JkwVwY58JZVPe7dz/6/jObGqDovIjcBvgQKMKs0mEfkCsE5VH8SIfnxHRD6G+WH8jWZhjaipItiRr4xptKq63rl/fLwnV9WHMEoyifs+m/B4M3DWeM8/XiZcsCMkIyZSNXFCsNVmcGSQtt42YpEY04umo2jmlQuj+Z9tW+GJW6BsLrQ8A50t7Nm9i1e7Szj72i9mHC1I9EHHyju47Z+uQoDVq1ePux2vqG0zuH07w/v3E6mpIVpfjwSQyTaljNb1ZdfvXU/PUA9P7XqKyuJK6srqaJzZ6H/lQuUCmLcM1t5pIgNSgIhiiW1ybH3MPUisDdvfNXC48vbbvpw9PWrbHPzpz+hdv/7wvlhjI9MuW+W74XrRPThZVUddNyHXcH3ZOfE5zI3P5UD/AVq6W3jHwndw5uwz/b8IsyxjmM/dC8XTIBpjRlWcGSODMNDle+6BWxs21oTBaGmLfuXTDm7fTu/69RTOno1YFmrb9K5fT2xpE0Xz5/vShouXkfbbIhLFrKHwI1U96GsPskiyL1tZUknvSC9FBUXBRQ269sJgrymvGR6A/oMgETMrFnC5TSJBiywP7zdBH3dUde9H2tog20arqstFZDFwLbBORNYA31PV3/nakyzgVuGO2CN0DXXRPdhN92A3lcWVwTRo23Cg2Yyo0TJjtKoweAiqjw+83AaO+LLb9nePmrboVyJ4pKYGMG6CO9ICFFT7P4HpyadV1VdF5F+AdcCtwGliKhw/pao/971XAVFXXsdptafxwNYH6OjvQFHihXF+9urPWLV4FQ0VDf6NuG7lwqYHwB6B4S6IFJmc2vgsOOHirJXbwNhpi37lH0Tr64k1Nr7Op43W12d87mS8+LSnYgSPL8FM575TVTeIyGzgaSBnjNYSi8aZjTy560nqy+vZ07OHA/0HeHTHo7QcauEt897i3wSDW7lQudCMtFYUBjphTpMx1vLsromSLJfkV/5BYsQg1tRISeMS7I4OCqqrJzR68F8YPa9PqWrf4c6qtjqjb07R3tdOabSUeCRO50AnFUUVWIMW8Wjc3wkGt3IhVgXT6qCzxfiSfe1w0qVZcQ0SSZRLSvRpM8k/yGbEIBEvRvt2oE9VRwCchJZiVe1V1btTv3Ty4fq13UPdR2a6BeLRON1D3f5NMLiVCyjMfCOUz4UD2+CMG2HheVl1DWDskFgmF2HZjBgk4uWT+z1QkrAdc/blJHXldSypXcL+vv109HfQ1tfGrNgsygrL/JVGSqxcOLjDRA9OunRCDNbFDYktW1DFgpp4xlGDlBGDAPEy0harare7oardIpK7Sg+AohQXFFMeLadnqIeO/g5ae1ppnOHjBINlmRmv+uVHEr+zqHXgEmRsNpsRg6Pa9XBMj4gsUdUNACLSCPSlec2kxZ1gmFE6g/JoOb955De0FbVxwzU3+D/BkFi5MAEEHZvNZsQgES9G+1HgJyLilsLMAq4IrkvBsrd3Lzu6dtA/0m9ittNHkGGh0CrMu7REt6R8dkUxh/qH6e4f4g8v76OpYTqLar3k8adGLItpl60itrSJkba2QCMGiXiZXFgrIidgNL0EeFlVhwLtVYAM28N09HdQW1KLZVkwCHapzZCds29pTNzY7Eu7u9h10Pw59g4Oc9+6nXxyxQm+jLZiWeaiK8ALr2S8JswsxUghRTATC6hqTqrMFFqFVBZXHokeRMHqs4gWRLPXiSytJVZbVkTXwDA7D/RSUVIImAm519q6aW7vycmiRvC+zOhCjGDHiLNb8SCNJCIrgP/AJIH/j6p+bZRjLscUOCrwnKomVzeMi1PuOsXbgVXm7trfXpv20BeufiGDHjlkaS0xMLHZhqoYL+8+ZJYbBeZOLyEejeRsJS54G2mbgJOOtaLA0T34Fgm6ByLyYJLuwWLgn4GzVPWAiNQeSxs5SRbWEnOxLOGKpfNobu+hrKiQeHGE8uIIrZ39OVuJC97itC8CM8dx7sO6B6o6CLi6B4l8CPiWqh4ASFXQmDekWkssABZUx7nghBkg0N0/TGtnf05X4oK3kbYa2OxkdyXqHqQrbPSie3AcgIg8hXEhblbVhz30KXfJwlpiiQQxEzbReDHam8d5bi+6BxFgMWa18rnAE07S+VE5u+MR60j2P5s7m7nj+TuOKmZs7W7liVueAODsT5z9uueuO/U6/2vGsrCWWDLpksNzDS8hr8dFpB5YrKq/d2bDCjyceycwL2F7Lq+XPdoJ/NkJob0mIlswRrw2qQ93AHeA0T3w0PbrGGtxkKcOPQWQvfVyJ8lMWS7jRYDuQ5hRrlJVFzoXT99W1QvSvC6CEaC7AKN7sBZ4r6puSjhmBUaU7moRqQY2Am9KJXaXiViHK9CRuDjI+eedD8Cjjz36uufybbIhB/Dks3gx2mcxF1XPqOppzr4XVDVtTMmR7vx3jugefDlR98BJJP8GsAITTvuyqt6b5py5o5oYcqz4ZrTPqOoyEdmoqqc5I+gGVT3Vj14eK6HR5jWejNbL/9/jIvIpoERE3gr8BPhlJj0LCckELyOtBXwQuAjzS/gtZnZrQka1cKTNa/xxDyYbeWG0Wco9yEH8UQIXkbMwsdp653gBVFWzW+SUL2Qx9yBf8TK58F3gY8B6jiTMhIyXLOYe5CtejLZTVX8TeE+mCqlyD0Kj9YQXo31MRG7B6Bsk5h5sCKxX+UyWcw/yES9G6ya5JF78KHC+/93JPr7L2KdjAnIP8o0pHT2YMBn7w9GDMPcgCX8mF0SkQkS+KSLrnNs3RCSAVTGyi602T7c+zSPNjxArjDGndA6z47PZsHcDOw4FvIKqW6XbcJa5Dw32mPDyad0JdAGXO7dDwPeC7FTQuCPsXZvvoqWrhY37NvJSx0sIcljGPtexbWXb/m7+vK2dbfu7se3c+kdNhRefdqGqrkrY/ryTRJOzuNoHdfE62nvbiUfi7O7ZzUxnZUPfVGYmiKD1DiYaLyNtn4gsdzecyYacFeuAI+LK04unMys+i+6hbnqHemk51BJcHm0WcfUO5kwrYV5ljDnTSljT3EFze89Ed80XvIy01wN3OX6sAB3A1YH2KmBcETpFObHyRGbGZtLS3cJVJ10VjIx9lglai3ai8VK58CzwRhEpd7YPBd6rgBmtiuGt9W/NC4OF4LRoJwtesryqgM8ByzEe0pPAF1JVFwSJ38uM5mOlQg77tL4lgf8O+CPwQ2fX+4BzVXXUFcOTXptWrMM57jJMnu5SVU1ZAxZUllfWJxkCxlVLzLEKXN+Mdr2qNibtS2U47jEFmBqxw2IdmHqwzUnHlQG/xqxyfuNEGG24Vm7m+LTwnT+piZjcg/cA9znbl2GMLB2HxToARMQV69icdNwXga8DH/fS4SAI18rNjGzL2Hs5498BPwIGndu9wE0i0iUiqS7KRhPrOGp1DBE5DZinqr9K1QERuc6dkcOj7sGxkGqt3JD0JMrYR+fOpXD2bHrXr2dw+/ZA2ktrtKpapqqWqkacm+XsK1PV8hQvTSnW4ZTx/BtmJfJ0fbhDVZsct8B3S3JDYLYaJWtb7byYZMgW2Zax9yT1KSIrORI9eEJVH/DwsnRiHWXAycBqU0nOTOBBEXlXOr/Wb8YS8sj1SYZskW0Zey8XYv8NLALucXZdAfxFVW9I87q0Yh1Jx68GPh7UhVi66EA+hsCCXG8hER99Wt+iB5uAk93qW+dv/QVVfUPak6cR60g6djUBGe1UjA5kO1brRg8ylLH3LXqwBagDXK96HvC8l5Or6kPAQ0n7PjvGsed6Oed4mIrRgcT8gyDWwk3GlbE/95prgCNr8gaBl59CFfCSiKx2RsPNQI2IPCgiD6Z+6eQgVXTAVpvmzmbW7llLc2fz4YuxXCdV/kGu42WkHXVkzCUSowPuSItCZUll3roN+Zx/4FXqcwZmsRCANbmm2D1WdAAlb92GINbCTYfaNrXDw1TYNgOvvRbY8kxeLsQuB24BVmO+8rOBT6jqT33vjQcyjR4kRgfW713P/VvvZ2587uHjdnbvZOWilTTNTDlLnRP4nX/w0gkn+tg7w4kvv5S46duF2KcxiSz7AESkBrM27oQY7XixxKKhouGoEXQstyFfJhXyTQHcxYvRWknuQDveLuAmPeGkQm7ixWgfFpHfcvTkwkMpjs8ZLLFYuXglTTOa8mpSId/xpHuQMI0rwB9V9f6gO5aiL7mvmpiHuLNiv/nSl5gzPMycRYuIzp5N0QknMLx7N9V//3dmOdLU+CaqjKr+XFVvUtWPTaTBhkxe3MWdN5YUsy9SQGlTI4WzZzO8Zw92d/fhpBo/8Lo2bkhIStxpXIBhhMFduxjesxdVRXt76V23nlhTky8hsNBoQzImMWHmxP4BKkdG6H/+BSI1NVgiRBYuZKC5mcHt2724CGkJjTYkYwa3b6dn3ToKSktZcfHFDB86RP+mTRRUV2HFSonU1qA9PSa/NhtGGyqBh6RjaO9ehrZvZ7C/HwB7YAC7q4vBHTsoKI2b+5ISCqqqfGkvVAIPyRgdGmK4vZ3IjBmIZWH39KBDQ4jqUeEAvxQ6QyXwkIyxolEiVVWMdHXB0BB2dzcUFVF0wolEysuxYjGGu7qwOzpg4cLM2/NwzGMicouInCkiS9ybl5OLyAoR2SIiW0Xkk6M8f5OIbBaR50XkD84avCE5RqS2lsK6OiLTKmB4GIlGEdvG7uoiMmsWVnk5IuJb+U1gSuCO7sG3SNA9EJEHk3QPNgJNqtorItdjSsmv8Nr5kMlBtL6eovnz6du4EYnFsACJxRhqbaV/yxYK4nFijY1E6/0Zk7ykJp43znOn1T1Q1ccSjv8z8P5xthUygYhlEWtcYiIG5WVILEZBeQUDr7xC6RnLKF2+3Nc0xSCVwNPqHiTxQWBU3zlo3YOQzInU1mLF40RmziIybToAVjxO6fLlFM2f72terRf34E7gRYwKOMAHMErgK9O8LqXuwVEHirwf436cM9rzqnoHcIdzbFbLy0PSo7aN2jaRqkoGXtmCFY8D4qtLkEiQSuDpdA8AEJELMTm756hq7hcwTTESZ8NUFUWIVFUz7bJVRBsaAqlc8GK0fSKyXFWfhGNSAl8LLBaR+Rjdg/cA7008wJFFuh1YkWslPCGGREkkV6hjqLUVRAIxWBi/EvjfpHuRqg6LyI2YVctd3YNNSboHtwBx4CeOyswOVX3XuN5JSNZR26bv+RcYbmujIB7Hqqg4WhLJhynb0fC8jthkUQIP82knB65b0P34agZeeRUpKaHw2PNnk8msRkxE3q+qPxSRm5L2m06rfvNYexSSP7huQXTxceiIzWBrKwNbt0JBAWXnnBPIBZhLKvfArTUuG+W5cESb4rhJ3VZBAcUnnUThnNkMbd9B2UVvJXbKKfSuW5eJuHJKxjRaVb3defh7VX0q8TnnYixkCpOslFhQXoFdWcnwrlba/vT04eOCEFf2cqb/9LgvZAoRra8n1tjIUGsrgzt3MtTaSrShgYHm5sDFlVP5tGcCb8bodiX6teWYaEDIFMatCYstbTqslDi0dy+D27ePLq7sYyQhlU8bxYSjIhzt1x7CrLsQMsVxlRKTDTJoceVUPu3jwOMi8n1V3e6sQqOq2u1rD0LyBtdlSBZX9juS4EXL62TgbqDS2dUGXK2qL/raE4+EcdrJTYbiyr4pgf8J+LSbRigi5wJfUdU3e+2Jn4RGm9f4JtZRmpj3qqqrORLDDQnJOl5yD7aJyGcwLgKYRO3XgutSSEhqvIy01wI1wM+B+53H1wTZqZCQVHhOmJkshD5tXpNxwkzKRUDCFMKQiSKVT3smpsbrHuAZPP4KQkKCJpXRzsSUf1+JqTj4NXDPWCsuhoRkizEvxFR1RFUfVtWrgTOArZh1bP/B68k9iHUUiciPneefEZGGcbyHkClGyuiBY1QrgR8CNwC3YqIIaUkQ67gYOAm4UkROSjrsg8ABVV2EWZH8/x1b90OmIqkuxO7CrBL+G+Dz45i2TSvW4Wzf7Dz+KfBfIiKaayGNkKySyqf9ANADHAd8xC2z4YjUZ3mac48m1rFsrGOcQshOzLKmbYkHich1wHXOZnGKNsOLxSlAqiyvTFPNvYh1eBL0SBTrCAkJcu0hL2Idh48RkQhQgSlRDwkZkyCN9rBYh4hEMWIdyRMWDwJXO48vAx4N/dmQdAS25oJHsY7vAneLyFbMCPueoPoTkj/kXO5BSEi4nmZIzpFXSzKJyMOMT7+2mqQwW8CE7Y1Om6quSHdQ6B6QNt0xbG+StRe6ByE5R2i0ITlHaLSGbM+2he1lQOjThuQc4UgbknOERhuSc0xpoxWRyvRHZQ8ncR5JyAMNeT1T1mhF5O3AfSKy0Mkwm+j+nAc8KSIVqqpBG66IXCAiXxeRj49SUZLpuQPt+5Q0WhFZgSnt+U9V/YuqDk90n4BFmCT5/04w3EC+H+cH+6+YpbLmcmRhQz/OfbjyREQ+ciw1hV6ZckYrIosxFcY3qeovRKRURKpEpElEUi2DGjS/Az4PDAEPOPtifjfirOt2A/ARVf0PzKLaM0VkmYgcl+n5Ewz2w5gq7j9kes5kppTRisg84GOYnItFIjIX+ApwF3AvZjXKt0xAvyygCFiIkaHaJCIbgY3OD8rP72kn5gf7hIhUAR/HlDhdgnGXGsdzUhEpdO4LnP4uA/5WVTcnPOeL2zCljBaYjvkbbgbqgZcwauffBt4F9AILst0pVbVVdQuwG1OCdCfQAPSqaruq2n604/x1DzltgXmv/6Kqf62qnwV+CSwfx3nPB74rIrNVdQRTx3cc4P4AXPfrjMzegWFKGa2qPo+pLv4rzKLUF6rqB1X1V6q6GbN86gLI7hV8Qlt9mBH/RxgXZpOI3O9GFTIluSpEVdeq6i8SdhVyDC5Jwj/AFcBZwNdFZL6q9mJW43y/iFzi+OfvA74jIjWZvQvT8by+YRTM4wnbZZiR9S2u++XcXwU8BxyXxb4VJN0vAP4EXJZwzJwM2ygaZZ81yr4rgPXA8cdw7gbn/u0Yf/x6jBRAzPmcL8FUW/+P89m+wZfPbaKNKmCjmAY8hrlSvtTZJxg/9n+d7ZnAh4AX/fpQU/TnIuCrTvt1iQaE+Vv+KFDhbBf60N75wGpgUcI+t703Ax91Hl8OrDmW948RYXkVI/16Csa1eAPwSeBR4HfOcccB84FZvn2OE21YQd8wFzdXYSqBv4LRJ4sAvwf+GlO/1gjMC7gflwAbMKo6twK3Jjy3zHnuEp/bfA/Q5bzvhQn7G51R9R3OdjEw+xjO+zaM6MoZznYc+I7z+B2YFZAeB8qD+CynTMKME865DKN8UwRsB3pU9R+z0PZc4DbgFlX9o4i8G7gQE9p6AagDpqvqb31u90Tg/2J+mAp8GLPU1kXAPlV9VEQiegxxahG5CKMK/wTwKVV9xdn/RYwrdjHm32Qmxt25TlWH/HtX5P9ImzRCuL7jlzC+YxvG95IA26x27hc595UYQ70HMwLuAOqd5zLuBzDDPRdm4uDrmB/p1zGx4I3AzPG0B1wAvAK8D7gJ+BpwjvPcKszofYmzPQ+oCeQznWhDyuYt8UsCZrhfcIDtzQM+R8LFELAUWJGw/U3gZh/b+wwQS9h3O8a3Pxc4gPFxK8Z5/qXAm53HxwNfcAy3EeNyNWTje5wy7oFLNgXuRCSG+TIXACeo6r3J/RCRzwJDqvpVn9t7g6r+r4j8H+CNmBjppzE+Zz/wMVUdHGc7lqrazuziBzA+8U9VdU3ie8v0/YzFhCeKZJtsGKwzq7YE2INR0TkbaBKRflV9wO2HiFwBvBPzxfvd3lIR2Y1Z3OXDwGdV9X4R+QXmH2ZcBuv03XbuXxWRuzHTte93jPWZoD/jKWe0QeMk43wDo+n7ZcwFyZ3eTSkmAAABzElEQVSYiYNLnC/2fhG5FLgRuEZVXw6gvX7Mj+FuzCivIlKkqgOYmTdfcAz3x8C7ydZSXdnwQabKDTgV4ze+29l+E7Df2V+DiQffhtHlLSRDn9pDe9dhJlKuzMJ7zziu7PU2paZxs0Ax8GPgJBGpUdVngYcxF2L7MQH4lzDhrkJV3Rtwew9i4qnLRCSeYVspUb/DWimYchdiQSMiZ2GSb/qAWswV/Sp1fEgRqQUGVfVgLrY3GQiN1idEpEBNhhMishzjT74NOEtVd4lIVDO4+Jno9iYToXuQASKyVES+B2Y1IDfrSVWfxFwM3QdcKyJz/TCgbLc3WQmNNjOeA04RkdvBhIISDOkZ4GeYGbAP+ZRemO32JiWh0Y4T5+95EJMtdbyIfBcOG5JrMEWY7Klb3b/yXGlvMhMa7Thx/p5dQ7oIWCgidyY8dz1mbbQnVbU919qbzIQXYhniXhCJWVfiEUyK4YuYKdNVThgqZ9ubjIRG6wNJhvQ4sBg4X015T863N9kIjdYnEgwpgpnp2pVP7U0mQqP1kcTYaT62N1kIjTYk5wijByE5R2i0ITlHaLQhOUdotCE5R2i0ITlHaLQhOUdotCE5R2i0ITnH/wcRax+Kbg//hQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Make Scatter Plot of MOI\n", "\n", "pos = [0,1,2,3]\n", "\n", "fig = plt.figure(figsize = (0.5*len(pos), 2.25))\n", "\n", "#plt.subplots_adjust(wspace = 0.5)\n", "ax = fig.add_subplot(111)\n", "w = 0.6\n", "wb = 0.6\n", "data = [R24E09_ArcLight_MOI, R20D01_ArcLight_MOI, R38G02_ArcLight_MOI, SplitT4c_ArcLight_MOI]\n", "colors = ['C2', 'C1', 'C0', 'C3']\n", "\n", "ax.axhline(y=0,xmin=0.0001,xmax=0.9999, color='0.05', linestyle=':', linewidth = 1)\n", "\n", "for j in range (len(data)):\n", " ax.plot([pos[j]-w/2, pos[j]+w/2],[np.mean(data[j]),np.mean(data[j])],\n", " color=colors[j], linestyle='-', linewidth = 4)\n", " ax.errorbar(pos[j],np.mean(data[j]),yerr=[stats.sem(data[j])], color='k',zorder = 0) \n", " for i in range(len(data[j])):\n", " ax.scatter(pos[j] + np.random.random(1) * w/2 - w/4, data[j][i], s=25, color=colors[j], alpha=0.5, zorder = 2) \n", " \n", "ax.spines[\"top\"].set_visible(False)\n", "ax.spines[\"right\"].set_visible(False)\n", "ax.spines['left'].set_position(('outward', 7))\n", "ax.spines['bottom'].set_position(('outward', 7))\n", "ax.set_xticks(pos)\n", "ax.set_xticklabels(('VS',' LPi3-4', ' LPi4-3', 'T4c'), rotation=45, size=10)\n", "ax.set_ylabel('Motion opponency index', size=10)\n", "plt.ylim([0,1])\n", "ax.tick_params(axis='both', which='major', labelsize=10)\n", "\n", "bbox_inches = 'tight'\n", "#plt.savefig('Fig1_R24E09_R20D01_R38G02_MOI.pdf',bbox_inches='tight', dpi=600, transparent=True)\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 25, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAK0AAAC5CAYAAABa6L/OAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAHVBJREFUeJztnXmcnFWV97+/6n3vJN3ZyNJZIOxiFgICQVC2wYFhkUWRGWRgXFnU+cjrKKK+KiiMCvo6oLKIo8jLGhgV2bdhSwggiTBg0tkgSS9Jek1vdeaP+1RTaTrVT3c91d3Vfb+fT32qnue5de+pqlP3Offcc8+VmeHxZBOxkRbA4xksXmk9WYdXWk/W4ZXWk3V4pfVkHV5pPVlHxpRW0s2Stkl6fQ/XJel6SW9Lek3SwkzJ4hlbDElpJeWHKHYrcGKK6ycBewePi4GfD0UWz/hjQKWV9ISkmqTjQ4GXBnqfmT0FNKYocirwa3M8D1RKmjagxJ5xT26IMt8H/iTpemAvXA95QQRt7wVsTDreFJx7N4K6PWOYAZXWzB6S9BngYaAe+KCZbYmgbfXXXL8FpYtxJgRAoZkduIc6/Zx0dtOfTryPMObBN4AbgGXAVcATkk5OSzTHJmBm0vEM4J3+CprZTWa22MwWA7siaNuTxYQZiFUBh5rZc2Z2I3ACcFkEbS8Hzg+8CIcBO83MmwaeAVHYKC9JJWbWGrpi6XfAh3FKvxX4JpAHYGb/IUnAT3EehjbgAjNbEaLeFUGP2x/ePMhuQpkHAyqtpMOBXwGlZjZL0geAfzGzz6Uv4+DxSjumicamBX6MMwkaAMzsVZx96/GMCKEmF8xsY59TPRmQxeMJRRg/7UZJHwIsmAm7BPhrZsXyePZMmJ72M8DncY7/TcAhwbHHMyKEmVyoBz45DLJ4PKHYo9JKuoEUo3EzuyQjEnk8A5DKPFgBrAQKgYXAW8HjEPxAzDOChPHTPg4cb2ZdwXEe8GczO2YY5OtPHu+nHbtE5qedDpQlHZcG5zyeESGMy+tqYFXQ4wIcjQuc8XhGhFCxB5KmAkuDwxciCk0cEt48GNNEZh4A5AB1wHZgH0l+GtczYgxoHki6BjgbWA3Eg9MGPJVBuTyePRLGpv0HYIGZdWRaGI8nDGGUdi0uDtYr7VCJx6FxLbRsgdKpMHEuxMZWygmLx+lcv57uujpyq6vJnz0bZegzhlHaNuAVSY+SpLh+Riwk8Tis+jVseA5MIINZh8MHzx8zimvxODvuupu2lSt7zxUvWkTlmWdkRHHDKO3y4DFoJJ0I/AQ3kPulmV3d5/os4DagMihzhZn9YShtDYa4xdnQtIG69jqqi6qZUTaDTc2beo9nlc8ipoi+7Ma1TmErZoFiYHF3PPtIqJofTRv9EI8btQ2tbGvuYHJZATWTSojFQg3OB03n+vW0rVxJ3vTpKBbD4nHaVq6keMliCubMiby9MAEztw2lYkk5wM+A43DRYS9JWm5ma5KKfR2408x+Lml/4A9AzVDaC0vc4tzz1j28vPVlECRcfhY32nra2NW9i0VTFnHhQReSGwvznx6Ali1BDxv8CRRzxy1bM6a08bhx54qNvFjb6IbMgkNrJnLW4pkZUdzuujqA3l418dxTXw/DqbSS7jSzsyT9hX78n2Z28AB1Hwq8bWZrg/ruwCXoSFZaA8qD1xXsYTVulGxo2sDLW19meul0YorR2N7Is5ufZWrJVJq6msDgvrfvA+Cigy9Kv8ctnepMAou/19Ni0LULap/JiI1b29DKi7WN7FVZREwibsaLtY0cOmcic6tLI2snQW51NeDMhERPC5BTVRV5W5C6p700eP7YEOvuLxnH0j5lrgL+LOmLQAnw0f4q6pP3IK1voq69DkSvMrb3tNNlXWxp28KU4ilIwjBe3vYyG5o2UFNRk05zTiFnHQ7rn4OOZuhqh7xCeO0OUE7kNm48bry6cQf1zR2UFuRSWZRHTAKDuuaOjCht/uzZFC9a9D6bNn/27MjbghRKm1jObWbrh1h3mGQc5wK3mtl1wQLK2yUdaGbx3d5kdhNwE7gZsSHKQ9zidPZ0UtdWR3dPNx09HbR2ttLR3UEsN0bjrkbMDMMoyCmgvr0+faWNxeCQ86BlWzAY64Gtq6F6X5h2CGCR2bgJs+DRN7ayrr6VrU27mDGhmP2mlYGguqwgvc+yBxSLUXH6aeROm0rnunXkz5lD8dKlGfMmRGC07ZEwyTguJEhSZ2bPSSrE9aTb0m38oNsOGlT5HZ07djte37yexzc+vtu5v/zjX4YmTONaeHcVlEyGzhbIKXbnYnlQPt1N2URg4ybMggVTyojHYdP2dt7a1kwsJo5ZUE3cjOfXNkQ+MLN4nJ333Nvb07a/vpqmBx50F+XaiNKbkEmlfQnYW9IcYDNwDvCJPmU2AB8BbpW0Hy52ty6DMg0/8Ti8/GvYsgbyS6F7F7QGH7GnG7avg/xiKKlOu6ktTbtobu/iXWB6ZSHTKwtZ39DG8ftPZvXmnfzuxY3kxURRfg6HzZ3IZ4+eT25u+krU13vQ3dhIy9NPU3LUkeRNnBS5NyHMNO7HgD/0vWUPhJl1S/oC8BDOnXWzma2W9G1ghZktB74M/ELS5TjT4Z9srO0R1bjWPfKLobAcOvOgabOzZ3NyiSrGJx43VtY28saWZnJzYvTE40yrKGJqWT73rtrMKxt2EDfojhtFeTHWN7QC4vPHzE+7x+3rPYi3tWHd3XRv3IRiMXIqKoHovAlhetpzgJ9Iuhu4xcxCr8QNfK5/6HPuyqTXa4AjwtY3qrmqYvDvaXjrvdc/W9JPnTtDV1Xb0MraulaK82Js3rmLnrixrbmDzillbN7ZTl5OjM6eOAW5oqM7DhIv1TZS29Ca9uAsZ9Ikelpa6Ny0iVhJCZ2bNxNvaqJj0ya66urInTqV3IqKyLwJYUMTy3GDpgtwXcMtwO/MrDkSKQbBUEMTE/7ZlVtWsrp+NbXNtQgRZ/cbiBB5yqOmooYfHv1D5lbODSfYUJR2wDpTK23NFf8VaXO1Vw8ur6DF43SsW8fOu+6m9bnn6Glpwbq76WlvJ2/yZBTYs/G2NipPP52qz/zLQDZtqC4/lE1rZk1BT1uESz53GvCvkq43sxvC1DHSxBTj9L1PZ/GUxbzT8g7fe+F7bG7eDNCruDnkkJ+Tz4ySGZy94Oz0PQdjmMTUbfOTT9Lx5puQk0OstBQVFcE771C8eDGxWIx4WxvdTU0UL1o4fN4DSafgeth5wO24DIrbJBXjknZkhdKCU9yEIh5YdSD7T9yfN7e/SUtXC509nUwqmsTBVQdz/gHnM7dybnRTuWOQxOArVlKCdXVhbW30NDaSU1GBdXQQb20lb+ZMYuXlWDxO7uTJkbUdpqc9E/hRkI6+FzNrk/TpyCQZRura64jFYpQXlJPblMvMwpk0dzWz74R96bZucmO5g1fYvrfy2mfg5dvdDNi6p6A7iDVKzIiZQeUs5xKa9SE38bDwU1CTHSZ+YvCFxYm3txMrLkaS62lbW+mp20ZnkrsryomGMEr7bl+FlXSNmX3VzB6NTJJhpLqoGgxaOlves4INSvJKaOluiWZSoWQyNP7NTSp0tkJPJygXiidCewP0dDmfbfV+UFgBnU1QOmVQTSTboDc/s45bnllLHNjW1IEEPXFjYkk+BXkx9ptaTmlhLnOrSlhUM5Gp5YWD9tX+dd/99nitp6nJPTc0AND23PO915ruv58tV17Z7/v2e2PwGbbCdCfH9XPupEG3NIqYVT6LhVMW0tLZQlt3G00dTUwrnUZZfhkYVBVFOGceywuCZWJgXa5HNYOCCqe4bfWwc6Obyp0YctDXD/OqSyAmppYXMKW8gILcGDk5wszYe3IZ+04tY2dbF3eu3MQtz6zjp4+/zZ0rNhKPZ5+HMVXAzGeBzwHzJL2WdKkMeDbTgmWSxKBs4eSF3PvWvWxo3kBJXgnvtr7LwikLmVU+K/1GWrfBhHkwpRT+9gi0NUCPQVcrVM6G/f4BOnbA9rVw0Fkw75i0Yg+OmFfF4lkTeG5dA/EeozAvRs2kYqaUF7L/9Ap2tHfxzs5dlOTnUlGUz7TKwowG0WSSVObBb4E/4na3uSLpfLOZpdpqKSuIKcbcyrlcvvhyNjRtoL69nqqiquhiaUunuvuYgPwSKKx0YYq5+e7R2QQlVc50yCtMO1gmFhOLayayrbmD5o5uygpyWTJnAusb2jCgtaMbM0MSxQU5GQ+iySSplNbMrFbS+zIkSpo4FhQX3vMoRO7eSkR3rbkfOtugoBiq9oa2IMa1swWKKl2U1yBt2f6obWhl5YbtLJw9oTcccX1DG3OrSlnX0ErLri7aunqYX11KZVEecbNBB9H0Z392rFvHtmuvpWPjJnLLy3n1lVcoMmPfD3+YKV/58rAHgf8WF5a4kt5Q4l4MGLoBNh6IxVy4YflMeP6nMGEuFE2ELa9C/ZuwaycQT9uWTbCtuQMM14PinoVYUjORkw+exramDl6qbWRtfQubtrf3BobXTCpJq9382bMpXbaM7nvupXvrVvafOpXcSZMoPeqojIUmht4oZLSQdck6+q4RIw6T5jtlLZsWWQD42roWfvr427sFfm/e0c4Xjpnfe/tPLMGpa+6gOsJIL4vH6aytpf311xFQcMABFMyZM5TJhMg2CrkfuAO438zaBitF1GSd0kLSatytzhTIwGrc4V5ikyEiU9qjcck6TgZeBH4PPGhmI7IJXVYq7TCRqZ40DBEtIY9GaXsLuoWKxwIXASeaWfkAb8kIXmlHF4mgmR133UXHG2+gWA4qLKB02TIqP/7xwSpudAEzkoqAv8f1uAtxy74945zeoJknnqB91SriXV3klJQQKy+n+557KV6yhIK50Y/Xw+yN+3tcYMyxuCXh88zsi2Eql3SipDclvS3pij2UOUvSGkmrJf12MMJ7RpZE0Iwk6O4mp6wM6+4mVpBPT2Mj7a+/npF2w/S0twCfMLNBpawPk/dA0t7A/wGOMLPtkqILBfJknN6gGal3LRjgor7Mwt3rh0CqadxjzewxoBg4VdpdBDO7Z4C6w+Q9uAj4mZltD+pMe0GjZ/hI5DvIqapCRUVYRwd0dWEdneROmkTBAQdkpt0U144GHsPZsn0xYCClDZP3YB8ASc/i1pFdZWZ/GqBezyghke+gdcUKcqdMoXvrVmKVleROnkzhvvu6HBJBAo8oCePymmNm6wY618/7Pg6cYGb/HBx/ChdA/sWkMg8CXcBZuCXmTwMHmtmOPnXtlqzDzGr20Kz3Hgwzva6ubdvo2bWL9hdeoHP9emKlpYAGu3Q8skzgd/dz7q4Q7wuT92ATbtKiK/gTvAns3bciM7vJzBYHrq76EG17hgnFYhTMmUPJ0qUUzJpFd0MjBfssIH/GTPKmT6dt5Uo61w8130v/pLJp9wUOACoknZ50qRyXn2AgwuQ9uI8gy4ykKpy5sDa8+J7RxHAloktl0y7ABcxUsrtd24wbQKUkZN6Dh4DjJa3Bbaj3r2bWMLSP4hlphisRXRib9nAzey7SVtPAz4iNXiJIrhxZ7MFtwKWJwZGkCcB1ZjYiixq90o5uEgOznvp6cqqqBhuDENk07sHJo/lgEuCDYaXwjC8SA7NMJFNOEOYvEAt6VyeUNJHMJq7zeFISRvmuA/5b0l242+9ZwHczKpXHk4Kwubz2xwXMCHi0z74Jw4q3acc0kW4zOhFoDfJ21QW+V49nRAgTmvhN4Ku4aCxwG+H9JpNCeTypCNPTngacArQCmNk7uIQdHs+IEEZpO4Ps3AYgKb01xx5PmoRR2jsl3QhUSroIeAT4RWbF8nj2TFjvwXHA8bjR3UNm9nCmBUshi/cejF2iXY07WvBKO6ZJbxpXUjP9K4Fweb5GZAm5x5Nqx0bvIfCMSsLFi0lHSrogeF3lJxc8I8lQJhfy8ZMLnhEko5MLYZJ1BOXOlGSS9jTA8nh6ydjkQlKyjpOA/YFzg8CbvuXKgEuAF8IK7RnfZHJyoTdZh5l14tKFntpPue8APwBGJAujJ/sYUGnN7FrckvG7cYsdrwy5S2N/yTr2Si4QrICYaWYPpqpI0sWSVkhaAUS7Ss6TdYTZsbEEeMzMHpa0AFggKc/MugZ6az/nev2+kmLAj4B/GkgGM7sJuCl434qBynvGNmHMg6eAAkl74UyDC4BbQ7xvoGQdZcCBwBOSaoHDgOV+MOYZiDBKqyBt/enADWZ2Gm5gNRC9yTok5eOSdSxPXDSznWZWZWY1QZqj54FTzMz3pJ6UhFJaSYcDnwQSe7UPaFaYWTeQSNbxV+DORLKOYJNoj2dIhMl7sAz4CvCsmV0jaS5wmZldMhwC9iOPD5gZu/goL0/WEenCRo9n1OCV1pN1eKX1ZB1hJheqcak9a5LLj1QCOo8nTFqk+3Fp5R/B5ZD1eEaUMEpbbGZfzbgkHk9Iwti0D0r6u4xL4vGEJMzkQjNQAnTidqKBEVzY6P20Y5pokir7BY6e0UbYDZ1PAZYFh08MFP/q8WSSMAsbrwYuxW0Puga4NDjn8YwIYWza14BDzCweHOcAq8zs4GGQrz95vE07dok09qAy6XXF4GXxeKIjjE37fWCVpMdx/4RlvJcDweMZdsIsbPwdbinMPcHjcDO7I0zlA+U9kPQlSWskvSbpUUmzB/sBPOOPPSptsDcukhYC03BrvjYC04NzKQmZ92AVsDiwj+/CLSX3eFKSyjz4Em67+uv6uWa43W5S0Zv3AEBSIu9B7844ZvZ4UvnngfNCyOwZ56TKmnhx8PIkM9stkYakMLuQ95f3YGmK8hcCf+zvgqSLcX8g8HkPxj1hvAf/HfJcX1LmPditoHQesBj4YX/XzewmM1scuLrqQ7TtGcOkSqo8FddbFgWZYBJKWA4Uh6h7oLwHiXY+CvwbcLSZdYSU2zOOSWXTnoDL/jIDZ9cmlLYJ+FqIunvzHgCbcXkPPpFcIPgz3AicaGbbBiW5Z9wSZkbsDDO7e0iVu5DGHwM5wM1m9l1J3wZWmNlySY8ABwHvBm/ZYGYpcyL4GbExTTRLyCV9D/iBme0IjicAXzazr6ct4hDwSjumiWwa96SEwgKY2XbAB4V7RowwSpsjqSBxIKkIKEhR3uPJKGFiD34DPCrpFtzt99PAbRmVyuNJQdgdG08EPoqzOf5sZg9lWrAUsnibduwSzXKbgL8C3Wb2iKRiSWVm1jx02TyeoRNm5cJFuGCWG4NTewH3ZVIojycVYQZinweOwE0qYGZvAZMzKZTHk4owStsR7E4DgKRcvO3oGUHCKO2Tkr6Gi0E4Dvj/wAOZFcvj2TNhZsRiuLDB43Gju4eAX9oIZWP23oMxTfreg2D1wW1mdh7hNrzzeDJOSvPAzHqA6mB3Go9nVBDGT1sLPCtpOcGmzgBm9u+ZEsrjSUUYpX0neMQIufu4x5NJ/O42ntFEegMxST82s8skPUA/yjBQsLbHkylSmQe3B8/XDrXyINDmJ7iVC780s6v7XC8Afg0sAhqAs82sdqjtecYHqZaQrwyenww2C8HM6sJWnJSs4zjcIseXJC03szVJxS4EtpvZfEnnANcAZw/+Y3jGE6kyzEjSVZLqgTeA/5FUJ+nKkHX3JusIpoETyTqSOZX3YnPvAj4iKZRd4xm/pPLTXoYLlFliZpPMbAIu2cYRki4PUXd/yTr22lOZYAPoncCkvhVJuljSCkkrgFSJQuQfWf0IRSqlPR8418zWJU4EKY7OC64NRJhkHaESeiQn6zCzA0O07RnDpFLaPDN7XzaXwK7NC1F3mGQdvWWC6LEKoDFE3Z5xTCql7RzitQS9yTqCaeBzgOV9yiwH/jF4fSbw2EgF4niyh1Qurw9IaurnvEhtVwLORpX0BVxUWCJZx+rkZB3Ar4DbJb2N62HPGfQn8Iw7sm5GzOPxu5B7so6wq3GzAkl/Ymj5a6sY3hSivr3+qTezEwcq5M0DBgzC8e2Nsva8eeDJOrzSerIOr7SOm3x72dOet2k9WYfvaT1Zh1daT9bhlXYUEQTO42OKUzOulVbSxJGWIYGkY4BnJFWYmWVacSV9RNIPJH2ln+1f0607o7KPW6UNdt65U9K8ICxypJmPC7L/f0mKm5HfJ/js1+K2ypoBnBVh3UpE6km6RNIXo6o7wbhU2mDB5TXADWb2t2DVxEjzMPAtoIv38v+G2WRwUAT7un0euMTMfoLbVHuqpKWS9km3/iSF/Rxu37hH062zL+NOaSXtDZwLfMnM7pdUImmSpMWS+i4HGi6ZYrjNV+bh9rRYLWkVsCqQLcrfaRPusz8taRLwFdwSp5Nxd55FQ6lUUl7wnBPIuxT4ZzNbk3QtErNhXCmtpJnA5bhAofmSZgDfwy2uvAP4lqRlwy2XmcXN7E3cJoACbgZqgDYzazCzeBTtBLfurqAtgLnA183s42Z2JS6F65FDqPdY4FeSpgf53wqBfXCpAQASd7LD0vsEjnGltMAEnO1YC8zG7SVRCvwHcArQhvshh5WkHqgd9+f5Le5usFrSvQmvQrr0XRViZi+Z2f1Jp/IYhEmSdAc4G7cI9geS5phZG25z7vMknRzY558EfpFIR5AWZjauHriedjWwAFja59o1wLcTZtkwyJLT53kubof3M5PK7JVmGwX9nIv1c+5sYCWwYBB11wTPf4ezxz+LSwVQjMv7djJutfUvgVeBAyL53kZaiYZBMSYCpUnHZbiedVlwnJjKPj/4YvfJoCzHA9/HmSSzkhUId1u+DKgIjvMiaO9Y4AlgftK5RHsfAi4LXp8FvDgYpQJOAt4CqnH7Gz8AHABcATwGPByU2weYA0yL7HscaaXK5AOoBB7HuXdODc4pUJr/DI6nAhcBr0fVE+xBlpOBl3FZda4Hrk+6tjS4dnLEbZ4DNAefd17S+UVBr/qx4LgQmD6Iek8A1gCHBcelwC+C1x/DbSrzJFCeie9yzAfMSJqHs7euBm7FKfHjwJ9w20zdAxwCbDOzjXuoJl0ZZgA/B35oZk9JOg23meB9wF+AWcAEi3hTQUn7AV/FLSw14HNAPq7H32Zmj0nKtUG4/CQdj8vz9jTwNTP7n+D8d3B3tZNwd5OpOHPnYjPriu5TjaMor8AHeSYuXVMBsB5oNbMvZ7jdKjOrlzTfzN4OZuGexPXs63DJT44ys/XJjvk02ptiZluDwd1ewCXAN4DvAB/ELYU5ycy2DLY9SR/B/fm+BUzBbc31R3P53s4AvgZcaWb/FXhqdtkg8r+FZjhu06PlwXsDnv+LG/DU42zcjAy6cIlIvknSYAhYApyYdPzvwFURtvcNoDjp3I04M+nDwHacjVsxxPqXAB8KXi8Avo27gy3CuRFrhuN3HDc9LbxvinEKgJltzWB7xbgfcy6wr5nd0VeWIKFfl5l9P+L2DjCz/5R0KfABnI/033A25y7gckvaH26Q7cTMLB5M1HwKZxPfZWYvJn+2dD/PnhgNc+7DRqAkMkcmlXUZsBDYgsuicxSwWNIuM7svSZazgb/H/fBRt7dE0rvAvThb9kozu1fS/cCUoSpsIHs8eH5L0u246drzgu/2hUwqLIwzpYX3O9ijJohruA74DfBd3IDkZtzEwcnBD3uvpFOBLwAXmNkbGWhvF+7PcDuulzdJBWbWgZt5i4RAcX8PnIaz0TPPcNgg4+UBHIyzG08Ljg8B6oLz1TjX2s9xeXnzcD1eJtu7GOeTPncYPnvafuWwj/E2jZtpCoHfA/tLqjazV3CutQJzo+gHcFPHH8X9yOmaKAO1txznT10qqTTNtlJiEbu1UjGuBmLDgaQjcHEM7TiX0EzgDAtsSEmTgU4z25GN7Y0GvNJGhKQccxFOSDoSZ0+eABxhZpsl5Vsag5+Rbm804c2DNJC0RNIt4LZkTUQ9mdkzuMHQncCnJc2IQoGGu73Rilfa9HgVOEjSjeBcQUmK9AJwN25q86KIwguHu71RiVfaIRLcnjtx0VILJP0KehUpoTAFuOip6xO38mxpbzTjlXaIBLfnhCIdD8yTdHPStc8CPwKeMbOGbGtvNOMHYmmSGBDJ7SvxZ1yI4eu4KdMzAjdU1rY3GvFKGwF9FOlJYG/gWDN7bSy0N9rwShsRSYqUi5vp2jyW2htNeKWNkGTf6Vhsb7TgldaTdXjvgSfr8ErryTq80nqyDq+0nqzDK60n6/BK68k6vNJ6sg6vtJ6s438BGXx+AX9IK3IAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Make Scatter Plot of LDir\n", "\n", "pos = [0,1,2,3]\n", "\n", "fig = plt.figure(figsize = (0.5*len(pos), 2.25))\n", "\n", "ax = fig.add_subplot(111)\n", "w = 0.6\n", "wb = 0.6\n", "data = [R24E09_ArcLight_LDir, R20D01_ArcLight_LDir, R38G02_ArcLight_LDir, SplitT4c_ArcLight_LDir]\n", "colors = ['C2', 'C1', 'C0', 'C3']\n", "\n", "ax.axhline(y=0,xmin=0.0001,xmax=0.9999, color='0.05', linestyle=':', linewidth = 1)\n", "\n", "for j in range (len(data)):\n", " ax.plot([pos[j]-w/2, pos[j]+w/2],[np.mean(data[j]),np.mean(data[j])],\n", " color=colors[j], linestyle='-', linewidth = 4)\n", " ax.errorbar(pos[j],np.mean(data[j]),yerr=[stats.sem(data[j])], color='k',zorder = 0) \n", " for i in range(len(data[j])):\n", " ax.scatter(pos[j] + np.random.random(1) * w/2 - w/4, data[j][i], s=25, color=colors[j], alpha=0.5, zorder = 2) \n", " \n", "ax.spines[\"top\"].set_visible(False)\n", "ax.spines[\"right\"].set_visible(False)\n", "ax.spines['left'].set_position(('outward', 7))\n", "ax.spines['bottom'].set_position(('outward', 7))\n", "ax.set_xticks(pos)\n", "ax.set_xticklabels(('VS',' LPi3-4', ' LPi4-3', 'T4c'), rotation=45, size=10)\n", "ax.set_ylabel('Direction selectivity index', size=10)\n", "ax.tick_params(axis='both', which='major', labelsize=10)\n", "plt.ylim([0,1])\n", "\n", "bbox_inches = 'tight'\n", "#plt.savefig('Fig1_R24E09_R20D01_R38G02_LDir.pdf',bbox_inches='tight', dpi=600, transparent=True)\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 26, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAALAAAAC5CAYAAACIuLWWAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAG5hJREFUeJztnXl03MWV7z/f1r7a2JJtjC3LNjYGQyBe2CFsk7AFxuwkE7Iw+IQwCeucSeblkXDeTAZImPfCm5nMkIEZQngJJOBANhgmCVsIOLbZbLOZWF5BkuVdq6W+74+qthsjq3+Sfi11S/U5p8+vf1vd++u+XX2r6tYtmRmBQL6SGG4FAoHBEAw4kNcEAw7kNcGAA3lNMOBAXhMMOJDXZM2AJZVKWirpVUmrJN3mj0vS30t6W9Ibkr6SdvxuSWskvSZpXrZ0C4wcCrNYdidwhpntllQEPC/p18DhwFRgjpklJU3w158DzPKv44Dv+W0gcECyZsDmRkh2+90i/zLgWuBTZpb01zX5ay4EfuDve1HSWEkHm9l72dIxkP9k1QeWVCDpFaAJeMrMXgJmApdLWibp15Jm+csPATak3b7RHwsEDkg2XQjMrAc4RtJYYImkI4ESoMPMFki6CLgPOAVQb0Xsf0DSYmCx3y01syMPJH7QDxAYTnqzhw8xJL0QZrYdeBo4G1ezPuJPLQE+4t9vxPnGKaYAm3sp6x4zW2BmC4CObOkcyA+y2QtR62teJJUBZwFvAj8DzvCXfQx4279/HLjK90YcD+wI/m8gE9l0IQ4G7pdUgPuhPGxmv5D0PPCgpBtxjby/9Nf/CjgXWAO0AZ/Pom6BEYLyOZxS0jLvSvRG/j5YAHLJBw4EskUw4EBe06cP7Lu5MtFhZr+KSZ9AoF/06QNLagEeo29/5FQzmxm3YlEIPvCIJpIPnKkX4tdm9oU+pUg/jKxSIBAzoRcikKvE0wshabakwwavTyAQP5kacbf7a5KSZGZ/PTRqBQLRyOQDnwmc7N+/kGVdAoF+k8mAbwHux/kjX8u+OoFA/wiNuECuMvhGnKRJGaVEuCYQyBaZeiGijLCFUbjAsJHJBz5a0s4+zgvo63wgkFX6NGAzKxgqRQKBgRCi0QJ5TTDgQF4TDDiQ10SeE+fntk1Mv8fM1mdDqUAgKpEMWNKXgW8AjUDSHzb2TYkPBIaFqDXw9cBhZtaSTWUCgf4S1QfeAOzIpiKBwECIWgP/CXha0i9xWScBMLN/PNANkqYCPwAm4dyOe8zsu5IeAlLxxWOB7WZ2jL/na8DVQA/wFTN7sp/PExhlRDXg9f5V7F9R6AZuNrMVkqqA5ZKeMrPLUxdIugtfs0s6ArgCmAtMBv5b0myfXy0Q6JVIBmxmqeTUVW7Xdme4BZ8W6j3/fpekN3DZJlf7sgRcxr40UxcCPzazTmCtpDXAscAf+vVEgVFFJB9Y0pGSXgZWAqskLZc0N6oQSfXAR4GX0g6fAjSa2Tt+P1J6VUmLfWrWZUBNVB0CI5OoLsQ9wE1m9jsASacB3wdOzHSjpEpcNsobzCw98OdK4Efpl/Zy+4dies3sHq8P3ogHRdKSrN+5nub2ZmrLaqmrriOhML6TL0Q14IqU8QKY2dOSKjLd5JcWeAR40MweTTteCFwEzE+7PFJ61ThJWpJH33mUFY0r3M/HYN7EeVw066JgxHlC1G/pT5L+p6R6//o6sLavG7yPey/wRi+9FWcBb5rZxrRjjwNXSCqRNB23VsbSiPoNiPU717OicQWTKyczpXIKkysns6JxBet3hgHGfCFqDfwF4DbgUVxd9SyZ05+eBHwGeN0vMwDwtz4N1RV80H3AzFZJehjXyOsGrst2D0RzezOIvbVtQgmQP+7P56tbkUwaDS2tNO3qZEJVCfXjK0gkIs3SyStG9Zy4hh0N3PPaPUyunExCCZKWZNPuTdRX1bN662o6ujsoLSzllENO4aLZ+eNWJJPGw8s2sLRhq/sUBMfWj+OyBVNJJo3fv7uFd5tbmVlbwUkzaygszMnnGnxqKUn/x8xukPRzem9QXTBA5XKCuuo65k2Yx3ObnttrrHPGzeGZjc/Q0d2BEsKSxpL2JcybNI8ZY2YMt8qRaGhpZWnDVg4ZW4aAltZOfrp8I9vaOlnesJ2Vm7fT49aIYs6kam49fy4zJ1TmZQ2dyYV4wG+/k21FhgvDSCaTtPe009bdxttb36alvYWDKw9GEmZGY1sjq7eszhsDbtrVCQbb27pYuWkHG7a109bZzTtNu9jd0U1xgXOX9iSNF9du5Zs/X8Unj57MZQum5p0RZ5pStNy/PcbMvpt+TtL1wDPZUmwoSDXikiRp3dPKmjVrSJYnqaiqoKKzgpLCEooTUQcec4eaymIatrSyra2Tpl1ddHX3UFyYYOf2bVhRBd2JBAmMAkFXj7GjvYuX1rZw7PRxzKitHG71+0UkH1jSCjObt9+xl83so1nTLAJRfeCj7j8qdtmvf/b12MscDPVf/WXsZTbcfl7sZfaDWHzgK4FPAdMlPZ52qgoIoZWBYSeTD/wCLp6hBrgr7fgu4LVsKRUIRCWTD7wOWCfp08BmM+uAveu+TQEasq5hINAHUX3gZcCJZtbl94uB35vZwizrl1GvgfYD//H9P/LNF77Jhl0bKFQhhtHd3e08rwQUJ4oxM0oLS5k3YR63LLyF+jH1sT9DtljTtIs7nniL6pICGnd1snzdNto7Okj07KG4vJI9PUZxYYLaqhIMuGz+FK47fVYu9ULEssTA3utSxgtgZl3eiPOW2rJaygvLAeixHoT2fmRlBWVIQogCFXBkzZHUVdcNo7b9Z0ZNJWfOmcDShq2UFRVQXJDASkopKiils9soKRTjKos5fFI1STMW1I/LJeONTFQDbpZ0gZk9DiDpQmBL9tTKPnXVdZwx9QzW7VxHV9L9NgsLCimkkNOmnLb3WNKSzJ80P29G4VIkEuKyBVM5dvo4mnZ2cuTkap5Y1UjTrg72JHuoKiti9oQqZk2sZPOODiZWlw63ygMiqgF/Ebc87D/j/po3AldlTashIKEE1xx9DQC/2/A7uq2b0kQpCKZWT6UgUUDSkmzevZnastph1nZgJBJiRm0lM2orOXb6OM47ejKvrN/Ob99spHFnJwmJTdvbOW76eOrHZwwuzEn6FQvhY3tlZruyp1J04sgPnIoHbmprorOnk5fee4n1O9dTUVyB0IgKr0zFSLy0toXWzm469iRZWH8Q137s0FyMh4jPB5Y0EfgWMNnMzvHz104ws3sHoWBOkFCCuuo6ljUuc3HBgCRqSmtYNGsR9WPqR4Txwr4YiSkHlZOQSJqxtqWN9dva8m4ELkXUb+Y/gSdxky0B3gZuyIZCw8EH4oKrpjDroFm0dLSQUGLEGC/si5FIyFVuCQkMmnd1Zrgzd4n67dSY2cP4rDxm1o2b+j4iOFBc8Jb2vG6nfogJVSUgSHq3MWkGgtqqkmHWbOBEbcS1ShqP9yslHc8ISnRSW1YL5vzhVFwwBjVlI2vOaP34Co6tH/ehOOF8bcBBdAO+CTflZ6ak3wO1wCVZ02qIqauuY97EeR+aG5dvfb+ZSO9aa97VSe0ImKmRsRdCUgI4Hjc/7TDcV/yWme3Jvnp9E+cqRaneiC3tW6gpq8nLaUQjjEi/qqhDyX8wsxMGrVLMhGW2RjTxrJXs+S9JF/uZxoFAztAfH7gC6JbUgfcUzaw6a5oFAhGIVAObWZWZJcys2Myq/X6fxivpPklNklamHTtG0ouSXvHpoY71xyXpbklrJL0mad6BSw4E9pFppc45fjuvt1eGsv8TOHu/Y3cCt/l0qrf6fYBzcIlMZgGLge/17zECo5VMLsTNwDV8cDZGCmNfZskPnzR71if12/+eVM09hn2poy4EfmCuRfmipLGSDvYZLgOBA5JpRsY1fnt6TPJuAJ6U9B1c7Z9KDnigzJTBgAN9kmlS50V9nU9P2BeRa4EbzewRSZfhcqedRcTMlF6nxTg3A0J61VFPJhfik347AVdb/tbvnw48jcuV1h8+i1swBuAnwL/795EzU8adXjWQ3/TZiDOzz5vZ53G14RFmdrGZXYxbBmAgbAY+5t+fAaSSWz8OXOV7I44HdgT/NxCFqP3A9fsZVCMwu68bJP0IOA2okbQRt87cNcB3fX7gDva5Ar8CzgXWAG1kznwZCADRh5L/CdfF9SNcbXwFsMbMvpxd9TLqFYaSRy7xxUIASFoEnOp3nzWzJQNULDaCAY9o4jXgXCQY8Igm1mCeQCAnCQYcyGuCAQfymkwjca/Thy9pZh+JXaNAoB9k6gc+32+v89vUkgOfxvXXBgLDStR+4N+b2UmZjg01oRdiRBNrL0SFpJP3liydiJuhEQgMK1GHkq8G7pM0Blez7cAtfhgIDCuRDNivVnS0pGqc2zFikpoE8ptILoSkiZLuBR4ysx2SjpB0dZZ1CwQyEpL7BfKakNwvkNeE5H6BIcGSSbrWraO7uZnC2lqKp01DicEPBIfkfoGsY8kk23/6CG3Ll+89Vj5/PmMvuXjQRpzRgH1yv1LcVKCcSu4XyA+61q2jbflyiiZPRokElkzStnw55QsXUDJ9+qDKzmj+ZpYE7jKzbjNbZWYrg/EG+kN3czPA3to2te3ZMvgE4iG5XyDrFNa6VZ4smfzAtqBm8FkRQnK/QNYpnjaN8vnzP+QDF0+bNuiyoyS4FjDVzNYPWlrMhGCe/CHVC9GzZQsFNTVReiFiTXC93MzmR9R1yAgGPKKJNRrtRUmxLewt6XpJKyWtknSDPzZO0lOS3vHbg+KSFxi5RDXg03FG/K7P3/u6pNcGIlDSkbgEJ8cCRwPnS5oFfBX4jZnNAn7j9wN5hCWTdK5dS+vSpXSuXbu3sZZNojbizolR5uHAi2bWBiDpGWARLsXqaf6a+3G51/4mRrmBLGLJJNt/8hN2P/ss1tmJSkqoPPVUxl56aSwjbgciaob2dbjke2f4921R7+2FlcCpksZLKsellJoKTEylr/LbCb3dLGmxz+6+jJCdMmfoamhg+6NL6NywkT3NW+jcsJHtjy6hq6Ehq3KjhlN+A1cbfs0fKgJ+OBCBZvYGcAfwFPAE8CrQ3Y/77zGzBb7xNrKW0sxj2l9fSc/WrRRWV1NQXU1hdTU9W7fSvnJl5psHQdRadBFwAdAKYGabgaqBCjWze81snpmdCmzFZalslHQwgN82DbT8nCWZhC1roOF5tx0CH3EoMbO9XT/m97M98hXVB+4yM5OUikYb1Hw4SRPMrElSHXARcAIwHZc/+Ha/fWwwMnKOZBJe/gGs/wOYQAZ1J8BHr4Is+ohDRemRcykcP57kzp1Iwsx48/33+fodd/DQBRdkTW5UA35Y0r8BYyVdg5sP9/1ByH3Eh2fuAa4zs22SbvdyrgbWA5cOovzcY+ufnPGOqQMlwJJuf9rJUHPocGs3aEqmT2fMokXsfu45rKMDlZby4nubaSqMamIDI1NikxIz6zSz70j6M2AnLiLtVjN7aqBCzeyUXo61AGcOtMycZ/f7vub1ta0Sbn9344gwYCUSHHTZpVQcd+ze0bYXPvc5LMvhM5l+Hn8A5kl6wMw+g2t4BQZC5STnNiR7oGOHe+1+H7avhy0TYdyMvHcllEi48Mi0EMkJ3d20Ll0aaxB7OpkMuFjSZ4ETe1vwZQCLvIxexs2AqSfAKz+EXU3QtQuSe6BlLYybDtNOhPoToWryiDBmSyY5sb2dmV172P6IM5O4gtjTyWTAX8SlkRrLvgVf9upI/xd5Gb0kElB3PKx+BKomQWMznR1tqG0Xxbvfh03LYPOrUF49Ihp3nWvXMrejk3aJgspKVFUVWxB7OlGDea42s3tjkxoTeRXMk0zCU7fC6sdhz25oa6GzxylZWlQMGNSfDFMWwo71cMpf57Rv/Macw2Mv8/A330jfjTUarQK4Eagzs8U+duEwM/vFQBSNi5w14G+OyUKZuTWHNlcMOOp/1H1AF/tW1twI/F1kzQKBLBHVgGea2Z24flvMrJ2Iv5BAIJtEHomTVMa+vBAzgc6saRXIefb7u6fj3XfZfPPNdG3aDEVFWGcnXa2tAJRUVJCorKRgzBgSZWVMuOlGKo47LhY9ohrwN3CBN1MlPQicBHwuFg1GIgfyV3u64cmvw+pHoaSKjRvWUZZIMr6qFKoPhpJq6O6EqolwzrdzuhG3P5KguASKiqCjA8z2RmhVTJtGQXm5O9bYSLKrKza5UfJCCHgTF7NwPM51uN7MQiRYfykohMPPdb0MpWNYvWoLR4/d5QY3dr4HaoLCIjjiAtcXnEfsaWpCgsLqaqy4GARtHR0AjGlvI9nTDYjC8eNJlJTEJjejAfsgnp/5OXG/jE3yaKVqMpRWwZgpfPz8S+GtX0BXG5SNhWQ3FBTD5I8Ot5b9xvbsobtlK4UTJyKJrqYmKs1oB8ygoKyckjlzsLa2vdPs42BY5sSNasbNcAMVO9bDtnchUQyTPgKVE6CoDJIGL/2ri1zLo3DLRHGxi0bbvZvulhaS27bRJtFSUECiuJielhZ6mpsoX7Aglun0KaL6wKcDX5TUgIsJTuWFCKsU9ZdEwo2yTTsZNq2AN5ZAxcGw8UUor3VDzAfNyLtItcIJE1ysQ0U5ezZsxMxobWtjQ1ERCxcupLOhgepPfILq888f0qHkFHHOiQskEs4wx82A7jZY/Rh0tgGC6ilQPh62t+ZVpFrxtGmUL1hA2/Llzsft7mbOySdzzNy5YEZhTQ2lRx01tME8kkpx8RCHAq8D9/rcwIE4SNXG1VPhxX9yNW+5z2Irg8qJw61hZJRIMPaSiylfuIDu5mbali2ns6GBro0bSba2UlJXB2ZYMhmrEfc5lCzpIdzgxXO4WnidmV0fm/RBkrNDyf1lBM7WSE2x3/Hoo3Q1NJCorATUn4i0SANlmVyII8zsKAC/RsbSKIUG+km6X7y70dW8eR5SqUQCJRJ0b2mhZPZhsadVTZHJgPemUTWz7pCcMouk/OI88Xmj0Gda1SEy4KMl7fTvBZT5/ZCdMpCR9LSqqRoY4kmruldGXyfNrCA2SYFRRzbTqqaIFA+cq4yYRtwIZgBpVVPEF9CeqwQDHtHEGtA+ZEg6W9JbktZIChkqA32SUwYsqQD4Z1yf8xHAlZKOGF6tArlMThkwLmfwGjP7k5l1AT/GpV0NBHolu3l/+s8hwIa0/Y3AB0L3JS0GFvvd0j7KCp3Wo4BcM+DejO4DjTEzuwe4Z2jUCeQ6ueZCbMQlu04xBdg8TLoE8oBcM+A/ArMkTZdUDFyBW6M5EOiVnHIhfLzFXwFPAgXAfWa2apjVCuQweT2QEQjkmgsRCPSLnHIh4kTSEwxsFaMahnbxmCCvd7aY2dmZLgouxH5kiK8I8nJMXnAhAnlNMOBAXhMM+MMM9ShfkDcIgg8cyGtCDRzIa4IBB/KaYMAeSeOGW4d0fHB/Kr1t4AAEAwYknYtb5nampGEf3JF0OvC8pDE+vW1WjVjSmZLulHRL3DNgsq37qDdgSWcDdwD/18zezZHcb4fiAvn/Jc2Is/Jd+R/vd4BNuPDVy2IsW+Z7CSR9RdKX4yo7xag2YL9c2JXATWb2mKQKSeMlLZB0yDCq9hRwGy4z0s/8sfK4hUiaDlwHfMXMvgu8DEySdJyk2YMtP814vwR8CvjNYMvcn1FrwJKm4ta+KwQOlTQF+BZwP24u3m2STh0GvRJACTAT+AKwStLLwMv+xxXnd7YR9+N9TtJ44BZgPHAezqWaP5BCJRX5bYHX9zjgL81sddq5WFyLUWvAwEG4v+oGYBrwBlAJ/CtwAdAGDPlCFWaWNLO3gPdwU6zuA+qBNjNrMbNY0rb7v/c9Xha4Z/26mV1qZrcCPwdOHkC5ZwD3SppsZj24eYuzgdSPIeWiHT+4J3CMWgM2s9eAXwN/DvwHcJaZXW1mvzCz1UA73oCHsicgTVY77p/g/+HcnFWSlqR6JwaL7TeCZWZ/NLPH0g4V0Q+3Je2f4XLcKlZ3SppuZm3At4G/kHSe9+c/DXxf0uAXyzCzUfMCxgGVaftVuBr31JS75rdXAa8Cs4dQt4L9tjOAF4BL0q45ZJAySno5lujl2OXActxywlHLrvfbc3H++7XAT3E/giqcW7IB+Hf/2c6N5XMbbqMaQgMZC/wO1+K+0B8Tzu990O9PAq4BVsb1Afehz8eBf/Dy69KNCffXfQMwxu8XxSDvDOBp4NC0Yyl5JwI3+PeX4fJAR35+XCKad4Ba4Cic+zEX+CrwW+Apf91sYDpwcGyf43Ab1lC+cA2jq3Aznb8F/BmuEfffwKW4eXjzgalZ1uM8YAVwNXA3cHfaueP8ufNilnkFsMs/98y04/N9bXu+3y8FJvej3E8Aq4Hj/X4l8H3//nxgJ/AMUJ2Nz3JUBvP4LqJLcJmASoB1QKuZ3TwEsqcA3wO+bWbPSloEnIXrLnsdqAMOMrMnY5Z7OPA3uB+pAV8CinH/BE1m9ltJhdaPfnBJHwcewC1B8bdm9rY//r9w7to5uH+ZSTiXaLGZ7TlAcQNjuGvF4Xqxz9f8O5yvuQXnqymLMmv89lC/HYcz2h/hasb1wDR/btB6ABNTZeEGKe7E/WDvxPU1vwxMGog84EzgbeDTwE3A7cDH/LmLcbX6eX5/KlCbjc90VNbA8KFRookAZtaYRXlTcf26t5tZpz+2EBhvZk/4/X8EdprZN2OS9zngLnM9AUj6N1wtfAywBNeYutDMDrC4c5/lL8T55i9IOgxnyMXAT3y5U8ysYbDPkVGP0WrA8EEjHgJZ5Th/ewYwx8x+vL8ekm4F9pjZP8Qsb66ZPSjpeuBoXB/s/8D5qB3AjeaSKQ5ETsLMkn5U8zM4H/qnZrY0/dkG+zwHYtgDV4aToTBeP5o3D3gfl2XoFGCBpA4z+1lKD0mXA5/EGUHc8hZKeg9X634JuNXMlkh6DOdmDHj5ePMDK2b2jqQHcEPGf+EN96Vsf8aj2oCzjQ8Uugv4IfD3uMbMfbhBivP8l7xE0oXAXwGfN7M3syCvA/fDeABX+5ukEu/KvDfgB9wPb8QPAYuAtXGVm0loeGXhBXwE2AYs8vvHAM3+eC2uv/l7uPzHRfgGVxblLcYN2lw5BM8+6H7rqK9RO5Q8BJQCDwFHSKo1s1eAJ3CjYc24zv43cF1oRTb4BmQmeY/j+muPk1Q5SFl9YnF3lfXBqG7EZRtJJ+ECg9qBCbjupIvN+5ySJgBdZrY9H+XlAsGAs4CkAnORWEg6Ged/fgI4ycw2SSq2QTSchlteLhFciJiQtFDSfwCYWU8qOsvMnsc1pB4GviBpShzGNNTycpVgwPHxKnCUHyzAXN9oyqheAh7BjbxdE1NI5FDLy0mCAceA/wvvwkV1HSbpXthrVCnjKcFFed2d+rvPF3m5TDDgGPB/4Smj+jgwU9J9aeeuBf438LyZteSbvFwmNOJiJNWYklvf479wYZErccO2F/uurbyVl4sEA46Z/YzqGWAWcIa5KUx5Ly/XCAacBdKMqhA3wrZpJMnLJYIBZ4n0vtmRKC9XCAYcyGtCL0QgrwkGHMhrggEH8ppgwIG8JhhwIK8JBhzIa4IBB/KaYMCBvOb/Az0MRRWennN2AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Make Scatter Plot of R_deg\n", "\n", "pos = [0,1,2,3]\n", "\n", "fig = plt.figure(figsize = (0.5*len(pos), 2.25))\n", "\n", "ax = fig.add_subplot(111)\n", "w = 0.6\n", "wb = 0.6\n", "data = [np.radians(R24E09_R_deg), np.radians(R20D01_R_deg), np.radians(R38G02_R_deg), np.radians(SplitT4c_R_deg)]\n", "colors = ['C2', 'C1', 'C0', 'C3']\n", "\n", "ax.axhline(y=0,xmin=0.0001,xmax=0.9999, color='0.05', linestyle=':', linewidth = 1)\n", "\n", "for j in range (len(data)):\n", " ax.plot([pos[j]-w/2, pos[j]+w/2],[stats.circmean(data[j]),stats.circmean(data[j])],\n", " color=colors[j], linestyle='-', linewidth = 4)\n", " ax.errorbar(pos[j],stats.circmean(data[j]),yerr=[stats.circstd(data[j])], color='k',zorder = 0) \n", " for i in range(len(data[j])):\n", " ax.scatter(pos[j] + np.random.random(1) * w/2 - w/4, data[j][i], s=25, color=colors[j], alpha=0.5, zorder = 2) \n", " \n", "ax.spines[\"top\"].set_visible(False)\n", "ax.spines[\"right\"].set_visible(False)\n", "ax.spines['left'].set_position(('outward', 7))\n", "ax.spines['bottom'].set_position(('outward', 7))\n", "ax.set_xticks(pos)\n", "ax.set_xticklabels(('VS',' LPi3-4', ' LPi4-3', 'T4c'), rotation=45, size=10)\n", "ax.set_yticks(np.radians(rotations[::2]))\n", "ax.set_yticklabels(rotations[::2])\n", "ax.set_ylabel(u'Preferred direction [°]', size=10)\n", "ax.tick_params(axis='both', which='major', labelsize=10)\n", "plt.ylim([0,np.pi*2])\n", "\n", "bbox_inches = 'tight'\n", "#plt.savefig('Fig1_R24E09_R20D01_R38G02_SplitT4c_R_deg.pdf',bbox_inches='tight', dpi=600, transparent=True)\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.15" } }, "nbformat": 4, "nbformat_minor": 2 }