Skip to content

Instantly share code, notes, and snippets.

@ethanwhite
Last active March 1, 2016 17:39
Show Gist options
  • Save ethanwhite/faac5420e4ae07844233 to your computer and use it in GitHub Desktop.
Save ethanwhite/faac5420e4ae07844233 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"import numpy as np\n",
"import pandas as pd\n",
"import matplotlib.pyplot as plt\n",
"\n",
"%matplotlib inline"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"cbcdata15 = pd.read_csv(\"opttestscipy15/cbc_likelihoods.csv\")\n",
"cbcdata16 = pd.read_csv(\"opttestscipy16/cbc_likelihoods.csv\")\n",
"cbcdata17 = pd.read_csv(\"opttestscipy17/cbc_likelihoods.csv\")\n",
"cbcdata17b = pd.read_csv(\"opttestscipy17b/cbc_likelihoods.csv\")"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/html": [
"<div style=\"max-height:1000px;max-width:1500px;overflow:auto;\">\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>site</th>\n",
" <th>S</th>\n",
" <th>N</th>\n",
" <th>likelihood_logseries</th>\n",
" <th>likelihood_pln</th>\n",
" <th>likelihood_negbin</th>\n",
" <th>likelihood_zipf</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td> L105249</td>\n",
" <td> 57</td>\n",
" <td> 7400</td>\n",
" <td>-315.187361</td>\n",
" <td>-315.811735</td>\n",
" <td>-311.991052</td>\n",
" <td>-334.021427</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td> L105252</td>\n",
" <td> 57</td>\n",
" <td> 3730</td>\n",
" <td>-262.599618</td>\n",
" <td>-264.710188</td>\n",
" <td>-262.482895</td>\n",
" <td>-274.556152</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2</th>\n",
" <td> L105254</td>\n",
" <td> 60</td>\n",
" <td> 17359</td>\n",
" <td>-328.224545</td>\n",
" <td>-324.407750</td>\n",
" <td>-328.224533</td>\n",
" <td>-338.805827</td>\n",
" </tr>\n",
" <tr>\n",
" <th>3</th>\n",
" <td> L105355</td>\n",
" <td> 64</td>\n",
" <td> 12490</td>\n",
" <td>-367.702703</td>\n",
" <td>-364.746419</td>\n",
" <td>-365.366943</td>\n",
" <td>-387.606697</td>\n",
" </tr>\n",
" <tr>\n",
" <th>4</th>\n",
" <td> L105358</td>\n",
" <td> 59</td>\n",
" <td> 5622</td>\n",
" <td>-298.761115</td>\n",
" <td>-293.324891</td>\n",
" <td>-297.554643</td>\n",
" <td>-314.489174</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" site S N likelihood_logseries likelihood_pln \\\n",
"0 L105249 57 7400 -315.187361 -315.811735 \n",
"1 L105252 57 3730 -262.599618 -264.710188 \n",
"2 L105254 60 17359 -328.224545 -324.407750 \n",
"3 L105355 64 12490 -367.702703 -364.746419 \n",
"4 L105358 59 5622 -298.761115 -293.324891 \n",
"\n",
" likelihood_negbin likelihood_zipf \n",
"0 -311.991052 -334.021427 \n",
"1 -262.482895 -274.556152 \n",
"2 -328.224533 -338.805827 \n",
"3 -365.366943 -387.606697 \n",
"4 -297.554643 -314.489174 "
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"cbcdata15.head()"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/html": [
"<div style=\"max-height:1000px;max-width:1500px;overflow:auto;\">\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>site</th>\n",
" <th>S</th>\n",
" <th>N</th>\n",
" <th>likelihood_logseries</th>\n",
" <th>likelihood_pln</th>\n",
" <th>likelihood_negbin</th>\n",
" <th>likelihood_zipf</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td> L105249</td>\n",
" <td> 57</td>\n",
" <td> 7400</td>\n",
" <td>-315.187361</td>\n",
" <td>-315.811735</td>\n",
" <td>-315.185351</td>\n",
" <td>-334.021427</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td> L105252</td>\n",
" <td> 57</td>\n",
" <td> 3730</td>\n",
" <td>-262.599618</td>\n",
" <td>-264.710188</td>\n",
" <td>-262.482895</td>\n",
" <td>-274.556152</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2</th>\n",
" <td> L105254</td>\n",
" <td> 60</td>\n",
" <td> 17359</td>\n",
" <td>-328.224545</td>\n",
" <td>-324.407750</td>\n",
" <td>-328.225260</td>\n",
" <td>-338.805827</td>\n",
" </tr>\n",
" <tr>\n",
" <th>3</th>\n",
" <td> L105355</td>\n",
" <td> 64</td>\n",
" <td> 12490</td>\n",
" <td>-367.702703</td>\n",
" <td>-364.746419</td>\n",
" <td>-367.701879</td>\n",
" <td>-387.606697</td>\n",
" </tr>\n",
" <tr>\n",
" <th>4</th>\n",
" <td> L105358</td>\n",
" <td> 59</td>\n",
" <td> 5622</td>\n",
" <td>-298.761115</td>\n",
" <td>-293.324891</td>\n",
" <td>-298.761042</td>\n",
" <td>-314.489174</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" site S N likelihood_logseries likelihood_pln \\\n",
"0 L105249 57 7400 -315.187361 -315.811735 \n",
"1 L105252 57 3730 -262.599618 -264.710188 \n",
"2 L105254 60 17359 -328.224545 -324.407750 \n",
"3 L105355 64 12490 -367.702703 -364.746419 \n",
"4 L105358 59 5622 -298.761115 -293.324891 \n",
"\n",
" likelihood_negbin likelihood_zipf \n",
"0 -315.185351 -334.021427 \n",
"1 -262.482895 -274.556152 \n",
"2 -328.225260 -338.805827 \n",
"3 -367.701879 -387.606697 \n",
"4 -298.761042 -314.489174 "
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"cbcdata16.head()"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/html": [
"<div style=\"max-height:1000px;max-width:1500px;overflow:auto;\">\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>site</th>\n",
" <th>S</th>\n",
" <th>N</th>\n",
" <th>likelihood_logseries</th>\n",
" <th>likelihood_pln</th>\n",
" <th>likelihood_negbin</th>\n",
" <th>likelihood_zipf</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td> L105249</td>\n",
" <td> 57</td>\n",
" <td> 7400</td>\n",
" <td>-315.187361</td>\n",
" <td>-315.811735</td>\n",
" <td>-315.185351</td>\n",
" <td>-334.021427</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td> L105252</td>\n",
" <td> 57</td>\n",
" <td> 3730</td>\n",
" <td>-262.599618</td>\n",
" <td>-264.710188</td>\n",
" <td>-262.482895</td>\n",
" <td>-274.556152</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2</th>\n",
" <td> L105254</td>\n",
" <td> 60</td>\n",
" <td> 17359</td>\n",
" <td>-328.224545</td>\n",
" <td>-324.407750</td>\n",
" <td>-328.225260</td>\n",
" <td>-338.805827</td>\n",
" </tr>\n",
" <tr>\n",
" <th>3</th>\n",
" <td> L105355</td>\n",
" <td> 64</td>\n",
" <td> 12490</td>\n",
" <td>-367.702703</td>\n",
" <td>-364.746419</td>\n",
" <td>-367.701879</td>\n",
" <td>-387.606697</td>\n",
" </tr>\n",
" <tr>\n",
" <th>4</th>\n",
" <td> L105358</td>\n",
" <td> 59</td>\n",
" <td> 5622</td>\n",
" <td>-298.761115</td>\n",
" <td>-293.324891</td>\n",
" <td>-298.761042</td>\n",
" <td>-314.489174</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" site S N likelihood_logseries likelihood_pln \\\n",
"0 L105249 57 7400 -315.187361 -315.811735 \n",
"1 L105252 57 3730 -262.599618 -264.710188 \n",
"2 L105254 60 17359 -328.224545 -324.407750 \n",
"3 L105355 64 12490 -367.702703 -364.746419 \n",
"4 L105358 59 5622 -298.761115 -293.324891 \n",
"\n",
" likelihood_negbin likelihood_zipf \n",
"0 -315.185351 -334.021427 \n",
"1 -262.482895 -274.556152 \n",
"2 -328.225260 -338.805827 \n",
"3 -367.701879 -387.606697 \n",
"4 -298.761042 -314.489174 "
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"cbcdata17.head()"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/html": [
"<div style=\"max-height:1000px;max-width:1500px;overflow:auto;\">\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>site</th>\n",
" <th>S</th>\n",
" <th>N</th>\n",
" <th>likelihood_logseries</th>\n",
" <th>likelihood_pln</th>\n",
" <th>likelihood_negbin</th>\n",
" <th>likelihood_zipf</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td> L105249</td>\n",
" <td> 57</td>\n",
" <td> 7400</td>\n",
" <td>-315.187361</td>\n",
" <td>-315.811735</td>\n",
" <td>-315.185351</td>\n",
" <td>-334.021427</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td> L105252</td>\n",
" <td> 57</td>\n",
" <td> 3730</td>\n",
" <td>-262.599618</td>\n",
" <td>-264.710188</td>\n",
" <td>-262.482895</td>\n",
" <td>-274.556152</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2</th>\n",
" <td> L105254</td>\n",
" <td> 60</td>\n",
" <td> 17359</td>\n",
" <td>-328.224545</td>\n",
" <td>-324.407750</td>\n",
" <td>-328.225260</td>\n",
" <td>-338.805827</td>\n",
" </tr>\n",
" <tr>\n",
" <th>3</th>\n",
" <td> L105355</td>\n",
" <td> 64</td>\n",
" <td> 12490</td>\n",
" <td>-367.702703</td>\n",
" <td>-364.746419</td>\n",
" <td>-367.701879</td>\n",
" <td>-387.606697</td>\n",
" </tr>\n",
" <tr>\n",
" <th>4</th>\n",
" <td> L105358</td>\n",
" <td> 59</td>\n",
" <td> 5622</td>\n",
" <td>-298.761115</td>\n",
" <td>-293.324891</td>\n",
" <td>-298.761042</td>\n",
" <td>-314.489174</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" site S N likelihood_logseries likelihood_pln \\\n",
"0 L105249 57 7400 -315.187361 -315.811735 \n",
"1 L105252 57 3730 -262.599618 -264.710188 \n",
"2 L105254 60 17359 -328.224545 -324.407750 \n",
"3 L105355 64 12490 -367.702703 -364.746419 \n",
"4 L105358 59 5622 -298.761115 -293.324891 \n",
"\n",
" likelihood_negbin likelihood_zipf \n",
"0 -315.185351 -334.021427 \n",
"1 -262.482895 -274.556152 \n",
"2 -328.225260 -338.805827 \n",
"3 -367.701879 -387.606697 \n",
"4 -298.761042 -314.489174 "
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"cbcdata17b.head()"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"'CBC Likelihood Difference. Min: -5.52335648243, Max: 104.347664841'"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYEAAAEACAYAAABVtcpZAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAD6pJREFUeJzt3X2onnd9x/H3Z8bOqp1ZcaQPCTR/pMNI1Vowbm54xC6L\nsjXdxmwLlnQGGaRT6R9jjWNrskEQhpuFkTJmW2OZGUFriNK5xM4DzoFRVzWaxibQ2ObMnLrNTR2U\npvjdH/eVentyHpJzTu6H/t4vuOnv+l1P36uc8/vc19NJqgpJUpt+btgFSJKGxxCQpIYZApLUMENA\nkhpmCEhSwwwBSWrYvCGQZE2SLyT5dpJvJXl/1395kkNJnkhyMMnKvnW2Jzme5FiSjX39NyQ50s27\n9+IdkiTpfC10JnAGuKuqXgu8GbgzyWuAu4FDVXUt8Gg3TZL1wC3AemATsDtJum3dB2ytqnXAuiSb\nlv1oJEkXZN4QqKrTVfX1rv1j4HHgauAmYE+32B7g5q69GdhbVWeq6iRwAtiQ5Ergsqo63C338b51\nJElDct73BJJcA1wPfBlYVVXT3axpYFXXvgo41bfaKXqhMbN/quuXJA3ReYVAklcCnwI+UFU/6p9X\nvb874d+ekKQxtGKhBZK8lF4APFRV+7vu6SRXVNXp7lLPM13/FLCmb/XV9M4Aprp2f//ULPsyTCRp\nEaoqCy91roWeDgpwP3C0qj7SN+sAsKVrbwH29/XfmuSSJGuBdcDhqjoN/DDJhm6bt/etM/NAxvZz\nzz33DL0G6x9+HS3WP861vxjqX4qFzgTeArwb+GaSx7q+7cCHgH1JtgIngXd1A/jRJPuAo8DzwLb6\naYXbgI8BlwKPVNXnllS5JGnJ5g2BqvpX5j5buHGOdXYBu2bp/xpw3YUWKEm6eHxjeBlNTEwMu4Ql\nsf7hGuf6x7l2GP/6lyJLvZ60nJLUKNUjSeMgCXUxbgxLkl7cDAFJapghIEkNMwQkqWGGgCQ1zBCQ\npIYZApLUMENAkhpmCEhSwwwBSWqYISBJDVvwH5Vp0Z133sVnP/voQPZ1xx2/z86dfzaQfUnSTIbA\nLI4d+y5PPXUHc/y17GW0jyefPOcfWJOkgTEE5nQN8LqLvI9/A/7zIu9DkubmPQFJapghIEkNMwQk\nqWGGgCQ1zBCQpIYZApLUMENAkhpmCEhSwwwBSWqYISBJDTMEJKlhhoAkNcwQkKSGGQKS1DBDQJIa\nZghIUsMMAUlqmCEgSQ0zBCSpYYaAJDXMEJCkhhkCktQwQ0CSGmYISFLDDAFJapghIEkNMwQkqWGG\ngCQ1zBCQpIYZApLUMENAkhq2YAgkeSDJdJIjfX07kpxK8lj3eUffvO1Jjic5lmRjX/8NSY508+5d\n/kORJF2o8zkTeBDYNKOvgL+uquu7zz8BJFkP3AKs79bZnSTdOvcBW6tqHbAuycxtSpIGbMEQqKov\nAj+YZVZm6dsM7K2qM1V1EjgBbEhyJXBZVR3ulvs4cPPiSpYkLZel3BN4X5JvJLk/ycqu7yrgVN8y\np4CrZ+mf6volSUO0YpHr3Qf8Rdf+S+DDwNblKGjHjh0vtCcmJpiYmFiOzUrSi8bk5CSTk5PLsq1F\nhUBVPXO2neSjwGe6ySlgTd+iq+mdAUx17f7+qdm23R8CkqRzzfyCvHPnzkVva1GXg7pr/Gf9DnD2\nyaEDwK1JLkmyFlgHHK6q08APk2zobhTfDuxfdNWSpGWx4JlAkr3AW4FXJ3kauAeYSPIGek8JPQn8\nIUBVHU2yDzgKPA9sq6rqNrUN+BhwKfBIVX1umY9FknSBFgyBqrptlu4H5ll+F7Brlv6vAdddUHWS\npIvKN4YlqWGGgCQ1zBCQpIYZApLUMENAkhpmCEhSwwwBSWqYISBJDTMEJKlhhoAkNcwQkKSGGQKS\n1DBDQJIaZghIUsMMAUlqmCEgSQ0zBCSpYYaAJDXMEJCkhhkCktQwQ0CSGmYISFLDDAFJapghIEkN\nMwQkqWGGgCQ1zBCQpIYZApLUMENAkhpmCEhSwwwBSWqYISBJDTMEJKlhhoAkNcwQkKSGGQKS1DBD\nQJIaZghIUsMMAUlqmCEgSQ0zBCSpYYaAJDXMEJCkhhkCktQwQ0CSGmYISFLDDAFJatiCIZDkgSTT\nSY709V2e5FCSJ5IcTLKyb972JMeTHEuysa//hiRHunn3Lv+hSJIu1PmcCTwIbJrRdzdwqKquBR7t\npkmyHrgFWN+tsztJunXuA7ZW1TpgXZKZ25QkDdiCIVBVXwR+MKP7JmBP194D3Ny1NwN7q+pMVZ0E\nTgAbklwJXFZVh7vlPt63jiRpSBZ7T2BVVU137WlgVde+CjjVt9wp4OpZ+qe6fknSEK1Y6gaqqpLU\nchQDsGPHjhfaExMTTExMLNemJelFYXJyksnJyWXZ1mJDYDrJFVV1urvU80zXPwWs6VtuNb0zgKmu\n3d8/NduG+0NAknSumV+Qd+7cuehtLfZy0AFgS9feAuzv6781ySVJ1gLrgMNVdRr4YZIN3Y3i2/vW\nkSQNyYJnAkn2Am8FXp3kaeDPgQ8B+5JsBU4C7wKoqqNJ9gFHgeeBbVV19lLRNuBjwKXAI1X1ueU9\nFEnShVowBKrqtjlm3TjH8ruAXbP0fw247oKqkyRdVL4xLEkNMwQkqWGGgCQ1zBCQpIYZApLUMENA\nkhpmCEhSwwwBSWqYISBJDTMEJKlhhoAkNcwQkKSGGQKS1DBDQJIaZghIUsMMAUlqmCEgSQ0zBCSp\nYYaAJDXMEJCkhhkCktQwQ0CSGmYISFLDDAFJapghIEkNMwQkqWGGgCQ1zBCQpIYZApLUMENAkhpm\nCEhSwwwBSWqYISBJDTMEJKlhhoAkNcwQkKSGGQKS1DBDQJIaZghIUsMMAUlqmCEgSQ0zBCSpYYaA\nJDXMEJCkhhkCktQwQ0CSGmYISFLDlhQCSU4m+WaSx5Ic7vouT3IoyRNJDiZZ2bf89iTHkxxLsnGp\nxUuSlmapZwIFTFTV9VX1pq7vbuBQVV0LPNpNk2Q9cAuwHtgE7E7imYgkDdFyDMKZMX0TsKdr7wFu\n7tqbgb1VdaaqTgIngDchSRqa5TgT+HySryZ5b9e3qqqmu/Y0sKprXwWc6lv3FHD1EvcvSVqCFUtc\n/y1V9b0kvwQcSnKsf2ZVVZKaZ/1z5u3YseOF9sTEBBMTE0ssUZJeXCYnJ5mcnFyWbS0pBKrqe91/\nv5/k0/Qu70wnuaKqTie5EnimW3wKWNO3+uqu72f0h4Ak6VwzvyDv3Llz0dta9OWgJC9PclnXfgWw\nETgCHAC2dIttAfZ37QPArUkuSbIWWAccXuz+JUlLt5QzgVXAp5Oc3c4/VNXBJF8F9iXZCpwE3gVQ\nVUeT7AOOAs8D26pqvktFkqSLbNEhUFVPAm+Ypf+/gRvnWGcXsGux+5QkLS+f05ekhhkCktQwQ0CS\nGmYISFLDDAFJapghIEkNMwQkqWGGgCQ1zBCQpIYZApLUMENAkhpmCEhSwwwBSWqYISBJDTMEJKlh\nhoAkNcwQkKSGGQKS1DBDQJIaZghIUsMMAUlqmCEgSQ0zBCSpYYaAJDXMEJCkhhkCktQwQ0CSGmYI\nSFLDDAFJapghIEkNMwQkqWGGgCQ1zBCQpIYZApLUMENAkhpmCEhSwwwBSWqYISBJDTMEJKlhhoAk\nNcwQkKSGGQKS1DBDYMgeeujvSDKQjyTNtGLYBQigBrAPQ0DSuTwTkKSGGQKS1LCBhkCSTUmOJTme\n5E8GuW9J0rkGFgJJXgL8LbAJWA/cluQ1g9r/YEwOu4AlmZycHHYJS2L9wzPOtcP4178UgzwTeBNw\noqpOVtUZ4B+BzQPc/wBMDruAeS309NDb3va2sX4Sadx/kce5/nGuHca//qUYZAhcDTzdN32q69PA\n1AKfe85jmfP5LBw44xw20otJqgbxeCIk+T1gU1W9t5t+N7Chqt7Xt0wNqp75bNz4u3zpS0+xYsWV\nF7Tes89+h5e97JfPe/nnnjvJs89+i8E9IrrQfnZ0n0HsazkYAhqOURin+iWhqhb1CzHI9wSmgDV9\n02vonQ38jHH/dvfcc8cXsdagjvl89rNzgPuSxtO4j1P9BnkmsAL4DvB24D+Aw8BtVfX4QAqQJJ1j\nYGcCVfV8kj8C/hl4CXC/ASBJwzWwMwFJ0ugZiTeGk/xVkseTfCPJw0le1Tdve/dy2bEkG4dZ53zG\n6UW4JGuSfCHJt5N8K8n7u/7LkxxK8kSSg0lWDrvW+SR5SZLHknymmx6b+pOsTPLJ7uf+aJINY1b/\nXd3PzpEkn0jy86Ncf5IHkkwnOdLXN2e9ozTuzFH7so2ZIxECwEHgtVX1euAJYDtAkvXALfReLtsE\n7E4yKjW/IOP3ItwZ4K6qei3wZuDOrt67gUNVdS3waDc9yj4AHOWnjyKNU/33Ao9U1WuA1wHHGJP6\nk1wNvA+4oaquo3d591ZGu/4H6f1+9pu13hEcd2arfdnGzJEYUKvqUFX9pJv8MrC6a28G9lbVmao6\nCZyg99LZqBmrF+Gq6nRVfb1r/xh4nN47GzcBe7rF9gA3D6fChSVZDbwT+Cg/fRRpLOrvvrX9elU9\nAL37ZVX1v4xJ/Z0VwMu7Bz5eTu9hj5Gtv6q+CPxgRvdc9Y7UuDNb7cs5Zo5ECMzwHuCRrn0VP/sY\n6ai+YDa2L8IluQa4nt4P0qqqmu5mTQOrhlTW+fgb4I+Bn/T1jUv9a4HvJ3kwyb8n+fskr2BM6q+q\nKeDDwFP0Bv//qapDjEn9feaqd1zGnbOWNGYO8m8HHequH878/HbfMn8KPFdVn5hnU6N4J3sUa1pQ\nklcCnwI+UFU/6p/XvbU3kseV5LeAZ6rqMeZ4IWGU66f3LfqNwO6qeiPwf8y4dDLK9Sf5RXrfoq+h\nN+i8snv58wWjXP9szqPekTyW5RgzB/mI6G/MNz/JHfRO79/e1z3zBbPVXd+oOa8X4UZJkpfSC4CH\nqmp/1z2d5IqqOp3kSuCZ4VU4r18FbkryTuBlwC8keYjxqf8UcKqqvtJNf5LeNd3TY1L/jcCTVfVf\nAEkeBn6F8an/rLl+XsZi3FmuMXMkLgcl2UTv1H5zVT3bN+sAcGuSS5KsBdbRe8ls1HwVWJfkmiSX\n0Lsxc2DINc0pSYD7gaNV9ZG+WQeALV17C7B/5rqjoKo+WFVrqmotvRuS/1JVtzM+9Z8Gnk5ybdd1\nI/Bt4DOMQf3Ad4E3J7m0+1m6kd4N+nGp/6y5fl5GftxZ1jGzqob+AY7T+8F6rPvs7pv3QXo3N44B\nvznsWuc5hnfQeyP6BLB92PUsUOuv0buW/vW+/+ebgMuBz9N72uAgsHLYtZ7HsbwVONC1x6Z+4PXA\nV4BvAA8Drxqz+nfQe6DgCL2bqi8d5fqBvfTuXzxH7/7dH8xX7yiNO7PU/p7lHDN9WUySGjYSl4Mk\nScNhCEhSwwwBSWqYISBJDTMEJKlhhoAkNcwQkKSGGQKS1LD/B1uSA7X31NmkAAAAAElFTkSuQmCC\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f88b8957c50>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"cbcdiff = cbcdata15[\"likelihood_negbin\"] - cbcdata16[\"likelihood_negbin\"]\n",
"cbcdiff = cbcdiff.replace([np.inf, -np.inf], np.nan)\n",
"cbcdiff = cbcdiff.dropna()\n",
"histogram = plt.hist(cbcdiff)\n",
"\"CBC Likelihood Difference. Min: {}, Max: {}\".format(min(cbcdiff), max(cbcdiff))"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"'CBC Likelihood Difference. Min: 0.0, Max: 0.0'"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAEACAYAAAC6d6FnAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAD35JREFUeJzt3G2MXNddx/HvjzgRDQ21oiA3TyhGctQYlTa1iAstyqCW\nyK1QEorUJNDUUq0qxfRBIKA2SLUFUhWQgBRVyYuSpg5UqawEggshxG2zoqA2biAPThyTWMJV1mCn\npYVSqS8c9c+LuQmT7XpnvDue2d3z/Ugrnzn33Dv/o5m5v5lzZ5yqQpLUnh+adgGSpOkwACSpUQaA\nJDXKAJCkRhkAktQoA0CSGrVgACS5NMnDSZ5O8lSSD3f95yfZn+TZJA8lWTuwz84kzyU5nOSagf5N\nSQ522z5x5qYkSRrFsE8AJ4HfqKqfBN4M/HqSK4AdwP6quhz4YnebJBuBG4CNwBbg9iTpjnUHsK2q\nNgAbkmwZ+2wkSSNbMACq6nhVPd61vws8A1wMXAvs6YbtAa7v2tcB91TVyao6ChwBNie5EDivqg50\n4+4e2EeSNAUjXwNIchlwJfAIsK6qTnSbTgDruvZFwOzAbrP0A2Nu/7GuX5I0JSMFQJJXA/cBH6mq\n/x3cVv3/S8L/T0KSVpg1wwYkOZv+yf8vqur+rvtEktdW1fFueeeFrv8YcOnA7pfQf+d/rGsP9h+b\n574MEklahKrK8FGvNOxbQAHuBA5V1W0Dm/YBW7v2VuD+gf4bk5yTZD2wAThQVceB7yTZ3B3z5oF9\n5k5i1f7t2rVr6jU4v8X/dc/QCf5N7vWw2h+71T6/xRr2CeAtwHuAJ5M81vXtBG4F9ibZBhwF3t09\nWQ8l2QscAl4Ettf/V7cd+AzwKuCBqnpw0VVLkpZswQCoqn/i1J8S3n6KfT4OfHye/n8BXn+6BUqS\nzgx/CTxBvV5v2iWcUat9fqvZan/sVvv8FitLWT8atyS1nOqRBvUvX03y+Zklre+qHUmocV8EliSt\nXgaAJDXKAJCkRhkAktQoA0CSGmUASFKjDABJapQBIEmNMgAkqVEGgCQ1ygCQpEYZAJLUKANAkhpl\nAEhSowwASWqUASBJjTIAJKlRBoAkNcoAkKRGGQCS1CgDQJIaZQBIUqMMAElqlAEgSY0yACSpUQaA\nJDXKAJCkRhkAktQoA0CSGmUASFKjDABJapQBIEmNMgAkqVEGgCQ1ygCQpEYZAJLUKANAkhplAEhS\nowwASWqUASBJjTIAJKlRBoAkNcoAkKRGGQCS1CgDQJIaNTQAknw6yYkkBwf6dieZTfJY9/eOgW07\nkzyX5HCSawb6NyU52G37xPinIkk6HaN8ArgL2DKnr4A/qaoru7+/B0iyEbgB2Njtc3uSdPvcAWyr\nqg3AhiRzjylJmqChAVBVXwa+Pc+mzNN3HXBPVZ2sqqPAEWBzkguB86rqQDfubuD6xZUsSRqHpVwD\n+FCSJ5LcmWRt13cRMDswZha4eJ7+Y12/JGlK1ixyvzuA3+/afwD8MbBtHAXt3r375Xav16PX643j\nsJK0aszMzDAzM7Pk46Sqhg9KLgM+X1WvX2hbkh0AVXVrt+1BYBfwdeDhqrqi678JuLqqPjDnWDVK\nPdI09C9nTfL5GXw9aBRJqKr5luUXtKgloG5N/yW/BLz0DaF9wI1JzkmyHtgAHKiq48B3kmzuLgrf\nDNy/mPuWJI3H0CWgJPcAVwMXJHme/jv6XpI30n879O/ALQBVdSjJXuAQ8CKwfeAt/XbgM8CrgAeq\n6sExz0WSdBpGWgKaFJeAtJy5BKTlaqJLQJKklc8AkKRGGQCS1CgDQJIaZQBIUqMMAElqlAEgSY0y\nACSpUQaAJDXKAJCkRhkAktQoA0CSGmUASFKjDABJapQBIEmNMgAkqVEGgCQ1ygCQpEYZAJLUKANA\nkhplAEhSowwASWqUASBJjTIAJKlRBoAkNcoAkKRGGQCS1CgDQJIaZQBIUqMMAElqlAEgSY0yACSp\nUQaAJDXKAJCkRhkAktQoA0CSGmUASFKjDABJapQBIEmNMgAkqVEGgCQ1ygCQpEYZAJLUKANAkhpl\nAEhSowwASWqUASBJjRoaAEk+neREkoMDfecn2Z/k2SQPJVk7sG1nkueSHE5yzUD/piQHu22fGP9U\nJEmnY5RPAHcBW+b07QD2V9XlwBe72yTZCNwAbOz2uT1Jun3uALZV1QZgQ5K5x5QkTdDQAKiqLwPf\nntN9LbCna+8Bru/a1wH3VNXJqjoKHAE2J7kQOK+qDnTj7h7YR5I0BYu9BrCuqk507RPAuq59ETA7\nMG4WuHie/mNdvyRpStYs9QBVVUlqHMUA7N69++V2r9ej1+uN69CStCrMzMwwMzOz5OOkavi5O8ll\nwOer6vXd7cNAr6qOd8s7D1fV65LsAKiqW7txDwK7gK93Y67o+m8Crq6qD8y5nxqlHmka+pezJvn8\nDL4eNIokVFWGj3ylxS4B7QO2du2twP0D/TcmOSfJemADcKCqjgPfSbK5uyh888A+kqQpGLoElOQe\n4GrggiTPAx8DbgX2JtkGHAXeDVBVh5LsBQ4BLwLbB97Sbwc+A7wKeKCqHhzvVCRJp2OkJaBJcQlI\ny5lLQFquJr0EJEla4QwASWqUASBJjTIAJKlRBoAkNcoAkKRGGQCS1CgDQJIaZQBIUqMMAElqlAEg\nSY0yACSpUQaAJDXKAJCkRhkAktQoA0CSGmUASFKjDABJapQBIEmNMgAkqVEGgCQ1ygCQpEYZAJLU\nKANAkhplAEhSowwASWqUASBJjTIAJKlRBoAkNcoAkKRGGQCS1CgDQJIaZQBIUqMMAElqlAEgSY0y\nACSpUQaAJDXKAJCkRhkAktQoA0CSGmUASFKjDABJapQBIEmNMgAkqVEGgCQ1ygCQpEYZAJLUqCUF\nQJKjSZ5M8liSA13f+Un2J3k2yUNJ1g6M35nkuSSHk1yz1OIlSYu31E8ABfSq6sqquqrr2wHsr6rL\ngS92t0myEbgB2AhsAW5P4icQSZqScZyAM+f2tcCerr0HuL5rXwfcU1Unq+oocAS4CknSVIzjE8AX\nkjya5P1d37qqOtG1TwDruvZFwOzAvrPAxUu8f0nSIq1Z4v5vqar/TPJjwP4khwc3VlUlqQX2/4Ft\nu3fvfrnd6/Xo9XpLLFGSVpeZmRlmZmaWfJxULXR+Po0DJbuA7wLvp39d4HiSC4GHq+p1SXYAVNWt\n3fgHgV1V9cjAMWpc9UjjloR53rOcyXvE14NGkYSqmrscP9Sil4CSnJvkvK79I8A1wEFgH7C1G7YV\nuL9r7wNuTHJOkvXABuDAYu9fkrQ0S1kCWgf8df9dEWuAz1bVQ0keBfYm2QYcBd4NUFWHkuwFDgEv\nAtt9uy9J0zO2JaBxcAlIy5lLQFquJr4EJEla2QwASWqUASBJjTIAJKlRBoAkNcoAkKRGGQCS1CgD\nQJIaZQBIUqMMAElqlAEgSY0yACSpUQaAJDXKAJCkRhkAktQoA0CSGmUASFKjDABJapQBIEmNMgAk\nqVEGgCQ1ygCQpEYZAJLUKANAkhplAEhSowwASWqUASBJjTIAJKlRBoAkNcoAkKRGGQCS1CgDQJIa\nZQBIUqMMAElqlAEgSY0yACSpUQaAJDXKAJCkRhkAktQoA0CSGmUASFKjDABJapQBIEmNMgAkqVEG\ngCQ1ygCQpEZNNACSbElyOMlzST46yfuWJL3SxAIgyVnAJ4EtwEbgpiRXTOr+l4OZmZlpl3BGrfb5\nrWar/bFb7fNbrEl+ArgKOFJVR6vqJPA54LoJ3v/UrfYn4Wqf32q22h+71T6/xZpkAFwMPD9we7br\nkyRNwZoJ3ldN8L5GsnXrVr71rW9N7P42bdo0sfuSpGFSNZnzcpI3A7urakt3eyfw/ar6w4Exyy4k\nJGklqKqc7j6TDIA1wL8BbwP+AzgA3FRVz0ykAEnSK0xsCaiqXkzyQeAfgLOAOz35S9L0TOwTgCRp\neZnqL4GTnJ9kf5JnkzyUZO0pxq1Ncm+SZ5Ic6q4nLHujzq8be1aSx5J8fpI1LsUo80tyaZKHkzyd\n5KkkH55GraMa5ceKSf6s2/5EkisnXeNSDJtfkl/t5vVkkn9O8lPTqHOxRv2xaZKfTvJikndNsr6l\nGvH52evOJU8lmVnwgFU1tT/gj4Df6dofBW49xbg9wPu69hrgNdOse9zz67b/JvBZYN+06x7n/IDX\nAm/s2q+mfx3oimnXfor5nAUcAS4DzgYen1sr8E7gga69GfjqtOse8/x+5qXXF/0fba6q+Q2M+xLw\nt8AvT7vuMT9+a4GngUu62xcsdMxp/19A19I/udP9e/3cAUleA/xcVX0a+tcSqup/JlfikgydH0CS\nS+ifWP4cOO0r+VM0dH5VdbyqHu/a3wWeAS6aWIWnZ5QfK74856p6BFibZN1ky1y0ofOrqq8MvL4e\nAS6ZcI1LMeqPTT8E3At8Y5LFjcEo8/sV4L6qmgWoqm8udMBpB8C6qjrRtU8A872Q1gPfSHJXkn9N\n8qkk506uxCUZZX4Afwr8NvD9iVQ1PqPOD4AklwFX0j+xLEej/FhxvjEr5SR5uj/G3AY8cEYrGq+h\n80tyMf2T5h1d10q6CDrK47cBOL9bdn00yc0LHfCMfwsoyX76ywBz/d7gjaqqU/wOYA3wJuCDVfW1\nJLcBO4CPjb3YRVjq/JL8IvBCVT2WpHdmqly8MTx+Lx3n1fTfdX2k+ySwHI16Mpj7KW2lnERGrjPJ\nzwPvA95y5soZu1Hmdxuwo3u+hpX1iXuU+Z1N/3z5NuBc4CtJvlpVz803+IwHQFX9wqm2JTmR5LVV\ndTzJhcAL8wybBWar6mvd7XvpB8CyMIb5/SxwbZJ3Aj8M/GiSu6vqvWeo5NMyhvmR5GzgPuAvq+r+\nM1TqOBwDLh24fSn9599CYy7p+laCUeZHd+H3U8CWqvr2hGobh1Hmtwn4XP/czwXAO5KcrKp9kylx\nSUaZ3/PAN6vqe8D3kvwj8AZg3gCY9hLQPmBr194K/MDJoaqOA88nubzrejv9ixwrwSjz+92qurSq\n1gM3Al9aLif/EQydX/cu607gUFXdNsHaFuNRYEOSy5KcA9xAf46D9gHvhZd/3f7fA8tgy93Q+SX5\nceCvgPdU1ZEp1LgUQ+dXVT9RVeu719u9wK+tkJM/jPb8/Bvgrd23Cs+l/0WFQ6c84pSvap8PfAF4\nFngIWNv1XwT83cC4NwBfA56g/+RcKd8CGml+A+OvZmV9C2jo/IC30r+28TjwWPe3Zdq1LzCnd9D/\nptIRYGfXdwtwy8CYT3bbnwDeNO2axzk/+l9E+K+Bx+rAtGse9+M3MPYu4F3Trnnc8wN+i/6b5IPA\nhxc6nj8Ek6RGTXsJSJI0JQaAJDXKAJCkRhkAktQoA0CSGmUASFKjDABJapQBIEmN+j8QuFNYW2JW\nIwAAAABJRU5ErkJggg==\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f88b8851e90>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"cbcdiff = cbcdata16[\"likelihood_negbin\"] - cbcdata17[\"likelihood_negbin\"]\n",
"cbcdiff = cbcdiff.replace(np.inf, np.nan)\n",
"cbcdiff = cbcdiff.dropna()\n",
"histogram = plt.hist(cbcdiff)\n",
"\"CBC Likelihood Difference. Min: {}, Max: {}\".format(min(cbcdiff), max(cbcdiff))"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 2",
"language": "python",
"name": "python2"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 2
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython2",
"version": "2.7.11"
}
},
"nbformat": 4,
"nbformat_minor": 0
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment