Skip to content

Instantly share code, notes, and snippets.

@fonnesbeck

fonnesbeck/ELSO Logistic Regression.ipynb

Last active August 29, 2015 14:22
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 fonnesbeck/ac400c96e1d4c9014662 to your computer and use it in GitHub Desktop.
Save fonnesbeck/ac400c96e1d4c9014662 to your computer and use it in GitHub Desktop.
Download ZIP
Display the source blob
Display the rendered blob
Raw
ELSO Logistic Regression.ipynb
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"'3.0'"
]
},
"execution_count": 1,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"%matplotlib inline\n",
"import numpy as np\n",
"import pandas as pd\n",
"import matplotlib.pyplot as plt\n",
"import pymc3 as pm\n",
"import pdb\n",
"\n",
"pm.__version__"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"flu = pd.read_csv('data/flu_organism.csv', index_col=0)"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"(2933, 63)"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"flu.shape"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"(1654,)"
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"flu.PatientID.unique().shape"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Turn into table of unique patients, with appropriate organism columns\n",
"\n",
"Create indicator for flu-only."
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"flu_only = flu.groupby('PatientID')['OrganismName'].apply(lambda s: \n",
" len([x for x in s if str(x).startswith('Influenza')])==len(s)).astype(bool)"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"PatientID\n",
"001B0EC6-66F5-48C9-8516-52D2ABFE28AA False\n",
"0034181F-59AB-4B2C-A06A-A83EE1DF1A17 True\n",
"004AE387-1096-410D-A7CB-F9A9086AA385 False\n",
"0084DCF0-8E3C-4F1C-8650-BE795A57EF2F False\n",
"0084E706-90D0-4314-8C99-496171AB639D True\n",
"Name: OrganismName, dtype: bool"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"flu_only.head()"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"flu_only.name = 'flu_only'"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"flu = flu.merge(pd.DataFrame(flu_only), left_on='PatientID', right_index=True)"
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"flu_staph_organism = flu.groupby('PatientID')['OrganismName'].apply(lambda s: \n",
" len([x for x in s if str(x).startswith('Influenza') or \n",
" str(x).startswith('Staphylococcus aureus')])>1).astype(bool)"
]
},
{
"cell_type": "code",
"execution_count": 14,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"flu_staph_organism.name = 'flu_staph_org'"
]
},
{
"cell_type": "code",
"execution_count": 15,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"flu = flu.merge(pd.DataFrame(flu_staph_organism), left_on='PatientID', right_index=True)"
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"0.13059250302297462"
]
},
"execution_count": 12,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"flu_staph_organism.mean()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"See if there are any flu-only or flu-other individuals with an ICD9 diagnosis of S. Aureus"
]
},
{
"cell_type": "code",
"execution_count": 18,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"flu_other = flu.groupby('PatientID')['OrganismName'].apply(lambda s: \n",
" len([x for x in s if str(x).startswith('Influenza') or \n",
" str(x).startswith('Other')])>1).astype(bool)"
]
},
{
"cell_type": "code",
"execution_count": 19,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"flu_other.name = 'flu_other'"
]
},
{
"cell_type": "code",
"execution_count": 20,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"flu = flu.merge(pd.DataFrame(flu_other), left_on='PatientID', right_index=True)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"There do not appear to be cases of flu-only or flu-other via organism that have a S.Aureus diagnosis via ICD9"
]
},
{
"cell_type": "code",
"execution_count": 21,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"0"
]
},
"execution_count": 21,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"(flu_other & flu.s_aureus_icd9).sum()"
]
},
{
"cell_type": "code",
"execution_count": 22,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"0"
]
},
"execution_count": 22,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"(flu_only & flu.s_aureus_icd9).sum()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Same for pneumo"
]
},
{
"cell_type": "code",
"execution_count": 23,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"0"
]
},
"execution_count": 23,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"(flu_other & flu.pneumo_icd9).sum()"
]
},
{
"cell_type": "code",
"execution_count": 24,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"PatientID \n",
"001B0EC6-66F5-48C9-8516-52D2ABFE28AA Klebsiella pneumoniae 1\n",
" Enterococcus 1\n",
"0034181F-59AB-4B2C-A06A-A83EE1DF1A17 Influenza A 1\n",
"004AE387-1096-410D-A7CB-F9A9086AA385 Candida albicans 1\n",
" Influenza A 1\n",
" Pseudomonas aerugenosa 1\n",
"0084DCF0-8E3C-4F1C-8650-BE795A57EF2F Influenza A 1\n",
" Staphylococcus aureus 1\n",
"0084E706-90D0-4314-8C99-496171AB639D Influenza B 1\n",
"00EFD11D-650A-4317-8C4F-4B7378AD7946 Other 2\n",
"0118B486-8ACE-4BFF-8ED6-753D96E6B441 Influenza A 1\n",
"011B5A5D-9967-44A7-92F9-02F82D31216D Candida albicans 2\n",
" Influenza A 1\n",
" Epstein-Barr virus (EBV) 1\n",
"01AC9644-7CDF-4E28-A5F4-FE5B7F9FA22A Influenza A 1\n",
" Streptococcus pneumoniae 1\n",
" Moraxella catarrhalis (Branhamella) 1\n",
" Staphylococcus, coag neg 1\n",
"01F84F52-92D3-403B-B8AC-21A72A4BC771 Influenza A 1\n",
"02373407-74E6-4A1E-9629-F4DCBF707B0C Influenza A 1\n",
"025D87C9-CF70-4B9A-94BE-B15EA2185472 Eschericha coli 2\n",
" Influenza A 1\n",
" Enterococcus 1\n",
"02A94D2E-768A-4ABB-86BF-C8A188D698CC Influenza A 1\n",
"02C9FF5A-A4D6-461A-A650-E3EB405D62AE Influenza A 2\n",
"02D33150-FDDB-4F12-AAAB-D6F4B8765EBF Influenza A 1\n",
"031C1231-B585-478F-A98F-35E017B788CD Other 2\n",
"035A053E-91B1-4B68-87D4-F92C675FE6E3 Other 2\n",
"03AA4AD0-1B0E-4CAD-BD76-03FE52AA580D Influenza A 1\n",
" Pneumocystis carinii 1\n",
" ..\n",
"FC9F84B4-B2CD-4C28-B824-D5CE1AF7BEF1 Gram positive, other 2\n",
" Aspergillus fumigatus 1\n",
" Candida albicans 1\n",
" Enterococcus 1\n",
" Pseudomonas fluorescens 1\n",
" Other 1\n",
" Yeast sp. 1\n",
" Stenotrophomonas maltophilia (Xanthomonas) 1\n",
"FCBAEC19-EFA1-4A95-A98D-1D8FE65B53CE Influenza A 1\n",
"FCD795A0-9FD1-4C20-961B-04BEF1B23E96 Influenza A 1\n",
" Staphylococcus aureus 1\n",
"FCDF1967-C307-4FF1-8CF1-7AA02FE0BE47 Influenza A 1\n",
"FD2FA339-7417-4F13-811B-2E1873E529C7 Influenza A 1\n",
"FDC32D1D-0D67-4752-A2E1-72E790FCB797 Pseudomonas aerugenosa 1\n",
" Influenza A 1\n",
"FDED68D1-1DFF-4899-9CF9-797E6AE70D18 Influenza A 1\n",
"FDFF23E3-C0CB-468D-98D9-998192CA7590 Eschericha coli 1\n",
" Candida albicans 1\n",
" Influenza A 1\n",
" Yeast sp. 1\n",
"FEA94F00-083A-479C-B862-05274A9842B5 Streptococcus, group A 1\n",
" Influenza A 1\n",
" Staphylococcus aureus 1\n",
"FF0D88BF-E131-4E56-9E78-CC6B60E461FD Influenza A 1\n",
" Staphylococcus, coag neg 1\n",
" Yeast sp. 1\n",
"FF0FB831-9804-4506-AE0B-B447AB54D4D5 Hemophilus influenza 1\n",
"FF528500-5274-4225-96C2-CF1153EE3AB0 Other 1\n",
"FF6CD7E8-8956-4C15-AF6C-3033A38C48D7 Influenza A 1\n",
" Staphylococcus epidermidis 1\n",
"dtype: int64"
]
},
"execution_count": 24,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"flu.groupby('PatientID')['OrganismName'].value_counts()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Create indictors for coinfection type"
]
},
{
"cell_type": "code",
"execution_count": 25,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"for org in flu.Type.unique():\n",
" flu[org] = (flu.Type==org).astype(int)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Create data frame of unique patients"
]
},
{
"cell_type": "code",
"execution_count": 26,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"flu_unique = flu.drop_duplicates(subset=['PatientID']).set_index('PatientID')"
]
},
{
"cell_type": "code",
"execution_count": 27,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"99"
]
},
"execution_count": 27,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"flu_unique.s_aureus_icd9.sum()"
]
},
{
"cell_type": "code",
"execution_count": 28,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"flu_unique['flu_only'] = flu_only"
]
},
{
"cell_type": "code",
"execution_count": 29,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"0.34159613059250304"
]
},
"execution_count": 29,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"flu_unique.flu_only.mean()"
]
},
{
"cell_type": "code",
"execution_count": 23,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"0.13059250302297462"
]
},
"execution_count": 23,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"flu_unique.flu_staph.mean()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Several missing values for admission to time on ECMO"
]
},
{
"cell_type": "code",
"execution_count": 30,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"0.16686819830713423"
]
},
"execution_count": 30,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"flu_unique.AdmitToTimeOnHours.isnull().mean()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Since we need this field to calculate event time, we will have to drop individuals with missing values from the survival analysis."
]
},
{
"cell_type": "code",
"execution_count": 31,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"flu_unique = flu_unique.dropna(subset=['AdmitToTimeOnHours', 'HoursECMO'])\n",
"assert not flu_unique.AdmitToTimeOnHours.isnull().sum()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Create variables for use in analysis"
]
},
{
"cell_type": "code",
"execution_count": 45,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"# ECMO Type\n",
"flu_unique['VA'] = flu_unique.Mode.isin(['VA', 'VV-VA'])\n",
"# Set \"Other\" type to NA (there are only a couple)\n",
"flu_unique.loc[flu_unique.Mode=='Other', 'VA'] = None"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Create oxygen index"
]
},
{
"cell_type": "code",
"execution_count": 46,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"flu_unique['OI'] = flu_unique.FiO2 * flu_unique.MAP / flu_unique.PO2"
]
},
{
"cell_type": "code",
"execution_count": 47,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x117075f60>"
]
},
"execution_count": 47,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": [
"iVBORw0KGgoAAAANSUhEUgAAAX4AAAEBCAYAAAB/rs7oAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\n",
"AAALEgAACxIB0t1+/AAAEtpJREFUeJzt3WGs3Xddx/H3924wZqS2mcgQp2Vr0fAIWQyRlazrWWdU\n",
"EkxMxJjAUB7IZixoQBZchDZgAuoSEsLwEYQHoCZKeEJcy+m2FhwjYSGLyZhclkZcnJimu92DdUD2\n",
"88E5xUt3Tu895/7+u/f3/b1fSZP+f/ufe36fnfZz/vd7zzmNUgqSpH6sbPcGJEkvLotfkjpj8UtS\n",
"Zyx+SeqMxS9JnbH4JakzGxZ/RHwmIh6KiJMR8Y7p2q0RcToiTkXEoXXnjmatS5J2jis3cU4Bfq+U\n",
"8j2AiAjgKDACArgPODldP3bp+hCbliQtbzPFH/zkdwb7gcdLKRcAImI1IvZNz3nBeilltfamJUnL\n",
"20zxPwN8PiLOAn8OXAOsRcQ9TJ4U1qZrK3PWLX5J2kE2LP5SyhGAiHg98DfAXwC7gTuYFPy9wFkm\n",
"xT9rXZK0g2zmiv+iC8APge8yGffApOD3lVJWI2Jl1vq8LzYej/2QIElawmg0iq3cfsPij4h/AF7F\n",
"ZOTzJ6WU5yPiKPAVJj/4PQowb/1ytrp5SepNjYvmzYx6fn/G2gngxGbXexQRB0spD2z3PoaQORuY\n",
"r3XZ89XgG7gkqTMW/0AyX3Fkzgbma132fDVY/JLUGYt/IBFxcLv3MJTM2cB8rcuerwaLX5I6Y/EP\n",
"JPOcMXM2MF/rsuerweKXpM5Y/APJPGfMnA3M17rs+Wqw+CWpMxb/QDLPGTNnA/O1Lnu+Gix+SeqM\n",
"xT+QzHPGzNnAfK3Lnq8Gi1+SOmPxDyTznDFzNjBf67Lnq8Hil6TOLPIvcFV37ty5pZ549uzZ83zt\n",
"vdSW+TPBM2cD87Uue74atrX4/+Xfv//Eorf51Z9/+Q/37Nmzf+MzJUmzbGvxf/rrT/7Sorc5ctN1\n",
"j73hhiF2U1fmK47M2cB8rcuerwZn/JLUGYt/IJlfS5w5G5ivddnz1WDxS1JnLP6BZJ4zZs4G5mtd\n",
"9nw1WPyS1BmLfyCZ54yZs4H5Wpc9Xw0WvyR1xuIfSOY5Y+ZsYL7WZc9Xg8UvSZ2x+AeSec6YORuY\n",
"r3XZ89Vg8UtSZyz+gWSeM2bOBuZrXfZ8NVj8ktQZi38gmeeMmbOB+VqXPV8NFr8kdcbiH0jmOWPm\n",
"bGC+1mXPV4PFL0mdsfgHknnOmDkbmK912fPVYPFLUmc2VfwR8dKIOBMRd06Pb42I0xFxKiIOrTtv\n",
"NGu9R5nnjJmzgflalz1fDZv9x9bfDXwTICICOAqMgADuA05O149dul57w5Kkrdnwij8irgZuA740\n",
"XdoPPF5KuVBKeRZYjYh9l1nvUuY5Y+ZsYL7WZc9Xw2au+I8AnwReOT2+BliLiHuYXNmvTddW5qyv\n",
"1t60JGl5l73ij4hdwJtLKf96cQk4C+wGPjj9tWe6Nm99EBFxcP0z+047vri2U/ZT87iU8sBO2o/5\n",
"zNdTvhqilDL/P0b8FvBnwP8C1wNXAO8CPgXcyuSJ43gp5UBErACnLl2f97XH43G565FYeMNHbrru\n",
"sbe/af/rFr6hJCUwHo/LaDRavDzXuewVfynly6WUw6WUP2BS9p8ppTzK5Ie4X2HyA9yj03Ofn/7+\n",
"J9Z7VfPZeafJnA3M17rs+WrY7Kt6KKV8bt3vjwPHZ5xzAjhRZ2uSpCH4Bq6BZH4tceZsYL7WZc9X\n",
"g8UvSZ2x+AeSec6YORuYr3XZ89Vg8UtSZyz+gWSeM2bOBuZrXfZ8NVj8ktQZi38gmeeMmbOB+VqX\n",
"PV8NFr8kdcbiH0jmOWPmbGC+1mXPV4PFL0mdsfgHknnOmDkbmK912fPVYPFLUmcs/oFknjNmzgbm\n",
"a132fDVY/JLUGYt/IJnnjJmzgflalz1fDRa/JHXG4h9I5jlj5mxgvtZlz1eDxS9JnbH4B5J5zpg5\n",
"G5ivddnz1WDxS1JnLP6BZJ4zZs4G5mtd9nw1WPyS1BmLfyCZ54yZs4H5Wpc9Xw0WvyR1xuIfSOY5\n",
"Y+ZsYL7WZc9Xg8UvSZ2x+AeSec6YORuYr3XZ89Vg8UtSZyz+gWSeM2bOBuZrXfZ8NVj8ktQZi38g\n",
"meeMmbOB+VqXPV8NFr8kdcbiH0jmOWPmbGC+1mXPV4PFL0mdsfgHknnOmDkbmK912fPVYPFLUmc2\n",
"LP6I+EhEnIyIExHxmunaKCJOR8SpiDi07tyZ6z3KPGfMnA3M17rs+Wq4cqMTSil3A0TETcAHIuIO\n",
"4BgwAgK4DzgZETFrfaB9S5KWtMio543At4H9wOOllAullGeB1YjYd5n1LmWeM2bOBuZrXfZ8NWx4\n",
"xQ8QEQ8CrwIOADcAaxFxD5Mr+zXgGiZPIrPWVwfYtyRpSZsq/lLKzRHxa8DngD8FdgN3MCn4e4Gz\n",
"TIp/1nqXMs8ZM2cD87Uue74aNlX8U/8DFOC7TMY6MCn4faWU1YhYmbVebaeXuPjt3MUH2WOPPfa4\n",
"l+OtilLK5U+I+EfgZ4FngfeUUr4bEYeBDzF5IjhWSjkxPXfm+izj8bjc9UgsvOEjN1332NvftP91\n",
"C9/wRRYRB7NeeWTOBuZrXfZ84/G4jEajxctznc28qudtM9ZOAC8o9XnrkqSdwzdwDSTzFUfmbGC+\n",
"1mXPV4PFL0mdsfgHkvm1xJmzgflalz1fDRa/JHXG4h9I5jlj5mxgvtZlz1eDxS9JnbH4B5J5zpg5\n",
"G5ivddnz1WDxS1JnLP6BZJ4zZs4G5mtd9nw1WPyS1BmLfyCZ54yZs4H5Wpc9Xw0WvyR1xuIfSOY5\n",
"Y+ZsYL7WZc9Xg8UvSZ2x+AeSec6YORuYr3XZ89Vg8UtSZyz+gWSeM2bOBuZrXfZ8NVj8ktQZi38g\n",
"meeMmbOB+VqXPV8NFr8kdcbiH0jmOWPmbGC+1mXPV4PFL0mdsfgHknnOmDkbmK912fPVYPFLUmcs\n",
"/oFknjNmzgbma132fDVY/JLUGYt/IJnnjJmzgflalz1fDRa/JHXG4h9I5jlj5mxgvtZlz1fDldu9\n",
"gWU8/J0nP/vU+ef2Lnq7a3dddeaN+1/9zvo7kqR2NFn8T51/bu9HTp65edHb3X1o7wC7mS0iDma9\n",
"8sicDczXuuz5anDUI0mdsfgHkvmKI3M2MF/rsuerweKXpM5Y/APJ/FrizNnAfK3Lnq8Gi1+SOmPx\n",
"DyTznDFzNjBf67Lnq2HD4o+IT0fE/RHxQES8Zro2iojTEXEqIg6tO3fmuiRp59iw+Esp7y6l3AIc\n",
"Bd4fEQEcAw4DvwF8GGDeeq8yzxkzZwPztS57vhoWeQPXM8APgP3A46WUCwARsRoR+5g8ibxgvZSy\n",
"WnvTkqTlLVL87wI+AVwDrEXEPUAAa9O1lTnrXRZ/5jlj5mxgvtZlz1fDpoo/It7C5Gr+2xHxWmA3\n",
"cAeTgr8XOMuk+GetS5J2kA2LPyJuBA6WUt43XVplMu6BScHvK6WsRsTKrPXaG96K80+f233x9xfn\n",
"gBevDgY4fi/wrQG//rYdr5+h7oT9mM98PeWjgiilXP6EiCeA7wHPA4+WUt4TEbcBfwUU4Fgp5cT0\n",
"3MPAhy5dn2U8Hpe7HomFN3zkpuse23XVFd9f8kPaHnzrjdcfXPhOl5D5g6IyZwPztS57vvF4XEaj\n",
"0eLluc6GV/yllOtnrB0Hjs9YPwHMLfueZP6DlzkbmK912fPV4Bu4JKkzFv9AMr+WOHM2MF/rsuer\n",
"weKXpM5Y/APJPGfMnA3M17rs+Wqw+CWpMxb/QDLPGTNnA/O1Lnu+Gix+SeqMxT+QzHPGzNnAfK3L\n",
"nq8Gi1+SOmPxDyTznDFzNjBf67Lnq8Hil6TOWPwDyTxnzJwNzNe67PlqsPglqTMW/0AyzxkzZwPz\n",
"tS57vhosfknqjMU/kMxzxszZwHyty56vBotfkjpj8Q8k85wxczYwX+uy56vB4pekzlj8A8k8Z8yc\n",
"DczXuuz5arD4JakzFv9AMs8ZM2cD87Uue74aLH5J6ozFP5DMc8bM2cB8rcuerwaLX5I6Y/EPJPOc\n",
"MXM2MF/rsuerweKXpM5Y/APJPGfMnA3M17rs+Wqw+CWpMxb/QDLPGTNnA/O1Lnu+Gix+SeqMxT+Q\n",
"zHPGzNnAfK3Lnq8Gi1+SOmPxDyTznDFzNjBf67Lnq8Hil6TOWPwDyTxnzJwNzNe67PlqsPglqTMb\n",
"Fn9EHIiIb0TEx9etjSLidESciohDG633KPOcMXM2MF/rsuer4cpNnHMV8NfAmwAiIoBjwAgI4D7g\n",
"5Lz1AfYsSdqCDa/4Sylj4Ny6pf3A46WUC6WUZ4HViNh3mfUuZZ4zZs4G5mtd9nw1bOaK/1LXAGsR\n",
"cQ+TK/u16drKnPXVSnuVJFWwTPGfBXYDdzAp+Hunaytz1neM80+f233x9xfngBevDgY4fi/wrQG/\n",
"/rYdr5+h7oT9mM98PeWjgiilbHxSxM3AW0op74+IFeAUcCuTsj9eSjkwb33e1xyPx+WuR2LhDR+5\n",
"6brHdl11xfc/cvLMzYve9u5Dex98643XH1z4TpcQEQezfsuZORuYr3XZ843H4zIajRYvz3U2vOKP\n",
"iA8Avwm8MiJ2lVL+OCKOAV8BCnAUoJTyfEQcvXS9V5n/4GXOBuZrXfZ8NWxY/KWUjwEfu2TtOHB8\n",
"xrkngBPVdidJqs43cA0k82uJM2cD87Uue74aLH5J6ozFP5DMc8bM2cB8rcuerwaLX5I6Y/EPJPOc\n",
"MXM2MF/rsuerweKXpM5Y/APJPGfMnA3M17rs+Wqw+CWpMxb/QDLPGTNnA/O1Lnu+Gix+SeqMxT+Q\n",
"zHPGzNnAfK3Lnq8Gi1+SOmPxDyTznDFzNjBf67Lnq8Hil6TOWPwDyTxnzJwNzNe67PlqsPglqTMW\n",
"/0AyzxkzZwPztS57vhosfknqjMU/kMxzxszZwHyty56vBotfkjpj8Q8k85wxczYwX+uy56vB4pek\n",
"zlj8A8k8Z8ycDczXuuz5arD4JakzFv9AMs8ZM2cD87Uue74aLH5J6ozFP5DMc8bM2cB8rcuer4Yr\n",
"t3sDL6bdV1/5K1/65hMPLHPba3dddeaN+1/9zro7kqQXX1fFf+7ZH73soyfP3LzMbe8+tHeh8yPi\n",
"YNYrj8zZwHyty56vBkc9ktQZi38gma84MmcD87Uue74aLH5J6ozFP5DMryXOnA3M17rs+Wqw+CWp\n",
"Mxb/QDLPGTNnA/O1Lnu+Gix+SepM9dfxR8QI+DBQgA+XUk7Wvo/tsOibv84/fW73rt17noZ8b/7K\n",
"/jpp87Ute74aqhZ/RARwDBgBAdwHpCj+5d78tQYs/uYvSRpS7VHPfuDxUsqFUsqzwGpE7Kt8H9pm\n",
"2a+mzNe27PlqqD3quQZYi4h7mFzxr03XVivfT1OW/Yyga1/+0l946pkf/Ncy97nsbbd0n8lGWlJW\n",
"tYv/LLAbuINJ8d87XZvp7W+4duGC+eVX/NRL/vv8c0tvcDss+xlBf3lo79pHT565YZn7XPa2W7nP\n",
"v/3tfct/CN52PFHNue36n8/MvN0WnuAe/s6Tn33q/HN7l7nt0v+PLtlv5hn4w9958rP/8Z9Pvf5y\n",
"j98svV20RCml3heLWAFOAbcyGSMdL6UcmHXueDyud8eS1JHRaBRbuX3V4geIiMPAh5i8qudYKeVE\n",
"1TuQJG1J9eKXJO1svoFLkjpj8UtSZyx+SerMthR/RIwi4nREnIqIQ9uxh62KiE9HxP0R8UBEvGa6\n",
"NjNXq3kj4qURcSYi7pwe35olX0S8OiJOTvf8d9O1TPn+MCIejoivRsQt07Vm80XEgYj4RkR8fN3a\n",
"Qn/fdnLOOfle0DHT9a3nK6W8qL+YvL7/a8DLgKuBUy/2HirnuQX41LxcLecFjgD/DNyZLR/wBeDX\n",
"1x1ny/cokwu7XcC/tZ6PycfA/A7w8WUer52e89J8l/y3W4B7a+bbjn9s/ccf6wAQEasRsa+U0uq7\n",
"e58BfsCcXEz+8jWXNyKuBm4D/gn4aRLlm77fZF8p5aF1y2nyTT3K5P00P8fkM7OazldKGUfE+jdB\n",
"LpRn3vpOyTkj33rPABfftVrlcdyO4s/2sQ7vAj7B/Fwrc9Z3et4jwCeBV06PM+V7BfCyiPgikyvi\n",
"TwJPkScfwGngdib7/wK5Hj9YPE+rOeH/OwYqPY7bUfwLfazDThYRb2HyLPvtiHgts3OtzFnfsSJi\n",
"F/DmUsrHIuJ2Jvue97g1l4/J/p4GfpfJ34GvAX9EknwRcQNwqJTytunx/UyeyFPkm1r0z2OTOdd3\n",
"zHSpyt/D7Sj+VSbfrsBkgzvm261FRMSNwMFSyvumSzNzTccKreU9AFwVEZ8HrgeuYHIFmSJfKeVH\n",
"EfE94FWllCcj4gK5Hr8V4GcAIuIlTAohS76LH1WwUJ6Gcv74oxhmdAzUehy36QcZh4GvMimTw9v9\n",
"g5UlMzwBPAjcD3xiunbbrFwt5wXeAdyZLR/wi8CXp/s+kjDfB4GHgK8Dt7eeD/gA8ADwGPD3y+TZ\n",
"yTnn5HtBx9TK50c2SFJnfAOXJHXG4pekzlj8ktQZi1+SOmPxS1JnLH5J6ozFL0mdsfglqTP/Bwws\n",
"I9sPIpCkAAAAAElFTkSuQmCC\n"
],
"text/plain": [
"<matplotlib.figure.Figure at 0x117068b00>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"flu_unique.OI.hist(bins=20)"
]
},
{
"cell_type": "code",
"execution_count": 48,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"covariates = ['AgeYears', 'pH', 'OI', 'VA']"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Get counts of missing values"
]
},
{
"cell_type": "code",
"execution_count": 49,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"AgeYears 0.000000\n",
"pH 0.077323\n",
"OI 0.442379\n",
"VA 0.011152\n",
"dtype: float64"
]
},
"execution_count": 49,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"flu_unique[covariates].isnull().mean()"
]
},
{
"cell_type": "code",
"execution_count": 50,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"image/png": [
"iVBORw0KGgoAAAANSUhEUgAAAYQAAAEKCAYAAAASByJ7AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\n",
"AAALEgAACxIB0t1+/AAAIABJREFUeJztnXu8XlV557+/k0gSC4fES4GJ4JGLk48zrbS0wyDRhBzE\n",
"G7VaW4WqXOp8KkobpOVW4IOg6Ahi1epU+umMUKpcSqukRSohJ4RbEUocZToWDQmBgBNaMbfagESe\n",
"+WPvN9nnzX7fs2/vvj7fz+f9nLPWu9fzrL332vt517PWepbMDMdxHMcZq7oCjuM4Tj1wg+A4juMA\n",
"bhAcx3GcEDcIjuM4DuAGwXEcxwlxg+A4juMAbhAyIelKSS9IOrJgub8g6d8kHdqXf4ekq4rU5Th1\n",
"RNJlkp6S9LSkuyT9QuS7t0jaFD4j2yQ9IelVVda3bcjXIaRDkoDHgaeBu8zsnILlfxY4zMzeHqZ/\n",
"A/gT4NVm9u9F6nKcOiHpXOBM4Dgze0zSKcCnCNr+v0WOuwbYZGaXVFTV1uI9hPQcBxjwB8DJoYEA\n",
"QNJBkm6T9KSk/yvpTkk/lHRE5JhF4S/+JyV9W9JxffI/CvyKpBMkzQGuAv4wagwk7S/py+EvpEck\n",
"ndVfSUnnSPpHSY+Hnw/GHPOCpLdKui/8Vbay7/vXhb/Sngg/f5PxmjlOEj4EXGlmjwGY2XXAJuBd\n",
"ldaqQ7hBSM97gZvM7B7gZ8Bk5LvPAhvM7BXAhcBRwC8BjwJI+jngDmBFeMxpwE2SFvYEmNl24Dzg\n",
"88BFwGNmdlNfHa4D5gCHAccAp0nqf2i+Q/BL65Whni9IOjDmfC4FTjezhez94P0V8BdmdoiZHQKc\n",
"PuS6OE5mJI0DE8CDfV+tJXiGnBJwg5CC8Bf7uwhelAB/DbwvcsgvAX8X/n87sC+A7fHLnQjsMLMv\n",
"hvkPA18D3h/VY2ZfAf6FwDCc2VeHA4BfA842s+fNbAtwBfC7fTJWRbrZ9wDPAEewNx8ysx+EZXb0\n",
"fbcJeJ2k14Tfb48p7zhFME7Q897Sl78l/M4pgdlVV6BhvJ3gJb8y9BTNApB0hpk9C9wPvEfSHcBJ\n",
"wHozezpS/mDgEEkbwrQIfunfHKPrGmDMzB7pyz8YeAH4VlgHhfX4l+hBkk4kMCb7As8D+wEvitHT\n",
"/wBGeTNBN/5/SvoPwGfM7AtDjnecrOwgaMsL+vIXhN85JeAGIR3vAz5uZpf1MiR9H3gHcCNwLoGr\n",
"5ofA94Bf7yu/Afi2mb0hgS4LP/1sIHjBvyY0Qnsh6bUEBmUy7IUgaWMCndMrYLaNYFDvU5IOA+6T\n",
"9F0zuzutLMcZhpltk/Q48F+AhyJfHQVcXU2tuoe7jBIi6aXAm4C/7Pvqq+xx+XwZ+JSZHWhmy8zs\n",
"e33HfgN4maRzJc0O5c6VtE/SepjZjwmMz7WSFoQyxiTtGzlsIfATYIOkfSRdAhxEfA9hIJJeGUm+\n",
"iKC9xBohxymAq4FzetOuJZ1KMK7gkxlKwg1Cct4NPGRmG/ryvwJMSnoZ8APginCu9CZJGyTdKOnF\n",
"AGb2HMEg9GuBdeGv9geA16SsyxnAOgK30RPAI8BvRr7/JrCKoDfxvwnGD/6BwFBEmWnO8cfD2VCP\n",
"AzcB55pZ/6Cf4xTFlQTP072S/gX4b8AJvbGr3joEgra+3NchFI+vQygQSV8jmJWzIkzvC3yX4EX6\n",
"tUor5ziOMwOZegiSrg7n2K/pWWhJ10i6X9LqcEFJ79hJSfdIulvSsqIqXjck7Ufw67/ns38RQa9i\n",
"AcHUOafBSLo8bNt3eJt32kquHkK4qOq3zOzD4erBS8xsU+R7AfcSvCgF3J5wQLWRSHovcAHBjB4D\n",
"7iMYU/inSivmFIakY4H3m9kZ3uadtpF3ltEO4KeRdH+P4wjg+73ZMJIelXS4mT2aU28tMbOvEgwy\n",
"O+3laIIxmx6dbvNOu8hrED5AsKIWAuNwvaRnCBZNrQdeCmyT9McEv5a2hXn+cDiNQ9JdBLO1FodZ\n",
"3uadVpHZIIQLn77fWzhlZsvD/CMJ4u+8k2B2y3yCxU0CvhTmxTI1NeUj3E4hTE5Oauaj0mFmSyT9\n",
"KkHokDd7m3fqRBFtPpNBkHQUsHRApM9nCRZOQfCrqBcuQcCMXec8JyVpqZmtyVq+CBlevvp7MOKX\n",
"7NPsPV23sjafhCLuicuvt46i2nzWHsLNwCZJdwIPm9lZkm4k6E7vIIy/Y2YvSLqMYE68AZcNElgE\n",
"RVzwvDK8fPX3YBRIugl4GbAT+L0wr5A2v/qfnpi28Gp87qynfuXwhcuLqvuor6fLr4eOIshkEMzs\n",
"0Ji8kwYcewdBhE/HaSxm9p6YvELa/Pl//+hvRNPnLXnlg7+SuoaOk59WrVSWtLRqGV6++nvgTGfU\n",
"19Pl10NHEbTKIDiO4zjZaZVBqIP/2stXfw+c6TTdB990+WXpKIJWGQTHcRwnO60yCHXwX3v56u+B\n",
"M52m++CbLr8sHUXQKoPgOI7jZKdVBqEO/msvX/09GAUDop3GRjWtW7TTpvvgmy6/LB1F4FtoOk4C\n",
"zOxi2B3t9HxJHwI+RiSqKbA6jHa6V34llXaclNSuh/Dguqc+F/1897H/976kZevgv/by1d+DEdOL\n",
"dro7qqmZ7QQelXT4kPzKaLoPvunyy9JRBLXrIfzhrevOiqYvnpz4+9cG2+o5TqX0RTs9jPiopmMD\n",
"8hNHO/3Jju3jEZ1LYY/LIWP6SCBPeZefM91jVPJXrVpFEdRqC82pqSm74NvT43x9/IRDb3vzayfe\n",
"VlGVnAYyNTVlowoYF0Y7/Tjw+8CFTI9qejmBQfij/vxBAe7i2vx5S1754CEL5v7z5u3PTUTzDxyf\n",
"s/HoIxaeVuT5OO2gqDafNdrp1cB/JGjwp5vZY5ImgUsJAnpdamarw2Nj8x2nofSina4nJqqppLG4\n",
"/LRKNm9/buLy1RuXRPMuXjaRudKOk4RMYwhmdoaZHUcQyfHcyEDaG4E3ERgABuWPijr4r7189fdg\n",
"FEi6SdIU8KfA75nZCwTtfxXBwPFlEEQ7jcuvkqb74JsuvywdRVDUFpqx2wYSGBzfTtBpPAOincZG\n",
"NfUIv05TKWoLzUHbBuYeYEtDHebAe/nq74EznabP42+6/LJ0FEHmaaeavoVmb9vAC8PPgjBvUH5i\n",
"tvzoX18S0bk02vXytKeHpR3HSUdRW2jGbhtYxADbgpe9/Me9//utbJzVjebNdHxcOvpC8fLll4fg\n",
"xW5ma7KWn5qawtlD73q6/Grkl6WjCIrcQvNj9G0baCVvoek4juNkp8gtNFcCK2PySxtgq4P/2stX\n",
"fw+c6TTdB990+WXpKILaha5wnDoi6WpJd0paoz3B7a6RdL+CoHenRI6tVXA7x0lKqwxCHebAe/nq\n",
"78Eo6F97E/nq3Wa2zMyug/LX3iSh6fP4my6/LB1F0CqD4Dgl0Ft706P/GapdcDvHSUrtgtvloQ7+\n",
"ay9f/T0YMb21NxAYh+slPQOcbWbrGbwmp7LFmE33wTddflk6isB7CI6TkL61N5jZcjM7FrgEuCo8\n",
"LPfam2i005g61Gqth6frlc5Lq3oIRcz1zSvDy1d/D0ZBzNqbKM8Cz4f/x67JSaPr5/Yb3z7ou6rW\n",
"hrj8fOm8a2tmShe19qZVBsFxRkjc2psbCfZH2AGcCb72xmk2rTIIdfBfe/nq78EoGLD25qQBx9Yq\n",
"uF3TffBNl1+WjiLwMQTHcRwHaJlBqMMceC9f/T1wptP0efxNl1+WjiJolUFwHMdxstMqg1AH/7WX\n",
"r/4eONNpug++6fLL0lEEmQyCpMWSHpR0ZSTP47o4juM0mKw9hDnAJ/vyjIrjutTBf+3lq78HznSa\n",
"7oNvuvyydBRBJoNgZlPAlr5sxcjzuC5OK1B8tNPY3q/3ip2mUuQ6hJHEdQm20JwA9ljZ6Oq/aLqX\n",
"N+j7pOmoLC9ffvm86VWrVlE0ZnZGqOM44FxJZxL0ficJ2vbtwOpIr3hafuEVSkHTffBNl1+WjiIo\n",
"zCCY2XIASUcSxHV5J3viunyI4OH4EinjuqTZQtPTnoaRb6HZi3a6u/cLIKnX+x2Ly7eU4Sscpwry\n",
"zjJSTF6hcV1SVaYG/msvX/09GDEfAK4m0vuV9Fn29H4H5VdG033wTZdflo4iyNRDkHQ+8BbgAEnj\n",
"ZvZBj+vitB1Fop1KejXxvd+xAfmJ+cmO7eP7zN9/54A6LIXUbrQjgTTHu/yC3Zg96u4mlZkVIqgI\n",
"pqam7IJvT+90fPyEQ29782sn3lZRlZwGMjU1ZZOTk3G918woiHZ6soXRTiWNAXcDxxMYgZVmtnhQ\n",
"/rC69rf585a88sF9Zmnn5as3LonmX7xs4q5fP+rQpcWdldMWimrzrQpu5zgjJC7a6cfo6/16r9hp\n",
"Mq1aqVwH/7WXr/4ejAIzO9TMlpjZcWZ2Vpi30swWm9nrLYhw2jv2jrj8qmi6D77p8svSUQStMgiO\n",
"4zhOdlplEIqY65tXhpev/h4402n6PP6myy9LRxG0yiA4juM42WmVQaiD/9rLV38PnOk03QffdPll\n",
"6SiCVhkEx3EcJzutMgh18F97+ervgTOdpvvgmy6/LB1F0CqD4DijQr4HiNMBWmUQ6uC/9vLV34MR\n",
"Ucs9QJLQdB980+WXpaMIWmUQHGdUmO8B4nSAIrfQrHyzkDr4r7189fegRHp7gPytpMPCvNpFO226\n",
"D77p8svSUQSFbKE5qJtcx+6z4xSFmS03s2OBSwj2AIE9e4BcGH4WkCHa6aDvJC2Nuh887eki3VGZ\n",
"gtuZ2ZSkaCTGWmwWIu3ZLa0qGV6++nswYka+B8jP7Te+fdB3WTYNir4wRrEpkcufOd1r06OSX9Sm\n",
"UEVFOx20VebYgHzfPcppFPI9QJwOUJRBGLRVZu7NQnxP5e6Vz5se0Z7KVwBX9OWdNODYO4DKo5z2\n",
"aLoPvunyy9JRBHkNQq/7HNtNVrBZSK7us++p7Om06RHvqew4rSXrLKPzCbrCJ0r6MzN7gWDweBVw\n",
"O5HNQsL/p+WPiiIGV/LK8PLV3wNnOqO+ni6/HjqKIOugclz3eSWwMubYWnWfHcdxnHhatTCtCD9d\n",
"Xhlevvp74Eyn6T74pssvS0cRtMogOI7jONlplUGI89M9sO6pa1es3bAm+nlg3VPXppGRtw5evnwZ\n",
"zh6a7oNvuvyydBRBUdNOa8vm7c9NXL56Y3QRHRcvm6ioNk5TkbQY+GNgjZmdF+ZNEqy+N+BSM1s9\n",
"LN9x6k6rDEId/Ndevvp7MCJ64VpeB9PCskwSTKm+HVg9KL+KCvdoug++6fLL0lEErXIZOc6osL2j\n",
"nQ6KaurRTp3G0iqDUAf/tZev/h6UxKCoprWLdtp0H3zT5Zelowhq7zLad59Zv7hi7YY1/fkHjs/Z\n",
"ePQRC08rv0bpeGDdU9du3v7cRDSvKXV3hjKycC0/2bF9fJ/5+++M+y5jOI8jgTTHu/yCw6n0qHu4\n",
"ltobhB//+/Mv7h8UhviB4bL817Ev+f32ecXmHT998paH1rNi7Ybd+bPGtOjy1RsPiB47bFC7ah9+\n",
"1eWLkjFCRh6upehop4QvuxzlXX7N03WLdtop4mYuXbRsYtsnVm88rP/Yi5ZNbEsiM87IwGh6E95r\n",
"SY/io532wrXsjmpqHu3UaTCtMgjRSKdVyshCnJEBuOpthy/qd5kNe3knqf+wqbh5z7/J92AYTQ7X\n",
"Murr6fLroaMICjUIkq4BFgE7gWvN7DpJxwMfxedkZ2LLzl1zP9H38k5qJAb1OmaNadEIquo4TsMp\n",
"uodgwLvNbBPsnqt9GSXNye6A/xqINxK9X/jR+g/qdQxzY/kYQvto+jz+pssvS0cRFD3tVH0yfU62\n",
"4zhOQyi6h7ADuF7SM8AfMHhrzdbuqZyE+fNm7+XyyevG6cncvnXL/PH5C7ZmleljCO2j6T74pssv\n",
"S0cRFGoQzGw5gKQjgU8D55FzTvauXbti6zh/3uxFfzm19jsAvRfgl//61hdLOi86R/e6VQ/NHyQ7\n",
"65zkWx5an7j+P4upf5zLZ5AbJ658HNNlBqIueMMrfpK4ohQ2Jzv3nO6++lQ2J9txusaoZhk9CzwP\n",
"rCfnnOzZs2fvisvfsnPX3D/59rbXBqngBXjV237h6VseWn9pbx3ALQ+tZ9aYDux930/WOcnRdQYz\n",
"MWtA/cson6bs/HmzF93y0PpLAVas3fAOgG/94MlpA9VNmdPtW2hOp+k++KbLL0tHERQ9y+hG4CAC\n",
"19GZZc/JTvPL25nOsIFqZzBlzqyLczX6+hGnSIp2GZ0Uk1fanOyk7pVBPLDuqWt/8MTmI3suKNiz\n",
"Ajl6nE/bHEwHxxBKm1mX1Wg33QffdPll6SiCVi1MS0rcLy0IXvR/8u1tB0RdTHErkL3X4UQYOLMO\n",
"QNKjklK7Sh2nClplEJL6zON+aYG/6Iugg+sQKp1Zl4Sm++CbLr8sHUXQqvDXjlM2ZrbczI4FLiGY\n",
"Wfcjgpl1F4afBWSIdpr02O1bt8yPzs6StNTT3U3npVU9hLxjCE5+OjiG0KOwmXXDop32Mz5/wdbo\n",
"tYqbeRV9YYxiZpfLnznda9N1n1nnL1BnID6rZWaqnlnnOEXSKoOQd86/M50ss1q6NoZQ9cy6JDTd\n",
"B990+WXpKAIfQ3Acx3GAlvUQfAxh9AyasttzJXV4DKG2NH0ef9Pll6WjCPwF6qRi0JRdX9XsOM2n\n",
"VS4jH0Oonq6NITSBpvvgmy6/LB1F0CqD4DiO42SnVS4jH0OoHh9DKJckU4Ob7oNvuvyydBRBaS9Q\n",
"SZPApfjeyk5HKKPNe5Rap0hKMQhhBMiPMeK9lX0MoXp8DCGgrDafhKb74JsuvywdRVBWD8EjQLac\n",
"WNdFTOjw2Lx2rn6urM37CnMnK2UZhMQRIM9efPD/iaZnjemApEp8DKE6Bm1OFBc6vD/vqrcdPu0F\n",
"1tsXuuEvscxtfmLB3J0/3P5cZsVx96L/GkMKgz1D3u77lUMexBut2D1KBhy3eftzE1l0v7Dt6Wff\n",
"ueyYN/fXp0i+vvr+b47tf8Dc/vy6tXGZ2eiVSK8G/ojpeytf3v9raWpqavSVcTrB5OSkqtTvbd4p\n",
"myLafFkGYQy4GzieYKrrSjNbPHLFjlMR3uadJlKKi8UjQDpdw9u800RK6SE4juM49cdXKjuO4ziA\n",
"GwTHcRwnpDYGQdKkpHsk3S1pWYLjF0t6UNKVWWRIulrSnZLWSHpVxjpcLmm1pDuyygjL7CNpo6QP\n",
"h+njk8iQdI2k+8M6nJKmbETGwrD83ZI+k1L/eHgNV4d/t6atg6TTJT0g6V5Jx2U8h98Nr8NKSYeH\n",
"eanvQ9kUVcc0bTmPzqTtNOMzkLgdZpSfuJ0llZ/mHZTlfgyQv9e9zio/FjOr/EMwLe8+YC4wD7g7\n",
"QZlJ4B3AlVllhOWOA/40a/lQxrHA1TnqsBz4G+DDaWQAXwYOznkdbwCOySMjLPeLwJ+nLQ88TPDD\n",
"ZBz4hwzl5wHfCv9/KXBznntZ5zafty3n1ZmkneZoP4naYQ75idpZyucv0Tsoq65++TH3+ktFX6u6\n",
"9BB2r+o0s53Ao71feoMwsylgSx4ZITuAn+YoD3A08EgWGZLmAScAKzKch5jey0ulX8HUyMPN7P6s\n",
"MiL8PvCFDOUfJpia+XaC8A5pywuYLWkfgsVfB+Q4hzIZRR1nasuZdaZop1megTTtMOs5JG1nieWn\n",
"eAdl0hUjP8oOoLd6sbBrVZeVvYlXdY5AxgeAz2ctL+kugk3WFwOHZZCxHPgiwYss7XnsAK6X9Azw\n",
"BxnO4eXAXElfJ/jl9EVgc9pzkPQSgp7Kw5KOSVn+HuBUAsN2Q9pzMLN/l/RJ4Jvh9VhAcC3ztqdR\n",
"U0Sb72emtjyWQ2fSdppFR5p2mPUckrazPNcorcw8unr3OovegfLrYhCeAeYzfVXnM6OWIelEAgv6\n",
"iIKVpanrYGZLJP0qcB3Br+TEMiSNA683sysknRqWSXweZrY8lHMk8GngvJTn8AywFXgXQVu4D/id\n",
"lDIAfpfAXdSTmai8pMOAZWb2njB9J8GLJ5V+M/sa8LVQxlqCl0ne9jRqimjzu0nYlsey6EzZTrPo\n",
"SNMOU8tP2c4yXaPIeaSRmfV+7L7XGfUOpC4G4VGC7g0EFU8TBKy3XDuVDElHAUvN7JwC6vA0weKj\n",
"9SllLAbmSLoeOBSYRfBLJm09ngWeT6vfzHZJ2gQcZGZPSXqW9NdxFvBrwOvDrDTlx4D9QzkvImi8\n",
"me+DpLcC3yX9faiCPO1tGknbcuiayaIzcTvNoiNNO8x4DonbWUb5Q99Bg2Sm0LU7JEXMvU6td+iZ\n",
"zDTIUNYHeCNwL0FDe2OC488H1gD/DPxZmHdCUhnABuAu4E7g82nLh8ffBEwBtwKHZTmPiKxTgA+n\n",
"qQdwY3gOtwKvzHgdDwFuC8ssz3AdfxO4IOu9BC4E7ge+BZya8T78r/D4bwIvy3Mf6tzmi2jLeXUm\n",
"aadZdKRphxnlJ25nSeWT4h00SOYwXQPk73Wvs8qP+/hKZcdxHAeo0ToEx3Ecp1rcINQMSWeHg179\n",
"+UdI2imp5wsdk/RDSd8pv5aOUy6SHutbcHWqpHuqrFMbcYNQP24AjpF0UF/+ycA3zGxbmJ4EdgIH\n",
"S/rPZVbQcWqC+7sLxg1CzTCzzQQDSSf1fXUy8JVI+n3A9cDXgfeXUjnHGRGSPirpZklfCHsDGySd\n",
"PqxIaZXrEG4Q6slXgd/uJcKpZj8PfCNMzwPeSdCbuIm9jYfjNJFlwI1m9irgvcD/kPSfIt9/JTQU\n",
"jxGsu3EKxg1CPfkasCiyzPxk4GYzez5MvxN4wsy+RzD9bK6kpeVX03EK5Rtmdh+ABWEsbiVYrNbj\n",
"vWZ2aGgwzokT4OTDDUINMbOfAH8LnCxJwHuY7i56L/AqST8G/pVguf/7Sq+o44yWx9kTKgPcTTRy\n",
"3CDUl68QuI2WAM+b2b0Akl5OEKTrNWb2EjNbQDDA/K4wwJvjNJX+yAmHAU9UUZGu4gahvqwkCNT2\n",
"KYLB4x4nAw+Y2eO9DDP7B+DHBJEcHaepvFPS2wEkTQJvIhgnc0rCDUJNMbOfEQwY/yrT3UW/DfxF\n",
"TJGv4G4jp9ncCvx6GNfoGuADZtbrIfgU0xLw0BWO41SOpI8SxAM7peq6dJmhPQSNeIs4x2kKKmmb\n",
"SsepkqE9hNCPtx/wOjM7L5zxci/BIKaA283sDWnzR3pGjjNCFOzH+1vAmXibLwzvIdSDofshmNmU\n",
"pCWRrN1bsgFI6m3JNpYm3+oXm95xkrLXNpXgbT4vZnZZ1XVw0m+QU+YWcY5TR0a5TaXjVEpagzDS\n",
"LeKmpqZ8hNsphMnJycIXMY1im0pv805RFNHmkxqEUW8Rt5tRPMjDkLTUzNa4zvboHcVLdpTbVBbV\n",
"5ou61i6nPFlFySmqzQ81CJLOB94CHCBp3Mw+KOljwCqCecGXAZjZC5IuS5pfJ6p4SXZFZ5V6R8DN\n",
"wKZwr4qHzeysNM9CGRR1rV1OebLq9nzMNKh8BXBFX95KglW0/cfeAdyRNN9xmoSZHRqTl+pZcJy6\n",
"U7uVyvc+sumr0c/933/yhi1btkyMSl8VUUK7orNKvV2kqGvtcsqTVbfnI+2g8sg5++/W/XY0fdTC\n",
"/Xb997cc/gVgYzU1chzH6Qa16yGUTVf8+T6G0H7q5tduq5wiZdXt+ei8QXAcx3ECaucyKpuuTAHt\n",
"0rTTrlK3qZBtkPPAuqeu3bz9uYlo3oHjczb+11e/4to6nVtRdN4gOI7TPQa96I8+YuFp0bzN25+b\n",
"uHz1xmj4Hi5eNq1Yq+i8QeiKP9/HENpP3fzadZaT90Vft3MrCh9DcBzHcQA3CJ1ZE+DrENpP3ebG\n",
"t1VOkbLq9nx03iA4juM4AZ03CF3x5/sYQvupm1+7rXKKlFW35yOzQZB0uqQHJN0b7iKFpON9S0HH\n",
"cZxmkqeHcDZwDPBW4BPh1oGXAW8E3gRcChDmf6w/vy50xZ/vYwjtp25+7bbKKVJW3Z6PPNNOHwaO\n",
"B34euB3fUtBxHKfR5DEI9wCnErzwb6ChWwp2xZ/vYwjtp25+7bbKKVJW3Z6PTAZB0mHAMjN7T5i+\n",
"E1hOAVsK9mPY7t2ket2r3kX0tKfj0qtWrcJxnPRk7SGMAfsDSHoRwQu/sC0Fowjt3hqu35oWkY76\n",
"8EYhPy7di19Slr6ozrL09dJVXN+pqSmKRtJi4I+BNWZ2Xph3DbAI2Alca2bXhfmTBGNlBlxqZqsL\n",
"r1B8HRsfO6gJcoqU1YpYRma2LpwxdD/BS/5zZrazblsKOk6BzAE+CbwukmfAu81sUy8jMolikuDZ\n",
"uB0oxSA4Tl4yjyGY2ScJHpBoXuO2FOyKP9/HEPJhZlOSlvRli71n6sVOrihjEkXd/Np1kdMLZHfL\n",
"Q+tZsXYDALPGtKj/uPnzZi9asXbDNF1xxxVRp6LlFEXng9s5Tg52ANdLegY428zWM3hyRS0mUXSR\n",
"uEB2Fy2b2NZ/3Jadu+Z+IsFxbabzK5W7sibA1yEUj5ktN7NjgUuAq8LsZwjG1C4MPwuYYRJF9BpJ\n",
"Wpo13fu/AHkfaVN9tm/dMp+M/GzXrr1+NG/fumV+/znmqN9Hirr/ReA9BMdJh2LyngWeD/+PnVwx\n",
"TGDUbVDFwH9/OhgGaU99VqzdsDXoqKVn1uzZu/rzxucv2JqnPn3ivlPE/S9qIkXnDUJX/Pk+hpAP\n",
"SecDbwEOkDRuZh+UdCNwEIHr6EyodhJF3fzadZNTJG09t84bBMdJgpldAVzRl3fSgGNrO4nCcYbh\n",
"Ywgd8ef7GEL7Kepat1VOkbT13DpvEBzHcZyAzruMuuLP9zGE9lM3v3YVcnprDqJ5g9YSlFWnMuQU\n",
"RecNguM47SHpmgMnns67jLriz/cxhPZTN7923eQUSVvPrfMGwXEcxwnovEHoij/fxxDaT9382nWT\n",
"UyRtPbc8eyovlLQ6jHr6mTDP91R2HMdpKHl6CFcBF5nZG8zsD8Owv76nsuushd4uUje/dt3kFElb\n",
"zy3rjmljBDFa7o9k+57KjuM4DSbrtNOXA3MlfR0YB74IbMb3VHadNdHbRerm166bnCJp67llNQjP\n",
"AFuBd4Uy7gN+B99T2dM1SPueyo6TjaxbaO6StAk4yMyekvQsvqdy4rTkeyqPMj2KPZWbQNw9djmj\n",
"oa3nlmdnM3IIAAALRUlEQVSl8gXAn0saB/7KfE9lx3GcRpNnT+UngLf25fmeyq6zFnq7SN382nWT\n",
"UyRtPbfOL0xzHMdxAjpvELqyJsDXIeRD0mJJD0q6MpIXu+CyqoWYdZsbXzc5RdLWc/Nop46TjDnA\n",
"J4HXwbQFl5MEkyVuB1YPyq+iwo6Tls73ELriz/cxhHyY2RSwJZK1eyGmme0EegsxB+WXUcc1Lqcc\n",
"2npu3kNwnGy8lAYuxHTyM3/e7EUr1m5YE807cHzOxqOPWHhaNTUqjs4bhCrmAXdFZ5V6S+AZClqI\n",
"Gb1GeRbnxa35yCjvSDP7XBPrs33rlvn08bNduzK/5+LK/mjHs/t+6u4np23Cc/GyCVKcTzT9EeA7\n",
"dVmM2XmD4Dgp6a2cL2whZtRgVrF4sD8dDIM0sz7j8xdsDTple5g1e/YuMhJXdpi8tOdHxBhkKV/0\n",
"YszOG4Su+PN9DCEfks4H3gIcIGnczD5Yt4WYdfNr101OHanbuXXeIDhOEszsCuCKvrzGLcR0nGF0\n",
"fpZRV9YE+DqE9lO3ufF1k1NH6nZunTcIjuM4TkAul5GkfYAfAFea2Z9KOh74KIHv9FIzWx0eN0mw\n",
"U9q0/DrQFX++jyG0n7r57Ect54F1T127eftzE9G8WWNaVITOsqjb85F3DOEMYC3sXrl5Gb5y03Gc\n",
"Eti8/bmJy1dvnDb986JlE9sGHe/MTGaXkaR5wAnAijCrdis3k9AVf76PIbSfuvns6yanjtTt3PL0\n",
"EJYTbJ15QJj2lZuO4zgNJlMPIdwU5/Vm9s1eFntWbl4YfhaEeYPyE9G/hWbUohaR7juvwuXHpaOr\n",
"DMvQB3tWjpalr5cmQtn6u0ZTfP9VyakjdTu3rD2ExcAcSdcDhwKzgHto4Baanm5fuqtbaDpOXrLu\n",
"qXwbcBuApFOAfc3s4bqt3ExC9Be762yW3rhZJgD7jkJZAyjqWrdVTh2p27nlXqlsZtdF/veVm05p\n",
"xM0yAfjUL1dRG8dpPp1fmNaVNQG+DqH91M1nXzc5daRu59Z5g+A4juMEdN4gdGVNgK9DaD91m/df\n",
"Nzl1pG7n1nmD4DiO4wR03iB0xZ/vYwjtp24++7rJqSN1OzffD8FxnNrThkB2TaDzPYSu+PN9DKH9\n",
"1M1nX6Sc3hTj6GfXCza3CPlVUrfnw3sIjpMDSdcAi4CdwLVmdt2gMPBOe5k/b/aiFWs3rInmHTg+\n",
"Z+PRRyw8rZoaZaPzBqEr/nwfQxgZBrzbzDbB4DDwpVSkZj77IuWsWLuhCFEjY8vOXXM/0bdI8uJl\n",
"EzOWq9vz0XmXkePkREx/jmod7t1xhtF5g9AVf76PIYyMHcD1kv42fPHvDgMv6bPsCfc+curo+6+T\n",
"nDpSt3PLGv76akl3Sloj6VVh3qSkeyTdLWlZ5NjYfMdpA2a23MyOBS4BPg38iJTh3ssOSZ4gZPmR\n",
"da3PTPxs16693OBxeWXJSxCi/cg6hXzPGu30jLAyxwHnSjqTmG0yQ39qrbfP7Io/38cQRs6zwPPA\n",
"elKGe49eozwhwOOudUZ5tavPirUb3tEvK45Zs2fvSpKXlLzyEpzf51IeH5suKuR73kHlHcBPifhN\n",
"AST1/KZjcflp9kNwnDoj6UbgIIJn4cy6h3t3nGHkHUP4AHA1g/2mlflTk9IVf76PIYwGMzvJzJaY\n",
"2Ylm9niYd4eZLTaz14eh30uhbj77usmpI3U7t8w9BEknEvz6f0TSqwn8ph8i6CZ/icBvOjYgPzH9\n",
"W2jCnm5SQekjCbvII5K/V3rE51O39Miu7/atW+bjOE5hZDIIko4ClprZOWHWo4xg+8yg0Mi30FyT\n",
"8vhGpgv056ZND6xD3vT4/AVbg06nA/VcP1CUnLqvQ8hK3cbYsvYQbgY2SboTeNjMzlIDt890HKd+\n",
"eNyi6sg6y+jQmLxGbp+pjuxvXIXOKvV2kaKuddVy4rZGvWjZRCu7gnV7Pjq/MM1xHMcJ6LxB6Mqa\n",
"gKp+hdTp10/bqaPvvwg5baZu16jzBsFxHMcJ6LxB6MqaAF+H0H7qNu/f7/3M1O0adT78teM41eEz\n",
"iupF5w1CV/z5PobQfurm+08ip0sziuKo2/PReZeR4ziOE9B5g9AVf76PIbSfuvn+/d7PTN2uUecN\n",
"guM4jhPgYwgd8ef7GEL7qfsYQtcGkOfPm71oxdoNa6J5B47P2Xj0EQtP66Xr9nx03iA4jlM8g17+\n",
"l6/eeEA0r80DyFt27pr7ib4B84uXTVRUm2R03iB0Ja6QxzJqP1XHIIrS9dlDSanb81HaGIJ8b2Wn\n",
"Y3ibd5pGKT2EOu+t3BV/vo8hlEsVbb6qMYSujQ3koX9c4ZaH1vPAuqemjStUSVkuo9g9l9NuluM4\n",
"DaLxbT7uRd8/KAruHkpD3LjCVW87fMbB57IoyyDs3luZ4NdSb2/lvR6O9//ygU9G0wfvP3efUVas\n",
"K/58H0MoncRtfhjrnnz6G/vPnX18NO/ZXS88cciBLz+i/9ik13qmF31PTtyLPu7l5b2BfNRp8Flm\n",
"NvNReZUEey7/EdP3Vr68/9fS1NTU6CvjdILJyUnNfNTo8DbvlE0Rbb4sgzAG3A0cTzCQvdLMFo9c\n",
"seNUhLd5p4mU4jLyvZWdruFt3mkipfQQHMdxnPrjsYwcx3EcwA2C4ziOE1K6QUizerOolZ4pdV4t\n",
"6U5JayS9qgyd4fH7SNoo6cNl6JS0UNLq8NjPlKTzdEkPSLpX0nE5dC6W9KCkK4us36hIc62Htb+0\n",
"92xYm0pZp4HHppQTe/8ljYfnvDr8u3UGObH3NIOc2GudVk5YJvZaZ6hT7PXMICfbs2ZmpX0Ipt/d\n",
"B8wF5gF3F3FsUTr7yh0HfKksncBy4G+AD5ehE7gBOKas+xke/zDBj5Bx4P4ceieBdwBXjuLeF/3J\n",
"cq3j2l9aOcPaVBpZw45NKWfG+w/8IvDnee/pTHJmutZp5SR5fpPISnI9E8rJ9KyV3UPYvXrTzHYC\n",
"j0o6vIBji9IZZQfwXAZ9qXVKmgecAKzIqC+VTgVTIg83s/tz6EulM+RhgmmYbwe+mVWpmU0BW0ZQ\n",
"v8LJca2ntb+0coa1qTSyhh2b4dyS3P/fB74wREbSezqTnCjDnvUZ5aR4fofKSnE9k5xbpmet7Gin\n",
"aVZvFrLSM4ecDwCfT6krq87lwBeBAwZ8X7TOlwNzJX2d4BfEF83s6yPWCXAPcCrBL5cbMugbdf1G\n",
"QdZr3d/+0soZ1qbSyBp2bNo6Db3/kl4CHGxmDw+RMeM9TSgnSuyznkLOjM9vQlkzXs8Udcr2rCXt\n",
"ShTxAV4NXMOe7t61BBYx17FF6YyUORH4SEnnOQ7cGv5/KnBmCTpnhw1mDNgH+Edg3oh1HgbcFEnf\n",
"Cbw4xzVewswuo0LaUJ5Plmsd1/7SyJmpTaWUNfDYlHJmvP/ABcC78t7TJHKGXeuU9Un0/CaUNeP1\n",
"TCgn87NWdg/hUYIuHwTWfViwrzTHFqUTSUcBS83snAy6suhcDMyRdD1wKDBL0p1m9r1R6TSzXZI2\n",
"AQeZ2VOSnk2pK7VOgka+P4CkFwHzgRcy6u0x01L9otpQZtJe60HtL6WcoW0qjaxhx6as09D7L2kW\n",
"8GvA64fIgBnuaQo5Q5/1FHJmfH6TyprpeqaoU/ZnLYnVKPIDvBG4l8ASvjGS/1vAW5McO2KdG4C7\n",
"CKzq58vQGfnuFDIOKmc4z0OA28LjzypJ54XA/cC3gFNz6DwfWAP8M/BnZbShPJ9B1zpt+0sjZ6Y2\n",
"lbJOA9tKSjkD7z/wm8AFMfVMdU9Tyhl2rRPLSXCt09Rp2LVOIyfTs+YrlR3HcRzAF6Y5juM4IW4Q\n",
"HMdxHMANguM4jhPiBsFxHMcB3CA4juM4IW4QHMdxHMANguM4jhPiBsFxHMcB4P8Dv//p3epcet4A\n",
"AAAASUVORK5CYII=\n"
],
"text/plain": [
"<matplotlib.figure.Figure at 0x1171529e8>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"flu_unique[covariates].hist(bins=25);"
]
},
{
"cell_type": "code",
"execution_count": 51,
"metadata": {
"collapsed": false,
"scrolled": true
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x116eb6f28>"
]
},
"execution_count": 51,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": [
"iVBORw0KGgoAAAANSUhEUgAAAXUAAAEBCAYAAACHTjUfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\n",
"AAALEgAACxIB0t1+/AAAE+lJREFUeJzt3VFsnfV5x/HfEzIC25IaUalhsMUQUmUX09BQw1TSxfWh\n",
"1SZF6nqTdZrKkHrRFjbaTc2K2mgbCCrBEGunaaM3A3GBqKZJTdUbAieE0Im1NKjqzWCEKBVUSqWx\n",
"xLZoEpr22cU5Rqe2j33O+T/H533+/n6kSDmv7dff/J08Pv779RtzdwEA6rBp0gEAgDgMdQCoCEMd\n",
"ACrCUAeAijDUAaAiDHUAqMiqQ93M7jezo2b2jJld3z3WMrMXzOy4mc2uTyYAYBA2yHXqZnarpE9K\n",
"+qyk70hqSTJJT7v7H4y1EAAwsEG3X26R9IqkXZJedfcL7n5e0kkzu3FsdQCAoWxe6xXM7HlJ10ja\n",
"K2mnpDkze0SdZ+pzkq6WdHKckQCAwaz5TN3d90n6M0lPSPpfSVOSvtT9dZWkt8YZCAAY3JrP1Lt+\n",
"Isklva7OFozUeaZ+o7v3fZbebre5sQwAjKDVatkob7fqUDezb0h6r6Tzkv7C3X9hZvdKeladIX/v\n",
"uMIAYKMqeUK86vaLu/+Ju7fcfb+7v9499oy773X3D7n7M6O+4yYxs5lJNwyCzlh0xsrQmaGxFD98\n",
"BAAVYahLcvdjk24YBJ2x6IyVoTNDYymGOgBUhKGuPPtsdMaiM1aGzgyNpRjqAFARhrry7LPRGYvO\n",
"WBk6MzSWYqgDQEUY6sqzz0ZnLDpjZejM0FiKoQ4AFWGoK88+G52x6IyVoTNDYymGOgBUhKGuPPts\n",
"dMaiM1aGzgyNpRjqAFARhrry7LPRGYvOWBk6MzSWYqgDQEUY6sqzz0ZnLDpjZejM0FiKoQ4AFWGo\n",
"K88+G52x6IyVoTNDYymGOgBUhKGuPPtsdMaiM1aGzgyNpRjqAFARhrry7LPRGYvOWBk6MzSWYqgD\n",
"QEUY6sqzz0ZnLDpjZejM0FiKoQ4AFWGoK88+G52x6IyVoTNDY6nN434HT7742g+jzrXrvb/6fx/Y\n",
"de1M1PkAoDZjH+r/+J03fifqXAf37XjpA1En62FmMxk+g9MZi85YGTozNJZi+wUAKrLmUDezR83s\n",
"OTM7ZmbXd489ZmYvmtlRM7t9/JnjleUzN52x6IyVoTNDY6k1t1/c/TOSZGYflnRQ0p3dFx1w9zfG\n",
"2AYAGNIw2y8Lkt4Z8W0bLcu1q3TGojNWhs4MjaWGGcyfkvRo9/cLkp40s2+Z2c74LADAKAa6+sXM\n",
"9kt61d1fkSR3v7t7/CZJD0v6+NgK10GWfTY6Y9EZK0NnhsZSg3yj9GZJM+7+1RVefEHSz8Kr+nh7\n",
"YW5r75dPZjbDYx7zmMc1Ph6Vufvqr2B2StIbkn4h6Yfu/jkze0rSNepsw9zl7j9a6W3b7bbf87KV\n",
"Nr7r4L4dLx3Ys3NP2Am7zHJcu0pnLDpjZejM0Ch1Zmer1RppeA5y9csNKxz7xCjvDAAwXtVcwVIi\n",
"w2duic5odMbK0JmhsRRDHQAqwlBXnmtX6YxFZ6wMnRkaSzHUAaAiDHXl2WejMxadsTJ0ZmgsxVAH\n",
"gIow1JVnn43OWHTGytCZobEUQx0AKsJQV559Njpj0RkrQ2eGxlIMdQCoCENdefbZ6IxFZ6wMnRka\n",
"SzHUAaAiDHXl2WejMxadsTJ0ZmgsxVAHgIow1JVnn43OWHTGytCZobEUQx0AKsJQV559Njpj0Rkr\n",
"Q2eGxlIMdQCoCENdefbZ6IxFZ6wMnRkaSzHUAaAiDHXl2WejMxadsTJ0ZmgsxVAHgIow1JVnn43O\n",
"WHTGytCZobEUQx0AKsJQV559Njpj0RkrQ2eGxlIMdQCoCENdefbZ6IxFZ6wMnRkaSzHUAaAiDHXl\n",
"2WejMxadsTJ0ZmgsxVAHgIqsOtTN7FEze87MjpnZ9d1jLTN7wcyOm9ns+mSOV5Z9Njpj0RkrQ2eG\n",
"xlKbV3uhu39Gkszsw5IOmtldku6T1JJkkp6WdHTckQCAwQy6/bIg6R1JuyS96u4X3P28pJNmduPY\n",
"6tZJln02OmPRGStDZ4bGUqs+U+/xKUlfk3S1pDkze0SdZ+pz3WMnx5MHABjGms/UzWy/Os/OX5H0\n",
"lqQpSV/q/rqqe2xdvL0wt7V3T8zMZiIeLx6LOt8YH3++YT2sJ+vZ+zjDen6+YT19H4/K3L3/C81u\n",
"lvSn7v6F7uNNko5Luk2dTwhH3H1vv7dvt9t+z8tW2viug/t2vHRgz849YSfsMrOZDF+W0RmLzlgZ\n",
"OjM0Sp3Z2Wq1Rhqea22//LukN8zsOUk/dPfPmdl9kp6V5JLuHeWdNk2GD7JEZzQ6Y2XozNBYaq2r\n",
"X25Y4dgRSUfGVgQAGBk/fKQ8167SGYvOWBk6MzSWYqgDQEUY6sqzz0ZnLDpjZejM0FiKoQ4AFWGo\n",
"K88+G52x6IyVoTNDYymGOgBUhKGuPPtsdMaiM1aGzgyNpRjqAFARhrry7LPRGYvOWBk6MzSWYqgD\n",
"QEUY6sqzz0ZnLDpjZejM0Fhq0PupY0Dffe3Hj5+Zvzgddb7t27acvmXXtXdEnQ9A3Rjqir0d55n5\n",
"i9P3Hz29L+JcknRodvrd32e5bSidseiMk6GxFNsvAFARhrry7LPRGYvOWBk6MzSWYqgDQEUY6spz\n",
"7SqdseiMlaEzQ2MphjoAVIShrjz7bHTGojNWhs4MjaUY6gBQEYa68uyz0RmLzlgZOjM0lmKoA0BF\n",
"GOrKs89GZyw6Y2XozNBYiqEOABVhqCvPPhudseiMlaEzQ2MphjoAVIShrjz7bHTGojNWhs4MjaUY\n",
"6gBQEYa68uyz0RmLzlgZOjM0llpzqJvZXjP7npk91HPsMTN70cyOmtnt400EAAxqkP/5aIukr0j6\n",
"YM8xl3TA3d8YS9U6y7LPRmcsOmNl6MzQWGrNZ+ru3pZ0dslhG+RtAQDra9TBvCDpSTP7lpntjAya\n",
"hCz7bHTGojNWhs4MjaVGGurufre73yrpbyU9HJvU39sLc1t7PyhmNtO0x/Pnzk6V/BmXWnK+myb9\n",
"56vpsVjPDbeekm5qUs8a6zkSc/e1X8lsn6T97n5wyfHdku5z9wMrvV273fZ7XrbSxncd3LfjpQN7\n",
"du4JO+EYHD5x6tj9R0/vizrfodnp5z928w0zUecD0HztdttbrdZIw3PNb5Sa2Rcl/ZGk95nZNnf/\n",
"tJk9JekadbZh7hrlHQMA4g3yjdIH3X3G3X/b3T/dPfYJd9/n7vvd/UfjzxyviC951gOdseiMlaEz\n",
"Q2MprmABgIow1JXn2lU6Y9EZK0NnhsZSDHUAqAhDXXn22eiMRWesDJ0ZGksx1AGgIgx15dlnozMW\n",
"nbEydGZoLMVQB4CKMNSVZ5+Nzlh0xsrQmaGxFEMdACrCUFeefTY6Y9EZK0NnhsZSDHUAqMgg//NR\n",
"9cxsJsNn8CZ2fve1Hz9+Zv7idO+x+XNnp7ZNXXVulPNt37bl9C27rr0jom0tTVzPldAZJ0NjKYY6\n",
"ipyZvzi98q2G50Y636HZ6bIgYINj+0V59tmydGaRZT3pjJOhsRRDHQAqwlBXnmtXs3RmkWU96YyT\n",
"obEUQx0AKsJQV559tiydWWRZTzrjZGgsxVAHgIow1JVnny1LZxZZ1pPOOBkaSzHUAaAiDHXl2WfL\n",
"0plFlvWkM06GxlIMdQCoCENdefbZsnRmkWU96YyTobEUQx0AKsJQV559tiydWWRZTzrjZGgsxV0a\n",
"G27qys27D584dSzqfNu3Xn7dmYV33ow632WbbHfUuQCUY6ir2fdYPnv+0hUPrHhr29F8eXZ67oGj\n",
"p3dGni/qXOutyR/3XnTGydBYiu0XAKjIqkPdzPaa2ffM7KGeYy0ze8HMjpvZ7PgTx6/2z9xYWZaP\n",
"O51xMjSWWmv7ZYukr0j6oCSZmUm6T1JLkkl6WtLRcQYCAAa36jN1d29LOttzaJekV939grufl3TS\n",
"zG4cZ+B62AjXrmK5LB93OuNkaCw17DdKr5Y0Z2aPqPNMfa577GR0GABgeMN+o/QtSVOSvtT9dVX3\n",
"2Lp4e2Fua+9nWjObiXi8uM8Wcb75c2enSv+cvX5+6VLoFUpNP1/v+kV9fPs9Xjy2Xu+vCX8/N/p6\n",
"Lm2ddM9a6zkKc/fVX8Fsn6T97n7QzDZJOi7pNnU+IRxx97393rbdbvs9L1tp47sO7tvx0oE9O/eE\n",
"nXAMDp84dez++EsQ37NRzndodvr5j918w0zU+YCM2u22t1qtkYbnWle/fFHSvZL2m9nX3f0X6nyj\n",
"9Fl1vkl67yjvtGk2wj4blsvycaczTobGUqt+6ezuD0p6cMmxI5KOjDMKADAafvhIG+PaVSyX5eNO\n",
"Z5wMjaUY6gBQEYa6NsY+G5bL8nGnM06GxlIMdQCoSLq7NH73tR8/fmb+4nTU+bZv23La3e+IOh/y\n",
"yLK/SmecDI2l0g31M/MXpyOvAz80Ox11KgCYOLZftDH22bBclo87nXEyNJZiqANARRjq2hj7bFgu\n",
"y8edzjgZGksx1AGgIgx1bYx9NiyX5eNOZ5wMjaXSXf0SberKzbufePb7Xz184tS5iPNdtsl2R5wH\n",
"GNU4Lvu9Zde1d0SdD+O14Yf62fOXrvinl+d+t/P/fZT78ux0zIkwdln2V4ftnNRlvxnWM0NjKbZf\n",
"AKAiDHVsWFn2V+mMk6GxFEMdACrCUMeGlWV/lc44GRpLMdQBoCIMdWxYWfZX6YyTobEUQx0AKsJQ\n",
"x4aVZX+VzjgZGksx1AGgIgx1bFhZ9lfpjJOhsRRDHQAqwlDHhpVlf5XOOBkaSzHUAaAiG/4ujWiW\n",
"qSs37z584tSxyHNu33r5dWcW3nlz6fH5c2entk1dNfQtl9f7VrRmNpPhGWaGzgyNpRjqaJSz5y9d\n",
"8UDgbWOlzu2QHzh6eufKLx3+TsmD3ooWmAS2X4CGy/LMMkNnhsZSIw91M3vMzF40s6NmdntkFABg\n",
"NCXP1F3SAXefdfcnooIA/LIs11Zn6MzQWKpkqFvh2wMAgpUM5QVJT5rZt8yszzehAJTKsg+coTND\n",
"Y6mRr35x97slycxukvSwpI9HRQEARhNxSeMFST8LOM+a3l6Y23pRfnnkOX9+6VLoZZ2cr1nnG4fF\n",
"a+nnz52dkqTFa91Hffz+39r+g1t2XXvH4n7v4rPJ3v1fdz/W7+VLH3/z+6+P44/d9/31PP68pB+s\n",
"1Rf9+L/+5807zsxfnB5kvc//9O1ff99vXPdmv5cvPt6+bcvp33//dY+vR3+/x6Ma+R+QmT0l6Rp1\n",
"tmHuKokY1K9tfc/Clsvsp6NcW9zPZZs3Xwo7Gedr3PnGYfm19Ev/Pg73+FB3qCz9x7zWP/Z+jw+f\n",
"OLVa/sjWev/qGeiDvH7U48MnTv39/cN8PF5Z+jMLy1//0Oz0uvX3ezyqku2XT0QEAFhdln3gLJ21\n",
"4+oVAKgIQx1ouCzXVmfprB1DHQAqwlAHGi7LXnWWzto1/vIxAJMVfTvk9b518UbDUAcabtL3AI++\n",
"HTK3Lh4vtl8AoCIMdaDh2KvGMBjqAFARhjrQcFz/jWEw1AGgIgx1oOHYU8cwGOoAUBGGOtBw7Klj\n",
"GAx1AKgIQx1oOPbUMQyGOgBUhKEONBx76hgGQx0AKsJdGoEJW+vWtt/8/utD/WfSl22y3RFd4xJ9\n",
"K9+m/3nXG0MdmLDoW9t+eXZ6Lupc47DR/rzrje0XAKgIQx0AKsJQB4CKMNQBoCIMdQCoCEMdACrC\n",
"UAeAijDUAaAiIw91M2uZ2QtmdtzMZiOjAACjGeknSs3MJN0nqSXJJD0t6WhgFwBgBKM+U98l6VV3\n",
"v+Du5yWdNLMbA7sAACMY9d4vV0uaM7NH1HmmPtc9djIqDAAwvFGH+luSpiR9Vp2h/q/dY8t88ve2\n",
"vzni+1hmx1VX/MqZ+YtRpwOA6pi7D/9GZpskHZd0mzpbOEfcfe/S12u328OfHACgVqtlo7zdSENd\n",
"kszsI5L+TpJLus/dnxnpRACAMCMPdQBA8/DDRwBQEYY6AFSEoQ4AFRnbUG/qbQTM7DEze9HMjprZ\n",
"7d1jtzWh1cz2mtn3zOyhnmMrruMk17dP57J1bUDno2b2nJkdM7PrV+tpYGej1tPM7u+2PNPwtVyp\n",
"s1FruaT3cjM7bWZ3dh+vOIuGanX38F/qXLv+n5KukHSlpOPjeD8jtv2bpN9sYqs6t134Y0kPrdY2\n",
"6ealnSutaxM6ezo+LOlfmrqeSzu7v3+siesp6VZJjyZYy1slPdrktex23C3pPyTdGbWm43qm3uTb\n",
"CJh++SuUxrS6e1vS2QHaJtq8Qqe0fF2l5qztgqR3VulpWueiJq7nLZJeWaWlCY29nYsat5ZmdqWk\n",
"j0o6vEbTUK2j/kTpWpp8G4EFSU+a2VuS/lrNbu3XtqnP8Uk2967rX7n762rO2n5K0tdW6WnKei52\n",
"Sg1cTzN7XtI1kvZK2tmnZeJruaRTauBadt0t6Z8lva/7OOTv57iG+sC3EVhv7n63JJnZTZL+QdLf\n",
"qKGt6r+Om/ocn5gl6/qwpI+rAX8PzGy/Os9yXjGz9/fpmfh69nZKzVxPd99nZh+Q9ISkv+zTMvG1\n",
"XNL5h01cSzPbJulD7v6gmf15tyHk3/u4hvpJdb5kUDfiRndvwjPfXhck/UzS62pe6+KPB6+4jta5\n",
"TUMTmlf6MebFdZUm/PfAzG6WNOPuX1itZ9LruUJnr8asZ9dP1Pkp8hX/3Ux6LXssdvZq0lrulbTF\n",
"zJ6UdIOkyyS9sFLT0Gs6xm8AfETSd7qhH5nENyH6dD0l6XlJ35a0o0mtkr4o6Zik/5b09e6xj67U\n",
"NsnmPp3L1rUBnae6Tc9J+lqD13Olzkatp6RvSGp3e3au1jLhtVyps1FruULz7ZLujPr7yW0CAKAi\n",
"/PARAFSEoQ4AFWGoA0BFGOoAUBGGOgBUhKEOABVhqANARRjqAFCR/wci6QM3xWJ/mwAAAABJRU5E\n",
"rkJggg==\n"
],
"text/plain": [
"<matplotlib.figure.Figure at 0x11692d3c8>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"flu_unique[flu_unique.AgeYears< 1].AgeDays.hist(bins=15)"
]
},
{
"cell_type": "code",
"execution_count": 52,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"172"
]
},
"execution_count": 52,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"flu_unique.flu_staph_org.sum()"
]
},
{
"cell_type": "code",
"execution_count": 53,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"84"
]
},
"execution_count": 53,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"flu_unique.s_aureus_icd9.sum()"
]
},
{
"cell_type": "code",
"execution_count": 55,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"218"
]
},
"execution_count": 55,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"(flu_unique.flu_staph_org | flu_unique.s_aureus_icd9).sum()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Logistic regression model"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Fate of each patient"
]
},
{
"cell_type": "code",
"execution_count": 35,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"died = (flu_complete.fate=='Died').astype(int).values"
]
},
{
"cell_type": "code",
"execution_count": 36,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"N_obs = len(died)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Functions for standarizing and centering"
]
},
{
"cell_type": "code",
"execution_count": 37,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"center = lambda x: (x - x[np.isnan(x) ^ True].mean())\n",
"standardize = lambda x: center(x) / x[np.isnan(x) ^ True].std()"
]
},
{
"cell_type": "code",
"execution_count": 38,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"['AgeYears', 'pH', 'OI', 'VA']"
]
},
"execution_count": 38,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"covariates"
]
},
{
"cell_type": "code",
"execution_count": 39,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"AgeYears, pH, OI, va = flu_complete[covariates].values.T\n",
"\n",
"AgeYears_center = center(AgeYears).round(1)\n",
"\n",
"# Center pH at 7.4\n",
"pH_center = pH - 7.4\n",
"\n",
"OI_std = standardize(OI)"
]
},
{
"cell_type": "code",
"execution_count": 63,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"no_coinfection = flu_complete.flu_only.values\n",
"with_staph = flu_complete.flu_staph.values"
]
},
{
"cell_type": "code",
"execution_count": 101,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"172"
]
},
"execution_count": 101,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"flu_unique.flu_staph.sum()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Model"
]
},
{
"cell_type": "code",
"execution_count": 40,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"# Values for missing elements, to be replaced in model\n",
"OI_std[np.isnan(OI_std)] = -999\n",
"pH_center[np.isnan(pH_center)] = -999\n",
"va[np.isnan(va)] = -999"
]
},
{
"cell_type": "code",
"execution_count": 90,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"from pymc3 import Model, find_MAP, NUTS, Metropolis, sample, forestplot, traceplot, Slice\n",
"from pymc3 import Normal, HalfCauchy, Beta, Bernoulli, MvNormal\n",
"from pymc3.distributions.timeseries import GaussianRandomWalk\n",
"\n",
"from numpy.ma import masked_values"
]
},
{
"cell_type": "code",
"execution_count": 56,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"def interpolate(x0, y0, x):\n",
"\n",
" idx = np.searchsorted(x0, x)\n",
" dl = np.array(x - x0[idx - 1])\n",
" dr = np.array(x0[idx] - x)\n",
" d = dl + dr\n",
" wl = dr / d\n",
"\n",
" return wl*y0[idx-1] + (1-wl)*y0[idx]"
]
},
{
"cell_type": "code",
"execution_count": 60,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"import theano.tensor as T\n",
"\n",
"def invlogit(x):\n",
" return 1. / (1 + T.exp(-x))"
]
},
{
"cell_type": "code",
"execution_count": 86,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"age_d = np.array([[xi - yi for yi in AgeYears_center] for xi in AgeYears_center])"
]
},
{
"cell_type": "code",
"execution_count": 97,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"with Model() as model:\n",
"\n",
" # Impute missing values\n",
" μ_pH = Normal('μ_pH', 0, 0.001, testval=0)\n",
" σ_pH = HalfCauchy('σ_pH', 5, testval=1)\n",
" pH_imputed = Normal('pH_imputed', μ_pH, σ_pH**-2, observed=masked_values(pH_center, value=-999))\n",
"\n",
" p_va = Beta('p_va', 1, 1)\n",
" va_imputed = Bernoulli('va_imputed', p_va, observed=masked_values(va, value=-999))\n",
"\n",
" # Non-linear age effect\n",
" '''\n",
" nknots = 10\n",
" knots = np.linspace(AgeYears_center.min(), AgeYears_center.max(), nknots)\n",
"\n",
" age_sd = HalfCauchy('age_sd', 5, testval=1)\n",
" age_knots = GaussianRandomWalk('age_knots', sd=age_sd, shape=nknots)\n",
" age_effects = interpolate(knots, age_knots, AgeYears_center)\n",
" '''\n",
" \n",
" # Parameters of kernel function\n",
" η2 = pm.HalfCauchy('η2', 5, testval=1)\n",
" ρ2 = pm.HalfCauchy('ρ2', 5, testval=1)\n",
" σ2 = pm.HalfCauchy('σ2', 5, testval=1)\n",
" \n",
" # Construct precision matrix\n",
" M = η2 * T.exp(-ρ2 * age_d**2)\n",
" S = M + T.identity_like(M)*σ2\n",
" τ = T.nlinalg.matrix_inverse(S)\n",
" \n",
" # GP of age effect\n",
" age_effects = MvNormal('age_effects', np.zeros(N_obs), τ, shape=N_obs)\n",
" \n",
" μ = Normal('μ', 0, 0.01, testval=0)\n",
" \n",
" # Linear covariates\n",
" β = Normal('β', 0, 0.01, shape=4)\n",
" \n",
" # Probability of mortality effect\n",
" π = invlogit(μ + β[0]*no_coinfection + β[1]*with_staph + β[2]*pH_imputed + β[3]*va_imputed + age_effects)\n",
" \n",
" deaths = Bernoulli('deaths', π, observed=died)"
]
},
{
"cell_type": "code",
"execution_count": 98,
"metadata": {
"collapsed": false
},
"outputs": [
{
"ename": "LinAlgError",
"evalue": "Singular matrix\nApply node that caused the error: MatrixInverse(Elemwise{Composite{((i0 * exp((i1 * i2 * i3))) + (i4 * i5))}}.0)\nInputs types: [TensorType(float64, matrix)]\nInputs shapes: [(1345, 1345)]\nInputs strides: [(10760, 8)]\nInputs values: ['not shown']\n\nBacktrace when the node is created:\n File \"<ipython-input-97-90b626f0aaef>\", line 29, in <module>\n τ = T.nlinalg.matrix_inverse(S)\n\nHINT: Use the Theano flag 'exception_verbosity=high' for a debugprint and storage map footprint of this apply node.",
"output_type": "error",
"traceback": [
"\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
"\u001b[0;31mLinAlgError\u001b[0m Traceback (most recent call last)",
"\u001b[0;32m/usr/local/lib/python3.4/site-packages/theano/compile/function_module.py\u001b[0m in \u001b[0;36m__call__\u001b[0;34m(self, *args, **kwargs)\u001b[0m\n\u001b[1;32m 594\u001b[0m \u001b[0;32mtry\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 595\u001b[0;31m \u001b[0moutputs\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfn\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 596\u001b[0m \u001b[0;32mexcept\u001b[0m \u001b[0mException\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/usr/local/lib/python3.4/site-packages/theano/gof/op.py\u001b[0m in \u001b[0;36mrval\u001b[0;34m(p, i, o, n)\u001b[0m\n\u001b[1;32m 767\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0mrval\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mp\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mp\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mi\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mnode_input_storage\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mo\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mnode_output_storage\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mn\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mnode\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 768\u001b[0;31m \u001b[0mr\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mp\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mn\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m[\u001b[0m\u001b[0mx\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mx\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mi\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mo\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 769\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mo\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mnode\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0moutputs\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/usr/local/lib/python3.4/site-packages/theano/tensor/nlinalg.py\u001b[0m in \u001b[0;36mperform\u001b[0;34m(self, node, xxx_todo_changeme2, xxx_todo_changeme3)\u001b[0m\n\u001b[1;32m 76\u001b[0m \u001b[0;34m(\u001b[0m\u001b[0mz\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m)\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mxxx_todo_changeme3\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 77\u001b[0;31m \u001b[0mz\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mnumpy\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mlinalg\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0minv\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mx\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mastype\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mx\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdtype\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 78\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/usr/local/lib/python3.4/site-packages/numpy/linalg/linalg.py\u001b[0m in \u001b[0;36minv\u001b[0;34m(a)\u001b[0m\n\u001b[1;32m 519\u001b[0m \u001b[0mextobj\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mget_linalg_error_extobj\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0m_raise_linalgerror_singular\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 520\u001b[0;31m \u001b[0mainv\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0m_umath_linalg\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0minv\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0ma\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0msignature\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0msignature\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mextobj\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mextobj\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 521\u001b[0m \u001b[0;32mreturn\u001b[0m \u001b[0mwrap\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mainv\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mastype\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mresult_t\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/usr/local/lib/python3.4/site-packages/numpy/linalg/linalg.py\u001b[0m in \u001b[0;36m_raise_linalgerror_singular\u001b[0;34m(err, flag)\u001b[0m\n\u001b[1;32m 89\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0m_raise_linalgerror_singular\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0merr\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mflag\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 90\u001b[0;31m \u001b[0;32mraise\u001b[0m \u001b[0mLinAlgError\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m\"Singular matrix\"\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 91\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;31mLinAlgError\u001b[0m: Singular matrix",
"\nDuring handling of the above exception, another exception occurred:\n",
"\u001b[0;31mLinAlgError\u001b[0m Traceback (most recent call last)",
"\u001b[0;32m<ipython-input-98-ffaed00ebadd>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[1;32m 1\u001b[0m \u001b[0;32mwith\u001b[0m \u001b[0mmodel\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 2\u001b[0;31m \u001b[0mstart\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mfind_MAP\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 3\u001b[0m \u001b[0;31m#step1 = NUTS([β, μ_pH, σ_pH, p_va, age_sd, age_knots, μ], scaling=start)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 4\u001b[0m \u001b[0mstep1\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mNUTS\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mβ\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mμ_pH\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mσ_pH\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mp_va\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mη2\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mρ2\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mσ2\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mμ\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mscaling\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mstart\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 5\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/Users/fonnescj/GitHub/pymc3/pymc3/tuning/starting.py\u001b[0m in \u001b[0;36mfind_MAP\u001b[0;34m(start, vars, fmin, return_raw, disp, model, *args, **kwargs)\u001b[0m\n\u001b[1;32m 79\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0;34m'fprime'\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mgetargspec\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfmin\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0margs\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 80\u001b[0m r = fmin(logp_o, bij.map(\n\u001b[0;32m---> 81\u001b[0;31m start), fprime=grad_logp_o, disp=disp, *args, **kwargs)\n\u001b[0m\u001b[1;32m 82\u001b[0m \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 83\u001b[0m \u001b[0mr\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mfmin\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mlogp_o\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mbij\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmap\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mstart\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdisp\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mdisp\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m*\u001b[0m\u001b[0margs\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/usr/local/lib/python3.4/site-packages/scipy/optimize/optimize.py\u001b[0m in \u001b[0;36mfmin_bfgs\u001b[0;34m(f, x0, fprime, args, gtol, norm, epsilon, maxiter, full_output, disp, retall, callback)\u001b[0m\n\u001b[1;32m 782\u001b[0m 'return_all': retall}\n\u001b[1;32m 783\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 784\u001b[0;31m \u001b[0mres\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0m_minimize_bfgs\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mf\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mx0\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0margs\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mfprime\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mcallback\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mcallback\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mopts\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 785\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 786\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mfull_output\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/usr/local/lib/python3.4/site-packages/scipy/optimize/optimize.py\u001b[0m in \u001b[0;36m_minimize_bfgs\u001b[0;34m(fun, x0, args, jac, callback, gtol, norm, eps, maxiter, disp, return_all, **unknown_options)\u001b[0m\n\u001b[1;32m 855\u001b[0m \u001b[0malpha_k\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mfc\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mgc\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mold_fval\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mold_old_fval\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mgfkp1\u001b[0m \u001b[0;34m=\u001b[0m\u001b[0;31m \u001b[0m\u001b[0;31m\\\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 856\u001b[0m _line_search_wolfe12(f, myfprime, xk, pk, gfk,\n\u001b[0;32m--> 857\u001b[0;31m old_fval, old_old_fval)\n\u001b[0m\u001b[1;32m 858\u001b[0m \u001b[0;32mexcept\u001b[0m \u001b[0m_LineSearchError\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 859\u001b[0m \u001b[0;31m# Line search failed to find a better solution.\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/usr/local/lib/python3.4/site-packages/scipy/optimize/optimize.py\u001b[0m in \u001b[0;36m_line_search_wolfe12\u001b[0;34m(f, fprime, xk, pk, gfk, old_fval, old_old_fval, **kwargs)\u001b[0m\n\u001b[1;32m 690\u001b[0m ret = line_search_wolfe1(f, fprime, xk, pk, gfk,\n\u001b[1;32m 691\u001b[0m \u001b[0mold_fval\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mold_old_fval\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 692\u001b[0;31m **kwargs)\n\u001b[0m\u001b[1;32m 693\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 694\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mret\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0;32mNone\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/usr/local/lib/python3.4/site-packages/scipy/optimize/linesearch.py\u001b[0m in \u001b[0;36mline_search_wolfe1\u001b[0;34m(f, fprime, xk, pk, gfk, old_fval, old_old_fval, args, c1, c2, amax, amin, xtol)\u001b[0m\n\u001b[1;32m 94\u001b[0m stp, fval, old_fval = scalar_search_wolfe1(\n\u001b[1;32m 95\u001b[0m \u001b[0mphi\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mderphi\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mold_fval\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mold_old_fval\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mderphi0\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 96\u001b[0;31m c1=c1, c2=c2, amax=amax, amin=amin, xtol=xtol)\n\u001b[0m\u001b[1;32m 97\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 98\u001b[0m \u001b[0;32mreturn\u001b[0m \u001b[0mstp\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mfc\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mgc\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mfval\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mold_fval\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mgval\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/usr/local/lib/python3.4/site-packages/scipy/optimize/linesearch.py\u001b[0m in \u001b[0;36mscalar_search_wolfe1\u001b[0;34m(phi, derphi, phi0, old_phi0, derphi0, c1, c2, amax, amin, xtol)\u001b[0m\n\u001b[1;32m 165\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mtask\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;36m2\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m==\u001b[0m \u001b[0;34mb'FG'\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 166\u001b[0m \u001b[0malpha1\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mstp\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 167\u001b[0;31m \u001b[0mphi1\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mphi\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mstp\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 168\u001b[0m \u001b[0mderphi1\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mderphi\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mstp\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 169\u001b[0m \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/usr/local/lib/python3.4/site-packages/scipy/optimize/linesearch.py\u001b[0m in \u001b[0;36mphi\u001b[0;34m(s)\u001b[0m\n\u001b[1;32m 80\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0mphi\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0ms\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 81\u001b[0m \u001b[0mfc\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m+=\u001b[0m \u001b[0;36m1\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 82\u001b[0;31m \u001b[0;32mreturn\u001b[0m \u001b[0mf\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mxk\u001b[0m \u001b[0;34m+\u001b[0m \u001b[0ms\u001b[0m\u001b[0;34m*\u001b[0m\u001b[0mpk\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m*\u001b[0m\u001b[0margs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 83\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 84\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0mderphi\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0ms\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/usr/local/lib/python3.4/site-packages/scipy/optimize/optimize.py\u001b[0m in \u001b[0;36mfunction_wrapper\u001b[0;34m(*wrapper_args)\u001b[0m\n\u001b[1;32m 280\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0mfunction_wrapper\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m*\u001b[0m\u001b[0mwrapper_args\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 281\u001b[0m \u001b[0mncalls\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m+=\u001b[0m \u001b[0;36m1\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 282\u001b[0;31m \u001b[0;32mreturn\u001b[0m \u001b[0mfunction\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m*\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mwrapper_args\u001b[0m \u001b[0;34m+\u001b[0m \u001b[0margs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 283\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 284\u001b[0m \u001b[0;32mreturn\u001b[0m \u001b[0mncalls\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mfunction_wrapper\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/Users/fonnescj/GitHub/pymc3/pymc3/tuning/starting.py\u001b[0m in \u001b[0;36mlogp_o\u001b[0;34m(point)\u001b[0m\n\u001b[1;32m 71\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 72\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0mlogp_o\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mpoint\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 73\u001b[0;31m \u001b[0;32mreturn\u001b[0m \u001b[0mnan_to_high\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m-\u001b[0m\u001b[0mlogp\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mpoint\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 74\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 75\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0mgrad_logp_o\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mpoint\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/Users/fonnescj/GitHub/pymc3/pymc3/blocking.py\u001b[0m in \u001b[0;36m__call__\u001b[0;34m(self, x)\u001b[0m\n\u001b[1;32m 117\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 118\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0m__call__\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mx\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 119\u001b[0;31m \u001b[0;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfa\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfb\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mx\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
"\u001b[0;32m/Users/fonnescj/GitHub/pymc3/pymc3/model.py\u001b[0m in \u001b[0;36m__call__\u001b[0;34m(self, state)\u001b[0m\n\u001b[1;32m 312\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 313\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0m__call__\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mstate\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 314\u001b[0;31m \u001b[0;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mf\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m**\u001b[0m\u001b[0mstate\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 315\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 316\u001b[0m \u001b[0;32mclass\u001b[0m \u001b[0mLoosePointFunc\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mobject\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/usr/local/lib/python3.4/site-packages/theano/compile/function_module.py\u001b[0m in \u001b[0;36m__call__\u001b[0;34m(self, *args, **kwargs)\u001b[0m\n\u001b[1;32m 604\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfn\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mnodes\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfn\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mposition_of_error\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 605\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfn\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mthunks\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfn\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mposition_of_error\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 606\u001b[0;31m storage_map=self.fn.storage_map)\n\u001b[0m\u001b[1;32m 607\u001b[0m \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 608\u001b[0m \u001b[0;31m# For the c linker We don't have access from\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/usr/local/lib/python3.4/site-packages/theano/gof/link.py\u001b[0m in \u001b[0;36mraise_with_op\u001b[0;34m(node, thunk, exc_info, storage_map)\u001b[0m\n\u001b[1;32m 204\u001b[0m exc_value = exc_type(str(exc_value) + detailed_err_msg +\n\u001b[1;32m 205\u001b[0m '\\n' + '\\n'.join(hints))\n\u001b[0;32m--> 206\u001b[0;31m \u001b[0;32mraise\u001b[0m \u001b[0mexc_type\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mexc_value\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mwith_traceback\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mexc_trace\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 207\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 208\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/usr/local/lib/python3.4/site-packages/theano/compile/function_module.py\u001b[0m in \u001b[0;36m__call__\u001b[0;34m(self, *args, **kwargs)\u001b[0m\n\u001b[1;32m 593\u001b[0m \u001b[0mt0_fn\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mtime\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mtime\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 594\u001b[0m \u001b[0;32mtry\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 595\u001b[0;31m \u001b[0moutputs\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfn\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 596\u001b[0m \u001b[0;32mexcept\u001b[0m \u001b[0mException\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 597\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mhasattr\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfn\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m'position_of_error'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/usr/local/lib/python3.4/site-packages/theano/gof/op.py\u001b[0m in \u001b[0;36mrval\u001b[0;34m(p, i, o, n)\u001b[0m\n\u001b[1;32m 766\u001b[0m \u001b[0;31m# default arguments are stored in the closure of `rval`\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 767\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0mrval\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mp\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mp\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mi\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mnode_input_storage\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mo\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mnode_output_storage\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mn\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mnode\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 768\u001b[0;31m \u001b[0mr\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mp\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mn\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m[\u001b[0m\u001b[0mx\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mx\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mi\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mo\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 769\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mo\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mnode\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0moutputs\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 770\u001b[0m \u001b[0mcompute_map\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mo\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;32mTrue\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/usr/local/lib/python3.4/site-packages/theano/tensor/nlinalg.py\u001b[0m in \u001b[0;36mperform\u001b[0;34m(self, node, xxx_todo_changeme2, xxx_todo_changeme3)\u001b[0m\n\u001b[1;32m 75\u001b[0m \u001b[0;34m(\u001b[0m\u001b[0mx\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mxxx_todo_changeme2\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 76\u001b[0m \u001b[0;34m(\u001b[0m\u001b[0mz\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m)\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mxxx_todo_changeme3\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 77\u001b[0;31m \u001b[0mz\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mnumpy\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mlinalg\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0minv\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mx\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mastype\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mx\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdtype\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 78\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 79\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0mgrad\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0minputs\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mg_outputs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/usr/local/lib/python3.4/site-packages/numpy/linalg/linalg.py\u001b[0m in \u001b[0;36minv\u001b[0;34m(a)\u001b[0m\n\u001b[1;32m 518\u001b[0m \u001b[0msignature\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;34m'D->D'\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0misComplexType\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mt\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;32melse\u001b[0m \u001b[0;34m'd->d'\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 519\u001b[0m \u001b[0mextobj\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mget_linalg_error_extobj\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0m_raise_linalgerror_singular\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 520\u001b[0;31m \u001b[0mainv\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0m_umath_linalg\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0minv\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0ma\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0msignature\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0msignature\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mextobj\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mextobj\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 521\u001b[0m \u001b[0;32mreturn\u001b[0m \u001b[0mwrap\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mainv\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mastype\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mresult_t\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 522\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/usr/local/lib/python3.4/site-packages/numpy/linalg/linalg.py\u001b[0m in \u001b[0;36m_raise_linalgerror_singular\u001b[0;34m(err, flag)\u001b[0m\n\u001b[1;32m 88\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 89\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0m_raise_linalgerror_singular\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0merr\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mflag\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 90\u001b[0;31m \u001b[0;32mraise\u001b[0m \u001b[0mLinAlgError\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m\"Singular matrix\"\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 91\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 92\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0m_raise_linalgerror_nonposdef\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0merr\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mflag\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;31mLinAlgError\u001b[0m: Singular matrix\nApply node that caused the error: MatrixInverse(Elemwise{Composite{((i0 * exp((i1 * i2 * i3))) + (i4 * i5))}}.0)\nInputs types: [TensorType(float64, matrix)]\nInputs shapes: [(1345, 1345)]\nInputs strides: [(10760, 8)]\nInputs values: ['not shown']\n\nBacktrace when the node is created:\n File \"<ipython-input-97-90b626f0aaef>\", line 29, in <module>\n τ = T.nlinalg.matrix_inverse(S)\n\nHINT: Use the Theano flag 'exception_verbosity=high' for a debugprint and storage map footprint of this apply node."
]
}
],
"source": [
"with model:\n",
" start = find_MAP()\n",
" #step1 = NUTS([β, μ_pH, σ_pH, p_va, age_sd, age_knots, μ], scaling=start)\n",
" step1 = NUTS([β, μ_pH, σ_pH, p_va, η2, ρ2, σ2, μ], scaling=start)\n",
"\n",
" step2 = Slice([va_imputed.missing_values, \n",
" pH_imputed.missing_values])\n",
" \n",
" trace = sample(5000, [step1, step2], start)"
]
},
{
"cell_type": "code",
"execution_count": 76,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.gridspec.GridSpec at 0x11c9389b0>"
]
},
"execution_count": 76,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": [
"iVBORw0KGgoAAAANSUhEUgAAAaUAAAEgCAYAAAAUr2zFAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\n",
"AAALEgAACxIB0t1+/AAAFjFJREFUeJzt3XuUrXV93/H3Bw4S6xVBvHBxgUpWophU5aCJxsupkmTR\n",
"CDaiXaahxYq3pJakJk01CHVpWkUUtGIqxUPIQi6N1WIilh6PcpPrUsLyUhWOAeWmBAXxEgPf/rGf\n",
"g9tn9p6ZPWdmnt+ceb/Wehazn+fZz+87v9lnf+b37N/8SFUhSVILdhm6AEmStjOUJEnNMJQkSc0w\n",
"lCRJzTCUJEnNMJQkSc0wlLRmJNklyZ8l2Zbk1iQfTvLg3jlHJ/l+kpuS3Nz9d8/uWJKc1e2/LMlj\n",
"u/0nJ/mDGWt5dJLNSW7r2vhikv+0fN/t1HbPSnL8lGMfTPJbE/ZvS/LCea75oK5Pbklyf5IVe19I\n",
"ctFYOzcled8Sr/PWJH+x3PUttyQf7/r2J/P9DPRTG4YuQJrBG4FfBQ4GfgicDrwXeE3vvPOq6pgJ\n",
"zz8M2Keq9ktyAvDGJGcCm4Cnz1jLJ4DrgAOr6gdJHgbsP+M1llVVvXaJz/sHYL8kTwBuXN6q5rT1\n",
"orF2nlA7+R9KVtVLAJKsaL/uTBwpaS15PnB+VX2/qu4DjgNemWSv3nmZ8vx/BHbtRgIP6h6fAryp\n",
"u96iJPlNYF/gtVX1A4Cquqeqvtg778PdyO593ajgm0meMnY8Sf5jkq91o5nTkuzeu8bLu+M3dSOD\n",
"h06o59rut/EfJZkUxgAHJ/lMktuTXJLkyYv9fsfaOSbJl5P8XZJzkzxq1mv0Lzmhja1J3tT13TeS\n",
"3JDkeb1zrmX0sz9ybDT8+N45hyf5QnfswiT7944/oRutPTfJ55N8K8npY8f/PMmpved8JMnbxh7v\n",
"keQ93Sj51iTXJXnOYr/Xsesck+T6rta/S/KmaeeuC1Xl5rYmNuAtwCXAnt3jjcBdwAvGzvmX3b5t\n",
"wGeB5/Wu8d+AG4BPAscCH1tCHccD/2sR530YuB14dff4Yb3jxwHXA49l9KZ1BvD+seO/wGhE+Lzu\n",
"8bOBO4Hjp7S3FThmwv5t3ff7yO7xu4FrJ5z3BOA+YJcJx44EvgUc1D0+EfjEEn+O87WzFfgS8E+7\n",
"x/8ZuHjCeW8F/mLK9Z8B3A08q3t8DHDdhBruBy4EHtP/+QDPBO4ANnSPHw58H9h/7JzdgRcBGeuT\n",
"L06paRvwwgn7n9T9jH9hbN9DV+PfU6ubIyWtJe8AtgAXJ/kq8Frgx8De20+oqo8weuM8gNGbxMeS\n",
"HDB2/A1V9UTgCEa3A49bQh37APcs8tzzq+pDXdv95xwLvK2qbqvRu9GbgVeNHX8ZcEFVfbZ7/ueA\n",
"v1lCvQDvqqrvdl8fD/xSdxttsY4FTqmqr3aP3wY8P8njlljPfE6qqs93X1/M7LdF/y3wl1V1BUBV\n",
"nQHsnuRZvfMKeGVV3d6d98DPp6quYRTCh3e7Xg5cWlU3jZ3z46q6qPvZAfxP4KAZa/0e8CPgxdtH\n",
"e1X1/RmvsVMxlLRmVNX9VXVCVT2lqg5i9ObzcEajkfHzvt3999OM3tQOn3Mx+A/AR4HXJvlckvOS\n",
"PHKRpdzNWBAu4O/nObYf8O4kN3afOVwO/GDsjf6xjH7DXlZVdS+j0eSjZ3jafsBxY7V+jdFv+Cv9\n",
"OdpPmP19aj/gZdtrTbINeFS3v++uea5zOvCvu69/F/jQ+MHu9utxSS5OcgmjUfgus0wU6V6rGxmN\n",
"3LZ2t2I3Lfb5OyNDSWvZbzH6LfOKec7ZrTvnAUn2BV4JXAkcUlXPZnTL6N8vst0rgUOS/JOZK/5Z\n",
"NwBHV9WB3XZAVe1ZVbd2x78JHNB7zq5LbOuBSU0ZzUbcA7hp+ukTa31zr9ZHV9WVS6xnJd0InNar\n",
"de+qOn/G6/wlo9HgsxndZvt47/jvAf8C+O2qei6j4JpZVX2tqv6gqn6e0ej+gv5ni+uJoaQ1I8lD\n",
"kjyi+/pA4J2MPl/50dg5+2//TTXJi4CnAR/rXepdjP7x9yc3PJjF+ThwC/DnSR7etbVnkkkjsvm8\n",
"B3hvkp8fq/8RY8fPBX4zyaHdsaMYBfFS/Nckj0mygVG/XVhVd0w4b9oH8u8F3ppk41itix1ZTmpj\n",
"6gf/i3QXo8kbuyXZtXcb8QPAa5Ic9kCDk2udt4aq+h7wv4GPAGdW1T/2TtkPuA34djfZ5l3d/t1m\n",
"+UbGby8zmoDzY+a+NtcNQ0lryVOB65J8E/gUcGpVvb93zh8C30zyDeCPgBdvv50H0M2OemxVndtd\n",
"47YkNwO/DpzKIlTVTxhNI78f+Fr3/EsYfcD+M6cucJ3NjALinG7W1Y3A28eO3wD8PqPPxW4Cnsf8\n",
"nylNa6+AsxmF8+3AgYw+I5p27jeSfKFX69buOad0s8RuBDbP9/1NkuQi4LKxdvp/p7TYKeJnAz8A\n",
"bmY0WeSIsVq/Avxz4I+7ft3GaPTxoCW09SFG4fM/Jhx7N6NboLcAfw38d0Zhss+Ecwv4SNd3D8x8\n",
"7H6pOX/7LELgDcCvTwjAdWP7rBFJkgbnSEmS1AxDSZLUDENJktQMQ0mS1AwXZN1JbNmyxRkrktaM\n",
"TZs2TZyS7+w7SepJUlW1o39LpSXw9p0kqRmGkiSpGYaSJKkZhpIkzXXi0AWsV050kCQ1w5GSJKkZ\n",
"hpIkqRmGkiSpGYaSJKkZhpIk9SQ5Yega1itn30lSj8sMDceRkiSpGYaSJKkZhpIkqRmGkiSpGYaS\n",
"JM3l2ncDcfadJKkZjpQkSc0wlCRJzTCUJEnNMJQkSc0wlCSpx7XvhuPsO0nqce274ThSkiQ1Y12E\n",
"UpKHJnnzDOcfluTzST62g+2+ejnqkaT1wtt3EyT5ALC5qq7awetcXVWHLFNZklaJt++Gs2ojpSTX\n",
"JTkpyRVJ3jG2/+gkl3fbxJFF7zqHJNmSZGuSzQtdp9u/NcnVi6znr4AjgPclOXts/2HdtS9N8tJF\n",
"1HMmcFCSTyd5yyLqmVb/eJ1vX6h/JM2Vkad1m2HTsqpalQ24Cdgb2BW4vtu3J3B5t29X4DPA4xe4\n",
"zrXA3r19C14HuGqhesaOnQH84tjjAF8AHtKdfzHwoGn1TGtz2rH56p+vTjc3t4W30b/fjefB5rtH\n",
"26Hn0t0lmuc5Jwxd93rdNuxoqM3gtqq6AyDJD7t9BwLXVNV93f4rgCcDt0y6QJK9xq8zZqbrzFPP\n",
"A031Hu8F7Atc0B17BLBPknum1DPtOtPMV/98dUrqSeh9JtH/iOLoo4CjJo2Xqkb/ZqvqhBUpTgta\n",
"zYkOmfD1jcAzk2xIsgH4FeCr0y5QVd8BHpdk396hxVyn/xKcVM+0dr8NfAl4SVW9oKp+qaq2zVPP\n",
"drvNc6tgfP989S+6TmlnlFCzbGulLU22miOl6n9dVXcm+SCj22EAZ1TVrQtc51jgrO7N/o6qOmqR\n",
"1+m/gObUM89jgD8BLkhSwLeq6nem1TP2nIuATybZVlWvm9bGAvXPV6e009s+elmq0b/NjefC6w8b\n",
"7TntQrjyFdXdp1NbnH0naafX/dJ4cPfwegOpXYaSJKkZ6+KPZyVpFq59NxxHSpLU4x/PDseRkiSp\n",
"GYaSJKkZhpIkqRmGkiSpGYaSJM114tAFrFfOvpMkNcORkiSpGYaSJKkZhpIkqRmGkiSpGYaSJPW4\n",
"9t1wnH0nST2ufTccR0qSpGYYSpKkZhhKkqRmGEqSpGYYSpI0l2vfDcTZd5KkZjhSkiQ1w1CSJDXD\n",
"UJIkNcNQkiQ1w1CSpB7XvhuOs+8kqce174bjSEmS1AxDSZLUDENJktQMQ0mS1AxDSZLmcu27gTj7\n",
"TpLUDEdKkqRmGEqSpGYYSpKkZhhKkqRmGEqS1OPad8Nx9p0k9bj23XAcKUmSmmEoSZKaYShJkpph\n",
"KEmSmmEoSdJcrn03EGffSZKa4UhJktQMQ0mS1AxDSZLUDENJktQMQ0mSelz7bjjOvpOkHte+G44j\n",
"JUlSMwwlSVIzdiiUkrx6RwtIcvWOXmMRbTw0yZtnOP+wJJ9P8rEdbHdi/8xajyStFzv0mVKSq6vq\n",
"kB0qILmqqjbuyDWWW5IPAJur6qodvM4O94+k1ednSsNZcKSU5PVJrklyyfhv/knOBA5K8ukkbxnb\n",
"/2+SfDzJdUl+f2z/l5OcnuTKJH861sTuSU5KckWSty+inkOSbEmyNcnmsf1HJ7m8217d27+1PyLr\n",
"6tve7jvG9v8VcATwviRnj+0/rLv2pUleuoh6pvXPtHqm1T9e54L9I2lZrMjadxl5WrcZepNU1bwb\n",
"8GlgjynHrpqwb0P3392B68b2bwMewigILwce0+2/Gdgb2BW4fhH1XAvs3du3Z3fNXbvtM8Dj56sV\n",
"uGlau8AZwC+OPQ7wha7+XYGLgQdNq2e+/pl0bL7656vTzc1t7Wyj95GN58Hmu0fboefS3a1y++m2\n",
"YRG59a+A1yR5FPDxqrps7NikpP+1JIcD9wIPHtv/naq6FyDJtcD+wO3ArVV1R7f/h/MVkmQv4Lbt\n",
"5485ELimqu7rzrsCeDJwyzyXu22edvvf117AvsAF3bFHAPskuWdKPdOuM8189c9Xp6RGJCzwWUj/\n",
"8NFHAUfNN16qWvR7yE5jwVCqqm8B/yXJ7ox+m3/G2OHd0t18Hdt3alU9Ncn+wMvH9j8+yR7A3cDT\n",
"geO7/eOdPu8PoKq+k+RxSfatqm+OHboReGaS7d/PrwCn9J7ev/Ys7X47yZeAl1TVPT9zkcn1bDep\n",
"fya1OV/9i65T0o5ZOFhW13LWs1YCbsFQSvIe4JeBhwHv7x2+CPhkkm1V9bpu36VJLmN0u+vvx879\n",
"LnAy8BTgrKq6q9s/3umL+QEcC5zV3Y+9o6qOqqo7k3yQ0W01gDOq6tbe8/rXnq/dSXX8CXBBkgK+\n",
"VVW/M62esedM6p85bSxQ/6z9I2mJVvKNe/QesfFceP1hoz2nXQhXvmLKL63r1qqt6OBMNEnrXffL\n",
"68Hdw+sNpLlW849n7XxJa8JKrX1XI3/bbb4nTuDad5LU498pDcdlhiRJzTCUJEnNMJQkSc0wlCRJ\n",
"zTCUJGmuFVn7Tgtz9p0kqRmOlCRJzTCUJEnNMJQkSc0wlCRJzTCUJKlnpda+08KcfSdJPa59NxxH\n",
"SpKkZhhKkqRmGEqSpGYYSpKkZhhKkjSXa98NxNl3kqRmOFKSJDXDUJIkNcNQkiQ1w1CSJDXDUJKk\n",
"Hte+G46z7ySpx7XvhuNISZLUDENJktQMQ0mS1AxDSZLUDENJkuZy7buBOPtOktQMR0qSpGYYSpKk\n",
"ZhhKkqRmGEqSpGYYSpLU49p3w3H2nST1uPbdcBwpSZKaYShJkpphKEmSmmEoSZKaYShJ0lyufTcQ\n",
"Z99JkprhSEmS1AxDSZLUDENJktQMQ0mS1AxDSZJ6XPtuOM6+k6Qe174bjiMlSVIzDCVJUjMMJUlS\n",
"MwwlSVIzDKXGJbl60teSVpRr3w1kw9AFaEE15WvNKEmAg7uH15dTTzVFVZ0wdA3rlSOlRiT5cpLT\n",
"k1yZ5E/HD035WjMYBdLGc2HzpaPt0HO6kJLUEP9OqRFJtgFPBX4IXAocWVW3J7kbuIZRID2xqvYf\n",
"sMxVl6yd0WGVvzRIO8rbd+34TlXdC5DkWmB/4HbgK1X1wm7/Vcvd6Fp6029dy31pYGqtMJTa8fgk\n",
"ewB3A08Hju/2r+jtu/XyZvXT23evP2y057QL4cpX+LmS1BZDqR3fBU4GngKcVVV3dfud6LAMqqqS\n",
"vByucqKDFpTkBCc7DMPPlBqR5OqqOmToOiS59t2QnH3XDn87kLTuOVKSpB5HSsNxpCRJaoahJElq\n",
"hqEkSXO59t1A/ExJktQMR0qSpGYYSpKkZhhKkqRmGEqSpGYYSpLUk+SEoWtYr5x9J0k9rugwHEdK\n",
"kqRmGEqSpGYYSpKkZhhKkqRmGEqSNJdr3w3E2XeSpGY4UpIkNcNQkiQ1w1CSJDXDUJIkNcNQkqQe\n",
"174bjrPvJKnHte+G40hJktQMQ0mS1AxDSZLUDENJktQMQ0mS5nLtu4E4+06S1AxHSpKkZhhKkqRm\n",
"GEqSpGYYSpKkZhhKktTj2nfDcfadJPW49t1wHClJkpphKEmSmmEoSZKaYShJkpphKEnSXK59NxBn\n",
"30mSmuFISZLUDENJktQMQ0mS1AxDSZLUDENJknpc+244zr6TpB7XvhuOIyVJUjMMJUlSMwwlSVIz\n",
"DCVJUjMMJUmay7XvBmIoNSDJqUmeO/Z4U5KTuq/PTHLOcNVJ609VnTB0DeuVodSGjwJHjj0+Evho\n",
"kg3AwcATk+w2SGWS1ryMPK3bmp7qbii14WLgV8ceb6yqy4EXAJ8DLgNePERhkta2UQhtPBc2Xzra\n",
"Dj2n5WDaMHQBgqq6P8l1SQ4Gfg64ujv0UuATwH3AUcBfD1SipMYlTFkJob/76KOAoybFUhWDh5Wh\n",
"1I7tt/B2Z3TrLsBvAHsBAZ6VZJequn/AGiWtgukB02a7yxlm3r5rxxZGt+t+DdgKPAf4bFW9rKp+\n",
"G/g/wPOHK09aP4Ze+66KLNcG2QUOPR/OvHu0Pes8yC7L28byce27hiQ5C/iHqnpVkpOBT1XVp7pj\n",
"/ww4oqp+b9AipXVgZ1v7rrvzcnD38Ppq+I3fUJKknp0tlNYSb99JkpphKEmSmmEoSZKaYShJ0lyu\n",
"fTcQJzpIkprhSEmS1AxDSZLUDENJktQMQ0mS1AxDSZJ6hl77bj1z9p0k9bjM0HAcKUmSmmEoSZKa\n",
"YShJkpphKEmSmmEoSdJcrn03EGffSZKa4UhJktQMQ0mS1AxDSZLUDENJktQMQ0mSelz7bjjOvpOk\n",
"Hte+G44jJUlSMwwlSVIzDCVJUjMMJUlSM5zosJPYsmWLP0hJa8amTZsmTiQxlCRJzfD2nSSpGYaS\n",
"JKkZhpIkqRmGkiSpGYbSOpPkg0m2JvlMkgMWcf6Hk3wuyaeT/G5jtW1KckmSi5O8cCVr69p7TpKr\n",
"krxzkeevZt/NWtuq9d2sba1Gv81S02q/zpZQX5OvsyX3W1W5rcMNeAFw2iLOOwPYr7XagACXAT8H\n",
"PBi4eBXq2gQcAbxzkeevWt/NUttq9t1S2lrpfpulpoFeZzO12eLrbEf6zZHS+nUP8ONFnBdWf0S9\n",
"mNqeDPy/qvpRVf0Q+HqSJ61kUVW1BbhrhqesWt/NWNtq9t1S2lrpfpulplV/nS2hzRZfZ0vutw07\n",
"UqDaleRFwB8BxehFW8AfVtX13SmvAk5ZxKXuAc5OcidwXFXd0EhtewLfS3Jyd43vdfu+vgr1LdYQ\n",
"fbcYK9J3U2p72xLaWvZ+65nl+1+x19ky1Qcr319LseR+M5R2UlV1EXDRpGNJDmf0W8xXFnGdf9c9\n",
"55eBk4AjG6ntTuCRwOsYvehP6/btsPnqm/E6q9p3M1iRvptUW5KDZm1rJfqtZ5bvf8VeZ8tU32r0\n",
"11Isud+8fbfOJHkG8Pyqeu+MT/0R8JMVKOkBM9b2dUa3CGD0on9SVa3kb6/jZv3/7Kx4341ZTG2r\n",
"2Xc70tZK9dssNQ3xOltqmy29zpbcb46U1p/zgZuTbAX+tqreuP1AkpcB91bV34ztOwd4HKNbBG9o\n",
"pbaquj/JicD/ZXSb6MQVro0kfwz8BvCYJA+vqtdMq6/bt2p9N0ttq9l3C7U1RL/NV1MLr7NZ6uv2\n",
"Df46W85+c+07SVIzvH0nSWqGoSRJaoahJElqhqEkSWqGoSRJaoahJElqhqEkSWqGoSRJasb/B3xD\n",
"9ylRkix1AAAAAElFTkSuQmCC\n"
],
"text/plain": [
"<matplotlib.figure.Figure at 0x11c938a58>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"forestplot(trace, vars=['β'], ylabels=['no coinfection', 'staph coinfection', 'pH', 'VA'])"
]
},
{
"cell_type": "code",
"execution_count": 78,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"array([[<matplotlib.axes._subplots.AxesSubplot object at 0x1157bbc18>,\n",
" <matplotlib.axes._subplots.AxesSubplot object at 0x11cb93b70>],\n",
" [<matplotlib.axes._subplots.AxesSubplot object at 0x11cbde5c0>,\n",
" <matplotlib.axes._subplots.AxesSubplot object at 0x11cc1f438>]], dtype=object)"
]
},
"execution_count": 78,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": [
"iVBORw0KGgoAAAANSUhEUgAAA1gAAAEXCAYAAACwDuU4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\n",
"AAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xm4JHV59//3PbLNAMMmICLIMixiRBRLXFiGPiCiqLiw\n",
"+MT4xBgVjeKStPpoEkH9oViCGo1oiBpJFBAVFWUbplBQEVsMi7I5ArKKijIMMAMMc//+qKo5NX16\n",
"qe6u7uqu/ryuq69zqrqW+67Tp6u+Vd/F3B0REREREREZ3LyyAxAREREREakKFbBEREREREQKogKW\n",
"iIiIiIhIQVTAEhERERERKYgKWCIiIiIiIgVRAUtERERERKQgKmCJTAkzu8PMDiw7DhERkWY6R0mV\n",
"qIAlIiIiIiJSEBWwRHpgZhuY2UfM7Bozu9PMbjazIzPvb2dm5yfv/drMLjWzu81st8wye5rZkmSZ\n",
"X5rZwX3EsWeyn9uT18+a3jczO8nM7jGzG83sA4NlLiIi407nKJHxoAKWSA/c/VHgp8Bz3P0pwMeA\n",
"r5pZ+r/0KeCW5L0PAPsCzwKWAZjZxsAS4LvJMn8LnG1m2/cYyn8CP3f3Hd19R+BFTe+/CTgmiXNP\n",
"4D5gux73ISIiE0TnKJHxoAKWSI/c/QJ3fyyZ/CawKZCefJ4FnJf8fhGwSbKOJ/OOAFa4++eS+dcC\n",
"3wb+pscw7gCeZWbPMTNz9wea3j8WOMXd70r280Xgnh73ISIiE0bnKJHyqYAl0iMze72ZRWZ2OfGJ\n",
"x4H1k7evAI5J7hYeC/zW3e/NrL4DsKOZ3ZK8bgVeTu937t4AXAx8HLjTzD7c9P6TgFt73KaIiEw4\n",
"naNEyrde2QGITBIzexnwL8Ch7n5bMm9NZpE6cDVwN3A98IqmTdwC/NLdB+opyd1XAf8O/LuZbQss\n",
"NbNl7n5GssidwM5Nqz1hkH2KiMh40zlKZDzoCZZIb3Ygrit+l5ltYmanAY8ze3fwy8DH3f1J7l5z\n",
"9+ub1v8B8EQzq5vZegBmtpGZbdBLEGa2g5mlJ6P1iE9MKzOLnAkcb2bbmNn6ZvZxYJte9iEiIhNH\n",
"5yiRMaAClkhv/gu4l7h++Y+BpcR34tL67TcDJyfjedyRVLE4y8wWALj7I8AM8EzgN2Z2G3AlsFeP\n",
"cbwN+J2Z3Z7E8CV3Pyfz/leJq2fcBPwqiffuHvchIiKT5b/QOUqkdDbbrlFEBmVm3wa+6u7fTaY3\n",
"Aa4B6u7+7VKDExGRqaZzlMhoDPUJlpl9IRlj4Ydm1lzXtnnZGTO73MwuM7PaMOMSGQYz25T4zt+1\n",
"yfT6wNHAFsBVJYYmIgUws4XJOS1Kft5fdkwieekcJTI6I3mClQxSd7S7v7XN+0b8KHsGMOCiQRtY\n",
"ipTBzP4aeD9xt7gO/IS4vvuvcqy7PfCzZL05byfzT3X3TxcXsYj0w8z2Bt7h7m8qOxaRvHSOEhmN\n",
"UfUiuAJ4pMP7uwE3Jb3OYGbLzGyRuy8bSXQiBXH3rwFf63Pdu4gbKIvI+HsH8NmygxDphc5RIqMx\n",
"qgLWG4HPdHh/K2C5mZ1KfBdkeTKvZQFr6dKlajgmIjKhZmZmrOwYBmFmWwI7JIOwzqFzlIjI5Cri\n",
"HDX0ApaZHUH8dOrGDovdB2wOvJW4gHVaMq+tST9BZ5nZYnf/YdlxFEk5jb+q5QPVy6lq+UBlCh9v\n",
"Bk7vtECVzlEiItOiqHPUUAtYZrYvsNjd/6nLosuIqwlCXMBS9UARERk7ydg+LwMOKDsWEREZT8N+\n",
"gnUOcIeZXQpc6+7vBDCzo4CH3P18AHdfY2YnApcQN5I8cchxjZWq3aEG5TQJqpYPVC+nquVTEa8E\n",
"znP3NWUHIiIi42moBSx336XN/HNazFsCLBlmPCIiIoNw92+WHYOIiIy3oY6DJfmY2eKyYyiachp/\n",
"VcsHqpdT1fIRERGZBipgiYiIiIiIFGRU3bSPtSCMDDgBuK1Rr31l1PuvYjsL5TT+qpYPVC+nquUj\n",
"IiIyDfQEK/YB4A3APwdhpGMiIlJRZraVmX3JzC5Kps3M3l52XCIiUh1TX5gIwug44O+A/YAHgENG\n",
"HUMV21kop/FXtXygejlVLZ8xcTpwPrAhgLs7cEypEYmISKVMdQErCKNjgH8BDm3Ua/cAXwCOKzcq\n",
"EREZoi3d/VvEQ4KkpvpcKCIixZrak0pSFfALwEsb9dotyeyvA4uDMHryKGOpYjsL5TT+qpYPVC+n\n",
"quUzJtaY2XYkBSwzeyWwPM+KZra9mUVmdpmZnTLMIEVEZHJNbQEL2AP4c6Neuzqd0ajXVgBnA28s\n",
"LSoRERmmfwIuAPYxs6uAjwDvyLnuJ4EPuvuB7v6PwwpQREQm2zQXsAKg0WL+F4A3BWE0sh4Wq9jO\n",
"QjmNv6rlA9XLqWr5jAN3/yXx9/9BwOuAvd39t93WM7N5wCJ3v2LIIYqIyISb5gLWc2hRwGrUa9cA\n",
"dwOHjzwiEREZOnd/zN2vc/cb3H1NztW2BjYys3PNbGlStVBERGSOaR4HKwC+1ea9tLOL80YRSBXb\n",
"WSin8Ve1fKB6OVUtn3FgZu9pNd/dT+2y6n3A/cCric+dPzGzC919ZcEhiojIhJvKJ1hBGK0P7A38\n",
"ss0i5wAHBmG06eiiEhGREdi06XUw8KxuK7n7auAOYDt3fxRYNcwgRURkck1lAQv4K+C2pFOLORr1\n",
"2kPAz4DaKIKpYjsL5TT+qpYPVC+nquUzDtz9xOwLeBXx06k83g+cbmY/Br6pp1ciItLKtFYRbNfB\n",
"RdaFwIuB7w4/HBERKYO7P2Zm2+Rc9nbgJUMOSUREJpwKWO1dCHw/CCNr1GveZdmBVLGdhXIaf1XL\n",
"B6qXU9XyGQdmdh7rDjK8DXBXSeGIiEgFdS1gmdlm7p5rEMYJEgCnd1nmeuLjsxtw89AjEhGRUfhk\n",
"0/Rf3P3aUiIREZFKytMG6xoz+7KZ7dfrxs1sfzP7uZl9IseyXzGzK8wsMrPX97qvvIIwWgDsDlzT\n",
"abnkqdVFxNUEh6qK7SyU0/irWj5QvZyqls84cPcfNb1UuBIRkULlKWAtIu6u/INm9gszO87MNsm5\n",
"/Q2Bk3Iu68DR7l5z9zNyrtOPfYAbGvXaIzmWTdthiYiIiIiIdNW1gOXuq939XHd/OfA24l6Ubjez\n",
"U7oVtNx9KfCXnLFYnngK0HKA4TYuAfYPwmijIcZTyXYWymn8VS0fqF5OVctnXJnZaWXHICIi1dG1\n",
"QGNmG5jZMWZ2IfBZ4idSTwV+QvuBevuxAvi6mX3PzHYtcLvN8nRwAUCjXrsfuBY4YIjxiIjIkJnZ\n",
"s9u9gJ6rwIuIiLSTpxfBZcAFwAfd/arM/G+b2duLCsTdjwcws32IGyG/stPyZrY4vbubtlPIOR3c\n",
"ef7pl9p7Z/Kuf+HK39/692b2WJ/7yzP9LuDqIW6/jOl93P3TYxTPwNPpvHGJR/nMnW7Orex4lM/c\n",
"6RJFxIPLW4v3cnXkZGZfAfYEVgL/NeTq7CIiMqHMvXMP5Ga2ibs/2Oa9vdz9+i7rHwQc4e71XAGZ\n",
"7Ql82N2PbrfM0qVLfWZmptVJsqMgjDYj7o5380a9tjrnOgHwX4167em97i+vbGGxKpTT+KtaPlC9\n",
"nKqWD/T//T0oM7vc3QeqjWBmXwY+5O53dFqurBxFRGQwRX1/52mD1bJwlbzXrXD1PuBE4Agz+2Jm\n",
"/lFm9pKmZc8ysx8RP73KVRjrw7OBa/IWrhJXAdsEYbTjkGIahzu7hVNO469q+UD1cqpaPiX7UAHb\n",
"GFVbYRERmWB52mC9scW8XFUD3f1kd1/s7k9z97dk5p/j7uc3LXusux/k7ke4++/ybL8PewK/6mWF\n",
"Rr22BjgfeMVQIhIRkaFz96iAzYyqrbCIiEywPHfi3tBi3lFFBzIiu9PfoMFnA8cUHMta2XYWVaGc\n",
"xl/V8oHq5VS1fMaVmW2cZzl3P97dXwj8K3MHLBYREQH6r+rwhEKjGJ1+C1iXAE8LwmiHguMREZER\n",
"MrMjzewmM1tuZg+Y2Qrghh43swp4bAjhiYhIBeQpYN1rZi9KJ8zslcCfhhfSUPVVwGrUa48C32FI\n",
"T+6q2M5COY2/quUD1cupavmMiY8T10g4EagBbwZyjYM1orbCIiIy4fJ00/5u4LtmdgLxk6sNmcD2\n",
"SEEYbQDsANza5ybOBj4CnFpYUCIiMmr3uPvVZrYI2NXdzzSzS4GPdVvR3Y8dfngiIjLp8vQieDuw\n",
"L/BG4G+BZw+xE4ph2hm4M3ka1Y8I2CUIo50LjAmoZjsL5TT+qpYPVC+nquUzJu43sw2AK4G3mtmh\n",
"wJNLjklERCokVxssd1/j7jckrzXDDmpI+m1/BUDStfu3mdwOPkREJO6gYv1kLKszgLcDb+m8ioiI\n",
"SH5dqwia2ZbAkcDm2fnuPmlV5QYqYCXOBkLgE4OHM6uK7SyU0/irWj5QvZyqls84cPfrMr9/Gfhy\n",
"ieGIiEgF5XmCdRFwILBp02vSFFHA+hGwfRBGiwqIR0REREREKiZPAetBd/9bdz8x+xp6ZMUbuIDV\n",
"qNceB75JwWNiVbGdhXIaf1XLB6qXU9XyGQdmdqWZvcnMNik7FhERqaY8BayrzGyPoUcyfEU8wQI4\n",
"E/jrIIysgG2JiMho/T2wG/BLM/svMzug7IBERKRa8hSwngksMbPvZV/DDqxIQRhtAmwB3FnA5n4K\n",
"zAeeVcC2gGq2s1BO469q+UD1cqpaPuPA3a9z9/cCewD/A5xsZrlvvpnZBmZ2m5m9bWhBiojIRMsz\n",
"DtZHhx7F8O0GLGvUawP3gNio1zwIo/8BXgf8cuDIRERkpMzsicBrib/H7wf+rYfVjwOuGkZcIiJS\n",
"DXnGwfpRq9cogitQUdUDU/8DvDYIozwF1K6q2M5COY2/quUD1cupavmMAzP7AfAT4p5xj3L3w9z9\n",
"rJzrzgdeBHx3iCGKiMiEyzUOlpntZGaHZ6Y3Hl5IQ1FoAatRr90E3AHMFLVNEREZiVPcfQ93/4i7\n",
"397juscDnwPUBldERNrqWsAys9cRd+zwsczsC4YW0XAU/QQL4L+BvyliQ1VsZ6Gcxl/V8oHq5VS1\n",
"fMaBu0f9rGdmC4ED3P3CdFZxUYmISJXkeYL1NuAg4C+ZeZN2YhlGAets4IikAw0REam2/YENzezr\n",
"wFuBvzWzvUqOSURExlCeAtZqd380nUjGDpk/vJCKlXSnXngBq1Gv/QH4MfDKQbdVxXYWymn8VS0f\n",
"qF5OVctnkrn7+e5+qLv/H+DzwFfc/fqy4xIRkfGTp4D1MzP7OLDQzF4GnA98Lc/GzWx/M/u5mX0i\n",
"x7IzZna5mV1mZrU8289pK8CB+wrcZirtTVBERKaEu5/h7p8vOw4RERlPeQpY7wduBW4jLkyc5u6f\n",
"yrn9DYGTui1kZgZ8GDgUOAw4Ief289gduLlRr3mB20x9D9gvCKNtBtlIFdtZKKfxV7V8oHo5VS0f\n",
"ERGRaZCnm/Y17v5Fd3+1ux/j7mfm3bi7L2Xdtlvt7Abc5O6r3H0lsMzMFuXdTxfDaH8FQKNeexi4\n",
"GHjZMLYvIiIiIiKTpZBxnAqwFbDczE4l7kBjeTJvWbsVzGxxenc3bafQZnr3h+648TGzmbzL9zp9\n",
"7mMP3He8mf12gO29C7h6SPGVNb2Pu396jOIZeDqdNy7xKJ+50825lR2P8pk7LSIiUnXm3rnmnJmt\n",
"IG7DtA53X5hrB2YHAS919/d2WGZ34P8R98xkwGnAR929ZQFr6dKlPjMzk6snwyCMvgV8o1GvnZ1n\n",
"+V4FYbQZ8ZhY2zfqtRX9bCNbWKwK5TT+qpYPVC+nquUDvX1/T6ppyFFEpIqK+v7OU0VwU3dfmL6A\n",
"FwMf7HE/3QJdRlxNMF12UbvCVR/2BG4saFtzNOq15cBPiY9LX6p2AQXKaRJULR+oXk5Vy2dcmNlO\n",
"ZnZ4ZnrjMuMREZFqydPJxTrc/afMFoY6MrP3AScCR5jZFzPzjzKzl2S2uSZZ7hLgouT3gQVhtB6w\n",
"C0Nqg5XxHeDIIe9DREQGZGavA84EPpaZfUFJ4YiISAV1LWCZ2bObXi8G9suzcXc/2d0Xu/vT3P0t\n",
"mfnnuPv5Tcsucff93f0Ad1/Scyat7Qzc06jXVha0vXa+CxwehNEG/axcxbFulNP4q1o+UL2cqpbP\n",
"mHgbcBDrdsCUqzqImX3UzCIzW2JmOw8lOhERmXh5Ork4pWn6z8TtpSbBUKsHphr12j1BGN0ELCbu\n",
"VVBERMbTand/NB4dBMxsE2B+nhXd/Z+TdV5IPITJWzqvISIi06hrAcvdDx5FIEPyNEZQwEqcS1xN\n",
"sOcCVhXbWSin8Ve1fKB6OVUtnzHxMzP7OLDQzF4G1IGv9biN/YAbCo9MREQqoec2WBNmJE+wEt8B\n",
"jgzCqOrHVERkkr0fuBW4DXgdcJq7fyrvymb2I+A44OtDiU5ERCZenjZYlyZ1zlu+RhHkAPZkRHcZ\n",
"G/XazcR1+oNe161iOwvlNP6qlg9UL6eq5TMO3H2Nu3/R3V/t7se4+5k9rn8Q8NfAGcOJUEREehGE\n",
"0fwgjMZqaIw8bbB+DtzPbNW3o5Kf5wwlooIkB3qUT7Agrib4GuDKEe5TRERG615ajA8pIiKlOBK4\n",
"Bri+7EBSeaqzPcvdP+buVyWv9wPPS6eHHeAAtiY+Af5phPv8BnBUr6XoKrazUE7jr2r5QPVyqlo+\n",
"48rM3pxzubPNbCnweeDtw41KRER6sGHZAWTleYL1FDPb2t3/CGBmmwFbDTesQjwNuLFRr43yLuN1\n",
"wCPAc4DGCPcrIiL9ex3wH90WcvdjRhCLiIhMuDxPsE4GrjazM8zsq8D/Ap8YbliFGHX1QJLC3DeA\n",
"o3tZr4rtLJTT+KtaPlC9nKqWT5nMbIWZPdDitQJ4ftnxiYhIdXQtYLn7V4mfyHwreT3P3f972IEV\n",
"YGQdXDTpq5qgiIgM1dXuvrDFa1N3X7/s4EREplEQRi8Mwmi/suMoWp4qgrj7PcB3hxxL0fYElpaw\n",
"318BK4l7E/x5nhWq2M5COY2/quUD1cupavmU7DtlByAiInPsCDxOxTqIyzVmk5m9zsxOTH43M3vB\n",
"cMMqxCgHGV4rqSZ4Dj1WExQRkeFx91PKjkFERFqqXK2vPONgnUI8av1hAO7ujHkbrCCMFgDbEg8k\n",
"WYZvAEfnrSZYxXYWymn8VS0fqF5OVctnnJjZ+mamgeFFRKRweU4uz3X3dxBXe0uN+/gfuwPLGvXa\n",
"6pL2/2vgQeKCqYiIjAkz28fMGsA9wO/N7BIz263suERkVhBGTwjCSDdAZGLl+fDOM7P1SApVZrYr\n",
"OdtulWjkPQhmZaoJHtVtWahmOwvlNP6qlg9UL6eq5TMmvgz8s7s/0d23Af4dOCPPimb2BTO71Mx+\n",
"aGY7DzVKken2CmASmqOItJSngPV54BLgqUl1wUuBDw01qsGV0v6qyVnA/wnCaIOS4xARkVmr3P2i\n",
"dMLdzwVy1XZw9+Pc/WDgROC9Q4pPROJBYzcvOwiRfuXppv1rwNuATwPLgIPc/eJhBzagUp9gATTq\n",
"tRuIexTsOjBlFdtZKKfxV7V8oHo5VS2fMXGhmR2QTpjZS4AlPW5jBfGg8iIiMh7GqvlSrvqt7n69\n",
"u3/W3U9z91t72YGZzZjZ5WZ2mZnVuiz7FTO7wswiM3t9L/tpUtYYWM0+DbxbY2KJiIyN9wJRZpDh\n",
"84B6OhBxzm28EfjC0CIUEahgz3KjFITRnkEYTdMg6mP1eenalsrMdnD3O/rZuJkZ8GFghjjxi4Co\n",
"wyoOHN3v/iBuGAnsBtzc7zYKdAFwCnAAcFm7harYzkI5jb+q5QPVy6lq+YwDd99kkPXN7AjgJncv\n",
"uxq6TLAgjDYFVpbYGZdU367AQuCKsgOZRnmeYJ03wPZ3Iz4RrXL3lcAyM1vUYXnLGVMnuwO/b9Rr\n",
"Dw64nYE16rU1wGeAd5cdi4iIDMbM9gUWu/uny46lmyCMtg7C6IllxzGugjDaMgijrUoM4QjgWSXu\n",
"fywFYTQThFGn60SZMEEYvXZENbkmrorgqgG2vxWw3MxONbNPAcuTee2sAL5uZt9Leivsx/MZr9L6\n",
"GcABQRjt0m6BKrazUE7jr2r5QPVyqlo+48DMTjOze5MqgSt6rBp4DhAkPQl+ZphxppLuqrfrY9VD\n",
"gEODMNq/6Jh6EYTRoiCMghHta70gjJ6Xc/EXJa8yqROsubYBti87CClcTwWsIIw2CcJoi2EFMwp5\n",
"Clj/aWanmNmW2VfO7d9H3AvMB5LXFsm8ltz9eHd/IfCvwCc7bTh74WFmizPTz3vgt1f/qcP7I51u\n",
"1GsPPXzPby9eee/vPtFueWCfcYm3qGlgn3GKR9Oa1vR4TI+BpwF7uPummdfCPCu6+y7ufpC7H+zu\n",
"7xxynCQ35vYAFg+wmR2Kiaa7IIzmBWG0ddPs3YGuTySCMNoiCKP1k9+f1Ocd700BdZ8/+QZ+EpE8\n",
"NSnzCaX0pvn//WDgxb1uJAijBUlTodLlGc/qg8nPV2XmOdD2iUzGMuJqghAfvEXuvizHequAxzot\n",
"kG2b0NRO4fkLd93nDe7+izbvj3x6wXa7vg+4JgijzRv12v0tlv9003Sp8RY0/cMel9e0pgeebtVm\n",
"aZzi63W6avmMiS8C55rZddmZ7n58SfF0sh/wx7KD6MHOwHOBM/OuEITRQmAl8cXUb4BfEF9cXQT8\n",
"ucf991Io8+blkwLe/Ea9lveJ5qAsCKMacHWjXus117GSPDlc0KjXOrWzzyut2jpotbJN6XBTH+Kb\n",
"AklzjrGWDLq8e6Nem9S2n93+lhaEkSXjyEL/zYVewez3SKm6FrDcve+7Qe6+xsxOJB5Hy4nHDgHA\n",
"zI4CHnL38zPzzgK2I64q+A+97i/5ot4ZuLbfmIehUa/dEYTRt4F/Af6x7HhERKbYR4DPAveXHUhW\n",
"EEbbA3dnLjCa398QoFGvFdY9fHKnd0tgm0a99usCNtnPRdFLiW/G9rx+8oTvtgIvkJ9NfPO4bQEx\n",
"qe54N7C8Ua89GITRRsC+jXrtJ33uc1vi656hF7CCMNpkiO3TtydnlccgjLYE/tLusw48RDwO1lAl\n",
"T1sPoYcbAiVaSNxmr5cC1kAF1CCMNgMe6PB3yrONXmLYHbgp+X1Bv/tkBJ+dPAbtUKIrd1/i7vu7\n",
"+wHuviQz/5xs4SqZd2xS/eIId/9dH7t7LvC/jXrt0UHjHoIPAK8PwujpzW+MWfWZQiin8Ve1fKB6\n",
"OVUtnzFxFvAAcF3Tq2wHEl9EtXMEfbYZCsLoSS3mPQM4mvi8uXc/2+0xhnkdqmx1vDAPwuilbdbd\n",
"j9mnHTD4E4/1O8SwYVKYWkT8t3pZ8tbOwI4D7rdvQRjVgjDaN8dymzEbc9kOIy5YtnNPQftZWzAI\n",
"wminFlXHBrmIL0UQRi/vYfFNB9zdSxhhFWMGj3estC1gmdnZmd/fMZpwBvY84GdlB9FKo177A/ET\n",
"vM9qXCwRkdK8EHg98RAa6atjm99+BWF0ZJ8dVLSyATC/z3UPDsJo46Z5eyU/2xYqUgW1adiZ9gXE\n",
"9JyYvVN+WLJvIy54blNADN10OjcfRnzB2WzP9JcgjDYIwui1fezPk/XnB2H0V9kFknmd4toWeEqO\n",
"fQ27XUqvHXbMg/ipWov30nZ8RV4rPR94cp4FgzDaKAijwwrcdxEWJz+b/48LEYRRu7/foH8Da/rZ\n",
"LY5DC/i+GYveBDs9wcp+EF857EAKMrYFrMQXiKtjHJ2dOYZtFAamnMZf1fKB6uVUtXzGQdJBRfOr\n",
"NqTdzWf2YjFPdZlctUpaFJYG0fFiJAijbWg6Z/WpU+6d8t6joH0Maj6tqx6lhSMDNutxm83xPhV4\n",
"RtO8I4l7Iu7UpXvhT2KSJ15FbWu9Dp/9lyVjgmX1XZgOwmiHIIxe0eH9XdsU6rI2J75WK0QQRts3\n",
"F5yb3n9tEEZP7bKZfm+u5PXqIIxaFUDXqX6bPMnt5SZCKs//phE/kV7n/yw5Pr1U+xv7AtZYBJhX\n",
"8s/7PMari/Z1JAMKvh34ZI5/cBERGQIz287MnmVmzzazfc3spSPa9bFBGO0A8TmrxZ3avZP3NgzC\n",
"qNVF8xOSi5u11YSaL1yDMAqCMGpVKGk+p6fn/27tl9bGEYTR5kEYHdxh2a5ttoMw6vVCMT1XDlx4\n",
"CsJo9yCMeu6ZrN3+m7qf35W4PU9aoMhTWG7+m7TLcXsyT8pG5CVJu8B1JO2nmud16/TsKKC5eUQ2\n",
"13YFqtx/8yCM9kqqcG7L7Ge21XXsc4l7Es29/eQCP/2/7efpyl+RFJyDMNo4CKP9k14zs5+RzfvY\n",
"bhqftfm+SN+fU4U1CKM9gjDaNvk9/ZvmKcSsfdIVhNHL054/O2j1f9Oce1tBGO2U/DrsAmbhOiX4\n",
"FDN7j5n9I7Bj8vva16gC7MEi4OFGvXZ32YF00qjXfgxcCnw8nVfFdhbKafxVLR+oXk5Vy2ccmNn/\n",
"R9zD1DeB04ALgDfkXHd/M/u5mX2i+9JrPb3p7mvazurZzH0ytHFy4XEIcW9Y0OE8HYTR3sQXr1mL\n",
"iBuL59XLzdQnA3Pac2V06hY73U+rqkgtqxH1UZ3+eV3WezLxcDHQ+rj2ur9s9/MbZX5/FXGVNNJ4\n",
"2lxQNl+c9jpWUNtxgtodgyCMDmzVJq+NA4Mw2rTp4v2w9MI8o+3TrkyMnW4qP7fN/I4X4U0F2WfS\n",
"vb1QekzWb/HUrJvNk6dMa/9ng3iMuqcnv2+Vs/D18iTOFxM/sSzCIma/L1ppleuzgfTJfVodstV3\n",
"QfPnKDu9Met+7jsKwmj95In4i+l8MyYbR/p/dHiO49vt8zKvzVO6oegUzFeJ/yibAP+d/J59jZtx\n",
"G2C4k3cS14l/V9mBiIhMmVcTX5CcCrwHmCHu9CKPDYGT+thnq7vLrTq02Iy4wJR9r1OhZXPip1rN\n",
"F5atzu3ZBv/ZeDbNzO92AdPqbvTGQRjt0zRvfhBGT2xaNE9BrtOFvyXb3q2psJKNKc0l90VfO8nT\n",
"rk4dMaS86SfE7Z2yBYoDgcNbrJtuv+caQ0nBpdPTuGODMNqrxfzt6a3jgiOI/0ey1lapTf4Wnaoo\n",
"poXQjoXHNk82Nww6j2f1ajIFWeLjmCe3pxLn1au1f9OkCuVOwN5J7C8C9mr1hI/2f9/1MoXdtp+B\n",
"Fk+pmwulgGrIAAAgAElEQVSrGyXz+34Klnh+m89Mut8Nmf0/PDTnNtPYN6T1Zykr7dm83d/w6CCM\n",
"DgvCaOs2T866VSfeDjioyzKFadtNu7uf2O69MTXu7a/WatRrfwnC6HDgp0EY3enu3yw7pqJVse1I\n",
"1XKqWj5QvZyqls+YuN3dV5rZbcAz3P0LZpar+pW7LzWzXCfo5gvGIIzSC4v0gmN1ulyjXluZWTRv\n",
"1Rljtnrf/kEYXdCo19Ku59sWlJI7yHMucpLC1dFBGP24Ua/d0Wb1VhfJOxFXubo6M+95wJOCMDo7\n",
"M69TVcR2Obfa33OS19qutYPZsbRSGyRPGx5v1Gu/ad5Asnw3+wL3BWG0Cmi0iSWr+QI5m++W9FHo\n",
"C8LomBbztgBeQNP1ThBGW7YYS6tdnu2ebu0JLEuaM2Q9IXm/1ZOqp9O5F8V2g0yv31RwODIIo7Pa\n",
"LLsXcHmL+fOYm2P2ODfn0Zx32tHGtsDjxDcsAuJaRs1jZG1MPIRQWqjMdniSFpR3A/4qCKOzmro2\n",
"b1d4WkD8WV4rKXA90KjXHk6mjwGah1DYAsh2t5/mdXgQRlcDtzfqtYfa7LObbYHr27z3KuDi5PeO\n",
"Y5UlTwhXMzum7bNY92/T6vtgp+TnPi3eS21J/IT/NuCKIO5EqHmss7Rd5AHAjY16reVYgkHc++Z1\n",
"jXrt0SCM1gP2btRrv+yw757kGWh4UjyP+EnbRGjUa7cHYXQEcHEQRr9Pqg6KiMhw3WVmWwKXAT82\n",
"sx0YTg9r2e6U12duO5O057cdgJv72H7zhc3hzBY6WralSJ4qtbuDnFZ96tqBRhBGT2nUa3cmk+lF\n",
"avauc1pgyrZRSi+oWl1sdrqLv3bg0TZPlIx4LK3fZualF79riAcdbd7W/umM5CL/WOAq5h7T9MnJ\n",
"rS32m0qf3sxpR5RciBvrPj18Iu3HYEufDixMBjtuVfDcitYFp8OYO57TmmR7XaseJjcEngU8EIRR\n",
"c3zp+s9tWmc78vceuFMQRjdlCoHPZ24e7eLst/1d85A9zftLb4Lsy7rVHLNPXdKC+87Mdh/f/H3x\n",
"zKbpdf7mqWBuT32tnhalbRzPDMJoZ+LPQPPnfmEQRs8hHpbo8ab39iH+m1wD/I74aV32KfV6zK1W\n",
"vE7syf/yqmT6mcDtmfebP5MvDeJxXjcEHib+jDyU5PYAswWy5ppvq5O2Yf122LNTUphc3GGZpwBP\n",
"CcLobuICYXOhc3fgDuAPxN9BewRhdMfHn91nRE2GPg7WKARxj0p7AIWVPEehUa9dDfz1mtWPfS8I\n",
"o/27rjBBqth2pGo5VS0fqF5OVctnTBzv7n929xXA64A/EffUVrTs+bVTtZh+zakO1OVCel6rdfJK\n",
"al2k2z8g81baO1r2HJYWNnZpsVx2m83xrmya3ojZ9mS7sm7VtLRtzw6Z95u1u8bJ3kk/Nvm5FTm7\n",
"8W6jVXf8zyP+bGUvtg8lvqDPar4Yb9npSlIISgutecZES7d7LOsWhHYNwugpyTZrQdxxyQaZdZq7\n",
"o9+ozVO/xTS19wvC6AVB+x7zDmuqQtdc1a1d9bTts22mgjDaMlNNbvPM9VPz57v589UcV7dOHV7Q\n",
"NJ3+jZvbTqWF7LSq8U6ZWNOe8SCu0tjOOrEHYfQC2nfGsTfx07L0GDQX+PZK/qbp53/nIO56fkc6\n",
"t4VL97U/SYct6fKZ/9VW/1ObEH9mn0l8bNK/Y7bqaPPTzzXJ8q2eVuW94dXue3vnpmqUTyb+fG8B\n",
"azstSWNLv5vT/5XCeiKtyhOs5xA/5itshPtRadRrS3Z69XtO2vp5R3w7CKO3N+q1b5Qdk4hIVSUF\n",
"q/T3fgcZHuQk/AzgV9kZbdoTdHMAc6vGHNtqwcTLgbtybNeDMNqd+ILkemZz3ZxMtZ4gjHZp1Gu3\n",
"5A937QXWvsClSXWz5gv53zVNZ99vbqOTtj9qVwUtjXNr4OmNeu2HzF64tbqwzravaW7L0c/fe0tm\n",
"L6zXS7abFuza9bq3dsDnFm3YoPNA1ARh9Dzg3ka9lj5xW0RcvRFmO/dIHRCE0Y3MPh1Jn5zMo/XY\n",
"aE9m3Sdx7a4fn5q8/whJlbom2Sedzce1uYv6rL2CMHqozXLpNndrmr8oCKM/0f5aN/07NxcaDkx+\n",
"tuu5LmgzP7UfkP5vbN1pwYwtgnV7+8sW4tp9/uYHYfQgc29MwNz2jM8hPk7N3a43V11tV7hJj2Gr\n",
"5ztpYb/5Cd16wGvabA+GO6xCp88SZG6yJDcb+m4P2U5VCliHAhNbxe62b536ySCMlgDfT/7BTmmq\n",
"vztxqth2pGo5VS0fqF5OVctn3JjZeu7e3E6j0/LvI66Kt62ZLXT3t/Sz32Du4MOdLkLaWUD89K0X\n",
"c7rdbiF7AfUE1n2qkH36sF/QovvnFprvXD8pecLRqprcc5oKN/0UPJulXae/mLmFjKxsgaaXXhjb\n",
"mdfi93ZjirY637d6mtaut73Uzsl6rao0tnpSkG17mBYm2l30No/D1e0z+5wu70NvVXO7dQXfyk7E\n",
"x6Nb7abm6mvp32vTFgWQTtYpTCWF0LxPrxcQD4LeSqvCNsT/zwcDzW3vWkkLoc2FyeYnv3M6FEna\n",
"Z6Z/705PwXsZO2w9hjRocmKnLu9nP+fZJ/LNVS77NvFVBJNeTf4e+FLZsQyiUa9dQ/w4+vXAaX3e\n",
"0RQRkRbM7O1N018C7jGzX5nZ09qstg53P9ndF7v70/otXCWyJ/Tm6mK9yFPAGUTzE4FmrarF5fEM\n",
"1j0GWcPqpbhV4arVnX+YW41qmHfaAQjCqPkz2K4qVjfrDAQb9D4obN6aQP0ek+y1Td7u4gexAUn3\n",
"/X1qVyjuKKm22KmtUxF+n/wcZFDkvbsvknug8Tydx6QKHxy7R+0Kzoe1md+ziS9gEVeJuKZRr91Q\n",
"diD9SttZJD037U/cMO+CoMMYF+Ouim1HqpZT1fKB6uVUtXxKtrZHNjP7G+IL+W2T+aeOOJZhdKox\n",
"DFW/0Ze3psiwj8Oz6dxz2ih1GhagCNnxjybl/6AfvXSF3692T7wmQdk16A7svshgJrqAlTS6eyfw\n",
"mbJjKUrSc9ArgGuBnyV14UVEZDDZE/q7gH9y9zXu/mvKv5sq5ch7DdStzc0wDP2pWRvNVQGLVvaF\n",
"dZVMcgG125hVE2/SP+j7E9fhvLDsQAbR3M4i6XbzPUEY3UBcyPoy8PFGvdZrffvSVLHtSNVyqlo+\n",
"UL2cqpZPyX5jZh8lbkNwvbtnux7OU/VKqmfgAYmHKNfYbCITqupPxyf7CRbx06vPZgaCq5RGvXY6\n",
"cbeiC4AbgzA6IQijfhp6iohMu+OIB738U/I7AGa2AAjLCkpERKpnYgtYySjtBwNfLTuWQXVqZ9Go\n",
"1+5u1GtvI+49aDvgp0EY3RCE0SeDMFrcoavUUlWx7UjVcqpaPlC9nKqWT5nc/WF3P9HdT3D3h5rm\n",
"n1VmbCNye/dFxla2x8HmwWPHyZ3Ar8sOYkIU1ltbG3d2X0RkeCa2gAX8A3BGo15rNc5C5TTqtVsa\n",
"9dpbiLsc/Rvi8SVOAX4fhNFXgzB6TZtR7kVEZPQup3X19Wu7rPejAvYLcG9m3jeBnzbqtTOBPyTz\n",
"0u7prxxwf1lXdHjvAeBs1h3ragWwJDO9vMN2/wxc1KjXvgX8tEscS7q834vmscaa/RyIiMcM+yXr\n",
"jjO0NPnZqZvw3zRNX9xTdIO7uo91eq011OoYtivwd/vbZv0k+fmXzLwlwJJGvXZ5i+WLlP5NbyYu\n",
"9H+/j200ui/Sk+uBP/a5XitX9RnHGmAVczuR+U7T9Ko+t9/JcuJxBptvdDyW/OxW8P55UYEM9emH\n",
"mc0AJxAf5BPcPSpi2cRewPHFRFquXtpZJNUhf5G8TkzGIXk58HfA6UEY/Zn4ZHQj8QfpTuAe4i+4\n",
"+xr12mMtN1ywKrYdqVpOVcsHqpdT1fKZdL2cpxr12p0AQRhBPJjxH4jHxLme+IJwNS3GyGnUa3cH\n",
"YfQD4jv8jxKPgfMgcVfML0iWObOp++0HSLpIbtRrdwZh9DDxxdv6wF+axlX8KXFHBjc16rX7khjv\n",
"Jh4PZ2/iAt6ezA68CfBoo177VhBGryG+Mfsd4jF4VgGLM7HfFoTR84HziQcJvpr4/HM4cGejXluT\n",
"nKPSQVSXN+q1PyXH6LZGvXZFMsDuauBh4rFsntSo1+4HLsrs53fA75LxxNbuP3EV6z7xepj4Yivt\n",
"qOKn6XFsshq4pFGv/SU5JnsBzyS+iH95i+XPBR7JHNt7k/VuTI5LkOybRr12UxBG2fHFrmnUa9cH\n",
"YfRSYBmz3eFf2ajX7kuOR2pVktO9wKsy83/H7HG8haZxopLPyN7A04kLA09PtrUpcDdxb8WXExdy\n",
"9yG+AL2dueMiQXw9sbZNWKNeOzvz+Xu8Ua99IwijLYk/ry9L5v+Z+Hg+noyjtH7y3nrE/wPpuFsX\n",
"EtfQ2TLZ9u+CMLoHeDVxYelA4kLMnIFjG/Xa7UEYvRD4caNee7BF3PczO27TfbTuIfEKZgt76xEP\n",
"zrwt8ENaD9x9WRLTPckNC0gKIsnf7WHi/93m4QZWJfN/BjzQqNdWJevck7z/xCT3jYAas+PH/Za4\n",
"WvN+xMfqxZltXkP8Gf09cZf3txL/Pbcm/ts+1qjX7s38rf5I64GPr23Ua9cEYbSA+HM0D1jWqNce\n",
"DsJo3ySnBcAlxN9Fa4g/iw8mx2kDZgcbBjiP+DvyceK/42pgvUa9trLpsx2x7mDi5yX7Snt8XUl8\n",
"w+ER4s/Ozck21xB/nuYx+7n4HnB4o147P91YEEa/AY5MJv+X+HP2QGZ/jzA76PjPgTsa9dqjS5cu\n",
"pQhDK2CZmQEfJj6BGPGXY8uTUS/Lphr12hFFxjupGvXa7cDngM8FYTSP+Evw+cRfkjXiL9HtiL9Y\n",
"tkhOvHczW/i6m/jEfy/xP9+fM6+Hqtq+TUSkF/2cpxI3ALc06rWHgfSC7G6AIIyuIP7+3ZP4ScYq\n",
"WNubbOqeZNmHk+nrkp9nk1yINOq1HyQXUWuS6e+2C6ZRr62k6SlBo15blXSq9FCjXrsbuDsIo4XE\n",
"FyCvYrZwswRY06jXHmV2MNszk3PPZsm2zkzi/VayHEEYfZvZO8g3EV80PpHZJ1Znp+eaps6cbk5e\n",
"7aTnp7OJL+DWVh8MwuhSYEWjXnsomX4QODi5iP0j8UXh4cSFk1uJC6LZG5A3JPFlC6ffB44AvpXd\n",
"V1ajXlsNLAvC6IFGvfZgEEbZp2lXNuq1WzLL/iCJ7XxgZWabVwP7ZC7g05xuJi6c7E/8ubiX+O9x\n",
"K8mTyKSAu1Oy/WuZfWJ6U3Lx7slx2zAtlARhdBNwczL98yCMNiMuLB0JrG7Ua/8bhNEtwPxGvZaO\n",
"vXQvsG2jXvtGsq8/JwWpe4hvAD+WdNiVdtz1eBBGlxAXEO4BHgLuTgq0FwVhdBTJ04VGvfZoEEbn\n",
"N+q15cC3kxjvTGI6EPhhWkABLm1TuIK4QLKIuAB1cDJvDfH/7SbE/2+/zxSSH02O7/XJPm8l/l+9\n",
"j7jb9Wc16rW7gjBKC6bNfgn8sVGv/TlZ/2Di/6HriQuic9ZJP5/J8SDJ6/zm5YIwui25QXEesE1y\n",
"rFYDz2zUa5dmFn0gCKPnAn/IfJ4uJi4ALW3Uax6E0cbJcbkZmJfmn3xH/app15c26rXfB2G0XvLZ\n",
"TuO5ML0ZkTgzuTmycbKddLnziD9D6d/reuIHJL9s1GvLMwWu7PcFwK+Tz2+q+elXuv07iI/tQ8RP\n",
"6tdKCnTfB2qNeu23QRjtRlwY3Yj4ptH/Jt+bNzXqtd+22v4gzD3vMBA9bthsd+D97v53yfSXgZPc\n",
"fdkgywIsXbrUZ2ZmyurCtHBmtngUd6qTbu03J65m+BTiL4wnEf+zbpv83IL4TtIWxHcsHia+S/Eg\n",
"8RfKA8nP5cR3h5YTfzGsIv4ieQRY/Zdf/2S3LZ7+wuuY/TJPP2jzklfavWj6d0zvdqxJfmbX88yy\n",
"lll/XtPLM+utyWwvW0hst410fhpHdt01gN931ZKnb7XvoemXz3D+cUboT1dd/FdP3PdFzV+mE63P\n",
"nNK/ffp7+rlI77qmr/Wafk858Ynu0eT1GPHnZzWz7QzSzxxN66Wf0/Tz5tlXjnyMdT/TadzrNeWR\n",
"7jv93DaSwc1HblK/v3s5Tw0rx+SJzb1pYSQpBHmjXlsRhNEriWspXFb0fsdVUrDbrlGv3dXn+gcD\n",
"V2Qu/tottw/xue93xBfZv+hnfz3GZk1PHnOvB1gRN0eTsTgfzRQCBhaE0bHAd5OC/sgkhcsNk0Jb\n",
"39vo9lkZV0EYrT+qGkxd4pgHbJQWwpLvtBXZQnLSt8Dj/Xz++4hnIfBg9v+lqO/vYRawnk88+rMz\n",
"e/FylrvPqe/dy7IQJ198xCIiMgoTWsDKfZ7SOUpEZHIVcY4aZhus+4iflryV+GR0Gu0bi/ay7ESe\n",
"nEVEZKLlPk/pHCUiMt2G2YtgtuGmAYvaVfnrcVkREZFR03lKRERyGdoTLHdfY2YnEvc64sCJ6Xtm\n",
"dhTwkLuf321ZERGRsuk8JSIieQ2tDZaIiIiIiMi0meSBhkVERERERMaKClgiIiIiIiIFmagClpnN\n",
"mNnlZnaZmdXKjmdQZra/mf3czD5RdixFMbMvmNmlZvZDM9u57HiKYGYfNbPIzJZUJScAM9vAzG4z\n",
"s7eVHcugzOwrZnZF8nd6fdnxFMHMtk/yuczMTik7nkGZ2cLkuyFKft5fdkxFq9o5qletzmntjkmv\n",
"86uk1XlSx2ldrc67OkbtNZ/PzewQHatZra4Rhn6M3H0iXsS9Nv2EeATm+cBlZcdUQE4zxCOlf6Ls\n",
"WIaQ28HAaWXHUXBOLwS+WHYcBeZzPPAt4G1lx1JALl8Gdig7joJzOhN4ftlxDCm3vYHTy46j4Jwq\n",
"d47q4xisc05rd0x6nV/VV3Ke/LyOU8dj9ELgCzpGXY/T2vO5jlXL47PONcIojtEwx8Eq2m7ATe6+\n",
"CsDMlpnZRHeT6+5LzeygsuMYkhXAI2UHUbD9gBvKDqIIZjYfeBHwDWCTksMpgjFhT+Q7MbN5xN2A\n",
"X1F2LEPyDuCzZQdRsMqdo3rV4pzW8pgQ/6/mnl/hY7gCeBQdp072A25Ex6itFudzHau5mq8Rhn6M\n",
"JqmAtRWw3MxOJT5Qy5N5VfoAVMkbgc+UHURRzOxHwHbA/mXHUpDjgc8B25YdSEFWAF83s/uAd7v7\n",
"b8sOaEBbAxuZ2bnAQuBz7n5uyTEVwsy2JL6TeG3ZsRRM56i52h2TeT3Or+oxTM+TOk4tNJ13d0XH\n",
"qJ3m87k+T3NlrxHewwiO0SQVsO4DNgfeSpzcack8GTNmdgRxSf/GsmMpirsfZGYBcAbw4rLjGYSZ\n",
"LQQOcPeTzez/Ev8/TTR3Px7AzPYBPgm8styIBnYfcD/wauLv6Z+Y2YXuvrLcsArxZuD0soMYAp2j\n",
"5mp3TOb1OL9ysudJM9sdHac5ms6770DHaI4253P93zVpukYIgfcy5GM0SQWsZcSP9CBOrkqPLyf+\n",
"AjdlZvsCi939n8qOZQjuJR5gdNLtD2xoZl8HdgGeYGaXuvv1JcdVhFXAY2UHMSh3X21mdwDbuftd\n",
"Zraq7JiKYGZPAF4GHFB2LENQ5XNUr9JzWstjklSBzT1/VEGPSovzpI5Te+l597foGLUy53wOXI6O\n",
"VTvpNcLQP08TU8By9zVmdiJwCfE/24klhzQwM3sfcDiwrZktdPe3lB1TAc4B7jCzS4Fr3f2dZQc0\n",
"KDM7G3gisBJ4e8nhDMzdzwfOB0h609lk0gtXZnYWcVWSFcA/lBxOUd4PnJ7coTynIk+vXgmc5+5r\n",
"yg6kaFU8R/Wq1TnNzD5M0zFpd6ym6BjOOU/qOK2r+bzb67GYhmMEbc/n1+rztK7ma4RRfJ4s6R1D\n",
"REREREREBlSZXrdERERERETKpgKWiIiIiIhIQVTAEhERERERKYgKWCIiIiIiIgVRAUtERERERKQg\n",
"KmCJiIiIiIgURAUsERERERGRgqiAJSIiIiIiUhAVsERERERERAqiApaIiIiIiEhBVMASEREREREp\n",
"iApYIiIiIiIiBVEBS0REREREpCAqYIlMKDP7ipl9ODP9VDNbY2b6vxYRkVLpHCXTTB9ykWrxsgMQ\n",
"ERFpQ+comQoqYIkUzGJ1M7vJzG43sz+a2crk92fm3MZBZnaHmR1nZr82s9+b2ae73PmzglIQEZGK\n",
"0jlKZPjWKzsAkQp6F3AscIC7/8HM5gM3Age7+y09bOfJgLn7081sW+Ay4Gbg85ll3mFmr0t+X7+A\n",
"2EVEpNp0jhIZMj3BEinem4EPufsfANx9JXAl8Noet3OXu5+WbONe4DPAa5qW+ay77+LuuwD7Dxa2\n",
"iIhMAZ2jRIZMBSyR4u0INN8FvB3YfsDt/g7YdsBtiIjIdNM5SmTIVMASKd5twG5N83YHbu1xO09o\n",
"mt6V+CQoIiLSr9vQOUpkqFTAEinefwAnmNl2AGb2auAA4H963M6TzOwkM5tnZjsC7wbO6LKOGhGL\n",
"iEgnOkeJDJk6uRAp3r8RN+b9kZktIK6Kcai739Pjdu4C/kTc+Hgr4D/d/czM+626u1UXuCIi0onO\n",
"USJDZu6j/6yb2RuA44DHgH9x90vN7BDgQ8T/fCe4ezTywETGhJkdBPy3u+9YdiwiVWVm+wOnAj90\n",
"9/d2WfYLwB7Ed+Df4O69VqcSqQydo0Q6K+sJ1ruBfYBNgAvN7IXAicAM8cnrIkAFLBERGaYNgZOA\n",
"F3Rb0N2PAzCzg4H3Am8dbmgiIjKpyipgXQscAmxDXJjaDbjJ3VcBmNkyM1vk7stKik9kKMzsNcCn\n",
"aF1NwpL5R400KJEp5e5LkzvxvVgBPDKMeETKpnOUSDHKqiL4FuBA4k42zgT+CBxN/I+bNoA8y92v\n",
"HHlwIiIyNZIC1ku7VRHMLH8a8Bl3v3G4kYmIyKQa+RMsM9sVqLn7Mcn0pcDxwObEVS4MOA24r902\n",
"li5dqkaSIiITamZmZiJ7EjOzI4hrW3QsXOkcJSIyuYo4R5VRRXAesBmAma1PXLBaxuyYDAZ0rR44\n",
"qSfofpnZYnf/YdlxjNI05gzTmfc05gzTmfeYFj66nk/MbF9gsbv/U54NTts5SkSkCoo6R428gOXu\n",
"vzGzy8zsCuKT2qfdfaWZfRi4hLia4ImjjktERKaLmb0POBzY1swWuvtbkvlHAQ+5+/mZxc8B7khq\n",
"XVzr7u8cfcQiIjIJSunkwt1PIu65KTvvYuDiMuKZBNN2lxumM2eYzrynMWeY3rzHhbufDJzcYv45\n",
"LebtMpKgRERk4s0rOwAREREREZGqUAFrQpjZ4rJjGLVpzBmmM+9pzBmmN28REZEqUwFLRERERESk\n",
"ICpgTYhpbKsxjTnDdOY9jTnD9OYtIiJSZSpgiYjIxDKzrczsS2Z2UTJtZvb2suMSEZHpNfIClpkt\n",
"NLNLzSxKft6fzD/EzC5PunCvjTqucTeNbTWmMWeYzrynMWeY3rwLdjpwPrAhgLs7cEypEYmIyFQr\n",
"YxysB4CDAcxsb+AdZmbEY1/NEI+NdREQjTq2UQjC6K3AHxv12jfLjkVEpAK2dPdvNT21Uu0MEREp\n",
"TdknoXcAnwV2A25y91XuvhJYZmaLyg1taJ4BPL3XlaaxrcY05gzTmfc05gzTm3fB1pjZdsSD1GNm\n",
"rwSWlxuSiIhMs1IGGgYwsy2BHdz9WjN7PrDczE4lfoK1HNgKWNZh/cXpxUlazWZCpuev/P2te01w\n",
"/JrWtKY13ff0EPwTcAGwo5ldRVxV8BVD2peIiEhXFldXL2HHZu8HfpNU7dgd+H/AW4kLWKcBH3X3\n",
"lgWspUuX+szMjI0u2uIEYXQ2sLxRr725l/WyBbJpMY05w3TmPY05w3TmPYzvbzNbH9gTWE1cG2JN\n",
"zvX2B04Ffuju7+2y7AxwAvGTshPcvW019kk+R4mITLOivr9LeYJlZk8AXgYckMxaRlxNEOIC1qJ2\n",
"hasKWADkOvmLiEh37v4YcF0fq24InAS8oNNCSTvhDzMF7YRFRGRwZVURfCVwXnqX0d3XmNmJwCXE\n",
"dwdPLCmuUZhPfILuybTd5YbpzBmmM+9pzBmmN+8imdl7Ws1391O7revuS83soBy7WdtOONnnMjOr\n",
"8o1AEREZQCkFLHef04Oeuy8BlpQQzqgtANYvOwgRkYrYtGn6OcD9Be9jK3psJywiItOrtE4uptgC\n",
"kt6uejGNbTWmMWeYzrynMWeY3ryL5O7r1HhI2mOFBe/mPmBz1m0nfF/B+xARkYpQAWv0FlB+9/gi\n",
"IpXk7o+Z2TY9rtat2vY0tRMWEZEB9V3AMrPN3F1jjfRuPrBBrytN413uacwZpjPvacwZpjfvIpnZ\n",
"eaxbK2Ab4K6c674POBzY1swWuvtbkvlHAQ+5+/kwde2ERURkQIM8wbrGzCLgi+5+ZVEBTYEFyUtE\n",
"RAb3yabpv7j7tXlWdPeTgZNbzD+nxbxpaScsIiIDGqSq2iLgPOCDZvYLMzvOzDbJs6KZbW9mkZld\n",
"ZmanJPMOMbPLk3m1AeIadwuAjYIw6qlwmw7WOU2mMWeYzrynMWeY3ryL5O4/anrlKlyJiIgMS99P\n",
"sNx9NXAucK6ZPRf4BnCSmX0F+JC7P9hh9U8CH3T3K2DtGCMnUvExRoIwmkc87soDxD1f/aXciERE\n",
"REREpEh9P8Eysw3M7BgzuxD4LPFgjU8FfgJ8q8N684gbCF+Rmb12jBF3XwksM7NF/cY2xjYCVjFb\n",
"wMptGttqTGPOMJ15T2POML15D5uZnVZ2DCIiMr0GqSK4DKgRP4naz93/w91XuPu36TzO09bARmZ2\n",
"rpktNbNXkhljxMw+xewYI21lq9aY2eIJmV4APPz4Iw+vvuP7X5gZg3g0rWlNa3qk00Uxs2e3ewH7\n",
"Fb0/ERGRvMy95yGZ4hXNNmlXDdDM9nL369u8tx5wKXAQcRXFnwB/B7yHdccY+Wi7bnCXLl3qMzMz\n",
"3brVHTtBGO0I/Bi4G3h3o167ossqa5lN33g505gzTGfe05gzTGfeRX1/m9n9wC9p08W6ux886D76\n",
"Nfg4f0oAABynSURBVKnnKBGRaVfU9/cgbbDatrFqV7hK3lttZncA27n7XWa2iukZY2Q+8DCwgh6r\n",
"CIqIyDquc/cqd4gkIiITapA2WG9sMe/tOVd/P3C6mf0YOCdpd/Vh4jFGLqK6Y4wsAFbSRwFr2u5y\n",
"w3TmDNOZ9zTmDNObd0E+VHYAIiIirQwyDtYbgC81zTsK+Fy3Fd39duAlTfMuBi4eIJ5JsID4CdYD\n",
"wMKSYxERmVjuXrmeZkVEpBoG6eSilScUvL2qSQtYPT/BGkYj8XE3jTnDdOY9jTnD9OY9bGa2cc7l\n",
"Zizn+Itm9gYzu9LMfmxmpbXvEhGR8TdIAeteM3tROmFxb4B/GjykSptPn1UERURkLjM70sxuMrPl\n",
"ZvaAma0AbsixnhFXTT8UOAw4ocsq7waeT1z74qTBohYRkSobpIrgu4HvmtkJxE+uNgReUURQFZZ9\n",
"grV5LytOY1uNacwZpjPvacwZpjfvgn0cOIZ42JDLiDtM2inHemvHXwQws2Vm1qmDpWuBQ4BtgAsH\n",
"DVpERKprkF4EbzezfYE9klk3ufuaYsKqrGwbrB1KjkVEpArucferLR6cfld3P9PMLgU+1mW9teMv\n",
"Evdem46/2K6AdTnwf4lrfpxZTOgiIlJFA7XBcvc17n5D8spVuDKzr5jZFWYWmdnrk3mH5K0HP+HU\n",
"BqsH05gzTGfe05gzTG/eBbvfzDYArgTeamaHAk/Osd59xDUJPpC8tkjmzWFmuwI1d/9rd38t8G4z\n",
"W1BI9CIiUjl9P8Eysy2BI2mq6ubup3ZZ1YGj3f2OZDtG3C37DPFdxIuAqvYOpTZYIiLF+ldgfXe/\n",
"w8zOAN4OvCXHer2MvzgP2AzAzNYnPu+pxoaIiLQ0SBusi4BfA7f2uJ6x7pOzXuvBT7LsE6yeummf\n",
"xrYa05gzTGfe05gzTG/eRXL36zK/fxn4cs711pjZicTjLzqZ8RfN7CjgIXc/P1n2N0kNiyuIz2Gf\n",
"Ts9ZIiIizQYpYD3o7n/bx3orgK+b2X3Ae+i9HvwkWwD8mbgNlp5giYiUyN2XAEtazD+nxbyTUO+B\n",
"IiKSwyBtsK4ysz26L7Yudz/e3V9IXK0jJO7aPVc9+Kxs2wUzWzwh0/OBlbd/93N7Pv7Iym17XP9d\n",
"YxD/SKfTeeMSzwin3zVm8Yxieuo+39nfxyWeUUwXzeKxqd5kZpsMax8iIiK9MHfvb0WzJcQ9CF6d\n",
"ne/uL8+5/p7EY5AcS9y17iHEBb6L3X3/TusuXbrUZ2ZmrJ+4yxSE0X8CPyPu4vfKRr22fd51zWzx\n",
"tFUnmsacYTrznsacYTrzLvr728yeAfwNcZvgnwJfcvfLi9p+Pyb1HCUiMu2K+v4epIrgR/tZyczO\n",
"ArYjrir4D53qwVeQ2mD1YBpzhunMexpzhunNu0hJG6z3mtn7iDtLOtnMnujuu5ccmoiITKlBxsH6\n",
"UZ/rHdtiXst68BWUFrAeBBYEYTSvUa+pJyoRkQGY2ROB1wKvA+4H/q3ciEREZJoNNA6Wme1kZodn\n",
"pjcePKRKmw+sbNRrjxN31577eA2zDcO4msacYTrznsacYXrzLpKZ/QD4CXFb3qPc/TB3P6vksERE\n",
"ZIoNMg7W64B/IC40XJDMvgA4sIC4qip9ggWz1QRXlBeOiMjEO8Xdqzp2ooiITKBBnmC9DTgI+Etm\n",
"nhr1dtZcwMrdVfs0ttWYxpxhOvOexpxhevMukgpXIiIybgYpYK1290fTiaSL3PmDh1Rp2QKWxsIS\n",
"EREREamYQQpYPzOzjwMLzexlwPnA1/KubGYbmNltZva2ZPoQM7vczC4zs9oAcY2z+cRtr6DHJ1jT\n",
"2FZjGnOG6cx7GnOG6c17XJjZTN7zjpltb2ZRsuwpo4pRREQmzyDdtL8feBNwG3HPTae5+5k9rH8c\n",
"cBWAmRlx9+wzxNUMLwKqWO2jVRssEREZseS882Hyn3c+CXzQ3a8YQXgiIjLBBummfQ3wxeTVEzOb\n",
"D7wI+AawCbAbcJO7r0reX2Zmi9x9Wb/xjSm1werBNOYM05n3NOYM05v3mMh93jGzecAiFa5ERCSP\n",
"QZ5gDeJ44HPAtsn0VsByMzuV+E7i8mReZQpYQRgZ61YRVBssEZHy9HLe2RrYyMzOJa558Dl3P3dk\n",
"kYqIyETpuw2Wma0wsweaXznWWwgc4O4XprOA+4jHMPlA8toimddpO4uzv4/79DUfOfoQYHWjXnvc\n",
"zBY/fM8tm5MUsHJu713jlM8optN54xLPCKffNWbxjGJ66j7f2d/HJZ5RTI+RXs479xEPYPxq4HDg\n",
"AxbXxBAREZnD3L2YDZm9ANjX3T/bZbmXAO8G/gjsAjwBeCPweeAQ4kLfxe6+f7ttLF261GdmZiaq\n",
"S/ggjLYAbmnUa1sk0/8CbNSo1z6YZ30zWzxt1YmmMWeYzrynMWeYzrzH5fvb4mp/l5HzvGNmXwfq\n",
"7n6XmV0OvMjdV7ZadlxyFBGR3hT1/V1YFUF3/6mZHZtjufOJexzEzF4PbOLu15rZh4FLACfu8KJq\n",
"su2vIK4iuHXelaftIgymM2eYzrynMWeY3ryLZmY7AU9z9wuS6Y3d/aFO67j7GjM7kRbnHTM7Cngo\n",
"OV+l3g+cbnEtjHPaFa5ERET6LmCZ2bObZm0D7NfLNtz9jMzvFwMX9xvPBMi2v4IeO7kQEZH/v737\n",
"D5KkrO84/v7s/eD2PLhDQgAVjXD4ozQpotkY9Agwh6RUErEMZVUkWmpiwB9I1EGwkhIMWsgWJFYs\n",
"gapECyxRQyxiNCcCu8cvxTga/BUC5oKKysmPA47j4Ni722/+6Gdgbndmd3qm58dOf15VUzv9THfP\n",
"8326t7uf6aefZz5JpwHvJjvGfj0lfx34w8WWjYjrgeubpF/dJO0e4LVdZdbMzEqhmztYc8cBeQg4\n",
"t4v1jbq5d7ByddNexqZEZYwZyhl3GWOG8sZdsHcBx5F1s17n5nlmZjYw3XTTfkKRGSmBZhUs38Ey\n",
"M+vOnoiYkbI6laQ1ZHezzMzMBqLjXgQtt9Xs20QwVzftZfyVu4wxQznjLmPMUN64C/ZtSRcCB0j6\n",
"Y7JnfD8/4DyZmVmJdfMM1mayB4ObiohKp+seUeP4DpaZWdHOAf4S+BlwGnBpRHxhoDkyM7NS6+YZ\n",
"rO+QjQtS75ji1PR33sPBjSRdALwS2Au8MyJ+KmkjcB5Zhe28iJjuIl/Dys9g5VTGmKGccZcxZihv\n",
"3EWKiFng8vQyMzMbuG4qWL8bESc1TH9P0o0Rcc5CC0XE3wBIehXwIUlnAB8FNpI9mPwNoCwVLN/B\n",
"MjMzMzMbId08g/UcSU+N4yRpLXBQjuVfAdwJHAXcFRG70rgiWySt7yJfw6ppN+0Tk9Nt9XZVxl+5\n",
"yxgzlDPuMsYM5Y271yS9c9B5MDOz8uqmgvUJ4PuSrpR0BXA7cFE7C0q6CTgduIqsUrZd0iWS/h7Y\n",
"Tr6K2lKxzx2sWrUyA+wBVg0sR2Zmo+m0QWfAzMzKq+MKVkRcAfwe8OX0+oOI+Fybyx4HvBm4EngQ\n",
"WAd8OL0OBLYttg5Jxze+H/bpnffc+RJSBavh8x1kPV+1s76zhimefkzX04YlP32cPmvI8tOP6dLt\n",
"343vhyU//ZguiqQdkh5t8toBHFP095mZmbVLES07AuztF0vPJXso+XXAzcCJZBW+6yJiw0LLTk1N\n",
"xcaNG5fUQJITk9MfB3bWqpWPNaTdDZxUq1a2LLa8VL6H4csYM5Qz7jLGDOWMu6jjt6RbIuLYLteR\n",
"q4MlSSuBnwAXRcSnW823FM9RZmZW3PG7m04ukHQacFREfESSgGMi4luLLPMl4DfInkd6T0TMSjof\n",
"uIHsJHd+N3kaYuPAA3PS2h4Lq2wXYVDOmKGccZcxZihv3AX5t24WTuesvB0snQ58r5vvNTOz0dfN\n",
"OFgXAyuBCeAjERGSLgIWvPsUEW9qknY9cH2neVki5vYiCDm7ajczs0xEXNzlKp7qYAlA0hZJ6yOi\n",
"aYsCSePAScC/AGu6/G4zMxth3XRy8fsR8V727RlvMO0Nl4ZWFay27mD14hmGYVfGmKGccZcxZihv\n",
"3L0gaYWkPOe0vB0snQl8iuxul5mZWUvdVLDGJC0nVaokHUmXTQ5H3Gr2rYyCx8IyM+uKpKMl1YCt\n",
"wK8l3SDpqDYW3UabHSxJOgA4NiKurSd1n3MzMxtV3VSIPk323NThqbngqcBfFJKr0TTO/DtY24G1\n",
"7Sxcxmc1yhgzlDPuMsYM5Y27YJ8Bzo2IbwBIegNZD7WL9SS4hayZIGQVppbNA8mavu8n6SrgCGCZ\n",
"pM0RcUfXuTczs5HTcQUrIj4v6XayB4T3AMdFxE8Ly9noadZE8D7g0AHkxcxsVOyqV64AIuIaSe9f\n",
"bKGFOliSdCqwMyI2pXk3AZvSZ28B1rhyZWZmrXTVpC+dYHySaU+zCta9wNHtLFzG7pzLGDOUM+4y\n",
"xgzljbtg10o6NiJuAZD0WtrsNKlVB0sRcfUCy1zZaUbNzKwcOn4GS9LhXSx7maTNkm6U9PyUtlHS\n",
"LZJullTpdN1DbJz5z2BtBZ41gLyYmY2Ks4HphkGGvwpU6wMRDzhvZmZWQt3cwfoqbd59mSsiTgeQ\n",
"dALZifDd5B+PZKlpdgdrK3BYOwuX8VfuMsYM5Yy7jDFDeeMuUkS4y3QzMxsq3fQiuKuA798BzNAw\n",
"HklEPAFskbS+gPUPk64qWGZmZmZmNvy6qWD9k6SLJT2z8ZVzHe8ALiP/eCT7jB8j6fhhn469e/cn\n",
"NRFs+PzXwMHL9huvtLG+s4Ypnn5M19OGJT99nD5ryPLTj+nS7d+N74clP/2YLpqkSyXdp6xJYP3l\n",
"poFmZjYwiuhsbGBJzXoMjIg4os3lTybrFvcfJL0AOBc4g6yJ4KXABa26zJ2amoqNGzcuqXFIJian\n",
"Z4A1tWplZk76A8BLa9XKfQstL5XvYfgyxgzljLuMMUM54y76+C3pRuCUiHikqHV2aymeo8zMrLjj\n",
"dzfdtD+/02UlvRw4PiI+mJLyjEey5ExMTq8gu1u4u8nH9Y4uFqxgle0iDMoZM5Qz7jLGDOWNu2CX\n",
"A9dI+lFjYkScOaD8mJlZyXXVTXsXrgZ+IWkz8MOIeJ+kj9JkPJIRMQ48XqtWmt0urD+HdXt/s2Rm\n",
"NhL+DvhHYGjuYJmZWbnlfgZL0pca3r+3ky+NiCMi4riIOCEi3pfSrouIDRFxbBqbZJQ066K9rq2O\n",
"Lnr5DMOwKmPMUM64yxgzlDfugn0ReBT40ZyXmZnZQHRyB6tx3KY3kP1yaAtr1oNg3b24J0Ezs069\n",
"Kr0aBdDWeIqSNgLnpWXOi4iWQ4RIugx4IVlT9rdFRLNnkc3MrOQ6qWB11itGuS1UwdoKvGixFZTx\n",
"WY0yxgzljLuMMUN54y5SRJzQ6bKSRI4xGOeM4Xg2WcdMZmZm++ikgvUcSe8nOxk9N71/SkRcUkjO\n",
"RstqFm4i2NYvrWZmNp+kw4BDyc5LAg6NiP9oY9GnxmBM69kiqZ1OlnYAT3aTZzMzG12djIN1BbA/\n",
"sAb4XHrf+LL5xln4DpafwWqijDFDOeMuY8xQ3riLJOljwHeBfyUb4uPrwNvaXDz3GIxJfQxHMzOz\n",
"eXLfwYqIrnv4k7QBuAS4MSLOTmltt4NfgvwMlplZb7wRWA+8Hfg+WYcXf93mstuAdew7BuO2hRZI\n",
"YzjeFRF3dpphMzMbbYPqpn0/4OPAKyF/O/glaLFnsA6dmJxWi27cgXI+q1HGmKGccZcxZihv3AW7\n",
"JyKekPQz4Lcj4jJJiz7XmuQag7HJGI6FmZic3g+YWeg8kHN9Lc8p6bsOqlUr985JX16rVvZMTE4v\n",
"I7s+EDADqFat7J2YnN6frGlk1KqV3ROT02PpfaT348CutNz+ZE3jVwMPAyvJmnE+CKwCdqZ1rwCe\n",
"QdYq5hfAWrJhTWZSnlbVqpVdc/JZb30TabkxgFq1sn1icnocWJbS9tSqlVbnXiYmp1emfO1M869O\n",
"8e1J5fNAWl+k+WbTvKvT+2XAXmCsVq3sbPU9Tb63PjZmvXxXpvWoVq3MTExOj9eqlSfSdphN5bsm\n",
"lcvsxOS0UlnX8/QMsruvs7VqZdfE5PThwC/JtsGONO842Y8Ps6nMHkvfX//utWTb4v4U/xjZPnBA\n",
"WsdMindZWmY12fidK9J8M2Q/Fj9Edr1TH6i1cR9cl9Lr+fgt4AHgpWQ/jiyrVSuPTUxOr0j7l1K+\n",
"9jZ8/7q0DQJ4rKE8xhpi279Wrdybym8NsDPt1yvSOup5W9FQDpHWtSpNrwMeTPl4Htk+/Diwt/5/\n",
"lf4fdtbXl/5HVjZs17GUdmAq0yfr+3XDfrA3fe/yWrWyp+GzVWlbrgZ21eNM+8KyVFaN6xpLZbBf\n",
"iqXpfp/iO4Ts2nP3Ysebicnp5WT7zsFk2/kBsv/rJxvKof6/vzvtvyvStjgQWFWrVn45MTl9SFp2\n",
"jGz/OQB4NvCTtH1m0/9afTsLeB7ZD14Pp7Q96e+hwK9SFtem7bKqVq08mvYF0ncIWF6rVh5L+ayX\n",
"30pgZT09fTYO/CYwc+HLFiqR9iliMH1WSDoOeF1EnC3pBcA5EfH29NlngI+3OtEVNcpyv0xMTp8G\n",
"vKZWrby5xecPA0fVqpUH+5szM7P+Kvr4LemzwAfIBnK/FfgacGJEvKLN5V8NfITsgu2j9WFCJJ0K\n",
"7IyITQ3z3k1WCZgljeHYbJ1TU1Nxzn/pT8kuHl5KVulYxdMX0914kuwiyszMCnbhy+KqIs5RnTyD\n",
"1QudtoNfKhZ6BgvaeA6rjM9qlDFmKGfcZYwZyht3wc6MiIciYgdwGtkdklPaXTgirm82BmNEXN1Y\n",
"uUpp88ZwXMBKssoVZJWrelq3uqlctepsqRsP9GCddXsX+Xx3weubbZL280U+b3RfvuwMtXbL9hGy\n",
"u1YLWazci/Bok7T7U/qvyO5+NPOT9Ldx7Lwgy/N9wD1kP6q4U5vuNfv/2Tpn+v4m8/x6zvQP2ljv\n",
"wA2qieBcnbSDP77evKZ+kTKs0zvu/sHvLH/GugMb8974+e7HHnnikR/f+kdQ+dEC6zsaGIp4+jXd\n",
"qrxGfRo4WtLQ5KdP06Xbv8s6XbRUsaq/H6ZBhq/h6Yuy5WQVoxmypjFryX5IfAg4gqw5V705yWyt\n",
"Wnnq4jY151sF7KhVK0N5IVG3WFP3JexbRa5shMupI0NQHt9Lf388wDwUamJyelmtWlm0YjsEZd+p\n",
"O3qx0tTE8Koi1jXoJoInR0RV0hhwM3Ai2V216yJiQ6tll2ATwXOBtbVq5ZwWn38OuKFWrVzR35yZ\n",
"mfVXL4/fkpZHRKtfqvtmqZ2jzMwsU9TxeyBNBCV9CDgfOFnS5RExS9bJxQ1kHVx03VPhkDmK7BfK\n",
"Vtrqqt3MzDKS3jNn+p+BrZJ+LOnFA8qWmZnZYCpYEfGJiDg+Il4cEX+V0q5r1g5+qUu9vZwMfHWB\n",
"2bYCz1poPWV8VqOMMUM54y5jzFDeuAvypvobSX9O1mPaISndA96bmdnADMszWKPsGODeWrXy8wXm\n",
"2Urqst7MzNrSeP46C3hDag3x35JWDyhPZmZmrmD1weuBrywyz6KDDffqAfFhVsaYoZxxlzFmKG/c\n",
"BflfSReQdZB0R0Tc0/DZmgHlyczMbGi6aR9JqTeSdipYfgbLzCyf08m6kn4wvQcg3b2aHFSmzMzM\n",
"XMHqrReRda97+yLzbQUOaxiBep4yPqtRxpihnHGXMWYob9xFiIjHI+L8iDgvInbOSf/iIPNmZmbl\n",
"NlQVLEkbJd0i6WZJlUHnpwCvB/59sTEGatXKY2SD2q3tS67MzCzXOWcEz09mZtYjQ/MMliSRddW+\n",
"kWyw4W8A0wPNVPdOAf62zXnrz2E90uzDMj6rUcaYoZxxlzFmKG/cwyDPOWdEz09mZtYjw3QH6yjg\n",
"rojYFRFPAFskrR90pjo1MTl9GPBC4KY2F1m0q3YzMytMnnPOSJ2fzMyst4bmDhZwELBd0iVkvxBu\n",
"T2lbBpqrHCYmp18DfABYBhwMbKpVKzNtLv5dYNPE5PT9wM+Bx4AAZoHY/ei2g1YccNC2lLYUKb3G\n",
"0itoiG/OfGOAdu94+Jkr9j9wW8O8zHk/kmYe3XbQymxbd0Nz3tfLHZ4u87nl2my5xu1WX3buNut6\n",
"exQU85LTw229PL0ma9XKpi7XP6rynHOW/PnJzMz6Z5gqWNvIuts9g+wEdmlKa2pqamroLrIvfNm8\n",
"pJdMTU39Wc5ln5NeczyTEa9XNLGO8sUM5dzWZYwZ+hD38VNTU71c/1KW55yT6/wEw3mOMjOz/lDE\n",
"cJwDJI0BNwMnkv0Ce11EbBhsrszMbBTlOef4/GRmZnkMzR2siJiVdD5wA9lPuucPOEtmZjaiFjrn\n",
"SDoV2BkRmxab18zMbK6huYNlZmZmZma21A1TL4JmZmZmZmZLmitYZmZmZmZmBXEFy8zMzMzMrCBD\n",
"WcGStFHSLZJullQpat5hlzPuyyRtlnSjpOf3K49Fy7v9JK2U9DNJ7+pH/nol57Z+tqTpNO/F/cpj\n",
"0XLG/DZJ/ynpVkkn9CuPRZO0QdJ3JF3UxrwjcSzLGfNIHMfqRmUbdqrZtm9VJnnTR0mz/d7ltC9J\n",
"F6Tz3vUuo8XNvTaSdKLL6mmSPivptrRPvSWl9baMImKoXmRjjHwTWAWMAzcXMe+wvzqNBTgBuHTQ\n",
"+e9XzMCZwJeBdw06//2KG/gCcMyg893nmH9I9gPQAcBtg85/F3FvBE4BLiqyfIb51W7Mc5ZZssex\n",
"UdyGRW37VmWSN31UX2m//7TLacEyehVwmcto0XJ66trIZdW0fD4DHN4w3fMyGsY7WEcBd0XEroh4\n",
"AtgiaX0B8w67TmPZATzZ26z1TK6YJY0DJwFf6VcGe6TtuJWNv7M+Im7raw6Ll3f//iHZmEN/Alzb\n",
"jwz2QkRMAQ+3MevIHMtyxNxoKR/H6kZmG3aqybZvVSZ500fVDmAGl9NCXgHcicuopSbXRi6r+cS+\n",
"rfZ6XkZDMw5Wg4OA7ZIuISuQ7SltS5fzDrtOY3kH8Mke561X8sZ8JvAp4JD+ZK9n8sR9MLBK0jVk\n",
"d3M+FRHX9C2nxcm7rW8B3kp2QPxCX3I4WKN0LOvEUj6O1ZV9GzbTqkzGcqaPahnW93uXUxOSbgIO\n",
"AzYAR+IyamXutZH3p/l2AFdJ2ga8nz6U0TBWsLYB64AzyIK4NKV1O++wyx2LpJPJatR39j57PdF2\n",
"zJIOAI6NiE9Iemuaf6nKu48/AryR7P/1m5KuTb+gLCV5tvWRQCUi3pSmN0u6ISIe71dmB2CUjmW5\n",
"jMBxrK6023ABrcpkLGf6yGnc7yW9AJfTPBFxnKQJ4ErgvbiM5mlxbeT/uzki4kwASUcDk8DZ9LiM\n",
"hrGCtYXsVhxkQayPiFY1xDzzDrtcsUh6OXB8RHywH5nrkTwxbwD2k3QVcASwTNLmiLijD/ksWttx\n",
"R8QeSb8ADouIX0na1a9MFizPth4D1gJIWkF2UJvteQ57a7EfBEbpWFa36I8gI3IcqxvFbdip+rZv\n",
"Wiap6XPb6f3KdL802e9dTq3dBwTwf7iMmpl3bUTWAsRl1dwuYDf92J/68XBZ3hfwauBWsp3k1Q3p\n",
"pwKvbWfepfjKGffdwE3AZuCTg857P2Ju+OwtLOFOLjrY1s8FNqX53zfovPcp5g8DtwHfBt466Lx3\n",
"EfOHgBuB/wEuXyTmkTiW5Yx5JI5jo7YNi9z2ZM+GNPu/b3U8GPkybLbfu5zmldGXgCnga8CRnZTF\n",
"qJdRkzJ76trI+9O8svli+p/7GvC8fuxPSguZmZmZmZlZl4axF0EzMzMzM7MlyRUsMzMzMzOzgriC\n",
"ZWZmZmZmVhBXsMzMzMzMzAriCpaZmZmZmVlBXMEyMzMzMzMriCtYZmZmZmZmBfl/fsYiK/72+40A\n",
"AAAASUVORK5CYII=\n"
],
"text/plain": [
"<matplotlib.figure.Figure at 0x11ca92860>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"traceplot(trace, vars=['age_sd', 'σ_pH'])"
]
},
{
"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.2"
}
},
"nbformat": 4,
"nbformat_minor": 0
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment