analysis.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Gender representation at specialized astronomy conferences"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### *Kyle Willett (Univ. of Minnesota)*"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Recently, members of the astronomy community have started to track participation at conferences as a function of gender, particularly relating to question/answer sessions after oral presentations. This is intended to address some basic questions about behavior at conferences, such as:\n",
    "\n",
    "* How equal are the allotments of talks among men and women?\n",
    "* Are men and women asking questions at the same rate?\n",
    "* Does it matter if the speaker/session chair is a man or a woman?\n",
    "* Are women/men more likely to ask the ***first*** question in a session? Does this affect the gender balance of remaining questions?\n",
    "\n",
    "In the broader sense, this is intended to measure if our collective behavior at conferences --- in choosing speakers, engaging with them, and interacting with the audience --- is as inclusive as it should be. If not, then hopefully the community can come up with some ideas of best practices to follow and alleviate some of these trends. \n",
    "\n",
    "Data-gathering efforts were first led by James Davenport at the [223rd American Astronomical Society (AAS) meeting in 2014](http://www.ifweassume.com/2014/03/report-gender-in-aas-talks.html), and repeated at the [225th AAS meeting in 2015](http://nbviewer.ipython.org/github/jradavenport/aas225-gender/blob/master/analysis.ipynb) and the [National Astronomy Meeting (NAM) in the UK in 2014](http://arxiv.org/abs/1412.4571). These focus on question/answer sessions since these are public, quantifiable, and one of the main ways in which professional interaction takes place. These studies found that men ask disproportionally more questions than women, and that the gender of the session chair has a strong impact on the gender ratio of questioners. However, there aren't any data (yet) on whether smaller, more specialized meetings follow the same trends as AAS and NAM, which are much larger and thematically broader overall.  \n",
    "\n",
    "I tracked this data for a recent conference on \"**Unveiling the AGN-Galaxy Evolution Connection**\". It was held in Puerto Varas, Chile from [9-13 March 2015](http://www.astro-udec.cl/agn15/). There were a total of 200 people on the list of participants. A one-person team [(me)](https://twitter.com/kwwillett) kept track of the gender of speakers, chairs, and questioners for all of the 72 talks given at the conference (hereafter referred to as <span style=\"color:red\">AGN2015</span>). Here are the results. "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## What are the overall demographics of the conference?"
   ]
  },
  {
   "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",
      "                       M         F\n",
      "First question  0.614286  0.385714\n",
      "Questioners     0.641509  0.358491\n",
      "Attendees       0.564246  0.435754\n",
      "Chairs          0.785714  0.214286\n",
      "Speakers        0.549296  0.450704\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "//anaconda/lib/python2.7/site-packages/matplotlib/collections.py:590: FutureWarning: elementwise comparison failed; returning scalar instead, but in the future will perform elementwise comparison\n",
      "  if self._edgecolors == str('face'):\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaoAAAGSCAYAAACypoU2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XfcHFXd/vHPRQhSQ1WpIQlKiVICCigtAQsIKFWlCeoj\nog9gQbGgJhFFH/X3KCKoIF2aSAQMRRATwdAeWug9AYQI0luCJPn+/jhnk8lm9773rnvu5Hq/Xvva\n2alnzszutWdmdkcRgZmZWamWaHcBzMzMOuKgMjOzojmozMysaA4qMzMrmoPKzMyK5qAyM7OiOais\nZZLGSTo7dw+TNFdSS/uQpDMkHdu3JWyvzupE0rckndLf5eqK6ja2gUXS3ZK2b3c5+oKDqnCSDpF0\nl6TXJM2QdJKkFdtUnJ786C56OP2AFxE/iojP9dfyJE2XtGMHw0dLeqKud69voxx+cyVt2WDYeyRN\nlPS8pBck3SPpB5JWysMPydN+vW66f9Y+lCUdLOkWSS9JekLS/0gaVBl3FUl/kvRqrpP9KsMGS/qj\npGl5OTs0KPubkl7Jj5clDevdGuq6Rl/8IuLdEXFtu8rUlxxUBZN0FPBj4ChgCLA1sC5wtaTBvbys\nJVsZraeL6eH0farFOhhIgq7Xea9uI0kCPgXclZ+rw94PTAKuAzaIiJWBnYHZwKaVUZ8Hjpa0fKVf\nNVCXAb4ErApsBewEfK0y/ERgFvA24ADg15JGVoZfCxwI/IuFgzqA8yJihfwYEhHTW1t76zUR4UeB\nD1IwvQLsU9d/OeAZ4NPAmsDrwMqV4aOAfwOD8uvPAPeS3uxXAkMr484Fvgg8BDyS+x0PPA68BNwC\nbFsZfxxwdu4elqdfokn5RwG3AS8D5wPnAcdWhu8G3AG8AEwBNq4Mm076oLkz18GpwNuBK3K5rgZW\nqoz/UeCePK9JwIaVYZsDt+dy/AG4oFYOYDTwT+BoYAZwJrASMDHX8fPAn4G1KvObDPwIuCmX5eJa\n/Vfq5FPAY3k7fLtR/eXX2wLX53I/Dhyc+38kr8/LuXxHNanj9YC/Ac/mZf0eWDEPOxuYk/ePV4Cv\nNdiPZuZxXsnLWgMYm+vozNzvbmCLynRrAhfl+nkUOKKT/Xj7XL7R+XlwZdg/gOM7mf5gUpBdAnyv\n0v8JYPsm03wFuLSynm8A76gMPxP4UYPpFppn/TbrpKy17X9I3p7PAYcB7yXtyy8AJ9RN09H7c0PS\nvv4ccD+wb+5/KPCfvF6vAJdU3jc7Vsr9hw62Y9P3RYmPthfAjyYbJn2zfJMGQQCcAZybu68B/qsy\n7KfASbn7Y6QQ2oDUej4GmFIZdy7wF9KH81tyvwOAlfP4XyV9gC+Vh81709JBUAFLkT6ovwQMAvbO\nb6zv5+GjgKfzG7j2jXsa+UMsd18PvJX0wfg0KfQ2Bd6S1/l7edz1gVdJ36IHAV/P67xkpRxH5GF7\n5jd3rRyjcx3/CBgMLA2sksdbGlg+v4n/VFm3yaTwGAksC/yxQZ38NpdzE9I3+Q3y8LGVcdfNHxKf\nyGVbBdgkD5sBbJO7VwRGNdlH1svrPRhYDfg78PPK8GnkD64m0+8APFHXbxwpwHbO2+Y44IY8bAng\nVuA7uX6HA48AH+pgGacCp+TuJ4C9cvdypJZTw7CpTH8IKag2JX2Yr1SZV7Oguhg4rrKvvVY3/Kvk\nIKvr3yioxgIvksLibuCwDspa2/4n5X3vg3l/+1PePrV9efvO3p+5fp4gBfUSwGakLyMb5eGnk/fj\nRtu7k+3Y4fuixEfbC+BHkw2TDkXMaDLsx8BVufuzwDW5W6Rvctvm11cAn6lMtwTwGrBOfj0XGN1J\nOZ4nt3ZoPai2B56s6zeF+QHx6wZvsvuB7XL3NGC/yrA/AidWXh9ODg/gu8D5lWEiBckOuRz/rFvO\ndSwYVG+Qg7jJ+m8GPF95PYn8IZhfb5TnoUqdrFkZfhPw8Qb19y3goibLfIz0rXlIF/eZPYDbKq87\nC6rRNA6qqyqvRwKv5+6tgMfqxv8WcFqT+S9LanV+KL/+BXBx7l4719X6lfF/Qmp1vAock/sdAlyX\nuy8Afpy7GwYVqYXyOLBKfr0dde8j4HPApAbTNgqqjYDV8/Z9H/AU8Mkm61vb/mtU+j1LbglV9uUj\nO3l/DiV9gbm2bv6/Zf4XtDOoawGxcFA1244dvi9KfPgcVbmeBVZrcgXZGqRvVwATgPdJWp20A86N\niH/kYesCx+eT1C+QvhUCrFWZ1wIn0yV9TdK9kl7M06xI+jbYFWsCT9b1e6zSvS5wVK1ceTlr5+lq\nnq50z6x7PYvU2qkt6/HagEjvuidI67hGg3LUXzzw74j4T+2FpGUl/TafdH+J1EpZMZ9raTSPx5nf\noqn5V6X79UpZq9YhHTprZG/S4b/pkiZL2rrRSJLeLun8fGHBS6TDfas2mWdXVOv6dWDpvB+uC6xZ\nt92+RTr308iepBbrNfn1hcAuklYlBdJc0jYCICKOjnSe6k+kb/r1vgd8QVLD5Unag9Ry2CUins+9\nXyUdRq9akXTIrFMRcV9E/CuSG0iHxvfpZLKO9t2ZzN8fOnp/rgtsVVfX+5MOgUNrF700246N3p9P\nUPA5ZAdVuW4gfVPfu9ozn1Demfzmj4gXgKtI38D2J50LqnkcODQiVq48louIGyvjzNvhJW1HOnS2\nb0SslD80XqLrO/AMFgxDSG+8arl+WFeu5SPigg7m2awMT1bnnQNlHVKrqlE5hta9rn/DH0U6nLhl\nRKxIapmpbvlD67rfJH2x6IrHSYfuFhIRt0TEHqRDnxeTDj82chzpHNO7c1kPYsH3dGcfZo2GdzTN\nE8C0uu02JCJ2azL+wcAKwD8lzSCd2xoMHBARr5Fam3s3mK6+vlPBIh4gfTH7zkITSDsDJwO7RcQ9\nlUEPAktKekel36akw3jt1uz9eUMe9ve6YStExH/naVsJqmaavS96Ms8+5aAqVES8BIwHTpD04XwZ\n7TDSh9YTpG/PNeeSPhT2zt01vwG+XbvCSdKKkvbtYLErkM4bPCtpKUnfY+Fvo624Hpgt6chc7r1I\n56NqTgEOk7SlkuUk7Vp3VVerLgR2lbRjvhLyKFKL63rgRmCOpMMlLSnpY3XlaGR50rfelyStQjpH\nUSXgQEkbSVoW+D5wYW7JdcW5wAck7ZvLtqqkTXN9HSBpxYioXegwp4Oyvga8LGkt0peMqqdpEoaV\n4atKqm7jjr6U3Ay8IuloSctIGiTp3ZLeUz9iLs+OwK6kYKg9/of5V/8dDXxG0jdqrSRJa5MOoTWr\nz/GkC4lWqixrR+Ac0vmvW6oj50CcAHw/t5a3BXan8v6R9BZJS+eX1W4kfUzSynk/3RI4knRhR0/U\n6rij9+dEYH1JB+Z9YrCk90raMA9/GhjRzeXfQNffF23loCpYRPwU+DbwM1LL5kbSIbSdIuLNyqiX\nAu8gHYu/qzL9xaQPhvPzoaG7gA9XF1G3yCvz40HSFUQzqRxWY+HfQjX8MMll24t0fuE54OOkb9O1\n4beSzhP8inQO7CHSh1dHH/b1y408rwdI5/NOIB0O3RXYPSJm50N6e5HO471AulBkIunCjmbr8AvS\n5c7PksLuigbLPpt0jmAG6cT0kR3Mr34dauV+nHR47yhSHd1OuviCvD7T8jY7NJe7kfGkq7deIl2d\neFHd8n8EfCcfOvrqQoWJuJ/UAn9U6XdMa7DwNp63Tjk4dyOdt3uUVN8n0/jLzEHA7RHx14h4Jj+e\nBn4JbCxpZERMIYXZ9sAD+fDWFaTzgCfU11kuw3TgLNL5r5rvkL5kXaH5v3e6rDL8i6Rt+gzpysjD\nIuK+yvAHSIfG1iRdXPSapFqr+ROk/fNl5l8t2NEPolv5wlKrz6bvz4h4FfgQ8EnSUYMZpO25VJ7H\nqcDIvG0nNFlGs+3YyvuiKOr6F0GzgUvSTaSrIs/s5vSTSBdEnNa7JTNrn56+L/qaW1S2SJO0vaTV\n8yGOg4F3k1qNPZptLxTNrG366H3RZxa1X+Kb1duAdF5vOdJvfvbJh6B6wochbKDri/dFn/GhPzMz\nK5oP/ZmZWdF86K+bJLkpambWDRHRpfO8blH1QLv/VqSUx9ixY9tehlIe7a6LkvbNdtdFSQ/XxfxH\ndziozMysaA4qMzMrmoPKemz06NHtLkIxXBfzuS7mc130jC9P7yZJ4bqz0tT+5N37ppVKEtHFiyl8\n1Z+ZWQEWvJPMoqG3vjA5qMzMCrEotYR7M3h9jsrMzIrmoDIzs6I5qMzMrGgOKjMzK5ovpjAzK9R4\nje/T+Y+NsS2NN2zYMGbMmMFTTz3FqquuOq//qFGjmDp1KtOnT2fo0KEdzKFn3KIyM7MOSWLEiBGc\nd9558/rdddddzJw5s18uq3dQmZlZpw488EDOOuusea/PPPNMPvWpT/XLJfUOKjMz69TWW2/Nyy+/\nzP3338+cOXO44IILOPDAA/tl2T5HZWZmLTnooIM466yz2H777Rk5ciRrrbVWvyzXQdUDfX2i06y7\nerpvjj1nXO8UxBYZkjjooIPYbrvtmDZtWmuH/c7tnfNXPvRnZmYtGTp0KCNGjOCKK65gr7326rfl\nukVlZmYtO/XUU3nxxRdZZpllmD17dr8s00FlZmYtGzFixAKv++PydAeVmVmhWv1Bbl+bNm1aw/5L\nLrkkc+bM6fPl+xyVmZkVzUFlZmZFc1CZmVnRHFRmZlY0B5WZmRXNQWVmZkVzUJmZWdHaFlSSjpF0\nt6Spkm6XtGUvzvsQSSf01vzMzKx92vKDX0nvA3YFRkXEm5JWAd7Si4vo9g1SJC0REXN7sSxmZtYD\n7fpnitWBZyPiTYCIeB5A0nTgAmAXYCawf0Q8IumtwK+B2r2OvxwR1+dW2C+ApfP4n46IB6sLkrQr\ncAywO7AFMI4Uio/k8V/Lyz0f+CDwE0lvBz4PzAbujYj9+qISzMw61Ev/Pt7U/q19px82bBjPPPMM\ngwYNAtLfJj344IOsvvrqfVm6edoVVFcB35P0APBX4IKIuJbUEnoxIjaRdBAphHYHjgd+HhFTJA0F\nrgRGAvcB20XEHEkfAI4D9gEEIGlP4Cuk4BtMCqydImKmpG8AXwWOzct9NiK2yNM9CQzLrb0h/VEh\nZmalksTEiRPZcccd27L8tgRVbsVsAWwHjAEukPStPPi8/Hw+8PPc/QFgo8qfH64gaVlgJeAsSe8g\nhU11fXYE3gN8MCJelbQbKdyuz/NZCri+Mv4Fle47gXMlXQxc3Gw9JjFpXvcwhjGc4S2svZnZ4mPy\nvTD5vp7No21/SpvPA/0d+Luku4BDGo2WnwVsFRH/qQ6UdBJwTUTsKWldYHJlukeA4cAGwK25/9UR\nsX+TIr1W6d4V2J7UmjtG0sYRsdA/L45hTIfraGa2qOj0JolNjB6ZHjXjJ3R9Hm256k/S+pLeWek1\nCpieuz9Rea61eK4CjqxMv2nuHAI8lbs/XV0E8BjpMOBZkkYCNwHbSFovz2O5ujLU5i1gaERMBr4J\nrAgs1/W1NDNbNEQEe+yxByuvvDIrr7xyv940EdrXoloeOEHSSqQLFh4iXbywG7CypKnALKB2EcOR\nwIm5/5KkltgXgZ8AZ0r6DnAZ81tgAUREPCDpAODCPO9DgPMk1a4wPCYvu2oQcLakFUmBd3xEvNyb\nK29mNpBI4pJLLlnszlHdBmxT3z+fO/pJRHyzbvzngE82mM+NpEN7Nd/N/c8EzszddwDvysOnAQv9\nXisihle6Z5POnZmZWQFK+2eKbv/+yczMFk1F3eE3IkZ0PpaZmS1OigoqMzOraPEHuX2t2a3o+0tp\nh/7MzMwW4KAyM7OiOajMzKxoDiozMyuag8rMzIrmoDIzs6I5qMzMrGj+HZWZWSEqtzKyCgeVmVkB\n4px2l6Bc6u49RhZ3ksJ1Z6WpfSP3vmmlkkREdKnp6HNUZmZWNAeVmZkVzUFlZmZFc1CZmVnRHFRm\nZlY0B5WZmRXNQWVmZkVzUJmZWdEcVGZmVjQHlZmZFc1BZWZmRXNQmZlZ0RxUZmZWNAeVmZkVzUFl\nZmZFc1CZmVnRHFRmZlY0B5WZmRXNQWVmZkVzUJmZWdEcVGZmVjQHlZmZFc1BZWZmRXNQmZlZ0RxU\nZmZWNAeVmZkVzUFlZmZFc1CZmVnRHFRmZlY0B5WZmRXNQWVmZkVzUJmZWdEcVGZmVrQl212AgWy8\nxre7CGYNed9sn7HnjGt3ERY5blGZmVnRHFRmZlY0B5WZmRXNQWVmZkVzUJmZWdEcVGZmVjQHlZmZ\nFc1BZWZmRXNQmZlZ0RxUZmZWNAeVmZkVbcAFlaTVJZ0v6WFJt0i6TNLnJP25i/MZL2mnviqnmZn1\njgH1p7SSBPwJOD0iPpn7bQJ8tKvzioixTZaxRETM7VFBzcys1wy0FtUY4D8RcXKtR0TcCVwHLC/p\nQkn3Sfp9bbik70q6WdJdkn5b6X+GpL1z93RJP5Z0K7CvpCMl3SNpqqTz+m/1zMys3oBqUQHvBm5t\n0F/AKGAkMAOYImmbiJgC/CoijgWQdJak3SJiIhD5QX5+NiK2yOM9CQyLiDclDenbVTIzs44MtKCK\nDobdHBFPAUi6AxgGTAF2lPR1YFlgFeBuYGKD6S+odN8JnCvpYuDiZgucxKR53cMYxnCGt7YWZmaL\nicn3wuT7ejaPgRZU9wD7NBn2RqV7DjBI0tLAicAWEfGkpLHA0k2mf63SvSuwPbA7cIykjSNiTv0E\nYxjT1fKbmS1WRo9Mj5rxE7o+jwF1jioi/ga8RdLnav3yxRTbNZmkFkrPSVoe2LezZeQLNoZGxGTg\nm8CKwHI9KbeZmXXfQGtRAewJ/ELSN4BZwDTgEhocFoyIFyWdQjrc9y/gphbmPwg4W9KKpHNfx0fE\ny71VeDMz65oBF1QRMQP4RINBv6uMc0Sl+7vAdxvM59OV7uGV7tk0b6GZmVk/G1CH/szMbPHjoDIz\ns6I5qMzMrGgOKjMzK5qDyszMiuagMjOzojmozMysaA4qMzMrmoPKzMyK5qAyM7OiOajMzKxoDioz\nMyuaIjq6F6E1Iylcd1aadJca8L5ppZJERKgr07hFZWZmRXNQmZlZ0RxUZmZWNAeVmZkVzUFlZmZF\nc1CZmVnRHFRmZlY0B5WZmRXNQWVmZkVzUJmZWdEcVGZmVjQHlZmZFc1BZWZmRXNQmZlZ0RxUZmZW\nNAeVmZkVzUFlZmZFc1CZmVnRHFRmZlY0B5WZmRXNQWVmZkVzUJmZWdEcVGZmVjQHlZmZFc1BZWZm\nRXNQmZlZ0RxUZmZWNAeVmZkVzUFlZmZFc1CZmVnRHFRmZlY0B5WZmRXNQWVmZkVbst0FGMjGa3y7\ni2DWUF/sm2PPGdfr8zRrhVtUZmZWNAeVmZkVzUFlZmZFc1CZmVnRHFRmZlY0B5WZmRXNQWVmZkVz\nUJmZWdEcVGZmVjQHlZmZFc1BZWZmReu3oJK0h6S5kjbIrzeVtEtl+A6S3teLyxsn6ajemp+ZmbVH\nf7ao9gMm5meAUcBHKsPHAO/vxeVFL87LzMzapF+CStLywFbA4cAnJA0Gvp+7b5d0NPB54Cv59TaS\n3irpj5Juzo/353mNk3SapEmSHpF0RGU5x0h6QNJ1wAaV/utJukLSLZKurbTqmi1jh1yO2yXdlstv\nZmZt0F+3+fgYcGVEPC7p38DGwHeBLSLiSABJywCvRMT/5tfnAj+PiCmShgJXAiPz/NYntcCGAA9I\nOgnYDPgEsCkwGLgNuCWPfzLw+Yh4WNJWwEnATsDxTZZxFPDFiLhB0rLAG31WM2Zm1qH+Cqr9gJ/n\n7gvz67sB1Y1Xff0BYCNpXq8VJC1HOqR3WUS8CTwn6RlgdWA7YEJEzAJmSboUIE/zfuDCyryW6mQZ\nU4CfSzonz/PJRis1iUnzuocxjOEMb6EqzMwWH5Pvhcn39WwefR5UklYhtX7eLSmAQaSwuaezSYGt\nIuI/dfMDqPabQ1qPYMGgq3UvAbwQEaNaXQbwP5ImArsCUyR9OCIeqJ94DGM6WQUzs8Xb6JHpUTN+\nQtfn0R/nqPYBzoqIYRExPCKGAtOAocAKlfFeqXt9FXBk7YWkTTtYRgDXAntIWlrSCsBuABHxCjBN\n0j55PpK0SZNlbJaf14uIeyLiJ8D/UTnfZWZm/as/guqTwJ/q+l1EOlw3Ml+wsC/wZ2DP2sUUpAB5\nj6Spku4hXWxRs9AVfRFxO3ABMBW4HLi5MvgA4LOS7iAdcvxo7l+/jENz/y9JukvSVFLr7YrurryZ\nmfVMnx/6i4gdG/Q7ocno9a2mTzaYdnzd640r3ccBxzWYZjqwS4P+zzVZxpH1/czMrD38zxRmZlY0\nB5WZmRXNQWVmZkVzUJmZWdEcVGZmVjQHlZmZFc1BZWZmRXNQmZlZ0RxUZmZWNAeVmZkVzUFlZmZF\nc1CZmVnRFLHQH5FbCySF685KU7sJqPdNK5UkIqL+prkdcovKzMyK5qAyM7OiOajMzKxoDiozMyua\ng8rMzIrmoDIzs6I5qMzMrGgOKjMzK5qDyszMiuagMjOzojmozMysaA4qMzMrmoPKzMyK5qAyM7Oi\nOajMzKxoDiozMyuag8rMzIrmoDIzs6I5qMzMrGgOKjMzK5qDyszMiuagMjOzojmozMysaA4qMzMr\nmoPKzMyK5qAyM7OiOajMzKxoDiozMyuag8rMzIrmoDIzs6I5qMzMrGgOKjMzK5qDyszMirZkuwsw\nkI3X+HYXwayhgbxvjj1nXLuLYIVxi8rMzIrmoDIzs6I5qMzMrGgOKjMzK5qDyszMiuagMjOzojmo\nzMysaA4qMzMrmoPKzMyK5qAyM7OiOajMzKxoPQ4qSWtLukTSg5IelvQLSYN7o3B5/h+TtFHl9XhJ\nO/XW/M3MrGw9CipJAiYAEyJifWB9YHngh71Qtpo9gZG1FxExNiKu6cX5L0CS/6jXzKwgPW1R7QjM\njIgzASJiLvAV4DOSviDphNqIkiZK2iF3f0jS9ZJulfQHScvl/j+WdI+kqZJ+Kul9wO7ATyXdJmmE\npDMk7Z3H3yn3v1PSqZKWyv2nSxqX53+npA1y/+UknSbppjzdR3P/QyRdKuka4GpJq0u6VtLtku6S\ntG0P68nMzLqpp0H1LuDWao+IeAV4HBhUN24AIWk14Bhgp4jYIk//VUmrAHtExLsiYlPg2Ii4AbgU\n+FpEbB4Rj1bmszRwOvDxiNiEdMuSL1SW9e88/18DX8v9jwGuiYitSCH7U0nL5mGjgL0jYgxwAHBl\nRIwCNgHu6GE9mZlZN/X0MFd0MKzReSoBW5MO5V2fjhyyFHA98BIwS9KpwMT8qE5XP58NgGkR8XDu\ndybw38Dx+fWE/HwbsFfu/hCwu6RacL0FGJrX4+qIeDH3vxk4LZ9ruzgipjZawUlMmtc9jGEMZ3ij\n0czMFluT74XJ9/VsHj0NqnuBfao9JA0B1gH+DbyjMmjpSvfVEbF//cwkbQnslOd5eO6GxoFY3091\n/d7Iz3NYcD33ioiH6pa7FfDavBlHXCdpO2A34AxJ/xsRZ9cXYAxjGhTLzMxqRo9Mj5rxE5qP20yP\nDv3lixqWlXQQgKRBwP8DzgWmAZspWQfYkhQkNwLbSFovT7OcpHfm81QrRcQVwFeBTfNiXgGG1C8a\neAAYVpsPcBDw906K/BfgyNoLSaNqndWRJA0lHTr8HfA70mFBMzNrg974HdWewD6SHgSeJYXK1yJi\nCims7iUdjrsVICKeBQ4BzpM0lXTYbwNgBeDPud91pIsyAM4Hvp4vjBhRW2hEvAF8GrhQ0p3AbOA3\ntcGV8kXl9bHA4HyBxd3A+AbjAIwG7pB0G/Bx5h9ONDOzfqaIjk4zdXFm6Sq9U4B9I6KHRyXLJinG\nMa7dxTBbQG2fHMj75thzxrW7CNaHdABERP11Bx3q1d8M5av03t2b8zQzs8Wb/0LJzMyK5qAyM7Oi\nOajMzKxoDiozMyuag8rMzIrmoDIzs6I5qMzMrGgOKjMzK5qDyszMiuagMjOzojmozMysaL36p7SL\nE0nhurPS5JuR4n3TSiWpy39K6xaVmZkVzUFlZmZFc1CZmVnRHFRmZlY0B5WZmRXNQWVmZkVzUJmZ\nWdEcVGZmVjQHlZmZFc1BZWZmRXNQmZlZ0RxUZmZWNAeVmZkVzUFlZmZFc1CZmVnRHFRmZlY0B5WZ\nmRXNQWVmZkVzUJmZWdEcVGZmVjQHlZmZFc1BZWZmRXNQmZlZ0RxUZmZWNAeVmZkVzUFlZmZFc1CZ\nmVnRHFRmZlY0B5WZmRXNQWVmZkVzUJmZWdEcVGZmVjQHlZmZFW3JdhdgIBuv8e0ugllDJe6bY88Z\n1+4i2ADlFpWZmRXNQWVmZkVzUJmZWdEcVGZmVjQHlZmZFc1BZWZmRXNQmZlZ0RxUZmZWNAeVmZkV\nzUFlZmZFc1CZmVnROg0qSXMk3Z4ft0laV9KUrixE0pclLdP9YnafpIMlrVF5fYqkjdpRFjMz67pW\n/pT29YgYVddvm/qRJC0ZEbObzONLwNnAzC6WrzccAtwNzACIiM+1oQxmZtZN3Tr0J+nV/Dxa0nWS\nLgHulrSspMsk3SHpLkkfl3QEsCYwSdI1Dea1s6T7JN0q6ZeS/pz7j5N0VGW8uyUNzd0HSropt/J+\nI2kJSYMknZGXe2duxe0NvAc4J7cGl5Y0WdIWeT775XHvkvTj6vpJ+kFejxskva079WRmZj3XSlAt\nUzn0d1HuF5Xho4AjI2JDYBfgyYjYLCI2Bq6IiBOAp4DREbFTdcaSlgZOBnaLiC2At1fmXV3GvNf5\nsN3Hgffnlt4c4ABgU2DNiNg4IjYBTouIi4BbgP0jYvOImJXnE5LWBH4MjAE2A94r6WN5WcsCN0TE\nZsC1gFthZmZt0sqhv5kNDv1V3RwRj+XuO4Gf5dbJxIj4Ryfz3hCYFhGP5Ne/Bw7tYHwBOwFbALdI\nAlgGeBr4MzBC0i+By4Cr6qarn897gckR8RyApHOA7YFLgP9ExGV53FuBDzYqzCQmzesexjCGM7yj\ndTUzW+xMvhcm39ezefTGjRNfq3VExEOSRgG7Aj+QdE1EHNvBtPWtpmqgzGbBFt/Sle4zI+Lb9TOT\ntAmwM3AEahYUAAAPN0lEQVQYqdX12SbLabbsWr83K/3n0qSexjCmUW8zM8tGj0yPmvETuj6PXr08\nPV9dNysizgF+RjosCPAKMKTBJA8AwySNyK/3Y35YTAc2z/PdHBieh10D7CPprXnYKpKGSloVWDIi\nJgDf7WTZAdwM7CBpVUmDgE8Cf+/uupuZWd9opUXVWWuk2r0x8FNJc0mtksNy/5OBKyU9WT1PFRGz\nJB0KXCbpdeA6YL08+CLgU5LuBm4ihRoRcZ+k7wBXSVoiL+eLwCzg9NwP4Jv5+QzgN3n+768s+1+S\nvglMIrWmJkbEn5usX6M6MDOzfqCIcj6DJe0AfC0idm93WTojKcYxrt3FMFtAbZ8scd8ce864dhfB\nCqADICLqrxvoUIn/TFFOcpqZWdv1xsUUvSYi/o7PE5mZWUWJLSozM7N5HFRmZlY0B5WZmRXNQWVm\nZkVzUJmZWdEcVGZmVjQHlZmZFc1BZWZmRXNQmZlZ0RxUZmZWtKL+lHYgkRSuOytNvpko3jetVJIW\niT+lNTMzm8dBZWZmRXNQmZlZ0RxUZmZWNAeVmZkVzUFlZmZFc1CZmVnRHFRmZlY0B5WZmRXNQWVm\nZkVzUJmZWdEcVGZmVjQHlfXY5MmT212EYrgu5nNdzOe66BkHlfWY34TzuS7mc13M57roGQeVmZkV\nzUFlZmZF840Tu0mSK87MrBu6euNEB5WZmRXNh/7MzKxoDiozMyuag6oTknaWdL+khyR9o8k4v8zD\np0oa1d9l7C+d1YWkA3Id3ClpiqRN2lHOvtbKPpHHe6+k2ZL26s/y9acW3x+jJd0u6W5Jk/u5iP2m\nhffHapKulHRHrotD2lDMfiHpNElPS7qrg3Fa/9yMCD+aPIBBwMPAMGAwcAewUd04HwEuz91bATe2\nu9xtrIv3ASvm7p0XxbpopR4q4/0NmAjs3e5yt3GfWAm4B1g7v16t3eVuY12MA35UqwfgOWDJdpe9\nj+pjO2AUcFeT4V363HSLqmNbAg9HxPSIeBM4H/hY3TgfBc4EiIibgJUkvb1/i9kvOq2LiLghIl7K\nL28C1u7nMvaHVvYJgCOAPwL/7s/C9bNW6mJ/4KKI+CdARDzbz2XsL63UxQxgSO4eAjwXEbP7sYz9\nJiKuA17oYJQufW46qDq2FvBE5fU/c7/OxlkUP6BbqYuqzwKX92mJ2qPTepC0FulD6te516J6aW0r\n+8Q7gVUkTZJ0i6SD+q10/auVujgFeJekp4CpwJf6qWwl6tLn5pJ9XpyBrdUPmPrfBCyKH0wtr5Ok\nMcBngG36rjht00o9/AL4ZkSEJLHw/rGoaKUuBgObAzsBywI3SLoxIh7q05L1v1bq4tvAHRExWtJ6\nwNWSNo2IV/q4bKVq+XPTQdWxJ4F1Kq/XISV/R+OsnfstalqpC/IFFKcAO0dER03/gaqVetgCOD9l\nFKsBu0h6MyIu7Z8i9ptW6uIJ4NmImAnMlHQtsCmwqAVVK3XxfuCHABHxiKRpwAbALf1SwrJ06XPT\nh/46dgvwTknDJC0FfAKo/7C5FPgUgKStgRcj4un+LWa/6LQuJA0FJgAHRsTDbShjf+i0HiJiREQM\nj4jhpPNUX1gEQwpae39cAmwraZCkZUknzu/t53L2h1bq4n7gAwD5fMwGwKP9WspydOlz0y2qDkTE\nbEmHA38hXdVzakTcJ+nzefhvI+JySR+R9DDwGvDpNha5z7RSF8D3gJWBX+fWxJsRsWW7ytwXWqyH\nxUKL74/7JV0J3AnMBU6JiEUuqFrcL44DTpc0ldRIODoinm9bofuQpPOAHYDVJD0BjCUdBu7W56b/\nQsnMzIrmQ39mZlY0B5WZmRXNQWVmZkVzUJmZWdEcVGZmVjQHlZmZFc1BZYscSeMkzZW0Q7vL0hOS\nhuRbIUzPtwuZuyjcOqU3tk+eflJvlsvK5aCyBeQPgI4eBxdQxkM6KUtUHgPZT4DDSX9g+kPSbSKK\n/9eTftw+A2L7Sjoj18fQdpdloPI/U1gjAYxvMuz2/ixIJ5p9UP0KOI8F/515INoNeCAiGt1GZCDo\ny+2zIfB6D6bvbwMiVEvloLKGIuL77S5DCxr+K3lEPEe6Kd1AtyYwud2F6IE+2z4R8WBPpm+DRflf\n9Pteu+8E6UdZD9L/sc3pZJxhebzTgfWBC4BngDnADnmcLYDjSYetngNmAg8CPwNW6mDenwCuAZ7P\n00wDzgW2yMMn52U3egzN44zLr7dvMP+dgCvz/GcBDwA/AoY0GLe2rEGkWzQ8lKd5HPgxMLiLdbsG\ncCIwHXgj19lFwOZNllv/mNTCMqbnOhtCark8mevxHuCIJtMcksvxKKmV8hLwD+CAJuPXyjeY9P+O\nD+R6Ob0Xts+GwGl5PWaRDnVeCxzWYD+dVNevNt8dgINJrf/X8zxOBd7eYHld2k9zXc3N8x+T1/fl\nXGcTgQ0blLPRY1plnBHAyaQ7BL+ey3En6X5mq7T7M6GEh1tU1hPrATeSPqjOBpYhvWEBPgfsQXoj\nX0U6H/oe4Kuk215sFRGv1maU79t0Oukflf/N/LvjrgOMJv3z9K15nBdINya8mHTL75qX6ED+g9Bf\nA68AF5KCYgzwDWB3SdvE/DsUV50HbEu6EeTLwK7A0cDbSPfd6pSk4aQP/zVIQXwOMBTYF9hV0t4R\ncVke/XRgEumPPKcDZ+T+01tYVABL5WUMIYX8W4C9geMlbRARh9dNcxJwN2lbzSDdmuQjwNl5/O81\nWdYE0ja9PHc/k+fR3e2zK2m7LAVcQaqjlYDNgK8Dv2mwro18BfgQ6S67l5Nui/5pYHTe76p3Ge7S\nflqxW17Hy0n71LtIdfZeSSMjtRohHULfg3Rrk18AL+b+L+Z1XgP4P2AF4LK8/kuTwutA4ATSl6rF\nW7uT0o+yHsz/xjeW9A21+jg4jzOsMt4PmsxnKPlPj+v6fyZPd3Rd/0Nz/xuBFeqGLQGsXnl9SB73\nU02WPY66b+zAuqRWzIvA+nXjn5jH/21d/8m5//9R+XZNugHgQ8BsGnxLb1Kmv+R5fauu//uAN4Fn\ngeUabIu/dXH7Tc/TXUulxUf6V/uH87Dt6qYZ3mA+g4G/Av8B1mxSL3fQ4Bt/N7fPaqQgm1Vfvjy8\nvgwL1U1lvrOATeuG/W8e9rse7qe1dfsPMKZu2HF52Nfr+p9BpUVZN+yIPGyh1i7pi9/SXdn+i+rD\nV/1ZM2NJh3Wqj/qruP5Fk4suIuLxyO+2OqeTWjQfqut/BOkb8uej7o6nETE3Iv7V5TVY0IGkD99f\nxcLnN44BXgUOzPcSqveNiKh9EyYiXid921+CdOioQ5LWBj4IPEa6km+eiLiB1GJbBdir5bXp3Lci\n4s3Kcl4Ajs0vF7ilQkRMq584T3sS6Tz2Tk2W8d3ovdtUHExqVfw6Iq5rUJ6nujCvsyNial2/caTW\n8H7VbdyN/bTm/Iiovzz+5Pz83i6UtbbsWQsNiJgZEQv1Xxw5qKyRiIglGjx2rBtvavXDsErSYEmH\nS/qHpOdrvwMincdaAVirMu5ypEMnTzf4gOktm+fnv9UPyCF0O+kb7Ib1g2l8B9ba3VtXbmHZo/Lz\ndRExp8HwWpk2a2FerXgTuL5B/8mNliNpqKQTJd0v6bXaTxFIh18hXdRRL4Cbe6m8AFvn5yt6YV5/\nr+8RES+TWoALbOOu7Kd1erpP1FxK+pJ0oqQ/SjpU0sguTL9Y8Dkq64mOWjkXkI7NPwL8KY/7BunK\npy+TzpvUrJSfm96KuhesmJ9nNBk+o268efKHXL3Z+XlQLyy7Vo8rNRneVc82aSXUfoM1bx0ljSAF\nzkqkw4VXkg7BzQGGk1o6b6GB6N07WffmPtCsXLV6rm7jruynVS/W94h080RobZ+oTfO4pC1JLb6d\nya3qfLPBn0XECa3Oa1HmoLKeaHgyW9J7SG/+q4FdImJuZZhIFy9U1d70zb699obaifw1gPsaDF+j\nbry+WPbqTYb39rJXk6QGYVVbfnU5XyUddjwkIs6qjixpPxY+3NtXavvA2qQrFHvi7U36L7D+3dhP\n+0RE3A98UtIg0kUXHyAdCj9e0msRcVp/lKNkPvRnfeEd+fnS6ps/24p0VdM8EfEa6aqz1SW1cvir\ndvis5W+uwG35eXT9AEm1K8tm0jjEeqq27G3zh1G9MXXj9dRgYJsG/Ufn5+qPtt9B+sJxUYPxu/sX\nR93ZPjfk5126ucyq0fU9JK3Iwtu4S/tpD7RUHxExJyJui4ifAPvl3gP1x969ykFlfaF2cn5Mtaek\nt5GusGvkl/n5t5KG1E23hKRqa6R26e+6XSjT70nnbo6QtF7dsGNJ5yN+3+ycW09ExJOkb+3DSYeT\n5pG0FbA/6RLkP/XiYn9UvWhA0irAd0ihdHplvGmkw1z12+rDwH91c9nd2T5nki52+IKk7eoH5gtS\nWnVQgy8840iX659X2cbd2U+7o2l9SNo8h2i92v4+kP59o8/40J/1hf8DpgB7SZqSu99OOgZ/P/AU\ndb/Sj4jf5Q+og4CHJF1K+h3VmqQPklOB2r9lXE96A39Z0qrMPyfxyybnk4iIxyR9mfQBdJukP5Au\nCd+BdCL/Phof6umtfxM4jFQPP5X0IdJvwtYh/Y5qNvDp3LLsDTNI51buzvW4FLAP6cPvxIj4R2Xc\nk0hXAV4o6Y952ncDHwb+QPoBdiMd1Ut3ts9zkvYnXcAxSdIVwF2kcNmEdEhwRIdrPd/lwJS8jf9F\n+g3cNqRg+mZlvC7vp930V+BrwCmSJpCuJnwhIk4k/W7wUEn/IP3g+gXS7xN3J10J+IteWP7A1+7r\n4/0o60HX/pnitA7GWZkUCtNIh1seAn5AuupqGvBok+n2J12d9mKe7hHSj4k3qxvvw6QPxFdqZWb+\nPx+Mza8b/fPBB0m/aar9M8WDpH+ZaPTPFJOa1QXp3M0cmvxWqMk0a5KCYTrz/5liAvlfN5psi+78\njupRUgvxV6Qr0WaRzvsc3mSa9zH/30BeJl1U8VFSiM8FvtdqvfTC9hlJal39M9fRjLy8/+qsbqj8\nPovW/5miS/tpZ9u92TYj/Qj53rwt5tbmC2yZ94k7SC2v1/M+eSowsq/f7wPloVxZZrYIkDQdmBsR\nrbY+FhmSxpF+7zc6Iq5tc3GsF/kclZmZFc1BZWZmRXNQmS1aFudj+cHivf6LLJ+jMjOzorlFZWZm\nRXNQmZlZ0RxUZmZWNAeVmZkVzUFlZmZFc1CZmVnR/j8UImpAGCBvKgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x107698d50>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%pylab inline\n",
    "\n",
    "'''\n",
    "Note: this notebook requires the following, fairly standard Python packages:\n",
    "\n",
    "    numpy\n",
    "    scipy\n",
    "    matplotlib\n",
    "    pandas\n",
    "    \n",
    "And one not-as-common package (also available through PyPI):\n",
    "\n",
    "    sexmachine\n",
    "'''\n",
    "\n",
    "from __future__ import division\n",
    "from matplotlib import pyplot as plt\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "\n",
    "gencolors =('purple','orange')\n",
    "urlbase = 'https://raw.githubusercontent.com/willettk/chile2015-gender/master'\n",
    "q = pd.read_csv('%s/question_data.csv' % urlbase)\n",
    "c = pd.read_csv('%s/chair_data.csv' % urlbase)\n",
    "\n",
    "fig = plt.figure(1,(6,6))\n",
    "\n",
    "# Speakers\n",
    "vc_speakers = q['speaker'].value_counts()\n",
    "\n",
    "# People asking questions of the speakers\n",
    "qa=list(q['questions'])\n",
    "qa.remove(' ')\n",
    "vc_questioners = pd.value_counts(list(''.join(qa)))\n",
    "\n",
    "# Chairs of the sessions\n",
    "vc_chairs = c['gender'].value_counts()\n",
    "\n",
    "# Attendees\n",
    "countrydata = pd.read_csv('%s/map/countries.csv' % urlbase)\n",
    "names = countrydata['name']\n",
    "firstnames = [x.split(' ')[0] for x in names]\n",
    "\n",
    "# Guess (based on first names) what gender the attendees are\n",
    "from sexmachine import detector as gender\n",
    "d = gender.Detector(case_sensitive=False)\n",
    "from collections import Counter\n",
    "genders = [d.get_gender(fn) for fn in firstnames]\n",
    "cg = Counter(genders)\n",
    "attendees = list('M'*(cg['male'] + cg['mostly_male'])+'F'*(cg['female'] + cg['mostly_female']))\n",
    "# Ignores users whose gender cannot be reliably determined from first name\n",
    "vc_attendees = pd.value_counts(attendees)\n",
    "\n",
    "# People who asked the first question of the speakers\n",
    "first = [x[1]['questions'][0] for x in q.iterrows()]\n",
    "first.remove(' ')\n",
    "vc_firstquestion = pd.value_counts(first)\n",
    "\n",
    "# Load everything into a single dataframe\n",
    "\n",
    "data = [vc_speakers,vc_chairs,vc_attendees,vc_questioners,vc_firstquestion][::-1]\n",
    "labels = ['Speakers','Chairs','Attendees','Questioners','First question'][::-1]\n",
    "\n",
    "normdata = [x/x.sum() for x in data]\n",
    "\n",
    "# Plot stacked bar chart\n",
    "\n",
    "ax1 = fig.add_subplot(111)\n",
    "df = pd.DataFrame(normdata,index=labels)\n",
    "dfplot = df.plot(kind='barh',stacked=True,ax=ax1,color=gencolors,legend=True)\n",
    "\n",
    "print df\n",
    "\n",
    "# Plot the 50-50 split for guidance\n",
    "\n",
    "ylims1 = ax1.get_ylim()\n",
    "ax1.vlines(0.5,ylims1[0],ylims1[1],color='k',linestyle='-')\n",
    "\n",
    "ax1.set_xlabel('Fraction of participants',fontsize=20)\n",
    "ax1.set_title('Overall demographics at the AGN2015 meeting');"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Analysis of overall demographics:"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "I tracked gender participation for five aspects of the question/answer sessions:\n",
    "\n",
    "***\n",
    "* speakers\n",
    "* chairs of the sessions\n",
    "* all attendees of the conference\n",
    "* people who asked questions of the speaker\n",
    "* people who asked the **first** question of the speaker for any given talk\n",
    "***\n",
    "\n",
    "Every one of these categories had more male than female participants. The gender ratio of the speakers closely matched that of the attendees as a whole (45% and 43%, respectively). \n",
    "\n",
    "Women were slightly less likely to ask questions than men; women constituted 43% of the attendees, but only 35% of the people asking questions. \n",
    "\n",
    "The gender of people who asked the first question in a given session is close to that of the overall gender of questioners (38% vs 35%). "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Overall demographics for <span style=\"color:red\">AGN2015</span> compared to other conferences"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "                       M         F\n",
      "First question  0.614286  0.385714\n",
      "Questioners     0.641509  0.358491\n",
      "Attendees       0.564246  0.435754\n",
      "Chairs          0.785714  0.214286\n",
      "Speakers        0.549296  0.450704\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABEgAAAGSCAYAAADw5vQ4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XmcXFWZ//HvNwQI+44IGpOwqAECAQEFWRJQQfZ9DYsK\nLqP8ZhQVBxnS4jAobogLKmEJExaRyCoKKAEkKAMhYd8TlEX2sIcleX5/nFOd25XqTnfnpqu67+f9\netWrq87dzrl1+6lbT517riNCAAAAAAAAVTao2RUAAAAAAABoNhIkAAAAAACg8kiQAAAAAACAyiNB\nAgAAAAAAKo8ECQAAAAAAqDwSJAAAAAAAoPJIkAC9ZPtI2zc3ux5YPGz/wfa4ZtcDAABgoLH9qu1h\nvVx2iu3P5ueH2v5TN5e7x/Z2nUzbwfY/e1OfrurXSmyfa/vkZtej1ZEgaSH5n+lF20s1mHaQ7b/b\nfs32M7b/ZvuLhenn2p5ne4tC2Xq25xVe/8D2Q7ZfsX1//Zc/25vavsP267Zvt71JYdpGtv9k+7ni\nOuvq/mYOdq/avr+MfYKe6W7gc/KY7Xu7mGeHfEx9o8G0z+Zj6BXb/7J9te3lO1lP048724Pz/86W\nhbJDc/vqy+6XpIj4dESc39n+aVCPlvsgBLrCZ07vFeLjz+vK/2r7iE7m/UZd+bBcPq2ufHXbb9ue\n2cm217B9oe0nbc/O29yybp5DbD+e37/f216lMO0A21Pzfr+hizYenutHbEO/RIzrPXeSLMj75R3b\nazUoP7murBbjGn7fjIgVImJWL6sY+aGImBQRn+rWQhEbRcRNvdxmT7TXr8W0ar1aCgmSFuGUQd1S\n0rOS9qib9jVJP5H0PUnviYj3SPqCpG1sL1mY9UVJ3+1iM69J2i0iVpR0hKTTbX8sb2MpSZdLmihp\nZUnnSbq8sP63JV0kqbMTpZD0bznYrRARH+5OuyHZXqIJm91O0tKS1rD9kU7mOULSPZIOLxba3l7S\nf0s6KB9LH1Y6NjrT9OMuIt6VNDW3u2Y7Sfc3KLuxi7Z0ZpE+bDo7eQAWFz5zSvG6pMNsf6CuXvXx\noGEsLVjG9oaF14dIeqzBemqWl/R3SZtJWkVp311tezlJyus6U9Khkt4j6Q1Jvygs/4KkH0k6tbOG\n5YTKf+Z6czKNfocYV74cY/aVdJ+kw+om88W7/3ApKxnA564DtmH90OGSrpd0vlKQlSTZXklSm6Qv\nRsTkiHhdkiJiekQcFhHv5FlDKfiOciddxyJifEQ8lJ/fJulmSR/Nk3eQtEREnB4R70TEGUr/QGPz\n/A9FxDlKQbEzC/2Hs710/sVrw0LZGrbfyL+arW77Ktsv2X7B9k22u7PeBZYrTJtl+3jb9+ZfEs62\nvXRh+m62p+dlb7G9cWHa8bYfyb8O3Gt7ry7qcJrtm22vYHsl2xNsP2X7Cdsn1wKJ06U5t9j+ke3n\nJZ2Uf5W4Me+b52x3mnCwfYntp/O8N9oemcuPUTqx/kb+teHyLnbZEZIuVfrwPqJ+YuFD8AuShtre\nvDB5C0m3RsQMSYqIlyLi/Ih4rdGGWuG4y25Sx2TIx5VOjopl2+b56rtvHun0K+1p+Rh6zPbOedp/\n5+V+lvf7T3P5h2xfl4/HB2zv317h9EvLL50u43kt7wegL/GZswifOdlsSedKOqmL7XcVS2s6vAeS\nxil9qWpYj4iYGRE/iYhnIvmNpKUkbZBnOVTSFRHx1/z+nShpn1oCJSL+HBG/k/R0F237H0mnKyVT\ngP6IGLfoMa7evpJmSvq+Gpw7dqe+dXWfZ3tEfn6u7Z/nur7i1KNnRGHeT+Rzqdm2a/uyNu1I50ve\n87nVaXXbudz2v+fns2zvmJ8vk7f7olOP6i3qlptXV4f2XjK2V8l1fTYvf6XtdbrZ7vH5XP783Na7\nbK9v+1tOvZket/2Jwvzd/U7xktN3lq1tH2X7H3l99cn51W1fm7c9xfbQwrZ6dO5q+9O278vresIp\n+djvkSBpHYdLuljSbyV9yvYaufxjSr/0d/Vlt+YNSaco/brfJdvLKAWC2iUWG0q6q262Gbm8u/7H\n6cv9X516GSwgIt5S+mJ+cKH4AElTIuJ5SV+T9E9Jq0taU9K3IqI7GekFlqubfoikT0paV+kk8tuS\nZHu0pAmSjpa0qqRfSbrC8zP8j0j6eP51oE3S/9p+T3HFTn4jaSNJn4iIV5VOmt/O2xudt/25wmJb\nSno01/UUSSdL+mNErCxpHUk/7aKtV0taT9IakqZJmiRJEfHr/Px7+deGPRstbHtZpQ+52vF2kDv+\nYiJJ+0h6JiKmSrpSHT8I/6Z0jI63vY0LyaaFadZxl90kaZtcj9UlLSfpEqX3olb2oTyftOCvIVtK\nekDSakonBxMkKSJOUDopqv3Sc2z+InKdpP9Vep8OkvQL28VfgA6WdHJELC/plh60FygDnzmL9plT\nc4qkfW1v0Mn0rmJpzSSlOGynhHeth0i32N5UKUHySC4aqbQvJUkR8ZiktzQ/gbKw9W2p1DvlzO7W\nAWhBxLhyYlzREUr79ApJ69nerJfr6cyBksYr9Yx7RHm/5/OzS5V6ta2mdP68TSfruCCvR3nZVSR9\nQvN7OhfP7U6SNFzSCEmfUmpfV/umuKyVzgOH5sebkn7WnUZmuyklwleRdKfSOaMkra30neBXhXnP\n1cK/U8xQ+h5zodIxv1me/zClH/CWLdT7UEnfUTompit/j+jFuevUvA+Ozt+TNpT0lx7sg5ZFgqQF\n2P640pfiKyLiYaVs8qF58uqSno+I4jWPU3OW8I28bE0o/UMNdf51uwtnSpoeEdfm18tLerlunlck\nrdDNZnxTKcisLenXkq4sZl3rXKD0T1dzSC6TUgB4r6RhETE3Irr7xbGr5ULSzyLiyYh4SSng1j5I\njpH0q4j4v/xL3ESlE8mPSVJE/C4i/pWf/1bSw5K2Kqx7SaWgu7Kk3SNiTk6g7CLpPyLizYh4Tqkr\nZ7HNT0XEzyNiXkTMyfUfZnudiHg7n0w3FBHnRsTr+VeONkmb2C6+TwvL4O8j6ZW8j2qBbNe6eY5Q\nSh4o/z3I9uC8/b/mdWwm6SpJz9v+obvX1a6Zx91tkpa1PUqpx8fNEfGmpJmFslkR8UQnyz8eERPy\nicVESe+1vWZhenG/7yZpZkScl9/j6ZImS9q/MM9lEXGr1H6CA/QJPnNK+cyRJEXEM0pt+04ns3Qa\nSwuekPSg0kn84UrxpVtsr6j0C/n4nJyXFmHfOl3y+XNJX16EL1FAUxHjyotxNbmXwQ6SLsmx5k/q\n/LLB3ghJkyPi9oiYq/SlfdM87dOS7sk9fuZGxE8k/auT9fxVUtjeNr/eT9LU2rl8nf0l/XdEzM7n\nfqdr4efQlqSIeDEifh8RcyL1oD5FUlc/0tW7KSKuy239nVLi59T8+mKl7wQrdvM7Re18M5SSI2tL\n+k7uuXSd0jGwXmH+q3IPw7clnSDpY7bfp56fu9a+v2xoe8WIeDki7uzBPmhZJEhawxGSri2c3Fyi\n+b8yvaDUFar9vYqIrSNilTytw3uYD/aT86PhyY1T17ORShnmmlclrVg360pKwXyhIuK22pf2nGS4\nRSmgNTJF6Yvqlk7XiG4i6fd52mlKWeNrbT9q+5vd2X43lisONPUPpeAhSR+Q9LX8wfiS7ZckvU/p\nw6Q2SN2dhWkbKQWxmvUk7a4UiN4trHNJSU8XljtTKRvbqD6S9A2loHub0wjbRzVqpO1Btk916kL3\nslJXRyl94HfXEUoBTzkQX6aO3U/fr/whmIv+KGmICkmUiPhjROyRj8M9JR2pjtnsRnVv6nGXA/lt\nSpfUbKvU60NKH6a1sq7GH2n/cI2IN/LT4sC0xf+3D0jaqu64OkRpPIDavKWMlA70Ap85i/6ZU/R9\npV+oRxULuxNLa81RSoocpXTSe76617V+GaVeKVMj4nuFSa8p7cuilZT2+cJ8SdJdkS4XaN9UN5YD\nWgkxrtwYJ6VL/+6JfEmR0j49xPPH0XtX6dy3aElJ84rJqIV4pvD8Tc0/x1pbKZFc1PAcKicJLtL8\nH0IPUe4h0cDaWvD7QbfYXtb2r5wu2XlZ6fxxJbvbly89W3j+plLSLgqvpdT+7nynqN9vyomUYllt\nX4YK+zLSJWYvKu2L3py77qt0XM5yulznoxoASJA0WT7BOUDSWKdxJZ5W6g63ST7ZulWpR0OnY18U\nV5f/nqvUo2HfBttrU+pG9snoOGbEvZJG1c0+SvO7CpYmfyn/rVLwOljSlTH/GtDXIuK4iFhXaVCt\nr9oe2411NlpuTGGWoXXPn8zP/6GUPV6l8Fg+Ii52Gnjv15L+TdKq+cPzHnU8Wbxf0mckXVPoYv1P\npfdstcI6V4qIjQvLdfiQjXQt+TERsY6kzyt1aWv0S8GhuX07RsRKSr8uqFCnLn/xyxnisZKOKBxv\nB0j6tO1a4mecUmz4Q54+U+mkvlHXcEXEX5R6onTabbQVjrusNg5JMUFys1LWv338kV6o3+//kHRj\n3XG1QkT8Wy/XD5SCz5xyPnPq1v+C0i969YM59iSWTlY6yXw0Ou/F1s7p0sbLJP0jIj5fN/lepS9I\ntXnXVboE56G6+Rp9XoyVtHfh2Nha0g+dx1YCWh0xrvwYlx0uaf3CPv2J0o9ztYTvPyQNq1tmuMr5\nMegpSe+vvchJiPd3PrsulLRfPo/fUunynEae1oLfD4rekLRs4fV7NT9ufk3pssUt8/n49krHS3cS\nJD3pnded7xQ90WHfOd2BclWl70U9PneN1ONnL6WEzWVKx2G/R4Kk+fZSyrp+WOmEZpP8/GZJh0fE\nbKXLKH5he1+nAUAHOV1zvFxhPe3/kLknw0lK3fPmz2B/SylwfiLSpSZFUyTNtX2s04BPx0qap8K1\nZLaHKJ1k1QaFWjo/X8n2p2wPcbqd6qFKXzb/2EW7a90Bi90AZXtXpwFLrZRln5sfXepkuVrG2pK+\nZHsd26sqdSe7OE/7jaQv5Ky7bS+X17W80v4NSc9LGpR7dWxUv+2IuEjpusjrbY+IiKclXSvpR4X3\na113MshXrv/+OXkhpYH/olD/ouWVAuWLTtcKnlI3/Rmlayk7M05pHI0NNP9420Apm1zLth+hdA3o\nJoXHvkpJlFVt72n7QKcBqux0vfr2SmOTNGpbKx13Nyl9AXhfRNRumXeL0q+8m6r3CZJnlK71rLlK\n0ga2D7O9ZH5sYftDtWb1cjvAouIzp4TPnAZ+pHRpZvFa7S5jaXHh/GVmjBbSEy/XeUmlLtlvKPXe\nqzdJ0u62P54/J06WdGntC1N+P4co/So5KO/b2i+/RyqNxbSJUky8PbfhhIXVC2gRxLiSY5zTnXlG\nKI2xUtunG+Xt1C6zmSxpV6fBVJewvbbSeH8XdnczXUz7g9JlHHs7XaJ4rKS1Ops5XxryvKSzlMb3\n66zXzm8lfcv2yvkc/Ct106dLOjS3Z2d1HNR/eaWeGS/neN5osO7O2tTtc8DefKfohk87jSG4lNLn\nw60R8aTSGIfdPnfN0w+1vVJO0r2q3n1+tp6I4NHEh6RrJJ3WoHx/pYzpoPz6EKVB215X6pb1N6UT\nqSXz9HOULvOoLW9Jd0uaWyibp/TP/GrhcXxheu1k6I38d5PCtGF5+Xman3x4LE9bQ+nShVckvaQ0\naM+O3Wj7w0oBbHCh7N+VfmV7TSlrekJh2h+K9a1bV1fLzVT6ULs31+8cSUMK0z+V6/9S3ucXS1o+\nT/uuUpfL5yT9UNINkj6Tpx2hdA1hbT2fkzRLKQO9otJtFf+plPCYJumARsvlsu8pJSleVeoK+blO\n2rmcUob2ldyucfn9GJGnr6c02NNLStdy1i9/v9JgovXlX8/7YKv8/q/WYJ57lLpfb6s0MvxzuR4P\nSjqui/e5ZY47pQ+0t5WuoSyW3yvpibqyTt/rXFbc7x/N++FFST/JZRsoJUqeVTrOr5c0qtH/Kw8e\nffUQnzllfebsoNR7o1j29VzXw3NMWFgsHZbnH9Rgnp1q7W0wbfu8P16r27fbFOY5WNLjeZ7fS1q5\nMO3Iwr6tPc7uZFvtcZAHj/7wIMaVH+Mk/VJp7JH6ebbI7V85v94tt3O20vnw9yQt3UV9i+dR9fu7\nffv59aeUzrNmSzpDCz9H+3Ze/7515TMljc3Pl1G6U9FLSnH5uLptbp7LX1G6DHJSrY5KvUluyO/5\nA0pjGrbHc3URO5WSKRMLrzvEe0mD87rWzq+7/Z1C6XvA3Lrt/VPS1oX9/AulpMurSom8DxTm7fa5\nq1KS/Rqlc9+Xlf6ftm52DCjj4dxAYMCyPVPSZyNdCgIAAAAAwAK4xAYAAAAAAFQeCRIAAAAAAFB5\nXGIDAAAAAAAqjx4kAAAAAACg8gY3uwL9lW263gBoSRFRidsIE4cBtDJiMQA0V2/iMD1IFkGzb0HU\nV4+TTjqp6XUYyG1t1vHE+zowH1XT7P3NMTww2trM44n3dmA+qqbZ+5tjeGC0tVnHE+/rwHz0FgkS\nAAAAAABQeSRIAAAAAABA5ZEgwULtsMMOza5CnxmobZ329DQNahukj5/98fay+rY2mqfoxlk3aseJ\nO2qN09bQcqcsp/V+up4Om3yYXn3r1cVZ9VIM1PcV1VGlY7hKbZV63t4jLztSg9oG6bs3fbdD+ZRZ\nUzSobZBefPPFDuV7XLiHBrUN0vWPXb/AusZPGa9BbYO008SdFpj2y//7pQa1DdLGv9y4R/XrStXe\nWww8VTqGaevA1FVbj7zsSO1+4e4dyp557RktdfJS+sBPPtDwspVBbYM0+f7J3VpXf0GCBAtF0Oj/\nzpp2lrZYZwv97Ym/6YHnH5C0YFsbzVNz33P3aedJO2vT92yqKUdM0b1fuldn7namVh6yst6a+1Zf\nNaPXBur7iuqo0jFcpbZKPW+vbQ0ZPESnTT1Nz7/xfJfzPv3q07rmkWu01fu20lnTzlpwXbLWWn4t\n3fLPW/T47Mc7TJtw5wQNXWmorPLGGa3ae4uBp0rHMG0dmLpqq+0FYv55M87TBqttoNlzZutPj/6p\n29tptK7+ggRJyWwPqAf6vzffeVMX3nOh2nZo09jhYzVh2oQez3Pto9dq9WVX1w8/9UNtuOaGGrby\nMO00Yif97NM/0+rLrt5XTQEASBozfIyGrTxMJ994cpfznTv9XG3+3s31w0/+UJc9cNkCvUtCodWW\nXU27rr+rzpl+Tnv5Xc/cpQdfeFD7jdxPoeoNOAoAVRQRC8T8s+88W1/e8ss6cMMDNeHOBb9D9GRd\n/QUJksWg2SP2NnvkX7SW3933O6209Eraeb2ddczmx2jiXRP17rx3ezTPe5d/r557/TlNmTWlj2sP\nACiKCA3yIJ2646k6844z9dhLj3U639nTz9bRmx2trd+/tdZfbX2dP+P8hvN+ZvRndN6M89pfT5g2\nQQdueKBWWGqFxdIGAEDru/nxm/XEK0/o0I0P1dGbHa0rHrxioT0XBwISJMAAN+HOCfrM6M9Ikvb6\n0F6yrMsfuLxH8+y/4f46ZONDNPa8sVrrB2tpjwv30I9v/XElgiQAtBrL2mX9XbTN+7fRCX85oeE8\nNz5+o5557RkdvPHBkqSjNzu64a9/lrXzejvrnbnv6PrHrtdb776lSXdP0mdGf6bf/voHAFh0E+6c\noAM2PEArLL2CtlhnC224xoY6b/p5C1+wnyNBAgxgj7z4iG755y06atOjJEmDBw3WEZsc0eEkuTvz\nDPIgnb3n2Xriq0/oB5/8gYauNFSnTT1NH/rZh3Tfc/f1baMAoOJqiYvv7fQ9XXLvJZr29LQF5jlr\n2lk6aKODtOySy0qSxo0ap0defES3PXnbAusa5EE6YpMjdPadZ+uyBy7Tmsutqa3fv/XibwgAoCW9\n8tYr+t19v9PnNvtce1lnifbO9NcxSAY3uwIAFp+zpp2lufPmasRPR7SX1S6feuKVJ/S+Fd/XrXlq\n1l5hbR026jAdNuowfXfsd7XBGRvotKmn6Zw951+7DgDoG1uss4X2HbmvvnHdN3Tidie2l8+eM1uX\n3n+p3pn7ToexRebOm6uzpp2lLdfZcoF1HTX6KG38y401a/as9oQ5AKBaakmNC+6+QG+884a2P3f7\nDtPnxTxN/efU9iT6CkuvoJfnvLzAembPma2Vhqy0+Cu8GJAgAQaod+e9q/NmnKdTdzpVu22wW3t5\nRGjc78fpnDvP0X9u+58LnefE7U9stHqtPGRlrbX8Wnr97dcXe1sAAI2dMvYUjfzFSF3zyDXtZZPu\nmqQ1l1tT1xx6TYd5p/5zqr527dd0+s6na5kll+kwbb1V19NW62ylW5+4VZcddFmf1B0A0Jom3DlB\nX9nyK/r8Rz7fXhYROv7Px2vCtAntCZIPrvZB3f7U7Tpq9PzE+tx5czXjmRn67OjP9nm9y0CCpA+0\nuW2xrv+kOGmh8wwbNkxPP/20nnrqKa222mrt5aNHj9aMGTM0a9YsDR06dHFWE33s6oeu1gtvvKCj\nNztaqyyzSodpB210kM68/UyNes+ohc5z4vYn6le3/0oznpmhvT+0t0asMkJz3p2jiTMm6p5n79G3\nPv6tvmwWAKBg3VXX1TGbHaPT/366bCsiNOHOCdp/5P4aucbIDvOut+p6+vp1X9dF91zU4WS25ppD\nr9Hbc9/ut7/6AQAWTSh01zN36Y6n7tA5e56zwOfIuFHj9NkrPqvTdzldyy+1vL76sa/qqMuP0oZr\nbqidRuykN955Q2f8/QzNnjNbx2x+TJNasWgYg6QibGvEiBG68MIL28vuvvtuvfnmm9zOd4A6e/rZ\nGjt87AKJD0naf+T+mjV7lg6ZfEin8+w3cj89/vLjuv6x67XV+7bSG++8oS9e/UVt/MuNtf252+vm\nf9ys8/c+v30AQADA4md7geu6/2v7/9KSg5aUJD320mOa/q/p2m/kfgssu9QSS2mPD+7Rfg251XFd\nyyy5TIfkSP10AMDAVft8mTBtgjZYbQNttOZGC8yz6/q7al7M00X3XCQp/aB6zp7n6Jzp52iL32yh\nXSbtomffeFY3H3Wz1lxuzb5uQinM7Vx7x3Y02ne1X2+KWqEHyfDhw/W5z31Ol19+uW67LQ3Qdtxx\nx2nVVVfVt7/97YY9SBq1BeWrJajY1yhD/r+txDeazuIw0FPEYZSNWAz0HLEYZeptHKYHSYV89KMf\n1SuvvKIHHnhAc+fO1cUXX6zDDjus2dUCAAAAAKDpGIOkYsaNG6eJEydqu+2208iRI7XOOus0u0oA\nAAAAADQdCZJFsLgvnSmbbY0bN07bbrutZs6cqcMPP3yhXdj6Wxv7s97s65MmjS+/IkA/QoxCmXp7\nPBGLUXXEYpSJc+LW4IfT31i/ufXoa1xiUzFDhw7ViBEjdM0112ifffZpdnUAAAAAAGgJ9CCpoAkT\nJmj27NlaZpll9O677za7OgAAAAAANB0JkgoaMWJEh9fc5hc1Ve1KBwADGbEdAIDuIUHSB7pzG97F\nbebMmQ3LBw8erLlz5/ZxbQAAAAAAaC2MQQIAAAAAACqPHiQAAAAAAKBdVS/LpAcJAAAAAACoPBIk\nAAAAAACg8rjEBkC7qnalA4CBjNgOAED30IMEAAAAAABUHgkSAAAAAABQeU1LkNg+wfY9tmfYvtP2\nliWu+0jbZ5S1PgAAAAAAqsIPp0fVNGUMEtsfk7SrpNER8Y7tVSUtXeImorcL2h4UEfNKrAsAAAAA\nAGhxzRqkdS1Jz0fEO5IUES9Kku1Zki6WtIukNyUdEhGP2l5D0i8lDc3L/3tETM29Tn4iaUie/6iI\neKi4Idu7SjpB0u6SNpc0XikZ82ie//W83YskfULS922/R9LnJb0r6b6IOHiRWnuBF2nxhTpk4fmg\nYcOG6dlnn9USSywhSbKthx56SGuttdbirRsAAAAAAP1AsxIk10r6L9sPSrpe0sURcZNSz4/ZETHK\n9jil5Mfukk6X9OOIuMX2UEl/lDRS0v2Sto2IubZ3knSKpP0kWZJs7y3pP5QSLksqJUp2jIg3bX9T\n0lclnZy3+3xEbJ6Xe1LSsNy7ZcW+2CGLm21dddVVGjt2bLOrghZW60bHHQ8AYOAgtgMA0D1NSZDk\nXhubS9pW0hhJF9v+Vp58Yf57kaQf5+c7Sfqw3d4TYwXby0paWdJE2+spJTmK7Rkr6SOSPhERr9ne\nTSmpMjWvZylJUwvzX1x4fpekC2xfJumyztpxg25ofz5MwzRcw7vRegAoz5T7pCn3N7sWzUMcBtAK\niMXEYgDNVVYcblYPEuVxPm6UdKPtuyUd2Wi2/NeStoqIt4sTbf9C0p8jYm/bH5A0pbDco5KGS/qg\npDty+XURcUgnVXq98HxXSdsp9V45wfbGETG3foExGtNlG1tNRK+HZgHQonYYmR41bZObV5dm6G9x\nGMDARCwmFgNorrLicFPuYmN7A9vFjp6jJc3Kzw8s/K318LhW0rGF5TfJT1eU9FR+flRxE5IeV7rc\nZqLtkZL+Lmkb2+vmdSxXV4faui1paERMkXS8pJUkLdfzVraWiNBee+2lVVZZRaussor22WefZlcJ\nAAAAANCCYv1qXprZrB4ky0s6w/bKSgOhPqw0KOpuklaxPUPSHEm1wVGPlfTzXD5YqefJlyR9X9J5\ntr8t6WrN73ESkiIiHrR9qKRL8rqPlHSh7dodc07I2y5aQtL5tldSSrScHhGvlNn4ZrCtyy+/nDFI\nAAAAAABooFljkEyTtE19eR4b5PsRcXzd/C9IOqjBev6mdAlNzYm5/DxJ5+Xn0yVtmKfPlLRlg/UM\nLzx/V2lsFAAAAAAAUBFNG4OkEwySATRRFbvRAcBAR2wHAKB7WipBEhEjml0HAAAAAABQPS2VIBmw\nDml+x5iZM2c2uwoAAAAAALSsptzFBgAAAAAAtCY/nB5VQ4IEAAAAAABUHgkSAAAAAABQeSRIALSr\nalc6ABjIiO0AAHQPCRIAAAAAAFB5JEgAAAAAAEDlcZtfAAAAAADQLtZvdg2agwTJYmC72VUAAAAA\nAAA9QIKkZOM1vtlVAAAAAAAAPeSIaHYd+iXbwb5DGWo9jjieUAbbiohKdGMjDqMsxGGUjVgM9Byx\nGGXqbRzOgpmTAAAgAElEQVRmkFYAAAAAAFB5JEgAAAAAAEDlkSABAAAAAADt3Ga5rRJXCnZAggQA\nAAAAAFQeCRIAAAAAAFB5JEgAAAAAAEDlkSABAAAAAACVR4IEAAAAAABUniOi2XXol2wH+w5lsNPo\n0BxPKINtRUQlhhwnDqMsxGGUjVgM9ByxGGXqbRymBwkAAAAAAKg8EiQAAAAAAKDySJAAAAAAAIDK\nI0ECAAAAAAAqjwQJAAAAAACoPBIkAAAAAACgndsst1XiZlwdkCABAAAAAACVR4IEAAAAAABUHgkS\nAAAAAABQeSRIAAAAAABA5ZEgAQAAAAAAleeIaHYd+iXbwb5DGew0OjTHE8pgWxFRiSHHicMoC3EY\nZSMWAz1HLEaZehuH6UECAAAAAAAqjwQJAAAAAACoPBIkAAAAAACg8kiQAAAAAACAyiNBAgAAAAAA\nKo8ECQAAAAAAaOc2y22VuBlXByRIAAAAAABA5ZEgAQAAAAAAlUeCBAAAAAAAVB4JEgAAAAAAUHkk\nSAAAAAAAQOU5Ippdh37JdrDvUAY7jQ7N8YQy2FZEVGLIceIwykIcRtmIxUDPEYtRpt7G4cGLozJV\n0ea2ZlcBAwjHU3lOmjS+2VVAH+H/BmXieCoPcbha+N9BmTieykMs7jkusQEAAAAAAJVHggQAAAAA\nAFQeCRIAAICKGj9+vMaPH9/sagAA0BJIkAAAAAAAgMojQQIAAAAAANr54fSoGhIkAAAAAACg8kiQ\nAAAAAACAyiNBAgAAAAAAKm9wsysAAACA5uAONgAAzEcPEgAAAAAAUHn0IAEAAAAAAO1i/WbXoDn6\nXQ8S22vZvsj2I7Zvt3217aNtX9nD9bTZ3nFx1RMAAAAAAPQf/aoHiW1L+r2kcyLioFw2StIePV1X\nRJzUyTYGRcS8RaooAAAAAADoV/pbD5Ixkt6OiF/XCiLiLkk3S1re9iW277f9v7Xptk+0fZvtu23/\nqlB+ru198/NZtk+1fYek/W0fa/te2zNsX9h3zQMAAAAAAM3Qr3qQSNpI0h0Nyi1ptKSRkp6WdIvt\nbSLiFkk/i4iTJcn2RNu7RcRVkiI/lP8+HxGb5/melDQsIt6xveLibRIAAEBz1O5iw91sAADofwmS\n6GLabRHxlCTZni5pmKRbJI21/XVJy0paVdI9kq5qsPzFhed3SbrA9mWSLutsgzfohvbnwzRMwzW8\ne60AgJJMuU+acn+za9E8xGEArYBYTCwG0FxlxeH+liC5V9J+nUx7q/B8rqQlbA+R9HNJm0fEk7ZP\nkjSkk+VfLzzfVdJ2knaXdILtjSNibv0CYzSmp/UHgFLtMDI9atomN68uzUAcBtAKiMXEYmCg8cPp\nb3+5m01ZcbhfjUESEX+RtLTto2tleZDWbTtZpJYMecH28pL2X9g28kCwQyNiiqTjJa0kablFqTcA\nAAAAAGht/a0HiSTtLekntr8paY6kmZIuV4PLbyJitu3fKF1W8y9Jf+/G+peQdL7tlZTGNjk9Il4p\nq/IAAAAAAKD19LsESUQ8LenABpPOKszzlcLzEyWd2GA9RxWeDy88f1ed90gBAAAAAAADUL9LkAAA\nAKAc3L0GAID5+tUYJAAAAAAAAIsDPUgAAAAAAEC7/nL3mrLRgwQAAAAAAFQeCRIAAAAAAFB5JEgA\nAAAAAEDlkSABAACoqPHjx3MnGwAAMhIkAAAAAACg8kiQAAAAAACAdn44PaqGBAkAAAAAAKg8EiQA\nAAAAAKDySJAAAAAAAIDKc0Q0uw79ku1g36EMtiVJHE8og21FhJtdj75AHEZZiMMoG7EY6DlicWtx\nW34/Tuqf70dv4zA9SAAAAAAAQOXRg6SXyJajLGTLUSZ+tQR6jjiMshGLgZ4jFqNM9CABAAAAAADo\nJRIkAAAAAACg8kiQAAAAAACAyiNBAgAAAAAAKo8ECQAAAAAAqDwSJAAAAAAAoJ3bLLdV4mZcHZAg\nAQAAAAAAlUeCBAAAAAAAVB4JEgAAAAAAUHkkSAAAAAAAQOWRIAEAAAAAAJXniGh2Hfol28G+Qxns\nNDo0xxPKYFsRUYkhx4nDKAtxGGUjFgM9RyxGmXobh+lBAgAAAAAAKo8ECQAAAAAAqDwSJAAAAAAA\noPJIkAAAAAAAgMojQQIAAAAAACqPBAkAAAAAAGjnNsttlbgZVwckSAAAAAAAQOWRIAEAAAAAAJVH\nggQAAAAAAFQeCRIAAAAAAFB5JEgAAAAAAEDlOSKaXYd+yXaw71AGO40OzfGEMthWRFRiyHHiMMpC\nHEbZiMVAzxGLUabexmF6kAAAAAAAgMojQQIAAAAAACqPBAkAAAAAAKg8EiQAAAAAAKDySJAAAAAA\nAIDKI0ECAAAAAADauc1yWyVuxtUBCRIAAAAAAFB5g5tdgf6szW3NrgIGkLKOp5MmjS9lPUB/QBxG\nmYjDQO8Qi1GmMo4n4nCJLqhWLxJ6kAAAAAAAgMojQQIAAAAAACqPBAkAAADkh9MDAICqYgwSAAAA\nAADQLtZvdg2agx4kAAAAAACg8kiQAAAAAACAyiNBAgAAAAAAKo8ECQAAAAAAqDwGaQUAAEBlB+QD\nAKCGHiQAAAAAAKBdVW/9ToIEAAAAAABUXp8lSGzvZXue7Q/m15vY3qUwfXvbHytxe+Ntf62s9QEA\nAAAAgIGrL3uQHCzpqvxXkkZL+nRh+hhJW5e4vShxXQAAAAAAYADrkwSJ7eUlbSXpy5IOtL2kpO/k\n53fa/oakz0v6j/x6G9tr2P6d7dvyY+u8rvG2z7Z9g+1HbX+lsJ0TbD9o+2ZJHyyUr2v7Gtu3276p\n0Iuls21sn+txp+1puf4AAAAAAGCA6qu72Owp6Y8R8Q/bz0naWNKJkjaPiGMlyfYykl6NiB/l1xdI\n+nFE3GJ7qKQ/ShqZ17eBUo+TFSU9aPsXkjaVdKCkTSQtKWmapNvz/L+W9PmIeMT2VpJ+IWlHSad3\nso2vSfpSRNxqe1lJby22PQMAANACaoPxcTcbAEBV9VWC5GBJP87PL8mv75HkuvmKr3eS9GG7vWgF\n28spXTpzdUS8I+kF289KWkvStpImR8QcSXNsXyFJeZmtJV1SWNdSC9nGLZJ+bHtSXueTjRp1g25o\nfz5MwzRcw7uxKwCgPFPuk6bc3+xaNA9xGEArIBYTi4GBpr8ly8uKw4s9QWJ7VaXeHhvZDklLKCU5\n7l3YopK2ioi369YnScWyuUrtCHVMsNSeD5L0UkSM7u42JH3P9lWSdpV0i+1PRcSD9QuP0ZiFNAEA\nFq8dRqZHTdvk5tWlGYjDAFoBsZhYDKC5yorDfTEGyX6SJkbEsIgYHhFDJc2UNFTSCoX5Xq17fa2k\nY2svbG/SxTZC0k2S9rI9xPYKknaTpIh4VdJM2/vl9dj2qE62sWn+u25E3BsR35f0fyqMZwIAAAAA\nAAaevkiQHCTp93VllypdFjMyD4S6v6QrJe1dG6RVKXHxEdszbN+rNIhrzQJ3qImIOyVdLGmGpD9I\nuq0w+VBJn7U9XenSnj1yef02jsnl/8/23bZnKPVWuaa3jQcAAAAAAK1vsV9iExFjG5Sd0cns9b1E\nDmqwbFvd640Lz0+RdEqDZWZJ2qVB+QudbOPY+jIAAAAAADBw9dUgrQAAAGhh/W1APgAAytYXl9gA\nAAAAAIB+wg/Pv/17lZAgAQAAAAAAlUeCBAAAAAAAVB4JEgAAAAAAUHkkSAAAAAAAQOWRIAEAAEBl\nB+QDAKCG2/wCAAAAAIB2Vb31Oz1IAAAAAABA5ZEgAQAAAAAAlUeCBAAAAAAAVB4JEgAAAAAAUHmO\niGbXoV+yHew7lMG2JInjCWWwrYhws+vRF4jDKAtxGGUjFgM9RyxGmXobh+lBAgAAAAAA2rnNclsl\n8rwdkCABAAAAAACVR4IEAAAAAABUHgkSAAAAAABQeSRIAAAAAABA5ZEgAQAAAAAAlcdtfnuJW5qh\nLNzSDGXi1pJAzxGHUTZiMdBzxGKUidv8AgAAAAAA9BIJEgAAAAAAUHkkSAAAAAAAQOWRIAEAAAAA\nAJVHggQAAAAAAFQeCRIAAAAAANDObZbbKnEzrg5IkAAAAAAAgMojQQIAAAAAACqPBAkAAAAAAKg8\nEiQAAAAAAKDySJAAAAAAAIDKc0Q0uw79ku1g36EMdhodmuMJZbCtiKjEkOPEYZSFOIyyEYuBniMW\no0y9jcP0IAEAAAAAAJVHggQAAAAAAFQeCRIAAAAAAFB5JEgAAAAAAEDlkSABAAAAAACVR4IEAAAA\nAAC0c5vltkrcjKsDEiQAAAAAAKDySJAAAAAAAIDKI0ECAAAAAAAqjwQJAAAAAACoPBIkAAAAAACg\n8hwRza5Dv2Q72Hcog51Gh+Z4QhlsKyIqMeQ4cRhlIQ6jbMRioOeIxShTb+MwPUgAAAAAAEDlkSAB\nAAAAAACVN7jZFejP2tzW7CpgAGn14+mkSeObXQVgAa3+f4P+pT8cT8RitKL+8L+D/qPVjyfi8MBG\nDxIAAAAAAFB5JEgAAAAAAEDlkSABUAl+OD0AAK2F+AwAraeqsZkECQAAAAAAqDwSJAAAAAAAoPJI\nkAAAAAAAgMojQQIAAAAAACqPBAkAAAAAAKi8wc2uAAD0hVi/2TUAADRCfAaA1lPV2EwPEgAAAAAA\nUHkkSAAAAAAAQOUtcoLE9vtsX277IduP2P6J7SXLqFxe/562P1x43WZ7x7LWDwAAAAAAsEgJEtuW\nNFnS5IjYQNIGkpaX9N8l1K1mb0kjay8i4qSI+HOJ6+/ANuOyAAAAAABQMYvag2SspDcj4jxJioh5\nkv5D0mdsf9H2GbUZbV9le/v8/JO2p9q+w/ZvbS+Xy0+1fa/tGbZPs/0xSbtLOs32NNsjbJ9re988\n/465/C7bE2wvlctn2R6f13+X7Q/m8uVsn23773m5PXL5kbavsP1nSdfZXsv2TbbvtH237Y8v4n4C\nAAAAAAAtbFETJBtKuqNYEBGvSvqHpCXq5g1JYXt1SSdI2jEiNs/Lf9X2qpL2iogNI2ITSSdHxK2S\nrpB0XERsFhGPFdYzRNI5kg6IiFFKd+T5YmFbz+X1/1LScbn8BEl/joitlJI7p9leNk8bLWnfiBgj\n6VBJf4yI0ZJGSZq+iPsJQJP54fQAALQW4jMAtJ6qxuZFvZwkupjWaBwSS/qo0iUzU9MVOlpK0lRJ\nL0uaY3uCpKvyo7hc/Xo+KGlmRDySy86T9G+STs+vJ+e/0yTtk59/UtLutmsJk6UlDc3tuC4iZufy\n2ySdncdSuSwiZjRq4A26of35MA3TcA1vNBsALDZT7pOm3N/sWjQPcRhAKyAWE4sBNFdZcXhREyT3\nSdqvWGB7RUnvl/ScpPUKk4YUnl8XEYfUr8z2lpJ2zOv8cn4uNU7E1Je5ruyt/HeuOrZzn4jokAuz\nvZWk19tXHHGz7W0l7SbpXNs/iojz6yswRmMaVAsA+s4OI9Ojpm1y5/MORMRhAK2AWEwsBtBcZcXh\nRbrEJg+WuqztcZJkewlJP5R0gaSZkjZ18n5JWyolMP4maRvb6+ZllrO9fh6HZOWIuEbSVyVtkjfz\nqqQV6zct6UFJw2rrkTRO0o0LqfKfJB1be2F7dO1pcSbbQ5Uu0TlL0llKl98AAAAAAIABapFv86t0\nl5n9bD8k6XmlZMZxEXGLUpLkPqXLXu6QpIh4XtKRki60PUPp8poPSlpB0pW57GalwV4l6SJJX88D\nro6obTQi3pJ0lKRLbN8l6V1JZ9YmF+oXhdcnS1oyD9x6j6S2BvNI0g6SptueJukAzb9sBwAAAAAA\nDECLfEvbiHhC0p6SlO868xulS2vuj4jDOlnmBqUeJfW2ajDvVKXBYGuOKkz7i6TNGixTTKTcoTQg\nqyJijqQvNJj/PKUxTGqvJ0qa2KjuAAAAAABg4FnkBElRvuvMRmWuEwDKEOs3uwYAgEaIzwDQeqoa\nm8u4xAYAAAAAAKBfI0ECAAAAAAAqjwQJAAAAAACoPBIkAAAAAACg8kiQAAAAAACAyiNBAqAS/HB6\nAABaC/EZAFpPVWMzCRIAAAAAAFB5JEgAAAAAAEDlkSABAAAAAACVR4IEAAAAAABUHgkSAAAAAABQ\neY6IZtehX7Id7DuUwbYkieMJZbCtiHCz69EXiMMoC3EYZSMWAz1HLEaZehuH6UECAAAAAAAqjwQJ\nAAAAAACoPBIkAAAAAACg8kiQAAAAAACAyiNBAgAAAAAAKo8ECQAAAAAAaOc2y22VuBlXByRIAAAA\nAABA5ZEgAQAAAAAAlUeCBAAAAAAAVB4JEgAAAAAAUHkkSAAAAAAAQOU5Ippdh37JdrDvUAY7jQ7N\n8YQy2FZEVGLIceIwykIcRtmIxUDPEYtRpt7GYXqQAAAAAACAyiNBAgAAAAAAKo8ECQAAAAAAqDwS\nJAAAAAAAoPJIkAAAAAAAgMojQQIAAAAAANq5zXJbJW7G1QEJEgAAAAAAUHkkSAAAAAAAQOWRIAEA\nAAAAAJVHggQAAAAAAFQeCRIAAAAAAFB5johm16Ffsh3sO5TBTqNDczyhDLYVEZUYcpw4jLIQh1E2\nYjHQc8RilKm3cZgeJAAAAAAAoPJIkAAAAAAAgMojQQIAAAAAACqPBAkAAAAAAKg8EiQAAAAAAKDy\nSJAAAAAAAIB2brPcVombcXVAggQAAAAAAFQeCRIAAAAAAFB5JEgAAAAAAEDlDW52BfqzNrc1uwoY\nQJp9PJ00aXxTtw/0RrP/bzCwtMLxRCxGf9QK/zsYOJp9PBGH61xQrXFI6EECAAAAAAAqjx4kAAAA\nAACgXazf7Bo0Bz1IAPQ7fjg9AADVQvwHACxOJEgAAAAAAEDlkSABAAAAAACVR4IEAAAAAABUHgkS\nAAAAAABQeSRIAAAAAABAu6oOis1tfgH0O1W97RgAVB3xHwCwONGDBAAAAAAAVB4JEgAAAAAAUHkL\nTZDYnmv7zvyYZvsDtm/pyUZs/7vtZXpfzd6zfYTt9xZe/8b2h5tRFwAAAAAA0Jq6MwbJGxExuq5s\nm/qZbA+OiHc7Wcf/k3S+pDd7WL8yHCnpHklPS1JEHN2EOgAAAAAAgBbWq0tsbL+W/+5g+2bbl0u6\nx/aytq+2Pd323bYPsP0VSWtLusH2nxusa2fb99u+w/ZPbV+Zy8fb/lphvntsD83PD7P999yr5Uzb\ng2wvYfvcvN27cq+VfSV9RNKk3PtliO0ptjfP6zk4z3u37VOL7bP93dyOW22v2Zv9BAAAAABAfxPr\nV3Ng7O4kSJYpXGJzaS6LwvTRko6NiA9J2kXSkxGxaURsLOmaiDhD0lOSdoiIHYsrtj1E0q8l7RYR\nm0t6T2HdxW20v86Xxxwgaevcs2WupEMlbSJp7YjYOCJGSTo7Ii6VdLukQyJis4iYk9cTtteWdKqk\nMZI2lbSF7T3ztpaVdGtEbCrpJkn0OgFaSFVvOwYAVUf8BwAsTt25xObNBpfYFN0WEY/n53dJ+kHu\njXFVRPx1Iev+kKSZEfFofv2/ko7pYn5L2lHS5pJuty1Jy0h6RtKVkkbY/qmkqyVdW7dc/Xq2kDQl\nIl6QJNuTJG0n6XJJb0fE1XneOyR9olFlbtAN7c+HaZiGa3hXbQWA0k25T5pyf7Nr0TzEYQCtgFhM\nLAbQXGXF4e4kSBbm9dqTiHjY9mhJu0r6ru0/R8TJXSxb30ukmMh4Vx17uAwpPD8vIv6zfmW2R0na\nWdIXlHqZfLaT7XS27VrZO4XyeepkP43RmEbFANBndhiZHjVtk5tXl2YgDgNoBcRiYjGA5iorDpd6\nm998t5g5ETFJ0g+ULr+RpFclrdhgkQclDbM9Ir8+WPOTFLMkbZbXu5mk4XnanyXtZ3uNPG1V20Nt\nryZpcERMlnTiQrYdkm6TtL3t1WwvIekgSTf2tu0AAAAAAKD/6k4PkoX1vig+31jSabbnKfXC+EIu\n/7WkP9p+sjgOSUTMsX2MpKttvyHpZknr5smXSjrc9j2S/q6UTFFE3G/725KutT0ob+dLkuZIOieX\nSdLx+e+5ks7M69+6sO1/2T5e0g1KvUeuiogrO2lfo30AAAAAAAAGiIUmSCJigZ4ftbKImCJpSqH8\nWnUc+6NW/jNJP+tk/X+S9GFJsr29pONy+RxJn+pkmd9K+m2DSZs3mHeypGIHmzGFaRdJuqjBMisW\nnl+qlKwBAAAAAGDAqw2IXbU72ZQxBknZ6K0BoEtVC9QAgIT4DwBYnFoqQRIRN4pxQAAAAAAAQB8r\ndZBWAAAAAACA/ogECQAAAAAAqDwSJAAAAAAAoPJaagwSAAAAAADQXFUdFJseJAD6HT88/9ZjAIDq\nIP4DABYnEiQAAAAAAKDySJAAAAAAAIDKI0ECAAAAAAAqjwQJAAAAAACoPBIkAAAAAACgXVUHxXZE\nNLsO/ZLtYN+hDLYlSRxPKINtRYSbXY++QBxGWYjDKBuxGOg5YnFrcVt+P07qn+9Hb+MwPUgAAAAA\nAEDlkSABAAAAAACVR4IEAAAAAABUHgkSAAAAAABQeQzS2ksMSIWyMCAVysTAgEDPEYdRNmIx0HPE\nYpSJQVoBAAAAAAB6iQQJAAAAAACoPBIkAAAAAACg8kiQAAAAAACAyiNBAgAAAAAAKo8ECQAAAAAA\naOc2y22VuBlXByRIAAAAAABA5ZEgwUJNmTKl2VXoM7R1YKpSWzEwVekYrlJbpWq1t0ptxcBUpWOY\ntg5MVWprb5EgwUJV6R+Jtg5MVWorBqYqHcNVaqtUrfZWqa0YmKp0DNPWgalKbe0tEiQAAAAAAKDy\nSJAAAAAAAIDKc0Q0uw79km12HICWFBGVGHKcOAyglRGLAaC5ehOHSZAAAAAAAIDK4xIbAAAAAABQ\neSRIAAAAAABA5ZEgWQjbO9t+wPbDtr/ZyTw/zdNn2B7d13Usy8LaavvQ3Ma7bN9ie1Qz6lmG7ryv\neb7/3979B9tR1nccf38JiRaHEGjGEAgxARWKHUEQkAJNUpSfKg60Q6FAArWgM8kUnYo/2obUXzjW\n6QASKCIEiy2ogG1GA1YlKT8EKyWJoEQISSRACCYkgCTRJPfbP57nJHv37p6759w9P/acz2tm59y7\n++zZ59lnz+fcee7+OMbMdpjZ2e2sX5kKHsPTzWypmT1hZkvaXMXSFDiGx5vZvWa2LLZ1VgeqOWJm\ndouZrTezx+uU6YlcAuVwarlyuIKUw4OW90QOg7I4p0xPtFdZnFlOWVwh/ZLFLclhd9eUMwGjgJXA\nFGA0sAz4o1SZM4BF8efjgEc6Xe8WtvV4YJ/482m93NZEufuA7wHndLreLezXccAvgEnx9/GdrncL\n2zoPuKrWTmAjsGen695EW08C3gU8nrO8J3KpgX7tifYqh5XDyuFqTcpiZXH8WVlcgUlZ3JtZ3Ioc\n1hkk9R0LrHT3Ne6+HbgDOCtV5oPANwDc/afAODOb0N5qlmLYtrr7w+7+Svz1p8CkNtexLEX6FWAO\ncCfwm3ZWrmRF2no+cJe7Pwfg7hvaXMeyFGnrOmBs/HkssNHdd7SxjqVw9weATXWK9EougXJYOawc\nrpK+yWFQFqMsVhZXi7J4sJ7I4lbksAZI6jsQWJv4/bk4b7gyVQzJIm1N+mtgUUtr1DrDttXMDiQE\nyQ1xVlUf91SkX98G7Gdmi83sUTO7sG21K1eRtt4EvMPMXgCWA3/bprq1W6/kEiiHlcPK4SpRDg/W\nK9kEymJlsbK4SpTFuzWcS3u2tDrVVzQA0s9XrmJwFK6zmc0ALgFOaF11WqpIW68GPuXubmbG0D6u\niiJtHQ0cBZwM7AU8bGaPuPvTLa1Z+Yq09TPAMnefbmaHAD80syPc/bUW160TeiGXQDmcSTlcKcrh\nwfoph6E3sgmUxZmUxZWiLB6sn7K4oVzSAEl9zwMHJX4/iDDqVK/MpDivaoq0lXgTqpuA09y93ulM\n3axIW48G7gjfA4wHTjez7e6+sD1VLE2Rtq4FNrj7VmCrmd0PHAFU7cugSFv/BPgCgLs/Y2argUOB\nR9tSw/bplVwC5bByWDlcJcrhwXolm0BZrCxWFleJsni3hnNJl9jU9yjwNjObYmZjgHOBdBgsBC4C\nMLP3AJvdfX17q1mKYdtqZpOBu4EL3H1lB+pYlmHb6u4Hu/tUd59KuObyoxX8IoBix/B/ASea2Sgz\n24twA6NftrmeZSjS1hXAewHi9YeHAqvaWsv26JVcAuWwclg5XCXK4cF6JZtAWawsVhZXibJ4t4Zz\nSWeQ1OHuO8xsNvADwt2Ab3b3J83ssrj8RndfZGZnmNlK4HXg4g5WuWlF2grMBfYFboijyNvd/dhO\n1blZBdvaEwoewyvM7F7g58AAcJO7V+7LoGC/fhFYYGbLCQPEV7j7yx2rdJPM7HZgGjDezNYCVxJO\nC+2pXALlsHK4+pTDvZnDoCxWFiuLq0RZ3JtZ3IocNvcqXhooIiIiIiIiIlIeXWIjIiIiIiIiIn1P\nAyQiIiIiIiIi0vc0QCIiIiIiIiIifU8DJCIiIiIiIiLS9zRAIiIiIiIiIiJ9TwMkIiIiIiIiItL3\nNEAiHWdm88xswMymdbouI2FmY83sWjNbY2Y7Ypve2el6jVQZ/RPXX1xmvUSkXMri7qYsFul9yuHu\nphzuDxogqbj4Ias3zeyCOs4api6emKrsy8BsYDnwBWAesL6TFSqijf1Tif41s1vj/pjc6bpIdSiL\nu4qyuL5K9K+yWBqlHO4qyuH6KtG//ZrDe3a6AlIKB/4pZ9nSdlZkGHlhcB1wO7C2jXVphfcDv3L3\nszpdkSa1sn8OA7aMYP12q8QXl3QdZXF3UBbnUxZLr1MOdwflcD7lcJfTAEmPcPfPdroOBVjWTHff\nCGxsc11a4QBgSacrMQIt6x93f2ok63eAkbM/ROpRFncFZXEOZbH0A+VwV1AO51AOV4C7a6rwBAwA\nO4cpMyWWWwC8HfgW8BKwE5gWyxwNXEM4FW4jsBV4CvgKMK7Oe58L/Bh4Oa6zGvgP4Oi4fEncdtY0\nOQr+ZfUAAAo4SURBVJaZF3//04z3Pxm4N77/NuBXwFXA2IyytW2NAj4DPB3XeRb4EjC6wX07EZgP\nrAF+F/fZXcBROdtNT4sLbGNN3GdjCaPSz8f9+AtgTs46s2I9VhFGoF8BHgT+Kqd8rX6jgblxH26L\nx8NI++cw4JbYjm2E0yfvBz6ScZwuTs2rve80YCbhPztb4nvcDEzI2F5Dx2ncVwPx/WfE9r4a99n3\ngMMy6pk1rU6UORj4GrAy1ncj8HPgBmC/TmeCps5MKIuTZWvbUhZn109ZrCzW1IIJ5XCybG1byuHs\n+imHlcO5k84g6S+HAI8QwuA24A8IHwqAvwE+RPiw/Dfh/jTvBj4OnG5mx7n7b2tvZGZGCJOLgN8A\nd8bXg4DpwArg/2KZTcBZwH8CyxL1eYU6zOwywgfsNeA7hDCeAXwS+ICZneDuWe9xO3AisIjwwT8T\nuAJ4M3BJvW0mtj2VELATCV92/w5MBv4CONPMznH378fiC4DFwJWEULw1zl9TYFMOjInbGEv4In0D\ncA5wjZkd6u6zU+tcDzxB6Kt1wHjgDOC2WH5uzrbuJvTpovjzS/E9mu2fMwn9Mga4h7CPxgFHAp8A\n/jWjrVk+BpwC3BHrdhJwMTA9HncbEmUbOk4T3h/buIhwTL2DsM+OMbPDPfxHAMJpuR8CjgCuBjbH\n+ZtjmycCPwP2Br4f2/9GwhfEBcBXCX+4iNSjLFYWK4uVxdJZymHlsHJYOZyt0yM0mkY2sXs070rC\n6GNymhnLTEmU+3zO+0wGLGP+JXG9K1LzL43zHwH2Ti3bA9g/8fusWPainG3PIzUaC7yFMEK9GXh7\nqvz8WP7G1Pwlcf7PSIycAnsRRs53kDECm1OnH8T3+nRq/vHAdmAD8KaMvrivwf5bE9e7n8RoPrAv\nYTR2ADgptc7UjPcZDfwI+D1wQM5+WUbGaG6T/TOe8GWxLV2/uDxdhyH7JvG+24AjUsv+JS77+giP\n01rbfg/MSC37Ylz2idT8W0n8tyC1bE5cNuQ/GYQ/rt7YSP9r6p0JZXFyfi1zlMXZ+0VZPHiZslhT\nKRPK4eT8Wt4oh7P3i3J48DLlcGLSU2x6x5WEU8WS08xUmRfJuXGVuz/r8YhOWUAYrT4lNX8OYfTz\nMnd/LfVeA+7+YsMtGOwCQsBd50Ov1ft74LfABWY2JmPdT7p7bZQTd99CGMndg3A6Wl1mNgl4H/Br\nwl24d3H3hwmj8fsBZxduzfA+7e7bE9vZBHwu/npxqg6r0yvHda8n3Ffo5Jxt/KO7lzWSO5MwYnyD\nuz+QUZ8XGniv29x9eWrePMJ/Os5L9nETx2nNHe6+ODXva/H1mAbqWtv2tiEL3Le6+5D50neUxbsp\ni7MpiwdTFkvZlMO7KYezKYcHUw4naICkN7i775Ex/Vmq3PJk4CSZ2Wgzm21mD5rZy7VnlhOuydwb\nODBR9k2E07HWZ3yIy3JUfL0vvSAG/VLC6ORh6cXAoxnv91x83bfAtt8VXx9w950Zy2t1OrLAexWx\nHfhJxvwlWdsxs8lmNt/MVpjZ67XH1xFO6YRwY6w0B/63pPoCvCe+3lPCe/1Peoa7v0oY3R/Ux40c\npykjPSZqFhL+EJlvZnea2aVmdngD60tvUxYnFqMsVhYPpSyWVlMOJxajHFYOD6UcHobuQdJf6o1g\nf4twndkzwHdj2d8R7lp8OeEawJpx8fX5FtSxZp/4ui5n+bpUuV1ikKTtiK+jSth2bT+Oy1neqA05\nI8C158XvaqOZHUwI9XGEUxDvJZzWtxOYShjFfgMZ3L3M58+XeQzk1au2n5N93MhxmrQ5PcPdd4TL\nhgsdE7V1njWzYwmj+acR/2NiZmuBr7j7V4u+l/Q1ZbGyuCzKYmWxNEc5rBwui3K4x3JYAyT9JStw\nMLN3Ez5gPwROd/eBxDIj3AAqqfbByhuZLEPtZkgTgSczlk9MlWvFtvfPWV72tsebmWV8IdS2n9zO\nxwmnMs5y939LFjaz8xh6Cmmr1I6BSYS7i4/EhJz5g9rfxHHaEu6+AvhLMxtFuHHVewmn115jZq+7\n+y3tqIdUmrK4sW0ri/Mpi5XF0hzlcGPbVg7nUw73WA7rEhsBeGt8XZj8gEXHEe5IvIu7v064Y/T+\nZlbklLraKXmFRyWBx+Lr9PQCM6vdFXor2V8UI1Xb9onxA582I1VupEYDJ2TMnx5flybmvZXwpX5X\nRvlpTW6/mf55OL6e3uQ2k6anZ5jZPgzt44aO0xEotD/cfae7P+buXwbOi7PPKqkO0p+UxdnbVhbn\nUxYri6VcyuHsbSuH8ymHeyyHNUAiEJ45DrtDDgAzezPh7thZro2vN5rZ2NR6e5hZcqS59riotzRQ\np28SrkOcY2aHpJZ9jnBt3Tfzrh8dCXd/njAiO5VwitouZnYccD7hsVXfLXGzV1nixktmth/wD4Tg\nX5Aot5pw6ly6r04FPtzktpvpn28Qbhj1UTM7Kb3Qwk29irow44+KeYRHvN2e6ONmjtNm5O4PMzsq\nflGl1Y73LSXWQ/qPsjhBWVyIsngwZbGMlHI4QTlciHJ4sMrnsC6xEQiPAHsIONvMHoo/TyBcT7YC\neIEQQLu4+9djCFwIPG1mCwnPfD+A8GG9GfhsLP4TwofkcjP7Q3ZfX3dtzrWRuPuvzexywof8MTP7\nNuExYtMIN0N6kuzTxyxjXjM+QtgP/2xmpxCeX38Q4ZnvO4CL438NyrCOcJ3gE3E/jgH+nBAw8939\nwUTZ6wl38P6Omd0Z1/1j4FTg28C5Oduot1+a6Z+NZnY+4SZYi83sHuBxQoC/k3Ca4cF1W73bIuCh\n2McvAicS/nuwGvhUolzDx2mTfgT8HXCTmd1NuBP4JnefD1wEXGpmDwKrgE3AIcAHCHfxvrqE7Uv/\nUhYPpSxWFiuLpZ2Uw0Mph5XD/ZXD3gXPGtbU/ER4/vTOYcpMieVuqVNmX0LwriacwvU08HnCHZNX\nA6ty1jufcGfpzXG9Z4DbgCNT5U4lhM5rtToTn6lNeBzbThLPFE+s9z7C89dfJnzYngK+BIzNKLs4\nb18QrkPcSc5zzXPWOYAQvmsINz16CbgbOLpOXzTzzPdVhNH/6wh3kd5GuIZxds46xwM/jvvkVcKN\nqT5I+KIcAOYW3S8l9M/hhJHz5+I+Whe39+Hh9g2JZ8nH/llK+FJaT/hjYsJIj9Ph+j2vz4CPAb+M\nfTFQe1/g2HhMLCOMqm+Jx+TNwOGt/rxr6t4JZXGyrLJYWaws1tT2CeVwsqxyWDmsHG5ysthAEekA\nM1sDDLh70ZHlnmFm84C5wHR3v7/D1RGRPqYsVhaLSGcph5XD3UL3IBERERERERGRvqcBEhERERER\nERHpexogEemsfr7Gzenv9otI9+jnLFIWi0g36OccUg53Ed2DRERERERERET6ns4gEREREREREZG+\npwESEREREREREel7GiARERERERERkb6nARIRERERERER6XsaIBERERERERGRvqcBEhERERERERHp\ne/8PvmxZdCrSNysAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x10a2fe890>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Plot stacked bar chart\n",
    "fig = plt.figure(1,(18,6))\n",
    "ax1 = fig.add_subplot(131)\n",
    "\n",
    "df = pd.DataFrame(normdata,index=labels)\n",
    "dfplot = df.plot(kind='barh',stacked=True,ax=ax1,color=gencolors,legend=False)\n",
    "\n",
    "print df\n",
    "\n",
    "# Find data positions of the plots\n",
    "\n",
    "patches = dfplot.patches\n",
    "yc = [p.get_y() for p in patches]\n",
    "yc = yc[:int(len(yc)/2)]\n",
    "height = p.get_height()\n",
    "\n",
    "ylims1 = ax1.get_ylim()\n",
    "ax1.vlines(0.5,ylims1[0],ylims1[1],color='k',linestyle='-')\n",
    "\n",
    "def getfrac(m,f):\n",
    "    return m/(f+m)\n",
    "\n",
    "# Speaker/questioner data from 225th AAS Meeting (Seattle, WA)\n",
    "# Data from https://github.com/jradavenport/aas225-gender\n",
    "\n",
    "aas225_speakers = getfrac(83,51)\n",
    "aas225_questionaskers = getfrac(305,73)\n",
    "aas225_firstquestion = getfrac(102,32)\n",
    "\n",
    "ax1.vlines(aas225_speakers,yc[-1],yc[-1]+height,color='g',linestyle='--')\n",
    "ax1.vlines(aas225_questionaskers,yc[-4],yc[-4]+height,color='g',linestyle='--')\n",
    "ax1.vlines(aas225_firstquestion,yc[-5],yc[-5]+height,color='g',linestyle='--')\n",
    "\n",
    "ax1.text(aas225_speakers,yc[-1]+height,'AAS',ha='center',va='bottom',fontsize=14,color='g')\n",
    "\n",
    "ax1.set_xlabel('Fraction of participants',fontsize=20)\n",
    "\n",
    "p,l = ax1.get_legend_handles_labels()\n",
    "ax1.legend(p,l,loc='upper left')\n",
    "\n",
    "ax1.set_title('AGN2015 vs. speakers at AAS 2015 Winter')\n",
    "\n",
    "# Speaker/questioner data from National Astronomy Meeting 2014 (Portsmouth, UK)\n",
    "# Data from Pritchard et al. (2014) http://arXiv.org/abs/1412.4571\n",
    "\n",
    "ax2 = fig.add_subplot(132)\n",
    "dfblank = pd.DataFrame(normdata,index=[' ']*5)\n",
    "dfplot2 = dfblank.plot(kind='barh',stacked=True,ax=ax2,color=gencolors,legend=False)\n",
    "\n",
    "nam_speakers = getfrac(181,81)\n",
    "nam_chairs = getfrac(188,75)\n",
    "nam_attendees = getfrac(452,172)\n",
    "nam_questionaskers = getfrac(476,101)\n",
    "nam_firstquestion = getfrac(216,35)\n",
    "\n",
    "ylims2 = ax2.get_ylim()\n",
    "ax2.vlines(0.5,ylims2[0],ylims2[1],color='k',linestyle='-')\n",
    "\n",
    "ax2.vlines(nam_speakers,yc[-1],yc[-1]+height,color='g',linestyle='--')\n",
    "ax2.vlines(nam_chairs,yc[-2],yc[-2]+height,color='g',linestyle='--')\n",
    "ax2.vlines(nam_attendees,yc[-3],yc[-3]+height,color='g',linestyle='--')\n",
    "ax2.vlines(nam_questionaskers,yc[-4],yc[-4]+height,color='g',linestyle='--')\n",
    "ax2.vlines(nam_firstquestion,yc[-5],yc[-5]+height,color='g',linestyle='--')\n",
    "\n",
    "ax2.text(nam_speakers,yc[-1]+height,'NAM',ha='center',va='bottom',fontsize=14,color='g')\n",
    "\n",
    "ax2.set_xlabel('Fraction of participants',fontsize=20)\n",
    "\n",
    "ax2.set_title('AGN2015 vs. NAM 2014')\n",
    "\n",
    "# IAU individual members (as of Apr 2015)\n",
    "# Data from http://www.iau.org/administration/membership/individual/distribution/\n",
    "\n",
    "ax3 = fig.add_subplot(133)\n",
    "dfplot3 = dfblank.plot(kind='barh',stacked=True,ax=ax3,color=gencolors,legend=False)\n",
    "iau_frac = getfrac(9546,1803)\n",
    "ylims3 = ax3.get_ylim()\n",
    "ax3.vlines(0.5,ylims3[0],ylims3[1],color='k',linestyle='-')\n",
    "ax3.vlines(iau_frac,ylims3[0],yc[-1]+height,color='g',linestyle='--')\n",
    "ax3.text(iau_frac*1.02,yc[-1]+height,'IAU',ha='center',va='bottom',fontsize=14,color='g')\n",
    "\n",
    "ax3.set_xlabel('Fraction of participants',fontsize=20)\n",
    "\n",
    "ax3.set_title('AGN2015 vs. IAU individual members');"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Analysis of the demographics compared to other meetings/organizations:"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The data from the **225th AAS meeting** did not track the gender of either the session chair or the total number of attendees. The <span style=\"color:red\">AGN2015</span> meeting had a higher proportion of women for both the speakers (45% to 38%) and the people asking questions (35% to 19%). \n",
    "\n",
    "The **NAM 2014** meeting tracked data in an identical format to <span style=\"color:red\">AGN2015</span>. The <span style=\"color:red\">AGN2015</span> meeting had a higher percentage of women for speakers, attendees, and question askers. In particular, the gender distribution of people asking questions is twice as high for <span style=\"color:red\">AGN2015</span> (35%) as for NAM (17%). This is driven by the large disparity in gender distribution between attendees and questioners at NAM. The gender of the session chairs is the only category in which NAM has a higher percentage of females than <span style=\"color:red\">AGN2015</span>. However, we note that there were only 14 sessions at <span style=\"color:red\">AGN2015</span>, compared to 64 sessions at NAM. Some (although likely not all) of the high fraction of males at <span style=\"color:red\">AGN2015</span> may be due to the smaller sample size.\n",
    "\n",
    "Finally, we look at the gender representation compared to the most recent data for members of the **International Astronomical Union** (IAU). The IAU is 85% male/15% female, a distribution significantly in excess of every demographic category measured for <span style=\"color:red\">AGN2015</span>. Membership in the IAU is likely not indicative of the astronomical community, especially since the requirements involve nominations by national professional organizations and thus skew toward more senior astronomers (which are more likely to be male). "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## How many questions per talk were there? Was this affected by the gender of the speaker?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAswAAAGECAYAAAAiMnaGAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xt8VfW57/vvA5jI/VZAbjEgFmXRpbi0raWUeGNZ2bpx\n7SqKClhrtWq1QO32Hqa69imnS+X0eLRe8IIYadXdSmUrcaGRutqKVaqUQhUE8QIY7iAogTznjzkS\nk5CMzIQ5MuZIPu/Xa7wy55hjjvHNJDx5MuZv/oa5uwAAAADUr13cAQAAAIBcRsMMAAAAhKBhBgAA\nAELQMAMAAAAhaJgBAACAEDTMAAAAQIjIGmYzG2xmr5jZCjP7m5ldG6yfaWYfmdmyYDkzqgwAgMxR\ntwGgfhbVPMxmdoSkI9z9r2bWRdKbkiZIOl/SLne/O5IDAwCahboNAPXrENWO3X2jpI3B7d1mtlLS\nwOBhi+q4AIDmoW4DQP1aZAyzmRVKGiXpz8GqH5vZ22Y2x8x6tEQGAEDmqNsA8KXIG+bgbb1nJF3n\n7rsl3S9piKTjJW2QdFfUGQAAmaNuA0BtkY1hliQzO0zS85JecPfZ9TxeKOn37v61OuujCwUALcDd\nEzmEgboNoC1qrGZHOUuGSZoj6e81i66Z9a+x2bmSltf3fHdPzFJcXBx7htaYNWl5k5Q1aXmTlNU9\nuX1jW6nbSft5SlLeJGVNWt4kZU1a3kxE9qE/SaMlXSzpHTNbFqy7SdKFZna8JJe0VtIVEWYAAGSO\nug0A9YhylozXVP8Z7BeiOiYAoPmo2wBQP670lwVFRUVxR8hYkrJKycqbpKxSsvImKStyX9J+npKU\nN0lZpWTlTVJWKXl5GxPph/6ay8w8F3OhZaWHUyrj8UVArjAzeUI/9Ndc1O3WK2UpSVKxF2d3xyXB\nf5FJ/NwgXpnU7CjHMAMAAESu6gQL0Jjm/mFPwwwAABKPdzjQmEP5w4oxzAAAAEAIGmYAAAAgBA0z\nAAAAEIJZMgAgy5glA2hZwf85SV/O6tESsj5zSAPWrVunoUOHav/+/WrXLvfPdU6dOlWDBw/WHXfc\nEXeUWmr+nNSzPp5LYwMAALR1hYWFys/P15YtW2qtHzVqlNq1a6f169fHlCw6ZtbqZi6hYQYAAIiI\nmWno0KF66qmnqtctX75ce/fubXVNZU3ZeMepsrIyC0myg2nlAABAqxTFkInmDPm4+OKLNXfuXF1z\nzTWSpMcff1yTJ0/WLbfcUr3NwoULdcstt+j9999X9+7dddlll6m4uP78O3bs0PTp0/XCCy+oXbt2\nuvTSS5VKpeodrrF06VJdddVVeu+999SxY0dddNFFuuuuu6qHeTzwwAOaOXOm3F0zZszQjBkzJKUb\n3lmzZunhhx/W9u3bddppp+lXv/qVevbsKUk677zz9Nprr2nv3r067rjjdP/992vEiBEHHX/Xrl06\n55xzdNxxx2n27NlatWqVfvzjH+utt95Snz59dMcdd+i8886TlB7K0bFjR33wwQdasmSJFixYoFNP\nPbXJr3cUOMMMAAAQoW9+85vauXOnVq1apQMHDujXv/61Lr744lrbdOnSRfPmzdOOHTu0cOFC3X//\n/Xruuefq3d/UqVOVl5enNWvWaNmyZSotLdXDDz9c77bXXXedpk2bph07duj999/X+eefX+vxsrIy\nrV69WqWlpZo1a5YWL14sSfrlL3+pBQsWaMmSJdqwYYN69uypq6++uvp548eP1+rVq1VeXq4TTjhB\nF110Ua39mpm2bNmi0047TWPGjNHs2bP12Wef6YwzztDFF1+s8vJyzZ8/X1dddZVWrlxZ/bynnnpK\nt956q3bv3q3Ro0dn/iJHjIYZAAAgYpdcconmzp2rl156SSNGjNDAgQNrPT527Fj90z/9kyTpa1/7\nmi644AK9+uqrB+1n06ZNeuGFF3TPPfeoY8eO6tOnj37yk59o/vz59R43Ly9P7733njZv3qxOnTrp\nG9/4Rq3Hi4uL1bFjR40cOVKXXnpp9dCRX/3qV7rzzjs1YMAAHXbYYSouLtYzzzxTPUxi6tSp6ty5\nc/Vjb7/9tnbt2lW9348//lhFRUWaOHGibr/9dknS888/ryFDhmjKlClq166djj/+eP3bv/2bnn76\n6ernTZgwQSeffLIkKT8/v0mvcZQYkoGcVTW2i0/eA0B8qoYgZH14Q0kwfndS66/xZqZLLrlEY8aM\n0dq1azV58uSDfre9/vrruuGGG7RixQrt27dPX3zxxUFngyXpgw8+UEVFhfr371+9rrKyUgUFBfUe\ne86cObrtttt07LHHasiQISouLtb48eOrHx88eHD17YKCAi1fvrz6OOeee26tYR4dOnTQpk2b1Ldv\nX91888165plnVF5eXr3N5s2b1bVrV7m7Fi5cqK5du+qKK66olf3111+vHtYhSfv379fkyZOrX6dB\ngwY1/oLGgIYZAAAgYgUFBRo6dKheeOEFPfLIIwc9PmnSJF177bVatGiR8vLyNG3aNG3evPmg7QYP\nHlw960YmU8wNGzZMJSUlkqRnn31W3/ve97R169bqx9evX6/hw4dX3646811QUKBHH320+mxvTU88\n8YQWLFigxYsX68gjj9T27dvVq1ev6j8CzEyXX365tm3bprPOOksvvviiOnXqpIKCAo0dO1alpaUZ\nvGK5hSEZAAAALWDOnDl6+eWX1bFjx4Me2717t3r27Km8vDwtXbpUJSUl9c6i0b9/f40bN07Tp0/X\nrl27VFlZqTVr1mjJkiX1HnPevHkqLy+XJHXv3l1mVqvRvvPOO7V3716tWLFCjz32mCZOnChJuvLK\nK3XTTTdVT3tXXl6uBQsWVGfNz89Xr1699Nlnn+mmm26qdcyqxvnee+/V8OHDdfbZZ+vzzz/X+PHj\n9e6772revHmqqKhQRUWF3njjDa1atarW83IRZ5gBAECr1JIXMcnE0KFDa92v2RDfd999mjFjhq65\n5hqNHTtWEydO1Pbt2+vddu7cubrhhhs0YsQI7dq1S0OHDtUNN9xQ7zEXLVqkGTNmaM+ePSosLNT8\n+fNrjQ0eO3ashg0bpsrKSl1//fU6/fTTJaU/LOjuGjdunD755BP17dtXF1xwgc455xxNnjxZixYt\n0sCBA9W7d2/dfvvteuCBB2plrcr74IMPaurUqZowYYKee+45lZaWavr06Zo+fboqKyt1/PHH6+67\n7z7oebmGK/0hZzGGGUnFlf7QmiRhDHNrv9JfFJJ29cBs4Ep/AAAAQEQ4wwwAWcYZZqBlNXTmEA1b\nt26djjrqKFVUVHCGOYOaTcMMAFlGwwy0LBpmZIIhGQAAAEBEaJgBAACAEDTMAAAAQAgaZgAAACAE\nDTNyVi5PYA4AbUXKUtHMbVxiX87FDOQ4rvQHAABaj5ZswrNw0ZWoJe0CJVOnTtXgwYN1xx13xB2l\nltx/5QAAABKqsLBQnTp1UteuXdW1a1d169ZNGzdujDtWzsrVd5dpmAEAACJiZnr++ee1a9cu7dq1\nSzt37tQRRxwRd6yclo05tSsrK7OQ5Es0zAAAoHWa5NlfsmTHjh267LLLNGDAAA0aNEi33nprdZP3\n2GOPafTo0Zo+fbp69uypYcOG6Y9//KMeffRRFRQUqF+/fpo7d271vhYuXKhRo0ape/fuKigoUCrV\n8JjzsOPWtXTpUp144onq3r27jjjiCM2YMUNSephHu3bt9NBDD2ngwIEaMGCA7rrrrurnubt+/vOf\na9iwYfrKV76iiRMnatu2bdWPn3feeerfv7969OihsWPH6u9//3u9x9+1a5dOOeUU/eQnP5EkrVq1\nSmeccYZ69+6tY445Rk8//XT1tlOnTtWPfvQjnXXWWerSpYvKysoa+RdoGhpmAACACNV3xnTq1KnK\ny8vTmjVrtGzZMpWWlurhhx+ufnzp0qU67rjjtHXrVl144YU6//zz9dZbb2nNmjWaN2+errnmGu3Z\ns0eS1KVLF82bN087duzQwoULdf/99+u5556rN0tjx63puuuu07Rp07Rjxw69//77Ov/882s9XlZW\nptWrV6u0tFSzZs3S4sWLJUm//OUvtWDBAi1ZskQbNmxQz549dfXVV1c/b/z48Vq9erXKy8t1wgkn\n6KKLLqq1XzPTli1bdNppp2nMmDGaPXu2PvvsM51xxhm6+OKLVV5ervnz5+uqq67SypUrq5/31FNP\n6dZbb9Xu3bs1evTosH+SpnP3nFvSsQAgmYIaFnstbcmFuo041fr5e1JfLlFo4v6PPPJI79Kli/fo\n0cN79Ojh5557rm/cuNHz8/N979691duVlJT4Kaec4u7ujz76qB999NHVj73zzjtuZv7pp59Wr+vd\nu7e//fbb9R7zuuuu82nTprm7+9q1a93M/MCBA40et67vfOc7Xlxc7OXl5bXWV+3zH//4R/W6n/3s\nZ37ZZZe5u/sxxxzjixcvrn7sk08+8cMOO8wPHDhw0DG2bdvmZuY7d+50d/epU6f697//fR85cqT/\nx3/8R/V28+fP9zFjxtR67g9/+ENPpVLu7j5lyhSfMmVKvd9HlYbqVCY1m1kyAAAAImJmeu6553Tq\nqadWr1u6dKkqKirUv3//6nWVlZUqKCiovt+vX7/q2x07dpQk9enTp9a63bt3S5Jef/113XDDDVqx\nYoX27dunL7744qCzwZL0wQcfNHrcmubMmaPbbrtNxx57rIYMGaLi4mKNHz+++vHBgwdX3y4oKNDy\n5curj3PuuefWmpWjQ4cO2rRpk/r27aubb75ZzzzzjMrLy6u32bx5s7p27Sp318KFC9W1a1ddccUV\ntbK//vrr6tmzZ/W6/fv3a/LkydWv86BBg+r9PrKBhhkAAKAFDR48WPn5+dqyZUtWpnqbNGmSrr32\nWi1atEh5eXmaNm2aNm/efMjHHTZsmEpKSiRJzz77rL73ve9p69at1Y+vX79ew4cPr749cOBASenm\n+dFHH9XJJ5980D6feOIJLViwQIsXL9aRRx6p7du3q1evXtXDVsxMl19+ubZt26azzjpLL774ojp1\n6qSCggKNHTtWpaWlTX+BsoAxzAAAAC2of//+GjdunKZPn65du3apsrJSa9as0ZIlS5q1v927d6tn\nz57Ky8vT0qVLVVJSUu/UbE097rx581ReXi5J6t69u8ysVqN95513au/evVqxYoUee+wxTZw4UZJ0\n5ZVX6qabbtL69eslSeXl5VqwYEF11vz8fPXq1UufffaZbrrpplrHrGqc7733Xg0fPlxnn322Pv/8\nc40fP17vvvuu5s2bp4qKClVUVOiNN97QqlWraj0vKjTMAACgdaq6mmA2lyyZO3eu9u3bpxEjRqhX\nr14677zzqudnrm8u4rC5ie+77z7ddttt6tatm+64447qxrW+54Ydt65FixZp5MiR6tq1q6ZNm6b5\n8+crPz+/+vGxY8dq2LBhOv3003X99dfr9NNPl5T+sOA555yjcePGqVu3bjr55JO1dOlSSdLkyZN1\n5JFHauDAgRo5cqROPvnkWvlqfu8PPvigBg0apAkTJigvL0+lpaWaP3++Bg4cqP79++vGG2/Uvn37\nGnzNssmi7sibw8w8F3MBQCbMTO6eezPvR4i6jTgF/+fSd7jSX+SSdvXAKrV+Tg5eH/qDk5zvEm1O\nrl7tBwDakpSllLKG5/VttiyfsQWixIf+AABA69FGz/q2tLZ2QouGGQAAABkrLCzUgQMH4o7RohiS\nAQAAAISgYQYAAABC0DADAAAAIZhWDgCyjGnlgJbV0HRhQE2HMq0cH/oDAACJ19ZmbUDLomEGAACJ\nxtllRI0xzAAAAEAIGmYAAAAgBEMykJOuvfZaPfHEE3HHaLKCggK9/fbbcccAAABZRMOMnLRnzx5t\n37497hhN1qNHj7gjAEBWpSwlSSr24uzuuCT4kB6XskYCMCQDAAAACMEZZuS82bNna/LkyXHHaNAH\nH3ygUaNGxR0DAABEhIYZOa9z587q2bNn3DEalMShIwAAIHMMyQAAAABC0DADAAAAIRiSAQAAGpT1\n2TGqMDsGEoQzzAAAAEAIGmYAAAAgRGQNs5kNNrNXzGyFmf3NzK4N1vcys5fM7F0zKzUzrvQAADGj\nZgNAw6I8w1whaZq7/5Okb0q62syOlXSDpJfc/auSFgf3AQDxomYDQAMia5jdfaO7/zW4vVvSSkkD\nJZ0j6fFgs8clTYgqAwAgM9RsAGhYi4xhNrNCSaMkvS6pn7tvCh7aJKlfS2QAAGSGmo2aUpZSylLZ\n33GJpRcgASJvmM2si6RnJV3n7rtqPubuLol5ZQAgR1CzAeBgkc7DbGaHKV14n3D33wWrN5nZEe6+\n0cz6S/q0vufOnDmz+nZRUZGKioqijAoAzVZWVqaysrK4YxyyQ6nZEnUbQDI0p2ZH1jCbmUmaI+nv\n7j67xkMLJE2RNCv4+rt6nl6r8AJALqvbHKZSEbx9HbFDrdkSdRtAMjSnZkd5hnm0pIslvWNmy4J1\nN0r6uaTfmNllktZJOj/CDACAzFCzAaABkTXM7v6aGh4jfXpUxwUANB01GwAaFukYZgAAkGzFXhzN\njifx+VEkB5fGBgAAAELQMAMAAAAhaJgBAACAEDTMAAAAQAgaZgAAACAEDTMAAGhQylJKWQQX4ymx\n9AIkAA0zAAAAEIKGGQAAAAhBwwwAAACEoGEGAAAAQtAwAwAAACE6xB0AAADkrmIvjmbHkzya/QIR\n4AwzAAAAEIKGGQAAAAhBwwwAAACEoGEGAAAAQtAwAwAAACFomAEAQINSllLKUtnfcYmlFyABaJgB\nAACAEDTMAAAAQAgaZgAAACAEDTMAAAAQgoYZAAAACNEh7gAAACB3FXtxNDue5NHsF4gAZ5gBAACA\nEDTMAAAAQAgaZgAAACAEDTMAAAAQgoYZAAAACEHDDAAAGpSylFKWyv6OSyy9AAlAwwwAAACEoGEG\nAAAAQtAwAwAAACFomAEAAIAQNMwAAABAiA5xBwAAALmr2Iuj2fEkj2a/QAQ4wwwAAACEoGEGAAAA\nQtAwAwAAACFomAEAAIAQNMwAAABACBpmAADQoJSllLJU9ndcYukFSAAaZgAAACAEDTMAAAAQgoYZ\nAAAACEHDDAAAAISgYQYAAABCdIg7AAAAyF3FXhzNjid5NPsFIsAZZgAAACAEDTMAAAAQgoYZAAAA\nCEHDDAAAAISgYQYAAABC0DADAIAGpSyllKWyv+MSSy9AAtAwAwAAACFomAEAAIAQNMwAAABAiEgb\nZjN7xMw2mdnyGutmmtlHZrYsWM6MMgMAIDPUbACoX9RnmB+VVLe4uqS73X1UsLwYcQYAQGao2QBQ\njw5R7tzd/2BmhfU8xMdiASDHULNRn2IvjmbHkzya/QIRiGsM84/N7G0zm2NmPWLKAADIDDUbQJsW\nR8N8v6Qhko6XtEHSXTFkAABkhpoNoM2LdEhGfdz906rbZvawpN/Xt93MmTOrbxcVFamoqCjqaADQ\nLGVlZSorK4s7RiQyrdkSdRtAMjSnZpt7tGOIgvFwv3f3rwX3+7v7huD2NEknufukOs/xqHMht/3g\nBz/QnDlzJEkPPfSQfvCDH8ScqGFr167V0KFDJUmFhYVau3ZtzIkQNzOTuydy3G9zanbwGHUbQCJl\nUrMjPcNsZk9JGivpK2b2oaRiSUVmdrzSn7xeK+mKKDMAADJDzQaA+kU9S8aF9ax+JMpjAgCah5qN\n+qQsJSmC2TJKghN6zJaBBOBKfwAAAEAIGmYAAAAgBA0zAAAAEIKGGQAAAAhBwwwAAACEaPELlwAA\ngOTI+uwYVZgdAwnCGWYAAAAgBA0zAAAAEIKGGQAAAAhBwwwAAACEoGEGAAAAQtAwAwCABqUspZSl\nsr/jEksvQALQMAMAAAAhaJgBAACAEDTMAAAAQAgaZgAAACBEow2zmX27nnWjo4kDAAAA5JYOGWzz\n/0oaVWfdvfWsAwAArUyxF0ez40kezX6BCDTYMJvZyZK+JamPmU2XVDX3S1cxlAMAAABtRNgZ5jyl\nm+P2wdcqOyV9L8pQAAAAQK5osGF291clvWpmj7n7upaLBAAAAOSOTMYw55vZQ5IKa2zv7n5qZKkS\n4he/+IV+/etfxx2jyXr16qXS0tK4Y7RKH3/8sU488cS4Y2Tszjvv1Jlnnhl3DAAJ9+7z76psZlnc\nMZpl0sJJ6tKvS9wxkOMyaZiflnS/pIclHQjWMVJf0ocffqg333wz7hhN1rdv37gjtFoVFRWJ+pnY\nunVr3BEAtAJ7tuzRhjc3xB2jWSorKuOOgATIpGGucPf7I08CAAByTspSkrI/W0bxkzPT+79oZlb3\nC0Qhk4b592Z2taT/LemLqpXuzqmpGmbMmKELLrgg7hgN2rJlC2+9R2TAgAF644034o6RsZtvvpkh\nOQAic/RZR6soVdT4hu/OlCRd/sblkeapz5NnPak95Xta/LhIrkwa5qlKD8H4aZ31Q7KeJsEKCgpy\neuzqpk2b4o7QauXn5+f0v31dvXv3jjsCgFas01c6acCJAxrf8N30l4y2zbL2ee1b/JhItkYbZncv\nbIEcAAAAQE5qtGE2symq50N+7j43kkQAAABADslkSMZJ+rJh7ijpVElvSaJhBgAAQKuXyZCMa2re\nN7MekpI3+TAAAGiybM+OUW0SM9QiOdo14zl7xAf+AAAA0EZkMob59zXutpM0QtJvIksEAAAA5JBM\nxjDfFXx1SfslrXf3D6OLBAAAAOSORodkuHuZpFWSuknqqRoXLwEAAABau0yGZJwv6ReSXg1W3Wtm\n17v705EmQ2T27NmjBx54IO4YoVauXBl3BAAAAEmZDcm4RdJJ7v6pJJlZH0mLJdEwJ9Tu3bt15ZVX\nxh0DAJAAKUtJimC2jBJLf2W2DCRAJrNkmKTyGve3BOsAAACAVi+TM8wvSlpkZiVKN8oTJb0QaSpk\nXceOHfXDH/4w7hhN8uCDD0qSjjnmmJiTAACAtqzBhtnMjpbUz92vN7P/IWl08NAfJZW0RDhkT7du\n3XJ+3HJdVQ3zt7/97ZiTAACAtixsSMZsSTslyd2fdffp7j5d0u8k3dMS4QAAAIC4hTXM/dz9nbor\ng3Vc6Q8AAABtQtgY5h4hjx2e7SBAXe58choA4pb12TGqMDsGEiTsDPNfzOygT4mZ2eWS3owuEgAA\nAJA7ws4w/0TSb83sIn3ZIP+LpHxJ50YdDAAAAMgFDTbM7r7RzL4l6RRJIyW5pOfd/eWWCgcAAADE\nLXQeZk8PIn05WAAAAIA2J5Mr/QEAAABtViZX+gNiYZa+AjuzZUSnsrJSBw4ciDtGxtq3bx93BKDN\nSVlKUgSzZZSkazyzZSAJaJiBNuySSy7RJZdcEneMjIwYMUIrVqyIOwYAoA1iSAYAAAAQgjPMQBtj\nZmrXLhl/K7s7Q3IAALFLxm9NAFnz5JNP6sCBA4lYli9fHvfLBQAADTMAAAAQhiEZyFm8FQ8A8cv6\n7BhVmB0DCcIZZgAAACAEDTMAAAAQgoYZAAAACEHDDAAAAISItGE2s0fMbJOZLa+xrpeZvWRm75pZ\nqZn1iDIDACAz1GwAqF/UZ5gflXRmnXU3SHrJ3b8qaXFwHziImcnM4o4BtCXUbBwkZSmlLJX9HZdY\negESINKG2d3/IGlbndXnSHo8uP24pAlRZgAAZIaaDQD1i2MMcz933xTc3iSpXwwZAACZoWYDaPNi\n/dCfp69MwczlAJAA1GwAbVUcV/rbZGZHuPtGM+sv6dP6Npo5c2b17aKiIhUVFbVMOgBoorKyMpWV\nlcUdIyoZ1WyJug0gGZpTs+NomBdImiJpVvD1d/VtVLPwAkAuq9scplIRfEAqPhnVbIm6DSAZmlOz\nI22YzewpSWMlfcXMPpR0m6SfS/qNmV0maZ2k86PMgORKv/sLoKVQs1GfYi+OZseTqPFIjkgbZne/\nsIGHTo/yuACApqNmA0D9uNIfAAAAEIKGGQAAAAhBwwwAAACEoGEGAAAAQtAwI2eZmcws7hgA0Kal\nLKWURTBVYomlFyABaJgBAACAEDTMAAAAQAgaZgAAACAEDTMAAAAQgoYZAAAACBHppbGBQ+HucUcA\ngDav2Iuj2fEkajySgzPMAAAAQAgaZgAAACAEDTMAAAAQgoYZAAAACEHDDAAAAISgYUbOMjOZWdwx\nAKBNS1lKKUtlf8clll6ABKBhBgAAAELQMAMAAAAhaJgBAACAEDTMAAAAQAgujQ0AaBNWv7hai6Yv\nijtGs1z0fy5Sj8IeccdoUHd7RVf93/dKkvK7Hy49f1PmT35+RESpGjZlxlZV7j8gSXr67A76/Ivc\nfW2HjfiLThrzoiQpv1u+OvfpFHOiJhj3Rykvd1/bpqBhRs5y97gjAGhFvtj1hTav3Bx3jGY5sO9A\nbMcu9uJGt2mnXeozsMZru7MJB9i5sumhDlHvfl/e3rF+i3Zt3d/iGTJVOKBcvfps/HJFU17buHll\n3AmyhiEZAAAAQAjOMAMA2pyjxh2lf73nX+OOEerJ7z6pHet3xB2jyTZ8epL6X/pY3DFCVS4qUrv9\n5ZKkKYunqDJvQMyJGrb1xX+vvv1J+bc1YOoDMabJQOnJUkWSToNnhoYZANDm5HfLV58RfeKOEap9\nfvu4IzRLRUVnqXvLj0tuinaH5UnBKIzeX+0tdcrdn4XP/pRffXvf/q45/9rKkvlz2xiGZAAAAAAh\naJgBAACAEDTMyFlmJjOLOwYAtGkpSyllqezvuMTSC5AANMwAAABACBpmAAAAIAQNMwAAABCChhkA\nAAAIQcMMAAAAhODCJchZ7h53BADICcseWabOfTvHcuxxd42TJP3p7j81uE2HDeU6clQTdzyJGo/k\noGEGACDH/des/4o7Qqjjxnykk5raMAMJwpAMAAAAIARnmAEAyEHHX3q89pTviTtGRgb03lV9u+vA\nrjEmAaJBwwwAQA4ac+OYuCNk7v2N0p/TN3sO6RlvFiACDMkAAAAAQtAwI2eZmcws7hgA0KalLKWU\npbK/4xJLL0AC0DADAAAAIWiYAQAAgBA0zAAAAEAIGmYAAAAgBA0zAAAAEIJ5mJGz3D3uCADQ5hV7\ncTQ7nkTDo6pZAAASmElEQVSNR3JwhhkAAAAIQcMMAAAAhKBhBgAAAELQMAMAAAAhaJgBAACAEDTM\nyFlmJjOLOwYAtGkpSyllqezvuMTSC5AANMwAAABACBpmAAAAIAQNMwAAABCChhkAAAAIQcMMAAAA\nhOgQ14HNbJ2knZIOSKpw96/HlQW5yd3jjgAgQM1uu4q9OJodT6LGIzlia5gluaQid98aYwYAQGao\n2QDarLiHZDABIwAkBzUbQJsU9xnm/zSzA5IecPeHYswCAAiX+JrdsfKv+h8/flqS1H1QN+m1Z2JO\n1IrsXht3guZb+iOpQ6e4UzSoT947cUeA4m2YR7v7BjPrI+klM1vl7n+IMQ8AoGGJr9mHaZNGfnPF\nlyvW/ym+MMgdnzwfd4JQndvHnQBSjA2zu28Ivpab2W8lfV1SdfGdOXNm9bZFRUUqKipq4YQAkJmy\nsjKVlZXFHSNSjdVsiboNIBmaU7NjaZjNrJOk9u6+y8w6SxonqdaF6msWXrRNZunhksyWgVxXtzlM\npVINb5xAmdRsKVl1e8vWo9X77DvijpEIqcJVkqTidcdk9oTOhZltVxIMiY9jtoyT7pMO7G354zbD\nmpfWaNmcZZKkft/+lgrjjdMqNKdmx3WGuZ+k3wYNUQdJT7p7aUxZAADhWl3N3rO3t3ofOTHuGAkR\nNBOt6fUadE7cCTK2teINrfhzhSTp8OOPjjlN2xVLw+zuayUdH8exAQBNQ80G0NbFPa0cAAAAkNNo\nmAEAAIAQNMwAAABAiDjnYQZCMTsGAMSv2Iuj2XEcs2MAzcQZZgAAACAEDTMAAAAQgoYZAAAACEHD\nDAAAAISgYQYAAABC0DAjZ5mZgkvxAgBikrKUUpbK/o5LLL0ACUDDDAAAAISgYQYAAABC0DADAAAA\nIWiYAQAAgBA0zAAAAECIDnEHABri7nFHAIA2r9iLo9nxJGo8koMzzAAAAEAIGmYAAAAgBA0zAAAA\nEIKGGQAAAAhBwwwAAACEYJYM5Cwzk8RsGQAQp5SlJEUwW0ZJusYzW0bmls9brtUvrI47Rqgrb/lC\nh3dK3/58x+c6vG+8ebKFhhkAACAB9u3ep32798UdI5RXVta43Xr+GGJIBgAAABCCM8wAAAA56p8v\n+mcd/d2j446RuZdmxZ0gEjTMAAAAOSq/W77yu+XHHSNje+MOEBGGZAAAAAAhOMOMnMXsGAAQv6zP\njlGF2TGQIJxhBgAAAELQMAMAAAAhaJgBAACAEDTMAAAAQAgaZgAAACAEDTNylpnJzOKOAQBtWspS\nSlkq+zsusfQCJAANMwAAABCChhkAAAAIQcMMAAAAhKBhBgAAAELQMAMAAAAhOsQdAGiIu8cdAQDa\nvGIvjmbHk6jxSA7OMAMAAAAhaJgBAACAEDTMAAAAQAgaZgAAACAEDTMAAAAQgoYZOcvMZGZxxwCA\nNi1lKaUslf0dl1h6ARKAhhkAAAAIwTzMABKhoqJCGzdujDsGACBDtu9TaW/7uGNkBQ0zgER47733\n1L9//7hjAAAydPhrJ8QdIWsYkgEAAACE4AwzgJzVoUMH9evXL+4YGdu2bZv27dsXdwwAiM2eXV10\noCJ9PrZz386ydjn8wc79n0n7d2e0KQ0zcpa7xx0BMRs+fHiixi1PmDBBzz33XNwxgKwq9uJodjyJ\nGt8aPXznDH2+7XNJ0s+2/Ewde3WMOVGIlXdJy36a0aYMyQAAAABC0DADAAAAIWiYAQAAgBA0zAAA\nAEAIGmYAAAAgBA0zcpaZySyHp6MBgDYgZSmlLJX9HZdYegESIJaG2czONLNVZvaemf3PODIAADJH\n3QbQlrV4w2xm7SXdK+lMSSMkXWhmx7Z0jraqrKws7gitVtJe2yTlTVLW1qi11e2yv8edoGmS9POf\npKxSsvImKaskrdXauCNkVRxnmL8uabW7r3P3CknzJf33GHK0SUn7D5ckSXttk5Q3SVlbqVZVt8tW\nxp2gaZL085+krFKy8iYpqySt07q4I2RVHFf6Gyjpwxr3P5L0jbobnX322S0WqLn+9re/xR0BAFpC\nRnX741+c0GKBmuPwDluk5FxpHUi8Zy98Vu3z2scdo0HHHLNMo0Zltm0cDXNG18J8/vnno84BAFl1\n1FFH6YQTTtBbb70Vd5Rsy6huDxy4LOocABJkTemauCOEOnz7HvXr0V/Shka3NfeWvZa7mX1T0kx3\nPzO4f6OkSnefVWMbLjAPINHcvdV8/J+6DaC1a6xmx9Ewd5D0D0mnSfpE0lJJF7p7wkaVAUDbQN0G\n0Na1+JAMd99vZtdIWiSpvaQ5FF0AyF3UbQBtXYufYQYAAACSJOeu9JekyfHN7BEz22Rmy+PO0hgz\nG2xmr5jZCjP7m5ldG3emMGZ2uJm9bmZ/NbO/m9n/FXemxphZezNbZma/jztLY8xsnZm9E+RdGnee\nMGbWw8yeMbOVwc/CN+POVB8zGx68nlXLjlz/f5YN1OxoULOjR82ORlJqttS0up1TZ5iDyfH/Iel0\nSR9LekM5PE7OzMZI2i1prrt/Le48YczsCElHuPtfzayLpDclTcjV11aSzKyTu+8Jxk++Jumn7v5a\n3LkaYmbTJf2LpK7ufk7cecKY2VpJ/+LuW+PO0hgze1zSq+7+SPCz0Nndd8SdK4yZtVO6hn3d3T9s\nbPukomZHh5odPWp2NJJYs6XG63aunWFO1OT47v4HSdvizpEJd9/o7n8Nbu+WtFLSgHhThXP3PcHN\nPKXHTeZsoTCzQZLOkvSwpKTMjpDzOc2su6Qx7v6IlB5Lm4TCq3QDuaY1N8sBanZEqNnRomZHI8E1\nW2qkbudaw1zf5PgDY8rSaplZoaRRkl6PN0k4M2tnZn+VtEnSK+6eyxe0vUfS9ZIq4w6SIZf0n2b2\nFzO7PO4wIYZIKjezR83sLTN7yMw6xR0qAxdIKok7RAugZrcAanYkqNnRSGrNlhqp27nWMOfO+JBW\nKnhr7xlJ1wVnLXKWu1e6+/GSBkn6jpkVxRypXmb23yR96u7LlIAzAIHR7j5K0nclXR28VZ2LOkg6\nQdJ97n6CpM8k3RBvpHBmlifpbElPx52lBVCzI0bNzj5qdqQSV7OlzOp2rjXMH0saXOP+YKXPWCAL\nzOwwSc9Kmufuv4s7T6aCt3MWSjox7iwN+Jakc4IxZk9JOtXM5sacKZS7bwi+lkv6rdJvreeijyR9\n5O5vBPefUboY57LvSnozeG1bO2p2hKjZkaFmRyeJNVvKoG7nWsP8F0lHm1lh0O1PlLQg5kytgpmZ\npDmS/u7us+PO0xgz+4qZ9Qhud5R0hqScvO6uu9/k7oPdfYjSb+m87O6T487VEDPrZGZdg9udJY2T\nlJOzBrj7RkkfmtlXg1WnS1oRY6RMXKj0L+G2gJodEWp2dKjZ0UlozZYyqNstfuGSMEmbHN/MnpI0\nVlJvM/tQ0m3u/mjMsRoyWtLFkt4xs6oidqO7vxhjpjD9JT0efGq1naQn3H1xzJkyletvU/eT9Nv0\n72N1kPSku5fGGynUjyU9GTRkayRdGnOeBgW/zE6XlMtjDLOGmh0panbLoWZnV2JqtpR53c6paeUA\nAACAXJNrQzIAAACAnELDDAAAAISgYQYAAABC0DADAAAAIWiYAQAAgBA0zAAAAEAIGuZWzMzGmdkL\nZrbZzPaa2T/M7OdVk8vHnG2qmR00N2OwvtLMCuLIlavMbIKZTTuE5x/0uprZOjN7IjsJARwqanbr\nQc1ufWiYWykzu0nSi5L2SLpM6SsD/UrSVElvmNmg+NJJQY7v17P+eUnflLSxRdPkvgmSpmd5n67c\nn7AfaBOo2a0ONbuVyakr/SE7zOwUSXdIusfdZ9R46A9m9ltJb0p6XNJpceQL4+6bJW2OO0cczCzf\n3b+IOweAlkXNTiZqdtvCGebW6WeStki6se4D7r5O0s8lnWJmJ0mSmRUGb/1MqbmtmRUF679TZ/2/\nmdmfzewzM9tmZr8xs8F1tplkZsvMbJeZ7TCzd8zsh8FjZZK+I2l0sP9KM3s5eKy+t6EOM7M7g7ej\nvjCztWZ2h5l1qLFN1ffwQzO73cw+CbItMLOBmWZriJnNDPY/0sxeCb73T8wsZcH1Smts28fMfmVm\nH5nZ52a20swur7NN1fc5xsyeNrNtkv7cwLEfkzRZ0sAar9fa4LF8M7vHzJYH38+G4HseHvb9NHCc\n9mb2YPCanNrU5wNoNmo2NZuaneM4w9zKBAVprKTfuvu+Bjb7vaRZkk6R9EaN9Y2+1WNmV0q6T9Ij\nkmZK6hZ8fdXM/tndd5vZtyU9Ien/kTRD6T/MjpXUPdjNjyTNC9ZfEazbGXLYxyWdJ+nfJb0mabSk\nmyUNlXRRnW1vlPRfSl+7vp+ku4JjnRLkbyxbY34naU6Q5UxJt0qqlJQK9t8tyJgvqVjS2mC7+y19\nNuLeOvt7UlKJpPvV8P/H2yV9RdJJks4O1lWd1ciX1FXS/5L0saSekq6W9CczO9bdN2XyTZlZR0lP\nSfqGpLHu/tdMngfg0FCzqdmiZieDu7O0okXpglMp6d9Dtjk82Obe4H5hcH9yne2KgvXfCe53kbRD\n0sN1titUuhhcF9z/qaQtjeQsk7SknvVTg2MWBPdHBvdvq7PdzcH6r9X5Hl6us92MYP0RmWZrIO/M\nYD8/q7P+QaV/cXQL7t8qaa+ko+rZrlxSuzrf510ZHv8xSR9msF07SZ2CTD9p6HUN1q2TNFfpgv2a\npPckDYn7Z5iFpS0t1GxqNjU7GQtDMtq2/U3c/mSl/zIuMbMOVYukjyT9Q+m37CRpqaSeZvaEmf03\nO7RPeFftc16d9fPqPF7l/9S5/7fga9XbhYea7Td17v9a6V9KI4P7Zyr9Nt26Oq9RqaTekkbUef5v\nm3j8g5jZ+Wb2evAW4X5Ju4NMX23kqS5poNKF93BJ33L3tYeaB0BkqNnUbGp2TGiYW58tSv+1XBiy\nTdVjHzVx332Dr/8paV+dZaSkXpLk7kuUfjtusKT/LelTM3vJzL7WxOOpap+SNtRZv6nO41W21rlf\n9TbY4VnKVvftsqr7VWPu+ir99mqFar8+v1G62PWu8/y631eTmNnZkuZLWiHpQklfV/ptwHIF33PY\n0yX9s9Jvb/7G3csPJQuAZqFm10bNDnm6qNmxYQxzK+Pu+81siaRx1vAneM8JvpYGXz8PvubV2a5u\nodgSfJ2i9H/2unbVyPGspGfNrJPSY9FmKT1l0sB6nhemqpj2l/R+jfVH1Hk8Yw1lM7NBHrzvFeII\npce4VekXfP04+LpZ6emVrmvg+e/WjZNx8PpdIOk9d6+e7snMDtPB/3b1cUkvSHpH0iwz+9zdf3mI\neQA0ATW7cdTsWsemZseEhrl1+oXSZxT+l9LjwaqZ2RBJ/1NSmbu/E6zepPRf9XX/Yh9f5/4flS6w\nR7t7RpOnu/seSQvN7ChJs82st7tvCY6XSYF4Nfh6QfD9VKn64EhZJjkyyab0mY8t4c/S+UoX6yoX\nKP2aLA/uvyjpx0qPXcvmX/9fSOpYz/pOkg7UWXeJmvDukbv/h5kdUPrfp527z25+TADNQM1uRjZR\ns6nZLYiGuRVy95fNrFhSyswKlf6E8TZJJ0i6QekPgUyqsb2b2a8lXWZm7yr9F/V4pd+mqrnfnWZ2\nvaT/z8z6KF1odih9BmKspFfc/Skzu13pt7leUfrtq0GSrpW0LCi8Uvpsx1Vmdr7SZyF2unvdv+Tl\n7ivM7ClJM4NxZX9SelzeLZJK3L2+syYNyjBbmB+YWTtJf5H0r0pfYKDY3avO1NwjaaLS86feo/Rr\n2VnSMZK+7e4TmpK3hhWSLg8+8f6mpM/dfbnSZxv+u5ndLWmhpBMlXSNpu9Jv34Wpftzd7wkK8D1B\nAb67mTkBNBE1u2HU7Fqo2XGK+1OHLNEtSheHF5V+C6wyWF6T1Lmebbsr/enbcqX/Yr9P0llK/yX8\nnTrbflfSy0oX3s+ULjAPSzomePys4LifKP3W4XpJDyn41HOwTT+li8VO1fiktNKfDD6g2p8MPkzp\nSf3XKT22bK3S0/a0r7FNYbCf79fJWlTze8gkWwOv5cxg/yOC731PsI9UPdv2kHS30r9UvlD6bNCr\nkq6tsU3V9zk0w3/LTkpPZVT1b/l+sN6C1+bj4N/iFUnHB6/RI/Ucr+brulbS3DrHuSrY7qdx//yy\nsLS1hZpNza7neNTsHFkseMHRBlj6GvQTJY1395fizpMkZjZT0m2SOrh7ZcxxALQB1Ozmo2Yj25gl\no235vtITxP/GzOpOlQMAyC3UbCBHMIa5DXH3CgVXT0KTuQ7909EAkDFq9iGhZiOrGJIBAAAAhGBI\nBgAAABCChhkAAAAIQcMMAAAAhKBhBgAAAELQMAMAAAAhaJgBAACAEP8/JZR0XLPmxiYAAAAASUVO\nRK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x10a82d150>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "There are 3.0 +- 1.4 total questions per talk\n",
      "There are 2.8 questions per talk when the speaker is male\n",
      "There are 3.3 questions per talk when the speaker is female \n",
      "There is a 41.2 percent chance that the questions are drawn from the same distribution for male and female speakers.\n"
     ]
    }
   ],
   "source": [
    "# How many questions were there per talk? Did the gender of the speaker affect it?\n",
    "\n",
    "fig2 = plt.figure(2,(12,6))\n",
    "ax4 = fig2.add_subplot(121)\n",
    "qpt = [len(x) for x in q['questions']]\n",
    "ax4.hist(qpt,bins=range(0,8),histtype='step',range=(0,8),linewidth=3, color='k')\n",
    "ylims4 = ax4.get_ylim()\n",
    "ax4.vlines(np.mean(qpt),ylims4[0],ylims4[1],linestyle='--',color='black')\n",
    "ax4.set_xlabel('Questions per talk',fontsize=16)\n",
    "ax4.set_ylabel('Count')\n",
    "\n",
    "ax5 = fig2.add_subplot(122)\n",
    "mq = [len(x[1]['questions']) for x in q.iterrows() if x[1]['speaker'] == 'M']\n",
    "fq = [len(x[1]['questions']) for x in q.iterrows() if x[1]['speaker'] == 'F']\n",
    "ax5.hist(mq,bins=range(0,8),histtype='step',range=(0,8),linewidth=3, color='purple',label='Male speaker')\n",
    "ax5.hist(fq,bins=range(0,8),histtype='step',range=(0,8),linewidth=3, color='orange',label='Female speaker')\n",
    "ax5.set_ylim(ax4.get_ylim())\n",
    "ylims5 = ax5.get_ylim()\n",
    "ax5.vlines(np.mean(mq),ylims5[0],ylims5[1],linestyle='--',color='purple')\n",
    "ax5.vlines(np.mean(fq),ylims5[0],ylims5[1],linestyle='--',color='orange')\n",
    "ax5.set_xlabel('Questions per talk',fontsize=16)\n",
    "ax5.legend(loc='upper right')\n",
    "\n",
    "plt.show()\n",
    "\n",
    "# Test to see if the distribution is different for male vs. female speakers\n",
    "\n",
    "from scipy.stats import ks_2samp\n",
    "\n",
    "D,p = ks_2samp(mq,fq)\n",
    "print 'There are %.1f +- %.1f total questions per talk' % (np.mean(qpt),np.std(qpt))\n",
    "print 'There are %.1f questions per talk when the speaker is male' % np.mean(mq)\n",
    "print 'There are %.1f questions per talk when the speaker is female ' % np.mean(fq)\n",
    "print 'There is a %.1f percent chance that the questions are drawn from the same distribution for male and female speakers.' % (p*100)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Engagement with the speakers was nicely high at the <span style=\"color:red\">AGN2015</span> conference, with an average of 3 questions per talk being asked. This was higher than the averages for both the NAM $(2.2\\pm1.1)$ and AAS $(2.8\\pm1.5)$ meetings. The number of questions per talk was slightly higher when the speaker was female instead of male, with 3.3 compared to 2.8. However, a two-sided Kolmogorov-Smirnov test gives a 41% chance that the two samples are drawn from the same distribution, with no indication of a statistically significant difference. \n",
    "\n",
    "A difference in the number of questions per talk depending on the speaker might be interpreted as either a positive or negative effect (preferentially ignoring speakers of certain genders, overaggressively questioning/harassing them, paying attention and engaging with them, etc). More analysis on this data set would be pure speculation from me as to which was responsible at this conference. "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Did the gender of people of asking questions depend on the speaker's gender?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "                                 M         F\n",
      "Male asks 1st question    0.694444  0.305556\n",
      "Female asks 1st question  0.586538  0.413462\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAeYAAAGDCAYAAAD+qrMmAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xu8XFV99/HPl4uCVJRLVURDiLcWioIo4gUMaq1WrQpe\nWgVEfarSavWpVlvtY0K1VqtPvT9oLaAooLViQQS8kioIiAIhKIhoUlFQEEFQLoXwe/7Y+5DJZM7J\nOUnOmZXweec1rzmz99p7r1mzM99Za+89k6pCkiS1YbNxV0CSJK1iMEuS1BCDWZKkhhjMkiQ1xGCW\nJKkhBrMkSQ3ZYtwVUPuSeE2dJK2DqspMl7HHrGmpKm9VLFq0aOx1aOVmW9gWtsXUt3VlMEuS1BCD\nWZKkhhjM0gwsXLhw3FVohm2xim2xim2x/rI+4+C6a0hS7ieSNDNJKE/+kiRp42YwS5LUEINZkqSG\nGMySJDXEYJYkqSEGsyRJDTGYJUlqiMEsSVJDDGZJkhpiMEuS1BCDWZKkhhjMkiQ1xGCWJKkhBrMk\nSQ0xmCVJaojBLElSQwxmSZIaYjBLktQQg1mSpIYYzJIkNcRgliSpIQazJEkNMZglSWqIwSxJUkMM\nZkmSGmIwS5LUEINZkqSGGMySJDXEYJYkqSEGsyRJDTGYJUlqiMEsSVJDDGZJkhpiMEuS1BCDWZKk\nhhjMkiQ1ZItxV0AbhyNyxLirIDVh0XGLx10FbeLsMUuS1BCDWZKkhhjMkiQ1xGCWJKkhBrMkSQ0x\nmCVJaojBLElSQwxmSZIaYjBLktQQg1mSpIYYzJIkNcRgliSpIQazJEkNMZglSWqIwSxJUkMMZkmS\nGmIwS5LUEINZkqSGGMySJDXEYJYkqSEGsyRJDTGYJUlqiMEsSVJDDGZJkhpiMEuS1BCDWZKkhhjM\nkiQ1xGCWJKkhBrMkSQ0xmCVJaojBLElSQwxmSZIaYjBLktQQg1mSpIYYzJIkNcRgliSpIWsN5iQr\nk1wwcJs3W5VJsiLJ9rO07sVJXr8Oy22f5IwkNyb54DTKvyTJTutWy/WT5NlJfn/g8RFJnjyOukiS\n1s0W0yhzU1XtNes16VSD674F+HvgD/rb2hwGXAxctY7bWx/PBb4AXAJQVYvGUAdJ0npYp6HsJHsn\nWZLkO0lOT3K/fvqSJP+S5LwklyR5dJLPJ7ksydsGlv98v+zFSf58km0cnOTcvpf+kSRr1DXJ/0ny\n7STLknx0YPpfJflekqVJjh9YpPr5f57k1CRbDZU9YXgbVXVTVZ0F3Dq07c2TfLzf9kVJXpfkIOBR\nwHFJzk+y1Yh2W5rkwiTvTrKsn37YYG88ySlJntj//dQk30ry3ST/nmSbfvo7B+r97iSPBZ4FvLvf\n9oK+fgf15Z/cT78oyVFJ7tZPX9GPJny3n/ewUa+HJGluTCeYtx4Yxv5cki2ADwIHVdWjgGOAf+zL\nFnBrVT0aOBI4CXgVXU/zsCTb9eVe1i/7aOCvBqYD0A/HvgB4XN9bvwN48Yi6faiq9qmqPfp6PrOf\n/iZgz6p6RL/9gVXn1cAfA8+uqluGyr5yinYY7nHvCdy/qvaoqocDR1fV54DvAC+qqkf26x90DPCX\nVbVnv77JevEFVJIdgbcAT66qvYHvAn/dD/c/p6p27+v9tqo6GzgZeEO/7R8PrGerftsv6Ou6BXD4\nwLau6dd/JPCGKdpAkjTLphPMN1fVXv3tIOD3gN2Brya5gC44dh4of3J/fzFwcVX9oqr+B/gx8MB+\n3muTXAic3U97yMDyAZ4M7A18p9/Gk4BdR9TtSUnOSXJRX2a3fvpFwPFJXgysHFjvocDTgOdV1W1T\nlJ2OHwELknwgyR8BNw49h9UkuTdwr6o6s5/0yVHlhtaxb/+cvtW3w6HAPODXwC19z/e5wM1TbDvA\nw4DlVXV5P+0TwP4DZU7s788H5k9RJ0nSLJvOMeZhAb5XVY+bZP7EkO8drD78ewewRZKFdMG7b1Xd\nkuQMYCvW9ImqevOkleh6gR8G9q6qnyVZBGzdz34GXfA8C3hLkj3oeobLgEfQfRhYMVnZqlprQFfV\n9UkeThf0r6Lr4b98Yvbalmf1AL2d1T8kDbbHV6rqRWssnOxD147PA17d/z3ZtoenZWjaxOu0kkn2\niTM4486/5zOfXUd+TpKku64l34cll6z/etYlmH8A/G6SfavqnCRbAg+pqu9PY9kA2wLX9aH8e3S9\nwkEFfA04Kcl7q+qafuj2d6rqJwPlJsLr2iS/Azwf+PckAeZV1ZIkZwF/CvxOX/YCuuHak/te7s9H\nlN0GuGGSuq96kOwA3FZVJya5DDi2n3Vj/xxXf1JdkF+f5PH9MevBofkVwOF93R8A7NO3wznAh5M8\nqKp+1B9fvj9wJbBNVZ2W5Ft0vffJtl10r9n8ifUAhwD/NeI5TuoADphJcUm6y1m4W3ebcMSJk5ed\nynSCebXeVlX9T5LnAR9Icq9+He8FhoN51DHUAk4HXpXk+3SBcfYaG6y6JMnfA1/uT/q6DfgL4CcD\nZa5P8jG6IfOfA+f2szYHPtnXLcD7q+rXXeZRVXVWkjcAXwSeOqLsGqGcZAVwT+BuSZ7dL3c34JiB\nk9L+tr//OPCRJDfRHSMfPM78UuDoJAV8eeC5nJlked+Gl9AdS6aqfpnkMOCEJHfvi7+FLoBP6kcN\nAvzvft6ngY8leQ3dB5WJ9d+a5KXAZ/tzBL4NfGRi9mDTM73eviRplqTK9+FxSLILcEp/4lrTktRi\nFo+7GlITFh23eNxV0EYiL4aqmupcopH85q/xGT7OK0nSOh1j1gZQVSuAh4+7HpKktthjliSpIQaz\nJEkNMZglSWqIwSxJUkMMZkmSGmIwS5LUEINZkqSGGMySJDXEYJYkqSEGsyRJDTGYJUlqiMEsSVJD\nDGZJkhpiMEuS1BCDWZKkhhjMkiQ1xGCWJKkhBrMkSQ0xmCVJaojBLElSQwxmSZIaYjBLktQQg1mS\npIYYzJIkNcRgliSpIQazJEkNMZglSWqIwSxJUkMMZkmSGmIwS5LUEINZkqSGGMySJDXEYJYkqSEG\nsyRJDTGYJUlqiMEsSVJDUlXjroMal6TcTyRpZpJQVZnpcvaYJUlqiMEsSVJDDGZJkhpiMEuS1BCD\nWZKkhhjMkiQ1xGCWJKkhBrMkSQ0xmCVJaojBLElSQwxmSZIaYjBLktQQg1mSpIYYzJIkNcRgliSp\nIQazJEkNMZglSWqIwSxJUkMMZkmSGmIwS5LUkC3GXQFJ0l1PknFXYYOqqg22LoNZkjQWGzLMxmlD\nf8hwKFuSpIYYzJIkNcRgliSpIQazJEkN8eQvSdLYHZEjZnX9i2rRtMrNnz+fq666iiuvvJIddtjh\nzul77bUXS5cuZcWKFcybN2+2qgnYY5Yk6U5JWLBgASeccMKd05YtW8bNN988Z5d4GcySJA04+OCD\nOfbYY+98/IlPfIJDDz10zi7vMpglSRqw7777csMNN3DppZeycuVKPvOZz3DwwQfP2fY9xixJ0pBD\nDjmEY489lv3335/ddtuNnXfeec62bTBLkjQgCYcccgj77bcfy5cvn9NhbHAoW5KkNcybN48FCxZw\n2mmnceCBB87ptu0xS5I0wlFHHcX111/P1ltvze233z5n2zWYJUljN93rjOfSggULVns8V5dLZVP5\ndQ/NniTlfiJpQ0qySf261Kjn0k+fcZp7jFmSpIYYzJIkNcRgliSpIQazJEkNMZglSWqIwSxJUkMM\nZkmSGmIwS5LUEINZkqSG+JWckqTxO36Wv+7yRdP7lrH58+dz9dVXs/nmmwPdt3dddtll3O9+95vN\n2q3GHrMkSb0knHLKKdx4443ceOON3HDDDXMaymAwS5LUFINZkqQB4/5xDYNZkqReVfGc5zyH7bbb\nju22244DDzxwzuvgyV+SJPWScNJJJ/GkJz1pbHWwxyxJUkMMZkmSGuJQtiRp/KZ5nfFdQcZ99pna\nl6TcTyRtSEnGfvbzhjLZc+mnz/ibUxzKliSpIQazJEkNMZglSWqIwSxJUkMMZkmSGmIwS5LUEK9j\nliSNRTLLv8G8kTKYJUlzblO5hnk2OJQtSVJDDGZJkhpiMEuS1BCDWZKkhhjMkiQ1xGCWJKkhBrMk\nSQ0xmCVJaojBLElSQwxmSZIaYjBLktQQg1mSpIYYzJIkNcRgliSpIQazJEkNMZglSWqIwSxJUkMM\nZkmSGmIwS5LUkC3GXQFtHI7IEeOugqQRFh23eNxV0AZmj1mSpIYYzJIkNcRgliSpIQazJEkNMZgl\nSWqIwSxJUkMMZkmSGmIwS5LUEINZkqSGGMySJDXEYJYkqSEGsyRJDTGYJUlqiMEsSVJDDGZJkhpi\nMEuS1BCDWZKkhhjMkiQ1xGCWJKkhBrMkSQ0xmCVJaojBLElSQwxmSZIaYjBLktQQg1mSpIYYzJIk\nNcRgliSpIQazJEkNMZglSWqIwSxJUkMMZkmSGmIwS5LUEINZkqSGGMySJDXEYJYkqSEGsyRJDZmV\nYE5yR5JPDjzeIsk1Sb6wluUWrq3MetZrRZLt12G55yf5XpKVSR65lrL3SnL4utdy/SR589Djs8ZV\nF0nSzM1Wj/m3wO5Jtuof/yHwU6BmaXvTta7bXwY8F/jGNMpuB/zFOm5nQ/i7wQdV9fhxVUSSNHOz\nOZR9KvCM/u8/A04AApBknyTfSnJ+krOSPHR44STbJDk6ybl9uT+ZpMxXk3w3yUUTZfrpX0xyYZJl\nSZ4/tNzWSU5L8vIk9xgq+4Lh7VTVpVV12Yjt797X74J++QcD7wQe1E9714hl3pLkB0m+meT4JK/v\npy9Jsnf/945Jlvd/b57k3Um+nWRpklf003dK8o1+O8uSPCHJO4Gt+2mf7Mv9pr9Pv55lfVu9oJ++\nsN/2Z5NckuRTI19NSdKc2GIW1/0Z4K1JTgH2AI4C9uvnXQLsV1UrkzwFeAfwvKHl3wJ8rapeluTe\nwLlJvlpVNw2UuQV4blXdmGRH4GzgZOBpwM+q6hkASe45sMw9+7p9oqo+leSgobLbzuA5vhJ4f1Ud\nn2QLuvZ8E7B7Ve01XLgP3hcCjwC2BM4HvtPPLkb36F8OXF9V+yS5O3Bmki8DBwKnV9U7kmwG3KOq\nzkzyl0Pbnljngf12Hw78LnBekokRgD2B3YCrgLOSPL6qHAKXpDGYtWCuqmVJ5tP1lr84NPvewLF9\nD7PoQmrYU4FnJXlD//juwAOBHwyU2Qz4pyT7AXcA909yH+Ai4D19D/KUqjqzLx/gJOBdVXVCP22y\nstNxNvCWJA8ATqyqy5NkivL79eVuAW5JcvI0tvFUYI8kEx9ctgUeDJwHHJ1kS+A/q2rpWtbzBOD4\nqirg6iT/BTwauAH4dlVdCZDkQmA+YDBL0hjMZo8Zut7re4An0vXSJryNrjf83CS7AEsmWf7Aqvrh\nFOt/MbAj8Mi+970c2KqqfphkL7qh9Lcn+VpVvY3uQ8CZwNPphtaZouxaVdUJSc4BngmcmuSVwPKp\nFqEfzu8N/n07qw4tbMXqXl1VXxleWf+B5JnAx5P8S1V9crjMFNuemAZw68C0lYzYL87gjDv/ns98\ndmXXKTYlSXc9S74PSy5Z//XM9uVSRwOLq+p7Q9O3Ba7s/37pJMt+CfiriQd9eA7bFri6D+UDgF36\nsjsBt1TVcXQfDAaXfStwXZIPT1J2yrOuGQi3JLtW1fKq+iBdT3wPuh7oPSdZ9hvAc5Js1Q+vP3Ng\n3grgUf3fg8P6XwL+oh8qJ8lD++Pi84Brqurf6A4TTDzH2ybKDvkm8MIkmyX5XWB/4NusGdYjHTDw\nz1CWpDUt3A0WH7Tqtq5mK5gLoKp+VlUfGpg20UP7Z7oh6POBzVn92OrE328DtuxPVLoYOGLEdo4D\nHpXkIuAQumPX0AXkuUkuoAvit69WuarX0p0k9a6hsv+n3+5qkjw3yRXAvsAXk5zWz3pBkov7ZXcH\njq2qX9Edp102fPJXVV1Ad3x7Kd3JceexKhjfAxzet8kOA+3wb8D3gfOTLAOOpOvRLgQu7Ms/H3h/\nX/5fgYuy6nK1idfi83TD9kuBrwF/U1VXM/rY9rjPnpeku6x0hxw1DkkWAb+pqv877rpMJUktZvG4\nqyFphEXHLR53FTSJvBiqalqjkoP85q/x85ORJOlOs33yl6ZQVaOG5yVJd2H2mCVJaojBLElSQwxm\nSZIaYjBLktQQg1mSpIYYzJIkNcRgliSpIQazJEkNMZglSWqIwSxJUkMMZkmSGmIwS5LUEINZkqSG\nGMySJDXEYJYkqSEGsyRJDTGYJUlqiMEsSVJDDGZJkhpiMEuS1BCDWZKkhhjMkiQ1xGCWJKkhBrMk\nSQ0xmCVJaojBLElSQwxmSZIaYjBLktQQg1mSpIYYzJIkNcRgliSpIQazJEkNMZglSWqIwSxJUkMM\nZkmSGmIwS5LUkFTVuOugxiUp9xNJmpkkVFVmupw9ZkmSGmIwS5LUEINZkqSGGMySJDXEYJYkqSEG\nsyRJDTGYJUlqiMEsSVJDDGZJkhpiMEuS1BCDWZKkhhjMkiQ1xGCWJKkhBrMkSQ0xmCVJaojBLElS\nQwxmSZIaYjBLktQQg1mSpIYYzJIkNcRgliSpIQazJEkNMZglSWqIwSxJUkMMZkmSGmIwS5LUEINZ\nkqSGGMySJDXEYJYkqSEGsyRJDTGYJUlqiMEsSVJDDGZJkhpiMEuS1BCDWZKkhhjMkiQ1xGCWJKkh\nBrM0A0uWLBl3FZphW6xiW6xiW6w/g1maAd90VrEtVrEtVrEt1p/BLElSQwxmSZIakqoadx3UuCTu\nJJK0DqoqM13GYJYkqSEOZUuS1BCDWZKkhhjMulOSpyW5NMkPk7xpkjIf6OcvTbLXXNdxrqytLZK8\nuG+Di5KcleTh46jnXJjOftGXe3SS25McOJf1myvT/P+xMMkFSS5OsmSOqzhnpvH/Y8ckpye5sG+L\nw8ZQzTmR5Ogkv0iybIoyM3vfrCpv3gA2By4H5gNbAhcCvz9U5o+BU/u/HwOcM+56j7EtHgvcq//7\naXflthgo93XgFOCgcdd7TPvEvYHvAQ/oH+847nqPsS0WA/800Q7AtcAW4677LLXHfsBewLJJ5s/4\nfdMesybsA1xeVSuq6jbg08Czh8r8CfAJgKo6F7h3kvvObTXnxFrboqrOrqpf9w/PBR4wx3WcK9PZ\nLwBeA/wHcM1cVm4OTacdXgR8rqp+ClBVv5zjOs6V6bTFVcC2/d/bAtdW1e1zWMc5U1XfBK6bosiM\n3zcNZk3YGbhi4PFP+2lrK7MpBtJ02mLQy4FTZ7VG47PWtkiyM90b85H9pE3xUo/p7BMPAbZPckaS\n7yQ5ZM5qN7em0xYfA3ZPciWwFHjtHNWtRTN+39xiVqujjcl030yHr8nbFN+Ep/2ckhwAvAx4/OxV\nZ6ym0xbvA/62qipJWHMf2RRMpx22BB4JPBm4B3B2knOq6oezWrO5N522eDNwYVUtTPIg4CtJHlFV\nN85y3Vo1o/dNg1kTfgY8cODxA+k+2U1V5gH9tE3NdNqC/oSvjwFPq6qphrI2ZtNpi72BT3eZzI7A\n05PcVlUnz00V58R02uEK4JdVdTNwc5JvAI8ANrVgnk5bPA74R4Cq+lGS5cDDgO/MSQ3bMuP3TYey\nNeE7wEOSzE9yN+CFwPAb68nAoQBJ9gWur6pfzG0158Ra2yLJPOBE4OCqunwMdZwra22LqlpQVbtW\n1a50x5kP38RCGab3/+Mk4AlJNk9yD7oTfb4/x/WcC9Npi0uBpwD0x1MfBvx4TmvZjhm/b9pjFgBV\ndXuSVwNfojvr8qiquiTJK/v5H62qU5P8cZLLgd8CLx1jlWfNdNoCeCuwHXBk31O8rar2GVedZ8s0\n22KTN83/H5cmOR24CLgD+FhVbXLBPM194h3AMUmW0nUA31hVvxpbpWdRkhOAJwI7JrkCWER3WGOd\n3zf9Sk5JkhriULYkSQ0xmCVJaojBLElSQwxmSZIaYjBLktQQg1mSpIYYzNKYJVmc5I4kTxx3XdZH\nkm37n7db0f/84x0b+89hDrw2+4+7Li3YVPbV1hnM2qj1bxJT3V7SQB0PW0tdauC2Mftn4NV0P1rw\nj3Q//TflNxz1Ib58Q1UgyZIkd2yo9d3V3IX21ab5zV/aFBRwxCTzLpjLiqzFZG9mHwJOYPVfoNkY\nPRP4QVWN+lnIqWzoN3lDY/1t6vtq0wxmbRKq6h/GXYdpGPmrS1V1Ld0PyW/s7g8sGXcltEFs6vtq\n0xzK1iav/7L9O5Ick+ShST6T5OokKyeOlSXZO8n7kyxNcm2Sm5NcluQ9Se49xbpfmORrSX7VL7M8\nyfFJ9u7nLwGO7osfMzTMPq8vM+lxzCRPTnJ6v/5bkvwgyT8l2XZE2SX9ejZP8uYkP+yX+UmSdybZ\ncobttlOSD/fDzbf2bfa5JI8ctd3+4cKB53fGFOte2C8zD5if1dvlmJm2wcRrDOzfPVxtfWcMlDsg\nyb8m+X6SXye5KcmyJG9NcveZtM8Uz+3BST6b5Lokv0lyVrrvSh45TDxVWyX5+OC+MjTvMUn+I8nP\n+9fnJ0k+kmSnEWUX9M/78v45X5vkoiRHJtm+L7OEhvfVJPsl+UKSn/Zlr0pydpK3jnwhNmL2mHVX\n8iDgHOAHwCeBrYFf9/P+HHgOXY/vy3QfWh8F/DXdzxg+pqp+M7GiJAGOofvVmGvoflXpGrqfd1tI\n9+s63+3LXAc8G/hP4MKB+vyaKaT7UYAjgRuBzwJXAwcAbwKeleTxVTVqHScATwBOBW4AngG8EbgP\n3W9Hr1WSXYEzgZ2ArwHH0YXo84FnJDmoqr7YFz8GOIPuy/tXAB/vp6+YYhPL6Q4/vK5//N6BeXe2\n0Qza4Lp+fYcBu9Ad354wWI830v3S0beALwBb0bXVYroPFU+pqnU+Rp3kIcDZwPZ07X8h8BC61/60\nvtioYeKpht/XmJfkZcC/AjfT/XrRFcBDgf9F1y77VtUVfdmdgPOAewJfpGvHrYAFwMHAB4Ff0fC+\nmuRpfd2v75/vz+jaeDfgcGBjGDGbvqry5m2jvdH9is8ddKGweOj2kr7M/IFyb59kPfPof9RlaPrL\n+uXeODT9Ff30c4B7Ds3bDLjfwOPD+rKHTrLtxf38/Qem7QLcSvdG9NCh8h/uy390aPqSfvp5wL0H\npt+D7jeBbwfuO812/VK/rr8bmv5Y4Dbgl8A2I16Lr8/w9VsB/HiSeevaBiun2N6uk0z/h359L1jb\na7OW5/Plvvxrhqb/ycA+eOjQvEnbje5Dzh3AvIFpDwX+B7gM2Gmo/JP61/nEgWmvGVWnft7WwFat\n76vA5/rye4yo0/Yz2ec2hptD2dpULKL7KcbB2/CZpT9nkpPEquon1f8vH3IMXS/gqUPTX0PXk3ll\nVd04tK47qurnM34GqzuY7qfjPlRVlw3NewvwG+DgdL+HO+xNVXX9QH1uouvxbgbsvbYNJ3kA8IfA\nf9OdaX2nqjqbrpezPXDgtJ/NulmfNhipqiY7A/x9/f3w6zxtfbs9he53hz80tN2Tgf9a13UPOZxu\ntPO1VXXV0Ha+TjcS8Kwk2wwtd8vwiqrq5qpaY/oMzeW+Ouo5bHI/J2kwa1NQVbXZiNuThsotrarb\nRq0gyZZJXp3kzP4Y2e39McuVdEOAOw+U3QbYHfhFVS2dpec0cRz368Mz+jeyC+h6O783PJvuh+yH\n/bS/324a296rv/9mVa0cMX+iTntOY13rY13bYFJJtumPaZ7XH2Ne2b/Ov+yL7DzV8msx0W5nTvIh\nb0MF82P7+4X9Md/VbnTDwJvTDdkDnEQXjh/uj0m/IsluG6guMDf76qf6+3P74+Iv7D8IbZI8xqy7\nkql6sZ+hO8b8I+Dzfdlb6c5OfR0weGLQxMlgP5uFOk64V39/1STzrxoqd6equmFE+dv7+803wLYn\n2nHSk+I2kHVug1H6E4q+DjwaWEbX87+Gbmg+dKMu63MC2EQ9Jrt2e31HUSbs0N//zRRlCtgGutGg\nJPvQDUM/jX6kI8kVwHuq6oPrWZ9Z31er6vNJngm8nu7w0isBknyX7nDLV9eh3s0ymHVXMvIEmySP\nogvlrwBPr4GTf/qTvN40tMjE0Nv69K7WZuJEmZ2AS0bM32mo3Gxs+36TzJ/NbY+qx4Zqg2fThfIx\nVfXywRn9CVKL1qWSAybqcd9J5k/WnjD5e/GoDz+/ptuX71UDJyROpaouBf40yebAI+iG3F8DvD/J\nb6vq6ClXMLU52Ver6lTg1CRbA/vSXTd/OHBKkr2qatS2N0oOZUvw4P7+5FrzjNzH0J3Beqeq+i1w\nMXC/JNMZzp0YDp5Ob3XC+f39wuEZ6S7f2pPujNzZeDOa2PYT+jfyYQcMlVsfK5m8XdalDVb280dd\nhzvxOp84Yt6G+IrJwXYb9d66cJLlrqM7m381fdvvyZofKM+m6+HP+GtCq2plVZ1fVf8M/Fk/efAL\nYZrfV/vj4mdU1euBdwB3A56+IdbdCoNZ6i7dgVWBA0CS+9CdVTrKB/r7jw5fp5lksySDvaOJL2TY\nZQZ1+hTdEOtrkjxoaN7b6I57f2qyY+bro6p+Rjd6sCurLmcCumtngRfRXV7z+Q2wuWuB+yTZasS8\ndWmDa+lCa1RbT/Y6LwDetQ51X81Auy2g+2rSwW08m8mD9FxglyR/ODT97+muFhj2Ibp2eW9/edZq\nktwtyX4Djx+ZZNRw/8Q+etPAtCb31ST7T/IhceI5/HZd190ih7Kl7pKNs4ADk5zV/31fuuNxlwJX\nMvRNSFX1b/2b3yHAD5OcTHe88v50b/xHseraym/Rvfm9LskOrDoG+YFJjrFRVf+d5HV0HwzOT/Lv\ndCcoPZFuGO8S1hxiZ7ie6+FVdO3w7iRPpbsm+4F01zHfDry0HzlYX1+lu1789CTfpDuuf2FVnbKO\nbfBV4HnAiUlOo+upraiqT9GdrXw58NdJ9qC7Tnce3bWzpwB/ugGez1/S9Wjf17fbRXQ99ef023/W\niGXeA/wRcFKSz9D1oB9Hd5nfEoZ6olX1g/465qOB7yU5ne4Soy3757Mf3T42cYLXocArkpxJd8b4\ndXTX9D9/0t4gAAABjklEQVSL7izn9w2svtV99QPA/fv/n/9Nd7nY3nT/11YAn57Buto37uu1vHlb\nnxvdtY2TXrfal5nflzt6ijLb0b2xLKd7M/8h8Ha6s0mXM/m1ti+ie/O8vl/uR3RfXrLnULk/onvT\nu3GizvTXptId21zJiGtl6S5b+hJdD/UWumtX3wlsO6LsGZO1Bd2lYyuZ5PrUSZa5P/D/6N74bqX7\n0ogTgb2neC1meh3zPfptXEHX61o5/DrNsA02o/sBjR/RvXmvVifgAXQ9vJ/SBdAy4A10Q7dr1H+q\n12aK5/Qgui/ZuI7ubOiz6IZaX8Ik1wjTheR5/T50DXA83QehYwb3laFl/qCfv6Jvl1/SfRA4Elg4\nUG6fvo0vpOsR39S34VHAbiPW29y+SveB8Ph+nTfSHa++iK5HvsP6vIe0eEv/pCVJsyjJYXS93MOq\n6tgxV0cN8xizJEkNMZglSWqIwSxJc6Pwt6I1DR5jliSpIfaYJUlqiMEsSVJDDGZJkhpiMEuS1BCD\nWZKkhhjMkiQ15P8D68rcidN/p9gAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x10a1e7d10>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig3 = plt.figure(3,(6,6))\n",
    "\n",
    "malefirst_maleafter = ['M'*x[1]['questions'].count('M') for x in q.iterrows() if x[1]['speaker'] == 'M']\n",
    "malefirst_femaleafter = ['F'*x[1]['questions'].count('F') for x in q.iterrows() if x[1]['speaker'] == 'M']\n",
    "femalefirst_maleafter = ['M'*x[1]['questions'].count('M') for x in q.iterrows() if x[1]['speaker'] == 'F']\n",
    "femalefirst_femaleafter = [\"F\"*x[1]['questions'].count('F') for x in q.iterrows() if x[1]['speaker'] == 'F']\n",
    "vc_malefirst = pd.value_counts(list(''.join(malefirst_maleafter+malefirst_femaleafter)))\n",
    "vc_femalefirst = pd.value_counts(list(''.join(femalefirst_maleafter+femalefirst_femaleafter)))\n",
    "\n",
    "# Load everything into a single dataframe\n",
    "\n",
    "firstdata = [vc_malefirst,vc_femalefirst]\n",
    "firstlabels = ['Male asks 1st question','Female asks 1st question']\n",
    "\n",
    "firstnormdata = [x/x.sum() for x in firstdata]\n",
    "df = pd.DataFrame(firstnormdata,index=firstlabels)\n",
    "\n",
    "print df\n",
    "\n",
    "# Plot stacked bar chart\n",
    "\n",
    "ax = fig3.add_subplot(111)\n",
    "dfplot = df.plot(kind='barh',stacked=True,ax=ax,color=gencolors,legend=True)\n",
    "\n",
    "ax.set_xlabel('Fraction of total questions',fontsize=20);"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "When women ask the **first** question in a session, women ask on average 40% of the total number of questions in such sessions. However, when men ask the first question in a session, women only ask 31% of the total questions in the session.\n",
    "\n",
    "But this is clearly affected by the fact that the gender of the first question is fixed. To isolate that effect, let's look at the *remaining* questions in the session. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "                                 M         F\n",
      "Male asks 1st question    0.666667  0.333333\n",
      "Female asks 1st question  0.637931  0.362069\n",
      "z = 0.354\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAeYAAAF/CAYAAACc+mpEAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xm4ZFV97vHvC62AKIhggqLYoMYEooiIwQHFMUYxAs4y\niHqdp5vovWqIdnONRqOJUTRqCMgQwBlBQBxpFQQUmZpJRGmjgorK0ArI9Lt/7H3o6qLOOdVNn1OL\n5vvhqadPrb323qtWbeqttfauqlQVkiSpDetMugGSJGkFg1mSpIYYzJIkNcRgliSpIQazJEkNMZgl\nSWrIgkk3QO1L4mfqJGk1VFVWdR1HzBpLVXmrYtGiRRNvQys3+8K+sC9mvq0ug1mSpIYYzJIkNcRg\nllbBLrvsMukmNMO+WMG+WMG+uP1ye+bBdeeQpDxOJGnVJKG8+EuSpDs2g1mSpIYYzJIkNcRgliSp\nIQazJEkNMZglSWqIwSxJUkMMZkmSGmIwS5LUEINZkqSGGMySJDXEYJYkqSEGsyRJDTGYJUlqiMEs\nSVJDDGZJkhpiMEuS1BCDWZKkhhjMkiQ1xGCWJKkhBrMkSQ0xmCVJaojBLElSQwxmSZIaYjBLktQQ\ng1mSpIYYzJIkNcRgliSpIQazJEkNMZglSWqIwSxJUkMMZkmSGmIwS5LUEINZkqSGGMySJDXEYJYk\nqSELJt0A3THsn/0n3QRprbLoiMWTboIa5YhZkqSGGMySJDXEYJYkqSEGsyRJDTGYJUlqiMEsSVJD\nDGZJkhpiMEuS1BCDWZKkhhjMkiQ1xGCWJKkhBrMkSQ0xmCVJaojBLElSQwxmSZIaYjBLktQQg1mS\npIYYzJIkNcRgliSpIQazJEkNMZglSWqIwSxJUkMMZkmSGmIwS5LUEINZkqSGGMySJDXEYJYkqSEG\nsyRJDTGYJUlqiMEsSVJDDGZJkhpiMEuS1BCDWZKkhhjMkiQ1xGCWJKkhBrMkSQ2ZNZiT3JzkrIHb\nlnPVmCTLktxrjra9OMmbV2O9eyU5KcnyJAeMUf8lSe6zeq28fZI8O8lfDNzfP8mTJ9EWSdLqWTBG\nnWuravs5b0mnGtz29cA/An/Z32azL3AecPlq7u/22B34EnAhQFUtmkAbJEm3w2pNZSfZIcmSJGck\nOTHJ5n35kiT/luT7SS5MsmOSo5NcnORdA+sf3a97XpJXTLOPvZKc3o/SP57kNm1N8o4k30uyNMkn\nBsrfmOT8JOckOXJgleqXvyLJCUnWH6p71PA+quraqjoF+OPQvtdNcki/73OT/O8kzwEeCRyR5Mwk\n64/ot3OSnJ3k/UmW9uX7Do7GkxyX5An9309L8t0kP0jymSQb9uXvHWj3+5M8GngW8P5+31v37XtO\nX//Jffm5SQ5Kcte+fFk/m/CDftlDRj0fkqT5MU4wbzAwjf35JAuAA4DnVNUjgU8C7+7rFvDHqtoR\n+BhwDPBqupHmvkk26eu9rF93R+CNA+UA9NOxzwce04/WbwH2HNG2j1TVo6rqoX07d+3L3wo8vKq2\n6/c/sOm8HngG8Oyqun6o7qtm6IfhEffDgftW1UOr6mHAwVX1eeAM4MVV9Yh++4M+Cbyuqh7eb2+6\nUXwBlWQzYD/gyVW1A/AD4O/76f7dqmrbvt3vqqpTgWOBt/T7/snAdtbv9/38vq0LgNcM7OuKfvsf\nA94yQx9IkubYOMF8XVVt39+eA/w5sC3w9SRn0QXHFgP1j+3/PQ84r6p+VVU3AD8B7t8ve1OSs4FT\n+7IHD6wf4MnADsAZ/T6eBGw1om1PSnJaknP7Otv05ecCRybZE7h5YLv7AE8HnltVN85Qdxw/BrZO\n8uEkfw0sH3oMK0lyT2Djqjq5Lzp8VL2hbezUP6bv9v2wD7AlcDVwfT/y3R24boZ9B3gIcGlVXdKX\nHQo8fqDOF/p/zwQWztAmSdIcG+cc87AA51fVY6ZZPjXlewsrT//eAixIsgtd8O5UVdcnOQlYn9s6\ntKr+YdpGdKPAjwI7VNUvkiwCNugXP5MueJ4F7JfkoXQjw6XAdnRvBpZNV7eqZg3oqroqycPogv7V\ndCP8l08tnm19Vg7Qm1j5TdJgf3ytql58m5WTR9H143OB1/d/T7fv4bIMlU09TzczzTFxEifd+vdC\nFrLVyPdJknTnteQCWHLh7d/O6gTzD4F7J9mpqk5LchfgwVV1wRjrBtgIuLIP5T+nGxUOKuAbwDFJ\nPlhVV/RTt3evqv8ZqDcVXr9NcnfgecBnkgTYsqqWJDkFeCFw977uWXTTtcf2o9xfjqi7IXDNNG1f\ncSfZFLixqr6Q5GLgsH7R8v4xrvyguiC/Kslj+3PWg1Pzy4DX9G2/H/Covh9OAz6a5IFV9eP+/PJ9\ngcuADavqy0m+Szd6n27fRfecLZzaDrA38K0Rj3FaT+SJq1Jdku50dtmmu03Z/wvT153JOMG80mir\nqm5I8lzgw0k27rfxQWA4mEedQy3gRODVSS6gC4xTb7PDqguT/CPw1f6irxuB1wL/M1DnqiQH0k2Z\n/xI4vV+0LnB437YAH6qqq7vMo6rqlCRvAY4Hnjai7m1COcky4B7AXZM8u1/vrsAnBy5Ke1v/7yHA\nx5NcS3eOfPA880uBg5MU8NWBx3Jykkv7PryQ7lwyVfWbJPsCRyVZr6++H10AH9PPGgT4u37Zp4AD\nk7yB7o3K1Pb/mOSlwGf7awS+B3x8avFg1zPeaF+SNEdS5evwJCR5AHBcf+Fa05LUYhZPuhnSWmXR\nEYsn3QTNsewJVTXTtUQj+c1fkzN8nleSpNU6x6w1oKqWAQ+bdDskSW1xxCxJUkMMZkmSGmIwS5LU\nEINZkqSGGMySJDXEYJYkqSEGsyRJDTGYJUlqiMEsSVJDDGZJkhpiMEuS1BCDWZKkhhjMkiQ1xGCW\nJKkhBrMkSQ0xmCVJaojBLElSQwxmSZIaYjBLktQQg1mSpIYYzJIkNcRgliSpIQazJEkNMZglSWqI\nwSxJUkMMZkmSGmIwS5LUEINZkqSGGMySJDXEYJYkqSEGsyRJDTGYJUlqiMEsSVJDDGZJkhpiMEuS\n1BCDWZKkhqSqJt0GNS5JeZxI0qpJQlVlVddzxCxJUkMMZkmSGmIwS5LUEINZkqSGGMySJDXEYJYk\nqSEGsyRJDTGYJUlqiMEsSVJDDGZJkhpiMEuS1BCDWZKkhhjMkiQ1xGCWJKkhBrMkSQ0xmCVJaojB\nLElSQwxmSZIaYjBLktQQg1mSpIYsmHQDJEl3Pkkm3YQ1qqrW2LYMZknSRKzJMJukNf0mw6lsSZIa\nYjBLktQQg1mSpIYYzJIkNcSLvyRJE7d/9p/T7S+qRWPVW7hwIZdffjmXXXYZm2666a3l22+/Peec\ncw7Lli1jyy23nKtmAo6YJUm6VRK23nprjjrqqFvLli5dynXXXTdvH/EymCVJGrDXXntx2GGH3Xr/\n0EMPZZ999pm3j3cZzJIkDdhpp5245ppruOiii7j55pv59Kc/zV577TVv+/ccsyRJQ/bee28OO+ww\nHv/4x7PNNtuwxRZbzNu+DWZJkgYkYe+992bnnXfm0ksvnddpbHAqW5Kk29hyyy3Zeuut+fKXv8we\ne+wxr/t2xCxJ0ggHHXQQV111FRtssAE33XTTvO3XYJYkTdy4nzOeT1tvvfVK9+fr41JZW37dQ3Mn\nSXmcSFqTkqxVvy416rH05auc5p5jliSpIQazJEkNMZglSWqIwSxJUkMMZkmSGmIwS5LUEINZkqSG\nGMySJDXEYJYkqSF+JackafKOnOOvu3zxeN8ytnDhQn7961+z7rrrAt23d1188cVsvvnmc9m6lThi\nliSpl4TjjjuO5cuXs3z5cq655pp5DWUwmCVJaorBLEnSgEn/uIbBLElSr6rYbbfd2GSTTdhkk03Y\nY4895r0NXvwlSVIvCccccwxPetKTJtYGR8ySJDXEYJYkqSFOZUuSJm/MzxnfGWTSV5+pfUnK40TS\nmpRk4lc/rynTPZa+fJW/OcWpbEmSGmIwS5LUEINZkqSGGMySJDXEYJYkqSEGsyRJDfFzzJKkiUjm\n+DeY76AMZknSvFtbPsM8F5zKliSpIQazJEkNMZglSWqIwSxJUkMMZkmSGmIwS5LUEINZkqSGGMyS\nJDXEYJYkqSEGsyRJDTGYJUlqiMEsSVJDDGZJkhpiMEuS1BCDWZKkhhjMkiQ1xGCWJKkhBrMkSQ0x\nmCVJasiCSTdAdwz7Z/9JN0G601l0xOJJN0ET4IhZkqSGGMySJDXEYJYkqSEGsyRJDTGYJUlqiMEs\nSVJDDGZJkhpiMEuS1BCDWZKkhhjMkiQ1xGCWJKkhBrMkSQ0xmCVJaojBLElSQwxmSZIaYjBLktQQ\ng1mSpIYYzJIkNcRgliSpIQazJEkNMZglSWqIwSxJUkMMZkmSGmIwS5LUEINZkqSGGMySJDXEYJYk\nqSEGsyRJDTGYJUlqiMEsSVJDDGZJkhpiMEuS1BCDWZKkhhjMkiQ1xGCWJKkhBrMkSQ2Zk2BOckuS\nwwfuL0hyRZIvzbLeLrPVuZ3tWpbkXqux3vOSnJ/k5iSPmKXuxkles/qtvH2S/MPQ/VMm1RZJ0qqb\nqxHzH4Btk6zf338q8HOg5mh/41rd/S8Fdge+PUbdTYDXruZ+1oS3D96pqsdOqiGSpFU3l1PZJwDP\n7P9+EXAUEIAkj0ry3SRnJjklyZ8Nr5xkwyQHJzm9r/e309T5epIfJDl3qk5ffnySs5MsTfK8ofU2\nSPLlJC9Pcrehus8f3k9VXVRVF4/Y/7Z9+87q138Q8F7ggX3Z+0ass1+SHyb5TpIjk7y5L1+SZIf+\n782SXNr/vW6S9yf5XpJzkryyL79Pkm/3+1ma5HFJ3gts0Jcd3tf7ff9v+u0s7fvq+X35Lv2+P5vk\nwiT/PfLZlCTNiwVzuO1PA+9MchzwUOAgYOd+2YXAzlV1c5KnAO8Bnju0/n7AN6rqZUnuCZye5OtV\nde1AneuB3atqeZLNgFOBY4GnA7+oqmcCJLnHwDr36Nt2aFX9d5LnDNXdaBUe46uAD1XVkUkW0PXn\nW4Ftq2r74cp98L4A2A64C3AmcEa/uBg9on85cFVVPSrJesDJSb4K7AGcWFXvSbIOcLeqOjnJ64b2\nPbXNPfr9Pgy4N/D9JFMzAA8HtgEuB05J8tiqcgpckiZgzoK5qpYmWUg3Wj5+aPE9gcP6EWbRhdSw\npwHPSvKW/v56wP2BHw7UWQf45yQ7A7cA903yJ8C5wAf6EeRxVXVyXz/AMcD7quqovmy6uuM4Fdgv\nyf2AL1TVJUkyQ/2d+3rXA9cnOXaMfTwNeGiSqTcuGwEPAr4PHJzkLsAXq+qcWbbzOODIqirg10m+\nBewIXAN8r6ouA0hyNrAQMJglaQLmcsQM3ej1A8AT6EZpU95FNxrePckDgCXTrL9HVf1ohu3vCWwG\nPKIffV8KrF9VP0qyPd1U+j8l+UZVvYvuTcDJwN/QTa0zQ91ZVdVRSU4DdgVOSPIq4NKZVqGfzu8N\n/n0TK04trM/KXl9VXxveWP+GZFfgkCT/VlWHD9eZYd9TZQB/HCi7mRHHxUmcdOvfC1nIVmw1w64k\n6c5nyQWw5MLbv525/rjUwcDiqjp/qHwj4LL+75dOs+5XgDdO3enDc9hGwK/7UH4i8IC+7n2A66vq\nCLo3BoPrvhO4MslHp6k741XXDIRbkq2q6tKqOoBuJP5QuhHoPaZZ99vAbknW76fXdx1Ytgx4ZP/3\n4LT+V4DX9lPlJPmz/rz4lsAVVfVfdKcJph7jjVN1h3wHeEGSdZLcG3g88D1uG9YjPXHgP0NZkm5r\nl21g8XNW3FbXXAVzAVTVL6rqIwNlUyO0f6Gbgj4TWJeVz61O/f0u4C79hUrnAfuP2M8RwCOTnAvs\nTXfuGrqAPD3JWXRB/E8rNa7qTXQXSb1vqO47+v2uJMnuSX4G7AQcn+TL/aLnJzmvX3db4LCq+h3d\nedqlwxd/VdVZdOe3z6G7OO77rAjGDwCv6ftk04F++C/gAuDMJEuBj9GNaHcBzu7rPw/4UF//P4Fz\ns+LjalPPxdF00/bnAN8A/k9V/ZrR57YnffW8JN1ppTvlqElIsgj4fVX966TbMpMktZjFk26GdKez\n6IjFk26CbofsCVU11qzkIL/5a/J8ZyRJutVcX/ylGVTVqOl5SdKdmCNmSZIaYjBLktQQg1mSpIYY\nzJIkNcRgliSpIQazJEkNMZglSWqIwSxJUkMMZkmSGmIwS5LUEINZkqSGGMySJDXEYJYkqSEGsyRJ\nDTGYJUlqiMEsSVJDDGZJkhpiMEuS1BCDWZKkhhjMkiQ1xGCWJKkhBrMkSQ0xmCVJaojBLElSQwxm\nSZIaYjBLktQQg1mSpIYYzJIkNcRgliSpIQazJEkNMZglSWqIwSxJUkMMZkmSGmIwS5LUEINZkqSG\nGMySJDUkVTXpNqhxScrjRJJWTRKqKqu6niNmSZIaYjBLktQQg1mSpIYYzJIkNcRgliSpIQazJEkN\nMZglSWqIwSxJUkMMZkmSGmIwS5LUEINZkqSGGMySJDXEYJYkqSEGsyRJDTGYJUlqiMEsSVJDDGZJ\nkhpiMEuS1BCDWZKkhhjMkiQ1xGCWJKkhBrMkSQ0xmCVJaojBLElSQwxmSZIaYjBLktQQg1mSpIYY\nzJIkNcRgliSpIQazJEkNMZglSWqIwSxJUkMMZkmSGmIwS5LUEINZkqSGGMySJDXEYJYkqSEGs7QK\nlixZMukmNMO+WMG+WMG+uP0MZmkV+KKzgn2xgn2xgn1x+xnMkiQ1xGCWJKkhqapJt0GNS+JBIkmr\noaqyqusYzJIkNcSpbEmSGmIwS5LUEINZt0ry9CQXJflRkrdOU+fD/fJzkmw/322cL7P1RZI9+z44\nN8kpSR42iXbOh3GOi77ejkluSrLHfLZvvoz5/8cuSc5Kcl6SJfPcxHkzxv8fmyU5McnZfV/sO4Fm\nzoskByf5VZKlM9RZtdfNqvLmDWBd4BJgIXAX4GzgL4bqPAM4of/7r4DTJt3uCfbFo4GN+7+ffmfu\ni4F63wSOA54z6XZP6Ji4J3A+cL/+/maTbvcE+2Ix8M9T/QD8Flgw6bbPUX/sDGwPLJ1m+Sq/bjpi\n1pRHAZdU1bKquhH4FPDsoTp/CxwKUFWnA/dM8qfz28x5MWtfVNWpVXV1f/d04H7z3Mb5Ms5xAfAG\n4HPAFfPZuHk0Tj+8GPh8Vf0coKp+M89tnC/j9MXlwEb93xsBv62qm+axjfOmqr4DXDlDlVV+3TSY\nNWUL4GcD93/el81WZ20MpHH6YtDLgRPmtEWTM2tfJNmC7oX5Y33R2vhRj3GOiQcD90pyUpIzkuw9\nb62bX+P0xYHAtkkuA84B3jRPbWvRKr9uLpjT5uiOZNwX0+HP5K2NL8JjP6YkTwReBjx27pozUeP0\nxb8Db6uqShJue4ysDcbph7sAjwCeDNwNODXJaVX1ozlt2fwbpy/+ATi7qnZJ8kDga0m2q6rlc9y2\nVq3S66bBrCm/AO4/cP/+dO/sZqpzv75sbTNOX9Bf8HUg8PSqmmkq645snL7YAfhUl8lsBvxNkhur\n6tj5aeK8GKcffgb8pqquA65L8m1gO2BtC+Zx+uIxwLsBqurHSS4FHgKcMS8tbMsqv246la0pZwAP\nTrIwyV2BFwDDL6zHAvsAJNkJuKqqfjW/zZwXs/ZFki2BLwB7VdUlE2jjfJm1L6pq66raqqq2ojvP\n/Jq1LJRhvP8/jgEel2TdJHeju9Dngnlu53wYpy8uAp4C0J9PfQjwk3ltZTtW+XXTEbMAqKqbkrwe\n+ArdVZcHVdWFSV7VL/9EVZ2Q5BlJLgH+ALx0gk2eM+P0BfBOYBPgY/1I8caqetSk2jxXxuyLtd6Y\n/39clORE4FzgFuDAqlrrgnnMY+I9wCeTnEM3APy/VfW7iTV6DiU5CngCsFmSnwGL6E5rrPbrpl/J\nKUlSQ5zKliSpIQazJEkNMZglSWqIwSxJUkMMZkmSGmIwS5LUEINZa40k+ya5ZZrbk+apDRsnWTzq\np92SLEly0ny0Y3UkeVaSpUmu6/tso9nXaleSZUkOXo31FvaPf5+5aNcdUf9zlov6r1wdLLev5oBf\nMKK10XO57VcEXjhP+96E7stH/gc4a2jZq+epDassyQLgCOBk4DXADcDvJ9qo2+/ZwDWrsd5lwE7A\nj9dsc+7QdqE7rt/Fyt/zbF/NAYNZa6Ozq2qsr/9Lsl5V/XEO2nCbH3KoqovmYD9ryhbA3YHPVtXJ\nq7pykrtW1Q1rvlmrr6rOWc31bgC+t4abs7ZY6bi2r+aGU9m60xiY6t45yWeTXAmc1i/bMcnnkvws\nybVJLkry7iTrj9jO7klOSbI8ydVJTu+ngRey4vuADxyYRp/6ntzbTGUneUiSo5Nc2e/31CR/PVRn\ncb+dByU5vt/vsiTvGJ5anOZx3yfJYUmuSHJ9knOS7Dm4feDS/u5B/b6mnXJPckjfT49O8t0k1wL/\n0i+7d5KPJ/l5v68Lk7ximufhMX2fX5Pkl0ne1i/ftW/jH5J8L8kjhtZ/WpITklzW11ma5O+TrDNU\nb1mST47Y718lOaJ/7n6R5ENJ1huoNzU9+5IRj/nhSb7T7/fi9F9DObTfpyQ5qz8l8KMkL+/Xv3S4\n7oh1753kyL5tVyY5NMlufXseP91jGyi/JcmiobLtkhyb5Hf9MXZykscN1dkxydeS/Kav8+MkH+2X\nLaYbLQPcOHVcT9dXffle/XN4XX/cHZZk86E6y5IcnuSF/XHy+yTfT/LYoXrTtm1t5YhZa6MF6aZm\np1RV3Txw/wjgSLrfD56qtyXd78YeClwF/CXdi9HWwIumVkzyBuBDwNHAB+ime3cAHkD33cF70P24\nxXtY8cX+g9N8NbCt+9JNHV8NvI5u2vV1wPFJdq2qE4ce19HAwcC/0v34+v50v2h0yHQdkWRD4FvA\nxsDb+/p7A4cnuVtVHUj3C1lLgc/STVUez+xTwBsDRwHvB95G92tKG/WPZz267wu+FHg63feJr1dV\nHxnaxiF0/f0fwPOB96T7wYOn9O34A13gfzHJA6vqxn69rYBvAh/p6+wILAbu3T/GKcXon9c7nO75\n353uV5AW0/3Q/eKhesPrbtSv98G+7sv6x/bDqloCkGQbuv47je7HHdYD3tH3183M7gvAQ/vH8SPg\nhcABI9oy3WNbqd39m5rvAD8A/hdwHd0pla8neUxVnZnk7nTH7mnAS4DldH386H4zB9LNqLyc7udN\nRz2OwX2+Evg48Cngrf267wH+KskjquoPA+vsDPwZsB/wR7rn/bgkC6vq6jHatnaqKm/e1oobsC/d\njwcM3749tPxfZ9lO6AJ7L7oXoU368o3oXhg+N8O6C/t9vGzEsiXANwfufwC4Edh6oGwdul/m+cFA\n2eJ+my8Z2t65wFdmeSyv79d9/FD514BfAev09x/U19tnjH4+pK/7rKHyd9C98D9wqPw/gSsG9jX1\nPPzjQJ11gV/Tndt+wED5s/q6O8/yXO0H/G5o2aXAwSOOj0VD9b4E/HDEc7jPiMf8hIGyuwK/AT4x\nUHZk36/rD5RtDlwP/GSWfn1qv4/nD5WfMPwcDj+2gfJbgHcO3P8GcD6wYOgYuwA4ur//yH69v5yh\nbVPH4DrTHO/7DDyPvwK+MVTvsX29NwyULQN+C2w8ULZDX+9F47Ztbbw5la210W50/0NP3V4+tPzo\n4RWSbJTkfUl+TPciegNwGN0L/4P7ao8BNqQLmjXh8cCpNXA+vKpuoRtpPLwfLQw6fuj++XQj/dn2\n8fOq+vZQ+RF0I8y/WOVWd24AjhsqezrdyGZZkgVTN+CrwKbANkP1vzz1R3UzGpfQBeRPB+r8sP/3\n1t+zTTc1/4kkP6UbZd1AN9LaOMmfjNH24X48j9n7EeAPVfWtgTbfAFzMyr+1uxNwQlVdP1Dvl8Ap\nY2z/0XRvBD8/VP6pMda9jSQb0D3/n+3vTz0f69AF9tTU+I/oZon+M8meSe4/antjegjdcXXEYGFV\nnQL8lO5XmAadWlVXD9w/r/93qg1rsm13GAaz1kbnVdWZA7fhH6q/fMQ6nwReBfw73VTqI+mmlQGm\nzjNv2v87fMX36rrXNG35Jd0bgk2Gyod/Nu+PA21bnX1MLV8dV1Q/pBnwJ3QvvDfSheXU7TN005ab\nDtW/cuj+DdOUQf84+/PIxwLPAP4f8ES65+rddH02W3/A6H5cb1TFWdo71b7BfW5ON/IfNqps2H2A\nK2vl0y7jrjvKvehGsO9k5efjBrpj+54AfTA+ke4K6/8Aftqft99jNfcJo4+5X7HyMV0MPRe14kLM\n9eegbXcYnmPWndFKgZLuAq+/pZviPGCgfLuh9X7T/3s/uqnA2+u3dC/Gwzbv2zgqCFbV7+jO4Y3a\nx9TyNeU3dIH/pmmWX7wG9vFAuunOvarqyKnCJM9eA9uezawX2tEF0p+OKB9VNmrdTZKsOxTOo9a9\nnm4qfUXjkuE3PlfRTQN/hG72Z1rVXcH+3P6Nz45057g/k2S7qjp/jLZPmTqepjuuv78K21rTbbvD\ncMQsdaOldYGbhsr3Hbr/XbqLvV45w7am3vFvMMZ+vwXslOQBUwVJ1qW7aOjMqhrnc8Sz/aD6EuB+\nSR4zVP5iuhHM6r7BGLXfE+mmxn82NGMxdVsTn4u+W//vrc9VkrsAe07TpjVpnO2fBjyjn0YGuql3\nunOss/ku3XH43KHyF46o+1O6i8QGPXPwTnUXWX0HeDhw1qjnZHijVXVLVZ1ON8peB/jzftHUcX23\n4XWG/JDuuFqpzf3xtyXd8bhaZmjbWscRs+70qrv68zTgzUkupxvJvgy471C95UneDhyQ5HN0F/os\np3vhu666q45/1a//oiRLgWvpLvqZGkkMjro+SBf+X+s/4rIceC3dhVgrvcjOYLZR3CF0I9gvJNkP\n+AVdiD0FeOWI6ehxjdrvB+neVHwnyQfpRsgb0r2APq6qdlvN7Q66gC6U3p3kZrqA/ju60Bxed5wR\n7qqYbnuD5f9EF6xfSfIBuinZd9DNJNwy08ar6utJTgY+kWQzunPuLwC2HVH9U8DBSf6N7pz5dnRX\nLQ/7e+Dz18p7AAAB+ElEQVTbfXsO6tuxGfAIugu53p5kV7o3m0fTXZC1IfBGuivzT+23MzUyfXOS\nE4Gbq+qMEY/h5iTv7B/D4XTnmregO9VwMd2nCqaM81G/cdq21nHErLXNbEEz3fIX0X2k5KN055sv\nowu0lepX1UeB59FNZ/838Dm6j0j9pF9+C93HUjYBvk735Qu7Duy7BrZ1OfA4uhe9j9FdpHNP4JlV\n9dWhNo9q90wfmZnax7V0532/CrwX+CLdSGuvqvqvmdadabOj9ltV19BdIHcC3cdkTgQOoruy+psj\ntjHWdof2cSPdxX2/pJuePYBuFPbeEetOt49V3u8MdYaf0wvp3lTdg+7c+nuAD9MdW1ePWH/YHnT9\n98904bsO3ZX1ww6l+0jaHnTn3J9K9/GvlRtXdRbd9O9v+3Z8he46im3pZmygC8xr6d5AnEAXnjcA\nT62qy/o6x9Gd430t3cj+9OkeQHUfwdub7jj7IvC+fr9PqKrrBqtO2wsrjNO2tU5W/w2zJGk2/dX1\nlwBfqqpXzFZ/xPq70L2x2WXE1fVaCzmVLUlrUJID6EaVl9GdDnkT3ReMfGiS7dIdh8EsSWvWenRT\n639KN+16OvCUqjpvxrVm5tTmnYhT2ZIkNcSLvyRJaojBLElSQwxmSZIaYjBLktQQg1mSpIYYzJIk\nNeT/A/YrsSj5FPT4AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x1076cd1d0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "malefirst_maleafter = ['M'*x[1]['questions'][1:].count('M') for x in q.iterrows() if x[1]['questions'][0] == 'M']\n",
    "malefirst_femaleafter = ['F'*x[1]['questions'][1:].count('F') for x in q.iterrows() if x[1]['questions'][0] == 'M']\n",
    "vc_malefirst_remaining = pd.value_counts(list(''.join(malefirst_maleafter+malefirst_femaleafter)))\n",
    "\n",
    "femalefirst_maleafter = ['M'*x[1]['questions'][1:].count('M') for x in q.iterrows() if x[1]['questions'][0] == 'F']\n",
    "femalefirst_femaleafter = [\"F\"*x[1]['questions'][1:].count('F') for x in q.iterrows() if x[1]['questions'][0] == 'F']\n",
    "vc_femalefirst_remaining = pd.value_counts(list(''.join(femalefirst_maleafter+femalefirst_femaleafter)))\n",
    "\n",
    "# Load everything into a single dataframe\n",
    "\n",
    "firstrdata = [vc_malefirst_remaining,vc_femalefirst_remaining]\n",
    "firstrlabels = ['Male asks 1st question','Female asks 1st question']\n",
    "\n",
    "firstrnormdata = [x/x.sum() for x in firstrdata]\n",
    "#print firstrnormdata\n",
    "dfr = pd.DataFrame(firstrnormdata,index=firstrlabels)\n",
    "\n",
    "print dfr\n",
    "\n",
    "# Plot stacked bar chart\n",
    "\n",
    "fig = plt.figure(4,(6,6))\n",
    "ax8 = fig.add_subplot(111)\n",
    "dfplot = dfr.plot(kind='barh',stacked=True,ax=ax8,color=gencolors,legend=True)\n",
    "\n",
    "ax8.set_xlabel('Fraction of remaining questions',fontsize=16)\n",
    "ax8.set_xlim(0,1)\n",
    "\n",
    "# Check statistical significance\n",
    "# http://stats.stackexchange.com/questions/113602/test-if-two-binomial-distributions-are-statistically-different-from-each-other\n",
    "\n",
    "def ztest(p1,p2,n1,n2):\n",
    "    pexp = (n1*p1exp + n2*p2exp) / (n1+n2)\n",
    "    z = (p1exp-p2exp) / np.sqrt(pexp*(1. - pexp)*(1/n1 + 1/n2))\n",
    "    return z\n",
    "\n",
    "p1exp,p2exp = dfr['M']\n",
    "n1 = len(list(''.join(malefirst_maleafter+malefirst_femaleafter)))\n",
    "n2 = len(list(''.join(femalefirst_maleafter+femalefirst_femaleafter)))\n",
    "z = ztest(p1exp,p2exp,n1,n2)\n",
    "print 'z = %.3f' % z"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "When the first question is eliminated from the data, the fraction of questions asked by females are nearly identical for females and males (33% and 36%, respectively). A two-proportion $z$-test indicates the distributions cannot reject the null hypothesis at $<1\\sigma$ level. \n",
    "\n",
    "This result differs from the AAS meeting data, who found an extremely strong difference between the two. If the first question was asked by a male, then only 10% of the remaining questions were asked by a female. If the first question was asked by a female, 49% of the remaining questions were asked by a female. "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Does the gender of the session chair affect the distribution of the questioners' gender?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "                     M         F\n",
      "Male chair    0.640000  0.360000\n",
      "Female chair  0.648649  0.351351\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAagAAAGDCAYAAAB+2YQ1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAG/1JREFUeJzt3XvcnGV95/HvlyQIUpCAlXNIomhFWQhnrUBAsFCLUtTS\nasBI19OuqFu7sqt9maRWbV1eqwWUui6kUBHZVigUARVMtoLAghAOQjhoIgjIyUACBCTJb/+4rieZ\nzDPzPDPPYeaXJ5/36zWvmbnv677va6655/7OfXZECACAbLbodwUAAGiFgAIApERAAQBSIqAAACkR\nUACAlAgoAEBKk/tdgU2JbY7JB4ARiAh3OwxrUF2KCB4RmjdvXt/rkOVBW9AWtMXQj5EioAAAKRFQ\nAICUCCiMyOzZs/tdhTRoiw1oiw1oi9HzaLYPbm5sB+0FAN2xreAgCQDAREFAAQBSIqAAACkRUACA\nlAgoAEBKBBQAICUCCgCQEgEFAEiJgAIApERAAQBSIqAAACkRUACAlAgoAEBKBBQAICUCCgCQEgEF\nAEiJgAIApERAAQBSIqAAACkRUACAlAgoAEBKBBQAICUCCgCQEgEFAEiJgAIApERAAQBSIqAAACkR\nUACAlAgoAEBKBBQAICUCCgCQEgEFAEiJgAIApERAAQBSIqAAACkRUACAlCb3uwKbmgVe0O8qABPa\nvAvn97sKSII1KABASgQUACAlAgoAkBIBBQBIiYACAKREQAEAUiKgAAApEVAAgJQIKABASgQUACAl\nAgoAkBIBBQBIiYACAKREQAEAUiKgAAApEVAAgJQIKABASgQUACAlAgoAkBIBBQBIiYACAKREQAEA\nUiKgAAApEVAAgJQIKABASgQUACAlAgoAkBIBBQBIiYACAKREQAEAUiKgAAApEVAAgJQIKABASgQU\nACAlAgoAkBIBBQBIadiAsr3W9m0Nj2njVRnby23vMMpxzLV9VpfDfNP260czXQDA2JrcQZnnI2LW\nuNekiH6MIyI+2Kq77S0iYt3oqwQA6NaINvHZPsD2Ytu32L7a9s61+2Lb/9P2zbbvsX2Q7Utt32f7\n8w3DX1qHvct2u3CYY/umutb2D7YH1bWO/3rbS2zfaPt3aq9dbV9Vp/t3DeW/Xut2l+35Dd0X296/\nvn7W9hm2l0g6dCTtAwAYvU4CauuGzXvftT1Z0lmS3hURB0paKOkLtWxIejEiDpJ0jqTLJH1E0hsl\nzbU9tZY7tQ57kKSPN3SXJNXNbX8i6c117W2dpPc1ldlS0nckfTwi9pN0tKTVkixpvzr8PpJOsr1b\nHeyztW77SjrC9j4N9R7wckk3RsR+EfGTDtoHADAOOtnEt7pxE5/tN0p6g6RrbEvSJEmPNJS/vD7f\nJemuiHisDvcLSXtIWiHpE7ZPqOX2kLSXpP83MAlJb5V0gKRb6jS2lvTrpnq9TtKjEfFTSYqIZ+t0\nQtK1EbGqvr9b0p6SHlYJqw/Wz72LpNdLurNpvGslfbeDdgEAjKNOAqqZJf0sIt7cpv+L9Xldw+uB\n95Ntz1YJoEMj4gXbiyRt1WI850fEZ0ZQv8Y6SCVwJtueIelTkg6MiGdsL2wz3Rciou1+rEVatP71\ndE3XDM0YYRUBYGJafLe0+J7Rj2ckAXWvpN+1fWhE3Gh7iqS9IuLuDoa1pO0krajh9HsavJ8nJF0r\n6TLbX4mIJ+qRfb8TEQ821WMX2wdGxC22t5X0fJ1Gq+luK+k5SStt7yTpOKkhbTp0pI7sdhAA2KzM\n3rs8Biy4ZGTj6SSgNlqbiIjf2n63pDNtv6KO4yuSmgMqmoet76+W9JG66e1eSTcMmmDEPbb/StIP\n6sERL0n6T5IebCjzW9snSTrL9tYq4XRMu+lGxB22b5O0VNJDkq7r5PMCAPrDQ2zNQhPbMV/z+10N\nYEKbd+H8flcBY8zvkyKi1datIXElCQBASgQUACAlAgoAkBIBBQBIiYACAKREQAEAUiKgAAApEVAA\ngJQIKABASgQUACAlAgoAkBIBBQBIiYACAKREQAEAUiKgAAApEVAAgJQIKABASgQUACAlAgoAkBIB\nBQBIiYACAKREQAEAUiKgAAApEVAAgJQIKABASgQUACAlAgoAkBIBBQBIiYACAKREQAEAUiKgAAAp\nEVAAgJQIKABASgQUACAlAgoAkBIBBQBIiYACAKREQAEAUnJE9LsOmwzbQXsBQHdsKyLc7XCsQQEA\nUiKgAAApEVAAgJQIKABASgQUACAlAgoAkBIBBQBIiYACAKREQAEAUiKgAAApEVAAgJQIKABASgQU\nACAlAgoAkBIBBQBIiYACAKREQAEAUiKgAAApEVAAgJQIKABASpP7XQEA2JzZ7ncVxlREjNm4CCgA\n6LOxXKj301iHLZv4AAApEVAAgJQIKABASgQUACAlDpIAgGQWeMG4jn9ezBu2zPTp0/Xoo4/qkUce\n0Y477ri++6xZs3T77bdr+fLlmjZt2nhWkzUoAMBgtjVz5kxddNFF67vdeeedWr16dc8OjSegAAAt\nzZkzRxdccMH69+eff75OOeWUnh0WT0ABAFo69NBDtXLlSi1dulRr167VxRdfrDlz5vRs+uyDAgC0\ndfLJJ+uCCy7Q4Ycfrr333lu77bZbz6ZNQAEAWrKtk08+WYcddpiWLVvW0817Epv4AABDmDZtmmbO\nnKmrrrpKJ554Yk+nzRoUAGBI5557rp5++mltvfXWWrNmTc+mS0ABQDKdnKfUSzNnztzofa8OM/dE\nuYpuL9gO2gvAWLI9oa5m3uqz1O5dpxr7oAAAKRFQAICUCCgAQEoEFAAgJQIKAJASAQUASImAAgCk\nREABAFIioAAAKXGpIwDI5tvjfCmh9w5/5Yrp06fr8ccf16RJkySVq0Hcd9992nnnnce3bg1YgwIA\nDGJbV1xxhVatWqVVq1Zp5cqVPQ0niYACACRFQAEAWur3RWwJKADAIBGhE044QVOnTtXUqVN7frNC\niYMkAAAt2NZll12mo446qm91YA0KAJASAQUASIlNfACQTQfnKW0OuOV7F7jlO4Cxxi3f22MTHwAg\nJQIKAJASAQUASImAAgCkREABAFIioAAAKXEeFAD0mT3O93/aRBFQANBHE+UcqPHAJj4AQEoEFAAg\nJQIKAJASAQUASImAAgCkREABAFIioAAAKRFQAICUCCgAQEoEFAAgJQIKAJASAQUASImAAgCkREAB\nAFIioAAAKRFQAICUCCgAQEoEFAAgJQIKAJDS5H5XYFOzwAv6XQVgwpl34fx+VwEJsQYFAEiJgAIA\npERAAQBSIqAAACkRUACAlAgoAEBKBBQAICUCCgCQEgEFAEiJgAIApERAAQBSIqAAACkRUACAlAgo\nAEBKBBQAICUCCgCQEgEFAEiJgAIApERAAQBSIqAAACkRUACAlAgoAEBKBBQAICUCCgCQEgEFAEiJ\ngAIApERAAQBSIqAAACkRUACAlAgoAEBKBBQAICUCCgCQEgEFAEiJgAIApERAAQBSIqAAACmNS0DZ\nXmf7nxreT7b9hO1/G2a42cOV6XD6i20f0EX5XW3/82inCwAYO+O1BvWcpDfY3qq+P0bSryTFOE2v\nWVfTiYhHIuI9zd1tTx67KgEAujGem/iulPT2+vrPJF0kyZJk+2DbP7F9q+3rbb+2eWDb29g+z/ZN\ntdw7Wk3E9um277C9xPYXG3q9pw57r+231LLTbf+77Z/Wx5saut9ZX8+1fbntayX9cKwaAwDQnfFc\nQ7hY0udsXyFpH0nnSjqs9rtH0mERsdb20ZK+KOndTcN/VtK1EXGq7e0l3WT7moh4fqCA7eMkvUPS\nwRHxQi03YFJEHFLLzFNZi3tM0jER8aLtvSR9W9JBLeo+S9I+EfH06JoAADBS4xZQEXGn7ekqa0/f\na+q9vaQLbL9GZXPclBajeJuk423/ZX3/Mkl7SLq3ocxbJZ0XES/UaTYGyiX1+VZJ0+vrLSWdbXtf\nSWslDVpzq35AOAFAf433PpbLJZ0h6QhJv9vQ/fMqa0d/bHtPSYvbDH9iRNw/zDTcpvuL9XmtNnzO\n/yLp0Yg42fYkSS+0Gfb5Nt21SIvWv56u6ZqhGcNUDwA2L4vvlhbfM/rxjHdAnSdpRUT8zPbshu7b\nSXqkvv5Am2G/L+njkk6TJNuzIuK2pjI/VNmMeGFErLY9NSJWDFGf7VQO1pCkUyRN6vyjFEfqyG4H\nAYDNyuy9y2PAgkvalx3KeB0kEZIUEQ9HxNkN3QaOrvuypC/ZvlUlJKJ5WJW1rCn1AIi7JC0YNJGI\n76uspd1i+zZJnxqqPpK+Lun9tpdIep2kZ1uUaawnAKBPHMGyuFO2Y77m97sawIQz78L5/a4CxpHf\nJ0VEu90xbXElCQBASgQUACAlAgoAkBIBBQBIiYACAKREQAEAUiKgAAApEVAAgJQIKABASgQUACAl\nAgoAkBIBBQBIiYACAKREQAEAUiKgAAApEVAAgJQIKABASgQUACAlAgoAkBIBBQBIiYACAKREQAEA\nUiKgAAApEVAAgJQIKABASgQUACAlAgoAkBIBBQBIiYACAKREQAEAUiKgAAApEVAAgJQIKABASgQU\nACAlAgoAkBIBBQBIiYACAKREQAEAUnJE9LsOmwzbQXsBQHdsKyLc7XCsQQEAUiKgAAApEVAAgJQI\nKABASgQUACAlAgoAkBIBBQBIiYACAKREQAEAUiKgAAApEVAAgJQIKABASgQUACAlAgoAkBIBBQBI\niYACAKREQAEAUiKgAAApEVAAgJQIKABASgQUACAlAgoAkBIBBQBIiYACAKREQAEAUiKgAAApEVAA\ngJQIKABASgQUACAlAgoAkBIBBQBIiYACAKREQAEAUiKgAAApEVAAgJQIKABASgQUACAlAgojsnjx\n4n5XIQ3aYgPaYgPaYvQIKIwIP74NaIsNaIsNaIvRI6AAACkRUACAlBwR/a7DJsM2jQUAIxAR7nYY\nAgoAkBKb+AAAKRFQAICUCKgWbB9re6nt+22f3qbMmbX/7bZn9bqOvTJcW9h+X22DO2xfb/s/9KOe\nvdDJfFHLHWR7je0Te1m/XurwNzLb9m2277K9uMdV7JkOfiOvtH217SW1Leb2oZrjzvZ5th+zfecQ\nZbpbbkYEj4aHpEmSHpA0XdIUSUskvb6pzB9KurK+PkTSjf2udx/b4k2SXlFfH7s5t0VDuR9JukLS\nu/pd7z7OF9tL+pmk3ev7V/a73n1si/mSvjTQDpKekjS533Ufh7Y4TNIsSXe26d/1cpM1qMEOlvRA\nRCyPiJckfUfSO5vKvEPS+ZIUETdJ2t72Tr2tZk8M2xYRcUNEPFPf3iRp9x7XsVc6mS8k6TRJ/yLp\niV5Wrsc6aYv3SvpuRPxKkiLiyR7XsVc6aYtHJW1XX28n6amIWNPDOvZERPxY0oohinS93CSgBttN\n0kMN739Vuw1XZiIumDtpi0Z/LunKca1R/wzbFrZ3U1k4nVM7TdRDZDuZL/aStIPtRbZvsX1yz2rX\nW520xTclvcH2I5Jul/SJHtUtm66Xm5PHtTqbpk4XKs3H9E/EhVHHn8n2kZJOlfT741edvuqkLb4q\n6b9FRNi2Bs8jE0UnbTFF0v6S3irp5ZJusH1jRNw/rjXrvU7a4jOSlkTEbNuvlvRD2/tGxKpxrltG\nXS03CajBHpa0R8P7PVSSfqgyu9duE00nbaF6YMQ3JR0bEUOt4m/KOmmLAyR9p2STXinpONsvRcTl\nvaliz3TSFg9JejIiVktabfvfJe0raaIFVCdt8WZJX5CkiPi57WWSXifplp7UMI+ul5ts4hvsFkl7\n2Z5ue0tJJ0lqXsBcLukUSbJ9qKSnI+Kx3lazJ4ZtC9vTJF0iaU5EPNCHOvbKsG0RETMjYkZEzFDZ\nD/XRCRhOUme/kcskvcX2JNsvV9kpfneP69kLnbTFUklHS1Ld5/I6Sb/oaS1z6Hq5yRpUk4hYY/tj\nkr6vcoTOuRFxj+0P1/7fiIgrbf+h7QckPSfpA32s8rjppC0kfU7SVEnn1DWHlyLi4H7Vebx02Bab\nhQ5/I0ttXy3pDknrJH0zIiZcQHU4X3xR0kLbt6usFHw6In7Tt0qPE9sXSTpC0ittPyRpnsqm3hEv\nN7nUEQAgJTbxAQBSIqAAACkRUACAlAgoAEBKBBQAICUCCgCQEgGFVGzPt73O9hH9rsto2N6u3lpg\neb31xrpN/VYkDd/N4f2uSwYTZV7NjICaQOqPZajH+xPUce4wdYmGx6bsy5I+pnJx0C+o3HJhyLPm\na5gtG6sK2F5se91YjW9zsxnNq2lxJYmJJyQtaNPvtl5WZBjtftRnS7pIG1/1eFP0R5LujYhWt+QY\nylgv7Fh4jt5En1fTIqAmoIj4637XoQMtr/QdEU+p3NBtU7erpMX9rgTGxESfV/Pq910YeYzpHS3X\nSVo7TJnptdxCSa+VdLGkxyWtlXRELXOApL9X2Tz1lKTVku6TdIak7YcY90mSrpX0mzrMMknflnRA\n7b+4TrvVY1otM7++P7zF+N8q6eo6/hck3SvpS5K2a1F2YFqTVG53cH8d5kFJfytpSpdtu4ukr0la\nLunF2mbflbR/m+k2PxYNMe7ZQ7TLwm7boOE7HrIeko6U9L9ULuL6jKTnJd2pcn3Fl7WoZ9vvZojP\n9hpJ/6xyI7tnJV2vcmfVuXVc728xD7dsK0n/2DivNPU7ROUCvb+u38+Dkv5B0i4tys6sn/uB+pmf\nUrlm4DmSdtgU5lWVu9f+m8qV019QuSniDZI+18tlzng/WIPafL1a0o0qP5x/krS1ykJKkj4o6QSV\nH84PVPZVHijpL1RuIXFIRDw7MKJ676OFKlcqfkIb7ii7h8rCd6mkn9YyK1Ru6vevKrfHHvCMhlAv\nvnmOpFUqC7zHVRawp0s63vbvx4Y7+za6SNJbVG6kuFLS2yV9WtKrVO5fNSzbMyRdpxJS10q6UNI0\nSe+R9Hbb74qI79XiCyUtUrlQ5nKVharq63aWqWyW/WR9/5WGfuvbqIs2WFHHN1fSnioL0gGN9fi0\nypW1f6KysNtKpa3mS5pt++iIGPE+LNt7qSw0d1Bp/yUqNzL8V0lX1WKtNp8NtVlyUD/bp6oEzmqV\nK2Y/pPLn6z+qtMuhEfFQLbuLpJslbSvpeyrtuJVKaM2RdJZKqKSdV20fW+v+dP28D6u08d6SPipp\nU9iC0pl+JySPsXtowz+8eSoLmcbH+2uZ6Q3l/qbNeKapXki4qfupdbhPN3X/UO1+o6Rtm/ptIWnn\nhvdza9lT2kx7vpr+laosZF9U+UG+tqn812r5bzR1X1y736yGtT6Vm+fdL2mNpJ06bNfv13H996bu\nb5L0kqQnJW3T4rv4UZff33JJv2jTb6Rt0HaNWtKMNt3/uo7vT4b7bob5PD+o5U9r6v6OhnnwlKZ+\nbdtNLdagVILotypr+Ls0lT+qfs+XNHQ7rVWdar+tJW2VfV5VWXNfJ2mfFnXaoZt5LvuDo/gmpnkq\nm2kaH81HIv1abQ6miIgHo87tTRaq/Ct8W1P301T+2X44mu4SGhHrIuLXXX+Cjc1RuWz/2RFxX1O/\nz6psOppT78fT7PSIeLqhPs+rrAFtobIpc0i2d5d0jKRfqhyZt15E3KDyr3cHSSd2/GlGZjRt0FJE\ntDti8Kv1ufl77lhtt6NV7nt0dtN0L5f0f0c67iYfVdmX/omIeLRpOj9SWTM83vY2TcO90DyiiFgd\nEYO6d6mX82qrzzChbuNBQE08ERFbtHgc1VTu9oh4qdUIbE+x/THb19n+zcB5PCr7qbaVtFtD2W0k\nvUHSYxFx+zh9pv3r84+ae9Qf9G0q/35/r7m3Wt+1dOCOp1M7mPas+vzjiFjbov9AnfbrYFyjMdI2\naMv2NrY/Y/tm28/YXlu/5ydrkd2GGn4YA+12XZs/O2MVUG+qz7PreUkbPVQ2j01S2ZQplRspPivp\na7b/xfaHbO89RnWRejOvfqs+32T7HNsn1T8EEw77oDZfQ63VXKyyD+rnki7Vhh3PVtlP8rKGstvX\n5/G85f0r6vOjbfo/2lRuvYhY2aL8mvo8aQymPdCO27fpP1ZG3Aat2J6ishA9SOXAiItU9hu+pPI9\nz9PG33O3BurR7tyv0a5VD9ixPv/XIcqEpG2ksnXA9sEqm+eOVV3zrTfYOyMizhplfcZ9Xo2IS23/\nkaRPqWx2/7Ak2f6pymboa0ZQ75QIqM1Xyx3Rtg9UCacfSjouGnaS14MhTm8aZGCTxGj+bQ9nYIfy\nLpLuadF/l6Zy4zHtndv0H89pt6rHWLXBO1XCaWFE/Hljj3ogwbyRVLLBQD12atO/XXtK7ZdLrf4E\nPKMyL78iGg7cGUpELJX0p7YnSdpXZVPkaZL+3vZzEXFeJ+NpoyfzakRcKelK21tLOlTlvLuPSrrC\n9qyIaDXtTQ6b+NDsNfX58hh8BNchKkc8rRcRz0m6S9LOtjvZzDWwmayTtZcBt9bn2c09bG+vsnlt\ntVovEEZrYNpvqQu0Zkc2lRuNtWrfLiNpg7W1f6vzeAa+50ta9BuLS/c0tlur5czsNsOtUDn6cyO1\n7ffT4D9WN6is8XV9+aWIWBsRt0bElyX9We3ceGJ1+nm17jdbFBGfUrm1/JaSjhuLcWdAQKHZwI7z\nIxs72n6VylFIrZxZn79he7um4baw3fhveeDExj27qNO3VDY9nWb71U39Pq+yX+xb7fapjUZEPKyy\nNjlDGw4DlyTZPkTSe1UOS750DCb3lKRX2d6qRb+RtMFTKgvvVm3d7nueKenvRlD3jTS020yVSz41\nTuOdah8oN0na0/YxTd3/SuXo0mZnq7TLV+ph7RuxvaXtwxre72+71WbQgXn0+YZuKedV24e3+bM0\n8BmeG+m4s2ETH5rdrHIy5Ym2r6+vd1LZXr9U0iNqOrM+Iv53XQicLOl+25er7M/YVWUBeK42nJvx\nE5WFwCdt76gN+yjObLMNXhHxS9ufVAnIW23/H5Ud+UeobN64R4M3Paq5nqPwEZV2+B+236ZyTtce\nKudBrZH0gbomOVrXqJxvdrXtH6vs91sSEVeMsA2ukfRuSZfYvkrln/vyiPiWytFtD0j6C9v7qJzn\nM03l3JsrJP3pGHye/6yyhvPV2m53qKy5nVCnf3yLYc6Q9AeSLrN9scoa1ZtVTo9YrKY1k4i4t54H\ndZ6kn9m+WuXQ7Cn18xymMo8NHAhxiqQP2b5O5QjDFSrnBB6vclTcVxtGn3VePVPSrvX3+UuVw+wP\nUPmtLZf0nS7GlVu/j3PnMXYPdXclifOGKDNV5Qe2TGWhdr+kv1E5+miZ2p+r816VhcjTdbifq5wE\nvF9TuT9Q+fGvGqizNpydP6++b3V2/jEq5yQNnJ1/n8qZ9q3Ozl/Uri1UDrlfqzbnt7QZZldJX9fG\nV5K4RPUqGW2+i27Pg3p5ncZDKv/C1zZ/T122wRYqF6r9ucpCbKM6Sdpd5R//r7ThKhJ/qbJJa1D9\nh/puhvhMr9bgK0kcV7+DlucYqYTFzXUeekLlaiR7qJzmsH5eaRrmjbX/8touT2rD1SFmN5Q7uLbx\nEpU1pOdrG54rae8W4003r6r8Mfp2Hecqlf1Zd6isoe04mmVItofrBwaAnrE9V2WtZ25EXNDn6iAp\n9kEBAFIioAAAKRFQAPohxL2qMAz2QQEAUmINCgCQEgEFAEiJgAIApERAAQBSIqAAACkRUACAlP4/\ns1/oUxX8whMAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x10a23e250>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "cdict = {}\n",
    "for k,v in zip(c['block'].values,c['gender'].values):\n",
    "    cdict[k]=v\n",
    "\n",
    "malechair_maleafter = ['M'*x[1]['questions'].count('M') for x in q.iterrows() if cdict[int(str(x[1]['session']).split('.')[0])] == 'M']\n",
    "malechair_femaleafter = ['F'*x[1]['questions'].count('F') for x in q.iterrows() if cdict[int(str(x[1]['session']).split('.')[0])] == 'M']\n",
    "femalechair_maleafter = ['M'*x[1]['questions'].count('M') for x in q.iterrows() if cdict[int(str(x[1]['session']).split('.')[0])] == 'F']\n",
    "femalechair_femaleafter = [\"F\"*x[1]['questions'].count('F') for x in q.iterrows() if cdict[int(str(x[1]['session']).split('.')[0])] == 'F']\n",
    "vc_malechair = pd.value_counts(list(''.join(malechair_maleafter+malechair_femaleafter)))\n",
    "vc_femalechair = pd.value_counts(list(''.join(femalechair_maleafter+femalechair_femaleafter)))\n",
    "\n",
    "# Load everything into a single dataframe\n",
    "\n",
    "chairdata = [vc_malechair,vc_femalechair]\n",
    "chairlabels = ['Male chair','Female chair']\n",
    "\n",
    "chairnormdata = [x/x.sum() for x in chairdata]\n",
    "df = pd.DataFrame(chairnormdata,index=chairlabels)\n",
    "\n",
    "print df\n",
    "\n",
    "# Plot stacked bar chart\n",
    "fig5 = plt.figure(3,(6,6))\n",
    "ax = fig5.add_subplot(111)\n",
    "dfplot = df.plot(kind='barh',stacked=True,ax=ax,color=gencolors,legend=True)\n",
    "\n",
    "ax.set_xlabel('Fraction of total questions',fontsize=20)\n",
    "\n",
    "# Check statistical significance\n",
    "\n",
    "p1exp,p2exp = df['M']\n",
    "n1 = len(list(''.join(malechair_maleafter+malechair_femaleafter)))\n",
    "n2 = len(list(''.join(femalechair_maleafter+femalechair_femaleafter)))\n",
    "z = ztest(p1exp,p2exp,n1,n2);"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "When looking at the gender of the session chair (who might have selection biases when there are more questions than time permits), there is no difference in the gender distribution of questioners. In this case the data are nearly identical for male vs. female chairs, at 36% and 35% females respectively. \n",
    "\n",
    "This null result differs from the NAM data, who saw a small but significant decrease in the fraction of females asking questions when the chair was male (16%) as opposed to a female chair (22%). "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## (Some) conclusions"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "I'm not sure what my expectations were in looking at data for a smaller, specialized meeting vs. a large general meeting like the AAS or NAM. I don't think I would have been hugely surprised by either significant differences or none at all --- I could easily have guessed that a smaller, specialized meeting is more subject to the whims of an SOC (who are more likely to be senior and thus might tend to invite senior colleagues, which would presumably skew the gender ratio towards more males). On the other hand, postdocs and early-career researchers might travel more (without teaching responsibilities mid-semester) and be more invested in giving talks to establish themselves for permanent jobs, which would skew toward younger people and a higher female-male ratio. The demographics of the AGN/galaxy community, the relatively remote conference location in South America, and the timing of the meeting are all priors that should be taken into account in a full analysis.\n",
    "\n",
    "One extremely important point is that this is a [sample size](https://youtu.be/dw9qqvm-LT8) of $N=1$. We need lots more data on other conferences, in different subfields, different parts of the world, and different conference structures. Unlike the large meetings, though, data collection can pretty easily be done by one person, and we're starting to put together some standard graphs and questions we can ask through open data sets [(all of this data is available on Github)](https://github.com/willettk/chile2015-gender) and joint analysis tools like the iPython notebook. I hope it's easy to clone this work for future conferences; I'm happy to help anyone who wants to try that. \n",
    "\n",
    "Given that limited sample size, I think the results for <span style=\"color:red\">AGN2015</span> were really encouraging. I found:\n",
    "\n",
    "* both the attendees and questioners had significantly higher female/male ratios than the NAM or AAS meetings\n",
    "* the male/female ratio was similar for both attendees and speakers, meaning that there was no apparent preference towards which gender got the opportunity to give a talk\n",
    "* both male and female speakers received about the same number of questions per talk\n",
    "* the gender of the first person to ask a question didn't affect subsequent questions\n",
    "* the gender of the session chair didn't affect subsequent questions\n",
    "\n",
    "The main (possible) negative was that the session chairs were overwhelmingly male (again, a possible reflection of a seniority bias). However, this is presumably simpler to correct for an SOC/LOC, is subject to small sample size variance, and doesn't appear to have affected the engagement with the audience, which is good. I also should note for <span style=\"color:red\">AGN2015</span> that the three end-of-day discussions were all moderated by women, and the conference summary was also delivered by a woman. My unofficial impression is that this conference probably had a younger-than-average age, especially among speakers. The targeted talks at this conference were unofficially designated as being reserved for \"*young people doing cool things*\" by the SOC. I'd love to see more of that in every conference, since I think it'd help with both scientific output and participant balance, but those are my own personal biases in play. \n",
    "\n",
    "So maybe smaller conferences are better at gender balance and engagement than larger meetings? As always, we need **lots more data** to see if the effect is real, and to think about what the community can do as a result. \n",
    "\n",
    "*Thanks to the conference organizers of <span style=\"color:red\">AGN2015</span> for inviting me to the meeting, to Meg Urry for presenting some of these results in her conference summary, and to James Davenport, Jonathan Pritchard, and Karen Masters for their previous hard work at the AAS and NAM meetings.*"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "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.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}