Skip to content

Instantly share code, notes, and snippets.

@aflaxman

aflaxman/2017_09_07a_smartva_bootstrap_for_csmf_uncertainty.ipynb

Created November 6, 2018 21:33
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save aflaxman/d5c523f7fe86f619b37aac35a3da6738 to your computer and use it in GitHub Desktop.
Save aflaxman/d5c523f7fe86f619b37aac35a3da6738 to your computer and use it in GitHub Desktop.
Download ZIP
Display the source blob
Display the rendered blob
Raw
2017_09_07a_smartva_bootstrap_for_csmf_uncertainty.ipynb
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Mon Sep 11 09:26:37 PDT 2017\r\n"
]
}
],
"source": [
"import numpy as np, pandas as pd, matplotlib.pyplot as plt, seaborn as sns\n",
"%matplotlib inline\n",
"!date"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# SmartVA user asked about CSMF uncertainty\n",
"\n",
"This is something I've always thought worth doing, but never got around to: bootstrap resample the predictions to get a UI around the CSMF:"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"import glob"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"results = []\n",
"\n",
"for module in ['adult', 'child', 'neonate']:\n",
" for fname in glob.glob(f'/homes/abie/projects/2017/smartva_testing/'\\\n",
" f'*/{module}-predictions.csv'):\n",
" df = pd.read_csv(fname)\n",
" results.append({'fname':fname,\n",
" 'module': module,\n",
" 'n': len(df),\n",
" 'undet': df.cause.isnull().sum()})\n",
"results = pd.DataFrame(results)"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"fname\n",
"/homes/abie/projects/2017/smartva_testing/phmrc_database_neonate/neonate-predictions.csv 2621\n",
"/homes/abie/projects/2017/smartva_testing/150817_neontae_only/neonate-predictions.csv 2823\n",
"/homes/abie/projects/2017/smartva_testing/150817_VA_Data_for_India/neonate-predictions.csv 2851\n",
"/homes/abie/projects/2017/smartva_testing/150817_VA_Data_for_India/adult-predictions.csv 6502\n",
"/homes/abie/projects/2017/smartva_testing/phmrc_database_adult/adult-predictions.csv 7839\n",
"Name: n, dtype: int64"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"results.groupby(['fname']).n.first().sort_values().tail()"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"df = pd.read_csv('/homes/abie/projects/2017/smartva_testing/phmrc_database_adult/adult-predictions.csv')"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<style>\n",
" .dataframe thead tr:only-child th {\n",
" text-align: right;\n",
" }\n",
"\n",
" .dataframe thead th {\n",
" text-align: left;\n",
" }\n",
"\n",
" .dataframe tbody tr th {\n",
" vertical-align: top;\n",
" }\n",
"</style>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>sid</th>\n",
" <th>cause</th>\n",
" <th>cause34</th>\n",
" <th>age</th>\n",
" <th>sex</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td>Adult1</td>\n",
" <td>6.0</td>\n",
" <td>Cirrhosis</td>\n",
" <td>51.0</td>\n",
" <td>1</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td>Adult2</td>\n",
" <td>22.0</td>\n",
" <td>Other Cardiovascular Diseases</td>\n",
" <td>24.0</td>\n",
" <td>1</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2</th>\n",
" <td>Adult3</td>\n",
" <td>9.0</td>\n",
" <td>Diabetes</td>\n",
" <td>62.0</td>\n",
" <td>2</td>\n",
" </tr>\n",
" <tr>\n",
" <th>3</th>\n",
" <td>Adult4</td>\n",
" <td>25.0</td>\n",
" <td>Other Non-communicable Diseases</td>\n",
" <td>80.0</td>\n",
" <td>2</td>\n",
" </tr>\n",
" <tr>\n",
" <th>4</th>\n",
" <td>Adult5</td>\n",
" <td>36.0</td>\n",
" <td>Ischemic Heart Disease</td>\n",
" <td>76.0</td>\n",
" <td>1</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" sid cause cause34 age sex\n",
"0 Adult1 6.0 Cirrhosis 51.0 1\n",
"1 Adult2 22.0 Other Cardiovascular Diseases 24.0 1\n",
"2 Adult3 9.0 Diabetes 62.0 2\n",
"3 Adult4 25.0 Other Non-communicable Diseases 80.0 2\n",
"4 Adult5 36.0 Ischemic Heart Disease 76.0 1"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"df.head()"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"0 Maternal\n",
"1 Stroke\n",
"2 Stroke\n",
"3 Maternal\n",
"4 Undetermined\n",
"Name: cause34, dtype: object"
]
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"def bootstrap_resample(s):\n",
" \"\"\" Resample pd.Series s *with* replacement\n",
" \"\"\"\n",
" \n",
" vals = np.random.choice(s.values, size=len(s), replace=True)\n",
" return pd.Series(vals, index=s.index, name=s.name)\n",
"\n",
"bootstrap_resample(df.cause34).head()"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"Maternal 0.050644\n",
"AIDS 0.054599\n",
"Undetermined 0.055237\n",
"Diabetes 0.057788\n",
"Stroke 0.072458\n",
"Name: cause34, dtype: float64"
]
},
"execution_count": 11,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"def csmf(s):\n",
" \"\"\" Find CSMF for pd.Series of cause-of-deaths\n",
" *without* redistribution\n",
" \"\"\"\n",
" return s.value_counts(normalize=True)\n",
"csmf(df.cause34).sort_values().tail()"
]
},
{
"cell_type": "code",
"execution_count": 20,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<style>\n",
" .dataframe thead tr:only-child th {\n",
" text-align: right;\n",
" }\n",
"\n",
" .dataframe thead th {\n",
" text-align: left;\n",
" }\n",
"\n",
" .dataframe tbody tr th {\n",
" vertical-align: top;\n",
" }\n",
"</style>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>lb</th>\n",
" <th>ub</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>Maternal</th>\n",
" <td>0.041873</td>\n",
" <td>0.044998</td>\n",
" </tr>\n",
" <tr>\n",
" <th>AIDS</th>\n",
" <td>0.047136</td>\n",
" <td>0.048953</td>\n",
" </tr>\n",
" <tr>\n",
" <th>Undetermined</th>\n",
" <td>0.047391</td>\n",
" <td>0.049241</td>\n",
" </tr>\n",
" <tr>\n",
" <th>Diabetes</th>\n",
" <td>0.049273</td>\n",
" <td>0.051537</td>\n",
" </tr>\n",
" <tr>\n",
" <th>Stroke</th>\n",
" <td>0.062348</td>\n",
" <td>0.065536</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" lb ub\n",
"Maternal 0.041873 0.044998\n",
"AIDS 0.047136 0.048953\n",
"Undetermined 0.047391 0.049241\n",
"Diabetes 0.049273 0.051537\n",
"Stroke 0.062348 0.065536"
]
},
"execution_count": 20,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"def csmf_ui(s):\n",
" \"\"\" Use `n_reps` bootstrap resamples to find the 95% UI\n",
" for CSMF values from cause-of-death pd.Series\"\"\"\n",
" n_reps = 1000\n",
" \n",
" df = pd.DataFrame(0, columns=range(n_reps), index=s.unique())\n",
" for rep in range(n_reps):\n",
" df[rep] = csmf(bootstrap_resample(s))\n",
" \n",
" df = df.fillna(0)\n",
" lb = np.percentile(df, .025, axis=1)\n",
" ub = np.percentile(df, .975, axis=1)\n",
" return pd.DataFrame({'lb': lb, 'ub':ub}, index=df.index)\n",
"\n",
"csmf_ui(df.cause34).sort_values('ub').tail()"
]
},
{
"cell_type": "code",
"execution_count": 21,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"df = pd.read_csv('/homes/abie/projects/2017/smartva_testing/phmrc_database_neonate/neonate-predictions.csv')"
]
},
{
"cell_type": "code",
"execution_count": 25,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<style>\n",
" .dataframe thead tr:only-child th {\n",
" text-align: right;\n",
" }\n",
"\n",
" .dataframe thead th {\n",
" text-align: left;\n",
" }\n",
"\n",
" .dataframe tbody tr th {\n",
" vertical-align: top;\n",
" }\n",
"</style>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>lb</th>\n",
" <th>ub</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>Undetermined</th>\n",
" <td>0.081744</td>\n",
" <td>0.084319</td>\n",
" </tr>\n",
" <tr>\n",
" <th>Meningitis/Sepsis</th>\n",
" <td>0.089755</td>\n",
" <td>0.093857</td>\n",
" </tr>\n",
" <tr>\n",
" <th>Birth asphyxia</th>\n",
" <td>0.090423</td>\n",
" <td>0.094140</td>\n",
" </tr>\n",
" <tr>\n",
" <th>Preterm Delivery</th>\n",
" <td>0.129530</td>\n",
" <td>0.138779</td>\n",
" </tr>\n",
" <tr>\n",
" <th>Stillbirth</th>\n",
" <td>0.359308</td>\n",
" <td>0.371797</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" lb ub\n",
"Undetermined 0.081744 0.084319\n",
"Meningitis/Sepsis 0.089755 0.093857\n",
"Birth asphyxia 0.090423 0.094140\n",
"Preterm Delivery 0.129530 0.138779\n",
"Stillbirth 0.359308 0.371797"
]
},
"execution_count": 25,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"csmf_ui(df.cause34).sort_values('ub').tail()"
]
},
{
"cell_type": "code",
"execution_count": 26,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<style>\n",
" .dataframe thead tr:only-child th {\n",
" text-align: right;\n",
" }\n",
"\n",
" .dataframe thead th {\n",
" text-align: left;\n",
" }\n",
"\n",
" .dataframe tbody tr th {\n",
" vertical-align: top;\n",
" }\n",
"</style>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>lb</th>\n",
" <th>ub</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>Meningitis/Sepsis</th>\n",
" <td>0.067999</td>\n",
" <td>0.076</td>\n",
" </tr>\n",
" <tr>\n",
" <th>Birth asphyxia</th>\n",
" <td>0.074250</td>\n",
" <td>0.082</td>\n",
" </tr>\n",
" <tr>\n",
" <th>Undetermined</th>\n",
" <td>0.073249</td>\n",
" <td>0.083</td>\n",
" </tr>\n",
" <tr>\n",
" <th>Preterm Delivery</th>\n",
" <td>0.126000</td>\n",
" <td>0.131</td>\n",
" </tr>\n",
" <tr>\n",
" <th>Stillbirth</th>\n",
" <td>0.355250</td>\n",
" <td>0.367</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" lb ub\n",
"Meningitis/Sepsis 0.067999 0.076\n",
"Birth asphyxia 0.074250 0.082\n",
"Undetermined 0.073249 0.083\n",
"Preterm Delivery 0.126000 0.131\n",
"Stillbirth 0.355250 0.367"
]
},
"execution_count": 26,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"csmf_ui(df.cause34.iloc[:1000]).sort_values('ub').tail()"
]
},
{
"cell_type": "code",
"execution_count": 27,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<style>\n",
" .dataframe thead tr:only-child th {\n",
" text-align: right;\n",
" }\n",
"\n",
" .dataframe thead th {\n",
" text-align: left;\n",
" }\n",
"\n",
" .dataframe tbody tr th {\n",
" vertical-align: top;\n",
" }\n",
"</style>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>lb</th>\n",
" <th>ub</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>Undetermined</th>\n",
" <td>0.022498</td>\n",
" <td>0.04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>Birth asphyxia</th>\n",
" <td>0.022498</td>\n",
" <td>0.04</td>\n",
" </tr>\n",
" <tr>\n",
" <th>Meningitis/Sepsis</th>\n",
" <td>0.030000</td>\n",
" <td>0.05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>Preterm Delivery</th>\n",
" <td>0.040000</td>\n",
" <td>0.06</td>\n",
" </tr>\n",
" <tr>\n",
" <th>Stillbirth</th>\n",
" <td>0.270000</td>\n",
" <td>0.32</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" lb ub\n",
"Undetermined 0.022498 0.04\n",
"Birth asphyxia 0.022498 0.04\n",
"Meningitis/Sepsis 0.030000 0.05\n",
"Preterm Delivery 0.040000 0.06\n",
"Stillbirth 0.270000 0.32"
]
},
"execution_count": 27,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"csmf_ui(df.cause34.iloc[:100]).sort_values('ub').tail()"
]
},
{
"cell_type": "code",
"execution_count": 42,
"metadata": {},
"outputs": [],
"source": [
"def csmf_hbar_w_error(df):\n",
" t = csmf_ui(df.cause34)\n",
" t['csmf'] = csmf(df.cause34)\n",
" \n",
" t *= 100\n",
" t = t.sort_values('csmf')\n",
" \n",
" t.csmf.plot(kind='barh', xerr=[(t.csmf-t.lb), (t.ub-t.csmf)])\n",
" return t\n"
]
},
{
"cell_type": "code",
"execution_count": 43,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAeYAAAD8CAYAAACiqQeGAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmYXHWB7vHvS4AQCCbjJCNNQNslyshiIBFlNTjIVWRE\nBxwU0ABzzTDjIOjlalyugmsUnxFBHScquBAwGhcYQJAlDQjEkEDIAkQEwkAA2YZAWIIk7/2jfg1F\n20t1p7vrdPf7eZ48deqc31qtvPU751SVbBMRERHVsFmzBxAREREvSDBHRERUSII5IiKiQhLMERER\nFZJgjoiIqJAEc0RERIUkmCMiIiokwRwREVEhCeaIiIgK2bzZA4ihYcKECW5tbW32MCIihpQlS5Y8\nbHtib+okmKMhra2tLF68uNnDiIgYUiTd3ds6OZUdERFRIQnmiIiICkkwR0REVEiCOSIiokISzBER\nERWSYI6IiKiQfFwqGrJ8zVpaZ13U7GFEjFirZ7+z2UOIQZIVc0RERIUkmCMiIiokwRwREVEhCeZB\nJOnTklZKWiZpqaQ3STpJ0tZ1ZS6WNL5sryuPrZJWlO1jJH2ri/av62L/eEn/Wvd8uqQL+3NuERHR\nPxLMg0TSXsAhwB62dwMOBO4BTgKeD2bbB9t+rC992N67k343B8YD//qXNSIiomoSzIOnBXjY9noA\n2w8DhwPbAwskLQCQtFrShB7a2lFSm6TbJX2ufWfdCnu6pGskXQDcAswGXl1W6aeV4mMlzZd0m6S5\nktSvs42IiD7Jx6UGz2+Bz0r6A3A5MM/2GZI+BhxQgrpRewK7AE8BN0i6yHbHn37aA9jF9l2SWsv2\nFKgFN7A7sDNwH3AtsA/wu75OLmIkeeDcWYPe5/SFp/VcqJ+1tbUNep+RFfOgsb0OmArMBB4C5kk6\npo/NXWb7EdtPA78E9u2kzCLbd3XTxiLb99reCCwFWjsWkDRT0mJJizc8tbaPQ42IiN7IinkQ2d4A\ntAFtkpYDM/raVA/PAZ7soY31ddsb6OR/C7bnAHMARrdM7qyPiBFpuyNnD3qfbfmCkREjK+ZBIul1\nkibX7ZoC3A08AWzby+beJumlksYA76Z2Kro7fekjIiKaICvmwTMWOLN8FOo54I/UTmu/H7hE0n22\nD2iwrUXAL4AdgHM6ub78IrYfkXRt+cjVb4B8t2ZEREXJzhnK6NnolslumXF6s4cRMWLlu7KHJklL\nbE/rTZ2cyo6IiKiQBHNERESFJJgjIiIqJDd/RUN2nTSOxbnGFREx4LJijoiIqJAEc0RERIUkmCMi\nIiokwRwREVEhCeaIiIgKSTBHRERUSII5IiKiQhLMERERFZJgjoiIqJAEc0RERIUkmCMiIiokwRwR\nEVEhCeaIiIgKSTBHRERUSH72MRqyfM1aWmdd1OxhDFur85OaEVFkxRwREVEhCeaIiIgKSTBHRERU\nyJAPZkkbJC2VtELSzyVt3cv6J/W2Tl9JWi1pefl3i6QvStqqgXrryuP2kuYP/EgjIqJZhnwwA0/b\nnmJ7F+BZ4Pj6g6rpbp4nAb0N81G9H+bzDrC9K7An8CrgPxutaPs+24dvQt8ASMpNfxERFTUcgrne\nNcBrJLVKWiXpx8AKYEdJB0m6XtKNZWU9VtJHgO2BBZIWAHRWruxfLemrkm4E3iupTdI3JC2WdKuk\nN0r6paTbJX2xp4HaXkftTcS7Jb209PF/Jd0gaZmkUzvWKfNaUbYXStq57libpGmStpF0lqRFkm6S\ndGg5foykCyRdCVwh6ceS3l1Xf2572YiIaJ5hs3Iqq8B3AJeUXZOBGbYXSpoAfAY40PaTkj4BfMz2\n5yV9jNoq9uGuygGfL20+YnuP0t/xwLO2p0k6ETgfmAo8Ctwh6Ru2H+luzLYfl3QXMFnSuDLmPQEB\nF0ja3/bVXVSfB/wj8DlJLUCL7cWSvgxcafs4SeOBRZIuL3X2AHaz/aiktwAfBX5d+t4bmNHASx3F\nA+fO6re2pi88rd/aamtr67e2ImLwDYdgHiNpadm+BvgBtVXw3bYXlv1vBl4PXCsJYEvg+k7a6qnc\nvA7lLyiPy4GVtu8HkHQnsCPQbTAXKo8HlX83ledjqQV1V8H8M+C3wOeoBXT7teeDgHdJOrk83wp4\nedm+zPajALavkvQdSROBw4Bf2H7uRQOTZgIzAUa9ZGIDU4mIiE01HIL5adtT6neUUH2yfhe1UHp/\nD231VO7JDs/Xl8eNddvtz3t8bSVtC7QCfyh9f8V2Q9ecba+R9Iik3YAjeOHauoDDbK/q0NebOhn/\nj4GjgfcBx3bSxxxgDsDolsluZFwjyXZHzu63ttryBSMRUQy3a8xdWQjsI+k1AOU67GvLsSeAbRso\n16/KtevvAL+2/T/ApcBxdde0J0n6mx6amQd8HBhne1nZdylwgsq7E0m7d1P/h9RufsP2LX2dS0RE\n9J8REcy2HwKOAc6TtIza6emdyuE5wCWSFvRQrr8sKDdwLQL+G/jnMsbfAucC10taTu3U9LZdtlIz\nn9pq92d1+74AbAEsk7SyPO+U7T8BtwJn920qERHR32TnDOVIVT6/vRzYw/ba7sqObpnslhmnD87A\nRqB8V3bE8CRpie1pvakzIlbM8ZckHUhttXxmT6EcERGDZzjc/BV9YPty4BXNHkdERLxYVswREREV\nkhVzNGTXSeNYnOugEREDLivmiIiICkkwR0REVEiCOSIiokISzBERERWSYI6IiKiQBHNERESFJJgj\nIiIqJMEcERFRIQnmiIiICkkwR0REVEiCOSIiokISzBERERWSYI6IiKiQBHNERESF5GcfoyHL16yl\nddZFzR7GgFudn7aMiCbLijkiIqJCEswREREVkmCOiIiokBEXzJI2SFoq6WZJN0rau+zfXtL8Luq0\nSjqy7vkxkr41SOOdLunCfmrruv5oJyIiBs6IC2bgadtTbL8B+CTwFQDb99k+vGNhSZsDrcCRHY8N\nNbb3bvYYIiKieyMxmOu9BPgfeH5VvKJsHyPpAklXAlcAs4H9ykr7o6Xu9pIukXS7pK911rikz0q6\nQdIKSXMkqez/iKRbJC2T9NOy7xRJP5F0fWnzQ3VNjZU0X9Jtkuaq5q2Sfl3X19sk/UrSK0r9CZI2\nk3SNpINKmXXlcaykK8oZg+WSDu3PFzUiIvpuJH5caoykpcBWQAvw1i7K7QHsZvtRSdOBk20fArXg\nBqYAuwPrgVWSzrR9T4c2vmX786XOT4BDgP8CZgGvtL1e0vi68rsBbwa2AW6S1P75pN2BnYH7gGuB\nfYAFwHckTbT9EHAscJbtuyV9FfgPYBFwi+3fdhjXM8B7bD8uaQKwUNIFtt3zy1dND5w7q1/amb7w\ntH5pp62trV/aiYiRZySumNtPZe8EvB34cftKtoPLbD/aTTtX2F5r+xngFuAVnZQ5QNLvJS2n9gZg\n57J/GTBX0tHAc3Xlz7f9tO2HqQXvnmX/Itv32t4ILAVaS4j+BDi6hPtewG8AbH+f2tmA44GTOxmX\ngC9LWgZcDkwCXvYXhaSZkhZLWrzhqbXdvBQREdFfRuKK+Xm2ry8rxomdHH6yh+rr67Y30OG1lLQV\n8B1gmu17JJ1CbZUO8E5gf+DvgU9L2rV9SB2H2ENfZ1NbgT8D/Nz2c6XvrYEdSpmxwBMd2j2K2pyn\n2v6zpNV1Y3uhc3sOMAdgdMvkSq+mtztydr+005YvGImIJhuJK+bnSdoJGAU80kPRJ4Bte9l8e9A9\nLGkscHjpczNgR9sLgE8A46iFJ8ChkraS9NfAdOCG7jqwfR+109ufoRbS7b4KzAU+C3yvk6rjgAdL\nKB9A56v9iIhogpG4Ym6/xgy1U7ozbG/o/Gz285YBGyTdDPyQcsNYd2w/Jul7wArgAV4I2VHAOZLG\nlf7PKGXb+1kATAC+YPs+Sa/toau5wETbtwJIegvwRmCfMq/DJB1r++wOdf6rnGJfDNzW03wiImJw\naAjf7zOslFPd62x/vZf1vgXcZPsHAzKwYnTLZLfMOH0gu6iEfFd2RPQnSUtsT+tNnZG4Yh42JC2h\ndi38/zR7LBER0T8SzBVh+5Q+1Jk6AEOJiIgmGtE3f0VERFRNVszRkF0njWNxrr9GRAy4rJgjIiIq\nJMEcERFRIQnmiIiICkkwR0REVEiCOSIiokISzBERERWSYI6IiKiQBHNERESFJJgjIiIqJMEcERFR\nIQnmiIiICkkwR0REVEiCOSIiokISzBERERWSn32Mhixfs5bWWRc1exgDbnV+2jIimiwr5oiIiApJ\nMEdERFRIgjkiIqJCRmQwS2qVtKLDvlMkndyLNtokTeuhzEmStu7rOBsYw3X91M5fvB4REdEcIzKY\nB9FJQK+CWdKoRsva3rvXI4qIiEpLMHdQVsJflbRI0h8k7Vf2j5H0U0m3SvoVMKauzkGSrpd0o6Sf\nSxor6SPA9sACSQu6Klf2ry593gi8t4zhG5IWl/7eKOmXkm6X9MW6fteVx+mlznxJt0maK0nl2FRJ\nV0laIulSSS11+2+WdDPw4UF5cSMiokf5uFTnNre9p6SDgc8BBwL/Ajxl+28l7QbcCCBpAvAZ4EDb\nT0r6BPAx25+X9DHgANsPd1UO+Hzp8xHbe5Q2jweetT1N0onA+cBU4FHgDknfsP1IhzHvDuwM3Adc\nC+wj6ffAmcChth+SdATwJeA44Gzg32xfLem0/n4Bm+WBc2dtUv3pC/v+UrS1tW1S3xERMHKD2T3s\n/2V5XAK0lu39gTMAbC+TtKzsfzPweuDaskjdEri+k7Z7KjevQ/kLyuNyYKXt+wEk3QnsCHQM5kW2\n7y1llpZxPwbsAlxW+hwF3C9pPDDe9tWl7k+Ad3QcsKSZwEyAUS+Z2MmUIiKiv43UYH4E+KsO+14K\n3FW215fHDfT8Ggm4zPb7N7Hckx2et49hY912+/POxlRfpn3cohbqe71oILVg7pHtOcAcgNEtk7t6\nM1Mp2x05e5Pqt+ULRiKiyUbkNWbb66itHN8KIOmlwNuB33VT7WrgyFJ+F2C3sn8htdPGrynHtpH0\n2nLsCWDbBsoNlFXAREl7lT63kLSz7ceAxyTtW8odNcDjiIiIBo3IYC4+CPy/ctr3SuBU23d0U/4/\ngLGSbqV2XXgJgO2HgGOA88rp7euBnUqdOcAlkhb0UG5A2H4WOBz4arnJaynQfif3scC3y/w1kOOI\niIjGyR4SZyijyUa3THbLjNObPYwBl+/Kjoj+JGmJ7W6/86KjkbxijoiIqJwEc0RERIUkmCMiIipk\npH5cKnpp10njWJzrrxERAy4r5oiIiApJMEdERFRIgjkiIqJCEswREREVkmCOiIiokARzREREhSSY\nIyIiKiTBHBERUSEJ5oiIiApJMEdERFRIgjkiIqJCEswREREVkmCOiIiokARzREREheRnH6Mhy9es\npXXWRc0exoBYnZ+zjIgKyYo5IiKiQhLMERERFZJgjoiIqJAhGcySLOmcuuebS3pI0oWb0ObFksb3\nse40SWeU7emS9q47drykD/ZQfwtJN5btT0taKWmZpKWS3tSXMXXRz/aS5vdXexER0f+G6s1fTwK7\nSBpj+2ngbcCaTWnQ9sGbUHcxsLg8nQ6sA64rx77bQBP7AtdK2gs4BNjD9npJE4At+zquTsZ5H3B4\nf7UXERH9b0iumIuLgfbbad8PnNd+QNI2ks6StEjSTZIOLfuPkfRLSZdIul3S1+rqrJY0QVKrpFsl\nfa+sXH8raUwp88a6lexpklaU/dMlXSipFTge+Ggps5+kUySdXMp9RNItpY2f1s3l7cBvgBbgYdvr\nAWw/XMIUSVMlXSVpiaRLJbWU/W2Svln6WyFpz7L/LWXf0vIabFvm1j7mncvrs7SMZ3K//nUiIqJP\nhuqKGeCnwGfL6evdgLOA/cqxTwNX2j6unJ5eJOnycmwKsDuwHlgl6Uzb93RoezLwftsfkvQz4DDg\nHOBs4EO2r5c0u+OAbK+W9F1gne2vA0j6u7ois4BXltVw/WnzA4BTqb1R+qykPwCXA/NsXyVpC+BM\n4FDbD0k6AvgScFypv7XtKZL2L6/DLsDJwIdtXytpLPBMh+EeD3zT9lxJWwKjunidh4wHzp3Vp3rT\nF57W5z7b2tr6XDciojNDdsVsexnQSm21fHGHwwcBsyQtBdqArYCXl2NX2F5r+xngFuAVnTR/l+2l\nZXsJ0FqCdFvb15f95/Zh2MuAuZKOBp4DkDQJeNT2U7bXAVOBmcBDwDxJxwCvoxa2l5U5fQbYoa7d\n8wBsXw28pIz1WuDfJX0EGG/7uQ5juR74lKRPAK8olwReRNJMSYslLd7w1No+TDciInprKK+YAS4A\nvk7tuu5f1+0XcJjtVfWFy41U6+t2baDz16BjmTH9MVhqp973B/4e+LSkXamdxr60vYDtDdTeTLRJ\nWg7MoPbmYKXtvbpo1x2f254t6SLgYGrXr/8Xdatm2+dK+n0Z08WS/tn2lR0amQPMARjdMrljH5Wz\n3ZF/cRKjIW35gpGIqJAhu2IuzgJOtb28w/5LgRMkCUDS7pvake3HgCfq7pJ+XxdFnwC27bhT0mbA\njrYXAJ8AxgFjeeH6MpJe1+Fa7xTgbmAVMLHcHNZ+F/fOdeWOKPv3BdbaXivp1baX2/4qcAOwU4fx\nvAq40/YZwPnULgdERESTDekVs+17gTM6OfQF4HRgWQnEu6jd7byp/gn4nqSNwFVAZ+d3/wuYX244\nO6Fu/yjgHEnjqK3oz6AW4q+xfVspMxY4s5yKfg74IzDT9rOSDgfOKPU3L/NbWeo9I+kmYAteuO58\nkqQDgI2lXPvNZe3+EfiApD8DDwBf7tMrEhER/Up25c9QVoakseU6MJJmAS22T9yE9vYFjrZ9/Ca0\n0QacXD6yNWBGt0x2y4zTB7KLpsl3ZUfEQJG0xPa03tQZ0ivmJninpE9Se93uBo7ZlMZs/w74XT+M\nKyIihokEcy/YngfMa/Y46tme3uwxRERE/xnqN39FREQMK1kxR0N2nTSOxbkWGxEx4LJijoiIqJAE\nc0RERIUkmCMiIiokwRwREVEhCeaIiIgKSTBHRERUSII5IiKiQhLMERERFZJgjoiIqJAEc0RERIUk\nmCMiIiokwRwREVEhCeaIiIgKSTBHRERUSH72MRqyfM1aWmdd1Oxh9Nrq/FRlRAwxWTFHRERUSII5\nIiKiQhLMERERFZJgjoiIqJAEcx9I2iBpqaQVkn4uaetmj6kRkraXNL/Z44iIiK4lmPvmadtTbO8C\nPAsc3+wBNcL2fbYPb/Y4IiKia/m41Ka7BthNUivwG+B3wN7AGuBQ209LejXwbWAi8BTwIdu3Sfoh\ncKHt+QCS1tkeK2k6cCrwGLAr8DNgOXAiMAZ4t+07Sp9nAROAh4Bjbf93afdxYBqwHfBx2/NL+Qtt\n71K2fwJsU+bxb7av6/+Xp388cO6sPtWbvvC0Xtdpa2vrU18REf0hK+ZNIGlz4B3UQhNgMvBt2ztT\nC9XDyv45wAm2pwInA99poPk3UFuJ/y3wAeC1tvcEvg+cUMqcCfzI9m7AXOCMuvotwL7AIcDsTtp/\nEHib7T2AIzrUbZ/fTEmLJS3e8NTaBoYcERGbKivmvhkjaWnZvgb4AbA9cJft9v1LgFZJY6mtoH8u\nqb3+6Ab6uMH2/QCS7gB+W/YvBw4o23sB/1C2fwJ8ra7+r21vBG6R9LJO2t8C+JakKcAG4LUdC9ie\nQ+1NBaNbJruBMQ+Y7Y7s7L1Fz9ryBSMRMcQkmPvmadtT6neU0F1ft2sDtdPOmwGPdSxfPFeOI2kz\nYMu6Y/Vtbax7vpHG/m719dXJ8Y8Cf6K2Mt8MeKaBNiMiYoDlVPYAs/04cJek9wKo5g3l8Gpgatl+\nF7VVbG9cB7yvbB9FbfXeqHHA/WVV/QFgVC/7joiIAZBgHhxHAf8k6WZgJXBo2f894C1l/17Ak71s\n9wTgWEnLqIXrib2o+x1gRul7pz70HRERA0B2Uy8dxhAxumWyW2ac3uxh9Fp+xCIimknSEtvTelMn\nK+aIiIgKSTBHRERUSO7KjobsOmkci3NaOCJiwGXFHBERUSEJ5oiIiApJMEdERFRIgjkiIqJCEswR\nEREVkmCOiIiokARzREREhSSYIyIiKiTBHBERUSEJ5oiIiApJMEdERFRIgjkiIqJCEswREREVkmCO\niIiokARzREREheT3mKMhy9espXXWRU3pe3V+BzoiRpCsmCMiIiokwRwREVEhCeaIiIgK6TGYJW0n\n6aeS7pC0RNLFkl47GIPrYVwXSxpf/v1rg3XWDcA4pku6sGyPlnS5pKWSjujvvroZw6c6PL9usPqO\niIj+1W0wSxLwK6DN9qttTwU+CbxsMAbXHdsH234MGA80FMyDYHcA21Nsz2ukgqRR/dDvi4LZ9t79\n0GZERDRBT3dlHwD82fZ323fYvhmeD+2vAe8ADHzR9jxJ04FTgIeBXYAlwNG2Lelg4N+BJ4FrgVfZ\nPkTSNsCZpfwWwCm2z5d0DPAuYGvg1cCvbH+89L8amAbMBl4taSlwGXAqcD7wV6Wtz9g+v6sJSmoF\nLgEWAnsDNwBnl3b+BjjK9iJJewLfBLYCngaOtb2qrp2/Ac4BJpaxHAa0Al8vr/MNwL/YXl/GPg94\nG/A1SccDNwH7AdsAH6T2BmhXYJ7tz5Q+fg3sWMbwTdtzJM0GxpQ+V9o+StI622P78jfq6nXqbw+c\nO6vhstMXntarttva2no5moiI6ugpmNv/o92ZfwCmAG8AJgA3SLq6HNsd2Bm4j1oA7yNpMfCfwP62\n75J0Xl1bnwautH2cpPHAIkmXl2NTSnvrgVWSzrR9T13dWcAutqcASNoceI/txyVNABZKuqCH0HkN\n8F7gOGoBeiSwL7U3BZ8C3g3cBuxn+zlJBwJfpha+ANh+UNL/Bk4ubza2AtqAv7P9B0k/Bv4FOL1U\necT2HmXMxwPP2p4m6URqbyymAo8Cd0j6hu1HgONsPyppTHm9f2F7lqR/a5//pvyNgN/VV5Y0E5gJ\nMOolE7t5+SIior9syueY9wXOs70B+JOkq4A3Ao8Di2zfC1BWcq3AOuBO23eV+udR/qMPHAS8S9LJ\n5flWwMvL9hW215a2bgFeAdQHc0cCvixpf2AjMInaqfcHuqlzl+3lpY+VpU9LWl7GDjAO+JGkydRW\nn1t00x7A60q7fyjPfwR8mBeCueOp7gvK43JqK9/7y3jupLZKfgT4iKT3lHI7ApPL/q709m/0omC2\nPQeYAzC6ZXK/rqa3O3J2w2Xb8jnmiBhBegrmlcDhfWh3fd32hgb6EXBY/alhAElv6kNbRwETgam2\n/1xOG2/Vi/FurHu+sa6/LwALbL+nnP5u66HNnjzZxRjq+39+DOX084HAXrafktRGz/PqTm9f14iI\nGAQ93ZV9JTC6nNIEQNJukvYDrgGOkDRK0kRgf2BRN22tAl5VQg2g/q7lS4ETyjVRJO3eizk8AWxb\n93wc8GAJ5QOorbD7wzhgTdk+poHyq4BWSa8pzz8AXLWJ/f9PCeWdgDfXHfuzpM5W8L39G0VERJN1\nG8zluux7gAPLx6VWAl+hdlr4V8Ay4GZqAf5x212eLrb9NLW7py+RtIRaoK4th79A7dTwstLHFxqd\nQLn2eq2kFZJOA+YC08pp6A9SuzbcH74GfEXSTTSwurT9DHAs8PMylo3Ad7uv1a1LqK2cb6V2w9vC\numNzqL12czvU6dXfKCIimk+DeCMuksbaXldWxt8Gbrf9jUEbQPTZ6JbJbplxes8FB0C+KzsihipJ\nS2xP602dwf7mrw+1f6yH2qnZ/xzk/iMiIiptUG/4KavjrJAjIiK6kDtxoyG7ThrH4pxSjogYcPkR\ni4iIiApJMEdERFRIgjkiIqJCEswREREVkmCOiIiokARzREREhSSYIyIiKmRQv5Izhi5JT1D7YY7h\nagLwcLMHMYAyv6FtOM9vOM8N4HW2t+252AvyBSPRqFW9/b7XoUTS4sxv6Mr8hq7hPDeoza+3dXIq\nOyIiokISzBERERWSYI5GzWn2AAZY5je0ZX5D13CeG/Rhfrn5KyIiokKyYo6IiKiQBHP0SNLbJa2S\n9EdJs5o9nk0l6SxJD0paUbfvpZIuk3R7efyrZo6xryTtKGmBpFskrZR0Ytk/XOa3laRFkm4u8zu1\n7B8W82snaZSkmyRdWJ4Pm/lJWi1puaSl7XcsD7P5jZc0X9Jtkm6VtFdv55dgjm5JGgV8G3gH8Hrg\n/ZJe39xRbbIfAm/vsG8WcIXtycAV5flQ9Bzwf2y/Hngz8OHy9xou81sPvNX2G4ApwNslvZnhM792\nJwK31j0fbvM7wPaUuo9JDaf5fRO4xPZOwBuo/R17Nb8Ec/RkT+CPtu+0/SzwU+DQJo9pk9i+Gni0\nw+5DgR+V7R8B7x7UQfUT2/fbvrFsP0HtPwqTGD7zs+115ekW5Z8ZJvMDkLQD8E7g+3W7h838ujAs\n5idpHLA/8AMA28/afoxezi/BHD2ZBNxT9/zesm+4eZnt+8v2A8DLmjmY/iCpFdgd+D3DaH7lNO9S\n4EHgMtvDan7A6cDHgY11+4bT/AxcLmmJpJll33CZ3yuBh4Czy6WI70vahl7OL8Ec0YFrH1UY0h9X\nkDQW+AVwku3H648N9fnZ3mB7CrADsKekXTocH7Lzk3QI8KDtJV2VGcrzK/Ytf793ULvUsn/9wSE+\nv82BPYD/sL078CQdTls3Mr8Ec/RkDbBj3fMdyr7h5k+SWgDK44NNHk+fSdqCWijPtf3LsnvYzK9d\nOUW4gNr9AsNlfvsA75K0mtplo7dKOofhMz9srymPDwK/ona5bLjM717g3nIWB2A+taDu1fwSzNGT\nG4DJkl4paUvgfcAFTR7TQLgAmFG2ZwDnN3EsfSZJ1K5v3Wr73+sODZf5TZQ0vmyPAd4G3MYwmZ/t\nT9rewXYrtf+vXWn7aIbJ/CRtI2nb9m3gIGAFw2R+th8A7pH0urLr74Bb6OX88gUj0SNJB1O77jUK\nOMv2l5o8pE0i6TxgOrVftfkT8Dng18DPgJcDdwP/aLvjDWKVJ2lf4BpgOS9co/wUtevMw2F+u1G7\neWYUtYUhxtz4AAAAd0lEQVTFz2x/XtJfMwzmV0/SdOBk24cMl/lJehW1VTLUTvuea/tLw2V+AJKm\nULtxb0vgTuBYyv9WaXB+CeaIiIgKyansiIiICkkwR0REVEiCOSIiokISzBERERWSYI6IiKiQBHNE\nRESFJJgjIiIqJMEcERFRIf8flMFO5OQ2Q/gAAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7ff22a988f60>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"csmf_hbar_w_error(df.iloc[:1000])\n",
"plt.axis(xmax=60);"
]
},
{
"cell_type": "code",
"execution_count": 44,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAeYAAAD8CAYAAACiqQeGAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmcXFWB9vHfQ4AQCCbjJCNNQNslyshiIBFlNSjyKjKi\nAw4KaIB5zTDjIOjwalxeBdcofkYE16jgQsBoXGAAQZY0IBBDB0IWICIQBgIIhCEQliDJM3/UaSja\n3re63Xm+n08+devcc849p1p56tx7q0q2iYiIiGrYrNEDiIiIiOclmCMiIiokwRwREVEhCeaIiIgK\nSTBHRERUSII5IiKiQhLMERERFZJgjoiIqJAEc0RERIVs3ugBxPAwYcIENzc3N3oYERHDyuLFix+2\nPbE3bRLM0SPNzc20trY2ehgREcOKpLt72yansiMiIiokwRwREVEhCeaIiIgKSTBHRERUSII5IiKi\nQhLMERERFZKPS0WPLFu9luZZFzV6GBHRB6tmv6PRQ4heyIo5IiKiQhLMERERFZJgjoiIqJAE8xCS\n9ClJKyQtlbRE0hsknSRp67o6F0saX7bXlcdmScvL9jGSvtlJ/9d1Uj5e0r/VPZ8u6cKBnFtERAyM\nBPMQkbQXcAiwh+3dgAOBe4CTgOeC2fbBth/tyzFs793BcTcHxgP/9tctIiKiahLMQ6cJeNj2egDb\nDwOHA9sDCyQtAJC0StKEbvraUVKLpNslfbatsG6FPV3SNZIuAG4BZgOvLKv000r1sZLmS7pN0lxJ\nGtDZRkREn+TjUkPnd8BnJP0RuByYZ/sMSR8FDihB3VN7ArsATwI3SLrIdvufftoD2MX2XZKay/YU\nqAU3sDuwM3AfcC2wD/D7vk4uopEeOHdWo4dQadMXntZ9pU1cS0tLo4fwnKyYh4jtdcBUYCbwEDBP\n0jF97O4y22tsPwX8Cti3gzqLbN/VRR+LbN9reyOwBGhuX0HSTEmtklo3PLm2j0ONiIjeyIp5CNne\nALQALZKWATP62lU3zwGe6KaP9XXbG+jgfwu25wBzAEY3Te7oGBGVsN2Rsxs9hEpryReMDCtZMQ8R\nSa+RNLmuaApwN/A4sG0vu3urpBdLGgO8i9qp6K705RgREdEAWTEPnbHAmeWjUM8Cf6J2Wvt9wCWS\n7rN9QA/7WgT8EtgBOKeD68svYHuNpGvLR65+C+S7NSMiKkp2zlBG90Y3TXbTjNMbPYyI6IN8V3bj\nSFpse1pv2uRUdkRERIUkmCMiIiokwRwREVEhufkremTXSeNozXWqiIhBlxVzREREhSSYIyIiKiTB\nHBERUSEJ5oiIiApJMEdERFRIgjkiIqJCEswREREVkmCOiIiokARzREREhSSYIyIiKiTBHBERUSEJ\n5oiIiApJMEdERFRIgjkiIqJC8rOP0SPLVq+ledZFjR5Gpa3Kz2JGxADIijkiIqJCEswREREVkmCO\niIiokGEfzJI2SFoiabmkX0jaupftT+ptm76StErSsvLvFklfkLRVD9qtK4/bS5o/+CONiIhGGfbB\nDDxle4rtXYBngOPrd6qmq3meBPQ2zEf1fpjPOcD2rsCewCuA7/W0oe37bB/ej2MDICk3/UVEVNRI\nCOZ61wCvktQsaaWknwDLgR0lHSTpekk3lpX1WEkfBrYHFkhaANBRvVK+StJXJN0IvEdSi6SvS2qV\ndKuk10v6laTbJX2hu4HaXkftTcS7JL24HOP/SbpB0lJJp7ZvU+a1vGwvlLRz3b4WSdMkbSPpLEmL\nJN0k6dCy/xhJF0i6ErhC0k8kvauu/dy2uhER0TgjZuVUVoFvBy4pRZOBGbYXSpoAfBo40PYTkj4O\nfNT25yR9lNoq9uHO6gGfK32usb1HOd7xwDO2p0k6ETgfmAo8Atwh6eu213Q1ZtuPSboLmCxpXBnz\nnoCACyTtb/vqTprPA/4J+KykJqDJdqukLwFX2j5O0nhgkaTLS5s9gN1sPyLpTcBHgN+UY+8NzOjB\nS91vD5w7aygOM+SmLzyt0UMYFC0tLY0eQsQmZSQE8xhJS8r2NcAPqa2C77a9sJS/EXgtcK0kgC2B\n6zvoq7t689rVv6A8LgNW2L4fQNKdwI5Al8FcqDweVP7dVJ6PpRbUnQXzz4HfAZ+lFtBt154PAt4p\n6eTyfCvgpWX7MtuPANi+StK3JU0EDgN+afvZFwxMmgnMBBj1ook9mEpERPTXSAjmp2xPqS8oofpE\nfRG1UHpfN311V++Jds/Xl8eNddttz7t9bSVtCzQDfyzH/rLtHl1ztr1a0hpJuwFH8Py1dQGH2V7Z\n7lhv6GD8PwGOBt4LHNvBMeYAcwBGN012T8bVE9sdOXuguqqUlnzBSEQMgJF2jbkzC4F9JL0KoFyH\nfXXZ9ziwbQ/qDahy7frbwG9s/w9wKXBc3TXtSZL+rptu5gEfA8bZXlrKLgVOUHl3Imn3Ltr/iNrN\nb9i+pa9ziYiIgbNJBLPth4BjgPMkLaV2enqnsnsOcImkBd3UGygLyg1ci4D/Bv6ljPF3wLnA9ZKW\nUTs1vW2nvdTMp7ba/Xld2eeBLYClklaU5x2y/WfgVuDsvk0lIiIGmuwBO0MZw0z5/PYyYA/ba7uq\nO7ppsptmnD40Axum8l3ZEdGepMW2p/WmzSaxYo6/JulAaqvlM7sL5YiIGDoj4eav6APblwMva/Q4\nIiLihbJijoiIqJCsmKNHdp00jtZcQ42IGHRZMUdERFRIgjkiIqJCEswREREVkmCOiIiokARzRERE\nhSSYIyIiKiTBHBERUSEJ5oiIiApJMEdERFRIgjkiIqJCEswREREVkmCOiIiokARzREREhSSYIyIi\nKiQ/+xg9smz1WppnXdToYQy5Vfmpy4gYYlkxR0REVEiCOSIiokISzBERERUyLINZkiWdU/d8c0kP\nSbqwH31eLGl8H9tOk3RG2Z4uae+6fcdL+kA37beQdGPZ/pSkFZKWSloi6Q19GVMnx9le0vyB6i8i\nIgbecL356wlgF0ljbD8FvBVY3Z8ObR/cj7atQGt5Oh1YB1xX9n23B13sC1wraS/gEGAP2+slTQC2\n7Ou4OhjnfcDhA9VfREQMvGG5Yi4uBtpumX0fcF7bDknbSDpL0iJJN0k6tJQfI+lXki6RdLukr9a1\nWSVpgqRmSbdK+n5Zuf5O0phS5/V1K9nTJC0v5dMlXSipGTge+Eips5+kUySdXOp9WNItpY+f1c3l\nbcBvgSbgYdvrAWw/XMIUSVMlXSVpsaRLJTWV8hZJ3yjHWy5pz1L+plK2pLwG25a5tY155/L6LCnj\nmTygf52IiOiT4bpiBvgZ8Jly+no34Cxgv7LvU8CVto8rp6cXSbq87JsC7A6sB1ZKOtP2Pe36ngy8\nz/YHJf0cOAw4Bzgb+KDt6yXNbj8g26skfRdYZ/trAJLeUldlFvDyshquP21+AHAqtTdKn5H0R+By\nYJ7tqyRtAZwJHGr7IUlHAF8Ejivtt7Y9RdL+5XXYBTgZ+JDtayWNBZ5uN9zjgW/YnitpS2BUJ6/z\ngHng3FmDfYgBN33haY0eQq+0tLQ0eggR0U/DdsVseynQTG21fHG73QcBsyQtAVqArYCXln1X2F5r\n+2ngFuBlHXR/l+0lZXsx0FyCdFvb15fyc/sw7KXAXElHA88CSJoEPGL7SdvrgKnATOAhYJ6kY4DX\nUAvby8qcPg3sUNfveQC2rwZeVMZ6LfCfkj4MjLf9bLuxXA98UtLHgZeVSwIvIGmmpFZJrRueXNuH\n6UZERG8N5xUzwAXA16hd1/3bunIBh9leWV+53Ei1vq5oAx2/Bu3rjBmIwVI79b4/8A/ApyTtSu00\n9qVtFWxvoPZmokXSMmAGtTcHK2zv1Um/bv/c9mxJFwEHU7t+/X+oWzXbPlfSH8qYLpb0L7avbNfJ\nHGAOwOimye2P0WvbHflXJxkqryVfMBIRQ2zYrpiLs4BTbS9rV34pcIIkAUjavb8Hsv0o8HjdXdLv\n7aTq48C27QslbQbsaHsB8HFgHDCW568vI+k17a71TgHuBlYCE8vNYW13ce9cV++IUr4vsNb2Wkmv\ntL3M9leAG4Cd2o3nFcCdts8Azqd2OSAiIhpsWK+Ybd8LnNHBrs8DpwNLSyDeRe1u5/76Z+D7kjYC\nVwEdnd/9L2B+ueHshLryUcA5ksZRW9GfQS3EX2X7tlJnLHBmORX9LPAnYKbtZyQdDpxR2m9e5rei\ntHta0k3AFjx/3fkkSQcAG0u9tpvL2vwT8H5JfwEeAL7Up1ckIiIGlOx+n6HcZEgaW64DI2kW0GT7\nxH70ty9wtO3j+9FHC3By+cjWoBndNNlNM04fzENUUr4rOyL6Q9Ji29N602ZYr5gb4B2SPkHtdbsb\nOKY/ndn+PfD7ARhXRESMEAnmXrA9D5jX6HHUsz290WOIiIiBM9xv/oqIiBhRsmKOHtl10jhac701\nImLQZcUcERFRIQnmiIiICkkwR0REVEiCOSIiokISzBERERWSYI6IiKiQBHNERESFJJgjIiIqJMEc\nERFRIQnmiIiICkkwR0REVEiCOSIiokISzBERERWSYI6IiKiQ/Oxj9Miy1WtpnnVRo4cx5Fblpy4j\nYohlxRwREVEhCeaIiIgKSTBHRERUyCYZzJKaJS1vV3aKpJN70UeLpGnd1DlJ0tZ9HWcPxnDdAPXz\nV69HREQ0xiYZzEPoJKBXwSxpVE/r2t671yOKiIhKSzC3U1bCX5G0SNIfJe1XysdI+pmkWyX9GhhT\n1+YgSddLulHSLySNlfRhYHtggaQFndUr5avKMW8E3lPG8HVJreV4r5f0K0m3S/pC3XHXlcfppc18\nSbdJmitJZd9USVdJWizpUklNdeU3S7oZ+NCQvLgREdGtfFyqY5vb3lPSwcBngQOBfwWetP33knYD\nbgSQNAH4NHCg7SckfRz4qO3PSfoocIDthzurB3yuHHON7T1Kn8cDz9ieJulE4HxgKvAIcIekr9te\n027MuwM7A/cB1wL7SPoDcCZwqO2HJB0BfBE4Djgb+HfbV0s6baBfwAfOnTXQXTbE9IUD/tI0REtL\nS6OHEBE9tKkGs7sp/1V5XAw0l+39gTMAbC+VtLSUvxF4LXBtWaRuCVzfQd/d1ZvXrv4F5XEZsML2\n/QCS7gR2BNoH8yLb95Y6S8q4HwV2AS4rxxwF3C9pPDDe9tWl7U+Bt7cfsKSZwEyAUS+a2MGUIiJi\noG2qwbwG+Jt2ZS8G7irb68vjBrp/jQRcZvt9/az3RLvnbWPYWLfd9ryjMdXXaRu3qIX6Xi8YSC2Y\nu2V7DjAHYHTT5M7ezHRouyNn96Z6ZbXkC0YiYohtkteYba+jtnJ8M4CkFwNvA37fRbOrgSNL/V2A\n3Ur5QmqnjV9V9m0j6dVl3+PAtj2oN1hWAhMl7VWOuYWknW0/Cjwqad9S76hBHkdERPTQJhnMxQeA\n/19O+14JnGr7ji7qfwcYK+lWateFFwPYfgg4BjivnN6+HtiptJkDXCJpQTf1BoXtZ4DDga+Um7yW\nAG13ch8LfKvMX4M5joiI6DnZvTpDGZuo0U2T3TTj9EYPY8jlu7Ijoj8kLbbd5XdetLcpr5gjIiIq\nJ8EcERFRIQnmiIiICtlUPy4VvbTrpHG05nprRMSgy4o5IiKiQhLMERERFZJgjoiIqJAEc0RERIUk\nmCMiIiokwRwREVEhCeaIiIgKSTBHRERUSII5IiKiQhLMERERFZJgjoiIqJAEc0RERIUkmCMiIiok\nwRwREVEh+dnH6JFlq9fSPOuiRg9j0KzKT1pGREVkxRwREVEhCeaIiIgKSTBHRERUyCYXzJI2SFoi\n6WZJN0rau5RvL2l+J22aJR1Z9/wYSd8covFOl3ThAPV13UD0ExERg2eTC2bgKdtTbL8O+ATwZQDb\n99k+vH1lSZsDzcCR7fcNN7b3bvQYIiKia5tiMNd7EfA/8NyqeHnZPkbSBZKuBK4AZgP7lZX2R0rb\n7SVdIul2SV/tqHNJn5F0g6TlkuZIUin/sKRbJC2V9LNSdoqkn0q6vvT5wbquxkqaL+k2SXNV82ZJ\nv6k71lsl/VrSy0r7CZI2k3SNpINKnXXlcaykK8oZg2WSDh3IFzUiIvpuU/y41BhJS4CtgCbgzZ3U\n2wPYzfYjkqYDJ9s+BGrBDUwBdgfWAyslnWn7nnZ9fNP250qbnwKHAP8FzAJebnu9pPF19XcD3ghs\nA9wkqe3zSbsDOwP3AdcC+wALgG9Lmmj7IeBY4Czbd0v6CvAdYBFwi+3ftRvX08C7bT8maQKwUNIF\ntt39y9d3D5w7azC775fpC09r9BA61dLS0ughRMQQ2hRXzG2nsncC3gb8pG0l285lth/pop8rbK+1\n/TRwC/CyDuocIOkPkpZRewOwcylfCsyVdDTwbF39820/ZfthasG7ZylfZPte2xuBJUBzCdGfAkeX\ncN8L+C2A7R9QOxtwPHByB+MS8CVJS4HLgUnAS/6qkjRTUquk1g1Pru3ipYiIiIGyKa6Yn2P7+rJi\nnNjB7ie6ab6+bnsD7V5LSVsB3wam2b5H0inUVukA7wD2B/4B+JSkXduG1H6I3RzrbGor8KeBX9h+\nthx7a2CHUmcs8Hi7fo+iNueptv8iaVXd2J4/uD0HmAMwumlyv1fT2x05u79dDJqWfMFIRFTEprhi\nfo6knYBRwJpuqj4ObNvL7tuC7mFJY4HDyzE3A3a0vQD4ODCOWngCHCppK0l/C0wHbujqALbvo3Z6\n+9PUQrrNV4C5wGeA73fQdBzwYAnlA+h4tR8REQ2wKa6Y264xQ+2U7gzbGzo+m/2cpcAGSTcDP6Lc\nMNYV249K+j6wHHiA50N2FHCOpHHl+GeUum3HWQBMAD5v+z5Jr+7mUHOBibZvBZD0JuD1wD5lXodJ\nOtb22e3a/Fc5xd4K3NbdfCIiYmhokO/3iR4qp7rX2f5aL9t9E7jJ9g8HZWDF6KbJbppx+mAeoqHy\nXdkRMRgkLbY9rTdtNsUV84ghaTG1a+H/0eixRETEwEgwV4TtU/rQZuogDCUiIhpok775KyIiomqy\nYo4e2XXSOFpzHTYiYtBlxRwREVEhCeaIiIgKSTBHRERUSII5IiKiQhLMERERFZJgjoiIqJAEc0RE\nRIUkmCMiIiokwRwREVEhCeaIiIgKSTBHRERUSII5IiKiQhLMERERFZJgjoiIqJD87GP0yLLVa2me\ndVGjh9Erq/IzlRExDGXFHBERUSEJ5oiIiApJMEdERFRIgjkiIqJCEsx9IGmDpCWSlkv6haStGz2m\nnpC0vaT5jR5HRER0LsHcN0/ZnmJ7F+AZ4PhGD6gnbN9n+/BGjyMiIjqXj0v13zXAbpKagd8Cvwf2\nBlYDh9p+StIrgW8BE4EngQ/avk3Sj4ALbc8HkLTO9lhJ04FTgUeBXYGfA8uAE4ExwLts31GOeRYw\nAXgIONb2f5d+HwOmAdsBH7M9v9S/0PYuZfunwDZlHv9u+7reTv6Bc2f1tsmQmb7wtEYPoVMtLS2N\nHkJEVFRWzP0gaXPg7dRCE2Ay8C3bO1ML1cNK+RzgBNtTgZOBb/eg+9dRW4n/PfB+4NW29wR+AJxQ\n6pwJ/Nj2bsBc4Iy69k3AvsAhwOwO+n8QeKvtPYAj2rVtm99MSa2SWjc8ubYHQ46IiP7Kirlvxkha\nUravAX4IbA/cZbutfDHQLGkstRX0LyS1tR/dg2PcYPt+AEl3AL8r5cuAA8r2XsA/lu2fAl+ta/8b\n2xuBWyS9pIP+twC+KWkKsAF4dfsKtudQe1PB6KbJ7miQ2x3ZUeZXQ0u+YCQihqEEc988ZXtKfUEJ\n3fV1RRuonXbeDHi0ff3i2bIfSZsBW9btq+9rY93zjfTs71bfXh3s/wjwZ2or882Ap3vQZ0REDLKc\nyh5kth8D7pL0HgDVvK7sXgVMLdvvpLaK7Y3rgPeW7aOord57ahxwf1lVvx8Y1ctjR0TEIEgwD42j\ngH+WdDOwAji0lH8feFMp3wt4opf9ngAcK2kptXA9sRdtvw3MKMfeqQ/HjoiIQSC7w0uHES8wummy\nm2ac3uhh9Ep+xCIiGk3SYtvTetMmK+aIiIgKSTBHRERUSO7Kjh7ZddI4WnNqOCJi0GXFHBERUSEJ\n5oiIiApJMEdERFRIgjkiIqJCEswREREVkmCOiIiokARzREREhSSYIyIiKiTBHBERUSEJ5oiIiApJ\nMEdERFRIgjkiIqJCEswREREVkmCOiIiokARzREREheT3mKNHlq1eS/Osiwal71X5neeIiOdkxRwR\nEVEhCeaIiIgKSTBHRERUSLfBLGk7ST+TdIekxZIulvTqoRhcN+O6WNL48u/fethm3SCMY7qkC8v2\naEmXS1oi6YiBPlYXY/hku+fXDdWxIyJiYHUZzJIE/Bposf1K21OBTwAvGYrBdcX2wbYfBcYDPQrm\nIbA7gO0ptuf1pIGkUQNw3BcEs+29B6DPiIhogO7uyj4A+Ivt77YV2L4ZngvtrwJvBwx8wfY8SdOB\nU4CHgV2AxcDRti3pYOA/gSeAa4FX2D5E0jbAmaX+FsApts+XdAzwTmBr4JXAr21/rBx/FTANmA28\nUtIS4DLgVOB84G9KX5+2fX5nE5TUDFwCLAT2Bm4Azi79/B1wlO1FkvYEvgFsBTwFHGt7ZV0/fwec\nA0wsYzkMaAa+Vl7nG4B/tb2+jH0e8Fbgq5KOB24C9gO2AT5A7Q3QrsA8258ux/gNsGMZwzdsz5E0\nGxhTjrnC9lGS1tke25e/UWev01/W3MsD587qbHe/TF942qD029LSMij9RkQMpu6Cue0/2h35R2AK\n8DpgAnCDpKvLvt2BnYH7qAXwPpJage8B+9u+S9J5dX19CrjS9nGSxgOLJF1e9k0p/a0HVko60/Y9\ndW1nAbvYngIgaXPg3bYfkzQBWCjpgq5CB3gV8B7gOGoBeiSwL7U3BZ8E3gXcBuxn+1lJBwJfoha+\nANh+UNL/BU4ubza2AlqAt9j+o6SfAP8KnF6arLG9Rxnz8cAztqdJOpHaG4upwCPAHZK+bnsNcJzt\nRySNKa/3L23PkvTvbfPvz98I+H19Y0kzgZkAjNqii5cvIiIGSn8+x7wvcJ7tDcCfJV0FvB54DFhk\n+16AspJrBtYBd9q+q7Q/j7b/6MNBwDslnVyebwW8tGxfYXtt6esW4GVAfTC3J+BLkvYHNgKTqJ16\nf6CLNnfZXlaOsaIc05KWlbEDjAN+LGkytdVnd0n1mtLvH8vzHwMf4vlgbn+q+4LyuIzayvf+Mp47\nqa2S1wAflvTuUm9HYHIp70xv/0YvCGbbc4A5AKObJnu7I2d3M+W+acnnmCMintNdMK8ADu9Dv+vr\ntjf04DgCDqs/NQwg6Q196OsoYCIw1fZfymnjrXox3o11zzfWHe/zwALb7y6nv1u66bM7T3Qyhvrj\nPzeGcvr5QGAv209KaqH7eXWlt69rREQMge7uyr4SGF1OaQIgaTdJ+wHXAEdIGiVpIrA/sKiLvlYC\nryihBlB/1/KlwAnlmiiSdu/FHB4Htq17Pg54sITyAdRW2ANhHLC6bB/Tg/orgWZJryrP3w9c1c/j\n/08J5Z2AN9bt+4ukjlbwvf0bRUREg3UZzOW67LuBA8vHpVYAX6Z2WvjXwFLgZmoB/jHbnZ4utv0U\ntbunL5G0mFqgri27P0/t1PDScozP93QC5drrtZKWSzoNmAtMK6ehP0Dt2vBA+CrwZUk30YPVpe2n\ngWOBX5SxbAS+23WrLl1CbeV8K7Ub3hbW7ZtD7bWb265Nr/5GERHReOr6nqgBPpg01va6sjL+FnC7\n7a8P2QCiz0Y3TXbTjNO7r9gH+a7siBipJC22Pa03bYb6m78+2PaxHmqnZr83xMePiIiotCG94aes\njrNCjoiI6ETuxI0e2XXSOFpzyjkiYtDlRywiIiIqJMEcERFRIQnmiIiICkkwR0REVEiCOSIiokIS\nzBERERWSYI6IiKiQIf1Kzhi+JD1O7Yc5RqoJwMONHsQgyvyGt5E8v5E8N4DX2N62+2rPyxeMRE+t\n7O33vQ4nklozv+Er8xu+RvLcoDa/3rbJqeyIiIgKSTBHRERUSII5empOowcwyDK/4S3zG75G8tyg\nD/PLzV8REREVkhVzREREhSSYo1uS3iZppaQ/SZrV6PH0l6SzJD0oaXld2YslXSbp9vL4N40cY19J\n2lHSAkm3SFoh6cRSPlLmt5WkRZJuLvM7tZSPiPm1kTRK0k2SLizPR8z8JK2StEzSkrY7lkfY/MZL\nmi/pNkm3Stqrt/NLMEeXJI0CvgW8HXgt8D5Jr23sqPrtR8Db2pXNAq6wPRm4ojwfjp4F/sP2a4E3\nAh8qf6+RMr/1wJttvw6YArxN0hsZOfNrcyJwa93zkTa/A2xPqfuY1Eia3zeAS2zvBLyO2t+xV/NL\nMEd39gT+ZPtO288APwMObfCY+sX21cAj7YoPBX5ctn8MvGtIBzVAbN9v+8ay/Ti1/yhMYuTMz7bX\nladblH9mhMwPQNIOwDuAH9QVj5j5dWJEzE/SOGB/4IcAtp+x/Si9nF+CObozCbin7vm9pWykeYnt\n+8v2A8BLGjmYgSCpGdgd+AMjaH7lNO8S4EHgMtsjan7A6cDHgI11ZSNpfgYul7RY0sxSNlLm93Lg\nIeDsciniB5K2oZfzSzBHtOPaRxWG9ccVJI0FfgmcZPux+n3DfX62N9ieAuwA7Clpl3b7h+38JB0C\nPGh7cWd1hvP8in3L3+/t1C617F+/c5jPb3NgD+A7tncHnqDdaeuezC/BHN1ZDexY93yHUjbS/FlS\nE0B5fLDB4+kzSVtQC+W5tn9VikfM/NqUU4QLqN0vMFLmtw/wTkmrqF02erOkcxg588P26vL4IPBr\napfLRsr87gXuLWdxAOZTC+pezS/BHN25AZgs6eWStgTeC1zQ4DENhguAGWV7BnB+A8fSZ5JE7frW\nrbb/s27XSJnfREnjy/YY4K3AbYyQ+dn+hO0dbDdT+//albaPZoTMT9I2krZt2wYOApYzQuZn+wHg\nHkmvKUVvAW6hl/PLF4xEtyQdTO261yjgLNtfbPCQ+kXSecB0ar9q82fgs8BvgJ8DLwXuBv7Jdvsb\nxCpP0r5P8ndIAAAAjklEQVTANcAynr9G+Ulq15lHwvx2o3bzzChqC4uf2/6cpL9lBMyvnqTpwMm2\nDxkp85P0CmqrZKid9j3X9hdHyvwAJE2hduPelsCdwLGU/63Sw/klmCMiIiokp7IjIiIqJMEcERFR\nIQnmiIiICkkwR0REVEiCOSIiokISzBERERWSYI6IiKiQBHNERESF/C89HlGSlTr9eQAAAABJRU5E\nrkJggg==\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7ff22a940c88>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"csmf_hbar_w_error(df.iloc[:100])\n",
"plt.axis(xmax=60);"
]
},
{
"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.6.1"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment