Predictive model.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "d136c3a7",
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import numpy as np\n",
    "import datetime\n",
    "from pyproj import Transformer\n",
    "\n",
    "df = pd.read_csv('./riista-utf8.csv')\n",
    "\n",
    "\n",
    "# Enrichment steps\n",
    "\n",
    "# convert the location data to longitude and latitude info\n",
    "# 3067 is the EPSG code used by the orignal data\n",
    "# 4326 is the Standard WGS84 longitude/latitude coords\n",
    "transformer = Transformer.from_crs(3067, 4326)\n",
    "df['lat'], df['lon'] = np.transpose(np.array([transformer.transform(x, y) for x, y in zip(df['x'], df['y'])]))\n",
    "\n",
    "\n",
    "# write the lat & long info to a file so we can use them in a seprate script to call https://sunrise-sunset.org/api\n",
    "# since we need to do 60k+ API calls, it will take some time. That's why it's being done on a separate server and not\n",
    "# on this Jupyter notebook\n",
    "df.to_csv('coordinates.csv', columns=['tapahtumaAika', 'lat', 'lon'], header=False)\n",
    "\n",
    "# read back the results and join with the original data set\n",
    "timeData = pd.read_csv('./daylight-api-output.txt', header=None, names=['id','tapahtumaAika', 'valoisuus'])\n",
    "df['valoisuus'] = timeData['valoisuus'] * -1 # Reverse the mapping codes used in the API fetcher script\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "8dbd5549",
   "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>id</th>\n",
       "      <th>tapahtumaAika</th>\n",
       "      <th>vuosi</th>\n",
       "      <th>kuukausi</th>\n",
       "      <th>x</th>\n",
       "      <th>y</th>\n",
       "      <th>kunta</th>\n",
       "      <th>kuntaNimi</th>\n",
       "      <th>maakunta</th>\n",
       "      <th>maakuntaNimi</th>\n",
       "      <th>tielaji</th>\n",
       "      <th>tielajis</th>\n",
       "      <th>tieNumero</th>\n",
       "      <th>tieYllapito</th>\n",
       "      <th>tieYllapitoNimi</th>\n",
       "      <th>riistalaji</th>\n",
       "      <th>riistalajiNimi</th>\n",
       "      <th>lat</th>\n",
       "      <th>lon</th>\n",
       "      <th>valoisuus</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1</td>\n",
       "      <td>2017-12-10T16:55:00.000+02:00</td>\n",
       "      <td>2017</td>\n",
       "      <td>12</td>\n",
       "      <td>278773.65900</td>\n",
       "      <td>6.980878e+06</td>\n",
       "      <td>743</td>\n",
       "      <td>Seinäjoki</td>\n",
       "      <td>14</td>\n",
       "      <td>Etelä-Pohjanmaa</td>\n",
       "      <td>3</td>\n",
       "      <td>Valtatie</td>\n",
       "      <td>18.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>Valtio</td>\n",
       "      <td>47507</td>\n",
       "      <td>Metsäkauris</td>\n",
       "      <td>62.890515</td>\n",
       "      <td>22.646311</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2</td>\n",
       "      <td>2017-09-03T13:00:00.000+03:00</td>\n",
       "      <td>2017</td>\n",
       "      <td>9</td>\n",
       "      <td>397501.00800</td>\n",
       "      <td>7.503084e+06</td>\n",
       "      <td>261</td>\n",
       "      <td>Kittilä</td>\n",
       "      <td>19</td>\n",
       "      <td>Lappi</td>\n",
       "      <td>4</td>\n",
       "      <td>Kantatie</td>\n",
       "      <td>80.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>Valtio</td>\n",
       "      <td>47503</td>\n",
       "      <td>Hirvi</td>\n",
       "      <td>67.625279</td>\n",
       "      <td>24.586614</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>3</td>\n",
       "      <td>2017-09-18T07:00:00.000+03:00</td>\n",
       "      <td>2017</td>\n",
       "      <td>9</td>\n",
       "      <td>418353.84100</td>\n",
       "      <td>7.514105e+06</td>\n",
       "      <td>261</td>\n",
       "      <td>Kittilä</td>\n",
       "      <td>19</td>\n",
       "      <td>Lappi</td>\n",
       "      <td>6</td>\n",
       "      <td>Muu maantie</td>\n",
       "      <td>9552.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>Valtio</td>\n",
       "      <td>47503</td>\n",
       "      <td>Hirvi</td>\n",
       "      <td>67.730620</td>\n",
       "      <td>25.069129</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>4</td>\n",
       "      <td>2017-12-01T17:00:00.000+02:00</td>\n",
       "      <td>2017</td>\n",
       "      <td>12</td>\n",
       "      <td>488970.82700</td>\n",
       "      <td>6.722504e+06</td>\n",
       "      <td>285</td>\n",
       "      <td>Kotka</td>\n",
       "      <td>8</td>\n",
       "      <td>Kymenlaakso</td>\n",
       "      <td>5</td>\n",
       "      <td>Seututie</td>\n",
       "      <td>357.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>Valtio</td>\n",
       "      <td>47629</td>\n",
       "      <td>Valkohäntäpeura</td>\n",
       "      <td>60.638182</td>\n",
       "      <td>26.798369</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>5</td>\n",
       "      <td>2017-09-18T07:00:00.000+03:00</td>\n",
       "      <td>2017</td>\n",
       "      <td>9</td>\n",
       "      <td>418766.58000</td>\n",
       "      <td>7.506931e+06</td>\n",
       "      <td>261</td>\n",
       "      <td>Kittilä</td>\n",
       "      <td>19</td>\n",
       "      <td>Lappi</td>\n",
       "      <td>6</td>\n",
       "      <td>Muu maantie</td>\n",
       "      <td>9552.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>Valtio</td>\n",
       "      <td>47503</td>\n",
       "      <td>Hirvi</td>\n",
       "      <td>67.666425</td>\n",
       "      <td>25.084131</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>63937</th>\n",
       "      <td>63938</td>\n",
       "      <td>2021-10-04T09:30:00.000+03:00</td>\n",
       "      <td>2021</td>\n",
       "      <td>10</td>\n",
       "      <td>236501.85584</td>\n",
       "      <td>6.691606e+06</td>\n",
       "      <td>445</td>\n",
       "      <td>Parainen</td>\n",
       "      <td>2</td>\n",
       "      <td>Varsinais-Suomi</td>\n",
       "      <td>5</td>\n",
       "      <td>Seututie</td>\n",
       "      <td>180.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>Valtio</td>\n",
       "      <td>47507</td>\n",
       "      <td>Metsäkauris</td>\n",
       "      <td>60.275376</td>\n",
       "      <td>22.233479</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>63938</th>\n",
       "      <td>63939</td>\n",
       "      <td>2021-09-20T19:51:48.017+03:00</td>\n",
       "      <td>2021</td>\n",
       "      <td>9</td>\n",
       "      <td>515750.42581</td>\n",
       "      <td>6.946907e+06</td>\n",
       "      <td>778</td>\n",
       "      <td>Suonenjoki</td>\n",
       "      <td>11</td>\n",
       "      <td>Pohjois-Savo</td>\n",
       "      <td>3</td>\n",
       "      <td>Valtatie</td>\n",
       "      <td>9.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>Valtio</td>\n",
       "      <td>47503</td>\n",
       "      <td>Hirvi</td>\n",
       "      <td>62.652474</td>\n",
       "      <td>27.307304</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>63939</th>\n",
       "      <td>63940</td>\n",
       "      <td>2021-05-29T06:00:00.000+03:00</td>\n",
       "      <td>2021</td>\n",
       "      <td>5</td>\n",
       "      <td>245406.95714</td>\n",
       "      <td>6.712757e+06</td>\n",
       "      <td>202</td>\n",
       "      <td>Kaarina</td>\n",
       "      <td>2</td>\n",
       "      <td>Varsinais-Suomi</td>\n",
       "      <td>4</td>\n",
       "      <td>Kantatie</td>\n",
       "      <td>40.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>Valtio</td>\n",
       "      <td>47507</td>\n",
       "      <td>Metsäkauris</td>\n",
       "      <td>60.470338</td>\n",
       "      <td>22.367095</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>63940</th>\n",
       "      <td>63941</td>\n",
       "      <td>2021-04-13T14:30:00.000+03:00</td>\n",
       "      <td>2021</td>\n",
       "      <td>4</td>\n",
       "      <td>238690.96803</td>\n",
       "      <td>6.705560e+06</td>\n",
       "      <td>853</td>\n",
       "      <td>Turku</td>\n",
       "      <td>2</td>\n",
       "      <td>Varsinais-Suomi</td>\n",
       "      <td>7</td>\n",
       "      <td>Puuttuu</td>\n",
       "      <td>NaN</td>\n",
       "      <td>99.0</td>\n",
       "      <td>Ei tietoa</td>\n",
       "      <td>47507</td>\n",
       "      <td>Metsäkauris</td>\n",
       "      <td>60.401653</td>\n",
       "      <td>22.254785</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>63941</th>\n",
       "      <td>63942</td>\n",
       "      <td>2021-06-12T14:15:00.000+03:00</td>\n",
       "      <td>2021</td>\n",
       "      <td>6</td>\n",
       "      <td>237601.43474</td>\n",
       "      <td>7.012767e+06</td>\n",
       "      <td>499</td>\n",
       "      <td>Mustasaari</td>\n",
       "      <td>15</td>\n",
       "      <td>Pohjanmaa</td>\n",
       "      <td>3</td>\n",
       "      <td>Valtatie</td>\n",
       "      <td>8.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>Valtio</td>\n",
       "      <td>47507</td>\n",
       "      <td>Metsäkauris</td>\n",
       "      <td>63.148286</td>\n",
       "      <td>21.788953</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>63942 rows × 20 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "          id                  tapahtumaAika  vuosi  kuukausi             x  \\\n",
       "0          1  2017-12-10T16:55:00.000+02:00   2017        12  278773.65900   \n",
       "1          2  2017-09-03T13:00:00.000+03:00   2017         9  397501.00800   \n",
       "2          3  2017-09-18T07:00:00.000+03:00   2017         9  418353.84100   \n",
       "3          4  2017-12-01T17:00:00.000+02:00   2017        12  488970.82700   \n",
       "4          5  2017-09-18T07:00:00.000+03:00   2017         9  418766.58000   \n",
       "...      ...                            ...    ...       ...           ...   \n",
       "63937  63938  2021-10-04T09:30:00.000+03:00   2021        10  236501.85584   \n",
       "63938  63939  2021-09-20T19:51:48.017+03:00   2021         9  515750.42581   \n",
       "63939  63940  2021-05-29T06:00:00.000+03:00   2021         5  245406.95714   \n",
       "63940  63941  2021-04-13T14:30:00.000+03:00   2021         4  238690.96803   \n",
       "63941  63942  2021-06-12T14:15:00.000+03:00   2021         6  237601.43474   \n",
       "\n",
       "                  y  kunta   kuntaNimi  maakunta     maakuntaNimi  tielaji  \\\n",
       "0      6.980878e+06    743   Seinäjoki        14  Etelä-Pohjanmaa        3   \n",
       "1      7.503084e+06    261     Kittilä        19            Lappi        4   \n",
       "2      7.514105e+06    261     Kittilä        19            Lappi        6   \n",
       "3      6.722504e+06    285       Kotka         8      Kymenlaakso        5   \n",
       "4      7.506931e+06    261     Kittilä        19            Lappi        6   \n",
       "...             ...    ...         ...       ...              ...      ...   \n",
       "63937  6.691606e+06    445    Parainen         2  Varsinais-Suomi        5   \n",
       "63938  6.946907e+06    778  Suonenjoki        11     Pohjois-Savo        3   \n",
       "63939  6.712757e+06    202     Kaarina         2  Varsinais-Suomi        4   \n",
       "63940  6.705560e+06    853       Turku         2  Varsinais-Suomi        7   \n",
       "63941  7.012767e+06    499  Mustasaari        15        Pohjanmaa        3   \n",
       "\n",
       "          tielajis  tieNumero  tieYllapito tieYllapitoNimi  riistalaji  \\\n",
       "0         Valtatie       18.0          1.0          Valtio       47507   \n",
       "1         Kantatie       80.0          1.0          Valtio       47503   \n",
       "2      Muu maantie     9552.0          1.0          Valtio       47503   \n",
       "3         Seututie      357.0          1.0          Valtio       47629   \n",
       "4      Muu maantie     9552.0          1.0          Valtio       47503   \n",
       "...            ...        ...          ...             ...         ...   \n",
       "63937     Seututie      180.0          1.0          Valtio       47507   \n",
       "63938     Valtatie        9.0          1.0          Valtio       47503   \n",
       "63939     Kantatie       40.0          1.0          Valtio       47507   \n",
       "63940      Puuttuu        NaN         99.0       Ei tietoa       47507   \n",
       "63941     Valtatie        8.0          1.0          Valtio       47507   \n",
       "\n",
       "        riistalajiNimi        lat        lon  valoisuus  \n",
       "0          Metsäkauris  62.890515  22.646311          3  \n",
       "1                Hirvi  67.625279  24.586614          0  \n",
       "2                Hirvi  67.730620  25.069129          0  \n",
       "3      Valkohäntäpeura  60.638182  26.798369          3  \n",
       "4                Hirvi  67.666425  25.084131          0  \n",
       "...                ...        ...        ...        ...  \n",
       "63937      Metsäkauris  60.275376  22.233479          0  \n",
       "63938            Hirvi  62.652474  27.307304          3  \n",
       "63939      Metsäkauris  60.470338  22.367095          0  \n",
       "63940      Metsäkauris  60.401653  22.254785          0  \n",
       "63941      Metsäkauris  63.148286  21.788953          0  \n",
       "\n",
       "[63942 rows x 20 columns]"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 360,
   "id": "0d16914c",
   "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>x</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>riistalajiNimi</th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>Hirvi</th>\n",
       "      <td>9192</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Kuusipeura</th>\n",
       "      <td>87</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Metsäkauris</th>\n",
       "      <td>22464</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Metsäpeura</th>\n",
       "      <td>325</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Valkohäntäpeura</th>\n",
       "      <td>31702</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Villisika</th>\n",
       "      <td>168</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                     x\n",
       "riistalajiNimi        \n",
       "Hirvi             9192\n",
       "Kuusipeura          87\n",
       "Metsäkauris      22464\n",
       "Metsäpeura         325\n",
       "Valkohäntäpeura  31702\n",
       "Villisika          168"
      ]
     },
     "execution_count": 360,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Check the distribution of different types of game animal crashes\n",
    "df[['riistalajiNimi', 'x']].groupby('riistalajiNimi').count()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 363,
   "id": "e760ea58",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "<class 'pandas.core.frame.DataFrame'>\n",
      "RangeIndex: 63942 entries, 0 to 63941\n",
      "Data columns (total 3 columns):\n",
      " #   Column     Non-Null Count  Dtype   \n",
      "---  ------     --------------  -----   \n",
      " 0   kuukausi   63942 non-null  category\n",
      " 1   kunta      63942 non-null  category\n",
      " 2   valoisuus  63942 non-null  category\n",
      "dtypes: category(3)\n",
      "memory usage: 261.0 KB\n",
      "LR  accuracy: 0.6162326999765423\n"
     ]
    }
   ],
   "source": [
    "from sklearn.model_selection import train_test_split\n",
    "\n",
    "from sklearn.metrics import accuracy_score\n",
    "from sklearn.pipeline import Pipeline\n",
    "from sklearn.ensemble import RandomForestClassifier\n",
    "\n",
    "#model to predict what animal is the most likely to be encountered in a crash on given month and location\n",
    "#drop unneeded dimensions, drop NaNs\n",
    "\n",
    "data = df.drop(['id', 'tielaji', 'tapahtumaAika','tieNumero', 'maakuntaNimi', 'vuosi', 'tielajis','x','y','lat','lon','maakunta', 'riistalajiNimi', 'kuntaNimi', 'tieYllapito', 'tieYllapitoNimi'], axis=1)\n",
    "# \n",
    "data['kuukausi'] = data['kuukausi'].astype('category')\n",
    "data['valoisuus'] = data['valoisuus'].astype('category')\n",
    "data['kunta'] = data['kunta'].astype('category')\n",
    "\n",
    "\n",
    "data = data.dropna()\n",
    "\n",
    "# extract the target category\n",
    "target = data['riistalaji'].astype('category')\n",
    "data = data.drop(['riistalaji'], axis=1)\n",
    "\n",
    "data.info()\n",
    "\n",
    "training_data, test_data, train_target, test_target = train_test_split(data, target, train_size=0.8)\n",
    "\n",
    "# Since the classes are highly imbalanced, let's use random forest classifier instead of logistic regression\n",
    "# With the RF we have bit more accuracy, at least\n",
    "rf = RandomForestClassifier(random_state=123, n_estimators=200)\n",
    "\n",
    "pipeline = Pipeline([('rf',rf)])\n",
    "\n",
    "model = pipeline.fit(training_data, train_target)\n",
    "predictedLr = model.predict(test_data)\n",
    "\n",
    "print('LR  accuracy:', accuracy_score(test_target, predictedLr))\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "9b637ec7",
   "metadata": {},
   "outputs": [],
   "source": [
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 356,
   "id": "228d6522",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "     dark             astr dawn        nautical dawn    civic dawn       daytime          civic twilight    nautical twilight    astr twilight    night\n",
      "---  ---------------  ---------------  ---------------  ---------------  ---------------  ----------------  -------------------  ---------------  ---------------\n",
      "Jan  Valkohäntäpeura  Valkohäntäpeura  Valkohäntäpeura  Valkohäntäpeura  Valkohäntäpeura  Valkohäntäpeura   Hirvi                Valkohäntäpeura  Valkohäntäpeura\n",
      "Feb  Valkohäntäpeura  Hirvi            Valkohäntäpeura  Valkohäntäpeura  Valkohäntäpeura  Valkohäntäpeura   Metsäkauris          Hirvi            Valkohäntäpeura\n",
      "Mar  Hirvi            Valkohäntäpeura  Valkohäntäpeura  Valkohäntäpeura  Valkohäntäpeura  Valkohäntäpeura   Metsäkauris          Metsäkauris      Valkohäntäpeura\n",
      "Apr  Valkohäntäpeura  Valkohäntäpeura  Valkohäntäpeura  Valkohäntäpeura  Valkohäntäpeura  Valkohäntäpeura   Metsäkauris          Valkohäntäpeura  Valkohäntäpeura\n",
      "May  Valkohäntäpeura  Valkohäntäpeura  Hirvi            Metsäkauris      Metsäkauris      Metsäkauris       Valkohäntäpeura      Valkohäntäpeura  Valkohäntäpeura\n",
      "Jun  Valkohäntäpeura  Valkohäntäpeura  Metsäkauris      Metsäkauris      Metsäkauris      Metsäkauris       Valkohäntäpeura      Valkohäntäpeura  Valkohäntäpeura\n",
      "Jul  Valkohäntäpeura  Hirvi            Hirvi            Villisika        Valkohäntäpeura  Valkohäntäpeura   Metsäkauris          Valkohäntäpeura  Metsäkauris\n",
      "Aug  Valkohäntäpeura  Hirvi            Metsäkauris      Hirvi            Metsäkauris      Metsäkauris       Metsäkauris          Valkohäntäpeura  Valkohäntäpeura\n",
      "Sep  Valkohäntäpeura  Hirvi            Metsäkauris      Valkohäntäpeura  Hirvi            Valkohäntäpeura   Metsäkauris          Valkohäntäpeura  Hirvi\n",
      "Oct  Valkohäntäpeura  Hirvi            Valkohäntäpeura  Valkohäntäpeura  Valkohäntäpeura  Valkohäntäpeura   Valkohäntäpeura      Metsäkauris      Valkohäntäpeura\n",
      "Nov  Valkohäntäpeura  Valkohäntäpeura  Valkohäntäpeura  Valkohäntäpeura  Valkohäntäpeura  Valkohäntäpeura   Valkohäntäpeura      Valkohäntäpeura  Valkohäntäpeura\n",
      "Dec  Valkohäntäpeura  Metsäkauris      Valkohäntäpeura  Valkohäntäpeura  Valkohäntäpeura  Valkohäntäpeura   Valkohäntäpeura      Valkohäntäpeura  Valkohäntäpeura\n"
     ]
    }
   ],
   "source": [
    "from tabulate import tabulate\n",
    "import datetime\n",
    "\n",
    "location = 'Sipoo'\n",
    "kunta = df['kunta'].where(df['kuntaNimi'] == location).dropna().head(1).values[0]\n",
    "\n",
    "output = [];\n",
    "for month in range(1,13):\n",
    "    monthData = []\n",
    "    datetime_object = datetime.datetime.strptime(str(month), \"%m\")\n",
    "    monthData.append(datetime_object.strftime(\"%b\"))\n",
    "\n",
    "    for light in range(-4,5):\n",
    "        data_to_predict = pd.DataFrame({'kuukausi': [month], 'kunta': [kunta], 'valoisuus': [light]})\n",
    "        p = lr.predict(data_to_predict)\n",
    "        game = df['riistalajiNimi'].where(df['riistalaji'] == p[0]).dropna().head(1).values[0]\n",
    "        monthData.append(game)\n",
    "    output.append(monthData)\n",
    "    \n",
    "    \n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 361,
   "id": "086b47c4",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "|     | dark            | astr dawn       | nautical dawn   | civic dawn      | daytime         | civic twilight   | nautical twilight   | astr twilight   | night           |\n",
      "|-----|-----------------|-----------------|-----------------|-----------------|-----------------|------------------|---------------------|-----------------|-----------------|\n",
      "| Jan | Valkohäntäpeura | Valkohäntäpeura | Valkohäntäpeura | Valkohäntäpeura | Valkohäntäpeura | Valkohäntäpeura  | Hirvi               | Valkohäntäpeura | Valkohäntäpeura |\n",
      "| Feb | Valkohäntäpeura | Hirvi           | Valkohäntäpeura | Valkohäntäpeura | Valkohäntäpeura | Valkohäntäpeura  | Metsäkauris         | Hirvi           | Valkohäntäpeura |\n",
      "| Mar | Hirvi           | Valkohäntäpeura | Valkohäntäpeura | Valkohäntäpeura | Valkohäntäpeura | Valkohäntäpeura  | Metsäkauris         | Metsäkauris     | Valkohäntäpeura |\n",
      "| Apr | Valkohäntäpeura | Valkohäntäpeura | Valkohäntäpeura | Valkohäntäpeura | Valkohäntäpeura | Valkohäntäpeura  | Metsäkauris         | Valkohäntäpeura | Valkohäntäpeura |\n",
      "| May | Valkohäntäpeura | Valkohäntäpeura | Hirvi           | Metsäkauris     | Metsäkauris     | Metsäkauris      | Valkohäntäpeura     | Valkohäntäpeura | Valkohäntäpeura |\n",
      "| Jun | Valkohäntäpeura | Valkohäntäpeura | Metsäkauris     | Metsäkauris     | Metsäkauris     | Metsäkauris      | Valkohäntäpeura     | Valkohäntäpeura | Valkohäntäpeura |\n",
      "| Jul | Valkohäntäpeura | Hirvi           | Hirvi           | Villisika       | Valkohäntäpeura | Valkohäntäpeura  | Metsäkauris         | Valkohäntäpeura | Metsäkauris     |\n",
      "| Aug | Valkohäntäpeura | Hirvi           | Metsäkauris     | Hirvi           | Metsäkauris     | Metsäkauris      | Metsäkauris         | Valkohäntäpeura | Valkohäntäpeura |\n",
      "| Sep | Valkohäntäpeura | Hirvi           | Metsäkauris     | Valkohäntäpeura | Hirvi           | Valkohäntäpeura  | Metsäkauris         | Valkohäntäpeura | Hirvi           |\n",
      "| Oct | Valkohäntäpeura | Hirvi           | Valkohäntäpeura | Valkohäntäpeura | Valkohäntäpeura | Valkohäntäpeura  | Valkohäntäpeura     | Metsäkauris     | Valkohäntäpeura |\n",
      "| Nov | Valkohäntäpeura | Valkohäntäpeura | Valkohäntäpeura | Valkohäntäpeura | Valkohäntäpeura | Valkohäntäpeura  | Valkohäntäpeura     | Valkohäntäpeura | Valkohäntäpeura |\n",
      "| Dec | Valkohäntäpeura | Metsäkauris     | Valkohäntäpeura | Valkohäntäpeura | Valkohäntäpeura | Valkohäntäpeura  | Valkohäntäpeura     | Valkohäntäpeura | Valkohäntäpeura |\n"
     ]
    }
   ],
   "source": [
    "print(tabulate(output, headers=['dark', 'astr dawn', 'nautical dawn', 'civic dawn', 'daytime','civic twilight', 'nautical twilight','astr twilight', 'night'], tablefmt=\"github\"))\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "d17389a1",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.9.13"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}