choppies_on_the_map.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "b44551c5",
   "metadata": {},
   "source": [
    "# Choppies Locations"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 224,
   "id": "60ebd8f0",
   "metadata": {},
   "outputs": [],
   "source": [
    "choppies <- read.csv(\"choppies.csv\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 225,
   "id": "2a259bca",
   "metadata": {},
   "outputs": [],
   "source": [
    "choppies$Date <- as.Date(choppies$Date, \"%d/%m/%Y\")\n",
    "choppies$year <- as.numeric(format(choppies$Date, \"%Y\"))\n",
    "choppies$District <- factor(choppies$District)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 226,
   "id": "6e4b0bdc",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "'data.frame':\t108 obs. of  11 variables:\n",
      " $ Name      : chr  \"Choppies Bobonong\" \"Choppies Letlhakane\" \"Choppies Mahalapye\" \"Choppies Mahalapye\" ...\n",
      " $ District  : Factor w/ 13 levels \"Central\",\"Chobe\",..: 1 1 1 1 1 1 1 1 1 1 ...\n",
      " $ Town      : chr  \"Bobonong\" \"Letlhakane\" \"Mahalapye\" \"Mahalapye\" ...\n",
      " $ Address   : chr  \"Moilamba Ward,Cash Bazaar building\" \"Nkosho Ward,Letlhakane\" \"Plot 6042, Main Mall, Mahalapye\" \"Watershed Mall, A1,  Mahalapye\" ...\n",
      " $ latitude  : num  -22 -21.4 -23.1 -23.1 NA ...\n",
      " $ longitude : num  28.4 25.6 26.8 26.8 NA ...\n",
      " $ Date      : Date, format: \"2012-09-12\" \"2007-10-12\" ...\n",
      " $ Cancelled : logi  NA NA NA NA NA NA ...\n",
      " $ Registered: chr  \"Y\" \"Y\" \"Y\" \"Y\" ...\n",
      " $ Store     : chr  \"choppies\" \"choppies\" \"choppies\" \"choppies\" ...\n",
      " $ year      : num  2012 2007 2012 NA 2012 ...\n"
     ]
    }
   ],
   "source": [
    "str(choppies)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 227,
   "id": "265321ba",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "     Name                 District      Town             Address         \n",
       " Length:108         South-East:42   Length:108         Length:108        \n",
       " Class :character   Kweneng   :18   Class :character   Class :character  \n",
       " Mode  :character   Central   :17   Mode  :character   Mode  :character  \n",
       "                    Southern  : 9                                        \n",
       "                    North-East: 7                                        \n",
       "                    North-West: 4                                        \n",
       "                    (Other)   :11                                        \n",
       "    latitude        longitude          Date            Cancelled     \n",
       " Min.   :-26.02   Min.   :21.64   Min.   :1999-11-02   Mode:logical  \n",
       " 1st Qu.:-24.66   1st Qu.:25.55   1st Qu.:2006-03-16   NA's:108      \n",
       " Median :-24.61   Median :25.87   Median :2012-06-22                 \n",
       " Mean   :-23.48   Mean   :25.81   Mean   :2011-09-14                 \n",
       " 3rd Qu.:-21.98   3rd Qu.:26.11   3rd Qu.:2014-10-28                 \n",
       " Max.   :-17.80   Max.   :28.42   Max.   :2021-10-28                 \n",
       " NA's   :29       NA's   :29      NA's   :10                         \n",
       "  Registered           Store                year     \n",
       " Length:108         Length:108         Min.   :1999  \n",
       " Class :character   Class :character   1st Qu.:2006  \n",
       " Mode  :character   Mode  :character   Median :2012  \n",
       "                                       Mean   :2011  \n",
       "                                       3rd Qu.:2014  \n",
       "                                       Max.   :2021  \n",
       "                                       NA's   :10    "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "summary(choppies)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 228,
   "id": "5aed08bb",
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "library(dplyr)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "3d86106f",
   "metadata": {},
   "source": [
    "Removing NAs..."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 230,
   "id": "8385e9c8",
   "metadata": {},
   "outputs": [],
   "source": [
    "choppies_map <- choppies %>% \n",
    "filter(!is.na(longitude),!is.na(latitude), !is.na(Date)) %>%\n",
    "data.frame"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 231,
   "id": "ae600eb1",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "'data.frame':\t74 obs. of  11 variables:\n",
      " $ Name      : chr  \"Choppies Bobonong\" \"Choppies Letlhakane\" \"Choppies Mahalapye\" \"Choppies Palapye\" ...\n",
      " $ District  : Factor w/ 13 levels \"Central\",\"Chobe\",..: 1 1 1 1 1 1 1 1 1 1 ...\n",
      " $ Town      : chr  \"Bobonong\" \"Letlhakane\" \"Mahalapye\" \"Palapye\" ...\n",
      " $ Address   : chr  \"Moilamba Ward,Cash Bazaar building\" \"Nkosho Ward,Letlhakane\" \"Plot 6042, Main Mall, Mahalapye\" \"Makwapa complex,Serorome Ward,Palapye\" ...\n",
      " $ latitude  : num  -22 -21.4 -23.1 -22.5 -22.5 ...\n",
      " $ longitude : num  28.4 25.6 26.8 27.1 27.1 ...\n",
      " $ Date      : Date, format: \"2012-09-12\" \"2007-10-12\" ...\n",
      " $ Cancelled : logi  NA NA NA NA NA NA ...\n",
      " $ Registered: chr  \"Y\" \"Y\" \"Y\" \"Y\" ...\n",
      " $ Store     : chr  \"choppies\" \"choppies\" \"choppies\" \"choppies\" ...\n",
      " $ year      : num  2012 2007 2012 2005 2020 ...\n"
     ]
    }
   ],
   "source": [
    "str(choppies_map)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 232,
   "id": "4b8880bb",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "     Name                 District      Town             Address         \n",
       " Length:74          South-East:24   Length:74          Length:74         \n",
       " Class :character   Central   :14   Class :character   Class :character  \n",
       " Mode  :character   Kweneng   :12   Mode  :character   Mode  :character  \n",
       "                    Southern  : 8                                        \n",
       "                    North-East: 4                                        \n",
       "                    North-West: 4                                        \n",
       "                    (Other)   : 8                                        \n",
       "    latitude        longitude          Date            Cancelled     \n",
       " Min.   :-26.02   Min.   :21.64   Min.   :1999-11-02   Mode:logical  \n",
       " 1st Qu.:-24.66   1st Qu.:25.52   1st Qu.:2005-09-13   NA's:74       \n",
       " Median :-24.61   Median :25.88   Median :2011-08-06                 \n",
       " Mean   :-23.42   Mean   :25.80   Mean   :2010-12-12                 \n",
       " 3rd Qu.:-21.97   3rd Qu.:26.13   3rd Qu.:2013-11-11                 \n",
       " Max.   :-17.80   Max.   :28.42   Max.   :2021-10-28                 \n",
       "                                                                     \n",
       "  Registered           Store                year     \n",
       " Length:74          Length:74          Min.   :1999  \n",
       " Class :character   Class :character   1st Qu.:2005  \n",
       " Mode  :character   Mode  :character   Median :2011  \n",
       "                                       Mean   :2010  \n",
       "                                       3rd Qu.:2013  \n",
       "                                       Max.   :2021  \n",
       "                                                     "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "summary(choppies_map)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 233,
   "id": "190334ac",
   "metadata": {},
   "outputs": [],
   "source": [
    "library(sf)\n",
    "library(rnaturalearth)\n",
    "library(rnaturalearthdata)\n",
    "library(ggplot2)\n",
    "library(plotly)\n",
    "library(viridis)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 234,
   "id": "bc219d45",
   "metadata": {},
   "outputs": [],
   "source": [
    "#install.packages(\"devtools\") # I guess you also need this\n",
    "#devtools::install_github(\"ropensci/rnaturalearthhires\")\n",
    "library(\"rnaturalearth\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 235,
   "id": "6621087c",
   "metadata": {},
   "outputs": [],
   "source": [
    "BW_states <- ne_states(country = 'botswana', returnclass = 'sf')\n",
    "#str(BW_states)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 240,
   "id": "bbfdbb04",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0gAAANICAMAAADKOT/pAAAB3VBMVEUAAAAem4ojnYglhY4m\nno4on4UqoIYroIMrsH8tcI4wopIxMTEzMzM0nnc0pIA4WYw7gpY7o3g9fZY9p3c9qHc+qHY+\nqHg/qXU/qXZAiW5BqX1DPoVDqYdDqodEAVRHqpxIE1VIIXNJqW5MLG5MmKBMpmxNKVZNTU1Q\nOmdRjaBRu5RRxWpSiKBTn2ZTr2ZVlGJXxm9ZWVlasmRatGZbdZ5fyHVht2litl9it2BkYJhl\nL3FluF1luGFmuV5nuV5oSYpoaGhoulxpSotrum1tbW1vu4Nvy4N4vLN8fHx8sLV/yK2ApbWA\nw1WF1UqGmLSH1U2J1VCLx1qLyFqMibCMjIyNaJWO1liPeaaSfKmTozuU06GY2Gqampqhz2Si\n0jmjzjak0jqmyzSnp6ep1Dyp1Tyry1utxC6w3JCysrKzpzyz2EC01yy02S211iu32i650yi7\n3DC9vb293Si+2iW/3SW/3iW/3inA3STA3iXB1GXB3iPB3iTB3yTB3ybC3yPC3yTC3yXC3yfD\n3yfD3yvE3yvF3zXF4CzG4DXHx8fJ4ErK4UzP5VDQ0NDS4nrU5H3Y2j3Z2dncyi/h4eHi3kLl\n5eXp6enr6+vs4TDw8PDy5nzy8vL450z75zf95yX////uxBpcAAAACXBIWXMAABJ0AAASdAHe\nZh94AAAgAElEQVR4nO2djZvltnWf6TppJVuu0jZNmirRSo0jupKbho3V1ls3mvhDupFrJ65H\nstcap/WoWa8U29JIiZq08UrKers76qwyns5sxlr7/q3lNwEQAAHygMTh/b3Ps/eDl5cHwD3v\nkgA5YLIFAEwmWboAAKwBiAQAARAJAAIgEgAEQCQACIBIABAAkQAgACIBQABEAoAAiAQAARAJ\nAAIgEgAEQCQACIBIABAAkQAgACIBQABEAoAAiAQAARAJAAIgEgAEQCQACIBIABAAkQAgACIB\nQABEAoAAiAQAARAJAAIgEgAEQCQACIBIABAAkQAgACIBQABEAoAAiAQAARAJAAIgEgAEQCQA\nCIBIABAAkQAgACIBQABEAoAAiAQAARAJAAIgEgAEQCQACIBIABAAkQAgACIBQABEAoAAiAQA\nARAJAAIgEgAEQCQACIBIABAAkQAgACIBQABEAoAAiAQAARAJAAIgEgAEQCQACIBIABAAkQAg\nACIBQAA3kRK/Al8ejPhS8bUk2RPe391PkmT/ZEwJ3NY/OUjziIeXVXC/APY4ScXejf5HIwMB\nHSsXqVrdW6TCm0377iyts3EzZmMu6+/XAZKTEQHscZotC/XxKRlwhFtbjhJpRJSz7s1ZvjM6\nzZ9PN8n++E3aOErSwqCLozIubYB6a2eb5NDwEaCAW1vOJZLwZi85ql9tij1GgPRLk4vqxVFy\nEEik7UWSmj4CBHBrS+nHP8t7FulBvfM42STJQdnJ2B7u5Ucyp9v6uKb7X7lbPV9ymCZ7J9qN\nJfXX6u2W+6FqleJ/9fyTo7TavLrJYvmJ5nVBvj+rF+RfSvIvneoqlb9qg0vbvtwri2HZSvml\nfr3abfcaoQnUbOg0KftMp0m5zYPiqW1IecPdYtDAWaQToWOxPSxflv/rpu1iSSRx9TwRuq+q\nG5NF2k9E3YpPD8rPT3ubvFG+Puy9Lr51N2kXnIrFLpEOu5rg8rb3y+/atlKLpNarqchlNXqi\nVlPcUJpUDVkWpmjKVCxBt+FUiQy2rEXK+y43LreXuUEXxZu8m3FZZuSNpBihulH2roXBBmn1\nXLnT7eV+1wOXPpV0TZNLpQRFpPKr6ibv5omaVt0c8XW5+eJ/8LPy2HAvyT/LneiGBS+KncTd\nszZAr0B5Gl8ObaUWSamX0Ec66W+1+KTb0I1ak+J/o5O8EeWGbDcsLgYNjEWq/+vMD0IOi4ei\nI1P+t7tXZX5zqNM8SqtX/59edluTPpVEUnsSwleVTVYdnbT3utx8NbJdHJ9puiYX1bDdvnAs\nqGz7dDu0lbq+Sr26Ubsbmq3K9bsoju0ukr2i8Af5g9yQ7YbFxaCBW2sIP99ena0Xgj01Zyc3\nNqpI0ur1ZrqtSZ8OiNQ+yZus/4feV19XwRuKY8V893Oh1OuyLLGwF9UV17oVob5bnUj7Z/qt\nChvaJMWe5m7hXF0DpSHb/VuzGDRwaw0pRYQX4s96lLa9HEEk3epSL194IW5tr3do1w9avtlr\n1ldeVxvsFLgoi7d3tO1xuilzWN5RSMW1bMUoUvl0mTfKqXarwobu5rudvMh56U/K471eQ1ZP\nwmLQwK01HEQ6yv9jP7x7QSXSgdirvpC+qhVpo75WN1hcx9Aea8mluCz6JzaRhrZiFKnsYBna\nrNtQclAcNu4nlwfF/x79hiyfxMWggVtrmA7thCGBveS0XZPg0E4Y/r5IN9JX5U3Wq+2rr4v1\nU/EMb8HZQXdeZyMPvZkP7axbGRJJU+JurXpDh8lR3tfM/5Uj4f2GrLdxKm8YbFmLdFid92gG\nG6rOcPs/+okqkrR6L+GkT+UkScvDnG15EHQifVXZZDUSkPZel2uUa571h9O24tDb3W6voSuu\ndSvDe6RUv1VhvdNkk7t6luwJ1yqdqBsWF4MGbq0h/Hz5wcVhNZR7Vvyu6Vk9/F1eiXCS1r9/\nO5otrd5LOOlTOcVOq0uEil6GcIlQ2U+RN1lc53NS7hrl19ty2PmwvGivHbg+7HZ0+S4pvZuv\neXFYn9+5MBTXupUBkU7T4rvKVotA0ob26q2k1RuxIdutiYtBA7fWELrbuhOyxX/VR80qp9U4\nl/B/qHB+sd5au2FpY0pnpB342ipf7Z001b4WtiGcSk27IbeLTROh6Kjs9WvX9mYsWzGK1LBR\nq1kFkjZUnR3arwbJxYaU+0jtYtDAWST5EqGjvSStTpIcFUtPy77N2V53sKdczyM8qZ+qx/+X\n5dD0gfRnFL0S5AtymbWvCy6KK2yqMbbT8pocaQD8ZL/8M4ryC1WZdcW1bmVApP27vWrWgcQN\nXZYqH9U9KaEhxQ0Li0EDN5HiRZQPhz07B35xKiDSToNfnAqItNPgF6cCIu00+MUBIAAiAUAA\nRAKAAIgEAAEQCQACIBIABEAkAAiASAAQAJEAIIBKpPs27J8GYneCrr6qRDkaFIjEP+jqq0qU\no0GBSPyDrr6qRDkaFIjEP+jqq0qUo0GBSPyDrr6qRDkaFIjEP+jqq0qUo0GBSPyDrr6qRDka\nFIjEP+jqq0qUo0GBSPyDrr6qRDkaFIjEP+jqq0qUo0GBSPyDrr6qRDkaFIjEP+jqq0qUo0GB\nSPyDrr6qRDkaFIjEP+jqq0qUo0GBSPyDrr6qRDkaFIjEP+jqq0qUo0GBSPyDrr6qRDkaFIjE\nP+jqq0qUo0GBSPyDrr6qRDkaFIjEP+jqq0qUo0GBSPyDrr6qRDkaFIjEP+jqq0qUo0GBSPyD\nrr6qRDkaFIjEP+jqq0qUo0GBSPyDrr6qRDkaFIjEP+jqq0qUo0GBSFEG/eIXv7hAVC8gkgRE\nijHoF7/97W97mMS5qk6hGACRIgxaeORjEuOquoViAESKMChEUkIxACJFETTLxNe//RmIJIVi\nAESKIqggUnZ+fv7cZz7zGfSRulAMgEjLB82y7Dz/V9qUFR7lJj3326GjTgUiSUCkhYNmtTrn\nlU3nNdnwN6dEJQAiSUCkRYN25ihAJDEUAyDSgkH1Gj333HN+HnGo6rRQDIBIiwU17I2e+9a3\nvvUcRJJCMQAiLRXU4tG3vuUz1BB/VSeHYgBEWiioqXMEkTShGACRlglq8ggi6UIxwEWktHrM\nERfKb+0NMVOD8wlq9Ah9JF0oBjiIVAnTPKTy2wZ7Q8zU4JEGzQrkBWaPqlE7DH9LoRgwLFK6\nFc1Jq5dp81mDvSFmavA4g1bWZOoCOxBJDMUA50M7rUgd9oaYqcGjDJqpZrh45GdSLFUNFooB\nY0USekiPFIQp3WJkKuO2Uj40ZjTbcBSJqipgFrwHGzqDpNEG+/8oM/3PRRG08mbs/kEMWl08\nl6nbcPMIeyQxFAN890iaT0rsDTFTg08MqlWo2Z14BzXJ6CiSj0lM2nd8KAZApBKLQ15p3QXV\nbC27PxQFIulDMWCkSGsatXPKbl+RtFt0lsjTpLjblyAUA8aKtI7zSIM7orEi+Skz1aRo25cq\nFAN29soGd4d8sroOSuLRuCGO2YBIErt4rZ2vQx5Zvb0v/M0rgUle+s4LRJLYMZHGOeRu0pZs\nb+QRNKL2DRSKAbsj0hSHnJN6S7c3cg8aR/uGDMWAnRBpukOuOU2tEUQqQzFg7SIROeSU04Sd\nI4+o9yFSDKxYJEqHqpTOrJ3/EBadQ6T7EKlriJkavIXcIYesDuSRi0kQaXlWJ1Iwh+qsNqY1\nRAoXigFrEimwQ3Vam/I6XORBkyDS8qxEpFkcqtNan9ch94JDtYdIy8NfpBkdqvNal9hBSzBk\nEkRaHtYize6QMa3DlgIixQ8bkXp//L2EQ6a8Dl2SAZMg0vLELdLy0uhQ03qGEY6BRBvbvlOA\nSBKRixQ6Q8eRzV1IiBQ9EGkM2dxltJsEkZYHIo1gdo8gUvRApBFApHmDEuVoUCDSCDynHyYO\nSdi+U4BIEhBpDNn8BbSZBJGWByKNoZnrcc7yWUyCSMsDkcaQLVE6s0kQaXkg0iiymXdHTVDi\n9p0CRJKASKNY5lILk0kQaXkgEicMJkGk5YFIrNCbBJGWByKxAiLFCkTihdYkiLQ8EIkXEClS\nIBIzdCZBpOWBSNzQTBkBkZYHIvGjpxJEWh6IxBHFJIi0PBCJJbJJEGl5IBJLIFJsQCSWQKTY\ngEgsgUixAZFYApFiAyLxRDIJIi0PROIJRGoSuCQ9vKjeCJ8cpcJK+uWEQCSeQKQmgRtOtoow\nwhvTckIgEk8gUpPAZQZfHCTppfYT0zeogUhMEU2CSNvtQXKjfnMjTfaOql1V8XiWbqrlF/v5\nEWCznByIxBSIVCdwncFnyaZ6c1ge6B21Im2Sg3L5ZVos2YdIQEYwCSLVL6p/F9vTJK0/SZLD\nbSPYQb48waEdUOhMgkjbTqQ0OTjpPim0ql7uJZfKN0iBSHxpTYJI206kk/wQbq+xp/5Y2hFB\nJKDQmASRtvlR237z5mwvSU8hkghEspNNbN8pRCbSfnLUvTlqzelEwqEdMJNNa98pRCVScR5p\n2/SRTrdnwmBDs9Jhcpgvx2AD0JJNat8pRCNSzelWHP4uTyqlokgX5fD3Xr2cHIjEnGxK+04h\nKpH2Di+rN/nDYZqkuUf58Z0k0vZskyQHF/VyciASc3ZdpFiASNzJIFIMQCTuQKQogEjsySBS\nBECkSMjGA5EiACJFwviqQqQYgEiRAJEsoRgAkSJhdFUx2BAFECkSIJIlFAMgUiRAJEsoBkCk\nSIBIllBmPrJBlNtOQKRIgEiWUGYgkhMQyeGLEAkiDQGRXL4JkSDSABDJ6asQCSLZgUhOX4VI\nEMkORHL77tj2nQJEkoBIkQCRLKHMQCQnIJLjl8c28AQgkgREigSIZAllBiI5QSDSK6+8Mn0j\nMzCtqguYxE0k58nzxymxdpFe+f73v8/CJIhkCWVGNufxxx83iFTPJuRiAETSUHjEw6SJVZ3f\npAhFevxx2aR+nmvSHSK5AJGCEZ9Ijz+umKRP8/oYL9nWN1DKH5Nt+xIiaYFIweAqUnOM10jT\nTmzcfjBKgJWLtCt9pAVMYi+S8JyIn0IkLbsxageR7H0kUaSkvilmvVwQSfjAm9WLxAWIZAll\nxnHUTt0jbXUiCcd6/kCkSIBIllBmPrLRz/PEIhL6SKtgclXnNomZSK0rxj5SM+oAkVgDkSyh\nzLiKJFzZII5yV4Pe2/YD7JHYA5Esocy4ixQYiBQJ06s6s0kQSQIiRQJEsoQyA5GcgEg+Wxjb\nyuOASBIQKRIIqjqvSRBJAiLFArddEkSSgEixwG2XBJEkIFIsUNR1TpMgkgREigWSus5oEkSS\ngEixQFPX6qaycwCRJCBSLBDWdRaXIJIERIoF2rqGNwkiSUCkWIBI5lBmIJITEGn01sY2uTMQ\nSQIixQJpXdFHgkgiEGnsxrBHgkgCEGn85kK7BJEkIFIs0Nc1rEkQSQIixQJEMocy4y5S86fm\nasYTGQCRYiFAXYOaFKNITz75pEEkZfITegMgUiyEqGvIflKEIj35pGxSP881c5tApJURpq7h\nLr2LT6Qnn1RM0qZ5PXm+Opm+xw2UtECkWAhX1zAmsRVJnFA16b8fCUSKhYB1DWISW5HqB9Pz\nSCBSLISsawiT4hPJ0kdyESmZdGwHkWIBIplDmXEctXPdI40HIsUCRDKHMvORjX6eWw7pINJK\nCFrXACYxE8k8iT4GG9ZF2LrSD4LzEkm5siHpFmD4e2WEriv1CSVuIgUGIrkS+h6aM9SV1CSI\nJAGRHAl+V+c56kppEkSSgEhuFB6FNWmWuhKaBJEkhkVKc8RnZXGNvSHGtiBEIo9CdvEdRJIY\nFCmtH9rnNBUXN9gbYmwLxiXSyy//ScAIM9aVRCWIJOErUiq8E02yN8TYFoxHpNykl195OeQu\nada6EqgEkSTc+kgmkTrsDTG2BSMS6fxP/viPXw55cDdzXSerBJEkxook9JAeKQhTuqhEeuWV\noLukbO7KZhVhfrm54CVSbY/YR9qKr7Y7sEd65eVXCpOCDjcsUt2xuybskSS8Rep9UGFviLEt\nGJNI3y9MCjnacF7uJMIG0Ecd9dtAJAkXkbTDdFv5vb0hxrZgVCJ9/+WXg+6QKhZQCSIR4CBS\n2j2m8rJdEin8pQ0N86s06twSRJJwOCErPHUi7dp5pPPwF9u18OgqQSSJ4fNIaT1Ct9tXNszJ\nQrX2NAkiSeBau/hYqtZ+JkEkCYgUH5NrffXq1VFx/bJ77I/qjyXvIJITEGkMV1966aXwJkEk\nCYgUHxNrXXg0zqQ1i5SYJtFXGWkERIoPiKSGMiOb89hjjxlEStpHqoxXgEhhmHIyaDGRvEyK\nUKTHHpNN6qf59Gm3TEAkeorzm1NMWqyPxFykxx5TTNKmuWkS/WamVRzaRUMmPY3+/nhGjtqd\ne5nEViTDJPrNahApDoTrTrORSixY61WLJOxx1EO9RFhrnAAQiZZMfjOqAhBJDWVmZB9JLxZE\niga1wBCJAkveOY7a6UTqTVWcSKt6ApGIyITH3uIxm1oE5iIp9NPcMok+RIqC8ihO93d5I6qw\nZK3dTWIlknoeKWmm/1b2UOONgEg0ZMZTR/51gEhqKDOuIilXNtTj3fLoHYa/I8BSUp9KlNOR\nTC7MBFYrUmggEg22kg7WImshLNE4INJIIBINY0SKR58OiDQSiESDtaTFh60wEe1+NECkkUAk\nGuwlra2JVp8OiDQSiEQDn5LagUgjgUg08CmplbWekA0PRKKBT0mtQKSxQCQa+JTUCkQaC0Si\ngU9JrWQaqH9Ufyx5B5Gc4JOefEpqRVcNg04QSQIi0cCnpAryX9MaqwGRBoBINPApqYwyv4O5\nGj2TIJIERKIhhpKOONurzjgEkcYCkWiIoaRZ059x/wpEogIi0RBDSdur+dy/ApGogEg0RFLS\n0iWfslwd6xFEkoFINMRTUs+ekjRqB5FGA5Fo4FNSG6sVyTwfg8Oc+k6OQCQa+JTURluL559/\nXvmE7kf1x5J3sjmPPvqoTSRdslPdnAIi0cCnpDaaWjz/ne98RzaJhUiPPiqbpKS5NtchUlzw\nKakF0SPFJA4iPfqoYpKc5s3hXTV5UDeHfvtS/UiaZ2gIiEQDn5JasIjEYfjbLpI4T748P6Rp\npjvp/RAQiQY+JbWwZpES3fzeSf+laU7WISASDXxKaqEba+h5xEEkWx9JnGNVM4e+Mp2+/B4i\nzQifklroKtEbtWMh0sConTLR6rY7ftvKeyTlPUSaET4lNWOfUozsR/XHkncf2VDTXJ36W9NH\nMj4PAZFo4FNSM6sXSR1BSBLTot7zIBCJBj4lNbN+kUptujHuzqLE9gyRZoRPSc2sWCT3rB8r\nBESigU9JzVjnL6f7Uf2x5J23SOZrHGwfDQORaOBTUjM7IVJ72YLfR4NAJBr4lNTE1atPWD5d\nj0iBgEg08CmpgasvvfTUVeOnmI5rCIhEA5+S6in+Uvapl4wmaWaJhEgSEIkGPiXVA5EmApFo\n4FNSPeUsKOZjO4g0BESigU9JDRQmmUcbINIQEIkGPiU1cfXqVZ8pGyCSDESigU9JLUCk8UAk\nGviU1IbHHPoQSQYi0cCnpDYg0mgiF4kPc+Z7OEzVgEhDRC7SrFk0CUZFtQCRxgKRiGBUVBuG\nakCkISASFZzKaqZXi2r2Bog0BEQig1VhjSi1qOcT4i5SgL/kC7MZiHS+EpPkSjQz3DER6eGH\nH/YSKTF/5AtEooNXaU1IteAl0sMPyyYNpjlEihNmxdVTVeLpp58uXv/uv/2O6dguPpEeflgx\nSUlzeeaTpH7s3hazscqTgHsAkSjhVl4tRSWefvHFF5/OX2W/W8+5ugKR1Am2EmlxpZm8mg8Q\niRR2BdaQVR69+OKnitfNnKsMZhGyiyTODqksM80M6QNEIoVdgXVkpUi/8RtPS/VRTYpPJGsf\nqZOjOXgTREo0k31DpEXhV+I+tUgvlr0kYTHRj+qPJe9cR+0SeY8kHsNJ+yiIFAsMi6xSHttV\nHgnHduoMKDGKpCCnueEYDiJFCscyKxQmfepT5ZW4mXSHF9EkriL1Zvc2+AWRFoZjmRWy5ikr\nR+30JrETSTfxd/82mNLk3x5AJHJYFlpGuFGScvO+ziRWIo1IeYi0ODxLLdLV4PeefXZFIjl2\nftBHigSmxe4Q7oH57PPPSnfBbE3iKJLr7N4jJgGHSAFgWuwO8a7Mzz77e/JnU39Ufyx5N+LQ\nLgwQKQRcy93gcntziCQBkYLAtuA1dfl1IjUmQSQJiBQGviUvEXdJ6v3Na5MgkgRECgPfkpd0\now3P9zyCSDogUiAYF/186Hay2aQf1R9L3kEkJzhnI+eyu9yYGSJJQKRgrLfwEKkPRAoG68IP\nmwSRJCBSONZbeojUAyIFhLD4V6+a75QchqFdEkSSgEghISt/cT+9mEyCSCoQKSRU5S/v8Dq/\nSdYbj0EkCYgUFKIKLCOSRSWIpAKRwkJTg6VEKlTS3v4JIqlApLAQ7pKW8Kgd41BuqZY/sRIp\nsWY7hQQQKTBUJs0+aldG1flbqxSjSB//+MdtIumSPTEs9wUihYZxHUxHlNWeaeyP6o8l71SP\nJJOUNNfmOkRq+fznP0+eQpTwNcnaNZvRJEve9TwSTZLTvDm8CzOPPn+RPv/Nb34zapNWKtKM\nOyZL3rmLJHSTQsyjz16kwiOYFAiXMQ7x1u5jf+cBLHnnLFIi9pOEqVUtc0Tu1h6JgUicTfIb\n4wikkiXvXPtI0vyP2+bgTRBp6jz6EGkO+JrkSxCVLHnnM2onP4jHcNI+akdFir+PdL5LItUq\n0R7sWfLuIxtqmou7JV0fSfOxM/xFin7UrmCXTDoXL4QgMcmSd14iiaMI0rPBrx0TiQXrqYkn\nFDslS975iVT3hsSJ8qnm0YdIM7Giqngy/QDPkneOIo3LeogUISuqij8TVbLknbdIjp2fXewj\nMWFNdfFnkkqWvPPfIwWaRx8izcaqKuPPBJUseTfu0C4AzEXiMGLXsuMmjVfJkncQyYmh1ONw\nDqlj10VqLyeCSGYWEYnFVQ0CMKnC1yRL3kGkXRQJJtV47pQseQeRdlIkmNTgZZIl7yAShUjM\n+kjnEKkFIhlYRiReo3YFMKkCIhlYSCR+rLBKI9iBPtLRfpJsN2e+24FIrqyxTi2ufwi4epEu\n96rZIJJTz+1AJFfWWKcG1+n3PMe/LXkXq0gHyWHxZxh3k43ndiCSM6usVInzhLDrFylJun9e\nQCR3WNeqP31xB0Tq3vZFSnPEZ2VxDUTygHO1nrj6hPEziNRSH9odJgfNkrR+aJ/TVFzcAJE8\n4Fut7ImnXnrKbAr6SA2XaTUvUXrRLFFESoV3okkQyQe29XoiN+UpiyuOo3azi5RITyHobfrG\nXpLsHV7KC9P2WRapAyJ5wbViTzyVM/2+GKFE+tjHPmYTKaBHjtvuiyT0kB4psHw5mwDFjx8j\nXCtGdaOmbHzGyqgeSSZ1q3n/5bg3ThtPhefOIGm0YRd3K1Pg2ihEN2oKckL2Yx9TTJLTvD9J\npGbifPW9s33ieomISaSt7gOI5A3XViG6UVOIa+3sIolzRHaP6sT56nvn80AuIhk9gkgTYNss\nNCr5mEQhkjCJflJbopsUUvfsLVLB/uZiu73Y7GtsSfvLINJ4uDYL1V04Pe5hYclf1z6SPFtx\n0h7ZSRPnT5lIX1lvP6nG65LOJO3+SDiZBJFGwrNdAtwXetAkV5EGRu3qh/a4LlE+1L4fKVJ9\nSHfZHdqlaTVC1zyLy4UvQiR/WDZM1CIpqGmudJBCirRJqkO7ZF+/uhGINAKOLRNApEGTLHnn\nJVI3eCcN49meRw025Fz0rmxwBCKNgGXLUPWRBGYTSR69002cr74fPdhweVhc2XDjUruyBYg0\nBpZNQzQALjJgkiXvHEUaz1iRxgKRRoG2qbAP3VnyLqRIk4a/xwKRxoHGqS8Es5lkybugeySf\nifR7o3a6KxscgEjj2PnGaW6UGeEeyQuItDA73DrO84Bb8i5WkSouNjd8twORxrKTzeM3l74l\n7+IWaXuZ+JoEkUYza/sEGHIbwfr/QrZdjEO72ZizfQKcBBqB741dLHkXuUh3k/6V3nYg0njm\na6AQlyX4swv3R2rHGg49twORJjBbC0Uhkv99+yx5F7lIqa9HEGkSczURRAoGTsjGwKy7JHZH\ndhAJIrkyn0nLj9rthEjNaF2KwYZ52aFGWr9IqWnyEwcg0kR2p5XWL9KR4NGR53Yg0kR2p5XW\nL9J2xInYBog0lV1pJn+POIo0Gog0mcXaad4BiB0QqbxVH/pIi7FQQ807JD7Co5AiSXk+Mvd7\nG4JIy7JMQ813ktZxFrsJIj300EMeInVToWhk6K86vDEKIBIBi7TUbCKNUchPpIcekk1yyH2I\ntEqWaKr1iPTQQ4pJYppr5smvJgkSjr4EGRJpfYdDNJyQjYmQTWUcUpijjzTymI5OpHrv081O\nvBXmrZNmhuxeqVPdWcEJ2agI11YWXYKP2k2RiEgknRjSHHbls5z76uSrVnBCNi5CNdaC131P\n1Yiij+QokihGfQyY9D7RgxOykRGotZYTabJGFKN2QyIJt6HovpB0SmCwgR3zihTyqM71pi2E\nIikoae6xR1IX+4t0iD7SwgQ0qe8R+Y0lunv/Uhi0qEi9WyVZUdY4xGDD4oQyqb/zoT7e63Y/\nNHsiepHkefKTRBGpy/3mM+GtHeXzNDnbJBeXm+TU1aAaiEQH1ykcKN0hFik8/cGGG8nJ9jLZ\neG4HItHBVKRgHjEV6aQY+sah3ZKwnMIhnEccRdpP7l4ke9tTiLQoDKdwCOgRR5EKgzZFB+vA\nczsQiRRmbUY7trAGkbYne9vtgf/8kBCJFl5tFtQipiKNBSLRwqrRIBJEihZGrRbaI3YiiRd/\nY7Bhadg0W3CPIBJEmgCfZsMeaYtDu4hh024QaQuRYoZHw4Ud+YZIEGk6HFouvEYQCSJNhEHL\nzeERRIJIE4m96WY4rINIEImAuNtuHo0gEkSaTtSNB5E6IFLcRN14EKkDIkVOzK3HR6VU/R0A\nABqcSURBVKTEnO3SsvHT6EOk2JnSfGGnfpzLI3eRnnnmGYNI+unpejNtmXXwnLNhNBApGOPb\nL+xkxLN55CzSM8/IJolpLs1yImY/RNoVRrdf0EkhA/8t3xiRnnlGMUlM8/48QrTT6EOk+Bnb\ngAFFyi2y/eTEWPLOWaRm79POVNfsjYim0YdI0ZONNSmcSMXOiLtITfabZrYT13MQBSJFz2iR\nAvaRshhFsvaRGitc51rdek6jD5FiJzufcHDXG7UjGseLUyTLqF31YBFp6jT6ECl2svaBAKp9\nVKQiKShpbhVJ1UFdDJGYk0lPEyHrNbEUqT8t/oBIPtPoQ6TIiVKkOAcbBkWSh7+r18bBBs9p\n9CFS3GS9F1MgEiljJ1J4IFLcZJpXE6DoI9VnYiGSBESKmszwejSTR+3aCxogkgREihpykSYi\nXBcEkSQgUsxklncLIN0SFiJJQKSYyaxvZ0e6TBUiSUCkiOm13MJNKV3wDZEkIFLEQKQ6lBmI\nBAbRNNxibVle8Sdnd1h5pFBmIBIYRNdwCzVmVvyKEMkCRIoWfbvN35r1WJ36J7EQSQIiRUs0\nIhmyO4gz+lBmIBIYwNRsszcnRHIBIsWKsdlmbk/THCcQSQIiRYql1WZtUONcQRBJAiJFCkQS\nQpnZKZHACKzJPZdF5xDJFeyR9MRd6AHRSEMZszuAMaZQZhxFSqQnowGJ7skNiKQn+kLP5RIr\nkb72ta8ZRNLPqaVOD9SsZZqBy2ILRNLCoswzFNIyL3GEIn3ta7JJYppLUzOI2Z+oC0wr6xc5\nfOSFpRlYJKUCjzIHLWXZU7NlN7UullBmeh6JJolp3s0aNDj3d39lhwnAIZIWJmUOU8zsfMCh\nKrsJTRkKZcZZpN6Uxc3eqDf3t/BPnbsLeyQfhsbMYiJEQfPaZ/fXL1KT/b157Uwi2W2BSCp8\nLDoPIlJukNMNW+ITydpHaoSwi6QYJ0xyN2QLRFLhVVz60rre9ihCkSyjdtWDRaSu67Rt1xXt\nGJoAHCKp7HZx3W8fFqNICkqaO8z9bRYJfSRfmBWXtrwed+HjJ5LD3N9KdwiDDRNgVly68roM\n1YnZPVUPj1BmvESSh7+r11aRlOFv2wTgEEmFWXHpCux5T1hWIoUHIqkwK+45SYmzodOvmuye\nJodXKDMQKVqYFbdgepG9NYJIChBJhVlxCyhE8s/u8WJ4hzIDkaJl+M+B4mNqcUd4BJFkIJIW\nZmWeaP4IjSCSAkTSwqzM2ZQrm8bsju5DJAWIpIdXoSfd+RwiUQCR9PAqdFXakTulcR5BJBmI\npIdVobs7NnsXO+vN6e2c3eOkGBXKDESKHFaFFgrrq5L/idg2u8dJMSqUGYgUOawKnRnfDH5z\nQnaP/6p3KDMQKXJYFVq51axH2SESFRBJD6tCq4V1VmmCRxBJBiLpYVXofmEdiw+RyIBIelgV\nWnePTJcKTPEIIslAJD2sCq0t7KBKY4frmuye8mXPUGYgUuSwKrShsFaVxp9AarJ7ypc9Q5mB\nSJHDqtBj7klWzl43KbsnfdsvlBlnkZxnSU2SdhIuHyCSHk6Ftuti/s76RLp586ZBJGVmBvkD\ndbIgA1ZXIJIeToW230tJ34EqH6Zl96Rv+4Uyo3okmdRLc022Q6SwcCr0QFmNY3prE+nmTcUk\nbZpLswhtDfPoS6s6TKEPkUxwKvTw8Jy6fn21+LTsnvRtv1BmXEVqJdBNVdebR9+8qhGIpIdT\noYfLKqnUeLRLIpWpLu6JZJFEE5TBhkSnmAaIpIdToZ3OvXZ/atEtm5bdk77tF8qMYx+py/Zh\nkcT1XafQh0gmOBXarayZ9HQ+1aMYRbKM2rXZbhGp7grJcnSTHUOkMXAqtOuFdcofpK9QJIVe\nmvvtkdTFEGkEnArtXFZ5nrGdEqmVQh1kGBCpN+++iWGR0pzmpX5xgaUZOOVkC6dC+/z9kfh6\nYnZP+7pXKDOOIgnjdsKY9jbZ9ufRFwYbms+Gp9B3ECntDCrNqf0RFkOkZRlZ1IkecRMpNB4i\npcVTWr1LhQ9LLM3AKCc7GBUaIrEQqaTTRxapw9IMjHKyg9GkxRCJs0hCD+mRAsuX+aSkQDZ1\nGuD5GCvStMSJBF4ipeK/bughXXsfiYdLs04LKewmpm7AI5QZdiKpgwvtBzWWZuCQjD2aQjNw\nCSJxEaneESkD3s0nFZZmiD4TdUhzLs5WgytXrvh/CSIxESk1v9oJkc7VM5nBuPLCCy/4mwSR\neIiU6l6m6qGepRlWIFK5JLhLhUcjTBpVrMkeQSSZ4fNIwiHdrl/ZENgliGQKZYaPSI5YmmE9\nIp2HdWmcSAsd2UEkGYikxzqhSKgajeojQSSIFDH2QodyacyoHUSCSBEzWOhoTjBBJIgUMW5/\nvh1D1SASRIoY5786Xbx2EAkiRYx7oUOPiQ91mxYa/YZIMhBJj1ehA174MDiQt9QOCSLJQCQ9\n3oUO49LwqSWI5CaS8yT6I4FIesYUOoBLgUQi8ChKkX7xi18YRJKmZVAFgEjhGPvXcsQuuYj0\n2c9+1nOrFB7FKNIvfiGb1EtzTbZDpLCMLzStS0N9pNyjr371q5/99Kc/7bzJaTfq67KbZCtu\nocz0PBJN0qa5bRL99jN1Zv1hIJKeSYWmdKkZtTOM3mWFR1/9lX/6R3/kaBKRRtxEcppEv10m\nLnQEIumZWugQx3g6k0qRfuVf5yI5mJRNvG2snN1kWxoOZcZRpDLVByfRV58h0nQICk26XzL1\nlfxEIs1uyo0NhDLj2Efqst0qUiJM/w2RSKApNJlLJpHyzX+2OrJzEInUoxhFsozatdk+uEfa\nQiRKyApN45JFpGLU7tPzexSlSAq9NIdIs0NZaIoLH6584hOf0B3Z1c8Oo3bEHjETqdVCHWRA\nHyks1IWe6pJeJJ/p88mzm3qDllBmHEXauk2iLw+Ne+kBkfQEKPQUlwyHdtePj11j02c3+RbN\nocw4ixQaiKQnTKFllx48eOD6Rb1Ix9fu3XMyiV4jiKQAkfQEK3Tn0oOf/exnriZpRTq+d+32\n++/fcogZJLtDbNQQygxEipyQha5cKjzyM6m3Q7r3xq3bt98f3CcF8QgiyUAkPUELXd7qwk8k\nzSVCx7fefe9/vnH71vu37CqF8QgiyUAkPYFFyvnpa6+95iFSj+N7t9979713/9et9961Ht4F\n8ggiyUAkPeFFKo7tXvvpuDj57in36PjWzZvv3bx58/3bt94zD99BpFmASHpmEKkYtRt3P7Oi\nw/RqIdL7776Vi/Tnt99/77Zp+C6UR5GIFA0QSc8cIjWvfF0qh/Be/0Yu0u2/+d83b/74zb94\nOxfJYBJEmgeIpGdOkc49ryKqxsLzXdKtW//tzTff/fGbf/X+2xBpYSCSnplFKl+bXFIH7OqT\nSsfHt97+0Xe/+86P3/mbW2+YTs0G8wgiyUAkPTOJ1L+hWX/t3Jv/93cPlCUvvHp8fH791bf+\n9PUXXn/r1tv37l27ro0UMLvDbboXigEQSc8yIp1rXCo8+uu//jPZpCuv3rt37+1vfPe7r/9S\nsbvKjnOtNEUO6BFEkoFIehYT6Vx16Urh0V//mXzG6Tj36I1X38x5/Uq7EbHDVXa6QnoEkWQg\nkp4lRTqXXDKJdO2Nt26+8/rr71wRDJKnZQjqEUSSgUh6wookYF2pfHElF+nvf/5zRaRr9/7i\nL//yzXf++w/+tPtjCsq5TQaBSBIQSU8cha5cuvJ3f//ggSzS+fV7P/rer//qr/7gr/7LO2+9\nUY7YzWrRfYikAJH0RFPowqUHP//5z+XL8rLjt979tf/wB//pn//gzR9/79Xb9+5dn9ei+xBJ\nASLpianQ2U/VC8Wz43vf+LUv/M4Xfv/3/8UPf/j6q8dvXPtgtqxugEgSEElPVIV+UFwnLoiU\ne/Sb//gffOELv/Mf//DX33zz9Vev3/kQIi0MRNITV6Ef/Oy11wSPzo9/8ytf+eV/+Af//j//\nq3/2vR/++fX/e2d+jyCSDETSE1mhiwvFhUHuf/mVnF/+R//u3/zWf/3R/7hz5/8s4BFEkplB\nJJYmxVjmbsD8ypdykb70S//kt+7k/OSDvPU/+GBumSCSBETSE2+ZS5Fe+NKXvvTClet3cn/y\nw7rt/Q8+/HDubhJEkphDJI4mxVvksmTVXCitSIVHc5sEkSRmEWkbb1qaiLbEdcHKP624fudv\n/7YYsINIyzOPSPz2SdEWWCzYrev5DgkiRcFMIrEzKZtA2IJ1L4/vvZEf2v3kPvpIETCXSOxM\nmkBIlTJJpGt37hSnkDBqtzyziXQ/8P/UURFuwmPxTS5StSOaMaeVH3WmUAyYT6RCpZ1xKVRF\n5e0eX6sO6CDS8swp0v0dOsALVFFliof6gA4iLc/MIu3MAd7Uahpu+aKI1G/f+YBIEnOLtDM7\npWnVNNzyRe8RRIqA+UWCScOY7lQBkaJlAZF25fBuQi2dRMpM7TsPEEliCZF2Zac0vpYGkQw7\nJIgUAcuIBJMG0PeR5MsmIFJMLCTSbhzeTTm4M92oWZg33Na+4YFIEkuJtBM7pYkjd61LolXt\nNrOB9g0NRJJYTqQdMGlSDbujO+k4DyLFCUQKCM0IuDTyoPMIIkUARAoIRCIKxYAFRVq/ScFE\nOj6+7tC+YYFIEhApIPR9pNqje/ekqVUh0vIsKdLqTaIftSu3WNwc6Zr4B7EQaXkWFWntJtFX\nDyLFCkQKSDiRrn0IkeJiWZFWbhJ57YoNZtePb8k7JIhUpXJifhceiBSQECJl16/du/e2PNUJ\nRCpTeYdFWrdHQUS6fu+aMmYHkepU3mGRVn7tKnXliu0VIl278xPskRpO95MkPdxW6uT/9pPN\nRfXusFrerRGShUVa904piEjXrt1744M76CPVnCQlh41IB/mb9LJ4tV8v79YICUQKSAiRzo/v\n3b5+50OM2tXsJXe327N6b5T/21xuN5VW+asbSSquERKIZGfSPMTEdas3d/z2nQ8hUsfFyY1N\nJ9JZviDZK17VB3jCGiFZXKSou0lZJhbTXyjaqjVbu65Omb/TIm2qI7e2j7RVX3VrhAQi6aiE\nyTK5oOWD33ZoS1U/9abM32WRDpK9o5MLm0jdGiGBSBqyQhhFo7a4XhuiLVZbug9wHqlN4CKD\nLzuRLopDu01Pqcv1ixRfL6na9VAUd3LVdH9k3i/Zbot0ur0U+kib4t0NWaRmjZBEIFJsJhkV\naorrXt7JExdr/shcU7xdFukwkftIRYco3YoidWuEBCKpDHl030OliTWTp7eDSFoOcnlOhUO7\nTXLQjdeVj+0aIYFICoaukVpktzJTimQ+stttkeR0nvnCICEy0XbsDTHQUPGY5KaRc5khEk0o\n93TebZHiGblz9sitzJR9JPORHUTq0nnHRYpmp+QhkkuZKUftanEhkjWdIdLUlCPBy6M5RJK3\nVrgEkeIkGpFiOLpz7yE1hR7a4DwlhEjLE4tIMZjkqZFDqQNM2jC2fcmBSBLxiESecyQ5Oq3c\n9Huk0e1LDkSSgEgDSTpYblvBIRJNKAZAJHuOupTcUnS6WlXjdxApViBSw2iRrGWnqlZ1RslQ\nxh0W6SMbRLntRDQiLW5SGJGITKqvcYBIKhCpz8ImTRApvEkQyQBE0rCsSVNECt5PqkT66bT2\npQUiSUQk0tInkiaZFLqf9MCyQ4JIEEmAcR+pKb/p9CyJSQ+MYw0QCSIJLLxD8r08SFeBuhL9\nGbxIqmYuIUSCSCKM+0htTQuB2vp0FRo7M55b+SASRBJg3UVqaqrM4CUFCFY8iASRBFawQ+rV\nlMok+4EnRIJIIouaFLVIA4WDSBBJYsHxBhqP+jUVazS2doPjIBBpUCRxPq4gk3NFJdKCO6VQ\nIk03yWE4ESLVfPKTn9SLlIi5HuTP0SFSla1Tf2tzTZtJS4pZ+K1F6I2ZV4yLGp4IRfrkl7/8\n5U9CpIKFRKLySFvT1gfrdUSZtLajQZaowYlPpMIjySQlzRPalJeJTKRlTCLzaHAqTFP1ppUA\nInmJFGT+Yoh0TnFVg3NN9SpNjA+RvPdIax9sWOjYbj6R+ir5HcWNjUpPfCK59pEo8558g/aG\ncG6ypeYSCjjYoK1kM7IwXSLnqNREKJLjqB1l3pNv0N4QPq22hElz7pGkui4SlV1QS959ZENJ\n8506tFtEpNkGGwIBkXxFWv9gwxIi0XkEkcKEMuMoUntlQyK8JiU+kYKb9LnPfQ4isQpqyTtX\nkYKzeyJ97utf/7psEkSKPKgl7yCSmbAiFR4pJtGdRoJIYUKZYSRSmiM+K4tr7A3h1Wyzi4TB\nhsiDWrKTj0hp/dA+p6m4uMHeEH7tFtQkiMQvqCU92YqUCu9Ek+wN4ddu4XdJoTyCSGFCmeEj\nUomsTirtiyrsDeHXbqFHGzBqxyyoJTM5ilQe1DUiCT2kRwooyzT3maSMsvBgN3ESqbSnU6hZ\nmnLsI/UsIhyzwx4pUCgGuIvUvFA/qLA3hF+7zSvShB94ck05R4VIEi4iaYfptvJ7e0P4tRtE\n4hAVIkk4iJR2j6m8DCJNrCnnqBBJwuGErPDUicT2PFJQjyBSmFAMGD6PlNYjdLNd2QCRWESF\nSBIRXms3p0nEHkGkMKEYsNsiUXsEkcKEYkCUIs1lErlHEClMKAbsskj0HkGkMKEYEKVIbD2C\nSGFCMSBGkdge2EGkQKEYEKFI4T0qZ88L4RFEChOKAbsoUjWhPURiE5QoR4MSn0izeBQKiBQk\nFAN2T6SQHkGkMKEYEJ1IrHdIEClMKAbsnEhBPYJIYUIxIDaReHsEkcKEYkBkIjH3CCKFCcWA\nuEQK7FGYIW8RiBQkFAOiEimkRwQ3xnMAIgUJxYBdEWkGiQogUpBQDIhJpKA7JILf0wGIFCQU\nA3ZEpJk8gkhhQjEgIpFW4BFEChOKAfGItAaPIFKYUAyIRqRVeASRwoRiAEQiBSIFCcWAWERa\nwYhdAUQKEooBOyDSjB5BpDChGACRSIFIQUIxACKRApGChGLA+kWa0yOIFCYUA1Yv0qweQaQw\noRgAkUiBSEFCMSAWkUKZNK9HEClMKAZEI1IYk2b2CCKFCcWAdYs0t0cQKUwoBqxapNk9gkhh\nQjEgGpHoPZrjT8tVIFKQUAxYq0izzNDQByIFCcWAWEQi8SjrWH92LR0VIkmsSSR5H7T27Fo6\nKkSSgEikQKQgoRgAkUiBSEFCMWBNIskmrT27lo4KkSRiEYl+l7T27Fo6KkSSgEikQKQgoRgA\nkUiBSEFCMSBakTJ/szDYsNKgRDkalGhEEk2qzqj6qqRcyrD27Fo6KkSSiEikxpvu4h6DSs3V\nC1aPVp9dS0eFSBLxiFTvk+Rr5ERfhMt/ug/bFXqX1q09u5aOCpEkIhNJd6mp6o9+hbFBaYFI\nQUIxICaR7pNesb327Fo6KkSSiEuk6W3uH5QUiBQkFAOiEomU3Qm6+qoS5WhQIBL/oKuvKlGO\nBgUi8Q+6+qoS5WhQIBL/oKuvKlGOBgUi8Q+6+qoS5WhQIBL/oKuvKlGOBgUi8Q+6+qoS5WhQ\nIBL/oKuvKlGOBgUi8Q+6+qoS5WhQIBL/oKuvKlGOBgUi8Q+6+qoS5WhQIBL/oKuvKlGOBgUi\n8Q+6+qoS5WhQIBL/oKuvKlGOBgUi8Q+6+qoS5WhQIBL/oKuvKlGOBgUi8Q+6+qoS5WhQIBL/\noKuvKlGOBgUi8Q+6+qoS5WhQIBL/oKuvKlGOBgUi8Q+6+qoS5WhQIBL/oKuvKlGOBgUi8Q+6\n+qoS5WhQIBL/oKuvKlGOBgUi8Q+6+qoS5WhQIBL/oKuvKlGOBgUi8Q+6+qoS5WhQIBL/oKuv\nKlGOBoVKJCuPzBFkd4PuUlWjBSLxD7pLVY0WiMQ/6C5VNVogEv+gu1TVaIFI/IPuUlWjZRaR\nAFg7EAkAAiASAARAJAAIgEgAEBBIpDRHfG4XSgvWFLR4OVvQHWpfLoQRKa0f2ueyoQM39qJB\ni8cy3BxBd6h92TCLSKnwLhyLBq0jzhN0h9qXDQH7SIu0eXeYNXPQdDujSE3UVI0/W1CIpACR\niIIuJFJ5fAWRliecSO0BdHc4Hb5fKvb7Zw3a/psraB1tu0D7pmk6b1U5MI9I3ZLApMrzXEHl\nis65G5y/fWcPyoJgIqXyg/QiGH2PZgqq/Lc8138Zi7QvRNIRSqS0e5yvzTUezRd03j3Sgu0L\nkXSEOiErPM3W5jqPZgu6nVWkJdsXIukIdB4p7TqkqbI0WL+02bhylDVP0PK1sjTsNQaLtS+u\nbNCBa+0AIAAiAUAARAKAAIgEAAEQCQACIBIABEAkAAiASAAQAJEAIAAiUZMMNmmxxpHuUoDh\nr4JYwU9HjZtI2rUgEl/w01HjZgNEWhn46agpbbg4SJKDi/LdxX6SHhYfXGySvZN6b5Qk3W6p\n+sIm2a9eXRZfvVyu/GAUEImawobLtDAlvSzelS8Pm2UGkcoP98v35Wp7C1cC+AKRqClsOEw2\n2+2m0CdJNpfboyTdbm/kyy43jT+lM51Ih+2HN4pvHSZHC9cCeAKRqCls2Evyw7qLYr+SJBfy\nMr1I3Yd71bL9JasA/IFI1MjHbKZXikjiq+YIEHACPxg1EGknwQ9GjXpopyzriXShO7QD3MDP\nRo062KAsE0VKk7vtEMPmsvrwsPjW3WJdwAmIRI06/C0v60RKC2lybsjD3/VqZ0tXA/gBkahR\nT8h2yzbJ3t1GpHJEfHuYJjfqD/ebE7LFVzeny5UfjAIizUuy8/NWrRSINBdJcloczB0sXQ4Q\nBIg0F4dVF+li6XKAIECk2Tjaq/tNYIVAJAAIgEgAEACRACAAIgFAAEQCgACIBAABEAkAAiAS\nAARAJAAIgEgAEACRACAAIgFAAEQCgACIBAABEAkAAiASAARAJAAIgEgAEACRACAAIgFAAEQC\ngACIBAABEAkAAiASAARAJAAIgEgAEACRACAAIgFAAEQCgACIBAABEAkAAiASAARAJAAIgEgA\nEACRACAAIgFAAEQCgACIBAABEAkAAiASAARAJAAIgEgAEACRACAAIgFAAEQCgACIBAABEAkA\nAiASAARAJAAIgEgAEACRACAAIgFAAEQCgACIBAABEAkAAiASAARAJAAIgEgAEACRACAAIgFA\nAEQCgACIBAAB/x/ZjNMe3+M2BgAAAABJRU5ErkJggg==",
      "text/plain": [
       "plot without title"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 420,
       "width": 420
      }
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "set.seed(3)\n",
    "ggplot(data = BW_states) + geom_sf() +\n",
    "geom_jitter(data = choppies_map, \n",
    "           mapping=aes(x = longitude, y = latitude, \n",
    "               colour = District, alpha = 0.5), \n",
    "            width = 0.1, height = 0.15) +\n",
    "ggtitle(\"Location of Choppies Stores in Botswana\" ) +\n",
    "theme(plot.title = element_text(hjust = 0.5)) +\n",
    "scale_color_viridis_d()\n",
    "\n",
    "#ggplotly(g)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 237,
   "id": "d38ac5e9",
   "metadata": {},
   "outputs": [],
   "source": [
    "library(gganimate)\n",
    "library(gifski)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 238,
   "id": "3c4c9dd9",
   "metadata": {},
   "outputs": [],
   "source": [
    "set.seed(3)\n",
    "map_with_data <- ggplot(data = BW_states) + geom_sf() +\n",
    "geom_jitter(data = choppies_map, \n",
    "           mapping=aes(x = longitude, y = latitude, \n",
    "               colour = District, alpha = 0.5, group = year), \n",
    "            width = 0.1, height = 0.15) +\n",
    "scale_color_viridis_d()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "id": "70d07fee",
   "metadata": {},
   "outputs": [],
   "source": [
    "num_years <- max(choppies_map$year) - min(choppies_map$year) + 1\n",
    "\n",
    "map_with_animation <- map_with_data +\n",
    "    transition_time(year) +\n",
    "ggtitle('The growth of choppies: {frame_time}')\n",
    "\n",
    "map_with_animation <- map_with_animation +\n",
    "shadow_mark()\n",
    "\n",
    "\n",
    "animate(map_with_animation, nframes = num_years, fps = 1)\n",
    "anim_save('choppies_growth.gif')"
   ]
  }
 ],
 "metadata": {
  "celltoolbar": "Raw Cell Format",
  "kernelspec": {
   "display_name": "R",
   "language": "R",
   "name": "ir"
  },
  "language_info": {
   "codemirror_mode": "r",
   "file_extension": ".r",
   "mimetype": "text/x-r-source",
   "name": "R",
   "pygments_lexer": "r",
   "version": "4.1.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}