model-extended_v4.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "from sklearn.metrics import mean_squared_error\n",
    "import pandas as pd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "N = 10000000"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Date</th>\n",
       "      <th>Uncovered</th>\n",
       "      <th>Du</th>\n",
       "      <th>Ru</th>\n",
       "      <th>Exposed</th>\n",
       "      <th>Infected</th>\n",
       "      <th>Di</th>\n",
       "      <th>Ri</th>\n",
       "      <th>Isolated</th>\n",
       "      <th>Apathetic</th>\n",
       "      <th>Da</th>\n",
       "      <th>Ra</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>16-06-21</td>\n",
       "      <td>1318</td>\n",
       "      <td>34</td>\n",
       "      <td>828</td>\n",
       "      <td>23807</td>\n",
       "      <td>3956</td>\n",
       "      <td>60</td>\n",
       "      <td>2679</td>\n",
       "      <td>1533</td>\n",
       "      <td>1102</td>\n",
       "      <td>32</td>\n",
       "      <td>686</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>17-06-21</td>\n",
       "      <td>1213</td>\n",
       "      <td>33</td>\n",
       "      <td>795</td>\n",
       "      <td>24871</td>\n",
       "      <td>3840</td>\n",
       "      <td>63</td>\n",
       "      <td>2714</td>\n",
       "      <td>1795</td>\n",
       "      <td>630</td>\n",
       "      <td>24</td>\n",
       "      <td>437</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>18-06-21</td>\n",
       "      <td>1408</td>\n",
       "      <td>40</td>\n",
       "      <td>865</td>\n",
       "      <td>20882</td>\n",
       "      <td>3883</td>\n",
       "      <td>54</td>\n",
       "      <td>1955</td>\n",
       "      <td>1951</td>\n",
       "      <td>864</td>\n",
       "      <td>25</td>\n",
       "      <td>592</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>19-06-21</td>\n",
       "      <td>798</td>\n",
       "      <td>27</td>\n",
       "      <td>577</td>\n",
       "      <td>16694</td>\n",
       "      <td>3057</td>\n",
       "      <td>67</td>\n",
       "      <td>1725</td>\n",
       "      <td>1234</td>\n",
       "      <td>361</td>\n",
       "      <td>15</td>\n",
       "      <td>245</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>20-06-21</td>\n",
       "      <td>1068</td>\n",
       "      <td>28</td>\n",
       "      <td>677</td>\n",
       "      <td>22231</td>\n",
       "      <td>3641</td>\n",
       "      <td>82</td>\n",
       "      <td>2509</td>\n",
       "      <td>1346</td>\n",
       "      <td>789</td>\n",
       "      <td>25</td>\n",
       "      <td>553</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "       Date  Uncovered  Du   Ru  Exposed  Infected  Di    Ri  Isolated  \\\n",
       "0  16-06-21       1318  34  828    23807      3956  60  2679      1533   \n",
       "1  17-06-21       1213  33  795    24871      3840  63  2714      1795   \n",
       "2  18-06-21       1408  40  865    20882      3883  54  1955      1951   \n",
       "3  19-06-21        798  27  577    16694      3057  67  1725      1234   \n",
       "4  20-06-21       1068  28  677    22231      3641  82  2509      1346   \n",
       "\n",
       "   Apathetic  Da   Ra  \n",
       "0       1102  32  686  \n",
       "1        630  24  437  \n",
       "2        864  25  592  \n",
       "3        361  15  245  \n",
       "4        789  25  553  "
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#data feeding\n",
    "df = pd.read_csv(\"Data_input_daily_n.csv\") # dataset\n",
    "df.head()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Date\n",
      "Uncovered\n",
      "Du\n",
      "Ru\n",
      "Exposed\n",
      "Infected\n",
      "Di\n",
      "Ri\n",
      "Isolated\n",
      "Apathetic\n",
      "Da\n",
      "Ra\n"
     ]
    }
   ],
   "source": [
    "for col in df.columns:\n",
    "    print(col)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(77, 12)"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# df.columns =['Date', 'Uncovered','Du', 'Ru','Exposed', 'Infected','Di', 'Ri','Isolated', 'Apathetic','Da', 'Ra']\n",
    "# for col in df.columns:\n",
    "#     print(col)\n",
    "df.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Available 77 days\n"
     ]
    }
   ],
   "source": [
    "available = len(df.index)\n",
    "print(\"Available\", available, \"days\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "start_date = df['Date'].iloc[1]   # first index take push korlam array te.\n",
    "uncovered = df['Uncovered'].values.tolist() \n",
    "Du = df['Du'].values.tolist() \n",
    "Ru = df['Ru'].values.tolist() \n",
    "exposed= df['Exposed'].values.tolist()\n",
    "infected = df['Infected'].values.tolist() \n",
    "Di = df['Di'].values.tolist()\n",
    "Ri = df['Ri'].values.tolist()\n",
    "isolated = df['Isolated'].values.tolist()\n",
    "apathetic = df['Apathetic'].values.tolist() \n",
    "Da = df['Da'].values.tolist() \n",
    "Ra = df['Ra'].values.tolist() "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "17-06-21\n"
     ]
    }
   ],
   "source": [
    "print(start_date)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "['2021-06-16', '2021-06-17', '2021-06-18', '2021-06-19', '2021-06-20', '2021-06-21', '2021-06-22', '2021-06-23', '2021-06-24', '2021-06-25', '2021-06-26', '2021-06-27', '2021-06-28', '2021-06-29', '2021-06-30', '2021-07-01', '2021-07-02', '2021-07-03', '2021-07-04', '2021-07-05', '2021-07-06', '2021-07-07', '2021-07-08', '2021-07-09', '2021-07-10', '2021-07-11', '2021-07-12', '2021-07-13', '2021-07-14', '2021-07-15', '2021-07-16', '2021-07-17', '2021-07-18', '2021-07-19', '2021-07-20', '2021-07-21', '2021-07-22', '2021-07-23', '2021-07-24', '2021-07-25', '2021-07-26', '2021-07-27', '2021-07-28', '2021-07-29', '2021-07-30', '2021-07-31', '2021-08-01', '2021-08-02', '2021-08-03', '2021-08-04', '2021-08-05', '2021-08-06', '2021-08-07', '2021-08-08', '2021-08-09', '2021-08-10', '2021-08-11', '2021-08-12', '2021-08-13', '2021-08-14', '2021-08-15', '2021-08-16', '2021-08-17', '2021-08-18', '2021-08-19', '2021-08-20', '2021-08-21', '2021-08-22', '2021-08-23', '2021-08-24', '2021-08-25', '2021-08-26', '2021-08-27', '2021-08-28', '2021-08-29', '2021-08-30', '2021-08-31', '2021-09-01', '2021-09-02', '2021-09-03', '2021-09-04', '2021-09-05', '2021-09-06', '2021-09-07', '2021-09-08', '2021-09-09', '2021-09-10', '2021-09-11', '2021-09-12', '2021-09-13', '2021-09-14', '2021-09-15']\n",
      "92\n"
     ]
    }
   ],
   "source": [
    "#Helper Function\n",
    "\n",
    "t_max = available + 15\n",
    "\n",
    "date = np.array('2021-06-16', dtype=np.datetime64) \n",
    "dates = date + np.arange(t_max)\n",
    "str_dates = []\n",
    "for i in dates:\n",
    "    str_dates.append(str(i))\n",
    "print(str_dates)\n",
    "print(len(str_dates))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "def seird_model(init_vals, params, t):\n",
    "\n",
    "    DATE_0, S_0, U_0, Du_0, Ru_0, E_0, I_0,  Di_0, Ri_0,  Is_0, A_0, Da_0, Ra_0 = init_vals  \n",
    "    DATE, S, U, Du, Ru, E, I, Di, Ri, Is, A, Da, Ra = [DATE_0] , [S_0] , [U_0] , [Du_0] , [Ru_0] , [E_0] , [I_0] ,  [Di_0] , [Ri_0] , [Is_0] , [A_0] , [Da_0] , [Ra_0]\n",
    "    \n",
    "    alpha, beta, lamda, gamma, rho, delta, miu, omega, eta, sigma, tao = params\n",
    "    \n",
    "    ##### Create next t days ######\n",
    "    date = np.array(DATE, dtype=np.datetime64)\n",
    "    dates = date + np.arange(len(t))\n",
    "    str_dates = []\n",
    "    for i in dates:\n",
    "        str_dates.append(str(i))\n",
    "        \n",
    "    ##### End creating t days #####\n",
    "    for tic in t[1:]:                # 1st din er initial value deya ache tai 1 theke shuru kore baki shob add dilam.\n",
    "\n",
    "        DATE.append(str_dates[tic])\n",
    "\n",
    "        \n",
    "        next_S = S[-1] - ( alpha * (S[-1]  * E[-1] / N ) ) - ( beta * (S[-1]  * U[-1] / N ) )  \n",
    "        next_U = U[-1] + ( beta * (S[-1]  * U[-1] / N ) ) - ( sigma * U[-1] ) - ( tao * U[-1] ) \n",
    "        next_Du = Du[-1] + ( sigma * U[-1] )\n",
    "        next_Ru = Ru[-1] + ( tao * U[-1]) \n",
    "        next_E = E[-1] + ( alpha * (S[-1]  * E[-1] / N ) ) - ( gamma * I[-1] ) - ( lamda * A[-1] ) - ( rho * Is[-1] ) \n",
    "        next_I = I[-1] + ( gamma * I[-1] ) - ( delta * I[-1] ) - ( miu * I[-1] ) # Infected\n",
    "        next_Di = Di[-1] + ( delta * I[-1] )  \n",
    "        next_Ri = Ri[-1] + ( miu * I[-1] )\n",
    "        next_Is = Is[-1] + (rho *  Is[-1] )  \n",
    "        next_A = A[-1] + ( lamda * A[-1] ) - ( omega * A[-1] ) - ( eta * A[-1] )\n",
    "        next_Da = Da[-1] + ( omega * A[-1] )\n",
    "        next_Ra = Ra[-1] + ( eta * A[-1] )\n",
    "        \n",
    "        \n",
    "        S.append(next_S)\n",
    "        U.append(next_U)\n",
    "        Du.append(next_Du)                   \n",
    "        Ru.append(next_Ru)\n",
    "        E.append(next_E) \n",
    "        I.append(next_I)\n",
    "        Di.append(next_Di)\n",
    "        Ri.append(next_Ri)\n",
    "        Is.append(next_Is)\n",
    "        A.append(next_A) \n",
    "        Da.append(next_Da)                   \n",
    "        Ra.append(next_Ra)\n",
    "                          \n",
    "\n",
    "    return np.stack([DATE, S, U, Du, Ru, E, I, Di, Ri, Is, A, Da, Ra]).T"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Spliting the data for one day rolling window approach\n",
    "\n",
    "train_min = 1  # diffrence ber korar jnno ami 1 din er diffrence k niye call kortesi\n",
    "train_max = available + 1\n",
    "\n",
    "uncovered_train = []\n",
    "Du_train = []\n",
    "Ru_train = []\n",
    "exposed_train = []\n",
    "infected_train = []\n",
    "Di_train = []\n",
    "Ri_train = []\n",
    "isolated_train = []\n",
    "apathetic_train = []\n",
    "Da_train = []\n",
    "Ra_train = []\n",
    "\n",
    "\n",
    "for i in range(train_min+1, train_max):\n",
    "    j = i - 2\n",
    "    uncovered_train.append(uncovered[j:i])\n",
    "    Du_train.append(Du[j:i])\n",
    "    Ru_train.append(Ru[j:i])\n",
    "    exposed_train.append(exposed[j:i])\n",
    "    infected_train.append(infected[j:i])\n",
    "    Di_train.append(Di[j:i])\n",
    "    Ri_train.append(Ri[j:i])\n",
    "    isolated_train.append(isolated[j:i])\n",
    "    apathetic_train.append(apathetic[j:i])\n",
    "    Da_train.append(Da[j:i])\n",
    "    Ra_train.append(Ra[j:i])\n",
    "    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0 [23807, 24871]\n",
      "1 [24871, 20882]\n",
      "2 [20882, 16694]\n",
      "3 [16694, 22231]\n",
      "4 [22231, 24057]\n",
      "5 [24057, 25028]\n",
      "6 [25028, 28256]\n",
      "7 [28256, 30391]\n",
      "8 [30391, 27653]\n",
      "9 [27653, 19262]\n",
      "10 [19262, 24400]\n",
      "11 [24400, 35059]\n",
      "12 [35059, 31982]\n",
      "13 [31982, 35105]\n",
      "14 [35105, 32055]\n",
      "15 [32055, 30012]\n",
      "16 [30012, 22687]\n",
      "17 [22687, 29879]\n",
      "18 [29879, 34002]\n",
      "19 [34002, 36631]\n",
      "20 [36631, 35639]\n",
      "21 [35639, 36850]\n",
      "22 [36850, 36586]\n",
      "23 [36586, 27884]\n",
      "24 [27884, 40015]\n",
      "25 [40015, 44067]\n",
      "26 [44067, 41755]\n",
      "27 [41755, 42490]\n",
      "28 [42490, 44941]\n",
      "29 [44941, 41945]\n",
      "30 [41945, 29214]\n",
      "31 [29214, 39806]\n",
      "32 [39806, 45012]\n",
      "33 [45012, 39510]\n",
      "34 [39510, 24979]\n",
      "35 [24979, 11486]\n",
      "36 [11486, 20493]\n",
      "37 [20493, 20827]\n",
      "38 [20827, 37587]\n",
      "39 [37587, 50952]\n",
      "40 [50952, 52478]\n",
      "41 [52478, 53877]\n",
      "42 [53877, 52282]\n",
      "43 [52282, 45044]\n",
      "44 [45044, 30980]\n",
      "45 [30980, 49529]\n",
      "46 [49529, 53462]\n",
      "47 [53462, 55284]\n",
      "48 [55284, 49514]\n",
      "49 [49514, 46995]\n",
      "50 [46995, 48015]\n",
      "51 [48015, 31714]\n",
      "52 [31714, 42003]\n",
      "53 [42003, 47207]\n",
      "54 [47207, 47424]\n",
      "55 [47424, 44430]\n",
      "56 [44430, 45087]\n",
      "57 [45087, 40641]\n",
      "58 [40641, 33330]\n",
      "59 [33330, 33001]\n",
      "60 [33001, 33015]\n",
      "61 [33015, 39278]\n",
      "62 [39278, 41014]\n",
      "63 [41014, 37226]\n",
      "64 [37226, 34892]\n",
      "65 [34892, 23882]\n",
      "66 [23882, 31689]\n",
      "67 [31689, 36789]\n",
      "68 [36789, 34708]\n",
      "69 [34708, 33640]\n",
      "70 [33640, 33111]\n",
      "71 [33111, 27578]\n",
      "72 [27578, 25129]\n",
      "73 [25129, 27921]\n",
      "74 [27921, 30855]\n",
      "75 [30855, 28097]\n"
     ]
    }
   ],
   "source": [
    "for i in range(len(exposed_train)): \n",
    "    print(i,exposed_train[i])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[ 0  1  2  3  4  5  6  7  8  9 10 11 12]\n"
     ]
    }
   ],
   "source": [
    "t2 = np.arange(0, 13, 1)   # 0 - 3 PORJNTO JABE.. R PROTTEKBAR 1 KORE KORE BARBE. 0 initial value tar mane 0 theke shuRa hoye aro 2 din porjnto simulation cholbe.\n",
    "print(t2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "2021-08-16\n",
      "2021-08-17\n",
      "2021-08-18\n",
      "2021-08-19\n",
      "2021-08-20\n",
      "2021-08-21\n",
      "2021-08-22\n",
      "2021-08-23\n",
      "2021-08-24\n",
      "2021-08-25\n",
      "2021-08-26\n",
      "2021-08-27\n",
      "2021-08-28\n",
      "2021-08-29\n",
      "2021-08-30\n"
     ]
    }
   ],
   "source": [
    "last5_vals = []\n",
    "last5_params = []\n",
    "\n",
    "\n",
    "t_incub = 15  # Assumption 30 days\n",
    "\n",
    "for sample in range(len(infected_train)):\n",
    "    \n",
    "    u_0 = uncovered_train[sample][0]\n",
    "    u_1 = uncovered_train[sample][1]\n",
    "    u_diff = u_1 - u_0\n",
    "    \n",
    "    Du_0 = Du_train[sample][0]\n",
    "    Du_1 = Du_train[sample][1]\n",
    "    Du_diff = Du_1 - Du_0\n",
    "\n",
    "    Ru_0 = Ru_train[sample][0]\n",
    "    Ru_1 = Ru_train[sample][1]\n",
    "    Ru_diff = Ru_1 - Ru_0\n",
    "    \n",
    "    e_0 = exposed_train[sample][0]\n",
    "    e_1 = exposed_train[sample][1]\n",
    "    e_diff = e_1 - e_0\n",
    "    \n",
    "    i_0 = infected_train[sample][0]\n",
    "    i_1 = infected_train[sample][1]\n",
    "    i_diff = i_1 - i_0\n",
    "    \n",
    "    Di_0 = Di_train[sample][0]\n",
    "    Di_1 = Di_train[sample][1]\n",
    "    Di_diff = Di_1 - Di_0\n",
    "    \n",
    "    Ri_0 = Ri_train[sample][0]\n",
    "    Ri_1 = Ri_train[sample][1]\n",
    "    Ri_diff = Ri_1 - Ri_0\n",
    "    \n",
    "    is_0 = isolated_train[sample][0]\n",
    "    is_1 = isolated_train[sample][1]\n",
    "    is_diff = is_1 - is_0\n",
    "    \n",
    "    a_0 = apathetic_train[sample][0]\n",
    "    a_1 = apathetic_train[sample][1]\n",
    "    a_diff = a_1 - a_0\n",
    "\n",
    "    Da_0 = Da_train[sample][0]\n",
    "    Da_1 = Da_train[sample][1]\n",
    "    Da_diff = Da_1 - Da_0\n",
    "\n",
    "    Ra_0 = Ra_train[sample][0]\n",
    "    Ra_1 = Ra_train[sample][1]\n",
    "    Ra_diff = Ra_1 - Ra_0\n",
    "    \n",
    "    s_0 = N - is_0 - Du_0 - Ru_0 - Di_0 - Ri_0 - Da_0 - Ra_0\n",
    "    \n",
    "    alpha = (( e_diff + i_diff + Di_diff + Ri_diff + a_diff + Da_diff + Ra_diff + is_diff) / (s_0 * e_0 )) * N  \n",
    "    beta =  ( ( u_diff + Du_diff + Ru_diff ) / (s_0 * u_0 ) )* N     \n",
    "    lamda =  ( ( a_diff + Da_diff + Ra_diff ) / a_0 )     \n",
    "    gamma = ( i_diff + Di_diff + Ri_diff) / i_0    \n",
    "    rho = ( is_diff / is_0 )      \n",
    "    delta = ( Di_diff / i_0 )  \n",
    "    miu =  ( Ri_diff / i_0 )     \n",
    "    omega =  ( Da_diff / a_0 )     \n",
    "    eta =  ( Ra_diff / a_0 )     \n",
    "    sigma = ( Du_diff / u_0 )     \n",
    "    tao = ( Ru_diff / u_0 )      \n",
    "    \n",
    "\n",
    "    init_vals = str_dates[sample], s_0, u_0, Du_0, Ru_0, e_0, i_0, Di_0, Ri_0, is_0, a_0, Da_0, Ra_0  \n",
    "    params = alpha, beta, lamda, gamma, rho, delta, miu, omega, eta, sigma, tao\n",
    "    pred = seird_model(init_vals, params, t2) # t2 cz etar upor porer day gula depend kortese.\n",
    "\n",
    "\n",
    "    # Store Last five cases for predecting future\n",
    "    if ((len(infected_train) - sample ) <= 15):\n",
    "        print(str_dates[sample])\n",
    "        last5_vals.append(init_vals)\n",
    "        last5_params.append(params)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[ 0  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 19 20]\n"
     ]
    }
   ],
   "source": [
    "t7 = np.arange(0, 21, 1)  #13 din er data ashbe\n",
    "print(t7)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Predictions using 2021-08-16 parameters\n",
      "\n",
      "------------------------------------------------------------------------------------------------------------------\n",
      "Date Uncovered Du Ru Exposed Infected Di Ri Isolated Apathetic Da Ra\n",
      "--------------------------------------------------------------------------------------------------------------------\n",
      "2021-08-22 2888 114 1923 105797 11214 351 36468 3504 2981 70 1925\n",
      "2021-08-23 2914 121 1963 131876 12142 389 42401 4017 2848 65 1854\n",
      "2021-08-24 2941 128 2004 165363 13147 431 48826 4605 2722 59 1786\n",
      "2021-08-25 2968 136 2045 208418 14235 477 55782 5279 2601 55 1722\n",
      "2021-08-26 2994 144 2086 263795 15413 526 63314 6052 2485 50 1660\n",
      "2021-08-27 3020 151 2128 334980 16689 579 71469 6937 2375 46 1601\n",
      "2021-08-28 3046 159 2170 426334 18071 637 80300 7953 2269 41 1545\n",
      "2021-08-29 3071 167 2212 543236 19566 699 89861 9117 2168 37 1491\n",
      "2021-08-30 3096 175 2255 692186 21186 766 100214 10452 2072 33 1439\n",
      "2021-08-31 3119 183 2298 880827 22939 839 111424 11982 1980 30 1390\n",
      "2021-09-01 3141 191 2341 1117802 24838 918 123561 13736 1892 26 1343\n",
      "2021-09-02 3161 199 2385 1412335 26894 1004 136703 15746 1808 23 1298\n",
      "2021-09-03 3178 207 2429 1773411 29120 1097 150933 18051 1727 19 1255\n",
      "2021-09-04 3193 215 2473 2208423 31531 1197 166341 20693 1651 16 1214\n",
      "2021-09-05 3203 223 2518 2721246 34140 1306 183024 23722 1577 13 1175\n",
      "---------------------------------------------------------------------------------------------------------------\n",
      "\n",
      "\n",
      "Predictions using 2021-08-17 parameters\n",
      "\n",
      "------------------------------------------------------------------------------------------------------------------\n",
      "Date Uncovered Du Ru Exposed Infected Di Ri Isolated Apathetic Da Ra\n",
      "--------------------------------------------------------------------------------------------------------------------\n",
      "2021-08-22 4763 96 2962 46859 6205 77 9067 2665 6898 152 4224\n",
      "2021-08-23 5314 102 3301 47944 5968 56 8377 2892 7796 167 4763\n",
      "2021-08-24 5927 107 3679 48825 5741 35 7713 3139 8811 184 5371\n",
      "2021-08-25 6612 114 4100 49471 5522 15 7074 3407 9958 203 6059\n",
      "2021-08-26 7374 121 4571 49844 5312 -3 6460 3697 11254 224 6836\n",
      "2021-08-27 8225 129 5095 49902 5110 -21 5869 4013 12718 248 7714\n",
      "2021-08-28 9173 138 5680 49595 4915 -39 5301 4355 14373 275 8706\n",
      "2021-08-29 10229 148 6332 48868 4728 -56 4754 4727 16244 306 9828\n",
      "2021-08-30 11407 159 7060 47659 4548 -72 4228 5130 18358 341 11096\n",
      "2021-08-31 12719 172 7871 45896 4374 -88 3723 5568 20748 380 12528\n",
      "2021-09-01 14182 185 8776 43498 4208 -103 3236 6043 23448 424 14147\n",
      "2021-09-02 15811 201 9784 40373 4048 -117 2768 6559 26500 474 15977\n",
      "2021-09-03 17626 218 10909 36416 3893 -131 2318 7118 29948 531 18045\n",
      "2021-09-04 19648 237 12162 31509 3745 -145 1885 7725 33846 595 20382\n",
      "2021-09-05 21901 259 13560 25518 3602 -158 1468 8384 38251 667 23023\n",
      "---------------------------------------------------------------------------------------------------------------\n",
      "\n",
      "\n",
      "Predictions using 2021-08-18 parameters\n",
      "\n",
      "------------------------------------------------------------------------------------------------------------------\n",
      "Date Uncovered Du Ru Exposed Infected Di Ri Isolated Apathetic Da Ra\n",
      "--------------------------------------------------------------------------------------------------------------------\n",
      "2021-08-22 1866 70 1160 27171 4881 126 5314 2302 1806 43 1197\n",
      "2021-08-23 1647 69 1021 24321 4422 118 3994 2409 1460 33 993\n",
      "2021-08-24 1453 67 899 21721 4006 110 2799 2520 1180 26 828\n",
      "2021-08-25 1282 66 791 19363 3629 103 1716 2637 954 20 695\n",
      "2021-08-26 1130 64 696 17232 3287 96 736 2759 771 15 588\n",
      "2021-08-27 997 63 612 15313 2978 90 -152 2887 624 11 501\n",
      "2021-08-28 879 62 538 13587 2698 85 -957 3021 504 7 430\n",
      "2021-08-29 775 61 473 12039 2444 80 -1686 3161 407 5 373\n",
      "2021-08-30 684 61 416 10652 2214 76 -2347 3307 329 3 327\n",
      "2021-08-31 603 60 365 9410 2005 72 -2945 3461 266 1 290\n",
      "2021-09-01 531 59 320 8298 1817 68 -3487 3621 215 0 260\n",
      "2021-09-02 469 59 281 7302 1646 65 -3979 3789 174 -1 236\n",
      "2021-09-03 413 58 246 6410 1491 62 -4423 3964 140 -2 216\n",
      "2021-09-04 364 58 215 5611 1350 59 -4827 4148 113 -2 200\n",
      "2021-09-05 321 58 188 4893 1223 57 -5192 4340 92 -3 187\n",
      "---------------------------------------------------------------------------------------------------------------\n",
      "\n",
      "\n",
      "Predictions using 2021-08-19 parameters\n",
      "\n",
      "------------------------------------------------------------------------------------------------------------------\n",
      "Date Uncovered Du Ru Exposed Infected Di Ri Isolated Apathetic Da Ra\n",
      "--------------------------------------------------------------------------------------------------------------------\n",
      "2021-08-22 1453 33 974 30132 4992 120 11308 1148 1761 65 1037\n",
      "2021-08-23 1180 23 817 27814 4556 109 11629 953 1412 60 803\n",
      "2021-08-24 957 16 691 25581 4159 100 11921 791 1132 56 615\n",
      "2021-08-25 777 9 588 23456 3796 91 12187 656 907 53 465\n",
      "2021-08-26 630 4 504 21451 3465 83 12431 544 727 51 344\n",
      "2021-08-27 511 0 436 19573 3162 75 12653 452 583 49 248\n",
      "2021-08-28 415 -3 381 17824 2886 69 12856 375 467 47 170\n",
      "2021-08-29 336 -5 336 16204 2634 62 13041 311 374 46 108\n",
      "2021-08-30 273 -8 300 14710 2404 57 13210 258 300 45 58\n",
      "2021-08-31 221 -9 271 13335 2194 52 13364 214 241 44 18\n",
      "2021-09-01 179 -11 247 12076 2003 47 13505 177 193 44 -13\n",
      "2021-09-02 145 -12 228 10924 1828 43 13633 147 154 43 -38\n",
      "2021-09-03 118 -13 212 9874 1668 39 13750 122 124 43 -59\n",
      "2021-09-04 95 -14 199 8917 1523 35 13857 101 99 42 -75\n",
      "2021-09-05 77 -14 189 8048 1390 32 13955 84 79 42 -88\n",
      "---------------------------------------------------------------------------------------------------------------\n",
      "\n",
      "\n",
      "Predictions using 2021-08-20 parameters\n",
      "\n",
      "------------------------------------------------------------------------------------------------------------------\n",
      "Date Uncovered Du Ru Exposed Infected Di Ri Isolated Apathetic Da Ra\n",
      "--------------------------------------------------------------------------------------------------------------------\n",
      "2021-08-22 1277 45 877 16332 2657 103 5729 770 1804 52 1200\n",
      "2021-08-23 971 41 704 11182 1769 92 4439 524 1464 43 1021\n",
      "2021-08-24 738 37 572 7676 1178 84 3580 356 1187 35 875\n",
      "2021-08-25 560 35 472 5289 784 79 3009 242 963 29 758\n",
      "2021-08-26 425 33 396 3663 522 76 2628 164 782 24 662\n",
      "2021-08-27 323 31 338 2551 348 74 2374 111 634 20 584\n",
      "2021-08-28 245 30 294 1789 231 73 2205 76 514 17 521\n",
      "2021-08-29 186 29 261 1264 154 72 2093 51 417 15 470\n",
      "2021-08-30 141 28 236 900 102 71 2018 35 338 12 429\n",
      "2021-08-31 107 28 216 647 68 71 1968 23 275 11 395\n",
      "2021-09-01 81 27 202 469 45 70 1935 16 223 9 368\n",
      "2021-09-02 61 27 191 343 30 70 1912 11 181 8 346\n",
      "2021-09-03 47 27 182 254 20 70 1898 7 146 7 328\n",
      "2021-09-04 35 27 176 189 13 70 1888 5 119 7 313\n",
      "2021-09-05 27 27 171 142 8 70 1881 3 96 6 301\n",
      "---------------------------------------------------------------------------------------------------------------\n",
      "\n",
      "\n",
      "Predictions using 2021-08-21 parameters\n",
      "\n",
      "------------------------------------------------------------------------------------------------------------------\n",
      "Date Uncovered Du Ru Exposed Infected Di Ri Isolated Apathetic Da Ra\n",
      "--------------------------------------------------------------------------------------------------------------------\n",
      "2021-08-22 2253 59 1411 31689 4804 139 8453 1217 3289 86 2064\n",
      "2021-08-23 3021 68 1821 41987 5782 161 9400 1306 4863 118 3014\n",
      "2021-08-24 4050 80 2372 55466 6960 189 10540 1401 7193 166 4421\n",
      "2021-08-25 5424 97 3110 72918 8378 222 11913 1504 10637 237 6500\n",
      "2021-08-26 7256 120 4099 95183 10085 262 13565 1614 15731 343 9576\n",
      "2021-08-27 9693 150 5421 123016 12139 310 15554 1732 23265 498 14124\n",
      "2021-08-28 12924 191 7188 156817 14612 368 17948 1859 34406 728 20851\n",
      "2021-08-29 17188 245 9543 196130 17589 437 20829 1995 50882 1069 30798\n",
      "2021-08-30 22789 316 12676 238765 21172 521 24298 2141 75248 1572 45509\n",
      "2021-08-31 30096 411 16829 279331 25485 622 28473 2298 111282 2316 67265\n",
      "2021-09-01 39555 537 22314 306857 30677 743 33498 2466 164571 3417 99439\n",
      "2021-09-02 51696 702 29524 301072 36926 889 39548 2646 243379 5045 147020\n",
      "2021-09-03 67141 917 38945 226703 44448 1065 46829 2840 359925 7453 217385\n",
      "2021-09-04 86649 1197 51182 24809 53503 1277 55594 3048 532282 11013 321446\n",
      "2021-09-05 111236 1558 66974 -400776 64402 1531 66145 3271 787174 16278 475339\n",
      "---------------------------------------------------------------------------------------------------------------\n",
      "\n",
      "\n",
      "Predictions using 2021-08-22 parameters\n",
      "\n",
      "------------------------------------------------------------------------------------------------------------------\n",
      "Date Uncovered Du Ru Exposed Infected Di Ri Isolated Apathetic Da Ra\n",
      "--------------------------------------------------------------------------------------------------------------------\n",
      "2021-08-22 2253 59 1411 31689 4804 139 8453 1217 3289 86 2064\n",
      "2021-08-23 1858 54 1119 36789 5717 117 8982 1620 2095 59 1307\n",
      "2021-08-24 1532 49 878 41561 6803 90 9611 2156 1334 41 824\n",
      "2021-08-25 1264 46 679 46117 8096 59 10360 2870 850 30 517\n",
      "2021-08-26 1043 43 515 50470 9635 22 11252 3821 541 23 322\n",
      "2021-08-27 861 41 380 54541 11466 -21 12313 5086 344 19 197\n",
      "2021-08-28 710 39 268 58160 13645 -74 13575 6770 219 16 118\n",
      "2021-08-29 587 37 176 61058 16239 -136 15078 9012 139 14 67\n",
      "2021-08-30 484 36 100 62841 19325 -210 16866 11997 89 13 35\n",
      "2021-08-31 400 35 37 62961 22997 -299 18994 15970 56 12 14\n",
      "2021-09-01 331 34 -14 60669 27368 -404 21527 21258 36 12 1\n",
      "2021-09-02 274 33 -57 54957 32570 -530 24540 28298 23 12 -6\n",
      "2021-09-03 226 33 -92 44476 38760 -679 28127 37668 14 11 -11\n",
      "2021-09-04 187 32 -121 27430 46126 -856 32395 50142 9 11 -15\n",
      "2021-09-05 155 32 -146 1434 54892 -1067 37474 66747 5 11 -17\n",
      "---------------------------------------------------------------------------------------------------------------\n",
      "\n",
      "\n",
      "Predictions using 2021-08-23 parameters\n",
      "\n",
      "------------------------------------------------------------------------------------------------------------------\n",
      "Date Uncovered Du Ru Exposed Infected Di Ri Isolated Apathetic Da Ra\n",
      "--------------------------------------------------------------------------------------------------------------------\n",
      "2021-09-12 9037 101 9861 -373686 1035 86 8231 39 277958 5766 181093\n",
      "2021-08-23 1858 54 1119 36789 5717 117 8982 1620 2095 59 1307\n",
      "2021-08-24 2011 55 1305 34708 5249 114 8907 1346 2675 71 1685\n",
      "2021-08-25 2176 56 1506 32376 4819 111 8838 1118 3415 86 2167\n",
      "2021-08-26 2356 57 1724 29744 4424 108 8774 929 4361 105 2783\n",
      "2021-08-27 2550 58 1960 26743 4062 106 8716 772 5568 130 3570\n",
      "2021-08-28 2760 59 2215 23278 3730 104 8663 641 7110 162 4575\n",
      "2021-08-29 2988 61 2491 19228 3424 102 8614 532 9078 203 5858\n",
      "2021-08-30 3234 62 2790 14432 3144 100 8569 442 11592 255 7496\n",
      "2021-08-31 3501 64 3114 8682 2886 98 8528 367 14801 321 9588\n",
      "2021-09-01 3790 66 3465 1711 2650 97 8490 305 18899 406 12258\n",
      "2021-09-02 4102 68 3844 -6824 2433 95 8455 253 24131 514 15668\n",
      "2021-09-03 4441 70 4255 -17363 2234 94 8423 211 30812 653 20022\n",
      "2021-09-04 4807 73 4699 -30468 2051 93 8394 175 39342 829 25582\n",
      "2021-09-05 5203 75 5181 -46859 1883 92 8367 145 50234 1054 32680\n",
      "---------------------------------------------------------------------------------------------------------------\n",
      "\n",
      "\n",
      "Predictions using 2021-08-24 parameters\n",
      "\n",
      "------------------------------------------------------------------------------------------------------------------\n",
      "Date Uncovered Du Ru Exposed Infected Di Ri Isolated Apathetic Da Ra\n",
      "--------------------------------------------------------------------------------------------------------------------\n",
      "2021-09-12 863 2 -414 13801 1831 114 -4365 5116 348 28 273\n",
      "2021-09-13 825 0 -470 12853 1732 114 -4748 5488 313 27 252\n",
      "2021-08-24 2011 55 1305 34708 5249 114 8907 1346 2675 71 1685\n",
      "2021-08-25 1924 51 1174 33640 4966 114 7808 1444 2403 66 1520\n",
      "2021-08-26 1840 47 1048 32530 4698 114 6768 1549 2158 61 1371\n",
      "2021-08-27 1760 43 928 31391 4444 114 5784 1661 1939 57 1238\n",
      "2021-08-28 1684 40 814 30231 4205 114 4853 1782 1741 53 1119\n",
      "2021-08-29 1611 36 704 29060 3978 114 3973 1912 1564 50 1011\n",
      "2021-08-30 1541 33 599 27884 3764 114 3140 2052 1405 47 915\n",
      "2021-08-31 1474 30 498 26710 3561 114 2352 2201 1262 45 828\n",
      "2021-09-01 1410 27 402 25541 3369 114 1606 2361 1134 42 750\n",
      "2021-09-02 1349 24 310 24383 3187 114 901 2533 1019 40 680\n",
      "2021-09-03 1290 21 223 23239 3015 114 233 2718 915 38 617\n",
      "2021-09-04 1234 19 139 22111 2853 114 -397 2916 822 36 561\n",
      "2021-09-05 1180 16 58 21000 2699 114 -994 3128 738 35 510\n",
      "---------------------------------------------------------------------------------------------------------------\n",
      "\n",
      "\n",
      "Predictions using 2021-08-25 parameters\n",
      "\n",
      "------------------------------------------------------------------------------------------------------------------\n",
      "Date Uncovered Du Ru Exposed Infected Di Ri Isolated Apathetic Da Ra\n",
      "--------------------------------------------------------------------------------------------------------------------\n",
      "2021-09-12 10915 558 6623 -36794 1829 -26 13729 717 51481 1056 22537\n",
      "2021-09-13 12018 621 7292 -51004 1730 -30 13916 690 61036 1249 26629\n",
      "2021-09-14 13232 690 8029 -67901 1637 -35 14092 664 72364 1477 31481\n",
      "2021-08-25 1924 51 1174 33640 4966 114 7808 1444 2403 66 1520\n",
      "2021-08-26 2119 62 1292 33111 4698 102 8314 1389 2849 75 1711\n",
      "2021-08-27 2333 74 1421 32467 4444 90 8792 1336 3377 85 1937\n",
      "2021-08-28 2570 87 1565 31686 4204 79 9245 1285 4004 98 2205\n",
      "2021-08-29 2830 102 1722 30738 3977 69 9673 1236 4747 113 2524\n",
      "2021-08-30 3117 118 1896 29591 3763 60 10079 1189 5629 131 2901\n",
      "2021-08-31 3433 136 2087 28205 3559 51 10462 1143 6673 152 3349\n",
      "2021-09-01 3780 155 2298 26533 3367 42 10825 1100 7912 177 3879\n",
      "2021-09-02 4163 177 2529 24522 3186 34 11168 1058 9381 206 4508\n",
      "2021-09-03 4585 201 2785 22105 3014 26 11493 1018 11122 241 5254\n",
      "2021-09-04 5049 227 3066 19206 2851 19 11800 979 13186 283 6138\n",
      "2021-09-05 5560 256 3376 15733 2697 12 12090 942 15634 332 7186\n",
      "---------------------------------------------------------------------------------------------------------------\n",
      "\n",
      "\n",
      "Predictions using 2021-08-26 parameters\n",
      "\n",
      "------------------------------------------------------------------------------------------------------------------\n",
      "Date Uncovered Du Ru Exposed Infected Di Ri Isolated Apathetic Da Ra\n",
      "--------------------------------------------------------------------------------------------------------------------\n",
      "2021-09-12 621 31 1109 1107 35 161 1044 234 1093 97 912\n",
      "2021-09-13 577 30 1103 941 26 161 1030 211 1033 98 884\n",
      "2021-09-14 537 29 1099 807 20 161 1019 190 977 98 859\n",
      "2021-09-15 499 29 1094 698 15 161 1012 171 923 99 834\n",
      "2021-08-26 2119 62 1292 33111 4698 102 8314 1389 2849 75 1711\n",
      "2021-08-27 1972 59 1274 27578 3525 117 6485 1251 2693 77 1640\n",
      "2021-08-28 1835 56 1257 22750 2644 128 5112 1126 2545 78 1572\n",
      "2021-08-29 1707 53 1241 18631 1984 136 4082 1014 2406 80 1509\n",
      "2021-08-30 1588 51 1227 15174 1489 143 3310 913 2274 82 1449\n",
      "2021-08-31 1478 48 1213 12309 1117 147 2730 823 2149 83 1392\n",
      "2021-09-01 1375 46 1201 9958 838 151 2295 741 2032 85 1339\n",
      "2021-09-02 1279 44 1189 8044 628 154 1969 667 1920 86 1288\n",
      "2021-09-03 1190 43 1178 6495 471 156 1724 601 1815 88 1240\n",
      "2021-09-04 1107 41 1168 5249 354 157 1540 541 1716 89 1195\n",
      "2021-09-05 1030 39 1159 4249 265 158 1402 487 1622 90 1152\n",
      "---------------------------------------------------------------------------------------------------------------\n",
      "\n",
      "\n",
      "Predictions using 2021-08-27 parameters\n",
      "\n",
      "------------------------------------------------------------------------------------------------------------------\n",
      "Date Uncovered Du Ru Exposed Infected Di Ri Isolated Apathetic Da Ra\n",
      "--------------------------------------------------------------------------------------------------------------------\n",
      "2021-09-12 1 13 -137 6244 2341 -375 -15113 3 0 -10 197\n",
      "2021-09-13 0 13 -137 5981 2282 -399 -16192 2 0 -10 196\n",
      "2021-09-14 0 13 -137 5748 2224 -423 -17243 1 0 -10 196\n",
      "2021-09-15 0 13 -137 5540 2168 -446 -18268 1 0 -11 196\n",
      "2021-09-16 0 13 -138 5350 2113 -469 -19267 0 0 -11 196\n",
      "2021-08-27 1972 59 1274 27578 3525 117 6485 1251 2693 77 1640\n",
      "2021-08-28 1243 42 751 25129 3436 80 4861 863 1622 42 1066\n",
      "2021-08-29 782 31 421 22399 3349 43 3278 595 976 20 720\n",
      "2021-08-30 492 24 213 19759 3264 8 1734 410 588 8 512\n",
      "2021-08-31 309 20 83 17375 3182 -25 230 283 354 0 386\n",
      "2021-09-01 194 17 0 15305 3101 -58 -1235 195 213 -4 311\n",
      "2021-09-02 122 15 -50 13551 3023 -91 -2664 134 128 -6 265\n",
      "2021-09-03 76 14 -83 12085 2947 -123 -4057 93 77 -8 238\n",
      "2021-09-04 48 14 -103 10869 2872 -154 -5415 64 46 -9 221\n",
      "2021-09-05 30 13 -116 9865 2800 -184 -6738 44 28 -10 211\n",
      "---------------------------------------------------------------------------------------------------------------\n",
      "\n",
      "\n",
      "Predictions using 2021-08-28 parameters\n",
      "\n",
      "------------------------------------------------------------------------------------------------------------------\n",
      "Date Uncovered Du Ru Exposed Infected Di Ri Isolated Apathetic Da Ra\n",
      "--------------------------------------------------------------------------------------------------------------------\n",
      "2021-09-12 455 -55 584 39796 27600 504 80610 3160 146 42 -286\n",
      "2021-09-13 426 -59 578 29576 31712 577 93502 3445 125 42 -306\n",
      "2021-09-14 398 -62 572 14953 36438 660 108316 3757 106 42 -323\n",
      "2021-09-15 373 -66 567 -5236 41868 755 125337 4096 90 42 -337\n",
      "2021-09-16 349 -69 562 -32420 48107 865 144894 4466 77 42 -349\n",
      "2021-09-17 326 -71 557 -68345 55275 991 167365 4870 65 42 -360\n",
      "2021-08-28 1243 42 751 25129 3436 80 4861 863 1622 42 1066\n",
      "2021-08-29 1162 32 734 27921 3948 89 6466 941 1382 42 846\n",
      "2021-08-30 1086 22 718 30822 4536 99 8310 1026 1177 42 658\n",
      "2021-08-31 1015 13 703 33815 5212 111 10429 1118 1003 42 498\n",
      "2021-09-01 949 5 689 36869 5988 124 12863 1219 854 42 362\n",
      "2021-09-02 887 -1 676 39937 6881 140 15661 1330 728 42 246\n",
      "2021-09-03 830 -9 664 42954 7906 158 18875 1450 620 42 148\n",
      "2021-09-04 776 -15 652 45834 9084 179 22569 1581 528 42 63\n",
      "2021-09-05 726 -21 642 48458 10438 203 26812 1724 450 42 -7\n",
      "---------------------------------------------------------------------------------------------------------------\n",
      "\n",
      "\n",
      "Predictions using 2021-08-29 parameters\n",
      "\n",
      "------------------------------------------------------------------------------------------------------------------\n",
      "Date Uncovered Du Ru Exposed Infected Di Ri Isolated Apathetic Da Ra\n",
      "--------------------------------------------------------------------------------------------------------------------\n",
      "2021-09-12 104109 2391 67496 -325191 1742 138 3709 247439 60849 1566 46728\n",
      "2021-09-13 142783 3287 92852 -521859 1643 140 3585 368398 79737 2051 61302\n",
      "2021-09-14 195615 4515 127626 -813438 1550 142 3469 548486 104489 2685 80400\n",
      "2021-09-15 267728 6199 175268 -1244196 1462 144 3359 816609 136925 3517 105426\n",
      "2021-09-16 366137 8503 240472 -1879212 1379 146 3255 1215801 179430 4607 138221\n",
      "2021-09-17 500571 11654 329643 -2814405 1301 148 3157 1810136 235128 6035 181196\n",
      "2021-09-18 684816 15962 451555 -4191641 1227 149 3065 2695006 308117 7907 237511\n",
      "2021-08-29 1162 32 734 27921 3948 89 6466 941 1382 42 846\n",
      "2021-08-30 1606 42 1017 30855 3724 94 6186 1401 1811 53 1177\n",
      "2021-08-31 2219 55 1408 33679 3512 98 5921 2085 2373 67 1610\n",
      "2021-09-01 3065 74 1948 36198 3313 103 5672 3105 3109 86 2179\n",
      "2021-09-02 4234 101 2695 38109 3125 107 5437 4623 4075 111 2923\n",
      "2021-09-03 5846 137 3726 38963 2948 111 5216 6883 5340 143 3900\n",
      "2021-09-04 8069 188 5150 38092 2780 115 5007 10248 6997 185 5179\n",
      "2021-09-05 11132 257 7115 34512 2623 118 4809 15259 9170 241 6855\n",
      "---------------------------------------------------------------------------------------------------------------\n",
      "\n",
      "\n",
      "Predictions using 2021-08-30 parameters\n",
      "\n",
      "------------------------------------------------------------------------------------------------------------------\n",
      "Date Uncovered Du Ru Exposed Infected Di Ri Isolated Apathetic Da Ra\n",
      "--------------------------------------------------------------------------------------------------------------------\n",
      "2021-09-12 18 18 2 6836 966 33 -14768 2 69 1 -244\n",
      "2021-09-13 12 18 -1 6084 871 31 -15492 1 53 0 -256\n",
      "2021-09-14 9 18 -3 5421 785 29 -16145 1 41 0 -266\n",
      "2021-09-15 6 17 -5 4835 708 28 -16733 0 32 0 -273\n",
      "2021-09-16 4 17 -6 4316 638 26 -17263 0 25 0 -279\n",
      "2021-09-17 3 17 -7 3857 575 25 -17741 0 19 0 -284\n",
      "2021-09-18 2 17 -8 3449 518 24 -18172 0 15 0 -288\n",
      "2021-09-19 1 17 -8 3087 467 23 -18561 0 11 0 -290\n",
      "2021-08-30 1606 42 1017 30855 3724 94 6186 1401 1811 53 1177\n",
      "2021-08-31 1141 35 719 28097 3357 86 3397 873 1409 41 849\n",
      "2021-09-01 810 30 507 25301 3026 78 882 543 1096 31 593\n",
      "2021-09-02 575 26 356 22620 2727 72 -1383 338 852 24 395\n",
      "2021-09-03 407 23 250 20131 2459 66 -3426 211 663 18 240\n",
      "2021-09-04 289 22 174 17866 2216 61 -5268 131 516 14 120\n",
      "2021-09-05 205 20 120 15833 1998 56 -6928 82 401 10 27\n",
      "---------------------------------------------------------------------------------------------------------------\n",
      "\n",
      "\n"
     ]
    }
   ],
   "source": [
    "for z in range(1, len(last5_params)+1):\n",
    "    results = seird_model(last5_vals[z-1], last5_params[z-1], t7)\n",
    "    t_z = np.arange(7 - z, 22 - z, 1)\n",
    "    print(\"Predictions using\", last5_vals[z-1][0], \"parameters\\n\")\n",
    "    print(\"------------------------------------------------------------------------------------------------------------------\")\n",
    "    print(\"Date\", \"Uncovered\", \"Du\", \"Ru\", \"Exposed\", \"Infected\", \"Di\", \"Ri\", \"Isolated\", \"Apathetic\", \"Da\", \"Ra\" )\n",
    "    print(\"--------------------------------------------------------------------------------------------------------------------\")\n",
    "    for day in t_z:\n",
    "        print(results[day][0],int(float(results[day][2])),int(float(results[day][3])),int(float(results[day][4])), int(float(results[day][5])), int(float(results[day][6])),int(float(results[day][7])),int(float(results[day][8])),int(float(results[day][9])), int(float(results[day][10])), int(float(results[day][11])),int(float(results[day][12])))\n",
    "    print(\"---------------------------------------------------------------------------------------------------------------\\n\\n\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "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.8.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}