DaskTest1000peopleonelocationMultipleSimulations.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This is DaskTest4 but assuming there's a distribution determining if an \"Ill\" person behaves an Ill or continues to follow their usual routine - so a new \"IA\", ill asymptomatic class"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "# various auxialiary routines\n",
    "#try json instead\n",
    "import csv\n",
    "import json\n",
    "def read_csv(file, json_file, format=','):\n",
    "    csv_rows = []\n",
    "    jsonfile = open(json_file, \"a\")\n",
    "    with open(file) as csvfile:\n",
    "        reader = csv.DictReader(csvfile)\n",
    "        title = reader.fieldnames\n",
    "        for row in reader:\n",
    "#            print(row)\n",
    "            json.dump(row,jsonfile)\n",
    "            jsonfile.write('\\n')\n",
    "        \n",
    "\n",
    "      \n",
    "\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "# read the data\n",
    "# json data is currently in .csv files \n",
    "# for some reason a read-in and write-out again is necessary for everythnig to go smoothly\n",
    "data_dir ='D:/Dask Tests/1000casesdata/start'\n",
    "\n",
    "def read_data(data_dir_start,first=True):\n",
    "#      import pdb;pdb.set_trace()\n",
    "      import os\n",
    "    \n",
    "      if first:\n",
    "# if this is the first time the data is being read then start by reading it from the csv files           \n",
    "         listoflocations = os.listdir(data_dir)\n",
    "         listofcsv= [file for file in listoflocations if 'csv' in file]\n",
    "         for file in listofcsv:\n",
    "             jsonfile = os.path.splitext(file)[0]\n",
    "             csvfile  = data_dir + '/' + file\n",
    "             jsonfile = data_dir + '/' + jsonfile + '.json'\n",
    "#    print(jsonfile)\n",
    "             read_csv(csvfile,jsonfile)\n",
    "\n",
    "      import dask.bag as db\n",
    "      dbl= db.read_text('D:/Dask Tests/1000casesdata/*.json')\n",
    "      dbl\n",
    "\n",
    "#convert it to json\n",
    "      import json\n",
    "      dblj = dbl.map(json.loads)\n",
    "      dblj\n",
    "      dblj.take(3)\n",
    "\n",
    "      return dblj"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "def  update_current_locations(dict, timeslice, stats):\n",
    "   \n",
    "    l=int(dict[timeslice])\n",
    "# stats is array of Status at location. At location l-1 entry [1] gives \"I\" count    \n",
    "    if stats[l-1][1]> 0 or dict[\"Status\"]==\"I\":\n",
    "       new_loc= dict[\"Iloc\"]\n",
    "#       dict.updat'e({'curloc':new_loc}) \n",
    "       dict[\"curloc\"]=new_loc\n",
    "    else:\n",
    "       new_loc =dict[timeslice] \n",
    "#       dict.update({'curloc':new_loc})\n",
    "#       import pdb;pdb.set_trace()\n",
    "       dict[\"curloc\"]=new_loc\n",
    "    return dict    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "def get_statistics(bag,locations,slice):\n",
    "# bag is dask bag with all the data in it    \n",
    "     timeslice =\"t\"+str(slice)\n",
    "#    print (timeslice)\n",
    "    \n",
    "#    import pdb;pdb.set_trace()\n",
    "     Count = 0\n",
    "     Total_I =0\n",
    "     stats=[]\n",
    "     while (Count < 2):\n",
    "        \n",
    "#        import pdb; pdb.set_trace()   \n",
    "        for l in range(locations):\n",
    "#        print(str(l+1))\n",
    "              l_list=[]\n",
    "#              import pdb;pdb.set_trace()\n",
    "# first scan based on individual's usual location or \"Ill\" location\n",
    "\n",
    "              if Count == 0:\n",
    "                  bgl=bag.filter(lambda d :d[timeslice]==str(l+1) and(d[\"Status\"] !=\"I\" or d[\"Iloc\"]==str(l+1))or (d[\"Iloc\"]==str(l+1) and d[\"Status\"]==\"I\") )\n",
    "              else:\n",
    "# second scan taking account of \"I\" status of everyone else at same location                    \n",
    "                  bgl=bag.filter(lambda d :d[\"curloc\"]==str(l+1))\n",
    "              Status_list = bgl.map(lambda d : d[\"Status\"]).frequencies().compute()\n",
    "              s_list = [s[0] for s in Status_list]\n",
    "              count_list = [s[1] for s in Status_list]\n",
    "              S_count = 0\n",
    "              I_count = 0\n",
    "              R_count = 0\n",
    "              if 'S' in s_list:\n",
    "                  S_count = (count_list[s_list.index('S')])\n",
    "              if 'I' in s_list:    \n",
    "                  I_count = (count_list[s_list.index('I')])\n",
    "              if 'IA' in s_list:\n",
    "                  I_count += (count_list[s_list.index('IA')])\n",
    "              if 'R' in s_list:  \n",
    "                  R_count = (count_list[s_list.index('R')])\n",
    "              Total_I += I_count   \n",
    "              l_list.append(S_count)      \n",
    "              l_list.append(I_count)\n",
    "              l_list.append(R_count)\n",
    "              l_list.append(sum(l_list))\n",
    "#        print(l_list)\n",
    "              stats.append(l_list)\n",
    "        if Count == 0:\n",
    "           \n",
    "            if Total_I >0:\n",
    "# if there are any ill then set up current_location field for use  in the second loop\n",
    "#               import pdb;pdb.set_trace()\n",
    "               stats0=stats\n",
    "#if anyone's Ill on the first scan decide where everyone should be allocated to \n",
    "#               import pdb;pdb.set_trace()\n",
    "               bagn =bag.map(lambda d: update_current_locations(d,timeslice,stats0))\n",
    "               bag = bagn\n",
    "               stats=[] \n",
    "#db               import pdb;pdb.set_trace()\n",
    "               Count =1\n",
    "            else:\n",
    "               Count=2\n",
    "        else:\n",
    "            Count =2\n",
    "#     import pdb;pdb.set_trace()      \n",
    "     return stats ,bag                 \n",
    "    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "def SIR_calc(ID_dict,timeslice,statistics):\n",
    "#def SIR_calc(ID_dict):\n",
    "# ID_dict is list of dictionnaries of individuals in current locations\n",
    "    beta = 1.2\n",
    "    gamma = 0.1\n",
    "    dt =0.25\n",
    "    import scipy\n",
    "    import numpy as np; \n",
    "#    import pdb; pdb.set_trace()\n",
    "    tslice =\"t\"+str(timeslice)\n",
    "    if int(ID_dict['curloc']) !=0:\n",
    "        l= int(ID_dict['curloc'])\n",
    "    else:    \n",
    "        l=int(ID_dict[tslice])\n",
    "      \n",
    "    totals = statistics[l-1]\n",
    "\n",
    "    N =totals[3]\n",
    "    S_tot=  totals[0]\n",
    "    I_tot = totals[1]\n",
    "    R_tot = totals[2] \n",
    "  \n",
    "    if I_tot==0:\n",
    "      return ID_dict\n",
    "# S-> I rate\n",
    "  \n",
    "    I_Rate = beta*I_tot/N*dt\n",
    "    \n",
    "  \n",
    "      \n",
    "   \n",
    "   \n",
    "\n",
    "# I-> R rate             \n",
    "    R_rate =gamma*dt\n",
    "# assume exponential distribution between events( if number of events in time period is Poisson)\n",
    "    I_prob = 1.0 - np.exp(-I_Rate)\n",
    "    R_prob = 1.0 - np.exp(-R_rate)\n",
    "   \n",
    "    outcomes = [0,1]     \n",
    "    if ID_dict[\"Status\"] ==\"S\":\n",
    "       SI_changes = np.random.choice(outcomes,1,p=[1.0-I_prob,I_prob])\n",
    "    \n",
    "       if SI_changes==1:\n",
    "# assume 30% asymptomatic            \n",
    "           IA =np.random.choice(outcomes,1,p=[0.7,0.3])\n",
    "           if IA==1:\n",
    "               ID_dict[\"Status\"]=\"IA\"\n",
    "           else:\n",
    "               ID_dict[\"Status\"]=\"I\"\n",
    "    elif ID_dict[\"Status\"]==\"I\" or ID_dict[\"Status\"]==\"IA\": \n",
    "# either I or IA can move to R        \n",
    "         IR_changes = np.random.choice(outcomes,1,p=[1.0-R_prob,R_prob])           \n",
    "         if IR_changes==1:\n",
    "              ID_dict[\"Status\"]=\"R\"\n",
    "    return ID_dict"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "from dask.distributed import Client, progress\n",
    "client = Client(n_workers=4, threads_per_worker=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "basic stats\n",
      "[[794, 206, 0, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 0 dbln00\n",
      "[('S', 753), ('I', 232), ('IA', 9), ('R', 6)]\n",
      "basic stats\n",
      "[[753, 241, 6, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 1 dbln01\n",
      "[('S', 701), ('I', 265), ('IA', 22), ('R', 12)]\n",
      "basic stats\n",
      "[[701, 287, 12, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 2 dbln02\n",
      "[('S', 655), ('I', 292), ('IA', 31), ('R', 22)]\n",
      "basic stats\n",
      "[[655, 323, 22, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 3 dbln03\n",
      "[('S', 617), ('I', 320), ('IA', 38), ('R', 25)]\n",
      "basic stats\n",
      "[[617, 358, 25, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 0 dbln10\n",
      "[('S', 585), ('I', 337), ('IA', 41), ('R', 37)]\n",
      "basic stats\n",
      "[[585, 378, 37, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 1 dbln11\n",
      "[('S', 555), ('I', 346), ('R', 51), ('IA', 48)]\n",
      "basic stats\n",
      "[[555, 394, 51, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 2 dbln12\n",
      "[('I', 366), ('S', 520), ('R', 59), ('IA', 55)]\n",
      "basic stats\n",
      "[[520, 421, 59, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 3 dbln13\n",
      "[('I', 376), ('S', 495), ('R', 70), ('IA', 59)]\n",
      "basic stats\n",
      "[[495, 435, 70, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 0 dbln20\n",
      "[('I', 389), ('S', 466), ('R', 78), ('IA', 67)]\n",
      "basic stats\n",
      "[[466, 456, 78, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 1 dbln21\n",
      "[('I', 398), ('S', 439), ('IA', 78), ('R', 85)]\n",
      "basic stats\n",
      "[[439, 476, 85, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 2 dbln22\n",
      "[('I', 410), ('S', 411), ('IA', 85), ('R', 94)]\n",
      "basic stats\n",
      "[[411, 495, 94, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 3 dbln23\n",
      "[('I', 420), ('S', 379), ('IA', 93), ('R', 108)]\n",
      "basic stats\n",
      "[[379, 513, 108, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 0 dbln30\n",
      "[('I', 425), ('S', 352), ('IA', 97), ('R', 126)]\n",
      "basic stats\n",
      "[[352, 522, 126, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 1 dbln31\n",
      "[('I', 420), ('S', 341), ('IA', 97), ('R', 142)]\n",
      "basic stats\n",
      "[[341, 517, 142, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 2 dbln32\n",
      "[('I', 423), ('S', 321), ('IA', 99), ('R', 157)]\n",
      "basic stats\n",
      "[[321, 522, 157, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 3 dbln33\n",
      "[('I', 419), ('S', 304), ('IA', 102), ('R', 175)]\n",
      "basic stats\n",
      "[[304, 521, 175, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 0 dbln40\n",
      "[('I', 419), ('S', 283), ('IA', 107), ('R', 191)]\n",
      "basic stats\n",
      "[[283, 526, 191, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 1 dbln41\n",
      "[('I', 419), ('S', 269), ('IA', 114), ('R', 198)]\n",
      "basic stats\n",
      "[[269, 533, 198, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 2 dbln42\n",
      "[('I', 418), ('R', 212), ('S', 255), ('IA', 115)]\n",
      "basic stats\n",
      "[[255, 533, 212, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 3 dbln43\n",
      "[('I', 414), ('R', 230), ('S', 240), ('IA', 116)]\n",
      "basic stats\n",
      "[[240, 530, 230, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 0 dbln50\n",
      "[('I', 416), ('R', 242), ('S', 226), ('IA', 116)]\n",
      "basic stats\n",
      "[[226, 532, 242, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 1 dbln51\n",
      "[('I', 416), ('R', 255), ('S', 209), ('IA', 120)]\n",
      "basic stats\n",
      "[[209, 536, 255, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 2 dbln52\n",
      "[('I', 418), ('R', 267), ('S', 195), ('IA', 120)]\n",
      "basic stats\n",
      "[[195, 538, 267, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 3 dbln53\n",
      "[('I', 415), ('R', 281), ('S', 183), ('IA', 121)]\n",
      "basic stats\n",
      "[[183, 536, 281, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 0 dbln60\n",
      "[('I', 418), ('R', 290), ('S', 171), ('IA', 121)]\n",
      "basic stats\n",
      "[[171, 539, 290, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 1 dbln61\n",
      "[('I', 418), ('R', 307), ('IA', 116), ('S', 159)]\n",
      "basic stats\n",
      "[[159, 534, 307, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 2 dbln62\n",
      "[('R', 322), ('I', 412), ('IA', 114), ('S', 152)]\n",
      "basic stats\n",
      "[[152, 526, 322, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 3 dbln63\n",
      "[('R', 331), ('I', 418), ('IA', 112), ('S', 139)]\n",
      "basic stats\n",
      "[[139, 530, 331, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 0 dbln70\n",
      "[('R', 353), ('I', 405), ('IA', 111), ('S', 131)]\n",
      "basic stats\n",
      "[[131, 516, 353, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 1 dbln71\n",
      "[('R', 359), ('I', 408), ('IA', 112), ('S', 121)]\n",
      "basic stats\n",
      "[[121, 520, 359, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 2 dbln72\n",
      "[('R', 372), ('I', 403), ('IA', 112), ('S', 113)]\n",
      "basic stats\n",
      "[[113, 515, 372, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 3 dbln73\n",
      "[('R', 384), ('I', 399), ('IA', 112), ('S', 105)]\n",
      "basic stats\n",
      "[[105, 511, 384, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 0 dbln80\n",
      "[('R', 398), ('I', 391), ('IA', 111), ('S', 100)]\n",
      "basic stats\n",
      "[[100, 502, 398, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 1 dbln81\n",
      "[('R', 409), ('I', 387), ('IA', 107), ('S', 97)]\n",
      "basic stats\n",
      "[[97, 494, 409, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 2 dbln82\n",
      "[('R', 415), ('I', 389), ('IA', 108), ('S', 88)]\n",
      "basic stats\n",
      "[[88, 497, 415, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 3 dbln83\n",
      "[('R', 430), ('I', 385), ('IA', 105), ('S', 80)]\n",
      "basic stats\n",
      "[[80, 490, 430, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 0 dbln90\n",
      "[('R', 452), ('I', 371), ('IA', 103), ('S', 74)]\n",
      "basic stats\n",
      "[[74, 474, 452, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 1 dbln91\n",
      "[('R', 462), ('I', 366), ('IA', 101), ('S', 71)]\n",
      "basic stats\n",
      "[[71, 467, 462, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 2 dbln92\n",
      "[('R', 474), ('I', 360), ('IA', 100), ('S', 66)]\n",
      "basic stats\n",
      "[[66, 460, 474, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 3 dbln93\n",
      "[('R', 480), ('I', 360), ('IA', 98), ('S', 62)]\n",
      "basic stats\n",
      "[[794, 206, 0, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 0 dbln00\n",
      "[('S', 742), ('I', 240), ('IA', 13), ('R', 5)]\n",
      "basic stats\n",
      "[[742, 253, 5, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 1 dbln01\n",
      "[('S', 692), ('I', 270), ('IA', 26), ('R', 12)]\n",
      "basic stats\n",
      "[[692, 296, 12, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 2 dbln02\n",
      "[('S', 658), ('I', 286), ('IA', 36), ('R', 20)]\n",
      "basic stats\n",
      "[[658, 322, 20, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 3 dbln03\n",
      "[('S', 620), ('I', 304), ('IA', 47), ('R', 29)]\n",
      "basic stats\n",
      "[[620, 351, 29, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 0 dbln10\n",
      "[('S', 582), ('R', 41), ('I', 322), ('IA', 55)]\n",
      "basic stats\n",
      "[[582, 377, 41, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 1 dbln11\n",
      "[('S', 547), ('R', 50), ('I', 338), ('IA', 65)]\n",
      "basic stats\n",
      "[[547, 403, 50, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 2 dbln12\n",
      "[('S', 511), ('R', 65), ('I', 355), ('IA', 69)]\n",
      "basic stats\n",
      "[[511, 424, 65, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 3 dbln13\n",
      "[('S', 486), ('R', 74), ('I', 367), ('IA', 73)]\n",
      "basic stats\n",
      "[[486, 440, 74, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 0 dbln20\n",
      "[('S', 452), ('R', 79), ('I', 388), ('IA', 81)]\n",
      "basic stats\n",
      "[[452, 469, 79, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 1 dbln21\n",
      "[('S', 431), ('R', 91), ('I', 393), ('IA', 85)]\n",
      "basic stats\n",
      "[[431, 478, 91, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 2 dbln22\n",
      "[('S', 403), ('R', 99), ('I', 411), ('IA', 87)]\n",
      "basic stats\n",
      "[[403, 498, 99, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 3 dbln23\n",
      "[('S', 379), ('R', 108), ('I', 421), ('IA', 92)]\n",
      "basic stats\n",
      "[[379, 513, 108, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 0 dbln30\n",
      "[('S', 358), ('R', 118), ('I', 426), ('IA', 98)]\n",
      "basic stats\n",
      "[[358, 524, 118, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 1 dbln31\n",
      "[('S', 335), ('R', 130), ('I', 435), ('IA', 100)]\n",
      "basic stats\n",
      "[[335, 535, 130, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 2 dbln32\n",
      "[('S', 320), ('R', 144), ('I', 429), ('IA', 107)]\n",
      "basic stats\n",
      "[[320, 536, 144, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 3 dbln33\n",
      "[('S', 293), ('R', 156), ('I', 441), ('IA', 110)]\n",
      "basic stats\n",
      "[[293, 551, 156, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 0 dbln40\n",
      "[('S', 274), ('R', 170), ('I', 443), ('IA', 113)]\n",
      "basic stats\n",
      "[[274, 556, 170, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 1 dbln41\n",
      "[('S', 262), ('R', 181), ('I', 440), ('IA', 117)]\n",
      "basic stats\n",
      "[[262, 557, 181, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 2 dbln42\n",
      "[('S', 242), ('I', 441), ('R', 194), ('IA', 123)]\n",
      "basic stats\n",
      "[[242, 564, 194, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 3 dbln43\n",
      "[('S', 225), ('I', 447), ('R', 203), ('IA', 125)]\n",
      "basic stats\n",
      "[[225, 572, 203, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 0 dbln50\n",
      "[('S', 211), ('I', 435), ('R', 227), ('IA', 127)]\n",
      "basic stats\n",
      "[[211, 562, 227, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 1 dbln51\n",
      "[('S', 199), ('I', 431), ('R', 241), ('IA', 129)]\n",
      "basic stats\n",
      "[[199, 560, 241, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 2 dbln52\n",
      "[('S', 186), ('I', 432), ('R', 251), ('IA', 131)]\n",
      "basic stats\n",
      "[[186, 563, 251, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 3 dbln53\n",
      "[('S', 174), ('I', 435), ('R', 263), ('IA', 128)]\n",
      "basic stats\n",
      "[[174, 563, 263, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 0 dbln60\n",
      "[('S', 165), ('I', 426), ('R', 283), ('IA', 126)]\n",
      "basic stats\n",
      "[[165, 552, 283, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 1 dbln61\n",
      "[('S', 153), ('IA', 127), ('I', 422), ('R', 298)]\n",
      "basic stats\n",
      "[[153, 549, 298, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 2 dbln62\n",
      "[('S', 146), ('IA', 128), ('I', 414), ('R', 312)]\n",
      "basic stats\n",
      "[[146, 542, 312, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 3 dbln63\n",
      "[('S', 135), ('IA', 129), ('I', 408), ('R', 328)]\n",
      "basic stats\n",
      "[[135, 537, 328, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 0 dbln70\n",
      "[('S', 122), ('IA', 133), ('I', 408), ('R', 337)]\n",
      "basic stats\n",
      "[[122, 541, 337, 1000]]\n",
      "SIR_calc\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "output results\n",
      "7 1 dbln71\n",
      "[('S', 115), ('IA', 132), ('I', 405), ('R', 348)]\n",
      "basic stats\n",
      "[[115, 537, 348, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 2 dbln72\n",
      "[('S', 110), ('IA', 129), ('I', 404), ('R', 357)]\n",
      "basic stats\n",
      "[[110, 533, 357, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 3 dbln73\n",
      "[('S', 103), ('IA', 124), ('I', 399), ('R', 374)]\n",
      "basic stats\n",
      "[[103, 523, 374, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 0 dbln80\n",
      "[('S', 100), ('IA', 121), ('I', 391), ('R', 388)]\n",
      "basic stats\n",
      "[[100, 512, 388, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 1 dbln81\n",
      "[('S', 93), ('IA', 122), ('R', 398), ('I', 387)]\n",
      "basic stats\n",
      "[[93, 509, 398, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 2 dbln82\n",
      "[('S', 90), ('IA', 120), ('R', 420), ('I', 370)]\n",
      "basic stats\n",
      "[[90, 490, 420, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 3 dbln83\n",
      "[('S', 81), ('IA', 119), ('R', 427), ('I', 373)]\n",
      "basic stats\n",
      "[[81, 492, 427, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 0 dbln90\n",
      "[('S', 77), ('IA', 116), ('R', 438), ('I', 369)]\n",
      "basic stats\n",
      "[[77, 485, 438, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 1 dbln91\n",
      "[('S', 72), ('IA', 115), ('R', 448), ('I', 365)]\n",
      "basic stats\n",
      "[[72, 480, 448, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 2 dbln92\n",
      "[('S', 67), ('IA', 112), ('R', 462), ('I', 359)]\n",
      "basic stats\n",
      "[[67, 471, 462, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 3 dbln93\n",
      "[('S', 64), ('IA', 113), ('R', 466), ('I', 357)]\n",
      "basic stats\n",
      "[[794, 206, 0, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 0 dbln00\n",
      "[('S', 750), ('I', 232), ('IA', 14), ('R', 4)]\n",
      "basic stats\n",
      "[[750, 246, 4, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 1 dbln01\n",
      "[('S', 706), ('I', 261), ('IA', 22), ('R', 11)]\n",
      "basic stats\n",
      "[[706, 283, 11, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 2 dbln02\n",
      "[('I', 283), ('S', 657), ('IA', 39), ('R', 21)]\n",
      "basic stats\n",
      "[[657, 322, 21, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 3 dbln03\n",
      "[('I', 309), ('S', 615), ('IA', 47), ('R', 29)]\n",
      "basic stats\n",
      "[[615, 356, 29, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 0 dbln10\n",
      "[('I', 323), ('S', 586), ('IA', 52), ('R', 39)]\n",
      "basic stats\n",
      "[[586, 375, 39, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 1 dbln11\n",
      "[('I', 330), ('S', 558), ('IA', 65), ('R', 47)]\n",
      "basic stats\n",
      "[[558, 395, 47, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 2 dbln12\n",
      "[('I', 344), ('S', 526), ('IA', 69), ('R', 61)]\n",
      "basic stats\n",
      "[[526, 413, 61, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 3 dbln13\n",
      "[('I', 354), ('S', 493), ('IA', 76), ('R', 77)]\n",
      "basic stats\n",
      "[[493, 430, 77, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 0 dbln20\n",
      "[('I', 371), ('S', 457), ('IA', 88), ('R', 84)]\n",
      "basic stats\n",
      "[[457, 459, 84, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 1 dbln21\n",
      "[('I', 380), ('S', 427), ('IA', 93), ('R', 100)]\n",
      "basic stats\n",
      "[[427, 473, 100, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 2 dbln22\n",
      "[('I', 385), ('S', 404), ('IA', 96), ('R', 115)]\n",
      "basic stats\n",
      "[[404, 481, 115, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 3 dbln23\n",
      "[('I', 401), ('S', 364), ('IA', 106), ('R', 129)]\n",
      "basic stats\n",
      "[[364, 507, 129, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 0 dbln30\n",
      "[('I', 414), ('S', 333), ('IA', 108), ('R', 145)]\n",
      "basic stats\n",
      "[[333, 522, 145, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 1 dbln31\n",
      "[('I', 405), ('S', 321), ('IA', 109), ('R', 165)]\n",
      "basic stats\n",
      "[[321, 514, 165, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 2 dbln32\n",
      "[('I', 411), ('S', 301), ('IA', 114), ('R', 174)]\n",
      "basic stats\n",
      "[[301, 525, 174, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 3 dbln33\n",
      "[('I', 416), ('S', 286), ('IA', 115), ('R', 183)]\n",
      "basic stats\n",
      "[[286, 531, 183, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 0 dbln40\n",
      "[('I', 413), ('S', 267), ('IA', 122), ('R', 198)]\n",
      "basic stats\n",
      "[[267, 535, 198, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 1 dbln41\n",
      "[('I', 413), ('S', 251), ('IA', 124), ('R', 212)]\n",
      "basic stats\n",
      "[[251, 537, 212, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 2 dbln42\n",
      "[('I', 413), ('S', 243), ('IA', 122), ('R', 222)]\n",
      "basic stats\n",
      "[[243, 535, 222, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 3 dbln43\n",
      "[('I', 410), ('S', 236), ('IA', 121), ('R', 233)]\n",
      "basic stats\n",
      "[[236, 531, 233, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 0 dbln50\n",
      "[('I', 405), ('S', 221), ('IA', 124), ('R', 250)]\n",
      "basic stats\n",
      "[[221, 529, 250, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 1 dbln51\n",
      "[('I', 399), ('S', 213), ('IA', 121), ('R', 267)]\n",
      "basic stats\n",
      "[[213, 520, 267, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 2 dbln52\n",
      "[('I', 396), ('S', 198), ('IA', 126), ('R', 280)]\n",
      "basic stats\n",
      "[[198, 522, 280, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 3 dbln53\n",
      "[('I', 393), ('S', 185), ('IA', 131), ('R', 291)]\n",
      "basic stats\n",
      "[[185, 524, 291, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 0 dbln60\n",
      "[('I', 392), ('S', 174), ('IA', 131), ('R', 303)]\n",
      "basic stats\n",
      "[[174, 523, 303, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 1 dbln61\n",
      "[('I', 396), ('S', 159), ('IA', 131), ('R', 314)]\n",
      "basic stats\n",
      "[[159, 527, 314, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 2 dbln62\n",
      "[('I', 394), ('S', 145), ('IA', 134), ('R', 327)]\n",
      "basic stats\n",
      "[[145, 528, 327, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 3 dbln63\n",
      "[('I', 396), ('S', 133), ('IA', 131), ('R', 340)]\n",
      "basic stats\n",
      "[[133, 527, 340, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 0 dbln70\n",
      "[('I', 392), ('S', 120), ('IA', 137), ('R', 351)]\n",
      "basic stats\n",
      "[[120, 529, 351, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 1 dbln71\n",
      "[('I', 391), ('S', 112), ('IA', 139), ('R', 358)]\n",
      "basic stats\n",
      "[[112, 530, 358, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 2 dbln72\n",
      "[('I', 383), ('S', 104), ('IA', 137), ('R', 376)]\n",
      "basic stats\n",
      "[[104, 520, 376, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 3 dbln73\n",
      "[('I', 376), ('S', 100), ('IA', 136), ('R', 388)]\n",
      "basic stats\n",
      "[[100, 512, 388, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 0 dbln80\n",
      "[('I', 366), ('S', 92), ('IA', 133), ('R', 409)]\n",
      "basic stats\n",
      "[[92, 499, 409, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 1 dbln81\n",
      "[('I', 360), ('S', 89), ('IA', 128), ('R', 423)]\n",
      "basic stats\n",
      "[[89, 488, 423, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 2 dbln82\n",
      "[('I', 356), ('S', 85), ('IA', 126), ('R', 433)]\n",
      "basic stats\n",
      "[[85, 482, 433, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 3 dbln83\n",
      "[('I', 351), ('S', 83), ('IA', 124), ('R', 442)]\n",
      "basic stats\n",
      "[[83, 475, 442, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 0 dbln90\n",
      "[('I', 345), ('S', 77), ('IA', 123), ('R', 455)]\n",
      "basic stats\n",
      "[[77, 468, 455, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 1 dbln91\n",
      "[('I', 346), ('S', 69), ('R', 463), ('IA', 122)]\n",
      "basic stats\n",
      "[[69, 468, 463, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 2 dbln92\n",
      "[('I', 344), ('S', 63), ('R', 472), ('IA', 121)]\n",
      "basic stats\n",
      "[[63, 465, 472, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 3 dbln93\n",
      "[('I', 337), ('S', 57), ('R', 487), ('IA', 119)]\n",
      "basic stats\n",
      "[[794, 206, 0, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 0 dbln00\n",
      "[('IA', 15), ('S', 746), ('I', 232), ('R', 7)]\n",
      "basic stats\n",
      "[[746, 247, 7, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 1 dbln01\n",
      "[('IA', 24), ('S', 704), ('I', 256), ('R', 16)]\n",
      "basic stats\n",
      "[[704, 280, 16, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 2 dbln02\n",
      "[('IA', 35), ('S', 663), ('I', 276), ('R', 26)]\n",
      "basic stats\n",
      "[[663, 311, 26, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 3 dbln03\n",
      "[('IA', 46), ('S', 629), ('I', 295), ('R', 30)]\n",
      "basic stats\n",
      "[[629, 341, 30, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 0 dbln10\n",
      "[('IA', 58), ('I', 306), ('S', 594), ('R', 42)]\n",
      "basic stats\n",
      "[[594, 364, 42, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 1 dbln11\n",
      "[('IA', 65), ('I', 317), ('S', 565), ('R', 53)]\n",
      "basic stats\n",
      "[[565, 382, 53, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 2 dbln12\n",
      "[('IA', 80), ('I', 333), ('S', 530), ('R', 57)]\n",
      "basic stats\n",
      "[[530, 413, 57, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 3 dbln13\n",
      "[('IA', 86), ('I', 347), ('S', 502), ('R', 65)]\n",
      "basic stats\n",
      "[[502, 433, 65, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 0 dbln20\n",
      "[('IA', 96), ('I', 364), ('S', 463), ('R', 77)]\n",
      "basic stats\n",
      "[[463, 460, 77, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 1 dbln21\n",
      "[('IA', 106), ('I', 377), ('S', 433), ('R', 84)]\n",
      "basic stats\n",
      "[[433, 483, 84, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 2 dbln22\n",
      "[('IA', 109), ('I', 389), ('S', 400), ('R', 102)]\n",
      "basic stats\n",
      "[[400, 498, 102, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 3 dbln23\n",
      "[('IA', 114), ('I', 394), ('S', 381), ('R', 111)]\n",
      "basic stats\n",
      "[[381, 508, 111, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 0 dbln30\n",
      "[('IA', 123), ('I', 397), ('S', 357), ('R', 123)]\n",
      "basic stats\n",
      "[[357, 520, 123, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 1 dbln31\n",
      "[('IA', 121), ('I', 393), ('S', 343), ('R', 143)]\n",
      "basic stats\n",
      "[[343, 514, 143, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 2 dbln32\n",
      "[('IA', 121), ('I', 391), ('S', 326), ('R', 162)]\n",
      "basic stats\n",
      "[[326, 512, 162, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 3 dbln33\n",
      "[('IA', 127), ('I', 395), ('R', 177), ('S', 301)]\n",
      "basic stats\n",
      "[[301, 522, 177, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 0 dbln40\n",
      "[('IA', 131), ('I', 397), ('R', 188), ('S', 284)]\n",
      "basic stats\n",
      "[[284, 528, 188, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 1 dbln41\n",
      "[('IA', 134), ('I', 395), ('R', 207), ('S', 264)]\n",
      "basic stats\n",
      "[[264, 529, 207, 1000]]\n",
      "SIR_calc\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "output results\n",
      "4 2 dbln42\n",
      "[('IA', 142), ('I', 390), ('R', 216), ('S', 252)]\n",
      "basic stats\n",
      "[[252, 532, 216, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 3 dbln43\n",
      "[('IA', 143), ('I', 390), ('R', 229), ('S', 238)]\n",
      "basic stats\n",
      "[[238, 533, 229, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 0 dbln50\n",
      "[('IA', 146), ('I', 392), ('R', 243), ('S', 219)]\n",
      "basic stats\n",
      "[[219, 538, 243, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 1 dbln51\n",
      "[('IA', 154), ('I', 400), ('R', 247), ('S', 199)]\n",
      "basic stats\n",
      "[[199, 554, 247, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 2 dbln52\n",
      "[('IA', 152), ('I', 402), ('R', 260), ('S', 186)]\n",
      "basic stats\n",
      "[[186, 554, 260, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 3 dbln53\n",
      "[('IA', 155), ('I', 399), ('R', 278), ('S', 168)]\n",
      "basic stats\n",
      "[[168, 554, 278, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 0 dbln60\n",
      "[('IA', 156), ('I', 392), ('R', 294), ('S', 158)]\n",
      "basic stats\n",
      "[[158, 548, 294, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 1 dbln61\n",
      "[('IA', 158), ('I', 386), ('R', 308), ('S', 148)]\n",
      "basic stats\n",
      "[[148, 544, 308, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 2 dbln62\n",
      "[('IA', 155), ('I', 389), ('R', 320), ('S', 136)]\n",
      "basic stats\n",
      "[[136, 544, 320, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 3 dbln63\n",
      "[('IA', 152), ('I', 388), ('R', 338), ('S', 122)]\n",
      "basic stats\n",
      "[[122, 540, 338, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 0 dbln70\n",
      "[('IA', 147), ('R', 357), ('I', 379), ('S', 117)]\n",
      "basic stats\n",
      "[[117, 526, 357, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 1 dbln71\n",
      "[('IA', 148), ('R', 366), ('I', 378), ('S', 108)]\n",
      "basic stats\n",
      "[[108, 526, 366, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 2 dbln72\n",
      "[('IA', 148), ('R', 378), ('I', 374), ('S', 100)]\n",
      "basic stats\n",
      "[[100, 522, 378, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 3 dbln73\n",
      "[('IA', 148), ('R', 392), ('I', 365), ('S', 95)]\n",
      "basic stats\n",
      "[[95, 513, 392, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 0 dbln80\n",
      "[('IA', 144), ('R', 412), ('I', 356), ('S', 88)]\n",
      "basic stats\n",
      "[[88, 500, 412, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 1 dbln81\n",
      "[('IA', 140), ('R', 423), ('I', 352), ('S', 85)]\n",
      "basic stats\n",
      "[[85, 492, 423, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 2 dbln82\n",
      "[('IA', 138), ('R', 435), ('I', 347), ('S', 80)]\n",
      "basic stats\n",
      "[[80, 485, 435, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 3 dbln83\n",
      "[('IA', 134), ('R', 447), ('I', 346), ('S', 73)]\n",
      "basic stats\n",
      "[[73, 480, 447, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 0 dbln90\n",
      "[('IA', 134), ('R', 450), ('I', 349), ('S', 67)]\n",
      "basic stats\n",
      "[[67, 483, 450, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 1 dbln91\n",
      "[('IA', 131), ('R', 464), ('I', 341), ('S', 64)]\n",
      "basic stats\n",
      "[[64, 472, 464, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 2 dbln92\n",
      "[('IA', 130), ('R', 475), ('I', 333), ('S', 62)]\n",
      "basic stats\n",
      "[[62, 463, 475, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 3 dbln93\n",
      "[('IA', 129), ('R', 492), ('I', 321), ('S', 58)]\n",
      "basic stats\n",
      "[[794, 206, 0, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 0 dbln00\n",
      "[('I', 237), ('S', 742), ('R', 5), ('IA', 16)]\n",
      "basic stats\n",
      "[[742, 253, 5, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 1 dbln01\n",
      "[('I', 261), ('S', 705), ('R', 7), ('IA', 27)]\n",
      "basic stats\n",
      "[[705, 288, 7, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 2 dbln02\n",
      "[('I', 294), ('IA', 39), ('S', 655), ('R', 12)]\n",
      "basic stats\n",
      "[[655, 333, 12, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 3 dbln03\n",
      "[('I', 322), ('IA', 50), ('S', 611), ('R', 17)]\n",
      "basic stats\n",
      "[[611, 372, 17, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 0 dbln10\n",
      "[('I', 342), ('IA', 60), ('S', 572), ('R', 26)]\n",
      "basic stats\n",
      "[[572, 402, 26, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 1 dbln11\n",
      "[('I', 357), ('IA', 72), ('S', 537), ('R', 34)]\n",
      "basic stats\n",
      "[[537, 429, 34, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 2 dbln12\n",
      "[('I', 366), ('IA', 79), ('S', 510), ('R', 45)]\n",
      "basic stats\n",
      "[[510, 445, 45, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 3 dbln13\n",
      "[('I', 371), ('IA', 81), ('S', 486), ('R', 62)]\n",
      "basic stats\n",
      "[[486, 452, 62, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 0 dbln20\n",
      "[('I', 386), ('IA', 86), ('S', 458), ('R', 70)]\n",
      "basic stats\n",
      "[[458, 472, 70, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 1 dbln21\n",
      "[('I', 395), ('IA', 91), ('S', 431), ('R', 83)]\n",
      "basic stats\n",
      "[[431, 486, 83, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 2 dbln22\n",
      "[('I', 405), ('R', 95), ('S', 400), ('IA', 100)]\n",
      "basic stats\n",
      "[[400, 505, 95, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 3 dbln23\n",
      "[('I', 412), ('R', 108), ('S', 376), ('IA', 104)]\n",
      "basic stats\n",
      "[[376, 516, 108, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 0 dbln30\n",
      "[('I', 418), ('R', 117), ('IA', 114), ('S', 351)]\n",
      "basic stats\n",
      "[[351, 532, 117, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 1 dbln31\n",
      "[('I', 426), ('R', 129), ('S', 325), ('IA', 120)]\n",
      "basic stats\n",
      "[[325, 546, 129, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 2 dbln32\n",
      "[('I', 423), ('R', 141), ('S', 312), ('IA', 124)]\n",
      "basic stats\n",
      "[[312, 547, 141, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 3 dbln33\n",
      "[('I', 439), ('R', 148), ('S', 286), ('IA', 127)]\n",
      "basic stats\n",
      "[[286, 566, 148, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 0 dbln40\n",
      "[('I', 444), ('R', 156), ('S', 271), ('IA', 129)]\n",
      "basic stats\n",
      "[[271, 573, 156, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 1 dbln41\n",
      "[('I', 445), ('R', 169), ('S', 257), ('IA', 129)]\n",
      "basic stats\n",
      "[[257, 574, 169, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 2 dbln42\n",
      "[('I', 446), ('R', 185), ('S', 239), ('IA', 130)]\n",
      "basic stats\n",
      "[[239, 576, 185, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 3 dbln43\n",
      "[('I', 447), ('R', 197), ('S', 222), ('IA', 134)]\n",
      "basic stats\n",
      "[[222, 581, 197, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 0 dbln50\n",
      "[('I', 447), ('R', 208), ('S', 215), ('IA', 130)]\n",
      "basic stats\n",
      "[[215, 577, 208, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 1 dbln51\n",
      "[('I', 451), ('R', 223), ('S', 198), ('IA', 128)]\n",
      "basic stats\n",
      "[[198, 579, 223, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 2 dbln52\n",
      "[('I', 454), ('R', 234), ('S', 186), ('IA', 126)]\n",
      "basic stats\n",
      "[[186, 580, 234, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 3 dbln53\n",
      "[('I', 457), ('R', 252), ('S', 167), ('IA', 124)]\n",
      "basic stats\n",
      "[[167, 581, 252, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 0 dbln60\n",
      "[('I', 456), ('R', 263), ('S', 154), ('IA', 127)]\n",
      "basic stats\n",
      "[[154, 583, 263, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 1 dbln61\n",
      "[('I', 448), ('R', 283), ('S', 146), ('IA', 123)]\n",
      "basic stats\n",
      "[[146, 571, 283, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 2 dbln62\n",
      "[('I', 437), ('R', 301), ('S', 141), ('IA', 121)]\n",
      "basic stats\n",
      "[[141, 558, 301, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 3 dbln63\n",
      "[('I', 433), ('R', 315), ('S', 134), ('IA', 118)]\n",
      "basic stats\n",
      "[[134, 551, 315, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 0 dbln70\n",
      "[('I', 428), ('R', 330), ('S', 123), ('IA', 119)]\n",
      "basic stats\n",
      "[[123, 547, 330, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 1 dbln71\n",
      "[('I', 421), ('R', 345), ('S', 119), ('IA', 115)]\n",
      "basic stats\n",
      "[[119, 536, 345, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 2 dbln72\n",
      "[('I', 408), ('R', 369), ('S', 112), ('IA', 111)]\n",
      "basic stats\n",
      "[[112, 519, 369, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 3 dbln73\n",
      "[('I', 401), ('R', 383), ('S', 102), ('IA', 114)]\n",
      "basic stats\n",
      "[[102, 515, 383, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 0 dbln80\n",
      "[('I', 399), ('R', 393), ('S', 95), ('IA', 113)]\n",
      "basic stats\n",
      "[[95, 512, 393, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 1 dbln81\n",
      "[('I', 401), ('R', 398), ('S', 88), ('IA', 113)]\n",
      "basic stats\n",
      "[[88, 514, 398, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 2 dbln82\n",
      "[('I', 381), ('R', 424), ('S', 82), ('IA', 113)]\n",
      "basic stats\n",
      "[[82, 494, 424, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 3 dbln83\n",
      "[('I', 383), ('R', 429), ('S', 75), ('IA', 113)]\n",
      "basic stats\n",
      "[[75, 496, 429, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 0 dbln90\n",
      "[('I', 377), ('R', 441), ('S', 72), ('IA', 110)]\n",
      "basic stats\n",
      "[[72, 487, 441, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 1 dbln91\n",
      "[('I', 374), ('R', 453), ('S', 69), ('IA', 104)]\n",
      "basic stats\n",
      "[[69, 478, 453, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 2 dbln92\n",
      "[('R', 463), ('I', 368), ('S', 68), ('IA', 101)]\n",
      "basic stats\n",
      "[[68, 469, 463, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 3 dbln93\n",
      "[('R', 469), ('I', 366), ('S', 63), ('IA', 102)]\n",
      "basic stats\n",
      "[[794, 206, 0, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 0 dbln00\n",
      "[('S', 748), ('I', 239), ('IA', 9), ('R', 4)]\n",
      "basic stats\n",
      "[[748, 248, 4, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 1 dbln01\n",
      "[('S', 702), ('I', 266), ('IA', 22), ('R', 10)]\n",
      "basic stats\n",
      "[[702, 288, 10, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 2 dbln02\n",
      "[('S', 660), ('I', 292), ('IA', 30), ('R', 18)]\n",
      "basic stats\n",
      "[[660, 322, 18, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 3 dbln03\n",
      "[('S', 611), ('I', 316), ('R', 25), ('IA', 48)]\n",
      "basic stats\n",
      "[[611, 364, 25, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 0 dbln10\n",
      "[('S', 580), ('I', 330), ('R', 35), ('IA', 55)]\n",
      "basic stats\n",
      "[[580, 385, 35, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 1 dbln11\n",
      "[('S', 554), ('I', 347), ('R', 38), ('IA', 61)]\n",
      "basic stats\n",
      "[[554, 408, 38, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 2 dbln12\n",
      "[('I', 354), ('S', 527), ('R', 49), ('IA', 70)]\n",
      "basic stats\n",
      "[[527, 424, 49, 1000]]\n",
      "SIR_calc\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "output results\n",
      "1 3 dbln13\n",
      "[('I', 366), ('IA', 81), ('S', 495), ('R', 58)]\n",
      "basic stats\n",
      "[[495, 447, 58, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 0 dbln20\n",
      "[('I', 376), ('IA', 90), ('S', 465), ('R', 69)]\n",
      "basic stats\n",
      "[[465, 466, 69, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 1 dbln21\n",
      "[('I', 377), ('IA', 94), ('S', 446), ('R', 83)]\n",
      "basic stats\n",
      "[[446, 471, 83, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 2 dbln22\n",
      "[('I', 390), ('IA', 98), ('S', 422), ('R', 90)]\n",
      "basic stats\n",
      "[[422, 488, 90, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 3 dbln23\n",
      "[('I', 396), ('IA', 106), ('S', 398), ('R', 100)]\n",
      "basic stats\n",
      "[[398, 502, 100, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 0 dbln30\n",
      "[('I', 406), ('IA', 111), ('S', 374), ('R', 109)]\n",
      "basic stats\n",
      "[[374, 517, 109, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 1 dbln31\n",
      "[('I', 410), ('IA', 117), ('R', 123), ('S', 350)]\n",
      "basic stats\n",
      "[[350, 527, 123, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 2 dbln32\n",
      "[('I', 417), ('IA', 115), ('R', 136), ('S', 332)]\n",
      "basic stats\n",
      "[[332, 532, 136, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 3 dbln33\n",
      "[('I', 423), ('IA', 120), ('R', 148), ('S', 309)]\n",
      "basic stats\n",
      "[[309, 543, 148, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 0 dbln40\n",
      "[('I', 427), ('IA', 120), ('R', 162), ('S', 291)]\n",
      "basic stats\n",
      "[[291, 547, 162, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 1 dbln41\n",
      "[('I', 426), ('IA', 122), ('R', 174), ('S', 278)]\n",
      "basic stats\n",
      "[[278, 548, 174, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 2 dbln42\n",
      "[('I', 427), ('IA', 126), ('R', 186), ('S', 261)]\n",
      "basic stats\n",
      "[[261, 553, 186, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 3 dbln43\n",
      "[('I', 426), ('IA', 138), ('R', 195), ('S', 241)]\n",
      "basic stats\n",
      "[[241, 564, 195, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 0 dbln50\n",
      "[('I', 430), ('IA', 140), ('R', 205), ('S', 225)]\n",
      "basic stats\n",
      "[[225, 570, 205, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 1 dbln51\n",
      "[('I', 431), ('IA', 140), ('R', 219), ('S', 210)]\n",
      "basic stats\n",
      "[[210, 571, 219, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 2 dbln52\n",
      "[('I', 432), ('IA', 142), ('R', 231), ('S', 195)]\n",
      "basic stats\n",
      "[[195, 574, 231, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 3 dbln53\n",
      "[('I', 426), ('IA', 143), ('R', 247), ('S', 184)]\n",
      "basic stats\n",
      "[[184, 569, 247, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 0 dbln60\n",
      "[('I', 429), ('IA', 146), ('R', 261), ('S', 164)]\n",
      "basic stats\n",
      "[[164, 575, 261, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 1 dbln61\n",
      "[('I', 423), ('IA', 144), ('R', 279), ('S', 154)]\n",
      "basic stats\n",
      "[[154, 567, 279, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 2 dbln62\n",
      "[('I', 422), ('IA', 147), ('R', 286), ('S', 145)]\n",
      "basic stats\n",
      "[[145, 569, 286, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 3 dbln63\n",
      "[('I', 419), ('IA', 147), ('R', 297), ('S', 137)]\n",
      "basic stats\n",
      "[[137, 566, 297, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 0 dbln70\n",
      "[('I', 414), ('IA', 149), ('R', 310), ('S', 127)]\n",
      "basic stats\n",
      "[[127, 563, 310, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 1 dbln71\n",
      "[('I', 412), ('IA', 147), ('R', 319), ('S', 122)]\n",
      "basic stats\n",
      "[[122, 559, 319, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 2 dbln72\n",
      "[('I', 408), ('IA', 142), ('R', 334), ('S', 116)]\n",
      "basic stats\n",
      "[[116, 550, 334, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 3 dbln73\n",
      "[('R', 342), ('IA', 143), ('I', 406), ('S', 109)]\n",
      "basic stats\n",
      "[[109, 549, 342, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 0 dbln80\n",
      "[('R', 360), ('I', 398), ('IA', 139), ('S', 103)]\n",
      "basic stats\n",
      "[[103, 537, 360, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 1 dbln81\n",
      "[('R', 376), ('I', 390), ('IA', 135), ('S', 99)]\n",
      "basic stats\n",
      "[[99, 525, 376, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 2 dbln82\n",
      "[('R', 389), ('I', 383), ('IA', 134), ('S', 94)]\n",
      "basic stats\n",
      "[[94, 517, 389, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 3 dbln83\n",
      "[('R', 402), ('I', 376), ('IA', 132), ('S', 90)]\n",
      "basic stats\n",
      "[[90, 508, 402, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 0 dbln90\n",
      "[('R', 412), ('I', 372), ('IA', 132), ('S', 84)]\n",
      "basic stats\n",
      "[[84, 504, 412, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 1 dbln91\n",
      "[('R', 424), ('I', 365), ('IA', 131), ('S', 80)]\n",
      "basic stats\n",
      "[[80, 496, 424, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 2 dbln92\n",
      "[('R', 440), ('I', 357), ('IA', 126), ('S', 77)]\n",
      "basic stats\n",
      "[[77, 483, 440, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 3 dbln93\n",
      "[('R', 450), ('I', 349), ('IA', 125), ('S', 76)]\n",
      "basic stats\n",
      "[[794, 206, 0, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 0 dbln00\n",
      "[('S', 759), ('I', 225), ('R', 9), ('IA', 7)]\n",
      "basic stats\n",
      "[[759, 232, 9, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 1 dbln01\n",
      "[('S', 722), ('I', 242), ('IA', 22), ('R', 14)]\n",
      "basic stats\n",
      "[[722, 264, 14, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 2 dbln02\n",
      "[('S', 675), ('I', 268), ('IA', 34), ('R', 23)]\n",
      "basic stats\n",
      "[[675, 302, 23, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 3 dbln03\n",
      "[('S', 631), ('I', 293), ('IA', 45), ('R', 31)]\n",
      "basic stats\n",
      "[[631, 338, 31, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 0 dbln10\n",
      "[('S', 599), ('I', 304), ('IA', 55), ('R', 42)]\n",
      "basic stats\n",
      "[[599, 359, 42, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 1 dbln11\n",
      "[('S', 565), ('I', 321), ('IA', 63), ('R', 51)]\n",
      "basic stats\n",
      "[[565, 384, 51, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 2 dbln12\n",
      "[('I', 341), ('S', 532), ('IA', 71), ('R', 56)]\n",
      "basic stats\n",
      "[[532, 412, 56, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 3 dbln13\n",
      "[('I', 356), ('S', 504), ('IA', 77), ('R', 63)]\n",
      "basic stats\n",
      "[[504, 433, 63, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 0 dbln20\n",
      "[('I', 365), ('S', 478), ('IA', 79), ('R', 78)]\n",
      "basic stats\n",
      "[[478, 444, 78, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 1 dbln21\n",
      "[('I', 373), ('R', 92), ('S', 443), ('IA', 92)]\n",
      "basic stats\n",
      "[[443, 465, 92, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 2 dbln22\n",
      "[('I', 385), ('R', 105), ('S', 415), ('IA', 95)]\n",
      "basic stats\n",
      "[[415, 480, 105, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 3 dbln23\n",
      "[('I', 391), ('R', 120), ('S', 386), ('IA', 103)]\n",
      "basic stats\n",
      "[[386, 494, 120, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 0 dbln30\n",
      "[('I', 394), ('R', 129), ('S', 365), ('IA', 112)]\n",
      "basic stats\n",
      "[[365, 506, 129, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 1 dbln31\n",
      "[('I', 396), ('R', 137), ('S', 343), ('IA', 124)]\n",
      "basic stats\n",
      "[[343, 520, 137, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 2 dbln32\n",
      "[('I', 394), ('R', 154), ('S', 323), ('IA', 129)]\n",
      "basic stats\n",
      "[[323, 523, 154, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 3 dbln33\n",
      "[('I', 392), ('R', 170), ('S', 308), ('IA', 130)]\n",
      "basic stats\n",
      "[[308, 522, 170, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 0 dbln40\n",
      "[('I', 404), ('R', 180), ('S', 287), ('IA', 129)]\n",
      "basic stats\n",
      "[[287, 533, 180, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 1 dbln41\n",
      "[('I', 410), ('R', 193), ('S', 269), ('IA', 128)]\n",
      "basic stats\n",
      "[[269, 538, 193, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 2 dbln42\n",
      "[('I', 413), ('R', 210), ('S', 249), ('IA', 128)]\n",
      "basic stats\n",
      "[[249, 541, 210, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 3 dbln43\n",
      "[('I', 412), ('R', 225), ('S', 234), ('IA', 129)]\n",
      "basic stats\n",
      "[[234, 541, 225, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 0 dbln50\n",
      "[('I', 413), ('R', 238), ('S', 219), ('IA', 130)]\n",
      "basic stats\n",
      "[[219, 543, 238, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 1 dbln51\n",
      "[('I', 419), ('R', 246), ('S', 206), ('IA', 129)]\n",
      "basic stats\n",
      "[[206, 548, 246, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 2 dbln52\n",
      "[('I', 423), ('R', 254), ('S', 191), ('IA', 132)]\n",
      "basic stats\n",
      "[[191, 555, 254, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 3 dbln53\n",
      "[('I', 419), ('R', 273), ('IA', 128), ('S', 180)]\n",
      "basic stats\n",
      "[[180, 547, 273, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 0 dbln60\n",
      "[('I', 419), ('R', 286), ('IA', 132), ('S', 163)]\n",
      "basic stats\n",
      "[[163, 551, 286, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 1 dbln61\n",
      "[('I', 415), ('R', 299), ('IA', 131), ('S', 155)]\n",
      "basic stats\n",
      "[[155, 546, 299, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 2 dbln62\n",
      "[('I', 418), ('R', 308), ('IA', 132), ('S', 142)]\n",
      "basic stats\n",
      "[[142, 550, 308, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 3 dbln63\n",
      "[('I', 417), ('R', 315), ('IA', 133), ('S', 135)]\n",
      "basic stats\n",
      "[[135, 550, 315, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 0 dbln70\n",
      "[('I', 416), ('R', 326), ('IA', 132), ('S', 126)]\n",
      "basic stats\n",
      "[[126, 548, 326, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 1 dbln71\n",
      "[('I', 411), ('R', 343), ('IA', 131), ('S', 115)]\n",
      "basic stats\n",
      "[[115, 542, 343, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 2 dbln72\n",
      "[('I', 404), ('R', 357), ('IA', 131), ('S', 108)]\n",
      "basic stats\n",
      "[[108, 535, 357, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 3 dbln73\n",
      "[('I', 403), ('R', 373), ('IA', 126), ('S', 98)]\n",
      "basic stats\n",
      "[[98, 529, 373, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 0 dbln80\n",
      "[('I', 396), ('R', 385), ('IA', 125), ('S', 94)]\n",
      "basic stats\n",
      "[[94, 521, 385, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 1 dbln81\n",
      "[('I', 388), ('R', 396), ('IA', 126), ('S', 90)]\n",
      "basic stats\n",
      "[[90, 514, 396, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 2 dbln82\n",
      "[('I', 380), ('R', 410), ('IA', 125), ('S', 85)]\n",
      "basic stats\n",
      "[[85, 505, 410, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 3 dbln83\n",
      "[('I', 373), ('R', 418), ('IA', 126), ('S', 83)]\n",
      "basic stats\n",
      "[[83, 499, 418, 1000]]\n",
      "SIR_calc\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "output results\n",
      "9 0 dbln90\n",
      "[('I', 368), ('R', 433), ('IA', 122), ('S', 77)]\n",
      "basic stats\n",
      "[[77, 490, 433, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 1 dbln91\n",
      "[('I', 370), ('R', 440), ('IA', 118), ('S', 72)]\n",
      "basic stats\n",
      "[[72, 488, 440, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 2 dbln92\n",
      "[('I', 365), ('R', 450), ('IA', 115), ('S', 70)]\n",
      "basic stats\n",
      "[[70, 480, 450, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 3 dbln93\n",
      "[('R', 463), ('I', 353), ('IA', 118), ('S', 66)]\n",
      "basic stats\n",
      "[[794, 206, 0, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 0 dbln00\n",
      "[('S', 734), ('I', 239), ('IA', 23), ('R', 4)]\n",
      "basic stats\n",
      "[[734, 262, 4, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 1 dbln01\n",
      "[('S', 685), ('IA', 38), ('I', 271), ('R', 6)]\n",
      "basic stats\n",
      "[[685, 309, 6, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 2 dbln02\n",
      "[('S', 651), ('IA', 46), ('I', 292), ('R', 11)]\n",
      "basic stats\n",
      "[[651, 338, 11, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 3 dbln03\n",
      "[('S', 610), ('IA', 56), ('I', 316), ('R', 18)]\n",
      "basic stats\n",
      "[[610, 372, 18, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 0 dbln10\n",
      "[('S', 575), ('IA', 67), ('I', 330), ('R', 28)]\n",
      "basic stats\n",
      "[[575, 397, 28, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 1 dbln11\n",
      "[('S', 544), ('IA', 74), ('I', 347), ('R', 35)]\n",
      "basic stats\n",
      "[[544, 421, 35, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 2 dbln12\n",
      "[('S', 518), ('IA', 81), ('I', 355), ('R', 46)]\n",
      "basic stats\n",
      "[[518, 436, 46, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 3 dbln13\n",
      "[('S', 484), ('IA', 88), ('I', 369), ('R', 59)]\n",
      "basic stats\n",
      "[[484, 457, 59, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 0 dbln20\n",
      "[('S', 454), ('IA', 96), ('I', 385), ('R', 65)]\n",
      "basic stats\n",
      "[[454, 481, 65, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 1 dbln21\n",
      "[('I', 394), ('IA', 99), ('S', 433), ('R', 74)]\n",
      "basic stats\n",
      "[[433, 493, 74, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 2 dbln22\n",
      "[('I', 398), ('IA', 105), ('S', 410), ('R', 87)]\n",
      "basic stats\n",
      "[[410, 503, 87, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 3 dbln23\n",
      "[('I', 410), ('IA', 105), ('S', 384), ('R', 101)]\n",
      "basic stats\n",
      "[[384, 515, 101, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 0 dbln30\n",
      "[('I', 421), ('IA', 107), ('S', 355), ('R', 117)]\n",
      "basic stats\n",
      "[[355, 528, 117, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 1 dbln31\n",
      "[('I', 424), ('IA', 108), ('S', 335), ('R', 133)]\n",
      "basic stats\n",
      "[[335, 532, 133, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 2 dbln32\n",
      "[('I', 432), ('IA', 114), ('S', 309), ('R', 145)]\n",
      "basic stats\n",
      "[[309, 546, 145, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 3 dbln33\n",
      "[('I', 434), ('IA', 122), ('S', 286), ('R', 158)]\n",
      "basic stats\n",
      "[[286, 556, 158, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 0 dbln40\n",
      "[('I', 439), ('IA', 126), ('S', 265), ('R', 170)]\n",
      "basic stats\n",
      "[[265, 565, 170, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 1 dbln41\n",
      "[('I', 439), ('IA', 128), ('S', 250), ('R', 183)]\n",
      "basic stats\n",
      "[[250, 567, 183, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 2 dbln42\n",
      "[('I', 441), ('IA', 129), ('S', 234), ('R', 196)]\n",
      "basic stats\n",
      "[[234, 570, 196, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 3 dbln43\n",
      "[('I', 438), ('IA', 131), ('S', 217), ('R', 214)]\n",
      "basic stats\n",
      "[[217, 569, 214, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 0 dbln50\n",
      "[('I', 434), ('IA', 133), ('S', 210), ('R', 223)]\n",
      "basic stats\n",
      "[[210, 567, 223, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 1 dbln51\n",
      "[('I', 435), ('IA', 133), ('S', 196), ('R', 236)]\n",
      "basic stats\n",
      "[[196, 568, 236, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 2 dbln52\n",
      "[('I', 428), ('IA', 136), ('S', 183), ('R', 253)]\n",
      "basic stats\n",
      "[[183, 564, 253, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 3 dbln53\n",
      "[('I', 424), ('IA', 134), ('S', 172), ('R', 270)]\n",
      "basic stats\n",
      "[[172, 558, 270, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 0 dbln60\n",
      "[('I', 427), ('IA', 136), ('S', 160), ('R', 277)]\n",
      "basic stats\n",
      "[[160, 563, 277, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 1 dbln61\n",
      "[('I', 425), ('IA', 135), ('R', 290), ('S', 150)]\n",
      "basic stats\n",
      "[[150, 560, 290, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 2 dbln62\n",
      "[('I', 423), ('IA', 131), ('R', 306), ('S', 140)]\n",
      "basic stats\n",
      "[[140, 554, 306, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 3 dbln63\n",
      "[('I', 423), ('IA', 127), ('R', 319), ('S', 131)]\n",
      "basic stats\n",
      "[[131, 550, 319, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 0 dbln70\n",
      "[('I', 412), ('IA', 121), ('R', 341), ('S', 126)]\n",
      "basic stats\n",
      "[[126, 533, 341, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 1 dbln71\n",
      "[('I', 402), ('IA', 121), ('R', 358), ('S', 119)]\n",
      "basic stats\n",
      "[[119, 523, 358, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 2 dbln72\n",
      "[('I', 397), ('IA', 123), ('R', 369), ('S', 111)]\n",
      "basic stats\n",
      "[[111, 520, 369, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 3 dbln73\n",
      "[('I', 390), ('IA', 123), ('R', 381), ('S', 106)]\n",
      "basic stats\n",
      "[[106, 513, 381, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 0 dbln80\n",
      "[('I', 389), ('IA', 126), ('R', 392), ('S', 93)]\n",
      "basic stats\n",
      "[[93, 515, 392, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 1 dbln81\n",
      "[('I', 387), ('IA', 122), ('R', 401), ('S', 90)]\n",
      "basic stats\n",
      "[[90, 509, 401, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 2 dbln82\n",
      "[('I', 383), ('IA', 119), ('R', 410), ('S', 88)]\n",
      "basic stats\n",
      "[[88, 502, 410, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 3 dbln83\n",
      "[('I', 380), ('IA', 117), ('R', 423), ('S', 80)]\n",
      "basic stats\n",
      "[[80, 497, 423, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 0 dbln90\n",
      "[('I', 375), ('IA', 119), ('R', 432), ('S', 74)]\n",
      "basic stats\n",
      "[[74, 494, 432, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 1 dbln91\n",
      "[('I', 364), ('IA', 118), ('R', 448), ('S', 70)]\n",
      "basic stats\n",
      "[[70, 482, 448, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 2 dbln92\n",
      "[('I', 360), ('IA', 118), ('R', 455), ('S', 67)]\n",
      "basic stats\n",
      "[[67, 478, 455, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 3 dbln93\n",
      "[('I', 359), ('IA', 121), ('R', 463), ('S', 57)]\n",
      "basic stats\n",
      "[[794, 206, 0, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 0 dbln00\n",
      "[('S', 735), ('I', 249), ('IA', 13), ('R', 3)]\n",
      "basic stats\n",
      "[[735, 262, 3, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 1 dbln01\n",
      "[('S', 693), ('I', 269), ('IA', 27), ('R', 11)]\n",
      "basic stats\n",
      "[[693, 296, 11, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 2 dbln02\n",
      "[('S', 653), ('I', 293), ('IA', 36), ('R', 18)]\n",
      "basic stats\n",
      "[[653, 329, 18, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 3 dbln03\n",
      "[('S', 618), ('I', 313), ('R', 25), ('IA', 44)]\n",
      "basic stats\n",
      "[[618, 357, 25, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 0 dbln10\n",
      "[('S', 574), ('I', 337), ('R', 36), ('IA', 53)]\n",
      "basic stats\n",
      "[[574, 390, 36, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 1 dbln11\n",
      "[('S', 544), ('I', 348), ('R', 46), ('IA', 62)]\n",
      "basic stats\n",
      "[[544, 410, 46, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 2 dbln12\n",
      "[('S', 511), ('I', 361), ('R', 60), ('IA', 68)]\n",
      "basic stats\n",
      "[[511, 429, 60, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 3 dbln13\n",
      "[('S', 475), ('I', 376), ('R', 71), ('IA', 78)]\n",
      "basic stats\n",
      "[[475, 454, 71, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 0 dbln20\n",
      "[('S', 452), ('I', 380), ('R', 86), ('IA', 82)]\n",
      "basic stats\n",
      "[[452, 462, 86, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 1 dbln21\n",
      "[('S', 420), ('I', 395), ('R', 100), ('IA', 85)]\n",
      "basic stats\n",
      "[[420, 480, 100, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 2 dbln22\n",
      "[('S', 393), ('I', 405), ('R', 110), ('IA', 92)]\n",
      "basic stats\n",
      "[[393, 497, 110, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 3 dbln23\n",
      "[('S', 365), ('I', 418), ('R', 122), ('IA', 95)]\n",
      "basic stats\n",
      "[[365, 513, 122, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 0 dbln30\n",
      "[('S', 345), ('I', 422), ('R', 134), ('IA', 99)]\n",
      "basic stats\n",
      "[[345, 521, 134, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 1 dbln31\n",
      "[('S', 329), ('I', 423), ('IA', 102), ('R', 146)]\n",
      "basic stats\n",
      "[[329, 525, 146, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 2 dbln32\n",
      "[('S', 297), ('I', 430), ('IA', 112), ('R', 161)]\n",
      "basic stats\n",
      "[[297, 542, 161, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 3 dbln33\n",
      "[('S', 285), ('I', 433), ('IA', 109), ('R', 173)]\n",
      "basic stats\n",
      "[[285, 542, 173, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 0 dbln40\n",
      "[('S', 268), ('I', 430), ('IA', 110), ('R', 192)]\n",
      "basic stats\n",
      "[[268, 540, 192, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 1 dbln41\n",
      "[('S', 252), ('I', 424), ('IA', 113), ('R', 211)]\n",
      "basic stats\n",
      "[[252, 537, 211, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 2 dbln42\n",
      "[('S', 233), ('I', 422), ('IA', 118), ('R', 227)]\n",
      "basic stats\n",
      "[[233, 540, 227, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 3 dbln43\n",
      "[('S', 221), ('I', 425), ('IA', 120), ('R', 234)]\n",
      "basic stats\n",
      "[[221, 545, 234, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 0 dbln50\n",
      "[('S', 208), ('I', 422), ('IA', 119), ('R', 251)]\n",
      "basic stats\n",
      "[[208, 541, 251, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 1 dbln51\n",
      "[('S', 194), ('I', 424), ('IA', 122), ('R', 260)]\n",
      "basic stats\n",
      "[[194, 546, 260, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 2 dbln52\n",
      "[('S', 180), ('I', 426), ('IA', 126), ('R', 268)]\n",
      "basic stats\n",
      "[[180, 552, 268, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 3 dbln53\n",
      "[('S', 166), ('I', 424), ('IA', 124), ('R', 286)]\n",
      "basic stats\n",
      "[[166, 548, 286, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 0 dbln60\n",
      "[('S', 159), ('I', 418), ('IA', 124), ('R', 299)]\n",
      "basic stats\n",
      "[[159, 542, 299, 1000]]\n",
      "SIR_calc\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "output results\n",
      "6 1 dbln61\n",
      "[('S', 150), ('I', 407), ('IA', 125), ('R', 318)]\n",
      "basic stats\n",
      "[[150, 532, 318, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 2 dbln62\n",
      "[('S', 144), ('I', 408), ('IA', 125), ('R', 323)]\n",
      "basic stats\n",
      "[[144, 533, 323, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 3 dbln63\n",
      "[('S', 132), ('I', 406), ('IA', 127), ('R', 335)]\n",
      "basic stats\n",
      "[[132, 533, 335, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 0 dbln70\n",
      "[('S', 127), ('I', 399), ('IA', 125), ('R', 349)]\n",
      "basic stats\n",
      "[[127, 524, 349, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 1 dbln71\n",
      "[('S', 121), ('I', 396), ('IA', 122), ('R', 361)]\n",
      "basic stats\n",
      "[[121, 518, 361, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 2 dbln72\n",
      "[('S', 116), ('I', 384), ('IA', 121), ('R', 379)]\n",
      "basic stats\n",
      "[[116, 505, 379, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 3 dbln73\n",
      "[('S', 108), ('I', 381), ('IA', 121), ('R', 390)]\n",
      "basic stats\n",
      "[[108, 502, 390, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 0 dbln80\n",
      "[('S', 103), ('I', 373), ('IA', 118), ('R', 406)]\n",
      "basic stats\n",
      "[[103, 491, 406, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 1 dbln81\n",
      "[('S', 96), ('I', 367), ('IA', 115), ('R', 422)]\n",
      "basic stats\n",
      "[[96, 482, 422, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 2 dbln82\n",
      "[('S', 89), ('I', 361), ('IA', 114), ('R', 436)]\n",
      "basic stats\n",
      "[[89, 475, 436, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 3 dbln83\n",
      "[('S', 86), ('I', 359), ('IA', 110), ('R', 445)]\n",
      "basic stats\n",
      "[[86, 469, 445, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 0 dbln90\n",
      "[('S', 80), ('I', 353), ('IA', 113), ('R', 454)]\n",
      "basic stats\n",
      "[[80, 466, 454, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 1 dbln91\n",
      "[('S', 71), ('I', 348), ('IA', 112), ('R', 469)]\n",
      "basic stats\n",
      "[[71, 460, 469, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 2 dbln92\n",
      "[('S', 64), ('I', 343), ('IA', 108), ('R', 485)]\n",
      "basic stats\n",
      "[[64, 451, 485, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 3 dbln93\n",
      "[('S', 62), ('I', 335), ('IA', 104), ('R', 499)]\n",
      "basic stats\n",
      "[[794, 206, 0, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 0 dbln00\n",
      "[('S', 756), ('I', 231), ('IA', 10), ('R', 3)]\n",
      "basic stats\n",
      "[[756, 241, 3, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 1 dbln01\n",
      "[('S', 707), ('I', 261), ('IA', 25), ('R', 7)]\n",
      "basic stats\n",
      "[[707, 286, 7, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 2 dbln02\n",
      "[('S', 670), ('I', 276), ('IA', 38), ('R', 16)]\n",
      "basic stats\n",
      "[[670, 314, 16, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "0 3 dbln03\n",
      "[('S', 624), ('I', 302), ('IA', 52), ('R', 22)]\n",
      "basic stats\n",
      "[[624, 354, 22, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 0 dbln10\n",
      "[('S', 595), ('IA', 57), ('I', 315), ('R', 33)]\n",
      "basic stats\n",
      "[[595, 372, 33, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 1 dbln11\n",
      "[('S', 555), ('IA', 68), ('I', 333), ('R', 44)]\n",
      "basic stats\n",
      "[[555, 401, 44, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 2 dbln12\n",
      "[('I', 351), ('S', 517), ('IA', 72), ('R', 60)]\n",
      "basic stats\n",
      "[[517, 423, 60, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "1 3 dbln13\n",
      "[('I', 357), ('S', 493), ('IA', 81), ('R', 69)]\n",
      "basic stats\n",
      "[[493, 438, 69, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 0 dbln20\n",
      "[('I', 366), ('S', 470), ('IA', 84), ('R', 80)]\n",
      "basic stats\n",
      "[[470, 450, 80, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 1 dbln21\n",
      "[('I', 381), ('S', 443), ('IA', 88), ('R', 88)]\n",
      "basic stats\n",
      "[[443, 469, 88, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 2 dbln22\n",
      "[('I', 388), ('S', 416), ('IA', 94), ('R', 102)]\n",
      "basic stats\n",
      "[[416, 482, 102, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "2 3 dbln23\n",
      "[('I', 393), ('S', 397), ('IA', 96), ('R', 114)]\n",
      "basic stats\n",
      "[[397, 489, 114, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 0 dbln30\n",
      "[('I', 401), ('S', 376), ('IA', 99), ('R', 124)]\n",
      "basic stats\n",
      "[[376, 500, 124, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 1 dbln31\n",
      "[('I', 404), ('S', 360), ('IA', 101), ('R', 135)]\n",
      "basic stats\n",
      "[[360, 505, 135, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 2 dbln32\n",
      "[('I', 403), ('S', 344), ('IA', 103), ('R', 150)]\n",
      "basic stats\n",
      "[[344, 506, 150, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "3 3 dbln33\n",
      "[('I', 409), ('S', 325), ('IA', 104), ('R', 162)]\n",
      "basic stats\n",
      "[[325, 513, 162, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 0 dbln40\n",
      "[('I', 413), ('S', 302), ('IA', 109), ('R', 176)]\n",
      "basic stats\n",
      "[[302, 522, 176, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 1 dbln41\n",
      "[('I', 420), ('S', 279), ('R', 189), ('IA', 112)]\n",
      "basic stats\n",
      "[[279, 532, 189, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 2 dbln42\n",
      "[('I', 420), ('S', 264), ('R', 199), ('IA', 117)]\n",
      "basic stats\n",
      "[[264, 537, 199, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "4 3 dbln43\n",
      "[('I', 416), ('S', 246), ('R', 215), ('IA', 123)]\n",
      "basic stats\n",
      "[[246, 539, 215, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 0 dbln50\n",
      "[('I', 417), ('S', 234), ('R', 228), ('IA', 121)]\n",
      "basic stats\n",
      "[[234, 538, 228, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 1 dbln51\n",
      "[('I', 416), ('S', 223), ('R', 238), ('IA', 123)]\n",
      "basic stats\n",
      "[[223, 539, 238, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 2 dbln52\n",
      "[('I', 416), ('S', 208), ('R', 248), ('IA', 128)]\n",
      "basic stats\n",
      "[[208, 544, 248, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "5 3 dbln53\n",
      "[('I', 416), ('S', 195), ('R', 258), ('IA', 131)]\n",
      "basic stats\n",
      "[[195, 547, 258, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 0 dbln60\n",
      "[('I', 420), ('S', 184), ('R', 268), ('IA', 128)]\n",
      "basic stats\n",
      "[[184, 548, 268, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 1 dbln61\n",
      "[('I', 412), ('S', 174), ('R', 284), ('IA', 130)]\n",
      "basic stats\n",
      "[[174, 542, 284, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 2 dbln62\n",
      "[('I', 407), ('S', 168), ('R', 298), ('IA', 127)]\n",
      "basic stats\n",
      "[[168, 534, 298, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "6 3 dbln63\n",
      "[('I', 403), ('S', 161), ('R', 310), ('IA', 126)]\n",
      "basic stats\n",
      "[[161, 529, 310, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 0 dbln70\n",
      "[('I', 404), ('S', 151), ('R', 317), ('IA', 128)]\n",
      "basic stats\n",
      "[[151, 532, 317, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 1 dbln71\n",
      "[('I', 400), ('S', 146), ('R', 330), ('IA', 124)]\n",
      "basic stats\n",
      "[[146, 524, 330, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 2 dbln72\n",
      "[('I', 389), ('S', 141), ('R', 345), ('IA', 125)]\n",
      "basic stats\n",
      "[[141, 514, 345, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "7 3 dbln73\n",
      "[('I', 381), ('S', 136), ('R', 357), ('IA', 126)]\n",
      "basic stats\n",
      "[[136, 507, 357, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 0 dbln80\n",
      "[('I', 379), ('S', 127), ('R', 365), ('IA', 129)]\n",
      "basic stats\n",
      "[[127, 508, 365, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 1 dbln81\n",
      "[('I', 367), ('S', 118), ('R', 387), ('IA', 128)]\n",
      "basic stats\n",
      "[[118, 495, 387, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 2 dbln82\n",
      "[('I', 365), ('S', 106), ('R', 401), ('IA', 128)]\n",
      "basic stats\n",
      "[[106, 493, 401, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "8 3 dbln83\n",
      "[('I', 361), ('S', 102), ('R', 412), ('IA', 125)]\n",
      "basic stats\n",
      "[[102, 486, 412, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 0 dbln90\n",
      "[('I', 355), ('S', 96), ('R', 428), ('IA', 121)]\n",
      "basic stats\n",
      "[[96, 476, 428, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 1 dbln91\n",
      "[('I', 347), ('S', 91), ('R', 440), ('IA', 122)]\n",
      "basic stats\n",
      "[[91, 469, 440, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 2 dbln92\n",
      "[('I', 337), ('S', 90), ('R', 453), ('IA', 120)]\n",
      "basic stats\n",
      "[[90, 457, 453, 1000]]\n",
      "SIR_calc\n",
      "output results\n",
      "9 3 dbln93\n",
      "[('I', 335), ('S', 85), ('R', 460), ('IA', 120)]\n",
      "--Return--\n",
      "> <ipython-input-8-edb864ae8f09>(92)<module>()->None\n",
      "-> import pdb;pdb.set_trace()\n",
      "(Pdb) continue\n"
     ]
    }
   ],
   "source": [
    "# try tying everything together \n",
    "# loop for 100 days\n",
    "from dask.distributed import Client, progress\n",
    "import dask.bag as db\n",
    "client = Client(n_workers=4, threads_per_worker=1)\n",
    "results=[]\n",
    "import os\n",
    "from copy import copy\n",
    "\n",
    "result_s =[]\n",
    "result_i =[]\n",
    "result_r =[]\n",
    "iterations =0\n",
    "first =True\n",
    "#outer loop over iterations\n",
    "for iterations in range(10):\n",
    "   \n",
    "# read the data from scratch at the start of every\n",
    "# iteration - assume this is in the 'start' directory\n",
    "    dblj = read_data('D:/Dask Tests/1000casesdata/start',first)\n",
    "#    import pdb;pdb.set_trace()\n",
    "    first = False\n",
    "#  copy the working .json files to the working directory\n",
    "\n",
    "    status_s=[]\n",
    "    status_i=[]\n",
    "    status_r=[]\n",
    "\n",
    "    for day in range(10):\n",
    "         for timeslice in range(4):\n",
    "# get statistics for current time slice\n",
    "             print(\"basic stats\")\n",
    "             test_dict = dblj.take(1)\n",
    "      \n",
    "             basic_stats, dblj =  get_statistics(dblj,1,timeslice+1)\n",
    "             print(basic_stats)\n",
    "#        print(dblj.take(3))\n",
    "#        import pdb;pdb.set_trace = lambda:None\n",
    "# apply the SIR calculation\n",
    "             print(\"SIR_calc\")\n",
    "#        import pdb;pdb.trace()\n",
    "             dbln =dblj.map(lambda d, basic_stats = basic_stats: SIR_calc(d,str(timeslice+1),basic_stats))\n",
    "# get a list of files currently in the directory\n",
    "             dblnoldfiles =os.listdir('D:/Dask Tests/1000casesdata/working/')\n",
    "\n",
    "# write out the data\n",
    "             file_name = 'dbln' + str(day)+str(timeslice)\n",
    "             dbln.map(json.dumps).to_textfiles('D:/Dask Tests/1000casesdata/working/'+file_name + '*.text')\n",
    "# tidy up the old files\n",
    "#        import pdb; pdb.set_trace()\n",
    "             if len(dblnoldfiles)>0:\n",
    "                   for file in dblnoldfiles:\n",
    "                         os.remove('D:/Dask Tests/1000casesdata/working/'+file)\n",
    "#read the latest data back\n",
    "             dbl= db.read_text('D:/Dask Tests/1000casesdata/working/'+file_name+'*.text')\n",
    "#        import pdb;pdb.set_trace()\n",
    "#convert to Json\n",
    "             dblj = dbl.map(json.loads)\n",
    "\n",
    "             results.append(dblj.map(lambda d : d[\"Status\"]).frequencies().compute())\n",
    "        \n",
    "             s_list = [s[0] for s in results[-1]]\n",
    "             count_list = [s[1] for s in results[-1]]\n",
    "             S_count = 0\n",
    "             I_count = 0\n",
    "             R_count = 0\n",
    "             if 'S' in s_list:\n",
    "                  S_count = (count_list[s_list.index('S')])\n",
    "             if 'I' in s_list:    \n",
    "                  I_count = (count_list[s_list.index('I')])\n",
    "             if 'IA' in s_list:\n",
    "                  I_count += (count_list[s_list.index('IA')])\n",
    "             if 'R' in s_list:  \n",
    "                  R_count = (count_list[s_list.index('R')])\n",
    "             status_s.append(S_count)\n",
    "             status_i.append(I_count)\n",
    "             status_r.append(R_count)\n",
    "        \n",
    "             dbljs= dblj.filter(lambda d:d[\"Status\"]==\"S\")\n",
    "             dblji= dblj.filter(lambda d: d[\"Status\"]==\"I\")\n",
    "             s =dbljs.map(lambda d :d[\"Id\"]).frequencies().compute()\n",
    "             i =dblji.map(lambda d :d[\"Id\"]).frequencies().compute()\n",
    "             print(\"output results\")\n",
    "             print(day,timeslice,file_name)\n",
    "#        print(s)\n",
    "#        print(i)\n",
    "             print(results[-1])\n",
    "    result_s.append(status_s)   \n",
    "    result_i.append(status_i)\n",
    "    result_r.append(status_r)  \n",
    "\n",
    "import pdb;pdb.set_trace()       \n",
    "csv_file = 'D:/Dask Tests/1000casesdata/SIR.csv'\n",
    "with open(csv_file, 'w') as csvfile:\n",
    "         \n",
    "     writer = csv.writer(csvfile,  lineterminator='\\n')\n",
    "     for irow in range(10):\n",
    "        writer.writerow(result_s[irow]) \n",
    "     for irow in range(10):  \n",
    "        writer.writerow(result_i[irow]) \n",
    "     for irow in range(10): \n",
    "        writer.writerow(result_r[irow])\n",
    "        "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "dfl.index()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "dfl.index()"
   ]
  }
 ],
 "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.7.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}