CouldYouLookAtThis.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "After comparing Keeling's Matlab code and output, and using multiple print statements throughout my code, I am still not sure why the size of the nearest infected farm is not significant. I found a couple of bugs, but they did not change the results. I also remembered that the reason I had previously decided not to use Keeling's kernel and use a Cauchy kernel was because Keeling's kernel makes the disease progress very slowly. So I have used the Cauchy kernel in the code below. The changes that I have made result in output that is consistent with Keeling's code, however I have noticed that in this new version the number of new infected is dramatically lower than it used to be."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "from random import randint\n",
    "Size = 20\n",
    "N = 75 \n",
    "np.random.seed(53080)\n",
    "x = Size*np.random.rand(N)\n",
    "np.random.seed(23003)\n",
    "y = Size*np.random.rand(N)  \n",
    "#np.random.seed(10)\n",
    "#Cows = np.array([randint(25,250) for p in range(N)])\n",
    "Cows = np.array([randint(25,51) for p in range(15)]+[randint(51,76) for p in range(30)]+[randint(76,95) for p in range(20)]+[randint(95,250) for p in range(10)])\n",
    "#np.random.seed(11)\n",
    "#Sheep = np.array([randint(25,250) for p in range(N)])\n",
    "Sheep = np.array([randint(25,51) for p in range(15)]+[randint(51,76) for p in range(30)]+[randint(76,95) for p in range(20)]+[randint(95,250) for p in range(10)])\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "#Calculates which grid square a particular location is in (turn a 2-d coordinate into a scalar)\n",
    "def WhichGrid(x,y,XRange,YRange,XNum,YNum):\n",
    "    #Essentially: floor(Unif[0,1)griddim)griddim+floor(Unif[0,1)griddim)+1\n",
    "    #Returns a number from 1 to griddim^2\n",
    "    return(np.floor(x*(XNum/XRange))*YNum+np.floor(y*(YNum/YRange))+1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def Kernel(dist_squared):\n",
    "    dist_squared = np.asarray(dist_squared)\n",
    "    is_scalar = False if dist_squared.ndim > 0 else True\n",
    "    dist_squared.shape = (1,)*(1-dist_squared.ndim) + dist_squared.shape\n",
    "    K = 1 / (pi * (1 + dist_squared**2))\n",
    "    K[(dist_squared < 0.0138)] = 0.3093\n",
    "    K[(dist_squared > 60*60)] = 0\n",
    "    return(K if not is_scalar else K[0])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "from math import pi\n",
    "def Iterate(Status, x, y, Suscept, Transmiss, grid, first_in_grid, last_in_grid, Num, MaxRate):\n",
    "    Event = 0*Status\n",
    "    INF = np.where(Status>5)[0]\n",
    "    NI = INF.size # Note reported farms still infectious\n",
    "    IGrids = grid[INF]-1\n",
    "        \n",
    "    for ii in range(NI):\n",
    "        INFi = INF[ii]\n",
    "        trans = np.multiply(-Transmiss[INFi],Num) #transmissibility of infected farm to all other grid squares \n",
    "        maxr = MaxRate[IGrids[ii],:] #max number of animals to be infected in infected grid square\n",
    "        # Elementwise multiplication\n",
    "        rate = np.multiply(trans, maxr) #max number of animals to be infected in each grid square based on infected grid square\n",
    "        MaxProb = 1 - np.exp(rate) #Max probability that infected farm infected noninfected farm\n",
    "        rng = np.random.rand(len(MaxProb))\n",
    "        m = np.where((MaxProb - rng)>0)[0]  #these grid squares need further consideration\n",
    "        for n in range(len(m)):\n",
    "            s = 1\n",
    "            M = m[n]\n",
    "            PAB = 1 - np.exp(-Transmiss[INFi]*MaxRate[IGrids[ii],M]) #Max probability that infected farm infects noninfected farms under consideration\n",
    "            if (PAB == 1):\n",
    "                # Calculate the infection probability for each farm in the susceptible grid\n",
    "                leng = last_in_grid[M]-first_in_grid[M]+1\n",
    "                R = np.random.rand(leng)\n",
    "                for j in range(leng):\n",
    "                    ind1 = first_in_grid[M]+j-1\n",
    "                    Q = 1 - np.exp(-Transmiss[INFi]*Suscept[ind1]*Kernel((x[INFi]-x[ind1])**2+(y[INFi]-y[ind1])**2))\n",
    "                    if ((R[j] < Q) & (Status[ind1] == 0)):\n",
    "                        Event[ind1] = 1\n",
    "            else:\n",
    "                R = np.random.rand(Num[M])\n",
    "                # Loop through all susceptible farms in the grids where an infection event occurred.  \n",
    "                for j in range(Num[M]):\n",
    "                    P = 1 - s*(1 - PAB)**(Num[M] - j)\n",
    "                    if (R[j] < (PAB / P)):\n",
    "                        s = 0\n",
    "                        ind1=first_in_grid[M]+j-1\n",
    "                        Q=1-np.exp(-Transmiss[INFi]*Suscept[ind1]*Kernel((x[INFi]-x[ind1])**2+(y[INFi]-y[ind1])**2))\n",
    "                        if ((R[j]< Q/P) & (Status[ind1] == 0)):\n",
    "                            Event[ind1] = 1\n",
    "    # Evolve the infection process of those farms which have been exposed and already infectious ones.  \n",
    "    Status[Status > 0] += 1\n",
    "    Status = Status + Event\n",
    "    #m=np.where(Status==13); # Initiate Ring Culling Around Reported Farm\n",
    "    #for i in range(len(m)):\n",
    "    #    Status[m[i]]=-1;\n",
    "    return {'Status':Status,'NI':NI}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/sandyalakkur/anaconda/envs/py3k/lib/python3.4/site-packages/IPython/kernel/__main__.py:56: DeprecationWarning: using a non-integer number instead of an integer will result in an error in the future\n",
      "/Users/sandyalakkur/anaconda/envs/py3k/lib/python3.4/site-packages/IPython/kernel/__main__.py:11: DeprecationWarning: using a non-integer number instead of an integer will result in an error in the future\n",
      "/Users/sandyalakkur/anaconda/envs/py3k/lib/python3.4/site-packages/IPython/kernel/__main__.py:20: DeprecationWarning: using a non-integer number instead of an integer will result in an error in the future\n",
      "/Users/sandyalakkur/anaconda/envs/py3k/lib/python3.4/site-packages/IPython/kernel/__main__.py:11: DeprecationWarning: using a non-integer number instead of an integer will result in an error in the future\n",
      "/Users/sandyalakkur/anaconda/envs/py3k/lib/python3.4/site-packages/IPython/kernel/__main__.py:20: DeprecationWarning: using a non-integer number instead of an integer will result in an error in the future\n",
      "/Users/sandyalakkur/anaconda/envs/py3k/lib/python3.4/site-packages/IPython/kernel/__main__.py:108: DeprecationWarning: converting an array with ndim > 0 to an index will result in an error in the future\n",
      "/Users/sandyalakkur/anaconda/envs/py3k/lib/python3.4/site-packages/IPython/kernel/__main__.py:160: DeprecationWarning: using a non-integer number instead of an integer will result in an error in the future\n"
     ]
    }
   ],
   "source": [
    "def Outbreaks(Size,N,Y0,farms,end,end2,x,y,Cows,Sheep,Maxtime=1000):\n",
    "    #This is an attempt of converting the Matlab Program 7.6 Code into Python\n",
    "    import numpy as np\n",
    "    import pandas as pd\n",
    "    from math import pi\n",
    "    \n",
    "    Status = np.array([0]*N)    #Initial Status of each farm\n",
    "    init_ind = np.random.randint(0,N)\n",
    "    for i in range(Y0):\n",
    "        Status[init_ind] = 6 #one farm is initially infected \n",
    "\n",
    "    #Cows are 10.5 times more susceptible to disease than sheep\n",
    "    Suscept = Sheep+10.5*Cows\n",
    "    Transmiss = 5.1e-7*Sheep + 7.7e-7*Cows \n",
    "\n",
    "    #Set up the grid\n",
    "    grid = WhichGrid(x,y,Size,Size,10.0,10.0)\n",
    "    tmp = sorted(grid) #Sort grid values\n",
    "    #i = np.argsort(grid) #get indexed values after sort\n",
    "    i = [i[0] for i in sorted(enumerate(grid), key=lambda x:x[1])]\n",
    "    x = x[i]\n",
    "    y = y[i]\n",
    "    Status = Status[i]\n",
    "    grid = grid[i]\n",
    "    Transmiss = Transmiss[i]\n",
    "    Suscept = Suscept[i]\n",
    "    Cows = Cows[i]\n",
    "    Sheep = Sheep[i]\n",
    "    Xgrid = []\n",
    "    Ygrid = []\n",
    "    Num = []\n",
    "    first_in_grid = []\n",
    "    last_in_grid = []\n",
    "    Max_Sus_grid = []\n",
    "    index_inf = np.where(Status==6)[0].astype(int)\n",
    "    m2 = np.array(np.where(grid==1))\n",
    "    \n",
    "    for i in range(1,int(max(grid))+1):\n",
    "        #turn the grid square number into an x-coordinate and y-coordinate (should not exceed XNum)\n",
    "        Xgrid.append(np.floor((i-1)/10))\n",
    "        Ygrid.append((i-1)%10)\n",
    "        m = np.array(np.where(grid==i))\n",
    "        Num.append(m.shape[1])\n",
    "    \n",
    "        if Num[i-1] > 0:\n",
    "            first_in_grid.append(m.min()) #Add the \"+1\" here so the indicies match those in the Keeling code\n",
    "            last_in_grid.append(m.max())\n",
    "            Max_Sus_grid.append(Suscept[m].max())\n",
    "        else:\n",
    "            first_in_grid.append(0)\n",
    "            last_in_grid.append(-1)\n",
    "            Max_Sus_grid.append(0)\n",
    "\n",
    "    #Work out grid to maximum grid transmission probabilities\n",
    "    from numpy import ndarray\n",
    "    MaxRate = ndarray((max(grid),max(grid)))\n",
    "\n",
    "    #Determine maximum number of animals to be infected in each grid square\n",
    "\n",
    "    for i in range (1,int(max(grid))):\n",
    "       for j in range(1,int(max(grid))):\n",
    "           if ((i-1)==(j-1)) | (Num[i-1]==0) | (Num[j-1] == 0):\n",
    "               MaxRate[i-1,j-1] = np.inf\n",
    "           else:\n",
    "               Dist2 = (Size*max([0,(abs(Xgrid[i-1]-Xgrid[j-1])-1)])/10)**2+(Size*max([0,(abs(Ygrid[i-1]-Ygrid[j-1])-1)])/10)**2\n",
    "               MaxRate[i-1,j-1] = Max_Sus_grid[j-1]*Kernel(Dist2)\n",
    "\n",
    "    #Susceptible, Exposed, Infectious, Reported.==> latent period is 4 days\n",
    "    i=1; S=len(np.where(Status==0)); E=len(np.where(np.logical_and(Status>0, Status<=5)));I=len(np.where(np.logical_and(Status>5, Status<=9))); R=len(np.where(Status==10)); R2=len(np.where(Status>9)); CullSheep=0; CullCattle=0;\n",
    "    i=i+1;  IterateFlag=1;\n",
    "\n",
    "    \n",
    "    S=[]\n",
    "    E=[]\n",
    "    I=[]\n",
    "    R=[]\n",
    "    R2=[]\n",
    "    CullSheep=[]\n",
    "    CullCattle=[]\n",
    "    t=[]\n",
    "    t.append(0)\n",
    "    results = np.c_[np.array([1]*N),np.arange(1,N+1),np.array([0]*N)]\n",
    "    \n",
    "    while(np.logical_and(t[-1]<end, IterateFlag)):\n",
    "        Status=Iterate(Status, x, y, Suscept, Transmiss, grid, first_in_grid, last_in_grid, Num, MaxRate)['Status']\n",
    "        Sus=np.where(Status==0)[0]; Exp=np.where(np.logical_and(Status>0, Status<=5))[0]; Inf=np.where(Status>5)[0]; \n",
    "        S.append(len(Sus)); E.append(len(Exp)); I.append(len(Inf)); \n",
    "        t.append(t[i-2]+1);i+=1;\n",
    "       \n",
    "        #This is how I stop the simulation (all farms are infected)\n",
    "        if t[-1]>5:\n",
    "            if np.logical_or((E[-4]+I[-4]==0),I == N):\n",
    "                IterateFlag=0\n",
    "        from scipy.stats import itemfreq\n",
    "        sim_num = np.array([i-1]*N)\n",
    "        seq = np.arange(1,N+1)\n",
    "        results_full = np.r_[results,np.c_[sim_num,seq,Status]]\n",
    "        results = results_full\n",
    "   \n",
    "    #Return information regarding only farm of interest\n",
    "    this = results_full[np.logical_or.reduce([results_full[:,1] == x for x in farms])]\n",
    "    #Extract rows relating to timepoint of interest\n",
    "    no_this = this[this[:,0]==end]\n",
    "    #turn status to an indicator\n",
    "    Status_ind = (no_this[:,2]>5).astype(int)\n",
    "    #Calculate distance to index farm - first infected is first in list of coords\n",
    "    coords = list(zip(x,y))\n",
    "    index = np.array((coords[index_inf][0],coords[index_inf][1]))\n",
    "    dist = []\n",
    "    for j in range(0,N):\n",
    "        b = np.array((coords[j][0],coords[j][1]))\n",
    "        dist.append(np.linalg.norm(b-index))\n",
    "    to_return = np.c_[no_this[:,1],Status_ind,dist,Cows,Sheep,x,y]\n",
    "    \n",
    "    from scipy import spatial\n",
    "    #Extract the infected farms\n",
    "    inf_farms = to_return[to_return[:,1]==1]\n",
    "    coords = list(zip(x,y))\n",
    "    #Create list of coordinates infected farms\n",
    "    inf_farm_coords = list(zip(inf_farms[:,5],inf_farms[:,6]))\n",
    "    list_of_inf_coords = [list(elem) for elem in inf_farm_coords]\n",
    "    #Create list of coordinates of all farms\n",
    "    list_of_coords = [list(elem) for elem in coords]\n",
    "    #Calculate Euclidean distance from each farm to all infected farms- each row in matrix represents\n",
    "        #distance of one farm to each infected farm\n",
    "    dist_to_inf = spatial.distance_matrix(list_of_coords,list_of_inf_coords)\n",
    "    #Find distance to closest infected farm\n",
    "    def minval(array):\n",
    "        #return(np.min(array[np.nonzero(array)]))\n",
    "        return(np.min(array))\n",
    "    closest_infected = np.apply_along_axis(minval,1,dist_to_inf)\n",
    "    average_infected = np.apply_along_axis(np.mean,1,dist_to_inf)\n",
    "    to_return = np.c_[to_return,closest_infected,average_infected]\n",
    "    \n",
    "    #Create list of number of Cows and number of sheep for infected farms\n",
    "    inf_farm_cows = list(inf_farms[:,3])\n",
    "    inf_farm_sheep = list(inf_farms[:,4])\n",
    "    #Create a function that extracts farm size based on closest infected farm\n",
    "    def where_minval(array):\n",
    "        #return(np.argmin(array[np.nonzero(array)]))\n",
    "        return(np.argmin(array))\n",
    "    closest_infected_size_ind = np.apply_along_axis(where_minval,1,dist_to_inf)\n",
    "    closest_infected_cows = [inf_farm_cows[i] for i in closest_infected_size_ind]\n",
    "    closest_infected_sheep = [inf_farm_sheep[i] for i in closest_infected_size_ind]\n",
    "\n",
    "    #Returns array: farmID, Status_ind, dist_to_index, num_Cows,num_Sheep,x,y,,disttoclosestinf,avgdisttoinf,#cowsinclosestiffarm,#sheepinclosestinffarm\n",
    "    to_return = np.c_[to_return,closest_infected_cows,closest_infected_sheep]\n",
    "    \n",
    "    #Now run the outbreak for the additional end - end2 steps\n",
    "    while(np.logical_and(t[-1]<end2, IterateFlag)):\n",
    "        Status=Iterate(Status, x, y, Suscept, Transmiss, grid, first_in_grid, last_in_grid, Num, MaxRate)['Status']\n",
    "        Sus=np.where(Status==0)[0]; Exp=np.where(np.logical_and(Status>0, Status<=5))[0]; Inf=np.where(Status>5)[0]; \n",
    "    \n",
    "        S.append(len(Sus)); E.append(len(Exp)); I.append(len(Inf)); \n",
    "        t.append(t[i-2]+1);i+=1;\n",
    "    \n",
    "    #Return only the data where a farm was not infected by the \"end\" day\n",
    "    not_inf = to_return[to_return[:,1]==0]\n",
    "    #append the status of these farms after the additional end2 - end days\n",
    "    newstatus = [Status[i] for i in not_inf[:,0]-1]\n",
    "    newstatus_ind = np.array([i>5 for i in newstatus]).astype(int)\n",
    "    final = np.c_[not_inf,newstatus_ind]\n",
    "    return(final)\n",
    "    \n",
    "test = Outbreaks(Size=Size,N=N,Y0=1,farms = np.arange(1,N+1),end=10,end2=50,x=x,y=y,Cows=Cows,Sheep=Sheep)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([ 0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,\n",
       "        0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,\n",
       "        0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,\n",
       "        0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,\n",
       "        0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  1.,  0.,  0.,  0.,\n",
       "        0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.])"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "test[:,11]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/sandyalakkur/anaconda/envs/py3k/lib/python3.4/site-packages/IPython/kernel/__main__.py:56: DeprecationWarning: using a non-integer number instead of an integer will result in an error in the future\n",
      "/Users/sandyalakkur/anaconda/envs/py3k/lib/python3.4/site-packages/IPython/kernel/__main__.py:11: DeprecationWarning: using a non-integer number instead of an integer will result in an error in the future\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0 3\n",
      "1 0\n",
      "2 0\n",
      "3 0\n",
      "4 0\n",
      "5 0\n",
      "6 0\n",
      "7 1\n",
      "8 0\n",
      "9 0\n",
      "10 1\n",
      "11 1\n",
      "12 0\n",
      "13 1\n",
      "14 1\n",
      "15 0\n",
      "16 1\n",
      "17 1\n",
      "18 0\n",
      "19 0\n",
      "20 1\n",
      "21 0\n",
      "22 2\n",
      "23 0\n",
      "24 0\n",
      "25 1\n",
      "26 0\n",
      "27 1\n",
      "28 1\n",
      "29 0\n",
      "30 0\n",
      "31 0\n",
      "32 0\n",
      "33 3\n",
      "34 3\n",
      "35 0\n",
      "36 0\n",
      "37 0\n",
      "38 0\n",
      "39 0\n",
      "40 2\n",
      "41 2\n",
      "42 0\n",
      "43 0\n",
      "44 0\n",
      "45 0\n",
      "46 0\n",
      "47 0\n",
      "48 0\n",
      "49 0\n",
      "50 0\n",
      "51 1\n",
      "52 0\n",
      "53 0\n",
      "54 1\n",
      "55 2\n",
      "56 3\n",
      "57 1\n",
      "58 2\n",
      "59 0\n",
      "60 0\n",
      "61 0\n",
      "62 2\n",
      "63 0\n",
      "64 2\n",
      "65 0\n",
      "66 1\n",
      "67 0\n",
      "68 1\n",
      "69 2\n",
      "70 0\n",
      "71 0\n",
      "72 0\n",
      "73 0\n",
      "74 1\n",
      "75 0\n",
      "76 1\n",
      "77 1\n",
      "78 0\n",
      "79 0\n",
      "80 0\n",
      "81 0\n",
      "82 0\n",
      "83 0\n",
      "84 0\n",
      "85 2\n",
      "86 0\n",
      "87 0\n",
      "88 0\n",
      "89 0\n",
      "90 0\n",
      "91 0\n",
      "92 0\n",
      "93 3\n",
      "94 0\n",
      "95 0\n",
      "96 0\n",
      "97 1\n",
      "98 1\n",
      "99 1\n",
      "100 0\n",
      "101 0\n",
      "102 0\n",
      "103 0\n",
      "104 1\n",
      "105 0\n",
      "106 0\n",
      "107 0\n",
      "108 0\n",
      "109 1\n",
      "110 2\n",
      "111 0\n",
      "112 1\n",
      "113 1\n",
      "114 0\n",
      "115 2\n",
      "116 0\n",
      "117 0\n",
      "118 1\n",
      "119 0\n",
      "120 0\n",
      "121 0\n",
      "122 0\n",
      "123 0\n",
      "124 1\n",
      "125 1\n",
      "126 2\n",
      "127 1\n",
      "128 0\n",
      "129 1\n",
      "130 0\n",
      "131 0\n",
      "132 0\n",
      "133 0\n",
      "134 0\n",
      "135 0\n",
      "136 0\n",
      "137 0\n",
      "138 2\n",
      "139 1\n",
      "140 2\n",
      "141 1\n",
      "142 1\n",
      "143 0\n",
      "144 0\n",
      "145 1\n",
      "146 0\n",
      "147 0\n",
      "148 2\n",
      "149 0\n",
      "150 1\n",
      "151 0\n",
      "152 0\n",
      "153 0\n",
      "154 1\n",
      "155 0\n",
      "156 2\n",
      "157 1\n",
      "158 2\n",
      "159 1\n",
      "160 0\n",
      "161 0\n",
      "162 0\n",
      "163 0\n",
      "164 0\n",
      "165 0\n",
      "166 0\n",
      "167 0\n",
      "168 0\n",
      "169 0\n",
      "170 3\n",
      "171 3\n",
      "172 0\n",
      "173 0\n",
      "174 0\n",
      "175 0\n",
      "176 0\n",
      "177 0\n",
      "178 1\n",
      "179 0\n",
      "180 1\n",
      "181 1\n",
      "182 1\n",
      "183 0\n",
      "184 0\n",
      "185 0\n",
      "186 0\n",
      "187 0\n",
      "188 0\n",
      "189 0\n",
      "190 0\n",
      "191 0\n",
      "192 1\n",
      "193 1\n",
      "194 0\n",
      "195 1\n",
      "196 2\n",
      "197 0\n",
      "198 0\n",
      "199 0\n",
      "200 1\n",
      "201 0\n",
      "202 0\n",
      "203 0\n",
      "204 4\n",
      "205 0\n",
      "206 0\n",
      "207 1\n",
      "208 1\n",
      "209 0\n",
      "210 0\n",
      "211 0\n",
      "212 1\n",
      "213 0\n",
      "214 1\n",
      "215 0\n",
      "216 0\n",
      "217 0\n",
      "218 0\n",
      "219 0\n",
      "220 0\n",
      "221 0\n",
      "222 0\n",
      "223 0\n",
      "224 0\n",
      "225 0\n",
      "226 0\n",
      "227 1\n",
      "228 4\n",
      "229 0\n",
      "230 0\n",
      "231 0\n",
      "232 1\n",
      "233 0\n",
      "234 3\n",
      "235 0\n",
      "236 0\n",
      "237 1\n",
      "238 0\n",
      "239 0\n",
      "240 2\n",
      "241 0\n",
      "242 0\n",
      "243 0\n",
      "244 1\n",
      "245 0\n",
      "246 0\n",
      "247 0\n",
      "248 1\n",
      "249 0\n",
      "250 1\n",
      "251 0\n",
      "252 0\n",
      "253 0\n",
      "254 4\n",
      "255 0\n",
      "256 1\n",
      "257 2\n",
      "258 0\n",
      "259 0\n",
      "260 0\n",
      "261 2\n",
      "262 0\n",
      "263 0\n",
      "264 0\n",
      "265 1\n",
      "266 0\n",
      "267 3\n",
      "268 0\n",
      "269 0\n",
      "270 1\n",
      "271 0\n",
      "272 0\n",
      "273 1\n",
      "274 0\n",
      "275 4\n",
      "276 0\n",
      "277 1\n",
      "278 0\n",
      "279 0\n",
      "280 0\n",
      "281 0\n",
      "282 0\n",
      "283 0\n",
      "284 0\n",
      "285 0\n",
      "286 0\n",
      "287 0\n",
      "288 0\n",
      "289 2\n",
      "290 0\n",
      "291 1\n",
      "292 0\n",
      "293 1\n",
      "294 1\n",
      "295 0\n",
      "296 1\n",
      "297 1\n",
      "298 0\n",
      "299 0\n",
      "300 0\n",
      "301 0\n",
      "302 0\n",
      "303 0\n",
      "304 2\n",
      "305 0\n",
      "306 0\n",
      "307 1\n",
      "308 1\n",
      "309 0\n",
      "310 2\n",
      "311 0\n",
      "312 1\n",
      "313 1\n",
      "314 0\n",
      "315 1\n",
      "316 0\n",
      "317 2\n",
      "318 1\n",
      "319 2\n",
      "320 2\n",
      "321 1\n",
      "322 0\n",
      "323 0\n",
      "324 2\n",
      "325 2\n",
      "326 0\n",
      "327 1\n",
      "328 0\n",
      "329 0\n",
      "330 2\n",
      "331 0\n",
      "332 1\n",
      "333 1\n",
      "334 0\n",
      "335 1\n",
      "336 0\n",
      "337 1\n",
      "338 2\n",
      "339 0\n",
      "340 1\n",
      "341 0\n",
      "342 0\n",
      "343 0\n",
      "344 0\n",
      "345 0\n",
      "346 2\n",
      "347 0\n",
      "348 0\n",
      "349 0\n",
      "350 1\n",
      "351 0\n",
      "352 1\n",
      "353 0\n",
      "354 1\n",
      "355 0\n",
      "356 0\n",
      "357 3\n",
      "358 0\n",
      "359 0\n",
      "360 1\n",
      "361 0\n",
      "362 0\n",
      "363 0\n",
      "364 0\n",
      "365 0\n",
      "366 0\n",
      "367 0\n",
      "368 1\n",
      "369 0\n",
      "370 0\n",
      "371 0\n",
      "372 0\n",
      "373 2\n",
      "374 0\n",
      "375 2\n",
      "376 2\n",
      "377 0\n",
      "378 0\n",
      "379 0\n",
      "380 1\n",
      "381 0\n",
      "382 0\n",
      "383 1\n",
      "384 0\n",
      "385 0\n",
      "386 0\n",
      "387 4\n",
      "388 1\n",
      "389 0\n",
      "390 0\n",
      "391 2\n",
      "392 0\n",
      "393 1\n",
      "394 0\n",
      "395 2\n",
      "396 0\n",
      "397 0\n",
      "398 1\n",
      "399 0\n",
      "400 1\n",
      "401 0\n",
      "402 0\n",
      "403 0\n",
      "404 1\n",
      "405 0\n",
      "406 0\n",
      "407 0\n",
      "408 0\n",
      "409 0\n",
      "410 0\n",
      "411 2\n",
      "412 0\n",
      "413 0\n",
      "414 2\n",
      "415 1\n",
      "416 0\n",
      "417 0\n",
      "418 0\n",
      "419 2\n",
      "420 0\n",
      "421 0\n",
      "422 0\n",
      "423 1\n",
      "424 0\n",
      "425 1\n",
      "426 1\n",
      "427 0\n",
      "428 2\n",
      "429 0\n",
      "430 1\n",
      "431 1\n",
      "432 2\n",
      "433 0\n",
      "434 2\n",
      "435 0\n",
      "436 1\n",
      "437 0\n",
      "438 0\n",
      "439 1\n",
      "440 0\n",
      "441 0\n",
      "442 0\n",
      "443 0\n",
      "444 1\n",
      "445 0\n",
      "446 0\n",
      "447 0\n",
      "448 0\n",
      "449 0\n",
      "450 0\n",
      "451 0\n",
      "452 1\n",
      "453 0\n",
      "454 0\n",
      "455 0\n",
      "456 0\n",
      "457 0\n",
      "458 0\n",
      "459 3\n",
      "460 0\n",
      "461 1\n",
      "462 1\n",
      "463 0\n",
      "464 0\n",
      "465 0\n",
      "466 0\n",
      "467 0\n",
      "468 1\n",
      "469 0\n",
      "470 1\n",
      "471 0\n",
      "472 0\n",
      "473 0\n",
      "474 0\n",
      "475 0\n",
      "476 1\n",
      "477 0\n",
      "478 0\n",
      "479 0\n",
      "480 0\n",
      "481 2\n",
      "482 0\n",
      "483 2\n",
      "484 0\n",
      "485 0\n",
      "486 1\n",
      "487 0\n",
      "488 0\n",
      "489 0\n",
      "490 0\n",
      "491 1\n",
      "492 0\n",
      "493 0\n",
      "494 0\n",
      "495 0\n",
      "496 0\n",
      "497 0\n",
      "498 1\n",
      "499 1\n",
      "500 0\n",
      "501 2\n",
      "502 0\n",
      "503 4\n",
      "504 1\n",
      "505 1\n",
      "506 0\n",
      "507 0\n",
      "508 0\n",
      "509 1\n",
      "510 1\n",
      "511 0\n",
      "512 2\n",
      "513 2\n",
      "514 0\n",
      "515 1\n",
      "516 1\n",
      "517 0\n",
      "518 0\n",
      "519 0\n",
      "520 1\n",
      "521 0\n",
      "522 1\n",
      "523 0\n",
      "524 0\n",
      "525 1\n",
      "526 2\n",
      "527 0\n",
      "528 0\n",
      "529 1\n",
      "530 0\n",
      "531 0\n",
      "532 0\n",
      "533 0\n",
      "534 0\n",
      "535 0\n",
      "536 1\n",
      "537 0\n",
      "538 0\n",
      "539 0\n",
      "540 0\n",
      "541 0\n",
      "542 1\n",
      "543 0\n",
      "544 1\n",
      "545 0\n",
      "546 0\n",
      "547 0\n",
      "548 0\n",
      "549 2\n",
      "550 4\n",
      "551 1\n",
      "552 1\n",
      "553 0\n",
      "554 1\n",
      "555 1\n",
      "556 2\n",
      "557 2\n",
      "558 1\n",
      "559 0\n",
      "560 0\n",
      "561 0\n",
      "562 0\n",
      "563 0\n",
      "564 1\n",
      "565 0\n",
      "566 0\n",
      "567 0\n",
      "568 0\n",
      "569 1\n",
      "570 0\n",
      "571 1\n",
      "572 1\n",
      "573 0\n",
      "574 0\n",
      "575 0\n",
      "576 0\n",
      "577 0\n",
      "578 0\n",
      "579 0\n",
      "580 0\n",
      "581 2\n",
      "582 0\n",
      "583 0\n",
      "584 0\n",
      "585 0\n",
      "586 0\n",
      "587 0\n",
      "588 1\n",
      "589 0\n",
      "590 0\n",
      "591 0\n",
      "592 0\n",
      "593 1\n",
      "594 0\n",
      "595 0\n",
      "596 0\n",
      "597 0\n",
      "598 0\n",
      "599 3\n",
      "600 1\n",
      "601 2\n",
      "602 1\n",
      "603 0\n",
      "604 0\n",
      "605 1\n",
      "606 0\n",
      "607 0\n",
      "608 0\n",
      "609 0\n",
      "610 0\n",
      "611 0\n",
      "612 0\n",
      "613 0\n",
      "614 0\n",
      "615 0\n",
      "616 0\n",
      "617 0\n",
      "618 1\n",
      "619 1\n",
      "620 0\n",
      "621 0\n",
      "622 1\n",
      "623 0\n",
      "624 0\n",
      "625 0\n",
      "626 2\n",
      "627 1\n",
      "628 0\n",
      "629 0\n",
      "630 1\n",
      "631 0\n",
      "632 0\n",
      "633 0\n",
      "634 0\n",
      "635 1\n",
      "636 0\n",
      "637 0\n",
      "638 0\n",
      "639 1\n",
      "640 0\n",
      "641 0\n",
      "642 1\n",
      "643 0\n",
      "644 0\n",
      "645 1\n",
      "646 1\n",
      "647 1\n",
      "648 0\n",
      "649 3\n",
      "650 1\n",
      "651 0\n",
      "652 0\n",
      "653 0\n",
      "654 0\n",
      "655 1\n",
      "656 1\n",
      "657 0\n",
      "658 0\n",
      "659 0\n",
      "660 1\n",
      "661 1\n",
      "662 0\n",
      "663 0\n",
      "664 0\n",
      "665 0\n",
      "666 1\n",
      "667 0\n",
      "668 2\n",
      "669 0\n",
      "670 2\n",
      "671 0\n",
      "672 0\n",
      "673 1\n",
      "674 1\n",
      "675 1\n",
      "676 0\n",
      "677 1\n",
      "678 2\n",
      "679 1\n",
      "680 0\n",
      "681 2\n",
      "682 0\n",
      "683 1\n",
      "684 0\n",
      "685 2\n",
      "686 0\n",
      "687 0\n",
      "688 0\n",
      "689 0\n",
      "690 1\n",
      "691 0\n",
      "692 0\n",
      "693 2\n",
      "694 0\n",
      "695 0\n",
      "696 1\n",
      "697 0\n",
      "698 0\n",
      "699 2\n",
      "700 0\n",
      "701 0\n",
      "702 0\n",
      "703 0\n",
      "704 0\n",
      "705 2\n",
      "706 2\n",
      "707 0\n",
      "708 0\n",
      "709 0\n",
      "710 0\n",
      "711 0\n",
      "712 0\n",
      "713 0\n",
      "714 0\n",
      "715 0\n",
      "716 0\n",
      "717 0\n",
      "718 0\n",
      "719 0\n",
      "720 1\n",
      "721 0\n",
      "722 1\n",
      "723 0\n",
      "724 3\n",
      "725 0\n",
      "726 0\n",
      "727 0\n",
      "728 0\n",
      "729 0\n",
      "730 2\n",
      "731 1\n",
      "732 0\n",
      "733 0\n",
      "734 1\n",
      "735 0\n",
      "736 0\n",
      "737 1\n",
      "738 1\n",
      "739 0\n",
      "740 0\n",
      "741 1\n",
      "742 1\n",
      "743 0\n",
      "744 0\n",
      "745 0\n",
      "746 0\n",
      "747 0\n",
      "748 2\n",
      "749 0\n",
      "750 2\n",
      "751 1\n",
      "752 2\n",
      "753 2\n",
      "754 0\n",
      "755 0\n",
      "756 0\n",
      "757 0\n",
      "758 0\n",
      "759 0\n",
      "760 1\n",
      "761 0\n",
      "762 1\n",
      "763 0\n",
      "764 0\n",
      "765 1\n",
      "766 0\n",
      "767 0\n",
      "768 0\n",
      "769 1\n",
      "770 0\n",
      "771 0\n",
      "772 0\n",
      "773 0\n",
      "774 0\n",
      "775 1\n",
      "776 2\n",
      "777 0\n",
      "778 0\n",
      "779 2\n",
      "780 0\n",
      "781 0\n",
      "782 0\n",
      "783 1\n",
      "784 0\n",
      "785 0\n",
      "786 0\n",
      "787 0\n",
      "788 0\n",
      "789 2\n",
      "790 0\n",
      "791 0\n",
      "792 0\n",
      "793 1\n",
      "794 1\n",
      "795 0\n",
      "796 0\n",
      "797 0\n",
      "798 0\n",
      "799 0\n",
      "800 0\n",
      "801 1\n",
      "802 2\n",
      "803 0\n",
      "804 0\n",
      "805 0\n",
      "806 0\n",
      "807 0\n",
      "808 1\n",
      "809 1\n",
      "810 1\n",
      "811 1\n",
      "812 1\n",
      "813 1\n",
      "814 1\n",
      "815 0\n",
      "816 2\n",
      "817 0\n",
      "818 1\n",
      "819 0\n",
      "820 1\n",
      "821 0\n",
      "822 0\n",
      "823 1\n",
      "824 0\n",
      "825 1\n",
      "826 3\n",
      "827 1\n",
      "828 0\n",
      "829 0\n",
      "830 0\n",
      "831 1\n",
      "832 0\n",
      "833 0\n",
      "834 0\n",
      "835 0\n",
      "836 0\n",
      "837 0\n",
      "838 0\n",
      "839 1\n",
      "840 0\n",
      "841 0\n",
      "842 0\n",
      "843 0\n",
      "844 0\n",
      "845 0\n",
      "846 0\n",
      "847 0\n",
      "848 0\n",
      "849 0\n",
      "850 1\n",
      "851 0\n",
      "852 0\n",
      "853 0\n",
      "854 1\n",
      "855 0\n",
      "856 0\n",
      "857 1\n",
      "858 0\n",
      "859 0\n",
      "860 2\n",
      "861 0\n",
      "862 2\n",
      "863 0\n",
      "864 0\n",
      "865 0\n",
      "866 0\n",
      "867 0\n",
      "868 0\n",
      "869 1\n",
      "870 0\n",
      "871 1\n",
      "872 3\n",
      "873 0\n",
      "874 0\n",
      "875 0\n",
      "876 0\n",
      "877 0\n",
      "878 0\n",
      "879 0\n",
      "880 1\n",
      "881 0\n",
      "882 0\n",
      "883 2\n",
      "884 0\n",
      "885 1\n",
      "886 1\n",
      "887 1\n",
      "888 0\n",
      "889 0\n",
      "890 1\n",
      "891 1\n",
      "892 0\n",
      "893 0\n",
      "894 2\n",
      "895 1\n",
      "896 0\n",
      "897 0\n",
      "898 3\n",
      "899 1\n",
      "900 1\n",
      "901 0\n",
      "902 0\n",
      "903 1\n",
      "904 1\n",
      "905 0\n",
      "906 2\n",
      "907 0\n",
      "908 1\n",
      "909 0\n",
      "910 3\n",
      "911 1\n",
      "912 0\n",
      "913 0\n",
      "914 1\n",
      "915 1\n",
      "916 0\n",
      "917 0\n",
      "918 2\n",
      "919 1\n",
      "920 0\n",
      "921 1\n",
      "922 0\n",
      "923 0\n",
      "924 1\n",
      "925 2\n",
      "926 0\n",
      "927 0\n",
      "928 0\n",
      "929 0\n",
      "930 0\n",
      "931 0\n",
      "932 1\n",
      "933 0\n",
      "934 1\n",
      "935 0\n",
      "936 0\n",
      "937 0\n",
      "938 0\n",
      "939 0\n",
      "940 0\n",
      "941 2\n",
      "942 0\n",
      "943 0\n",
      "944 0\n",
      "945 0\n",
      "946 0\n",
      "947 0\n",
      "948 0\n",
      "949 0\n",
      "950 0\n",
      "951 0\n",
      "952 1\n",
      "953 0\n",
      "954 1\n",
      "955 0\n",
      "956 0\n",
      "957 1\n",
      "958 0\n",
      "959 1\n",
      "960 0\n",
      "961 0\n",
      "962 2\n",
      "963 0\n",
      "964 0\n",
      "965 0\n",
      "966 1\n",
      "967 1\n",
      "968 1\n",
      "969 2\n",
      "970 0\n",
      "971 1\n",
      "972 0\n",
      "973 0\n",
      "974 0\n",
      "975 2\n",
      "976 0\n",
      "977 0\n",
      "978 0\n",
      "979 1\n",
      "980 0\n",
      "981 1\n",
      "982 0\n",
      "983 2\n",
      "984 0\n",
      "985 0\n",
      "986 0\n",
      "987 0\n",
      "988 0\n",
      "989 0\n",
      "990 1\n",
      "991 0\n",
      "992 0\n",
      "993 1\n",
      "994 0\n",
      "995 1\n",
      "996 0\n",
      "997 0\n",
      "998 2\n",
      "999 0\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/sandyalakkur/anaconda/envs/py3k/lib/python3.4/site-packages/IPython/kernel/__main__.py:20: DeprecationWarning: using a non-integer number instead of an integer will result in an error in the future\n",
      "/Users/sandyalakkur/anaconda/envs/py3k/lib/python3.4/site-packages/IPython/kernel/__main__.py:108: DeprecationWarning: converting an array with ndim > 0 to an index will result in an error in the future\n",
      "/Users/sandyalakkur/anaconda/envs/py3k/lib/python3.4/site-packages/IPython/kernel/__main__.py:160: DeprecationWarning: using a non-integer number instead of an integer will result in an error in the future\n"
     ]
    }
   ],
   "source": [
    "import numpy as np\n",
    "from random import randint\n",
    "Outbreak = Outbreaks(Size=Size,N=N,Y0=1,farms = np.arange(1,N+1),end=10,end2=50,x=x,y=y,Cows=Cows,Sheep=Sheep)\n",
    "Outbreak = np.c_[np.array([1]*Outbreak.shape[0]),Outbreak]\n",
    "Num_outbreaks = 1000\n",
    "for i in range(Num_outbreaks):\n",
    "    #Cows = np.array([randint(25,250) for p in range(N)])\n",
    "    Cows = np.array([randint(25,51) for p in range(15)]+[randint(51,76) for p in range(30)]+[randint(76,95) for p in range(20)]+[randint(95,250) for p in range(10)])\n",
    "    #Sheep = np.array([randint(25,250) for p in range(N)])\n",
    "    Sheep = np.array([randint(25,51) for p in range(15)]+[randint(51,76) for p in range(30)]+[randint(76,95) for p in range(20)]+[randint(95,250) for p in range(10)])\n",
    "    add = Outbreaks(Size=Size,N=N,Y0=1,farms = np.arange(1,N+1),end=10,end2=50,x=x,y=y,Cows=Cows,Sheep=Sheep)\n",
    "    add = np.c_[np.array([i+2]*add.shape[0]),add]\n",
    "    new_Outbreak = np.r_[Outbreak,add]\n",
    "    Outbreak = new_Outbreak \n",
    "    print(i,np.sum(add[:,12]==1))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>run</th>\n",
       "      <th>farmID</th>\n",
       "      <th>Status</th>\n",
       "      <th>DistToIndex</th>\n",
       "      <th>NumCows</th>\n",
       "      <th>NumSheep</th>\n",
       "      <th>x-coord</th>\n",
       "      <th>y-coord</th>\n",
       "      <th>DistToNearestInfected</th>\n",
       "      <th>AvgDistToInfected</th>\n",
       "      <th>CowsNearestInfected</th>\n",
       "      <th>SheepNearestInfected</th>\n",
       "      <th>Status2</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>20.757396</td>\n",
       "      <td>190</td>\n",
       "      <td>206</td>\n",
       "      <td>0.428463</td>\n",
       "      <td>1.864803</td>\n",
       "      <td>20.757396</td>\n",
       "      <td>20.757396</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>2</td>\n",
       "      <td>0</td>\n",
       "      <td>18.509917</td>\n",
       "      <td>91</td>\n",
       "      <td>87</td>\n",
       "      <td>0.391203</td>\n",
       "      <td>6.562081</td>\n",
       "      <td>18.509917</td>\n",
       "      <td>18.509917</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>1</td>\n",
       "      <td>3</td>\n",
       "      <td>0</td>\n",
       "      <td>17.417512</td>\n",
       "      <td>61</td>\n",
       "      <td>68</td>\n",
       "      <td>0.414427</td>\n",
       "      <td>10.122575</td>\n",
       "      <td>17.417512</td>\n",
       "      <td>17.417512</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>1</td>\n",
       "      <td>4</td>\n",
       "      <td>0</td>\n",
       "      <td>16.556293</td>\n",
       "      <td>43</td>\n",
       "      <td>31</td>\n",
       "      <td>0.996829</td>\n",
       "      <td>15.573354</td>\n",
       "      <td>16.556293</td>\n",
       "      <td>16.556293</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>1</td>\n",
       "      <td>5</td>\n",
       "      <td>0</td>\n",
       "      <td>16.035718</td>\n",
       "      <td>66</td>\n",
       "      <td>66</td>\n",
       "      <td>1.869069</td>\n",
       "      <td>17.528015</td>\n",
       "      <td>16.035718</td>\n",
       "      <td>16.035718</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>1</td>\n",
       "      <td>6</td>\n",
       "      <td>0</td>\n",
       "      <td>19.040339</td>\n",
       "      <td>27</td>\n",
       "      <td>36</td>\n",
       "      <td>2.842180</td>\n",
       "      <td>1.531988</td>\n",
       "      <td>19.040339</td>\n",
       "      <td>19.040339</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>1</td>\n",
       "      <td>7</td>\n",
       "      <td>0</td>\n",
       "      <td>19.803052</td>\n",
       "      <td>75</td>\n",
       "      <td>65</td>\n",
       "      <td>2.762186</td>\n",
       "      <td>0.461724</td>\n",
       "      <td>19.803052</td>\n",
       "      <td>19.803052</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>1</td>\n",
       "      <td>8</td>\n",
       "      <td>0</td>\n",
       "      <td>16.900562</td>\n",
       "      <td>65</td>\n",
       "      <td>69</td>\n",
       "      <td>2.449791</td>\n",
       "      <td>5.959629</td>\n",
       "      <td>16.900562</td>\n",
       "      <td>16.900562</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>1</td>\n",
       "      <td>9</td>\n",
       "      <td>0</td>\n",
       "      <td>16.450377</td>\n",
       "      <td>30</td>\n",
       "      <td>44</td>\n",
       "      <td>2.138491</td>\n",
       "      <td>7.732310</td>\n",
       "      <td>16.450377</td>\n",
       "      <td>16.450377</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>1</td>\n",
       "      <td>10</td>\n",
       "      <td>0</td>\n",
       "      <td>14.890135</td>\n",
       "      <td>69</td>\n",
       "      <td>64</td>\n",
       "      <td>2.561501</td>\n",
       "      <td>13.778480</td>\n",
       "      <td>14.890135</td>\n",
       "      <td>14.890135</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>1</td>\n",
       "      <td>11</td>\n",
       "      <td>0</td>\n",
       "      <td>13.484130</td>\n",
       "      <td>92</td>\n",
       "      <td>90</td>\n",
       "      <td>3.969138</td>\n",
       "      <td>13.955119</td>\n",
       "      <td>13.484130</td>\n",
       "      <td>13.484130</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>1</td>\n",
       "      <td>12</td>\n",
       "      <td>0</td>\n",
       "      <td>13.840768</td>\n",
       "      <td>67</td>\n",
       "      <td>64</td>\n",
       "      <td>3.767688</td>\n",
       "      <td>15.820503</td>\n",
       "      <td>13.840768</td>\n",
       "      <td>13.840768</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>1</td>\n",
       "      <td>13</td>\n",
       "      <td>0</td>\n",
       "      <td>16.091276</td>\n",
       "      <td>61</td>\n",
       "      <td>71</td>\n",
       "      <td>4.705106</td>\n",
       "      <td>3.921407</td>\n",
       "      <td>16.091276</td>\n",
       "      <td>16.091276</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>1</td>\n",
       "      <td>14</td>\n",
       "      <td>0</td>\n",
       "      <td>17.070804</td>\n",
       "      <td>63</td>\n",
       "      <td>76</td>\n",
       "      <td>4.246506</td>\n",
       "      <td>2.924207</td>\n",
       "      <td>17.070804</td>\n",
       "      <td>17.070804</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>1</td>\n",
       "      <td>15</td>\n",
       "      <td>0</td>\n",
       "      <td>15.546098</td>\n",
       "      <td>58</td>\n",
       "      <td>70</td>\n",
       "      <td>4.321338</td>\n",
       "      <td>5.419305</td>\n",
       "      <td>15.546098</td>\n",
       "      <td>15.546098</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>1</td>\n",
       "      <td>16</td>\n",
       "      <td>0</td>\n",
       "      <td>15.024453</td>\n",
       "      <td>70</td>\n",
       "      <td>54</td>\n",
       "      <td>5.452718</td>\n",
       "      <td>4.700431</td>\n",
       "      <td>15.024453</td>\n",
       "      <td>15.024453</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>1</td>\n",
       "      <td>17</td>\n",
       "      <td>0</td>\n",
       "      <td>16.036512</td>\n",
       "      <td>78</td>\n",
       "      <td>83</td>\n",
       "      <td>4.179607</td>\n",
       "      <td>4.741281</td>\n",
       "      <td>16.036512</td>\n",
       "      <td>16.036512</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>1</td>\n",
       "      <td>18</td>\n",
       "      <td>0</td>\n",
       "      <td>12.207746</td>\n",
       "      <td>53</td>\n",
       "      <td>74</td>\n",
       "      <td>5.860130</td>\n",
       "      <td>9.912863</td>\n",
       "      <td>12.207746</td>\n",
       "      <td>12.207746</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>1</td>\n",
       "      <td>19</td>\n",
       "      <td>0</td>\n",
       "      <td>13.352243</td>\n",
       "      <td>73</td>\n",
       "      <td>71</td>\n",
       "      <td>4.278867</td>\n",
       "      <td>11.560496</td>\n",
       "      <td>13.352243</td>\n",
       "      <td>13.352243</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>1</td>\n",
       "      <td>20</td>\n",
       "      <td>0</td>\n",
       "      <td>12.058908</td>\n",
       "      <td>34</td>\n",
       "      <td>28</td>\n",
       "      <td>5.472275</td>\n",
       "      <td>12.359426</td>\n",
       "      <td>12.058908</td>\n",
       "      <td>12.058908</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>1</td>\n",
       "      <td>21</td>\n",
       "      <td>0</td>\n",
       "      <td>11.943705</td>\n",
       "      <td>77</td>\n",
       "      <td>83</td>\n",
       "      <td>5.701048</td>\n",
       "      <td>15.881959</td>\n",
       "      <td>11.943705</td>\n",
       "      <td>11.943705</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>1</td>\n",
       "      <td>22</td>\n",
       "      <td>0</td>\n",
       "      <td>16.102311</td>\n",
       "      <td>61</td>\n",
       "      <td>56</td>\n",
       "      <td>7.859043</td>\n",
       "      <td>0.809430</td>\n",
       "      <td>16.102311</td>\n",
       "      <td>16.102311</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>1</td>\n",
       "      <td>23</td>\n",
       "      <td>0</td>\n",
       "      <td>14.762513</td>\n",
       "      <td>43</td>\n",
       "      <td>25</td>\n",
       "      <td>7.052037</td>\n",
       "      <td>3.264976</td>\n",
       "      <td>14.762513</td>\n",
       "      <td>14.762513</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>1</td>\n",
       "      <td>24</td>\n",
       "      <td>0</td>\n",
       "      <td>13.917072</td>\n",
       "      <td>91</td>\n",
       "      <td>76</td>\n",
       "      <td>6.062825</td>\n",
       "      <td>5.743816</td>\n",
       "      <td>13.917072</td>\n",
       "      <td>13.917072</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>1</td>\n",
       "      <td>25</td>\n",
       "      <td>0</td>\n",
       "      <td>11.940109</td>\n",
       "      <td>92</td>\n",
       "      <td>89</td>\n",
       "      <td>6.155299</td>\n",
       "      <td>9.874813</td>\n",
       "      <td>11.940109</td>\n",
       "      <td>11.940109</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>1</td>\n",
       "      <td>26</td>\n",
       "      <td>0</td>\n",
       "      <td>10.288324</td>\n",
       "      <td>179</td>\n",
       "      <td>210</td>\n",
       "      <td>7.325950</td>\n",
       "      <td>11.920267</td>\n",
       "      <td>10.288324</td>\n",
       "      <td>10.288324</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>1</td>\n",
       "      <td>27</td>\n",
       "      <td>0</td>\n",
       "      <td>14.432479</td>\n",
       "      <td>45</td>\n",
       "      <td>38</td>\n",
       "      <td>8.559388</td>\n",
       "      <td>2.374897</td>\n",
       "      <td>14.432479</td>\n",
       "      <td>14.432479</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>1</td>\n",
       "      <td>28</td>\n",
       "      <td>0</td>\n",
       "      <td>11.597777</td>\n",
       "      <td>42</td>\n",
       "      <td>36</td>\n",
       "      <td>9.859762</td>\n",
       "      <td>4.974922</td>\n",
       "      <td>11.597777</td>\n",
       "      <td>11.597777</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>1</td>\n",
       "      <td>29</td>\n",
       "      <td>0</td>\n",
       "      <td>11.520832</td>\n",
       "      <td>148</td>\n",
       "      <td>106</td>\n",
       "      <td>8.987554</td>\n",
       "      <td>5.926636</td>\n",
       "      <td>11.520832</td>\n",
       "      <td>11.520832</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>1</td>\n",
       "      <td>30</td>\n",
       "      <td>0</td>\n",
       "      <td>10.719772</td>\n",
       "      <td>83</td>\n",
       "      <td>92</td>\n",
       "      <td>8.402271</td>\n",
       "      <td>7.996031</td>\n",
       "      <td>10.719772</td>\n",
       "      <td>10.719772</td>\n",
       "      <td>222</td>\n",
       "      <td>200</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73976</th>\n",
       "      <td>1001</td>\n",
       "      <td>45</td>\n",
       "      <td>0</td>\n",
       "      <td>2.123537</td>\n",
       "      <td>89</td>\n",
       "      <td>86</td>\n",
       "      <td>12.698309</td>\n",
       "      <td>7.118964</td>\n",
       "      <td>2.123537</td>\n",
       "      <td>2.123537</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73977</th>\n",
       "      <td>1001</td>\n",
       "      <td>46</td>\n",
       "      <td>0</td>\n",
       "      <td>3.459675</td>\n",
       "      <td>28</td>\n",
       "      <td>30</td>\n",
       "      <td>12.213462</td>\n",
       "      <td>8.499905</td>\n",
       "      <td>3.459675</td>\n",
       "      <td>3.459675</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73978</th>\n",
       "      <td>1001</td>\n",
       "      <td>47</td>\n",
       "      <td>0</td>\n",
       "      <td>2.582251</td>\n",
       "      <td>131</td>\n",
       "      <td>209</td>\n",
       "      <td>13.633796</td>\n",
       "      <td>8.253553</td>\n",
       "      <td>2.582251</td>\n",
       "      <td>2.582251</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73979</th>\n",
       "      <td>1001</td>\n",
       "      <td>48</td>\n",
       "      <td>0</td>\n",
       "      <td>4.631741</td>\n",
       "      <td>73</td>\n",
       "      <td>66</td>\n",
       "      <td>13.963750</td>\n",
       "      <td>10.385687</td>\n",
       "      <td>4.631741</td>\n",
       "      <td>4.631741</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73980</th>\n",
       "      <td>1001</td>\n",
       "      <td>49</td>\n",
       "      <td>0</td>\n",
       "      <td>7.746719</td>\n",
       "      <td>107</td>\n",
       "      <td>106</td>\n",
       "      <td>13.311252</td>\n",
       "      <td>13.447489</td>\n",
       "      <td>7.746719</td>\n",
       "      <td>7.746719</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73981</th>\n",
       "      <td>1001</td>\n",
       "      <td>50</td>\n",
       "      <td>0</td>\n",
       "      <td>2.718216</td>\n",
       "      <td>121</td>\n",
       "      <td>214</td>\n",
       "      <td>14.277536</td>\n",
       "      <td>3.051502</td>\n",
       "      <td>2.718216</td>\n",
       "      <td>2.718216</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73982</th>\n",
       "      <td>1001</td>\n",
       "      <td>52</td>\n",
       "      <td>0</td>\n",
       "      <td>0.677720</td>\n",
       "      <td>83</td>\n",
       "      <td>85</td>\n",
       "      <td>14.668796</td>\n",
       "      <td>5.177794</td>\n",
       "      <td>0.677720</td>\n",
       "      <td>0.677720</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73983</th>\n",
       "      <td>1001</td>\n",
       "      <td>53</td>\n",
       "      <td>0</td>\n",
       "      <td>3.235225</td>\n",
       "      <td>26</td>\n",
       "      <td>44</td>\n",
       "      <td>14.581574</td>\n",
       "      <td>8.995065</td>\n",
       "      <td>3.235225</td>\n",
       "      <td>3.235225</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73984</th>\n",
       "      <td>1001</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "      <td>3.958684</td>\n",
       "      <td>75</td>\n",
       "      <td>60</td>\n",
       "      <td>15.396214</td>\n",
       "      <td>9.583559</td>\n",
       "      <td>3.958684</td>\n",
       "      <td>3.958684</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73985</th>\n",
       "      <td>1001</td>\n",
       "      <td>55</td>\n",
       "      <td>0</td>\n",
       "      <td>4.084082</td>\n",
       "      <td>54</td>\n",
       "      <td>70</td>\n",
       "      <td>14.764359</td>\n",
       "      <td>9.830775</td>\n",
       "      <td>4.084082</td>\n",
       "      <td>4.084082</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73986</th>\n",
       "      <td>1001</td>\n",
       "      <td>56</td>\n",
       "      <td>0</td>\n",
       "      <td>9.836703</td>\n",
       "      <td>64</td>\n",
       "      <td>54</td>\n",
       "      <td>15.416086</td>\n",
       "      <td>15.546453</td>\n",
       "      <td>9.836703</td>\n",
       "      <td>9.836703</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73987</th>\n",
       "      <td>1001</td>\n",
       "      <td>57</td>\n",
       "      <td>0</td>\n",
       "      <td>14.066189</td>\n",
       "      <td>195</td>\n",
       "      <td>169</td>\n",
       "      <td>14.067645</td>\n",
       "      <td>19.832654</td>\n",
       "      <td>14.066189</td>\n",
       "      <td>14.066189</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73988</th>\n",
       "      <td>1001</td>\n",
       "      <td>58</td>\n",
       "      <td>0</td>\n",
       "      <td>5.925386</td>\n",
       "      <td>76</td>\n",
       "      <td>73</td>\n",
       "      <td>17.187998</td>\n",
       "      <td>0.574325</td>\n",
       "      <td>5.925386</td>\n",
       "      <td>5.925386</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73989</th>\n",
       "      <td>1001</td>\n",
       "      <td>59</td>\n",
       "      <td>0</td>\n",
       "      <td>3.610333</td>\n",
       "      <td>80</td>\n",
       "      <td>89</td>\n",
       "      <td>16.627714</td>\n",
       "      <td>2.978049</td>\n",
       "      <td>3.610333</td>\n",
       "      <td>3.610333</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73990</th>\n",
       "      <td>1001</td>\n",
       "      <td>60</td>\n",
       "      <td>0</td>\n",
       "      <td>3.936865</td>\n",
       "      <td>83</td>\n",
       "      <td>84</td>\n",
       "      <td>16.576724</td>\n",
       "      <td>2.530989</td>\n",
       "      <td>3.936865</td>\n",
       "      <td>3.936865</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73991</th>\n",
       "      <td>1001</td>\n",
       "      <td>61</td>\n",
       "      <td>0</td>\n",
       "      <td>3.837140</td>\n",
       "      <td>87</td>\n",
       "      <td>91</td>\n",
       "      <td>17.180531</td>\n",
       "      <td>8.346103</td>\n",
       "      <td>3.837140</td>\n",
       "      <td>3.837140</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73992</th>\n",
       "      <td>1001</td>\n",
       "      <td>62</td>\n",
       "      <td>0</td>\n",
       "      <td>7.381665</td>\n",
       "      <td>69</td>\n",
       "      <td>65</td>\n",
       "      <td>17.041228</td>\n",
       "      <td>12.637765</td>\n",
       "      <td>7.381665</td>\n",
       "      <td>7.381665</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73993</th>\n",
       "      <td>1001</td>\n",
       "      <td>63</td>\n",
       "      <td>0</td>\n",
       "      <td>8.545854</td>\n",
       "      <td>87</td>\n",
       "      <td>76</td>\n",
       "      <td>17.048799</td>\n",
       "      <td>13.873426</td>\n",
       "      <td>8.545854</td>\n",
       "      <td>8.545854</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73994</th>\n",
       "      <td>1001</td>\n",
       "      <td>64</td>\n",
       "      <td>0</td>\n",
       "      <td>8.560077</td>\n",
       "      <td>245</td>\n",
       "      <td>108</td>\n",
       "      <td>17.451593</td>\n",
       "      <td>13.742620</td>\n",
       "      <td>8.560077</td>\n",
       "      <td>8.560077</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73995</th>\n",
       "      <td>1001</td>\n",
       "      <td>65</td>\n",
       "      <td>0</td>\n",
       "      <td>14.054256</td>\n",
       "      <td>90</td>\n",
       "      <td>82</td>\n",
       "      <td>16.378905</td>\n",
       "      <td>19.674271</td>\n",
       "      <td>14.054256</td>\n",
       "      <td>14.054256</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73996</th>\n",
       "      <td>1001</td>\n",
       "      <td>66</td>\n",
       "      <td>0</td>\n",
       "      <td>4.886189</td>\n",
       "      <td>38</td>\n",
       "      <td>51</td>\n",
       "      <td>18.330410</td>\n",
       "      <td>2.952671</td>\n",
       "      <td>4.886189</td>\n",
       "      <td>4.886189</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73997</th>\n",
       "      <td>1001</td>\n",
       "      <td>67</td>\n",
       "      <td>0</td>\n",
       "      <td>5.822442</td>\n",
       "      <td>80</td>\n",
       "      <td>90</td>\n",
       "      <td>18.864681</td>\n",
       "      <td>2.107629</td>\n",
       "      <td>5.822442</td>\n",
       "      <td>5.822442</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73998</th>\n",
       "      <td>1001</td>\n",
       "      <td>68</td>\n",
       "      <td>0</td>\n",
       "      <td>4.509184</td>\n",
       "      <td>63</td>\n",
       "      <td>54</td>\n",
       "      <td>18.775315</td>\n",
       "      <td>6.568489</td>\n",
       "      <td>4.509184</td>\n",
       "      <td>4.509184</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73999</th>\n",
       "      <td>1001</td>\n",
       "      <td>69</td>\n",
       "      <td>0</td>\n",
       "      <td>4.066472</td>\n",
       "      <td>54</td>\n",
       "      <td>73</td>\n",
       "      <td>18.186413</td>\n",
       "      <td>7.081517</td>\n",
       "      <td>4.066472</td>\n",
       "      <td>4.066472</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>74000</th>\n",
       "      <td>1001</td>\n",
       "      <td>70</td>\n",
       "      <td>0</td>\n",
       "      <td>6.158406</td>\n",
       "      <td>69</td>\n",
       "      <td>56</td>\n",
       "      <td>18.492750</td>\n",
       "      <td>10.314425</td>\n",
       "      <td>6.158406</td>\n",
       "      <td>6.158406</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>74001</th>\n",
       "      <td>1001</td>\n",
       "      <td>71</td>\n",
       "      <td>0</td>\n",
       "      <td>6.034976</td>\n",
       "      <td>52</td>\n",
       "      <td>74</td>\n",
       "      <td>18.078301</td>\n",
       "      <td>10.504857</td>\n",
       "      <td>6.034976</td>\n",
       "      <td>6.034976</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>74002</th>\n",
       "      <td>1001</td>\n",
       "      <td>72</td>\n",
       "      <td>0</td>\n",
       "      <td>9.298684</td>\n",
       "      <td>72</td>\n",
       "      <td>69</td>\n",
       "      <td>19.985391</td>\n",
       "      <td>13.156041</td>\n",
       "      <td>9.298684</td>\n",
       "      <td>9.298684</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>74003</th>\n",
       "      <td>1001</td>\n",
       "      <td>73</td>\n",
       "      <td>0</td>\n",
       "      <td>10.615504</td>\n",
       "      <td>93</td>\n",
       "      <td>81</td>\n",
       "      <td>18.989877</td>\n",
       "      <td>15.310801</td>\n",
       "      <td>10.615504</td>\n",
       "      <td>10.615504</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>74004</th>\n",
       "      <td>1001</td>\n",
       "      <td>74</td>\n",
       "      <td>0</td>\n",
       "      <td>9.998750</td>\n",
       "      <td>83</td>\n",
       "      <td>93</td>\n",
       "      <td>19.700055</td>\n",
       "      <td>14.208172</td>\n",
       "      <td>9.998750</td>\n",
       "      <td>9.998750</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>74005</th>\n",
       "      <td>1001</td>\n",
       "      <td>75</td>\n",
       "      <td>0</td>\n",
       "      <td>12.295354</td>\n",
       "      <td>76</td>\n",
       "      <td>92</td>\n",
       "      <td>18.040890</td>\n",
       "      <td>17.493436</td>\n",
       "      <td>12.295354</td>\n",
       "      <td>12.295354</td>\n",
       "      <td>67</td>\n",
       "      <td>54</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>74006 rows × 13 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "        run  farmID  Status  DistToIndex  NumCows  NumSheep    x-coord  \\\n",
       "0         1       1       0    20.757396      190       206   0.428463   \n",
       "1         1       2       0    18.509917       91        87   0.391203   \n",
       "2         1       3       0    17.417512       61        68   0.414427   \n",
       "3         1       4       0    16.556293       43        31   0.996829   \n",
       "4         1       5       0    16.035718       66        66   1.869069   \n",
       "5         1       6       0    19.040339       27        36   2.842180   \n",
       "6         1       7       0    19.803052       75        65   2.762186   \n",
       "7         1       8       0    16.900562       65        69   2.449791   \n",
       "8         1       9       0    16.450377       30        44   2.138491   \n",
       "9         1      10       0    14.890135       69        64   2.561501   \n",
       "10        1      11       0    13.484130       92        90   3.969138   \n",
       "11        1      12       0    13.840768       67        64   3.767688   \n",
       "12        1      13       0    16.091276       61        71   4.705106   \n",
       "13        1      14       0    17.070804       63        76   4.246506   \n",
       "14        1      15       0    15.546098       58        70   4.321338   \n",
       "15        1      16       0    15.024453       70        54   5.452718   \n",
       "16        1      17       0    16.036512       78        83   4.179607   \n",
       "17        1      18       0    12.207746       53        74   5.860130   \n",
       "18        1      19       0    13.352243       73        71   4.278867   \n",
       "19        1      20       0    12.058908       34        28   5.472275   \n",
       "20        1      21       0    11.943705       77        83   5.701048   \n",
       "21        1      22       0    16.102311       61        56   7.859043   \n",
       "22        1      23       0    14.762513       43        25   7.052037   \n",
       "23        1      24       0    13.917072       91        76   6.062825   \n",
       "24        1      25       0    11.940109       92        89   6.155299   \n",
       "25        1      26       0    10.288324      179       210   7.325950   \n",
       "26        1      27       0    14.432479       45        38   8.559388   \n",
       "27        1      28       0    11.597777       42        36   9.859762   \n",
       "28        1      29       0    11.520832      148       106   8.987554   \n",
       "29        1      30       0    10.719772       83        92   8.402271   \n",
       "...     ...     ...     ...          ...      ...       ...        ...   \n",
       "73976  1001      45       0     2.123537       89        86  12.698309   \n",
       "73977  1001      46       0     3.459675       28        30  12.213462   \n",
       "73978  1001      47       0     2.582251      131       209  13.633796   \n",
       "73979  1001      48       0     4.631741       73        66  13.963750   \n",
       "73980  1001      49       0     7.746719      107       106  13.311252   \n",
       "73981  1001      50       0     2.718216      121       214  14.277536   \n",
       "73982  1001      52       0     0.677720       83        85  14.668796   \n",
       "73983  1001      53       0     3.235225       26        44  14.581574   \n",
       "73984  1001      54       0     3.958684       75        60  15.396214   \n",
       "73985  1001      55       0     4.084082       54        70  14.764359   \n",
       "73986  1001      56       0     9.836703       64        54  15.416086   \n",
       "73987  1001      57       0    14.066189      195       169  14.067645   \n",
       "73988  1001      58       0     5.925386       76        73  17.187998   \n",
       "73989  1001      59       0     3.610333       80        89  16.627714   \n",
       "73990  1001      60       0     3.936865       83        84  16.576724   \n",
       "73991  1001      61       0     3.837140       87        91  17.180531   \n",
       "73992  1001      62       0     7.381665       69        65  17.041228   \n",
       "73993  1001      63       0     8.545854       87        76  17.048799   \n",
       "73994  1001      64       0     8.560077      245       108  17.451593   \n",
       "73995  1001      65       0    14.054256       90        82  16.378905   \n",
       "73996  1001      66       0     4.886189       38        51  18.330410   \n",
       "73997  1001      67       0     5.822442       80        90  18.864681   \n",
       "73998  1001      68       0     4.509184       63        54  18.775315   \n",
       "73999  1001      69       0     4.066472       54        73  18.186413   \n",
       "74000  1001      70       0     6.158406       69        56  18.492750   \n",
       "74001  1001      71       0     6.034976       52        74  18.078301   \n",
       "74002  1001      72       0     9.298684       72        69  19.985391   \n",
       "74003  1001      73       0    10.615504       93        81  18.989877   \n",
       "74004  1001      74       0     9.998750       83        93  19.700055   \n",
       "74005  1001      75       0    12.295354       76        92  18.040890   \n",
       "\n",
       "         y-coord  DistToNearestInfected  AvgDistToInfected  \\\n",
       "0       1.864803              20.757396          20.757396   \n",
       "1       6.562081              18.509917          18.509917   \n",
       "2      10.122575              17.417512          17.417512   \n",
       "3      15.573354              16.556293          16.556293   \n",
       "4      17.528015              16.035718          16.035718   \n",
       "5       1.531988              19.040339          19.040339   \n",
       "6       0.461724              19.803052          19.803052   \n",
       "7       5.959629              16.900562          16.900562   \n",
       "8       7.732310              16.450377          16.450377   \n",
       "9      13.778480              14.890135          14.890135   \n",
       "10     13.955119              13.484130          13.484130   \n",
       "11     15.820503              13.840768          13.840768   \n",
       "12      3.921407              16.091276          16.091276   \n",
       "13      2.924207              17.070804          17.070804   \n",
       "14      5.419305              15.546098          15.546098   \n",
       "15      4.700431              15.024453          15.024453   \n",
       "16      4.741281              16.036512          16.036512   \n",
       "17      9.912863              12.207746          12.207746   \n",
       "18     11.560496              13.352243          13.352243   \n",
       "19     12.359426              12.058908          12.058908   \n",
       "20     15.881959              11.943705          11.943705   \n",
       "21      0.809430              16.102311          16.102311   \n",
       "22      3.264976              14.762513          14.762513   \n",
       "23      5.743816              13.917072          13.917072   \n",
       "24      9.874813              11.940109          11.940109   \n",
       "25     11.920267              10.288324          10.288324   \n",
       "26      2.374897              14.432479          14.432479   \n",
       "27      4.974922              11.597777          11.597777   \n",
       "28      5.926636              11.520832          11.520832   \n",
       "29      7.996031              10.719772          10.719772   \n",
       "...          ...                    ...                ...   \n",
       "73976   7.118964               2.123537           2.123537   \n",
       "73977   8.499905               3.459675           3.459675   \n",
       "73978   8.253553               2.582251           2.582251   \n",
       "73979  10.385687               4.631741           4.631741   \n",
       "73980  13.447489               7.746719           7.746719   \n",
       "73981   3.051502               2.718216           2.718216   \n",
       "73982   5.177794               0.677720           0.677720   \n",
       "73983   8.995065               3.235225           3.235225   \n",
       "73984   9.583559               3.958684           3.958684   \n",
       "73985   9.830775               4.084082           4.084082   \n",
       "73986  15.546453               9.836703           9.836703   \n",
       "73987  19.832654              14.066189          14.066189   \n",
       "73988   0.574325               5.925386           5.925386   \n",
       "73989   2.978049               3.610333           3.610333   \n",
       "73990   2.530989               3.936865           3.936865   \n",
       "73991   8.346103               3.837140           3.837140   \n",
       "73992  12.637765               7.381665           7.381665   \n",
       "73993  13.873426               8.545854           8.545854   \n",
       "73994  13.742620               8.560077           8.560077   \n",
       "73995  19.674271              14.054256          14.054256   \n",
       "73996   2.952671               4.886189           4.886189   \n",
       "73997   2.107629               5.822442           5.822442   \n",
       "73998   6.568489               4.509184           4.509184   \n",
       "73999   7.081517               4.066472           4.066472   \n",
       "74000  10.314425               6.158406           6.158406   \n",
       "74001  10.504857               6.034976           6.034976   \n",
       "74002  13.156041               9.298684           9.298684   \n",
       "74003  15.310801              10.615504          10.615504   \n",
       "74004  14.208172               9.998750           9.998750   \n",
       "74005  17.493436              12.295354          12.295354   \n",
       "\n",
       "       CowsNearestInfected  SheepNearestInfected  Status2  \n",
       "0                      222                   200        0  \n",
       "1                      222                   200        0  \n",
       "2                      222                   200        0  \n",
       "3                      222                   200        0  \n",
       "4                      222                   200        0  \n",
       "5                      222                   200        0  \n",
       "6                      222                   200        0  \n",
       "7                      222                   200        0  \n",
       "8                      222                   200        0  \n",
       "9                      222                   200        0  \n",
       "10                     222                   200        0  \n",
       "11                     222                   200        0  \n",
       "12                     222                   200        0  \n",
       "13                     222                   200        0  \n",
       "14                     222                   200        0  \n",
       "15                     222                   200        0  \n",
       "16                     222                   200        0  \n",
       "17                     222                   200        0  \n",
       "18                     222                   200        0  \n",
       "19                     222                   200        0  \n",
       "20                     222                   200        0  \n",
       "21                     222                   200        0  \n",
       "22                     222                   200        0  \n",
       "23                     222                   200        0  \n",
       "24                     222                   200        0  \n",
       "25                     222                   200        0  \n",
       "26                     222                   200        0  \n",
       "27                     222                   200        0  \n",
       "28                     222                   200        0  \n",
       "29                     222                   200        0  \n",
       "...                    ...                   ...      ...  \n",
       "73976                   67                    54        0  \n",
       "73977                   67                    54        0  \n",
       "73978                   67                    54        0  \n",
       "73979                   67                    54        0  \n",
       "73980                   67                    54        0  \n",
       "73981                   67                    54        0  \n",
       "73982                   67                    54        0  \n",
       "73983                   67                    54        0  \n",
       "73984                   67                    54        0  \n",
       "73985                   67                    54        0  \n",
       "73986                   67                    54        0  \n",
       "73987                   67                    54        0  \n",
       "73988                   67                    54        0  \n",
       "73989                   67                    54        0  \n",
       "73990                   67                    54        0  \n",
       "73991                   67                    54        0  \n",
       "73992                   67                    54        0  \n",
       "73993                   67                    54        0  \n",
       "73994                   67                    54        0  \n",
       "73995                   67                    54        0  \n",
       "73996                   67                    54        0  \n",
       "73997                   67                    54        0  \n",
       "73998                   67                    54        0  \n",
       "73999                   67                    54        0  \n",
       "74000                   67                    54        0  \n",
       "74001                   67                    54        0  \n",
       "74002                   67                    54        0  \n",
       "74003                   67                    54        0  \n",
       "74004                   67                    54        0  \n",
       "74005                   67                    54        0  \n",
       "\n",
       "[74006 rows x 13 columns]"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import pandas as pd\n",
    "df = pd.DataFrame(Outbreak)\n",
    "df.columns = ['run','farmID','Status','DistToIndex','NumCows','NumSheep','x-coord','y-coord','DistToNearestInfected','AvgDistToInfected','CowsNearestInfected','SheepNearestInfected','Status2']\n",
    "long_random_everything_andind_anddist2 = df\n",
    "long_random_everything_andind_anddist2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "491"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.sum(long_random_everything_andind_anddist2['Status2']==1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "dist_diff = (long_random_everything_andind_anddist2['DistToNearestInfected'] - long_random_everything_andind_anddist2['DistToNearestInfected'].mean())/5\n",
    "closest_cow_diff = (long_random_everything_andind_anddist2['CowsNearestInfected'] - long_random_everything_andind_anddist2['CowsNearestInfected'].mean())/20\n",
    "closest_sheep_diff = (long_random_everything_andind_anddist2['SheepNearestInfected'] - long_random_everything_andind_anddist2['SheepNearestInfected'].mean())/20\n",
    "status = long_random_everything_andind_anddist2['Status2']\n",
    "from pymc import Normal, Binomial, Gamma, Lambda, invlogit, MCMC, Matplot, Bernoulli, MAP, AdaptiveMetropolis\n",
    "#N = df.shape[0]\n",
    "\n",
    "def pooled_model():\n",
    "    \n",
    "    # Common slope & intercept prior\n",
    "    intercept = Normal('intercept', mu=0., tau=0.001, value = 0) \n",
    "    \n",
    "    first_coef = Normal('first_coef', mu=0., tau=0.001, value = 0)\n",
    "    \n",
    "    size_coef = Normal ('size_coef', mu=0., tau= 0.001, value=[0]*2)\n",
    "    \n",
    "    #likelihood model\n",
    "    prob = Lambda('prob', lambda intercept=intercept,first_coef=first_coef, size_coef= size_coef: \n",
    "                  invlogit(intercept + first_coef*dist_diff + size_coef[0]*closest_cow_diff + \n",
    "                           size_coef[1]*closest_sheep_diff))\n",
    "    \n",
    "    y = Bernoulli('y', p=prob, value=status, observed=True)\n",
    "    \n",
    "    return locals()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      " [-----------------100%-----------------] 10000 of 10000 complete in 128.1 sec"
     ]
    }
   ],
   "source": [
    "chains = 2\n",
    "iterations = 10000\n",
    "burn = 4000\n",
    "M_pooled = MCMC(pooled_model())\n",
    "M_map = MAP(pooled_model())\n",
    "M_pooled = MCMC(M_map)\n",
    "#M_pooled.use_step_method(AdaptiveMetropolis, M_pooled.b)\n",
    "for i in range(chains):\n",
    "    M_pooled.sample(iterations, burn)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ4AAAEgCAYAAACEvBUKAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAFC1JREFUeJzt3Xu4JVV95vHvi7eoNCLioGgQJTFRbGMrISQYbYeB6Tw6\nXqOOPBPFSWImcWA0OpoEEzCAMfNMCPE6GjHiLQrBoCQBFKVJlAkYQAFj4oWLyKVBGhQQQeGXP2od\n2Jycffp09znrnL3P9/M8+6F616qqVWuvXW+tqjqbVBWSJPWyw3JXQJK0uhg8kqSuDB5JUlcGjySp\nK4NHktSVwSNJ6srg0YqR5ANJjmrTv5jkX0bmXZ7kgDHLrU9yZa96LpckG5P86nLXQ8tjmj5/g2fC\nJHlCks8luSnJ15M8f2TenknuSnLzyOvwkfkHJ7k6yWVJ1o+8v1eSLyTJFrb9yCTHt3V8L8lXkxyZ\n5EGLtHvVXlTVP1TVT881bzGNht0Cyh6Z5EOLXYetsCRtoOXTTqi+376r1yb5UJKdxhTf5s9/a/p5\nDwbPBElyX+CTwKeAhwKvAj6c5CdnFd2pqta01zEjy/4RsA74n8DbR8q/DXhNzfPXxEl2Af4/8ABg\nv6raCTgQeAiw15i6bot5w28JdDuYJ7lPj+1oohTwnKpaA/wMsBZ40/JWaekZPJPlp4FHVtVxNTgL\n+ALwK7PKzfW5Pgy4qqo2AZ8FHgeQ5JeBK6vqi1vY9m8D362q/1ZV3wKoqm9X1Wur6uK2rruS/FaS\nrwP/2t57TpIvJbmxjarWzqwwybokF7TR08eAHxuZN9fls32TfCXJ5iTvT/KAuSqaZPckJye5Lsml\nSQ7dwr6lLTczYnx5kiuSXJ/k99q8DcDvAi9tZ6cXtvcfMjIK/HaSo5Ls0OYd0vb52CTfAY5q7bD3\nSF0f3s54d03y0CR/0+q9OcmpSR41Zh9/IsnZbeR7fWs/TbD23fw0sPc8xfZM8vn2nTkjycNmZiQ5\nKck1rU+cneSJ7f1XAQcDb2h995NLuiMLYPBMvh2AJ81674okV7aD80zHvB54WDuQHQhckmRH4HCG\nA+qW/CfgEwso9zzgZ4EnJlkHHA/8OrAL8B7gU0nul+T+wCnACQyjt5OAFzF+9BGGL89BDCOsxzPH\nmWE76J8KXAjsDhwAvCbJQQuo+4z92/oPAP4gyU9V1enAW4CPtZHkulb2A8AdrU7rWv1+bWRd+wLf\nBP4D8IcMbfiykfkvATZW1XfaPh4P7NFetwHvGFPHo4DTq2pn4FEMo1ZNppkTn0cDG4Bz5yl3MHAI\nQ3+6P/D6kfl/C/wE8HDgAuAjAFX13jb9x63vPm/xd2HrGDyT5V+B65L873bwPgh4BvDANv96YB+G\ng9bTgDXc0/nuAn4T+CuG0cuvMxwI3wY8JcN9o9NHz8Zn2QW4ZgF1/KOquqmqbme4FPieqvpiG6F9\nELgd+HlgP+C+VfVnVXVnVZ0MzDfqKuAdVXVVVd0IHMO9D+AzfhbYtaqOrqofVdVlwPuA/7qAus94\nc1XdXlUXAV9muAQCwxf/7kuBSXYDfgl4bVXdVlXXA8fN2tbVVfXOqrqrqn4AfHTW/IPbe1TV5qr6\n66r6QVXdwhB0zxxTxzsYzn4fVVV3VNU5W7F/WjkCnJLke8C3GE5Sjh5TtoD3V9U3Wl86EXjK3TOr\nPlBVt1bVD4E3Az+TZM2sba0IBs8EaR3q+cCzGULgtQyd79tt/q1VdUE7yF3HcC/noCQPbvM/V1U/\nX1XPYuiET2UYcZwAvILhLPp9YzZ/A8MIYktGL489Bnhdu7x0Y5IbgUcDj2zrumrWsldsxbq/NaY+\njwF2n7XN32U4Q1yoa0emvw/sOKbcY4D7AdeMbOv/MZxxzlVngI3Ag5Lsm2RPhlD7a4AkD0rynnbD\n+bvA2cBDkjkf+ngDw2d4XpJLkrxyK/ZPK0cBz2v3TNcD/xHYJ8m7c88DQr8zUn60b95G65tJ7pPk\nrUm+0frOZa3Mrku/C1tvW28Aa5m0+ynrZ/6d5BzgL7aw2L1OMNqB7O0MwfRw4D5VdWWS64Anj1nH\nmcALkrx5vocQuPelsm8Bx1TVW2YXSvJMhktEox4DfGOede8xa/rqOcpcCVxWVY+fZz3bavZ+X8kw\ngntYG1FucZmqujPJiQyjteuAU6vq1jb7dQyX+PatquuSPIXhkknmWM8mhhElSfYHzkxydlVdus17\np2VVVX+f5O0Ml8SexXCFYqEOBp4LHFBVVyTZGdjMPaOcFfU0pCOeCZNkbZIfa2fHrwd2Y7jPQDuL\n/qkkO7R7O28Dzqqqm2et5teA89ulpBuAByZ5AvAshqH+XI4FdgJOSLJH296jkvxJktn3mGb8OfA/\nWr2S5MFJnt3uLZ0D/CjJYe2y4QsZLpON3XXg1W2buzDcm5rrhvp5wM1J3pDkge1M8ElJ9plnvQt1\nLcPlrQBU1TUMN4OPTbKmtfteSZ6xhfXMXG67+zJbsyPDWex32z4eMW4FSV7c7gkA3MRwYBkXfpoc\nxzE8RPNzY+aP6687MpwEbW5XOGaf7G2iPVC0Ehg8k+dXGM70NzEExYHtEhwMHes04HvAxQwHsXvd\nB0myK3AY8PsAVfUjhpHP54B3AXM+Adbuq/wC8EPg3HZN+kyGg97MKGX2Wfn5DPeS3sFw9vV14OVt\n3g+BFzLcKL2B4Sb7ybM3O2v6IwwH+m+2dR09u2xV3Qk8h+Ha96UM973eyxCac+7aHNsZ56T23xuS\n/FObfjnDTd5/bvt4EvCIMeum1fE84BaGS46njcw6juF+3XcYgvm0eeqzD/CPSW5meMT+sKq6fJ66\nawK0h0xOAN44rsis6Zl/f5DhUvVVwCUMf/owWvZ4hgd+bkyykIeEllT8H8FJknpyxCNJ6srgkSR1\nZfBIkroyeCRJXU313/Ek8ckJaYJV1bL9tb3Hj+037vOb6hFPVc35OuKII8bOW40v28P2WKFtsaw/\n8bKC22VS6jD285vq4JEkrTwGjySpq1UZPOvXr1/uKqwotse92R73sC3mthLaZZLrMO2/XDDVOydN\nueX+GX+PH9vHezySpJXB4JEkdWXwSJK6MngkSV0ZPJKkrgweSVJXBo8kqSuDR5LUlcEjSerK4JEk\ndWXwSJK6MngkSV0ZPJKkrgweSVJXBo8kqSuDR5LUlcEjSerK4JEkdWXwSJK6MngkSV0ZPJKkrgwe\nSVJXBo8kqSuDR5LUlcEjSerK4JEkdWXwSJK6MngkSV0ZPJKkrgweSVJXBo8kqSuDR5LUlcEjSerK\n4JEkdWXwSJK6MngkSV0ZPJKkrgweSVJXBo8kqSuDR5LUlcEjSerK4JEkdWXwSJK6MngkSV0ZPJKk\nrgweSVJXBo8kqSuDR5LUlcEjSerK4JEkdWXwSJK6MngkSV0ZPJKkrgweSVJXBo8kqSuDR5LUlcEj\nSerK4JEkdWXwSJK6MngkSV0ZPJKkrgweSVJXBo8kqSuDR5LUlcEjSerK4JEkdWXwSJK6MngkSV0Z\nPJKkrgweSVJXBo8kqSuDR5LUlcEjSerK4JEkdWXwSJK6MngkSV0ZPJKkrgweSVJXBo8kqSuDR5LU\nlcEjSerK4JEkdWXwSJK6MngkSV0ZPJKkrgweSVJXBo8kqSuDR5LUlcEjSerK4JEkdWXwSJK6Mngk\nSV0ZPJKkrgweSVJXBo8kqSuDR5LUlcEjSerK4JEkdWXwSJK6MngkSV0ZPJKkru673BXQZKgqLr74\nYgDWrl1LkmWukaRJ5YhHW1RVvOQlh7P//pew//6X8NKXvomqWu5qSZpQmfIDyFTsnIOLpTHdXX8q\nLHfPt4dsn7Gfn5faGg/uq8+kfOYGpKaNwdP45R5v5lLbGWc8CYANG77Cxz9+tPd5JG0TL7VpQXy4\nQMtguTuZx4/tM/bzM3gkrVQGz2Qb+/n5VJskqSuDR5LUlcEjSerK4JEkdTVv8CT5wpZWkOQ1SR64\neFVauCQPSfKby7FtSdK22e6n2pJcBuxTVTdsxTI7VNVd27XhYT17AqdW1doxRXwqRSuOj6Yv2HI3\njMePMRbYh7ftqbYkt7T/rk+yMclJSb6a5MPt/cOA3YGzkny2vXdQknOSnJ/kxCQPbu9fnuStSc4H\nXpxkQyvzpSRntjIPTvL+JOcmuSDJc9v7hyT5ZJKzknwtyR+0Kr4V2CvJhUn+eKGNJi0Xf/dOk24x\n+vC8I54kN1fVmiTrgVOAJwLXAF8AXl9V57QRz9OqanOSXYGTgQ1VdVuSNwL3r6qjWrl3VtX/TfJw\n4HzgF6vqiiQ7V9VNSd4CfKWqPpJkZ+BcYB3wEuAtwN7AbcAXgUOA7wB/44hH22K1DDQmONeW+xOa\n3JZbQhdddBH7738Jt9xyMABr1nyUz3/+STz5yU+eXXRRfqvtvKq6GiDJl4A9gXNmldmPIZzOaUOv\n+88q8/GRcmdX1RUAVXVTe/8g4L8keX379wOAPRg6wKer6sa2/U8AT2cIQy2R1XJgnnYr7XOc4CDU\nItma4Ll9ZPrOeZb9TFUdPGbere2/xfg0fGFVfX30jSQ/N6tMgO2+R6T5eYBYfP7unSbd2rVr2bDh\nY5xxxkeBoQ+vXfuyrVrHYvxI6M3ATsBmhktj70yyV1V9s93f2X12kLRy70qyZ1VdnmSXqtoMnAEc\nBhwKkGRdVV3IEDQHJnko8APgecArgVuANYuwD1IXSTjxxGNGbsy+zNDRRFmMPryl4Kkx06PeC5ye\n5KqqOiDJIcBfJnlAm384cK/gqarrk7wK+ESSHYBNwH8GjgKOS3IRw4MPlwLPbds+j+H+0aOBD1XV\nBTA88p3kYuDvquqNC9lpaTklmet6uDQxtrcPT8SPhLYwe1pVHbqVi678nZM0znIPBT1+bJ+J/5HQ\nwk4gSVNhIkY822Gqd06aco54JtvEj3gkSVPC4JEkdWXwSJK6MngkSV0ZPJKkrgweSVJXBo8kqSuD\nR5LUlcEjSerK4JEkdWXwSJK6MngkSV0ZPJKkrgweSVJXBo8kqSuDR5LUlcEjSerK4JEkdWXwSJK6\nMngkSV0ZPJKkrgweSVJXBo8kqSuDR5LUlcEjSerK4JEkdWXwSJK6MngkSV0ZPJKkrgweSVJXBo8k\nqSuDR5LUlcEjSerK4JEkdWXwSJK6MngkSV0ZPJKkrgweSVJXBo8kqSuDR5LUlcEjSerK4JEkdWXw\nSJK6MngkSV0ZPJKkrgweSVJXBo8kqSuDR5LUlcEjSerK4JEkdWXwSJK6MngkSV0ZPJKkrgweSVJX\nBo8kqSuDR5LUlcEjSerK4JEkdWXwSJK6MngkSV0ZPJKkrgweSVJXBo8kqSuDR5LUlcEjSerK4JEk\ndWXwSJK6MngkSV0ZPJKkrgweSVJXBo8kqSuDR5LUlcEjSerK4JEkdWXwSJK6MngkSV0ZPJKkrgwe\nSVJXBo8kqSuDR5LUlcEjSerK4JEkdWXwSJK6MngkSV0ZPJKkrgweSVJXBo8kqSuDR5LUlcEjSerK\n4JEkdWXwSJK6MngkSV0ZPJKkrgweSVJXBo8kqSuDR5LUlcEjSerK4JEkdWXwSJK6MngkSV0ZPJKk\nrgweSVJXBo8kqSuDR5LUlcEjSerK4JEkdWXwSJK6MngkSV0ZPJKkrgweSVJXBo8kqSuDR5LU1aoM\nno0bNy53FVYU2+PebI972BZzWwntMsl1MHhke8xie9zDtpjbSmiXSa7DqgweSdLyMXgkSV2lqpa7\nDksmyfTunLQKVFWWa9seP7bfuM9vqoNHkrTyeKlNktSVwSNJ6mpqgyfJi5N8JcmdSZ468v6eSW5L\ncmF7vWvM8kcm+fZIuQ39ar/4FqE9dknymSRfS/LpJDv3q/3iGtcWI/P3SHJLkteNWX5V9I2R+Vtq\nj6npG/NJ8v4km5JcvFK2mWR9ku+O9MU3LVE9fjzJWa2fXJLksO1Z39QGD3Ax8ALg7+eY942qWtde\nvzVm+QKOHSl3+pLVtI/tbY/fAT5TVY8HPtv+PanmawuAY4G/nWf51dQ3YMvtMU19Yz5/AfQ+yVjI\nNs8e6YtHL1E9fgi8tqr2BvYDXp3kCaMFkly+0JXdd3HrtnJU1b8AJNv1UMyyPVGz2BahPZ4LPLNN\nnwBsZEIPMPO1RZLnA5cCt25hNauibyywPaamb8ynqv4hyZ4rcJtL3her6lrg2jZ9S5KvArsDXx0t\nttD1TfOIZz6PbcPSjUmePk+5Q5N8Ocnx03r5oFlIe+xWVZva9CZgt0516ybJjsAbgCMXUHzq+8ZW\ntMfU940VrIBfaH3x75I8cak32IJwHXDutq5jokc8ST4DPGKOWb9XVaeOWexq4Mer6sZ2PfuUJHtX\n1c2zyr0b+MM2fRTwJ8CvLka9l8oSt8fdqqpW+t84bGNbHAn8aVV9P/MPDVdL3ziShbXH3Sahb0yZ\nCxi+v99P8kvAKcDjl2pj7WTkr4D/1UY+hwO/3GbvnuTCNv35qjp03HomOniq6sBtWOYO4I42fUGS\nbwI/yfABjpa7bmY6yfuAcV/OFWMp2wPYlOQRVXVtkkcC17GCbUtbAPsCL0ryf4CdgbuS3FZV93rg\nYrX0DRbYHkxY35gmoyeIVXVakncl2aWqNi/2tpLcDzgZ+HBVndK2eQxwTJt/WVWtW8i6VsultrvP\n1pLsmuQ+bfpxDAfZS//dAsMXaMYLGG7AToutbg/gU8Ar2vQrGM6spsHdbVFVz6iqx1bVY4HjgGPm\nOMiumr6x0PZgevvGipdkt5nRaJJ9GX4UYClCJ8DxwD9X1XHbvcKqmsoXwwHhSuA2hptip7X3XwRc\nAlwInA88e2SZPwee2qY/CFwEfJnhi7Tbcu/TMrXH09r0LsCZwNeATwM7L/c+LXZbzCpzBPDbq7lv\nLKA9pq5vbKGd/pLh0vTtrb1e2XGbd7Rt/nfgN4DfaPNf3b6/XwLOAfZbono8HbirbefC9towq8yl\nC12fP5kjSepqtVxqkyStEAaPJKkrg0eS1JXBI0nqyuCRJHVl8EiSujJ4JEldGTySpK7+DbyBusgz\nGj5fAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x1096f1be0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "import matplotlib.pyplot as plt\n",
    "Matplot.summary_plot(M_pooled.intercept)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZkAAAEgCAYAAABmYA5zAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAFGxJREFUeJzt3Xe4JGWZhvH7YTAhI4q6KqiMopgYFEUWHcMgymLEVVaU\nVVfXtK5KWHNY0VUMexkQMIOCLohixACsLII5EYdgQEcZMpIEAxLe/aPqSNOeNOHrPmfm/l1XX9On\nq07VW9Vf11PfV3V6UlVIktTCeuMuQJK09jJkJEnNGDKSpGYMGUlSM4aMJKkZQ0aS1IwhozkjycFJ\n3tY/f2SSnw1M+02SHab4vaVJVoyqznFJcnySF4y7Do3G2vJ+GzLzTJL7JTkuyRVJfpnkqQPTFiW5\nIclVA483DkzfLcn5SZYnWTrw+uZJvpckM6z7LkkO6pfx+yRnJXlLkg3W0OZV/6CqvlNV951s2po0\nGGyzmPctST69pmtYCU32gUanP1n6Y//ZvDDJp5PcZorZV/n9Xpl23ZohM48kWR/4CnAkcDvgxcD/\nJLn30Ky3qaqF/WOfgd99J7A18HJg/4H59wP2rGn+MjfJxsAPgFsA21XVbYDHARsBm09R66qYNuga\nGNmBO8mCUaxHc1oBT6qqhcADgcXAm8ZbUluGzPxyX+AuVbVvdb4FfA94ztB8k72vtwfOq6qLgP8D\n7gmQZBdgRVX9ZIZ1/wdwZVU9u6rOAaiqc6tqr6pa1i/rhiT/nuSXwM/7156U5JQkl/e9pcUTC0yy\ndZKT+l7R4cAtB6ZNNgS2bZIzklyW5BNJbjFZoUk2SfKFJBcn+XWSV8ywbel/b6In+Nwkv01ySZI3\n9NN2Al4P7NqfhZ7cv77RQO/u3CRvS7JeP+15/Ta/L8nvgLf1++EBA7XesT+zvUOS2yX5Wl/3ZUm+\nmmTTKbbxXklO6Hu0l/T7T/NI/1n8X+AB08y2KMl3+8/IMUluPzEhyRFJLujbwAlJ7t+//mJgN+A1\nfVv9StMNmYEhM/+tB2w59Npvk6zoD8QTjfIS4Pb9QetxwOlJNgTeSHfwnMljgS/OYr6dgYcC90+y\nNXAQ8CJgY+CjwJFJbpbk5sCXgUPoemVHAE9n6l5F6D44O9L1nLZgkjPA/gD/VeBkYBNgB2DPJDvO\novYJS/rl7wC8Ocl9qupo4B3A4X0Pcet+3oOBv/Q1bd3X98KBZW0L/Ar4O+C/6PbhswamPwM4vqp+\n12/jQcDd+8efgAOmqPFtwNFVdVtgU7reqOaHiZOauwI7AT+aZr7dgOfRtZ+bA68amP514F7AHYGT\ngEMBqupj/fN392115zW/CbNnyMwvPwcuTvLq/kC9I/Ao4Fb99EuAbegOUA8BFnJjw7sBeCnwebpe\nyYvoDnr7AQ9Kd53n6MGz7CEbAxfMosZ3VtUVVXUN3XDeR6vqJ33P61PANcDDgO2A9avqA1V1fVV9\nAZiuN1XAAVV1XlVdDuzDTQ/WEx4K3KGq3l5V11XVcuBA4JmzqH3CW6vqmqo6DTiVblgDug/9X4fz\nktwJeDywV1X9qaouAfYdWtf5VfXBqrqhqv4MHDY0fbf+Narqsqr6UlX9uaqupgu1R09R41/oznI3\nraq/VNX3V2L7ND4Bvpzk98A5dCcgb59i3gI+UVVn923nc8CD/jqx6uCq+kNVXQu8FXhgkoVD6xo7\nQ2Ye6RvTU4En0h3w96JreOf20/9QVSf1B7SL6a697Jjk1v3046rqYVW1PV0DfDBdT+IQ4F/ozo4P\nnGL1l9L1DGYyOMS1GfDKfojo8iSXA3cF7tIv67yh3/3tSiz7nCnq2QzYZGidr6c7E5ytCwee/xHY\ncIr5NgNuBlwwsK6P0J1ZTlYzwPHABkm2TbKILsC+BJBkgyQf7S8OXwmcAGyUTHpDxmvo3sMfJzk9\nyfNXYvs0PgXs3F/TXAo8BtgmyYdz4806rxuYf7At/om+LSZZkORdSc7u28ryfp47tN+ElbOqF2c1\nJv31j6UTPyf5PvDJGX7tJicT/UFrf7oQuiOwoKpWJLkY2GqKZRwL/GOSt053gwA3He46B9inqt4x\nPFOSR9MN8wzaDDh7mmXffej5+ZPMswJYXlVbTLOcVTW83Svoema373uKM/5OVV2f5HN0vbCLga9W\n1R/6ya+kG6bbtqouTvIgumGQTLKci+h6iiRZAhyb5ISq+vUqb51Gqqq+nWR/umGt7elGGmZrN+Ap\nwA5V9dsktwUu48bey5y5C9GezDyTZHGSW/Znva8C7kR3XYD+7Pg+Sdbrr8XsB3yrqq4aWswLgRP7\n4aBLgVsluR+wPV33fTLvA24DHJLk7v36Nk3y3iTD14QmfBz4t76uJLl1kif214K+D1yXZPd+6O9p\ndENdU2468LJ+nRvTXUua7GL3j4Grkrwmya36M74tk2wzzXJn60K6IaoAVNUFdBdu35dkYb/fN0/y\nqBmWMzFk9tehst6GdGerV/bbuPdUC0jyT/2YPsAVdAeVqYJOc9e+dDe0/P0U06dqnxvSneBc1o9U\nDJ/IXUR/c8+4GTLzz3PozuAvoguFx/XDaNA1qqOA3wPL6A5YN7lukeQOwO7AfwJU1XV0PZrjgA8B\nk96J1V8HeThwLfCjfkz5WLoD3ETvY/hs+0S6az8H0J1l/RJ4bj/tWuBpdBc1L6W7AP6F4dUOPT+U\n7qD+q35Zbx+et6quB55EN3b9a7rrVB+jC8hJN22S9UzliP7fS5P8tH/+XLoLsmf223gEcOcplk1f\n44+Bq+mGDY8amLQv3fW139GF8FHT1LMN8MMkV9Hd1r57Vf1mmto1B/U3fBwCvHaqWYaeT/z8Kbrh\n5fOA0+n+vGBw3oPobr65PMlsbthpJv6nZZKkVuzJSJKaMWQkSc0YMpKkZgwZSVIz68TfySTx7gZp\nHqmqsfy1useK1TPZ+7ZO9GSqqtlj7733brp8a56fNc+3eudYzWP7OpS5vH/Gvf5Z1DDp+7ZOhIwk\naTwMGUlSM4bMalq6dOm4S1hp1tzefKsX5mfNozTu/TPu9a9qDevKX/yvExsprSXG+RX1HitWnddk\nJEmjZchIkpoxZCRJzRgykqRmDBlJUjOGjCSpGUNGktSMISNJasaQkSQ1Y8hIkpoxZCRJzRgykqRm\nDBlJUjOGjCSpGUNGktSMISNJasaQkSQ1Y8hIkpoxZCRJzRgykqRmDBlJUjOGjCSpGUNGktSMISNJ\nasaQkSQ1Y8hIkpoxZCRJzRgykqRmDBlJUjOGjCSpGUNGktSMISNJasaQkSQ1Y8hIkpoxZCRJzRgy\nkqRmDBlJUjOGjCSpGUNGktSMISNJasaQkSQ1Y8hIkpoxZCRJzRgykqRmDBlJUjOGjCSpGUNGktSM\nISNJasaQkSQ1Y8hIkpoxZCRJzRgykqRmDBlJUjOGjCSpGUNGktSMISNJasaQkSQ1Y8hIkpoxZCRJ\nzRgykqRmDBlJUjOGjCSpGUNGktSMISNJasaQkSQ1Y8hIkpoxZCRJzRgykqRmDBlJUjOGjCSpGUNG\nktSMISNJasaQkSQ1Y8hIkpoxZCRJzRgykqRmDBlJUjOGjCSpGUNGktSMISNJasaQkSQ1Y8hIkpox\nZCRJzRgykqRmDBlJUjOGjCSpGUNGktSMISNJasaQkSQ1Y8hIkpoxZCRJzRgykqRmDBlJUjOGjCSp\nGUNGktSMISNJasaQkSQ1Y8hIkpoxZCRJzRgykqRmDBlJUjPrj7sAzX1VxbJlywBYvHgxScZckaT5\nwp6MplVVPOMZb2TJktNZsuR0dt31TVTVuMuSNE9kHTlgzMuNtMOwZq0bTX2tMM6WbytZdZO+bw6X\nDfHAvvaaj++twaj5zpAZ4of6piaGy445ZksAdtrpDD772bd7XUbSrDhcphl54V8j5nDZ/DTp+2bI\nSJprDJn5adL3zbvLJEnNGDKSpGYMGUlSM4aMJKmZGUMmye5JzkxyWZLXzHbBSTZL8qzVK2/W67pv\nklOSnJjkHqNYpyRpZjPeXZbkLGCHqjp/iukLqur6SV5fCryyqp68JgqdTpLXAQuqap8pZvGOEY2d\nt4LPmneXzTGzbLsrf3dZko8A9wSOTrJnkv371w9O8pEkPwT+O8mjk5zcP05MsiHwLuCR/Wt7TLH8\nBUnek2RZklOTvLx/fYckJyU5LclBSW7ev/6QJMcn+WmSo5PcOckTgD2AlyY5bubdJY2e3wGn+Wp1\n2+5sejLLgYcATwa2qapXJDkY2BjYuaoqyZHAO6vqB0k2AK4BHgG8arqeTJKXAtsDz6yqG5LcDvgT\n8AvgMVV1dpJDgJOADwLfBp5cVZcm2RXYsapekGRv4Kqqet8Uq/LTvI5blzoNa0F22ZOZQ0477TSW\nLDmdq6/eDYCFCw/ju9/dkq222mp41tX67rIMLaCAI+rGhPoe8P4khwJfrKrzMruxgB2AD1fVDQBV\ndXmSBwLLq+rsfp5DgJcBxwIPAI7tF70AGBzCW4cOI6OzLh2c1xZz9T1bC8JPq2B1vrvsjxNPqurd\nSb4GPBH4XpJ/WInlDH8khptiBv49o6oevtKVapV5YFgz/A44zVeLFy9mp50O55hjDgO6trt48ezv\n6VojX5CZZPOqOgM4I8lDgfsA5wILZ/jVbwIvSfKtqrq+Hy77BbCoX+avgOcAxwM/B+6YZLuq+mGS\nmwH3rqoz18Q2SC0l4XOf22fg4umzDBjNC6vbdmcTMjX0GHx9wh5JtgduAE4HjuqnX5/kFOCTVfWB\nSZZ9ILAFcFqSa4GPVdWHkjwfOCLJ+sCPgY9U1bVJdgH2S7JRX/v7gYmQ8Zxbc1qSycaxpTlvddqu\nX5Apaa7xwv/85BdkSpJGayT/aVl/I8C7hl7+dVU9fRTrlySNh8NlkuYah8vmJ4fLJEmjZchIkpox\nZCRJzRgykqRmDBlJUjOGjCSpGUNGktSMISNJasaQkSQ1Y8hIkpoxZCRJzRgykqRmDBlJUjOGjCSp\nGUNGktSMISNJasaQkSQ1Y8hIkpoxZCRJzRgykqRmDBlJUjOGjCSpGUNGktSMISNJasaQkSQ1Y8hI\nkpoxZCRJzRgykqRmDBlJUjOGjCSpGUNGktSMISNJasaQkSQ1Y8hIkpoxZCRJzRgykqRmDBlJUjOG\njCSpGUNGktSMISNJasaQkSQ1Y8hIkpoxZCRJzRgykqRmDBlJUjOGjCSpGUNGktSMISNJasaQkSQ1\nY8hIkpoxZCRJzRgykqRmDBlJUjOGjCSpGUNGktSMISNJasaQkSQ1Y8hIkpoxZCRJzRgykqRmDBlJ\nUjOGjCSpGUNGktSMISNJasaQkSQ1Y8hIkpoxZCRJzRgykqRmDBlJUjOGjCSpGUNGktSMISNJasaQ\nkSQ1Y8hIkpoxZCRJzRgykqRmDBlJUjOGjCSpGUNGktSMISNJasaQkSQ1Y8hIkpoxZCRJzRgykqRm\nDBlJUjOGjCSpGUNGktSMISNJasaQkSQ1Y8hIkpoxZCRJzRgykqRmDBlJUjOGjCSpGUNGktSMISNJ\nasaQkSQ1Y8hIkpoxZCRJzRgykqRmDBlJUjOGjCSpGUNGktSMISNJasaQkSQ1Y8hIkpoxZCRJzRgy\nkqRmDBlJUjOGjCSpGUNGktSMIbOajj/++HGXsNKsub35Vi/Mz5pHadz7Z9zrX9UaDJnVNBfe+JVl\nze3Nt3phftY8SuPeP+Ne/6rWYMhIkpoxZCRJzaSqxl1Dc0nW/o2U1iJVlXGs12PF6pnsfVsnQkaS\nNB4Ol0mSmjFkJEnNGDIrKcnhSU7uH8uTnDzFfK9PckaSZUkOS3KLUdc6UMtsa75tks8nOSvJmUm2\nG3WtfR2zqrefd0E/31dHWeMkdcxYc5K7JflW3y5OT7L7OGodqGe27WKnJD9L8sskrx11neOW5BNJ\nLkqybNzrSrI0yZUD79ubGtWxxtrq+muysHVBVT1z4nmS9wBXDM+TZBHwIuB+VXVNks8CzwQOGVGZ\nNzGbmnsfAL5RVbskWR+49SjqG7YS9QLsAZwJLGxd13RmWfO1wF5VdUqSDYETk3yzqs4aVZ2DZtmW\nFwAHAI8FzgN+kuTIcdU8Jp8E9gc+NUfWdUJVPaVxHTO21SS/qapFMy3InswqShLgGcBnJpn8e7o3\naYP+YL0B3Qd0rKarOclGwCOr6hMAVXVdVV054hKHa5puH5PkrsATgAOBsdyNNGy6mqvqwqo6pX9+\nNXAWsMloK/xbM+znbYGzq+o3VXUtcDiw8yjrG7eq+g5w+RxaV/O2Psu2Oqu7xgyZVfdI4KKq+tXw\nhKq6DHgvcA5wPnBFVR074vomM2XNwD2AS5J8MslJST6eZIMR1zdsunoB3g+8GrhhdCXNaKaagb/2\ndrcGfjSCmmYyXc2bAisGfj63f03jUcDDk5ya5BtJ7t96havbVh0um0SSbwJ3nmTSG6pqYuz/WcBh\nU/z+5sCewCLgSuCIJP9cVYc2KHdinatVM11beDDw8qr6SZJ9gdcBb17jxbJG9vGTgIur6uQkS1vU\nOMk6V3cfTyxnQ+DzwB79WWIza6Bm/8ZhbjkJuFtV/THJ44EvA1u0WtlwW03yRmCXfvImA9fxvltV\nr5h0IVXlYyUfdAfkC4FNppi+K3DgwM/PAT44x2u+M7B84OdHAF+bw/W+g+4MezlwAfAH4FNzeR/3\n89wMOAbYc5y1rsR+3g44euDn1wOvHXfdY9hPi4Blc21dffvfuFEd07bVwePFdA+Hy1bNY4Gzqur8\nKab/DNguya368e7H0l2cHqdpa66qC4EVSbYYmP+MURU3iZnqfUNV3a2q7kF3U8VxVfXckVb4t6at\nuW8LBwFnVtW+I61sajO15Z8C906yKMnN6U6gjhxZdbqJJHfq2xFJtqX7g/rLGqxnjbVVQ2bV7MrQ\nRdIkmyT5OkBVnUp3d8hPgdP6WT420gr/1rQ1914BHJrkVGArut7CuMym3kFzYVhnppqXAM8Gth+4\nBXWnURc5ZKa2fB3wcroz2jOBz9a6dWcZST4DfB/YIsmKJM8fwbru06/rX5O8JMlL+ll2AZYlOQXY\nl+4Eq4XZtNVZfeb8WhlJUjP2ZCRJzRgykqRmDBlJUjOGjCSpGUNGktSMISNJasaQkSQ1Y8hIkpr5\nfz/w8h0y6n1nAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x1096f1c18>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "import matplotlib.pyplot as plt\n",
    "Matplot.summary_plot(M_pooled.first_coef)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZcAAAEgCAYAAAB4qT7AAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGRNJREFUeJzt3Xm4ZHV95/H3R0Wj0mJABMGFuEWUJqJI1I6KY+RpRyMm\nMTDiGhM1mYxE4zYGJ90KqDNGh0ETI4ojKo6gRNxBicCIC0SBplEkEkHZQUAEN1C++eOcC8W16nbd\nW+dWne77fj1PPV1V59Q5v/Or5XOW7/11qgpJkrp0h1k3QJK05TFcJEmdM1wkSZ0zXCRJnTNcJEmd\nM1wkSZ0zXLRZSfKBJAe395+Q5DsD0y5K8pQRr9s7ycXTauesJDklyZ/Nuh2azJbwPhouW6Akuyb5\nUpIfJflukmcNTNslyS1Jbhi4HTQw/YAklyW5MMneA88/KMlXkmQT675PkiPbZfw4yXlJ1ie5W0eb\nV+2NqvpyVT1s2LQuDQbaGPOuT/KhrtuwCMvSB5pcu/Pz0/Y7d0WSDyW5x4jZl/w+LubzupwMly1M\nkjsBnwQ+Bfwm8FLgw0keMm/We1TVqvZ26MBr3wLsAfw34J0D8x8OvKIW+KvbJNsCXwPuAjy2qu4B\nPBXYBnjQiLYuxYIBtwym9oOd5I7TWI9mooBnVNUq4HeA1cAbZtuk5WO4bHkeBtynqg6rxsnAV4Dn\nz5tv2Hu/HXBpVV0J/AvwQIAkzwYurqp/3cS6/wa4vqqeV1U/AKiqS6rqlVW1sV3WLUn+a5LvAue3\nzz0jydlJrmuPjlbPLTDJHknObI+CPgr8xsC0Yae69kryrSTXJnl/krsMa2iSnZIcl+SqJN9L8vJN\nbFva180d+b0gyfeTXJ3kb9tpa4HXA/u3e6dntc9vM3A0d0mSg5PcoZ32onab35Hkh8DBbT88YqCt\n27d7vPdK8ptJPtO2+9okn06y84htfHCSU9sj2Kvb/lMPtN+xLwCPWGC2XZKc1n72T0yy3dyEJB9L\ncnn73p6a5OHt8y8FDgBe234GP7msG7IAw2VluAOw27znvp/k4vYHeO5DezWwXftj9VTg3CRbAwfR\n/Ghuyu8D/zzGfPsCjwEenmQP4EjgJcC2wHuATyXZKsmdgeOBo2iOwj4G/DGjjyJC88Xah+ZI6aEM\n2TNsf9g/DZwF7AQ8BXhFkn3GaPucNe3ynwL8XZLfrqoTgDcDH22PCPdo5/0AcFPbpj3a9v35wLL2\nAv4duDfwJpo+fM7A9P2AU6rqh+02Hgncv739DHjXiDYeDJxQVfcEdqY5+tRsze2k3BdYC5y+wHwH\nAC+i+VzcGXj1wPTPAg8GtgfOBI4GqKoj2vv/s/0M7tv9JozHcNnynA9cleQ17Q/0PsATgbu2068G\n9qT5YXo0sIrbPpi3AH8JfJzmKOQlND92hwOPTHMd54TBvep5tgUuH6ONb6mqH1XVL2hO272nqv61\nPdL6IPAL4HHAY4E7VdX/qapfVdVxwEJHTwW8q6ourarrgEO5/Y/0nMcA96qqQ6rql1V1IfA+4L+M\n0fY5b6yqX1TVOcAGmtMc0Pwo3HraLskOwNOAV1bVz6rqauCweeu6rKr+oapuqaqfAx+ZN/2A9jmq\n6tqq+kRV/byqbqQJsyeNaONNNHu/O1fVTVX11UVsn7oX4PgkPwZ+QLNDcciIeQt4f1Vd0H4mjgUe\neevEqg9U1U+q6mbgjcDvJFk1b10zZbhsYdoP27OAp9P80L+S5oN5STv9J1V1ZvtDdhXNtZV9kty9\nnf6lqnpcVT2Z5gP6KJojh6OAF9LsDb9vxOqvoTkS2JTBU1kPAF7Vngq6Lsl1wH2B+7TLunTea7+/\niGX/YER7HgDsNG+dr6fZQxzXFQP3fwpsPWK+BwBbAZcPrOufaPY4h7UZ4BTgbkn2SrILTXB9AiDJ\n3ZK8p704fD1wKrBNMrTQ4rU07+EZSc5N8qeL2D51r4B922uRewP/CdgzybtzW3HNfx+Yf/Az9jPa\nz1iSOyZ5a5IL2s/Ahe0891r+TRjfUi+oqsfa6xt7zz1O8lXg/27iZbfb0Wh/rN5JEz7bA3esqouT\nXAXsPmIZJwF/mOSNC1345/antX4AHFpVb54/U5In0ZzOGfQA4IIFln3/efcvGzLPxcCFVfXQBZaz\nVPO3+2KaI7Ht2iPDTb6mqn6V5Fiao66rgE9X1U/aya+iOR23V1VdleSRNKdFMmQ5V9IcGZJkDXBS\nklOr6ntL3jp1oqr+f5J30py+ejLNGYNxHQA8E3hKVX0/yT2Ba7ntaKUX1YIeuWyBkqxO8hvtXu6r\ngR1ozvvT7g3/dpI7tNdaDgdOrqob5i3mz4Fvtqd9rgHummRX4Mk0h/PDvAO4B3BUkvu369s5yduT\nzL/mM+e9wF+07UqSuyd5enut56vAL5Mc2J7i+yOaU1ojNx34q3ad29JcKxp2EfsM4IYkr01y13ZP\ncLckey6w3HFdQXMqKgBVdTnNhdt3JFnV9vuDkjxxE8uZOzV26ymx1tY0e7HXt9u4btQCkvxJe24f\n4Ec0PzqjAk7TdxhNAcrvjpg+6nO3Nc0Oy7XtGYf5O2ZX0hbjzJLhsmV6Ps0e+5U0YfDU9nQZNB+6\nzwM/BjbS/FDd7rpEknsBBwL/A6CqfklzBPMl4B+BoZVV7XWOxwM3A6e355ZPovlhmzvamL93/U2a\nazvvotn7+i7wgnbazcAf0VzUvIbmwvZx81c77/7RND/m/94u65D581bVr4Bn0JzD/h7NdagjaIJx\n6KYNWc8oH2v/vSbJN9r7L6C5IPvtdhs/Buw4Ytm0bTwDuJHm9ODnByYdRnP97Ic04fv5BdqzJ/D1\nJDfQlKcfWFUXLdB2TVFboHEU8LpRs8y7P/f4gzSnhy8FzqUp/x+c90iaYpnrkoxTYLMs4n8WJknq\nmkcukqTOGS6SpM4ZLpKkzhkukqTO+XcuY0hi1YM0I1U11b829/u+OKPeH49cxlBVU7mtW7duauva\n1G2a27259Etf2rHC2jL1YUz62O+zfL83se6R74/hIknqnOEiSeqc4dIje++996yb0Et96Ze+tANs\ny6zMaltn2cdLXbd/oT+eFddJ69evZ/369bNuhjSLoeNX3Pd9AiPfH8NlPHaSNBuGS795QV+SND2G\niySpc4aLJKlzhoskqXOGi4ayUkzSJKwWG8+K66Qk+NlQD1gt1m9Wi0mSpsdwkSR1znCRJHXOcJEk\ndc5w0VDr1q2bdRMkbcasFhuPnSTNhtVi/Wa1mCRpegwXSVLnDBdJUucMF0lS5wwXDeXYYpImYbXY\neFZcJzm2mHrCarF+s1pMkjQ9hoskqXOGiySpc4aLJKlzhouGcmwxSZOwWmw8dpI0G1aL9ZvVYpKk\n6TFcJEmdM1wkSZ0zXCRJnTNcNJRji0mahNVi41lxneTYYuoJq8X6zWoxSdL0GC6SpM4ZLpKkzhku\nkqTOGS4ayrHFJE3CarHx2EnSbFgt1m9Wi0mSpsdwkSR1znCRJHXOcJEkdc5w0VCOLSZpElaLjWfF\ndZJji6knrBbrN6vFJEnTY7hIkjpnuEiSOnenWTdA0mSqio0bNwKwevVqkllcppBuzyMXDeXYYpuH\nqmK//Q5izZpzWbPmXPbf/w0WYqgXrBYbj52kZTfrA46e/hRYLdZvI98fT4tpxZn1j3hfTbtfehpm\n6ojhohVnS/pRmzstduKJuwGwdu23OOaYQ7zuoplb0jWXJO9NsmvXjZlEkgOTfDvJh+Y9/6IkVyc5\non28XZKTk9yQ5J3z5p17/tHTbLu0VEk49thDOe203TjttN0MFvXGFnPNJcl5wFOq6rJ5z78QeHRV\nHdg+vhuwB7AbsFtVvXze/CcDr6qqMwee3jI6SVusLbhirBfXXLbg/p3U0v5CP8ndk3w2ydlJNib5\nk/b5U5I8OskfJDmrvZ2f5Hvt9Ee383wjyQlJdlxgHQ9OclK7jm8m+a32+be16zwnyX4D878myRlJ\nNiRZ3z73T8ADgROSvGKhDqiqn1bVV4BfLLTtK51ji20+rBhbXvbvElXVyBvwx8ARA4/v0f57MvCo\nefMeA/wlzXWcrwLbtc/vDxy5wDpOB/Zt798ZuGu73i/QhMK9ge8DOwL7AO9p570D8BngCe3jC4Ft\nhyz/hcA7F/H8r21brUDNR0PT1FwNmv2tZxb8jVqm2+1s2LChtt766Fv7Z9Wqo2vDhg3Ltb2bm5H9\nuKkL+ucAf5/krcBnquq0YTMleS3w06p6d5LdgEcAJ7WHjncELhvxulXATlX1yTbobmqfXwN8pKoK\nuCrJqcBjgCcB+yQ5q13E3YEHA1/exHZIC/Isx22Wuy/c6V8ZFgyXqvpukj2ApwOHJPmXqjp4cJ4k\nv09zpPHEuaeAb1XV4yds26iP+Fuq6ogJly3dzub8g1dWjC2r1atXs3btRznxxI8ATf+uXv2cGbeq\n/xYMlyT3Aa6rqqOTXA+8eN70BwD/AOxTVXPXMM4Htk/y2Kr6epKtgIdU1bfnL7+qbkhySZJ9q+qT\nSe5Cc7rry8DLkhwFbEcTXK8Gfg4cnOToqvpJkp2Bm6rq6iVsu988bRHmKsZuu+D8HIOlQ/bv0mzq\ntNhq4G1JbgFuBv5iYFporltsCxzfdvalVfWMJM8GDk+yTbuO/w38Wri0ng+8J8mb2nU8u6o+keRx\nwAaayo3XVNVVwBfbEuivteu7EXgucDWLqOhKchGwCrhzkn1pwvE7475e6psk7L777rNuxhbL/l28\nLaYUeZS2FHnPmldyvMD8liLTVItZMaYe6EUpskYa+f6shHDZD3gz8KWqeukm5j0Z+C3gD6pq48Ck\nLbuTpP4yXPpt9uGS5F3AmnlPH1ZVR02lAZPxwybNhuHSb7MPl82cnSTNhuHSb0v7C31JkpbCcJEk\ndc5w0VBWikmahNdcxrPiOimJg/OpD7zm0m9ec5EkTY/hIknqnOEiSeqc4SJJ6pzhoqHWrVs36yZI\n2oxZLTYeO0maDavF+s1qMUnS9BgukqTOGS6SpM4ZLpKkzhkuGsqxxSRNwmqx8ay4TnJsMfWE1WL9\nZrWYJGl6DBdJUucMF0lS5wwXSVLnDBcN5dhikiZhtdh47CRpNqwW6zerxSRJ02O4SJI6Z7hIkjpn\nuEiSOme4aCjHFpM0CavFxrPiOsmxxdQTVov1m9VikqTpMVwkSZ0zXCRJnTNcJEmdM1w0lGOLSZqE\n1WLjsZOk2bBarN+sFpMkTY/hIknqnOEiSeqc4SJJ6pzhoqEcW0zSJKwWG8+K6yTHFlNPWC3Wb1aL\nSZKmx3CRJHXOcJEkdc5wkSR1znDRUI4tJmkSVouNx06SZsNqsX6zWkySND2GiySpc4aLJKlzhosk\nqXOGi4ZybDFJk7BabDwrrpMcW0w9YbVYv1ktJkmaHsNFktQ5w0WS1DnDRZLUOcNFQzm2mKRJWC02\nHjtJmg2rxfrNajFJ0vQYLpKkzhkukqTOGS6SpM4ZLhrKscUkTcJqsfGsuE5ybDH1hNVi/Wa1mCRp\negwXSVLnDBdJUucMF0lS5wwXDeXYYpImYbXYeOwkaTasFus3q8UkSdNjuEiSOme4SJI6Z7hIkjpn\nuGgoxxaTNAmrxcaz4jrJscXUE1aL9ZvVYpKk6TFcJEmdM1wkSZ0zXCRJnTNcNJRji0mahNVi47GT\npNmwWqzfrBaTJE2P4SJJ6pzhIknqnOEiSeqc4aKhHFtM0iSsFhvPiuskxxZTT1gt1m9Wi0mSpsdw\nkSR1znCRJHXOcJEkdc5w0VCOLSZpElaLjcdOkmbDarF+s1pMkjQ9hoskqXOGiySpc4aLJKlzhouG\ncmwxSZOwWmw8K66THFtMPWG1WL9ZLSZJmh7DRZLUOcNFktQ5w0WS1DnDRUM5tpikSVgtNp4V0UlV\nxcaNGwFYvXo1ySwKdaTbsVqs36wW08Kqiv32O4g1a85lzZpz2X//N1iKLGnJPHIZz2bRSX080PDj\npQl55NJvI9+fO02zFStJH3/oZ6GP/WDgScvPcFkmm9sP2NxpsRNP3A2AtWu/xTHHHOJ1F0lLYrgI\naIZ7OfbYQ2+9oH/ccecbLJKWbEVcc0myC3Bee3sm8CHg3jTnVo+oqsPb+d4GPA/4+6p6+8AitvxO\nmsexxdQTXnPpN6+5ABdU1aOS7Ai8sqrOTrI18M0kX6yq86rqNUlunHVDZ2F+GbKk21imv3grKVwA\nqKorgCva+zcmOQ/YieaoZkWau95ywgnN9ZanPe2YGbdI6o9h3w+vR27aSjot9umqWj3k+VOBR1TV\nje1z64AbN5fTYsv3+Q6TbvYK+Ghp+c38tNg555zDmjXncuONBwCwatVHOO203dh9991n0LTe8bTY\nfO0psY8Dfz0XLN0st6slbf5m2RcGmzRbKzJckmwFHAd8uKqO73LZm+OP2rAy5F13/Tve+MYZN0zq\ngdWrV7N27Uc58cSPAM33Y/Xq58y4Vf234k6LpTlRehRwTVW9csi8m9Vpsa54wVI9NfPTYuD3YwGe\nFhuwhqbc+JwkZ7XP/W1VfX6GbZq5JJ5Dlkbw+7F4Ky5cquo0Fh6w010SSZrQShkV+ZfANknOXGim\n9o8onwusyL91kaSurIhrLh2wk6TZ6MU1F43k/+eixVm/fv2smyBpM+aRy3hWXCc5tph6wiOXfvPI\nRZI0PYaLJKlzhoskqXOGiySpc4aLhlq3bt2smyBpM2a12HjsJGk2rBbrN6vFJEnTY7hIkjpnuEiS\nOme4SJI6Z7hoKMcWkzQJq8XGs+I6ybHF1BNWi/Wb1WKSpOkxXCRJnTNcJEmdM1wkSZ0zXDSUY4tJ\nmoTVYuOxk6TZsFqs36wWkyRNj+EiSeqc4SJJ6pzhIknqnOGioRxbTNIkrBYbz4rrJMcWU09YLdZv\nVotJkqbHcJEkdc5wkSR1znCRJHXOcNFQji0maRJWi43HTpJmw2qxfrNaTJI0PYaLJKlzhoskqXOG\niySpc4aLhnJsMUmTsFpsPCuukxxbTD1htVi/WS0mSZoew0WS1DnDRZLUOcNFktQ5w0VDObaYpElY\nLTYeO0maDavF+s1qMUnS9BgukqTOGS6SpM4ZLpKkzhkuGsqxxSRNwmqx8ay4TnJsMfWE1WL9ZrXY\n5uCUU06ZdRN6qS/90pd2gG2ZlVlt6yz7eKnrNlx6ZCV9SRejL/3Sl3aAbZkVw2V8hoskqXOGiySp\nc17QH0MSO0makaqa6kV9v++LM+r9MVwkSZ3ztJgkqXOGiySpc4bLlCXZNskXk/xbki8kueeI+dYm\n+U6S7yZ53cDz65NckuSs9rZ2kesfutx58xzeTt+QZI/FvHaKbbkoyTltH5yx3G1J8rAkX0vy8ySv\nWux2TKkd0+6T57bvyzlJvpJk93FfuzlJ8v4kVybZOMv1JNk7yfUD3/03dLz++yU5Ocm3kpyb5MCJ\nFlhV3qZ4A/4X8Nr2/uuAtw6Z547ABcAuwFbA2cCu7bR1wN8scd0jlzswz38GPtfe/13g6+O+dlpt\naR9fCGzb0XsyTlu2B/YEDgFetZjXTqMdM+qTxwHbtPfXLtdnZdY34AnAHsDGWa4H2Bv41DKuf0fg\nke39rYHzh7znF427PI9cpu+ZwFHt/aOAZw2ZZy/ggqq6qKpuBj4K7DswfanVM5ta7u3aV1WnA/dM\nsuOYr51GW3YYmN5VFdEm21JVV1fVN4Cbl7Ad02jHnGn2ydeq6vr24enAfcd97eakqr4MXNeT9Sxb\n5VxVXVFVZ7f3bwTOA3aaP9u4yzNcpm+HqrqyvX8lsMOQeXYGLh54fEn73JyXt6cjjhx1Wm2ETS13\noXl2GuO1izFJW6D5kJ+U5BtJXjJBO8Zty3K8tutlzbJP/gz43BJfq/EU8Pj2u/+5JA9frhUl2YXm\nKOr0pS7jTl01RrdJ8kWaQ8z5Dhp8UFU1oqZ+ob2DdwNvau8fDLyd5os9jnH3OqbxdwWTtuX3quqy\nJNsDX0zynXbPbznb0vVru17Wmqq6fNp9kuTJwIuBNYt9rRblTOB+VfXTJE8Djgce2vVKkmwNfBz4\n66q6MclBwLPbyTslOau9f1pVvXzUcgyXZVBVTx01rb1gt2NVXZHkPsBVQ2a7FLjfwOP70ez9UVW3\nzp/kfcCnF9G0kctdYJ77tvNsNcZrF2OpbbkUoKoua/+9OsknaE7FLPWHdJy2LMdrO11WVV3e/ju1\nPmkv4r8XWFtV1y3mtVqcqrph4P7nk/xjkm2r6tqu1pFkK+A44MNVdXy7rkOBQ9vpF1bVHgss4lae\nFpu+TwEvbO+/kGbvY75vAA9JskuSOwP7t6+jDaQ5fwgspoJl5HLnte8F7boeC/yoPY03zmsXY8lt\nSXK3JKva5+8O7MPi+mEpbZkz/0iqy35Zcjtm0SdJ7g/8M/C8qrpgiduhMSXZIUna+3vR/BF8l8ES\n4Ejg21V12MQLXM7qB29DKzK2BU4C/g34AnDP9vmdgM8OzPc0mmqNC4DXDzz/QeAcYANNMO2wyPX/\n2nKBlwEvG5jnXe30DcCjNtWmCfpiSW0BHkhTgXQ2cO402kJzmvNi4Hqai64/ALbuul+W2o4Z9cn7\ngGuAs9rbGcv1WZnlDfh/wGXAL9q+/9NlXs9N7XpePK+//6p9b88Gvgo8tuP1/x5wS7v8ufd07bx5\nvjfu8hz+RZLUOU+LSZI6Z7hIkjpnuEiSOme4SJI6Z7hIkjpnuEiSOme4SJI6Z7hIkjr3H/WSLWwT\nfvCOAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x10bd0cda0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "import matplotlib.pyplot as plt\n",
    "Matplot.summary_plot(M_pooled.size_coef)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.4.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}