"
]
},
"metadata": {
"tags": []
}
}
]
},
{
"cell_type": "markdown",
"metadata": {
"id": "iu66LbfuMGKR",
"colab_type": "text"
},
"source": [
"#K-means\n",
"###Before applying k-means, we should point out that it is good practise to disregard highly correlated variables when clustering. This is because their importance will be hightened if we include all of them.\n",
"###The same applies to categorical data, hence we will omit Region and Channel varaiables. \n",
"###The correlations plot below, shows that Milk is highly correlated with both Grocery and Detergents_Paper. We will exclude the last two from our clusterin data. "
]
},
{
"cell_type": "code",
"metadata": {
"id": "x1OJ9K2AN9mV",
"colab_type": "code",
"outputId": "439d6ffb-5e4a-41ae-f9e8-173626239436",
"colab": {
"base_uri": "https://localhost:8080/",
"height": 367
}
},
"source": [
"corr = dat.iloc[:,2:8].corr()\n",
"#corr.style.background_gradient(cmap='coolwarm').set_precision(2)\n",
"sns.heatmap(corr, annot = True)"
],
"execution_count": 0,
"outputs": [
{
"output_type": "execute_result",
"data": {
"text/plain": [
""
]
},
"metadata": {
"tags": []
},
"execution_count": 13
},
{
"output_type": "display_data",
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAbMAAAFMCAYAAABBKVJLAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzs3Xd4VNXWwOHfSuglIYF0QHqX3kQg\n1FCkCiIgoIIfiNhAQUUFRAUEL16vDVARRAURpah0KQLSIx1EmkI6BAglkSSzvz9mCBOSwEQymYyu\n93nOw5y995mzJgmzZpc5R4wxKKWUUu7Mw9UBKKWUUndKk5lSSim3p8lMKaWU29NkppRSyu1pMlNK\nKeX2NJkppZRye5rMlFJK5TgRmS0isSJyIIt6EZH/icgxEdknIvXt6h4Wkd9t28OOnE+TmVJKKWeY\nA3S8RX0noLJtGwp8BCAivsB4oAnQGBgvIj63O5kmM6WUUjnOGPMzEH+LJt2Bz43VNqCEiAQBHYA1\nxph4Y8x5YA23ToqAJjOllFKuEQKctts/YyvLqvyW8uVoaCrHJZ894VbXG6tarZerQ8i2S8lXXR1C\ntnxSuIGrQ8i2XvEbXR1Cto0LauXqELJt3B9fyp0cn533mwJ+FYdhHR68bpYxZtadnP9OaDJTSill\nZUl1uKktcd1J8ooAytjtl7aVRQCtbirfcLsn02FGpZRSVsbi+HbnlgGDbKsamwIXjTFRwCogTER8\nbAs/wmxlt6Q9M6WUUlaWHElSAIjIfKw9rFIicgbrCsX8AMaYGcByoDNwDLgKPGqrixeR14Gdtqea\naIy51UISQJOZUkopG5OaknPPZUy/29QbYEQWdbOB2dk5nyYzpZRSVjkzfOgSmsyUUkpZZWMBSF6j\nyUwppZSV9syUUkq5vRxcAJLbNJkppZQCwGjPTCmllNvLwdWMuU2TmVJKKStdAKKUUsrtufEwo17O\n6l/mlUnTaXlfX3oMeNzVoaQzbtIY1u1YyvKNX1OzdrVM2zw3dgSb965g/6kt6cr7P9KbFT8v5If1\nC1j4w2wqVamQGyEz6a2X2fHrajZsWUbtOjUybTP21WfZc3ADpyLC05WHlA5i8fefs27TYjZsWUa7\n9i1zI+Q0/q1r03bz27TdOp3KT3bNUF9uUFtar59Cq7WTaL50PMWr3Pai5U7xzvSJHDm0mfDda6hX\nt1aG+sKFC7Fsyecc2L+RvXvWMenNlzK06dmzMynXImhQv7bT460YWpsn1k3jyY3/4d7hGX+uZRtX\n4/9+fINXjn9O9c6N09V5BZfkoXkvMvynqQxfOxXv0qWcHm8GFovjWx6jySwLIpIqInvstnI58Jyn\nRMQFf6E39OjcnhnT33BlCBm0atecchXK0qZxd8aOeoPXp43NtN1Pq36mZ9jADOXLFq2gU8s+dGnd\nl5nvz+Xl10c5O2TatW9JhYrlaFwvjOeeeZWp0ydk2m7VivV0aPNAhvJRo4ezdMkK2rToydDBI3nr\nP+OdHLEdD6H25EfZ2n8q61qOJqRnswzJ6sx3v7C+9YtsaDeWYx98T80JA3IvPptOHdtQuVJ5qtVo\nzvDhL/DB+5MzbTf9nRnUujuUho060OyeRnTs0Dqtrlixojz95BC2bw/P9NicJB5Cp9cf4auHp/Jh\nuzHU7HYPpSqn/7lejDzL0udmsn/pLxmO7zH9cbbO/IGP2o7hk26vcuVsgtNjziB3r82YozSZZS3R\nGFPXbjtlXykibjlE27Du3Xh7FXd1GOm06xTK4oU/ALBn9368vIvjF5Ax5+/ZvZ+4mLMZyi9fvpL2\nuEiRwphcuGlOx/va8vX8JQDs3rUXb28vAgL8MrTbvWsvMTFxGcqNMRQvXgwAL6/iREfHOjdgOz71\nKnHlZAxX/4zFJKcSsWQrgR3S31Ym5XJi2mPPIgWB3L8TUdeuHZj35SIAtu8Ix7uEN4GB/unaJCYm\nsWGjNTEkJycT/ut+QkKC0upfmzCGaW9/SFJSktPjDalbkfOnYrhwOg5LcioHv99G1fbpf64Xz5wl\n9shpjCX9z7NU5RA88nlyYvMB62u5+hcpSdecHnMGbtwzc8s3ZFcRkUeA+4FigCcQKiKjgT5AQWCx\nMWa8iBQFFmK9dYEn8Lox5mvb0zwlIl2xXnDzAWPMkVx+GXlOYJA/URHRafvRkTEEBvlnmriyMnBw\nHwYPH0D+AvkZ0HOYM8JMJygogEi7mCMjowkMDsg0cWVm2uT3Wbj4Ux4bOoAiRQvTq/ujzgo1g0JB\nPiRGnkvbT4yKx6d+pQztyj/anorDOuORPx9ber+Za/FdFxIcyJnTkWn7EWeiCAkOzDLxe3t70eW+\n9rz3/qcA1KtbizJlgli+4ieeG+X8YfXigb5cjLrxc02IiiekXkWHji1ZPpCkhKs8MPNZSpTx4+Tm\nA/w0ZUGGpOdsxpKcq+fLSdozy1phuyHGxXbl9YHexphQEQkDKgONgbpAAxFpifUW35HGmDrGmFrA\nSrvjzxpj6gMfAc/nzkv555s3eyGtG3Vj6sR3GTHqMVeHc1s9e9/Hgq8WU6dGKP16D+XDmVMRuaP7\nKua4k5+tYW3TkRx8Yz5VRvZwdTi35OnpyZfzPuD9D2Zz8uSfiAhvTxvP6DETXR2aQzzyeVK2UVXW\nvPEln3R9FZ+y/tR5IHfnUQG37plpMsua/TBjT7vyNXa3Iwizbb8C4UA1rMltP9BeRN4SkRbGmIt2\nx39n+3c3UC6zE4vIUBHZJSK7Pvl8fg6+pLxj4OA+/LB+AT+sX0BszFmCQgLT6gKDA4iO+nvDbt9/\nt4qwzq1yKMr0Bj/Wn/WblrB+0xJiYuIItos5ODiQ6MgYh5/roYG9Wbp4BQC7du6hYKGClCzpk+Mx\nZyYp6jyFg0um7RcO8iUpKus7bEQs2UpQx4a5ERrDH3+YXTtXs2vnaqKiYyhdJjitLqR0EBGR0Zke\nN+Ojqfx+7CT/e+8TAIoXL0bNmtX4ac0ijh3dRpMm9Vn83WdOXQRyKToe76AbP1evIF8uRZ936NiE\nqHhiDv3BhdNxmFQLR1btJqhWeWeFmjWdM/tXuWL3WIDJdkmvkjHmU2PMUaw9uP3AGyIyzu6Yv2z/\nppLFMK8xZpYxpqExpuFjg255FwW3NW/2Qrq07kuX1n1Zs3w9Pft0AaBug7u5lHA5W0OM5SqUTXvc\nOqwFp06czvF4AWZ/8hWtW/SgdYserPhhLQ/2s/ZWGjSsQ0LCJYeHGME6ZNYy9B4AKlepQKGCBTl7\n9ra3bMoRF/Ycp2iFQIqU9UPyexLS4x6iV+9O16Zo+RuJOqBdPa6czDyJ5LSPZsylYaMwGjYKY9my\nVQx8qDcATRrXJ+FiQqZDjBNfG4O3d3FGPXdjEU1CwiUCg++mUpWmVKrSlO3bw+l5/6PsDt/ntNgj\n9p7At3wgJcr44ZHfk5pdm3J0ze7bHwhE7j1OQa8iFPG1zmeXb1aDuN8jnBZrliypjm95jM6Z3ZlV\nwOsi8qUx5rKIhADJWH+u8caYL0TkApBnxr1Gj5/Czl/3ceFCAm17DOCJIQPp1bWDS2Nav2Yzrdo1\nZ/3OZSQlJjHm6QlpdT+sX0CX1n0BeGH8M3Tr1YnCRQqxZd9KFn6xmHenzmTgkAe5N7QJKckpXLyY\nwPMjXnV6zGtWb6RdWCg79qwh8WoiT4+4sQJz/aYltG5hTXTjJo6mV+8uFC5SmL2HNvLF598wbcr7\njHt5Cu/87w2GPfEIGMNTT7zo9JivM6kW9o2dwz3zX0Q8Pfhz/gYu/RZBtTG9ubDnBNGrwyk/OAy/\nlrUwySlcu3iF8Kc/yrX4rlu+4ic6dmzDb4e3cDUxkcceu7FKddfO1TRsFEZISBBjX3qGw0d+Z+cO\n682IP/zwM2Z/lvsjGibVwopxc3jo8xcQTw/2LNxI3O8RtBrVi8h9Jzm6Npzg2hXoM2skhbyLUKVd\nPUJH9mJG+xcwFsPaN79i4FdjQYSo/ScJn78u119DXuxxOUpMbiz9ckMictkYU+ymskeAhsaYJ+3K\nnuFGsroMDAAqAdMAC9bkNtwYs0tETtmOPysiDYG3jTGtbhVH8tkTbvULqlqtl6tDyLZLyVddHUK2\nfFK4we0b5TG94je6OoRsGxfUytUhZNu4P768o4nXpK3zHX6/KXRPvzw1yas9syzcnMhsZXOAOTeV\nvQu8e1PT41h7bTcfX87u8S6stxRXSqm8IQ8u7HCUJjOllFJWmsyUUkq5O2Py3sIOR2kyU0opZaU9\nM6WUUm7PjVczajJTSillpTfnVEop5fZ0mFEppZTb02FGpZRSbk97ZkoppdyeJjOllFJuT4cZlVJK\nuT1dzaiUUsrt6TCjchZ3uwr9b0e+dXUI2Zay/GNXh5AtYS/tcHUI2ZbwpmtvM/R3HPyvYzfW/EfR\nYUallFJuT3tmSiml3J4mM6WUUm4vVa+ar5RSyt25cc/Mw9UBKKWUyiOMxfHNASLSUUR+E5FjIvJi\nJvXviMge23ZURC7Y1aXa1S273bm0Z6aUUsoqB3tmIuIJfAC0B84AO0VkmTHm0PU2xpiRdu2fAurZ\nPUWiMaauo+fTnplSSikrYxzfbq8xcMwYc8IYcw1YAHS/Rft+wPy/G7omM6WUUlYWi+Pb7YUAp+32\nz9jKMhCRu4DywDq74kIisktEtolIj9udTIcZlVJKWWXjclYiMhQYalc0yxgz62+euS+wyBhjv5zy\nLmNMhIhUANaJyH5jzPGsnkCTmVJKKQCMxaHhQ2tba+K6VfKKAMrY7Ze2lWWmLzDipuePsP17QkQ2\nYJ1PyzKZ6TCjUkopq5wdZtwJVBaR8iJSAGvCyrAqUUSqAT7AVrsyHxEpaHtcCrgXOHTzsfa0Z6aU\nUsoqB6/NaIxJEZEngVWAJzDbGHNQRCYCu4wx1xNbX2CBMelWlVQHZoqIBWuna4r9KsjMaDJTSill\nlY1hRkcYY5YDy28qG3fT/oRMjvsFuDs759JkppRSyipF72em8phxk8bQqt29JCUmMfqp8RzcdyRD\nm+fGjqDng13w9vbi7nL3ppX3f6Q3Awf3ITXVwtUrVxk76g2OHT2Rm+Gn88qk6fy8ZQe+PiVY8sUM\nl8Vhb8vxaKau3ofFGHrWLcfgZlXT1U9bs4+dp+IASEpJJf7KX2x+visAURev8tqP4cQkJCIC7z3Y\njJASRZ0W6zMTR9C0TRP+SvyLSSOncvTA7xnaVLm7MmPfGUPBQgXZtm477477IF39g8Me4Mlxj9Ol\nVk8unk+gfc+2PPREXxC4eiWR/7z0X44fct7fiEe5WhRo2x9ESNm3iZQd6T7sk791XzzLVrPu5CuA\nFPEi8b0nnRZPVrxa1aPsxCGIhwdx89cS/cF36er9BnbA/+FOYLGQeiWJU2M+JOn3M0g+T8q9PYIi\ntSog+Tw5t2g9Ue9/l8VZnMix74/lSZrM/iYRMcCXxpgBtv18QBSw3RjTRUS6ATWMMVNEZAJw2Rjz\ntm1VzvPGmF3Oiq1Vu+aUq1CWNo27U7fB3bw+bSz3dxiUod1Pq37m80+/Zt32penKly1awVdzFgHQ\ntmMoL78+ikcfzP03hut6dG5P/17dGPv62y6LwV6qxTB55V5m9G9OgFdhHpq9ntDKQVT080prM7p9\n7bTH83ce50h02lV6eGXZLh67tyr3VAjg6rUURJwXa9M2jSldvjT9mg+iRv3qPDf5GYZ1zfi7fG7y\ns0wdM51D4YeZNm8yTVo3Zvt6633T/IP9aNyyAdFnYtLaR52O4sneI7l88TJNWjdmzFujMn3eHCFC\ngfYD+GvhfzCX4ik0cBypx/dgzkWmNUlev4Bk2+N89driEVDWObHciocHd705lKP9JnAt6hw1lk/l\nwuodJP1+Jq3JucU/EzdvFQAl2jei7PhHOTrgdXy6NEMK5ONgu2fxKFSAWhve49ySTVw7E5e7r0Gv\nzfivdAWoJSKFbfvtsVt2aoxZZoyZ4orA2nUKZfHCHwDYs3s/Xt7F8QsolaHdnt37iYs5m6H88uUr\naY+LFCns8g9rDevejbdXcdcGYedAZDxlfItS2qco+T096FCjNBuORmXZfsXB03SsWRqA43EJpFoM\n91QIAKBIgXwUzu+8z5TNO9zLykWrATgUfphi3sUo6e+brk1Jf1+KFi/CofDDAKxctJoWHW/01J+a\n8AQfvjkL+/n5A7sOcfniZQAOhh/CL8jPaa/BI6gC5nws5mIcWFJJObIdz0pZX+XIs3oTUg5vd1o8\nWSlarzJ/nYrirz9jMMkpxC/djE+HxunaWC4npj32KFLwxv8tY/AsUgg8PZDCBTHJKaTatc01FuP4\nlsdoz+zOLAfuAxZx41IsLQBE5BGgoTEm04+rIuIBzAbOGGNeycmgAoP8iYqITtuPjowhMMg/08SV\nlYGD+zB4+ADyF8jPgJ7DcjI8txd7KYnA4oXT9gO8CrM/Ij7TtpEXrxJ54QqNy/kD8Ef8ZYoXys+o\nRduIuHCFJuX9eaZ1LTw9nNM98wssRWzkjU/3cVFxlAosxbnYG/GWCixFXJR9m7P4BVo//DQPa0Zc\n1NlbDiF26dsprRfnDFKsBObSjXjNpfN4BFXIvK1XSTy8S2H587DT4slKgUBfrkXe+D92LeocRetV\nydDO/+FOBAzthkeBfBzpY10Lcf7HrZTo0Ji6v87Go3BBTk+YTeqFy7kWexo3vtO09szuzAKgr4gU\nAmoDjn4czAd8CfyeWSITkaG2y7jsSkhyPAHlpHmzF9K6UTemTnyXEaMec0kM/wSrDp6mXfWQtGSV\najH8evoso9rezZeDWxNx/grL9v3h4igzV7BQQQY+1Z9P356TZZt6zepyX79OfDTp49wL7BY8qzUm\n5eiuPD33Ezt3BfvvHc7pNz8n+JkHAChatzKkWthbfwj7mj5OwLDuFCwbkPvBuXHPTJPZHTDG7APK\nYe2VLb9163RmAgeMMW9m8byzjDENjTENvQplHB7MzMDBffhh/QJ+WL+A2JizBIUEptUFBgcQHRWb\njfBu+P67VYR1bvW3jv2n8i9eiOhLN4aAYhIS8bfrqdlbeegMHWveuAhCgFdhqgaUoLRPUfJ5eNC6\najCH7ebTckLPh7sze/VMZq+eybmYc/gH3xgC9Avy42x0+g9IZ6PPphsm9AsqRVz0WULKBRNUNpDP\n1sxi4bYv8Qvy49NVM/D18wGgYvUKvDDtOcYOHkfC+YQcfQ32zOULSPEbQ6NS3Adz+XymbfNVa0yq\nC4YYAa5Fx1Mg+Mb/1wJBJUmOPpdl+/ilmylhG4b07dmSixt+xaSkknLuIpd3HqFInYpOj/lmJiXV\n4S2v0WR255YBb5O9qz3/ArS29ehyxLzZC+nSui9dWvdlzfL19OzTBYC6De7mUsLlbA0xlqtwY/K8\ndVgLTp04fYvW/z41g334M/4yEReukJxqYdWhM4RWCcrQ7uTZSyQkJVMn5MYbcc0gHy4lJRN/5S8A\ndpyKpUKpnJ0PXDx3KYPDhjE4bBibVm2hY+8wAGrUr87lhCvphhgBzsXGc+XSVWrUrw5Ax95hbF61\nhRNHTtKtTm/6NH2IPk0fIi4qjiEdHic+7jz+wf688fEE3nhmMqdPnMkQQ06yRJ1EfAIQ71Lg4Um+\nak1IPbYnQzvxDYRCRbFEZnnFI6e6sud3CpYPokAZfyR/Pny7N+f86p3p2hQsf+PvxLtdA/46aZ1r\nvRYRR/F7rV+r8ihckGL1q5B0LKsrPzlRDt/PLDfpnNmdmw1cMMbsF5FWDh7zKdASWCgi9xtjcvTL\nHevXbKZVu+as37mMpMQkxjw9Ia3uh/UL6NK6LwAvjH+Gbr06UbhIIbbsW8nCLxbz7tSZDBzyIPeG\nNiElOYWLFxN4fsSrORleto0eP4Wdv+7jwoUE2vYYwBNDBtKraweXxZPPw4MXO9Rl+PwtWCyG7nXu\nopKfFx9uPESNoBK0qhIMwMpDp+lYozRit1zR00MY2bYWw77ahDFQPagEveqVd1qsW3/aTtM2TViw\nZR5JiUlMHjUtrW726pkMDrPOh04f++6Npfnrd7Bt3a3nwB4dORBvHy9GTXoGgNSUVP6v8xPOeRHG\nwrW1X1Cw9yjw8CBl/2bMuUjy39sDS/QpUo9bE1u+ak1IPeK8ubvbSrXw5ysfU/Wr8eDhwdmvfyLp\n6GmCn+/H1b3HuLBmJwGPdMarRW1rD+ziZU48+z8AYuesoPw7T1Fr3bsgwtmv15F42AXDz3lw+NBR\nYvLw2HJeJiKXjTHFbiprhXXZfRf7BSBZLc0XkdeAKsBDxmT+UadCqXpu9Qv67ci3rg4h21KW5435\nHkeFveTCN+y/adWzd7k6hGw7+N/MhzLzskYRi+9oJdGVCf0cfr8pOmG+E79Ukn3aM/ubbk5ktrIN\nwAbb4znAHNvjCXZtWtk9Hu/MGJVSKlvcuGemyUwppZRVHpwLc5QmM6WUUgB5cpWiozSZKaWUstJh\nRqWUUm5Pk5lSSim3p3NmSiml3J72zJRSSrk7k6I9M6WUUu7Oje9npslMKaWUlQ4zKqWUcnuazJRS\nSrk7d75WryYzpZRSVtozU85yKfmqq0PIFne7Aj1Avs7/5+oQsuXks9+7OoRs8wh93NUhZNvh/y51\ndQjZ1ugOj9fVjEoppdyf9syUUkq5PfftmGkyU0opZWW0Z6aUUsrtaTJTSinl9nSYUSmllLszKdoz\nU0op5eZ0zkwppZT702FGpZRS7s6N782pyUwppZSNJjOllFLuzqS4OoK/z8PVASillMobjMXxzREi\n0lFEfhORYyLyYib1j4hInIjssW2P2dU9LCK/27aHb3cu7ZkppZQCcnbOTEQ8gQ+A9sAZYKeILDPG\nHLqp6dfGmCdvOtYXGA80BAyw23bs+azOp8nsH2rSWy/TLiyUq1eTePqJF9m39+a/Hxj76rP06duD\nEiW8KBdSP608pHQQ73/0Ft4liuPh4ckbE95m7ZqfnRbrluPRTF29D4sx9KxbjsHNqqarn7ZmHztP\nxQGQlJJK/JW/2Px8VwCiLl7ltR/DiUlIRATee7AZISWKOi1WR7wyaTo/b9mBr08Jlnwxw6Wx3Gzi\n5Jdo074FiYlJjBzxMgf2Hc7QZszLT9O7bze8vb2oWrZxWvkD/brzymvPER0VC8CcT+Yzf963Tot1\ny76jvDVvORaLhZ6tGjCka2i6+qizF3hl1rdcupqExWLhmT5htKhbleSUFCbOXsqhk5F4iDBmYGca\nVa/gtDjtBbeqTaOJAxEPD47N38CBD9Lf4aD60E5U7tcKk5JKUvwlfhk1iysR5wAY8OfnXDhyGoAr\nEedY/+j0XInZXg4vAGkMHDPGnAAQkQVAdyDjm1FGHYA1xph427FrgI7A/KwO+FclMxEJAN4BmgLn\ngWvAVGPMYpcGlsPatW9JhYrlaFwvjAYN6zB1+gQ6tu2Tod2qFev5dNaXbA9fla581OjhLF2ygjmf\nzqdK1YrM/2YWDWq3dUqsqRbD5JV7mdG/OQFehXlo9npCKwdR0c8rrc3o9rXTHs/feZwj0RfS9l9Z\ntovH7q3KPRUCuHotBRGnhJktPTq3p3+vbox9/W1Xh5JOm3YtKF+xLM0bdqZ+w9pM/s+rdG3fP0O7\ntas2MOeTr9i0c3mGuu8Xr+SVFyY5PdZUi4VJc79n5guPEuDrRf9xM2hVvzoVQ/zT2ny8dAMdGtei\nT7smHI+I5cm3P2dF3ap8u34XAN9OfopzFy8z4u3P+eq1x/HwcO6singITd58mDX9pnA1Kp7Oyydy\nevVuLv4emdYm/sApfuz0KqlJ16gyqC0NXunHz8Pft77mpGv8EPayU2O8LeP4fyARGQoMtSuaZYyZ\nZbcfApy22z8DNMnkqXqJSEvgKDDSGHM6i2NDbhXPv2bOTEQEWAL8bIypYIxpAPQFSt/UzikJ3lnP\nm5mO97Xl6/lLANi9ay/e3l4EBPhlaLd7115iYuIylBtjKF68GABeXsWJjo51WqwHIuMp41uU0j5F\nye/pQYcapdlwNCrL9isOnqZjTeuv7HhcAqkWwz0VAgAoUiAfhfO7/vNZw7p34+1V3NVhZBDWuTWL\nFiwDIHzXPry8iuMfUCpDu/Bd+4iNOZvb4aVz4PgZygSUpLS/L/nz5aNj07vZsPumXqTA5aS/ALh8\nNQm/Etaf+YmIOBrXsPbESnoXo3iRQhw8GYmzlaxXkUunYrj8ZxyW5FROLd1GmQ4N0rWJ+eUwqUnX\nADi7+xhFgnydHld2ZGfOzBgzyxjT0G6bdfszZPA9UM4YUxtYA8z9u7H/a5IZ0Aa4ZoxJG/cxxvxh\njHnPNgm5TETWAT+J1TQROSAi+0XkwevHiMgLtrK9IjLFVlZRRFaKyG4R2SQi1Wzlc0RkhohsB6ba\nJjL9bHUetknRjFnmDgUFBRAZEZ22HxkZTWBwgMPHT5v8Pr37dGXvoY3MXzSLl8a8kdMhpom9lERg\n8cJp+wFehYm9lJhp28iLV4m8cIXG5ayfzv+Iv0zxQvkZtWgbD37yE9N/2k+qG1/BwNkCb/q7iIqM\nITDI8b8LgE5d27Nm03fMnDOdoJDAnA4xTez5BAJ9vdP2/X29iDmfkK7N8Pvb8uOWvbR/eioj3v6c\nFwd1AaBK2UA2hh8hJTWVM7HxHD4VSUz8RafFel2RQB+uRMan7V+NiqdIoE+W7Sv1CyVi/d60fc+C\n+em8fCKdvp+QIQnmFkuKOLw5IAIoY7df2laWxhhzzhjzl233E6CBo8fezPUfY3NPTSD8FvX1gdrG\nmHgR6QXUBeoApbBOXP5sK+sONDHGXLVNUgLMAh43xvwuIk2AD7EmT7D+EpoZY1JF5CLwEPBfoB2w\n1xiToWtk330vVsifQgVK3NELz66eve9jwVeL+ej9z2jYqC4fzpxKi6ZdMMa1iWLVwdO0qx6Cp4f1\nP1KqxfDr6bMsGNKWQO/CvPDdDpbt+4Oedcu5NM5/qjUrN7D02+Vcu5bMQw8/wH8/eJMHewxxWTwr\ntu6jW4t6PNy5OXt//5OXZyzi28lP0SO0Picj4+g/7iOCSpWgTqWyeOSF8Wc75e+/l5J1KrCq140P\nit82eZbE6PMUK+tH2MKxnD9ymst/OG9UJDMmG8OMDtgJVBaR8lgTUV8g3bi2iAQZY64PxXQDrne/\nVwGTROT6p4Ew4KVbnezf1DO0FfMjAAAgAElEQVRLR0Q+sPWudtqK0iYbgebAfGNMqjEmBtiI9Y7k\n7YDPjDFXAWyJrxjQDPhGRPYAM4Egu1N9Y4xJtT2eDQyyPR4MfJZZbPbdd0cT2eDH+rN+0xLWb1pC\nTEwcwXafmoODA4mOjHHoeQAeGtibpYtXALBr5x4KFipIyZJZf8K8E/7FCxFt1xOLSUjE366nZm/l\noTN0rHnjw1qAV2GqBpSgtE9R8nl40LpqMIft5tMUPDykL6s2LmLVxkXE3vR3ERQcQHSU438XF85f\n5Nq1ZADmz/uWu+vWyPF4r/P38SLarjcVG59AgI9XujaLN+6mQ5NaANSpXJa/klM4f+kq+Tw9GT2g\nMwvffJJ3Rw7g0tVE7grKOJya065Gn6do8I1hwyJBvlyNzrj4LqhFTe5+uhvrH5mO5dqNL3Yl2tpe\n/jOO6K2H8a11l9NjvllOLs03xqQAT2JNTIeBhcaYgyIyUUS62Zo9LSIHRWQv8DTwiO3YeOB1rAlx\nJzDR7v05U/+mZHYQa+8LAGPMCKAtcH2Y78rffF4P4IIxpq7dVt2uPu15bRObMSLSButKnxV/85wZ\nzP7kK1q36EHrFj1Y8cNaHuzXA4AGDeuQkHAp07mxrESciaJl6D0AVK5SgUIFC3L27C3/jv62msE+\n/Bl/mYgLV0hOtbDq0BlCqwRlaHfy7CUSkpKpE3LjzaJmkA+XkpKJv2IdpdhxKpYKpfLeXJUrzf10\nAR1Ce9MhtDcrf1xH777W95D6DWtzKeFytubG7OfXwjq15tjREzke73U1K4TwZ/Q5zsTGk5ySwspt\n+wmtXy1dm6CS3mw/aI3hREQs15JT8PUqSuJf17hqm5fauv8Ynp4e6RaOOMu5PScoXj6QYmX88Mjv\nSbnuTTm9Ov1gkG/Nu2g6ZTDrH51O0rkbw6YFvIvgUcA6UFbQpxj+japw8egtR9WcwljE4c2h5zNm\nuTGmijGmojHmTVvZOGPMMtvjl4wxNY0xdYwxrY0xR+yOnW2MqWTbMv3gb+/fNMy4Dmu3dbgx5iNb\nWZEs2m4ChonIXMAXaAmMxrr6cZyIfHl9mNHWOzspIg8YY76xLTSpbYzZm8VzfwJ8Acyz67HlqDWr\nN9IuLJQde9aQeDWRp0eMTatbv2kJrVtYE924iaPp1bsLhYsUZu+hjXzx+TdMm/I+416ewjv/e4Nh\nTzwCxvDUExm+65hj8nl48GKHugyfvwWLxdC9zl1U8vPiw42HqBFUglZVggFYeeg0HWuURuyGizw9\nhJFtazHsq00YA9WDStCrXnmnxeqo0eOnsPPXfVy4kEDbHgN4YshAenXt4OqwWLfmZ9q0b8Hm3StI\nSkxk1JOvptWt2riIDqG9AXh5wih69O5M4SKF2HlgLfPnfcf0tz5k8NABtO/UitSUVC6cv8jIEa84\nLdZ8np68NKgLw6fNxWKx0KNlAyqVDuCDb9dSs3wIrepX57n+nZj46RK+WPkLIjBx6P2ICPEJVxg+\ndS4eHoK/T3HefLy30+K0Z1It7HhlLu2+GmNdmv/1Ri4ejaDO8704t/ckZ9aE0+DVfuQrWojQmU8D\nN5bge1cOoemUwRhjQcSDA+9/n24VZG5x8UzCHRFXz4PkJhEJwro0vwkQh7XXNAMoDDS8/sU9W0Ka\nCnTC+oW9N4wxX9vqXsQ6VHgNWG6MGWsbE/4I6/BifmCBMWaiiMwBfjDGLLKLIT9wDmhs/ykkK37e\nVd3qF/Tne/e7OoRsy9f5/1wdQraUr9Lt9o3ymGMrx7s6hGxb2HOpq0PItkERX9zRpNcf9ds5/H5z\nV/jaPDUR+W/qmWGbaOybRfUcu3YGa09sdCbPMQWYclPZSaxf6Lu57SOZnKcO1oUft01kSimVmyyp\neSo/Zcu/Kpm5mq1XNxzrikallMpTHJ0Ly4s0meWizHp1SimVV+Tw0vxcpclMKaUUoDfnVEop9Q9g\n0Z6ZUkopd2dJdd+vHmsyU0opBbj398w0mSmllAJ0NaNSSql/AJ0zU0op5fZ0ab5SSim3p3NmSiml\n3F6qRVczKqWUcnPaM1NO80lh19w+/e8Ke2mHq0PItpPPfu/qELLl5NFlrg4h2/zLhbk6hGx7pmQT\nV4eQ63QBiFJKKbenC0CUUkq5Pe2ZKaWUcntuPGWmyUwppZSVrmZUSinl9tz4DjCazJRSSlkZdM5M\nKaWUm7O48aSZJjOllFIAWLRnppRSyt2lajJTSinl7nTOTCmllNvT1YxKKaXcniYzpZRSbs+dhxnd\n9+veyiH+rWvTdvPbtN06ncpPds1QX25QW1qvn0KrtZNovnQ8xauE5Fpsz0wcwfzNnzNnzcdUqVU5\n0zZV7q7MnLUfM3/z5zwzcUSG+geHPcCmiJ/w9vECoH3PtsxZ8zFz1n7Mh0v/R8UaFZwW/8TJL7F5\n13LWbPqOWrWrZ9pmzMtPs2P/Wn77M/3dBB7o1529R39m1cZFrNq4iH4DezktTke8Mmk6Le/rS48B\nj7s0jptNmfYqu/f+xOZtP1C7Ts1M27wyfhQHjmzidPTedOVlygSz5IfP2bztB75f8SXBwYFOj7dS\naG2eXDeNpzf+h+bDM/5/u6txNYb9+Abjjn9Ojc6N09V5B5dk4LwXGfHTVEasnUqJ0qWcHu/NLOL4\nltf865OZiKSKyB67rZyrY8oxHkLtyY+ytf9U1rUcTUjPZhmS1ZnvfmF96xfZ0G4sxz74npoTBuRK\naE3bNKZ0+dL0az6IqS9M57nJz2Ta7rnJzzJ1zHT6NR9E6fKladL6xhuAf7AfjVs2IPpMTFpZ1Oko\nnuw9kkfa/R9z//sFY94a5ZT427RrQfmKZWnesDMvjJzA5P+8mmm7tas20KVd30zrvl+8kg6hvekQ\n2pv58751SpyO6tG5PTOmv+HSGG7WPiyUihXL0aBOW5596hX+89/XMm23cvk62oben6F84qSXWPDV\nYpo37cLUKe8x7rXnnRqveAidX3+ELx+eygftxlCr2z34VU7//+1i5FmWPDeT/Ut/yXB8z+mPs2Xm\nD3zQdgwfd3uVK2cTnBpvZlIRh7e85l+fzIBEY0xdu+2UfaWIuO1QrE+9Slw5GcPVP2MxyalELNlK\nYIf090dLuZyY9tizSEFy61KjzTvcy8pFqwE4FH6YYt7FKOnvm65NSX9fihYvwqHwwwCsXLSaFh3v\nTat/asITfPjmLIzdHQUP7DrE5YuXATgYfgi/ID+nxB/WuTWLFljvKxa+ax9eXsXxD8j4STp81z5i\nY846JYac1LDu3Xh7FXd1GOl07tKOBfMXA7Br5x68vb0ICMj4+9y1cw8xMXEZyqtWq8SmjdsA2LRx\nG53ua+fUeEPqViT+VAznT8eRmpzKge+3UbV9+v9vF86cJebIacxN3072qxyCRz5PTmw+AMC1q3+R\nnHTNqfFmxpKNLa/RZJYJEXlERJaJyDrgJ7GaJiIHRGS/iDxoazfRrkcXISKf2coHiMgOW/lMEfG0\nlV8WkTdFZK+IbBORAGe+jkJBPiRGnkvbT4yKp1CQb4Z25R9tT7tt71Dz1f7sf/lzZ4aUxi+wFLGR\nN96A4qLiKBWYPhmUCixFXJR9m7P42do0D2tGXNRZjh86keU5uvTtxPb1zrlZaGBQAJER0Wn7UZEx\nBAZl79fZqWt71mz6jplzphMU4vwhMHcTFBRAxJmotP3IyGiCgh3/GR/cf5gu3a03Be3SLQwvr2L4\n+JbI8Tiv8wr0JSHqxv+3hKh4vAJ9HDq2ZPlAkhKu8uDMZxm2/E3aj+2HeOR+78ci4vCW12gyg8J2\nCWmxXXl9oLcxJhS4H6gL1AHaAdNEJMgYM84YUxdoBcQD74tIdeBB4F5bXSrwkO05iwLbjDF1gJ+B\n/8uF13dbJz9bw9qmIzn4xnyqjOzh6nBuq2Chggx8qj+fvj0nyzb1mtXlvn6d+GjSx7kXWDasWbmB\ne+qG0b7F/fy8fiv//eBNV4f0j/Pq2Cnc27wxG7cs497mjYmIiCY1NdXVYWXKI58nZRtVZfUbX/Jx\n11fxKetP3Qda5nocJhtbXuO2Q2g5KNGWdG62xhgTb3vcHJhvjEkFYkRkI9AIWCYiAnwBTDfG7BaR\nJ4EGwE5rFYWBWNvzXAN+sD3eDbTPLCARGQoMBRhevBEdilT6Wy8sKeo8hYNLpu0XDvIlKSo+y/YR\nS7ZS563B/Pq3znZ7PR/uTteHOgNwZM9v+AffGDLyC/LjbHT64biz0WfTDRP6BZUiLvosIeWCCSob\nyGdrZqUd++mqGQy9bwTxceepWL0CL0x7jtEDXyLhfM7NOzw8pC/9B/UGYO+vBwi2600FBQcQHRWT\n1aEZXDh/Me3x/Hnf8vJrzpnbczePDR3AoEf6ABC+ez8hpYPS6oKDA4mKdPxnHB0dy6D+1kVDRYsW\noWv3jiRcvJSzAdtJiI7HK+jG/zevIF8Sos87dmxUPNGH/uD8aetIxJFVuyldvxK/fr3RKbFmJS8O\nHzpKe2ZZu+JguwnAGWPMZ7Z9AebazcFVNcZMsNUlmxsTPKlk8WHCGDPLGNPQGNPw7yYygAt7jlO0\nQiBFyvoh+T0J6XEP0at3p2tTtPyNN+SAdvW4cjL65qfJMYvnLmVw2DAGhw1j06otdOxtHQKqUb86\nlxOucC42faI9FxvPlUtXqVHfulKwY+8wNq/awokjJ+lWpzd9mj5En6YPERcVx5AOjxMfdx7/YH/e\n+HgCbzwzmdMnzuRo/HM/XZC2YGPlj+vo3bcbAPUb1uZSwuVszY3Zz6+FdWrNsaNZD5f+m3wy6wta\nNutGy2bdWP7DGvr26wlAw0Z1SUi4lOncWFZ8S/pg+0DJyOcf58t53zgl5usi956gZPlASpTxwzO/\nJ7W6NuW3NbtvfyAQsfc4hbyKUMTXOm9ZvlkN4n6PcGa4mUoRcXhzhIh0FJHfROSYiLyYSf0oETkk\nIvtE5CcRucuuzn5x3rLbnUt7Zo7ZBAwTkbmAL9ASGC0iXbEOO7a2a/sTsFRE3jHGxIqIL1DcGPNH\nbgdtUi3sGzuHe+a/iHh68Of8DVz6LYJqY3pzYc8JoleHU35wGH4ta2GSU7h28QrhT3+UK7Ft/Wk7\nTds0YcGWeSQlJjF51LS0utmrZzI4bBgA08e+y9h3xlCwUEG2rd/BtnW3ngN7dORAvH28GDXJujoy\nNSWV/+v8RI7Hv27Nz7Rp34LNu1eQlJjIqCdvrGZctXERHUKtPbiXJ4yiR+/OFC5SiJ0H1jJ/3ndM\nf+tDBg8dQPtOrUhNSeXC+YuMHPFKjseYHaPHT2Hnr/u4cCGBtj0G8MSQgfTq2sGlMa1etYH2HVoR\nvm8diYmJjHj8hbS6n39ZRstm1g8Tr70+hl59ulGkSGEO/LaZeXMX8tak/9G8RRPGTXgeYwy/bNnJ\n6FETnBqvJdXC8nFzGPj5C4inB78u3Ejc7xG0HtWLyH0n+W1tOMG1K9B31kgKeRehSrt6tBrZiw/b\nv4CxGFa/+RUPfzUWRIjaf5Lw+eucGm9mcnL40LZW4AOsI1BnsI5WLTPGHLJr9ivQ0BhzVUSGA1Ox\nTtNA1qNmmZ/PfiXYv5GIXDbGFLup7BGsP+AnbfuC9YfcCevv+w1jzNcish4oD1ywHbrMGDPOtkDk\nJaw932RghDFmm/25RKQ30MUY88it4lsa2N+tfkFvezo+DJRXnEx0r5hPHr3th9Q8x79cmKtDyLZn\nSjZxdQjZNuGPL+9oZcbnIQMcfr8ZFPHFLc8lIvcAE4wxHWz7LwEYYyZn0b4e8L4x5l7bfob35lv5\n1/fMMvthGWPmAHPs9g0w2rbZt2tNJowxXwNf3+pcxphFwKK/GbZSSuW47MyZ2c/t28wyxsyy2w8B\nTtvtnwFu9QlhCLDCbr+QiOwCUoApxpglt4rnX5/MlFJKWWVnGMiWuGbdtqEDRGQA0BAItSu+yxgT\nISIVgHUist8Yczyr59BkppRSCsjxy1RFAGXs9kvbytIRkXbAy0CoMeav6+XGmAjbvydEZANQD8gy\nmelqRqWUUoB1PM/RzQE7gcoiUl5ECgB9gXQTvrZ5splAN2NMrF25j4gUtD0uBdwL2C8cyUB7Zkop\npQAwOdgzM8ak2L53uwrwBGYbYw6KyERglzFmGTANKAZ8Y/saxZ/GmG5AdWCmiFiwdrqm3LQKMgNN\nZkoppYCc/9K0MWY5sPymsnF2jzO9YKYx5hfg7uycS5OZUkopwL2vAKLJTCmlFJA3r7noKE1mSiml\ngLx5001HaTJTSikFOLxKMU/SZKaUUgrQYUallFL/ADrMqJRSyu3pakblNL3ic/fmfHcq4U3X3jbk\n7/AIfdzVIWSLO16BPvbUaleHkG1TGrx6+0b/MDrMqJRSyu2luHE602SmlFIK0J6ZUkqpfwCdM1NK\nKeX2dDWjUkopt2dx44FGTWZKKaUAnTNTSin1D6CrGZVSSrk9901lmsyUUkrZ6GpGpZRSbk8XgCil\nlHJ77pvKNJkppZSy0WFGpZRSbi/Vjftmmsz+od6ZPpFOHdtwNTGRIUNG8uueA+nqCxcuxNfzZ1Gh\n4l2kpqby449rGPvy5HRtevbszDdff0yTpp3YHb4vV+L2KFeLAm37gwgp+zaRsmN5uvr8rfviWbaa\ndSdfAaSIF4nvPZkrsdnbsu8ob81bjsVioWerBgzpGpquPursBV6Z9S2XriZhsVh4pk8YLepWJTkl\nhYmzl3LoZCQeIowZ2JlG1SvkWtxTpr1K+7BWJCYm8sSwF9i392CGNq+MH0Xffj3xLuFFmcA6aeVl\nygTz3kdTKFXKl/PnLzJsyHNERkbnWuwZ4pw0nZ+37MDXpwRLvpjhsjjsVQytTYfxAxFPD35dsIFf\nPvo+XX3ZxtUIGz+AgGpl+e6p9zm8fAcAd91Tg7BXB6S1K1UxiO+eep/fVu/O1fjdec7Mw9UBqJzX\nqWMbKlcqT7UazRk+/AU+eH9ypu2mvzODWneH0rBRB5rd04iOHVqn1RUrVpSnnxzC9u3huRU2iFCg\n/QD+WvQOSbNfIV/1JkjJ4HRNktcvIGnuBJLmTiAl/CdSf8/d/+wAqRYLk+Z+z4ejB7H4radZuXU/\nxyNi07X5eOkGOjSuxcI3RvDWiAeZNNf6pvbt+l3Wfyc/xYwXHuE/X63EYsmdwZ32YaFUrFiOBnXa\n8uxTr/Cf/76WabuVy9fRNvT+DOUTJ73Egq8W07xpF6ZOeY9xrz3v7JBvqUfn9syY/oZLY7AnHkLH\n1x/hq4en8lG7MdTqdg+lKoeka3Mx8izLnpvJgaW/pCv/Y+shPu48lo87j2VevzdJTrrG8Z/352L0\nViYbW15z22QmIqkiskdEDorIXhF5TkRueZyIlBOR/jkXZvbdSQx2r/mAiHwjIkVyOj5n6tq1A/O+\nXATA9h3heJfwJjDQP12bxMQkNmy0/odKTk4m/Nf9hIQEpdW/NmEM097+kKSkpFyL2yOoAuZ8LOZi\nHFhSSTmyHc9KdbNs71m9CSmHt+dafNcdOH6GMgElKe3vS/58+ejY9G427D6cvpHA5aS/ALh8NQm/\nEsUBOBERR+Ma1p5YSe9iFC9SiIMnI3Ml7s5d2rFg/mIAdu3cg7e3FwEBfhna7dq5h5iYuAzlVatV\nYtPGbQBs2riNTve1c27At9Gw7t14exV3aQz2gutW5PypGC6cjsOSnMrB77dRtX2DdG0unjlL7JHT\nGEvW6aB658Yc27CXlKRrzg45AwvG4S2vcaRnlmiMqWuMqQm0BzoB429zTDkgW4lERHJ6yDPbMdi5\n/pprAdcAp9290Qmvm5DgQM6cvvEGGXEmipDgwCzbe3t70eW+9qxbvxmAenVrUaZMEMtX/JTTod2S\nFCuBuRSftm8unUeK+WTe1qskHt6lsPx5ONN6Z4o9n0Cgr3favr+vFzHnE9K1GX5/W37cspf2T09l\nxNuf8+KgLgBUKRvIxvAjpKSmciY2nsOnIomJv5grcQcFBRBxJiptPzIymqDgAIePP7j/MF26W28M\n2qVbGF5exfDxLZHjcborr0BfEqLOpe0nRMVTPDDzv99bqdntHg4u3ZqToTnMko0tr8nWMKMxJhYY\nCjwpVp4iMk1EdorIPhEZZms6BWhh692MzKqdiLQSkU0isgw4ZCt7VUR+E5HNIjJfRJ63lVcUkZUi\nstt2TDVb+RwR+Z+I/CIiJ0SkdxYx1BSRHbb9fSJS2cGXvQmoZDvXEtv5D4rI0OsNROSyiLxjK/9J\nRPwciHmGiGwHpmbnd5DTPD09+XLeB7z/wWxOnvwTEeHtaeMZPWaiK8O6Lc9qjUk5ugtM3vuECLBi\n6z66tajHmv+N4YPnB/HyjEVYLBZ6hNYnwNeb/uM+YtqXy6lTqSwe4h6XKn917BTubd6YjVuWcW/z\nxkRERJOamurqsP5RivmXwL9qGY7/nDtz1DdLxTi85TXZ7hUYY06IiCfgD3QHLhpjGolIQWCLiKwG\nXgSeN8Z0AbC98WfWDqA+UMsYc1JEGgG9gDpAfiAcuD4pMgt43Bjzu4g0AT4E2tjqgoDmQDVgGbAo\nkxjeA941xnwpIgUAz9u9VluvqROw0lY02BgTLyKFgZ0i8q0x5hxQFNhljBkpIuOw9lyfvE3MpYFm\nxpgM7wa2n9dQAPH0xsOj6O1CZfjjDzNkyEMA7Nq1h9Jlbsw1hZQOIiKLifoZH03l92Mn+d97nwBQ\nvHgxatasxk9rrMOUgYF+LP7uM3re/6jTF4GYyxeQ4r5p+1LcB3P5fKZt81VrzLW1Xzg1nqz4+3gR\nbdebio1PIMDHK12bxRt389HoQQDUqVyWv5JTOH/pKiW9izF6QOe0doNem8ldQaWcFutjQwcw6JE+\nAITv3k9I6RtDycHBgURFxjj8XNHRsQzqPwKAokWL0LV7RxIuXsrZgN1YQnQ8XkEl0/a9gny5FJ35\n329WatzXhN9W7cKS4poPCSYPJilH3ekCkDBgkIjsAbYDJYHMejy3arfDGHPS9vheYKkxJskYcwn4\nHkBEigHNgG9szzETawK7bokxxmKMOQRkNW6yFRgrIi8AdxljEm/xugrbzrML+BP41Fb+tIjsBbYB\nZexegwX42vb4C6C5AzF/k1kiAzDGzDLGNDTGNHQkkQF8NGMuDRuF0bBRGMuWrWLgQ9YOapPG9Um4\nmEB0dGyGYya+NgZv7+KMeu7GqHFCwiUCg++mUpWmVKrSlO3bw3MlkQFYok4iPgGIdynw8CRftSak\nHtuToZ34BkKholgijzs9pszUrBDCn9HnOBMbT3JKCiu37Se0frV0bYJKerP94AkATkTEci05BV+v\noiT+dY2rtrmQrfuP4enpQcUQ/wznyCmfzPqCls260bJZN5b/sIa+/XoC0LBRXRISLmU6N5YV35I+\niK0XOfL5x/ly3jdOidldRe49gW/5QEqU8cMjvyc1uzbl6JrsLVCq2a0ZB5a5ZogR3HuYMds9MxGp\nAKQCsYAATxljVt3UptXNh92i3RUHTusBXDDGZLUa4K+bzpWBMeYr27DefcByERlmjFmXxfMl3nwu\nW6ztgHuMMVdFZANQKIvjjQMxO/K6/5blK36iY8c2/HZ4C1cTE3nssVFpdbt2rqZhozBCQoIY+9Iz\nHD7yOzt3WH8tH374GbM/m++ssG7PWLi29gsK9h4FHh6k7N+MORdJ/nt7YIk+Repxa2LLV60JqUd2\nuCzMfJ6evDSoC8OnzbUOHbZsQKXSAXzw7Vpqlg+hVf3qPNe/ExM/XcIXK39BBCYOvR8RIT7hCsOn\nzsXDQ/D3Kc6bj/e+/QlzyOpVG2jfoRXh+9aRmJjIiMdfSKv7+ZdltGzWDYDXXh9Drz7dKFKkMAd+\n28y8uQt5a9L/aN6iCeMmPI8xhl+27GT0qAm5FntmRo+fws5f93HhQgJtewzgiSED6dW1g8viMakW\nVo6bQ//PX0A8Pdi7cCNxv0cQOqoXUftOcnRtOEG1K9Bn1kgKeRehcrt6hI7sxYz21t+Dd+lSeAX7\n8se23J8Hvs6SR4ftHSHmNsGLyGVjTDHbYz/gS2CrMWa8bTisM/CAMSZZRKoAEViH+6YbY0Jtx2XV\nrhHphwIbYe3BNMOaaMOBWcaYt0XkF+AdY8w3Yv14WNsYs1dE5gA/GGMW2ccrIg1uiqECcNIYY0Tk\nbeCMMea/t3vNdmXdgceMMV1tc197gI7GmA0iYoB+xpgFIvIKEGCMecrRmG8lX4EQt/rrSnjTdW8m\nf5dHaCdXh5AtQW1ecnUI2RZ7avXtG+UxUxq86uoQsu3VP768ownYAXfd7/D7zRd/fJenJnsd6Zld\nH3LLD6QA84DptrpPsK4aDLe9WccBPYB9QKptSG4O8G4W7dIxxuy0LQbZB8QA+4HrkxMPAR/ZkkV+\nYAGw9xZx3xxDQWCgiCQD0cAkB167vZXA4yJyGPgN61DjdVeAxrbYYoEH/2bMSinlMnlxyb2jbpvM\njDFZLpQwxliAsbbtZm1u2s+s3QbbZu9tY8wEsX6362dsC0Bs82odM4nhkZv2i9n+Tc4khimZv5IM\nz1ksk7K/sC4GyeqYUZmUORSzUkrlBXlxlaKj8uLlrGaJSA2s81FzjTG5eAkKpZT69/pH98xymzEm\nV64cIiIlgcy+FdzWttzeYZn15JRSyt38m5fmuy1jzDnbVT5u3rKVyJRS6p8ip5fmi0hH20UwjonI\ni5nUFxSRr23120WknF3dS7by30TktivL/rXJTCmlVHrGGIe327FdXOMDrGsNagD9bFNI9oYA540x\nlYB3gLdsx9YA+gI1sa47+ND2fFnSZKaUUgrI8QsNNwaOGWNOGGOuYV3N3f2mNt2BubbHi4C2thXv\n3YEFxpi/bAvpjtmeL0uazJRSSgE5fm3GEOC03f4ZW1mmbYwxKVi/ilXSwWPT0WSmlFIKyF7PTESG\nisguu23o7c/gPHluNaNSSinXcGQuzK7tLKwXU89KBNZr2F5X2laWWZsztgu7ewPnHDw2He2ZKaWU\nAnJ8NeNOoLKIlP//9lZw/KEAABtySURBVO48vqrqXOP470mYDHMQBBwQELSIDIqzrXgVKFioinrF\noWq92PY629b2UhUL0qteS69TFa3XqWodagUVBVGBglaUSVAUFFBGGQIBRKbkvX/snXASM2qSvfc5\n77effMzZZx/ycFjNe/Zaa68V7lRyHsGuJqkmABeH358NvGlBRZ0AnBfOduxIsKh7hQuy+pWZc845\noGbvMzOzPZKuBCYRbLn1f2b2oaRRBFtmTSDYkeQJSZ8CeQQFj/C8Zwn2udwDXFHeLiNFvJg555wD\noMBqdnMXM5sITCx17OaU73cA55Tz2jHAmKr+rEpXzXfRGtXhgkT9Aw3csz3qCNW2iKrtGRcXS+sl\nqkkAVdgJN4Z+O3t01BGqrf6+nb7TSvanHNCvyo3rrZWvJ27VfOeccxkgyctZeTFzzjkHJHtzTi9m\nzjnnABJ8XebFzDnnXMi3gHHOOZd4NT2bsS55MXPOOQf4lZlzzrk04LMZnXPOJV6S7zv2Yuaccw7w\nbkbnnHNpwCeAOOecSzwfM3POOZd4vgKIc865xPMrM+ecc4nnV2YuVjqf3IMBIy8iKzuLuX+bysz7\nXyrx/EHHHMaAkRey32EH8fer7mXRxL0buDZr34rBtw+nWftcMHjqkjvIX7mhVvM269ubg0ZdhrKy\nWP/0FNbe90KJ51tfNIA2Fw+EwkIKvtrB8hv+zI4lK1G9bA6+8wpyundC9bLZ+PxbrLn3hXJ+Ss1q\n37cHR4+6CGVl8enTU1l4X8n3+HuXD6TLsL7YngJ25G3l7esf5KtVGwG48IvH2fzxCgC+WrWRty4d\nWyeZDzm5Bz8M28Wcv01lRql20eGYw/hh2C6ev+pePkppF83bt2JISrt48pI72FzL7aKoHStsx2+X\n0Y77h3lfSGnHHY7vRv+bLiw+b9/O7Xjhqnv5ZPLsWs1bmRv/MJbpM2eR27IFL/71gUizlMevzOqQ\npAJgAVCfYAfSx4E/mZU/DUfSwcDLZtZdUh/gJ2Z2dR3ErXPKEgNHX8JfL/hvtqzN4z8mjOaTKXPY\nsGRV8Tn5qzcw/pfjOP7y07/x+jPG/pwZ945n6YyF1M9piBXWcuPOyqLDmMtZPOwWdq3ZSLeJd7B5\n8ix2LFlZfMrGf0xn/ROTAGjR72gOGnkpiy8cTcsfnYAa1OPD064lq1EDuk+9h40v/pNdK9fXamRl\niWPHXMzrw25j+5o8Bk0cxYrJs8lfsrr4nLyFy3ll4E0U7NhF15+cylE3DmP6L+4FoGDHLl7u/7ta\nzVhW5kGjL+GJsF0MD9vF+lLt4sVfjuOEMtrFmWN/zvSwXTSog3ahLPHD0ZfwZEo7XlxGO55QRjv+\n/J2PeGjQCAAaNW/MldPH8tn0BbWatyrOGNSP84cOYcToO6OOUq4kz2bMijrAt/C1mfUys8OBfsBA\nYGRVX2xm76drIQPYv1dnNi3/ks0r1lO4u4APX/oXh/Y7qsQ5+Ss3sO7jFd/4hbRvl/3JqpfN0hkL\nAdi9fSd7duyq1byNe3dh5/I17PziS2z3HvLGz6DlgGNKnFO47evi77NyGlLcE2JGdk4jyM5C+zTE\ndu+hIOXc2tKqd2e2Lv+SbV8E7/Hy8f/iwAEl3+Mv315EQfjebZj9KTntcms9V0X279WZvOVfsmnF\negp2F7CwjHaxeeUGviyjXbQu1S52bd/J7lpuF+2/QztO9b1Bx/Dp1Pm13o6rok+vI2jerGnUMSpU\naFblr7hJYjErZmbrgMuBKxXIlvQ/kt6T9IGkn5V+jaS+kl4Ov28i6RFJC8Lzh4bH75f0vqQPJf0+\n5bW3SfooPPfO8Ng5khZKmi9peniszBzhz54q6XlJH0t6UlKN7tbatG0u+Ws2Fj/esiaPpm1bVum1\nrTq2ZceW7Zwz7lqGTxzDaSOGoaza3Uy2Qdtcdq3e2121a81G6rdt9Y3z2lw8kCNm3s+BN17MFzf/\nBYBNr7xDwfYd9Jr7f/Sc9SBrH3iRgs3bajUvQE7blny1Oq/48fY1eeRU8B4fMuxkVr01v/hxdsP6\nDJo4ioEv3fKNIlhbmrXNZUupdtGsmu3i38ddy88mjqFfHbSLsvJWtR2nOnzI8Xw4/p2ajJbWrBr/\ni5vEdTOWZmZLJWUDbYAfA/lmdrSkhsBMSZMpf5uem8LzjwCQVPT/lt+ZWV74574hqQewCjgTOMzM\nTFKL8NybgQFmtirl2GXl5ADoDRwOrAZmAicCM1JDSbqcoEgzOPcY+jQ55Nu+PdWSVS+bg44+lAcH\njSB/9UbOvu8qep7zA+Y9M61Ofn5F1j32Kusee5XcM75P+2vOYdm1d9O4VxcoKGT+kZeR3bwJh/1j\nDFv++QE7v/gy6rjFOp51Iq16dmLS0FuLj/392Gv5eu0mmhzUmv7PjmDTxyvY9vm6CFNWrKhdjEtp\nF73O+QFzY9AuKtKkTQvaHHogn03/IOooiVHBaE3sJfrKrAz9gZ9Imge8C7QCulRw/mnAfUUPzGxT\n+O25kuYAcwkKTzcgH9gBPCzpLGB7eO5M4FFJw4HsKuSYZWYrwzG+ecDBpUOZ2YNm1sfM+lS3kG1d\nm0fzdnuvbJq1y2Xr2k0VvGKvLWvy+PKjz9m8Yj1WUMjHk2bTrnvHav386tq1No8G7fctftygXSt2\nr91Y7vl542fQIuyGzD3zB+RPnYvtKWDPxny2vfcxOT0712pegO1rN9G4/d5uw5x2uWwv4z1u9/3D\nOeLqIbx1yVgKd+0pPv51eO62L9az9p1F5HbvUOuZt6zNo1mpdrGlGu1i7Uefs2nFegrrqF2Ulbeq\n7bhIt9OP5ZNJ71O4p6Cm46WtQqzKX3GT+GImqRNQAKwDBFwVjqn1MrOOZja54j/hG39eR+BXwKlm\n1gN4BWhkZnuAY4DngR8BrwGY2c+BG4EDgdmSWlWSY2fKjyughq+OV81fSm7HtrQ4sDVZ9bM5fPBx\nLH69arO4Vs//jIbNcsjJDfr1O57QrcQEgdrw1bwlNOzYjgYHtkH165H745PYNPm9Euc07Niu+Pvm\npx3FzmVrANi1aj1NTzwCgKx9GtLkyK7s+LR28wJsnLeUph3b0iR8jw/+8XGsmDynxDm5h3fguNt+\nyluXjmXHxi3Fxxs0zyGrQfBP3rBlE9oc3ZX8xbWfefX8pbQK20V2/Wy6Dz6OT6rYLlbN/4xGddwu\nVn+Hdlzk8CEnsHCCdzFWh5lV+StuEt3NKKk18ABwb9j1Nwn4haQ3zWy3pK4E3YPleR24Arg2/PNa\nAs2Ar4B8SfsRTDCZKqkJkGNmEyXNBJaGr+lsZu8C70oaSFDUqpujxlhBIa/e/CgXPP4blJ3FvGen\nsX7JKvpeP5TVHyxj8ZQ5tO/RiXMfvI5GzXPoelpvTr5uKA/0+w1WaEwZ8xQXPTUCJNYsWMacp9+s\n3cAFhXxx40Mc+tRIyMpiwzNvsGPxCtr/ahjb53/K5tffY79LBtHs+z2CK7D8bSy99m4A1j36Kh3/\ndBXd37wLJDY88yZfL/q8dvMSvMezbnyM0566IZia/8w08hevouevhrJx/jJWvj6Ho24aRr3GjTh5\nXDDXqGgKfvMu+3PcbT/FrBApi4X3vlRiFmRtKSwoZOLNj3JR2C7mhu3ilLBdfBK2i/NS2kXf64by\n57BdTB7zFBfXYbuwgkJeu/lRzg/zzg/znnz9UNaE7bhdSjvuktKOAZofsC/N2ufy+b8W1WrO6vj1\nyNt4b+4HbN68hVPPuJD/vOwihg4eEHWsEpI8m1FxrLAVKWNq/hPAWDMrlJQF3AoMJrg6Wg+cAbRk\n79T8vsCvzOxHYYG6DziK4Crp92b2gqRHgROAFQTdixMICtR4oFH4Z99pZo9JeoGgC1HAGwSFUeXk\n6F30s8O/y73A+2b2aHl/31EdLkjUP9DAPdsrPylmFtE46gjVsrReopoEsLf/PUl+O3t01BGqrf6+\nnb7TzJx2LbpVuXGt2fxR7c4CqqbEFbNM48Ws9nkxq31ezOrGdy1mbVt8r8qNa+3mRbEqZonuZnTO\nOVdzknxx48XMOecc4JtzOuecSwMFhcmdAOLFzDnnHODdjM4559KAdzM655xLPL8yc845l3hxXA2/\nqryYOeecA3xzTuecc2nAZzM655xLPL8yc845l3g+AcQ551ziJbmY+ULDGUrS5Wb2YNQ5qiNpmZOW\nFzxzXUha3qRI/Oac7lu7POoA30LSMictL3jmupC0vIngxcw551zieTFzzjmXeF7MMlcS++yTljlp\necEz14Wk5U0EnwDinHMu8fzKzDnnXOJ5MXPOOZd4Xsycc84lnhcz52qQpFZRZ0h3krIl3Rl1Dhcv\nvpxVBpHUFfg10IGUf3sz+7fIQlVC0igzuznlcTbwuJldEGGsivxL0jzgEeBVS8AMK0mtgeHAwZRs\nFz+NKlNFzKxA0klR56guSSfwzff48cgCpRkvZpnlOeAB4CGgIOIsVXWgpP8ys/+W1BB4FpgbdagK\ndAVOA34K3C3pWeBRM1scbawKjQf+CUwhOe1irqQJBG36q6KDZvZCdJHKJ+kJoDMwj73vsQFezGqI\nT83PIJJmm9lRUeeoDkkCngQWAKcAE83sf6NNVTWSTgH+CjQG5gO/NbN3ok31TZLmmVmvqHNUh6RH\nyjhscb2alLQI6JaEK/Wk8mKWASTlht9eDawD/gHsLHrezPKiyFURSUemPKwPjANmAg8DmNmcKHJV\nJhwzuxC4CPiSIO8EoBfwnJl1jDBemSTdCrxtZhOjzpKuJD0HXG1ma6LOkq68mGUAScsIujRUxtNm\nZp3qOFKlJL1VwdMW13E+SYuBJ4BHzGxlqed+Y2a3R5OsfJK2Elw97gq/RPAeN4s0WAXC8d/7gf3M\nrLukHsAQM7s14mhlCttzL2AWJT9IDoksVJrxYuZcDQknp9xhZr+MOku6kzSNYDLTODPrHR5baGbd\no01WNkknl3XczKbVdZZ05RNAMoikc4DXzGyrpBuBI4HRZha7CRWSrq/oeTMbW1dZqiqcZXdC1Dmq\nKxyXvADoaGajJR0ItDOzWRFHq0iOmc0KohfbE1WYypjZNEkdgC5mNkVSDpAdda504veZZZabwkJ2\nEsGMu4cJZjfGUdNKvuJqnqQJki6SdFbRV9ShKvFn4Hjg/PDxNuC+6OJUyQZJnQm6z5F0NhDb8ShJ\nw4HnCcZ+AfYHXowuUfrxK7PMUjQl+HTgQTN7JRz8jx0z+33UGb6lRsBGIHVMz4BYThkPHWtmR0qa\nC2BmmyQ1iDpUJa4gWH3+MEmrgGUEV5dxdQVwDPAugJktkdQm2kjpxYtZZlklaRzQD7g9vG8rllfn\nkm4wszsk3UP46TuVmV0dQaxKmdmlUWf4FnaH431FVzmtgcJoI1XMzJYCp0lqDGSZ2daoM1Vip5nt\nKuoWlVSPMtq1+/Zi+YvM1ZpzgUnAADPbDOQSDKLH0aLwv+8Ds8v4iiVJXSW9IWlh+LhHOD4ZZ3cT\n3K7RRtIYYAbwh2gjVUxSK0l3E9zsPVXSXTFfSmyapBHAPpL6Edzs/VLEmdKKz2bMMOF4WRczeyT8\nBN7EzJZFnStdJG2WXRFJhwGnEkzLf8PMFlXykkhJeh2YTnBTOgRdjH3N7LToUpVPUhZwGdCf4D2e\nBPzFb6KuOV7MMoikkUAf4FAz6yqpPcGNvCdGHO0bwqWKyhXX+3MkvWdmR0uam1LMYr3ChqSHgXvM\nbF7KsVvM7JboUlWsrA8IkhaY2RFRZaqqcBGDA8zsg6izpBMfM8ssZwK9gTkAZrZaUlxnBh4PrACe\nJhg0L+uG7zhK1Cy70ACgj6Q/pix8OwS4JbpIlZos6TyCtToBzia42oklSVMJ3tN6BN3k6yS9bWbX\nRRosjfiYWWbZFXZrFP2ibRxxnoq0BUYA3YG7CCatbDCzaTG/0fQKgunXRbPsrgV+EW2kSq0DfgCc\nI+m+cHJC3D88DAeeYu+qJX8DfiZpq6QtkSYrW3Mz2wKcRbDrw7EE3bquhngxyyzPhrMZW4T3vUwh\nWEE/dsyswMxeM7OLgeOATwkG+q+MOFqFzGxpOG7TGjjMzE4ys+URx6qMzCzfzAYD64GpQPNoI1XM\nzJqaWZaZ1Qu/ssJjTWO6DFc9Se0IJmG9HHWYdOTdjBnEzO4MZ1JtAQ4Fbjaz1yOOVa7w1oHTgWEE\n+0AVzbqLLUl/IFjSanP4uCXwSzOL84zG4vFJM7tF0mwg9t1f4XvbheDePgDMbHp0iSo0iqAbdIaZ\nvSepE7Ak4kxpxSeAZIjwPqIpZnZK1FmqQtLjBF2ME4G/mdnCiCNVSerEj5Rjc8zsyPJeEweS9gOO\nDh/OMrN1UeapjKT/AK4BDiDYI+w44J24LkDtap93M2YIMysACiXFuvsoxYUEn7qvAd6WtCX8iuuY\nSJHs8IoSAEn7AA0rOD9yks4lWM39HIJusHfDiStxdg1B8f08/IDWG9gcbaTySbpDUjNJ9cP7ENdL\nujDqXOnEuxkzyzZgQXiPTuruvLFbTcPMkvpB60ngDe3dPPJS4LEI81TF74Cji67GwvsPpxCsJRhX\nO8xshyQkNTSzjyUdGnWoCvQ3sxsknQksJ5gIknqfnPuOvJhllheI9xqBiWdmt0uaT7CQMwS7EsR2\nyngoq1S34kbi32uzUlILgsV6X5e0Cfg84kwVKfpdezrBvZ35pVb8d9+RF7MMIOkgM/vCzOJ+hZAu\n5hLsjm3h93H3mqRJBPf0Afw7wVhlbJnZmeG3t4QbXzYHXoswUmVelvQx8DXwi/Dqd0fEmdKKTwDJ\nAKkTECT93cyGRp0pXYXjT/9DML1dwPeBX5tZnLvsCLepOSl8+E8zi+WsUUmNgJ8DhwALgIfNLLb7\nmKUKV/7ID/e9ywGamdnaqHOlCy9mGaDU0krfmG3nak7Yxdiv9PiTmfWMNlnZEjjL9RlgN8ECwwMJ\nJoBcE22qqpHUHehGyVsJHi//Fa46vJsxM1g537ual6jxp/AqoVBSczPLjzpPFXQrWn8xXFMyzrth\nFwvXRe1LUMwmEhTiGYAXsxrixSwz9Ayns4tgC4qiqe0CLKYrJiRV4safSNAsV4KrMgDMbE+CJlGc\nDfQE5prZpeF9fT6TsQZ5McsAZpYddYZMYWa/LjX+9GBcx59SJGmWa89SH8b2SfmgFucPZl+bWaGk\nPZKaEayHeWDUodKJFzPnakip8afYF4ckznKt6gczSS3NbFNt56mG98NbCR4iWDV/G/BOtJHSi08A\nca4GSXoDOCsJ40/pPMs1zkuISTqYYCaj72dWg/zKzLmalaTxp9QBp06RpagdsRpMC1f+eDPcnWC5\npBaSzjCzF6POli68mDlXs1LHn4q6PWL1izVFOs9yjdvfZ2Tq2KmZbQ5nOHoxqyFezJyrAZJ+DBxg\nZveFj2cR7GlmwG+izFYBn+Vad8q6PcN//9YgfzOdqxk3AOelPG4AHAU0AR4BnosiVEXSfJZr3K6G\n35c0FrgvfHwFwUQQV0NiezOncwnTwMxWpDyeYWZ5ZvYF0DiqUOlKUueirXYk9ZV0dThbsMipEUUr\nz1XALuCZ8GsnQUFzNcRnMzpXAyR9amaHlPPcZ2bWua4zpTNJ84A+BDuQTwTGA4eb2aAoc7noeDej\nczXjXUnDzeyh1IOSfkZCllxKmMJwBZAzgXvM7B5JsduhQNL/mtm1kl6ijEkpZjYkglhpyYuZczXj\nOuBFSecDc8JjRxHsMn1GZKnS125Jw4CLgcHhsfoR5inPE+F/74w0RQbwbkbnapCkfwMODx9+aGZv\nRpknXUnqRrAVzDtm9rSkjsC5ZnZ7xNFcRPzKzLkaFBYvL2C1r1/qjehmtkxS7Da7lLSAsu95K7r9\noUcdR0pbfmXmnEucspariuNefZI6VPS8mX1eV1nSnV+ZOecSIxwnOx/oKGlCylNNgbxoUpUvtViF\nha2LmU2RtA/++7dG+ZvpnEuSt4E1wL7AH1OObwViu3CvpOHA5UAu0Bk4AHiA+N0Pl1jezeicc7Us\nvC/uGODdoq5QSQuKds12352vAOKcSxxJZ0laIilf0hZJW1PWloyjnWa2q+iBpHrEbzHkRPNuRudc\nEt0BDDazRVEHqaJpkkYQLOjcD/hP4KWIM6UV72Z0ziWOpJlmdmLUOapKUhZwGdCfYFr+JOAv5r+A\na4wXM+dc4ki6C2hLsB/YzqLjZvZCuS+KmKTWAGa2Puos6ci7GZ1zSdQM2E5wpVPE2LsxaixIEjAS\nuJJwjoKkAoL1JEdFmS3d+JWZc87VEknXAwOBy81sWXisE3A/8JqZ/SnKfOnEi5lzLnEkdSUoCPuZ\nWXdJPYAhZnZrxNFKCFfy72dmG0odbw1MjtuKJUnmU/Odc0n0EPBfwG4AM/uAkjt9x0X90oUMisfN\n4rjKf2J5MXPOJVGOmZXeJ25PJEkqtutbPueqySeAOOeSaIOkzoQ3Hks6m2CZq7jpWc7N3AIa1XWY\ndOZjZs65xAknUTwInABsApYBF5rZ8ihzueh4MXPOJZakxkCWmW2NOouLlhcz51zihFPeS8sHZpvZ\nvLrO46Lnxcw5lziSngL6sHd9wx8RbAFzMPCcmd0RUTQXES9mzrnEkTQdGGRm28LHTYBXgB8SXJ11\nizKfq3s+Nd85l0RtSFmTkeB+s/3M7OtSx12G8Kn5zrkkehJ4V9L48PFg4KlwQshH0cVyUfFuRudc\nIknqAxRtAzPTzN6PMo+Llhcz51yiSMoGPjSzw6LO4uLDx8ycc4liZgXAJ5IOijqLiw8fM3POJVFL\n4ENJs4Cvig6a2ZDoIrkoeTFzziXRTVEHcPHiY2bOuUSS1AHoYmZTJOUA2b6sVebyMTPnXOJIGg48\nD4wLD+0PvBhdIhc1L2bOuSS6gmBa/hYAM1tCcCO1y1BezJxzSbTTzIo3t5RUj3BvM5eZvJg555Jo\nmqQRwD6S+gHPsXfRYZeBfAKIcy5xJGUBlwH9CXZtnmRmD0WbykXJi5lzLnEkXWNmd1V2zGUO72Z0\nziXRxWUcu6SuQ7j48JumnXOJIWkYcD7QUdKElKeaAnnRpHJx4MXMOZckbwNrgH2BP6Yc30qw07TL\nUD5m5pxLpFIrgOwD1PMVQDKXj5k55xKnjBVADsBXAMloXsycc0nkK4C4EryYOeeSyFcAcSV4MXPO\nJZGvAOJK8AkgzrnEKWsFEOAv5r/QMpYXM+dcIklqDWBm66PO4qLn3YzOucRQ4BZJG4BPgE8krZd0\nc9TZXLS8mDnnkuQ6glmMR5tZrpnlAscCJ0q6LtpoLkrezeicSwxJc4F+Zrah1PHWwGQz6x1NMhc1\nvzJzziVJ/dKFDIrHzepHkMfFhBcz51yS7PqWz7k0592MzrnEkFQAfFXWU0AjM/Orswzlxcw551zi\neTejc865xPNi5pxzLvG8mDnnnEs8L2bOOecSz4uZc865xPt/TedtkK28x8MAAAAASUVORK5CYII=\n",
"text/plain": [
"
"
]
},
"metadata": {
"tags": []
}
}
]
},
{
"cell_type": "code",
"metadata": {
"id": "brI-mbPaSTfE",
"colab_type": "code",
"colab": {}
},
"source": [
"#Exclude Region, Channel, Groceries and Detergents_Paper variables from df dataframe\n",
"df2=df.iloc[:,[2,3,5,7]]"
],
"execution_count": 0,
"outputs": []
},
{
"cell_type": "markdown",
"metadata": {
"id": "7WkgLZrRF2E_",
"colab_type": "text"
},
"source": [
"###To perform k-means clusterring, firstly we load the necessary libraries"
]
},
{
"cell_type": "code",
"metadata": {
"id": "c-5yDX1A0iy3",
"colab_type": "code",
"colab": {}
},
"source": [
"from sklearn.cluster import KMeans"
],
"execution_count": 0,
"outputs": []
},
{
"cell_type": "markdown",
"metadata": {
"id": "9sQKNqAFGhAg",
"colab_type": "text"
},
"source": [
"###One critical decision is the number of clusters that k-means will create. For this we use the **elbow method**.\n",
"###For this, we apply k-means and create clustering solutions for 2,3,..,10 clusters.\n",
"###Then we plot the sum of squared distances of data to their closest cluster center \n",
"###According to the elbow method the best number of clusters is where an elbow appears at the graph."
]
},
{
"cell_type": "code",
"metadata": {
"id": "GUyIWKUCIigY",
"colab_type": "code",
"outputId": "c63c742b-24f4-45c3-c707-d40a48da7097",
"colab": {
"base_uri": "https://localhost:8080/",
"height": 296
}
},
"source": [
"wcss = []\n",
"for i in range(2, 10):\n",
" kmeans = KMeans(n_clusters=i, init='k-means++', max_iter=300, n_init=10, random_state=0)\n",
" kmeans.fit(df2)\n",
" wcss.append(kmeans.inertia_)\n",
"plt.plot(range(2, 10), wcss)\n",
"plt.title('Elbow Method')\n",
"plt.xlabel('Number of clusters')\n",
"plt.ylabel('Sum of sq distances of data to their closest center')\n",
"plt.show()"
],
"execution_count": 0,
"outputs": [
{
"output_type": "display_data",
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEXCAYAAAC6baP3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzt3Xd4VGX2wPHvSSMESCAhtICERHrH\nKEWKYsF1UbG3ta1l8aeuZXVX3V3brruuuvaGYsG6spYV14YiRboBpffQEUIHCZB2fn/cGxxZktwk\nM7kzk/N5nvvM3Dt37pxB4cx9y3lFVTHGGGO8ivE7AGOMMZHFEocxxpgqscRhjDGmSixxGGOMqRJL\nHMYYY6rEEocxxpgqscRhjDGmSipMHCISKyJv1VYwxoSSiFwpIlMD9lVEjvYzpmAJ5ncRkTUicnIw\nrmWiU4WJQ1VLgLYiklBL8RhTI+4/evtF5MeA7Rm/44JDiUtF5PHDjp/lHn/N43Umicg1IQnSGA/i\nPJyTB0wTkXHAvrKDqvpYyKIypmbOUNWv/A6iHKuAC0TkDlUtdo9dASz3MSZjqsRLH8cq4L/uuY0C\nNmOiwekikici20TkERGJARCRGBH5k4isFZF8EXldRFLc18aIyO/c5xnu3cIN7n62iOwou84RbAYW\nAMPc81OBAcC4wJNEpJ+ITBeRXSIyT0ROcI8/CAwCnjnC3dTJIrLCfc+zIiKVfRf39cvc17aLyB9r\n+Odp6oBKE4eq3q+q9wOPlD13942JBmcDOUAf4Czg1+7xK93tRCALaAiU/SM9GTjBfT4E5658cMD+\nN6paWsFnvg5c7j6/CPgIOFj2oohkAJ8AfwVSgduB90UkXVX/CHwD3KiqDVX1xoDrDgeOBXoAF+Am\np4q+i4h0AZ4HLgNaAWlA6wpiN6byxCEi/UVkMbDU3e8pIs+FPDJjqu8/7q/usu3aCs79h6ruUNV1\nwBPAxe7xS4HHVDVPVX8E7gIuEpE4nMQx0L2rGAw8DBzvvm+I+3pFPgROcH/1X46TSAL9CvhUVT9V\n1VJV/RLIBU6v5LoPqeou97tMBHp5+C7nAf9V1SmqehD4M1BR0jPGU1PVEzi/XLYDqOo8fvp1ZUw4\nGqGqjQO2lyo4d33A87U4v7pxH9ce9loc0FxVV+H09/XCaTb6L7BJRDriIXGo6n6cO4o/AWmqOu2w\nU9oC5wcmP2Ag0LKi6+I0g5UpwLmzqPC7uK8d+jNQ1X24f9eNKY+XznFUdb3bXFqmJDThGFPr2gCL\n3OdHAZvc55tw/gEn4LViYIu7Pxnn13qCqm4Ukck4ndxNgO89fO7rwNfAkZp91wNvqGp5d0pVXQuh\nou/yA9C57AURScJprjKmXF7uONaLyABARSReRG4HloQ4LmNqyx0i0kRE2gA3A++6x98BbhWRdiLS\nEPgb8G7ASKjJwI3AFHd/krs/1R3GXpnJwCnA00d47U3gDBEZ5s6lShSRE0SkrO9hC05fhVcVfZf3\ngOEiMtAddv8ANjHYVMLL/yAjgRuADGAjzu35/4UyKGNq6OPD5nF8WMG5HwFzcO4SPgFedo+/AryB\nkxhWAweAmwLeNxlndGFZ4pgKJAXsV0gdE1R1xxFeW4/TUX83sBXnDuQOfvr7+iRwnojsFJGnPHxc\nud9FVRfh/P1+G+fuYyewwct3MHWXVLYCoIgcf3gb7JGOGWOMqRu8JI65qtqnsmPGGGPqhnI7x0Wk\nP87EpHQRuS3gpWQgNtSBGWOMCU8VjapKwBnOF8fPZ4rvwRlNYowxpg7y0lTVVlXXVniSMcaYOsPL\nPI56IvIikBl4vqoODVVQNdW0aVPNzMz0OwxjjIkoc+bM2aaq6ZWd5yVx/Bt4ARhNhEz8y8zMJDc3\n1+8wjDEmooiIp9YlL4mjWFWfr2E8xhhjooSXCYAfi8j/iUhLEUkt20IemTHGmLDk5Y7jCvfxjoBj\nStVKHhhjjIkSlSYOVW1XG4EYY4yJDF7W40hyVw970d1vLyLDQx+aMcaYcOSlj+NVoBBnFjk4hQ7/\nGrKIjDHGhDUviSNbVR8GigBUtQCQit9ijDEmWnlJHIUiUh938RgRySZgfeRoUlBYzEOfLWX9jgK/\nQzHGmLDlJXHcC3wOtBGRt4AJwO9DGpVPdhUU8caMNfzxPwuprBSLMcbUVZUmDlX9EjgHuBJnJbEc\nVZ0U2rD80apxfe4Y1pEpy7fy0febKn+DMcbUQV5GVZ2NM3v8E1X9L1AsIiNCH5o/LuufSa82jXng\nv4vZsa/Q73CMMSbseGqqUtXdZTuqugun+SoqxcYID53bnT37i/jrJ4v9DscYY8KOl8RxpHO8zDiP\nWJ1aJDNySDYfzN3I1BXb/A7HGGPCipfEkSsij4lItrs9BswJdWB+u3Ho0bRr2oC7P1zA/sKIKAps\njDG1wkviuAlnAuC7wL+AA8ANoQwqHCTGx/K3s7uzbkcBT0xY7nc4xhgTNrzUqtoH3FkLsYSd/tlp\nXJjThtHfrOaMHq3olpHid0jGGOM7L3ccddrdp3emSVICd32wgJJSm9thjDGWOCqRkhTPvWd0YcHG\n3bw6bbXf4RhjjO+8zOM43suxaDa8R0uGdmrGP8cvt3Ikxpg6z8sdx9Mej0UtEeEvI7oRI/AnK0di\njKnjyu0cF5H+OKXU00XktoCXkoHYUAcWbjIa1+f2YR25/+PFjJu3ibN6ZfgdkjHG+KKiO44EoCFO\ncmkUsO0Bzgt9aOHn8v6Z9GzTmAc+XsxOK0dijKmjpLJmFxFpq6pr3ecxQENV3VMbwVVXTk6O5ubm\nhuTaS37YwxlPT2VE7wwePb9nSD7DGGP8ICJzVDWnsvO89HH8XUSSRaQBsBBYLCJ31DjCCNW5ZTLX\nDc7ivTkbmLbSypEYY+oeL4mji3uHMQL4DGgHXBbSqMLcb09qT2ZaEnd/uIADRVaOxBhTt3hJHPEi\nEo+TOMapahHuaoB1VWJ8LH87pztrtxfw5IQVfodjjDG1ykviGAWsARoAU0SkLU4HeZ02ILsp5x/T\nmhen5LF4U53/4zDG1CFeVgB8SlUzVPV0dawFTqzsfSLyiojki8jCgGOPiMhSEZkvIh+KSOOA1+4S\nkZUiskxEhgUcP809tlJEwqpm1h9/2ZkmSfHc9cF8K0dijKkzvMwcby4iL4vIZ+5+F+AKD9d+DTjt\nsGNfAt1UtQewHLgr4JoXAV3d9zwnIrEiEgs8C/wC6AJc7J4bFhonJXDPGV2Zt2E3r01f43c4xhhT\nK7w0Vb0GfAG0cveXA7dU9iZVnQLsOOzYeFUtdndnAq3d52cB/1LVg6q6GlgJHOduK1U1T1ULccq6\nn+Uh5lpzRo+WnNgxnX+OX8aGnVaOxBgT/bwkjqaqOhYoBXD/4Q/GUKJf44zSAsgA1ge8tsE9Vt7x\n/yEi14lIrojkbt26NQjheVNWjgTgz1aOxBhTB3hJHPtEJA13JJWI9AN2V/yWionIH4Fi4K2aXCeQ\nqr6oqjmqmpOenh6sy3rSukkSvzu1IxOXbeXj+T/U6mcbY0xt85I4bgPGAdkiMg14HWdVwGoRkSuB\n4cCl+tPP841Am4DTWrvHyjsedq4ckEnP1ik88PEidhVYORJjTPTyMqpqLjAEp+Dhb4Cuqjq/Oh8m\nIqcBvwfOVNXADoFxwEUiUk9E2gHtgdnAt0B7EWknIgk4HejjqvPZoRYbI/z9nB7sLCjib58u8Tsc\nY4wJGS+jqs4H6qvqIpxJgO+KSB8P73sHmAF0FJENInI18AxOocQvReR7EXkBwL32WGAx8Dlwg6qW\nuP0pN+J0zi8BxrrnhqUurZK5dlAWY3M3MH2VlSMxxkQnL0UO56tqDxEZCPwFeBS4R1X71kaA1RHK\nIoeVOVBUwrAnpiDA57cMJjG+zlWgN8ZEqGAWOSwbQfVL4CVV/QSn5Lo5gsT4WP52dnfWbC/gKStH\nYoyJQl4Sx0YRGQVcCHwqIvU8vq/OOv7oppznliNZ8oOVIzHGRBcvCeACnD6GYaq6C0gF6mxZda/+\neHpnUurHc+cHC6wciTEmqngZVVUArAKGiciNQDNVHR/yyCJckwYJ3HNGF+at38XrM9b4HY4xxgSN\nl1FVN+NM1Gvmbm+KSLXncdQlZ/ZsxZAO6TzyxTI27trvdzjGGBMUXpqqrgb6quo9qnoP0A+4NrRh\nRQcR4a8juqFq5UiMMdHDS+IQfl6bqsQ9Zjxok5rE707twNdL8/lkgZUjMcZEPi+J41VglojcJyL3\n4VS1fTmkUUWZKwdk0j0jhfvGLWZ3QZHf4RhjTI146Rx/DLgKp0T6DuAqVX0i1IFFk7jYGB46tzs7\nCwqtHIkxJuLFlfeCiKQG7K5xt0OvqeqOw99jyte1VQrXDGrHqMl5jOidQf/sNL9DMsaYaqnojmMO\nkOs+lj3PDXhuquiWkzpwVGoSd3+4gANFwVjSxBhjal+5iUNV26lqlvtY9rxsP6s2g4wW9RNiefDs\nbqzeto9nvl7pdzjGGFMtXuZxnC0iKQH7jUVkRGjDil6D2qdzTp8MXpi8iqWbrRyJMSbyeBlVda+q\nHlrxzy07cm/oQop+f/plF5Lrx3Pn+1aOxBgTebwkjiOdU26nuqlcaoME7hnehe/X7+LNmWv9DscY\nY6rES+LIFZHHRCTb3R7D6SA3NXBWr1YM7pDOw58vZZOVIzHGRBAvieMmoBB4F/gXcAC4IZRB1QUi\nwoMjulGqcM9HVo7EGBM5vEwA3Keqd6pqjqoeq6p3q+q+2ggu2rVJTeK2Uzrw1ZJ8Pl2w2e9wjDHG\nE1uQyWdXHe+UI7l33CIrR2KMiQiWOHwWFxvD389xypE89LmVIzHGhL8KE4eIxIrIrbUVTF3VLSOF\nqwe2453Z65mZt93vcIwxpkIVJg5VLQEurqVY6rRbT+5Am9T63P2BlSMxxoQ3L01V00TkGREZJCJ9\nyraQR1bH1E+I5W9ndydv2z6enWjlSIwx4cvLRL5e7uMDAccUGBr8cOq2Qe3TOad3Bs9PWsXwHq3o\n2KKR3yEZY8z/8DIc98QjbJY0QuSPv+xMo8Q47vpgPqVWjsQYE4YqWo/jV6r6pojcdqTX3QWeTJCl\nNazHn4d34bax83hz1lou75/pd0jGGPMzFd1xNHAfG5WzmRA5u3cGg9o35eHPl/HDbitHYowJLxKN\npS5ycnI0Nzey15pat72AU5+YzMCj03np8mMQEb9DMsZEORGZo6o5lZ3nZT2ODiIyQUQWuvs9RORP\nwQjSlO+otCRuPbkDXy3ZwucLrRyJMSZ8eBmO+xJwF1AEoKrzgYtCGZRxXD2wHV1bJXPPuEXs3m/l\nSIwx4cFL4khS1dmHHSsORTDm5+JiY3jonB5s//EgD3221O9wjDEG8JY4tolINs7cDUTkPOCHkEZl\nDuneuqwcyTpmr97hdzjGGOMpcdwAjAI6ichG4Bbg+pBGZX7m1lM60LpJfe76YD4Hi60ciTHGX14m\nAOap6slAOtBJVQeq6pqQR2YOSUqI48Gzu7Nq6z6enbjK73CMMXVcpSVHRKQecC6QCcSVDQtV1Qcq\neJsJsiEd0hnRqxXPT1rJ8B4t6dDcptIYY/zhpanqI+AsnA7xfQGbqWV/Ht6FhvXiuOuDBVaOxBjj\nGy9FDlur6mkhj8RUKq1hPf70yy787t/zeGv2Oi7r19bvkIwxdZCXO47pItI95JEYT87pk8HAo5vy\nj8+Wsnn3Ab/DMcbUQeUmDhFZICLzgYHAXBFZJiLzA45XSEReEZH8shnn7rFUEflSRFa4j03c4yIi\nT4nISvcz+gS85wr3/BUickXNvm7kExEePLsbxaWl3PPRwsrfYIwxQVbRHcdw4AzgF8DRwKnuftnx\nyrwGHN7EdScwQVXbAxPcfdzPaO9u1wHPg5NogHuBvsBxwL1lyaYua5vWgFtO7sD4xVv4fKFNqTHG\n1K5yE4eqrlXVtcBfy54HHqvswqo6BTh8xtpZwBj3+RhgRMDx19UxE2gsIi2BYcCXqrpDVXcCX/K/\nyahOumZgO7q0TOauDxbw/fpdfodjjKlDvPRxdA3cEZFY4Jhqfl5zVS37ibwZaO4+zwDWB5y3wT1W\n3vH/ISLXiUiuiORu3bq1muFFjrjYGJ7/VR8aJsZx8Yszmbgs3++QjDF1REV9HHeJyF6gh4jscbe9\nQD7OEN0aUaeee9DGlKrqi6qao6o56enpwbpsWGub1oD3rx9AVnoDrhmTy3tzNvgdkjGmDqioqerv\nqtoIeERVk92tkaqmqepd1fy8LW4TFO5j2c/kjUCbgPNau8fKO25czRol8u5v+tM/K43b/z2P5yat\nJBrXWDHGhA8vJUeqmySOZBxQNjLqCn66cxkHXO6OruoH7HabtL4AThWRJm6n+KnuMROgYb04Xrny\nWM7q1YqHP1/G/R8vpsQmCBpjQsTLBMBqEZF3gBOApiKyAWd01EPAWBG5GlgLXOCe/ilwOrASKACu\nAlDVHSLyF+Bb97wHVNVKxB5BQlwMj1/Qi/SG9Rg9dTVb9x7knxf0JDE+1u/QjDFRxpaOjUIvTcnj\nwU+X0C8rlRcvzyE5Md7vkIwxESBoS8e6F+spIje6W8+ah2dC6drBWTxxYS/mrN3JBS/MYMsem2Fu\njAkeL2uO3wy8BTRztzdF5KZQB2ZqZkTvDF658ljW7yjgnOemszL/R79DMsZECS93HFcDfVX1HlW9\nB+gHXBvasEwwDGqfzru/6c/B4hLOf2E6c9ft9DskY0wU8JI4BAhcdq7EPWYiQLeMFN6/fgAp9eO5\n5KWZTFiyxe+QjDERzkvieBWYJSL3ich9wEzglZBGZYKqbVoD3rt+AB2aN+K6N+Yw9tv1lb/JGGPK\n4WUex2M4w2N3uNtVqvp4qAMzwdW0YT3eubYfxx/dlN+/P5+nJ6ywiYLGmGrx0jn+hqrOVdWn3O07\nEXmjNoIzwdWgXhwvX5HDOb0z+OeXy/nzRwttoqAxpsq8TAAMZpFD47P42BgePb8n6cn1GDU5j217\nC3niol42UdAY45lvRQ6Nf2JihLt+0Zk/D+/C54s2c/nLs9ldUOR3WMaYCFHbRQ5NGLl6YDueurg3\n363fyfmjpvPD7v1+h2SMiQC1XeTQhJkze7ZizFXHsWnXAc59bjortuz1OyRjTJjzVHLERLcBRzfl\n3d/0o6hUOe+FGeSusTqSxpjyWeIwAHRtlcIH1w8gtUECl46exfhFm/0OyRgTpjwnDhFpJiJHlW2h\nDMr4o01qEu+N7E+nlsmMfHMOb89a53dIxpgw5GUex5kisgJYDUwG1gCfhTgu45O0hvV459q+DO6Q\nzt0fLuCJr5bbREFjzM94ueP4C05hw+Wq2g44CafsiIlSSQlxvHR5Ducd05onvlrB3R8upLik1O+w\njDFhwssEwCJV3S4iMSISo6oTReSJkEdmfBUfG8Mj5/WgeXI9np24im0/HuTpi3vbREFjjKc7jl0i\n0hCYArwlIk8C+0IblgkHIsIdwzpx/5ld+WrJFi4dPYtdBYV+h2WM8ZmXxHEWzjrgtwKfA6uA4aEM\nyoSXKwZk8uwlfViwYTfnvTCDTbtsoqAxdZmXxHGPqpaqarGqjlHVp4A/hDowE15O796SMb8+ji27\nD3DOc9NZttkmChpTV3lJHKcc4dgvgh2ICX/9s9MYO7I/inL+C9OZlbfd75CMMT6oqMjh9SKyAOgo\nIvMDttXA/NoL0YSTzi2Tef/6ATRtVI/LXpnN5wt/8DskY0wtq+iO423gDGCc+1i2HaOqv6qF2EyY\nat0kifdHDqBrq2Suf2sub8xc63dIxphaVFF13N2qukZVL1bVtcB+QIGGNnPcNGmQwNvX9GNox2b8\n+T8L+ef4ZTZR0Jg6wsvM8TNs5rg5kvoJsYy67BguzGnD01+v5A/vz7eJgsbUAV4mAP4VZ+b4V6ra\nW0ROBKypygAQFxvDQ+d2p3lyPZ76eiXbfyzkmUv6UD/BJgoaE628jKoqUtXtwKGZ40BOiOMyEURE\nuO3UjvxlRDe+XpbPJaNnsmOfTRQ0JlrZzHETNJf1a8vzl/Zh0aY9nPfCdDbsLPA7JGNMCHidOb6f\nn88cPyOUQZnIdVq3lrx5dV+27T3IOc9NZ8kPe/wOyRgTZF6Wjt2nqiWBM8fdpitjjui4dqn8e+QA\nYkS44IUZzFhl/7sYE00qmgC4V0T2lLfVZpAm8nRs0YgP/m8AzVMSueKV2Xwy3yYKGhMtKprH0UhV\nk4EngTuBDKA1Tp0qK6tuKtWqcX3eG9mfHq1TuPGdubw2bbXfIRljgsBLH8eZqvqcqu5V1T2q+jxO\nv4cxlWqclMCb1/Tl5M7Nue/jxfzxwwUcKCrxOyxjTA14SRz7RORSEYl1F3O6FBtVZaogMT6W5y/t\nw3WDs3hr1jrOfm46eVt/9DssY0w1eUkclwAXAFvc7Xz3mDGexcXGcPfpnXnlyhw2797P8Ken8p/v\nNvodljGmGryMqlqjqmepalNVTVfVEaq6phZiM1FoaKfmfHrzILq2SuaWd7/nD+/NZ3+hNV0ZE0m8\n3HEYE1QtU+rzzrX9uOHEbMbOWc9Zz05lxRZbGMqYSOFL4hCRW0VkkYgsFJF3RCRRRNqJyCwRWSki\n74pIgntuPXd/pft6ph8xm+CKi43hjmGdGHPVcWz/sZAznpnK2Nz1VmHXmAhQ0TyOm93H44P5gSKS\nAfwWyFHVbkAscBHwD+BxVT0a2Alc7b7lamCne/xx9zwTJQZ3SOezmwfRu00Tfv/efH43dh77Dhb7\nHZYxpgIV3XFc5T4+HYLPjQPqi0gckAT8AAwF3nNfHwOMcJ+f5e7jvn6SiEgIYjI+aZacyJvX9OWW\nk9vz4fcbOeOZqVaqxJgwVlHiWOKuw3H40rELRKTaS8eq6kbgUWAdTsLYDcwBdqlq2U/NDTgTDnEf\n17vvLXbPTzv8uiJynYjkikju1q1bqxue8UlsjHDLyR1465q+7D1QzFnPTuOtWWut6cqYMFTRzPGL\ngUHASn6+dOxwalDkUESa4NxFtANaAQ2A06p7vYB4X1TVHFXNSU9Pr+nljE8GZDfls5sH0bddKn/8\ncCE3vfMdew8U+R2WMSZAhZ3jqrpZVXvi3Bk0crdN7lKy1XUysFpVt6pqEfABcDzQ2G26Aqe0Sdkg\n/41AGwD39RTAquZFsaYN6zHmquO4Y1hHPlu4meFPT2Xhxt1+h2WMcXlZOnYIsAJ4FngOWC4ig2vw\nmeuAfiKS5PZVnAQsBiYC57nnXAF85D4f5+7jvv61WvtF1IuJEW448Wj+dV0/DhaVcs5z03lt2mpr\nujImDEhlfxFFZA5wiaouc/c7AO+o6jHV/lCR+4ELgWLgO+AanL6MfwGp7rFfqepBEUkE3gB6AzuA\ni1Q1r6Lr5+TkaG5ubnXDM2Fmx75Cbv/3PL5ems+wrs15+NyepCTF+x2WMVFHROaoaqUrvHpJHPNV\ntUdlx8KJJY7oU1qqvDx1Nf/4fCktUhJ5+uLe9D6qid9hGRNVvCYOLxMAc0VktIic4G4vAfavsqlV\nMTHCtYOzGDuyP6pw/gszeGlKnjVdGeMDL4njepw+iN+622L3mDG1rs9RTfj0t4MY2qkZD366hGvG\n5LJzX6HfYRlTp1TaVBWJrKkq+qkqY6av4W+fLiWtYQJPX9ybnMxUv8MyJqIFs6nKmLAjIlx5fDve\nv34A8bExXPjiTJ6duJLS0uj7IWRMuLHEYSJa99Yp/Pe3AzmtWwse+WIZV7w6m20/HvQ7LGOiWpUS\nh7sCYHKogjGmOpIT43nm4t48eHY3Zq3ewelPfsOMVTZH1JhQ8TIB8G0RSRaRBsBCYLGI3BH60Izx\nTkS4tG9b/vN/x9OwXhyXjp7JE18tp8SarowJOi93HF1UdQ9OtdrPcGpMXRbSqIyppi6tkvn4poGc\n1SuDJ75awWUvzyJ/zwG/wzImqnhJHPEiEo+TOMa59aXsZ5wJWw3qxfHYBT15+LwezF23k9Of+oZv\nVljFZGOCxUviGAWswaliO0VE2gK2WIIJayLCBTltGHfjQJokJXD5K7N59ItlFJeU+h2aMRGvWvM4\nRCQuYO2MsGPzOEyg/YUl3DtuIWNzN3BcZipPXtyLlin1/Q7LmLATtHkcItJcRF4Wkc/c/S78VK3W\nmLBXPyGWh8/ryeMX9mThpt2c/uQ3TFya73dYxkQsL01VrwFf4Cy6BLAcuCVUARkTKmf3bs3HNw2k\neXIiV732LX//dAlF1nRlTJV5SRxNVXUsUAqHlm8tCWlUxoRIdnpD/nPD8Vza9yhGTcnjglEz2LCz\nwO+wjIkoXhLHPhFJwx1JJSL9cNb9NiYiJcbH8uDZ3Xnmkt6s2PIjpz/5DV8s2ux3WMZEDC+J4zac\nVfiyRWQa8DpwU0ijMqYWDO/Rik9+O5C2aQ34zRtzuP/jRRwstptpYyoTV9kJqjrXXT62IyDAMncu\nhzERr21aA967vj9//3Qpr05bQ+6anTxzSW/apjXwOzRjwpaXUVU3AA1VdZGqLgQaisj/hT40Y2pH\nvbhY7juzK6MuO4a12/cx/KmpfDL/B7/DMiZseWmqulZVd5XtqOpO4NrQhWSMP4Z1bcEnvx1EdrOG\n3PD2XH43dh5rt+/zOyxjwo6XxBErIlK2IyKxQELoQjLGP21Sk/j3yP6MHJLNx/M2ceKjk7jh7bks\n3GjjQYwpU+nMcRF5BGiLU3oE4DfAelX9XYhjqzabOW6CIX/PAV6etpq3Z65j78FiBh7dlOtPyGZA\ndhoBv6WMiRpeZ457SRwxOMniJPfQl8BoVQ3b4SeWOEww7TlQxFsz1/HKtNVs3XuQ7hkpjBySzWnd\nWhAbYwnERI+gJY5IZInDhMKBohI+/G4joyavYs32AjLTkrh2cBbn9mlNYnys3+EZU2PBvOM4HrgP\np7kqDmdIrqpqVhDiDAlLHCaUSkqVLxZt5oXJq5i/YTdNG9bjquMz+VW/tqTUj/c7PGOqLZiJYylw\nKzCHgFIjqhq2a3Na4jC1QVWZsWo7z09exTcrtjkrD/Y9il8PbEfz5ES/wzOmyoKZOGapat+gRVYL\nLHGY2rZw425GTcnjk/mbiIuJ4ezeGVw3JIvs9IZ+h2aMZ8FMHA8BscAHwMGy46o6t6ZBhoolDuOX\ntdv38dI3efw7dwOFJaUM69ImywmtAAAWHUlEQVSCkSdk06tNY79DM6ZSwUwcE49wWFV1aHWDCzVL\nHMZv2348yGvT1vD6jDXsOVBMv6xURg7JZkiHdBvKa8KWjaqyxGHCwI8Hi3ln1jpenrqazXsO0Lll\nMiOHZPHL7i2Ji/Uy/9aY2hPUxCEivwS6Aod6/FT1gRpFGEKWOEy4KSwu5T/fO0N5V23dR5vU+lw7\nKIvzj2lD/QQbymvCQzCbql4AkoATgdHAecBsVb06GIGGgiUOE65KS5Uvl2zhhcmr+G7dLlIbJHDl\ngEwu79+WxklWycf4K5iJY76q9gh4bAh8pqqDghVssFniMOFOVZm9egcvTF7FxGVbSUqI5eLjjuLq\nge1o1bi+3+GZOspr4qh0PQ5gv/tYICKtgO1Ay5oEZ0xdJyL0zUqjb1YaSzfvYdTkPF6bvoYx09dw\nVq8MRg7Jon3zRn6HacwReUkc/xWRxsAjwFycJWRHhzQqY+qQTi2SefzCXtx2Sgdenrqaf327jvfn\nbuDkzs0YOSSbnMxUv0M05me8NFXVU9WDZc9xOsgPlB0LR9ZUZSLZjn2FjJm+hjEz1rCroIhjM5sw\nckg2J3ZsRowVVTQhFMw+jrmq2qeyY+HEEoeJBgWFxbz77XpGf7Oajbv206F5Q34zOJsze7Ui3oby\nmhCoceIQkRZABvAmcAlOcUOAZOAFVe0UpFiDzhKHiSZFJaV8PG8ToybnsWzLXlqlJHL1oCwuOrYN\nDep5aW02xptgJI4rgCuBHOBbfkoce4HXVPWD4IQafJY4TDRSVSYuy+eFSXnMXrODxknxXN4/kysH\nZJLawIbympoLZlPVuar6ftAic67ZGKeDvRtOZ/uvgWXAu0AmsAa4QFV3usvWPgmcDhQAV1ZWJ8sS\nh4l2c9bu4PlJeXy1ZAuJ8TFcmNOGc49pTfeMFCtpYqotmMNxW4tIMs6dxktAH+BOVR1fg/ieBD5X\n1fNEJAFnguHdwARVfUhE7gTuBP4A/AJo7259gefdR2PqrGPapjL6ilRWbNnLqCl5vDVrHWNmrKVl\nSiKndmnOsK4tOK5dqpU1MSHh5Y5jnqr2FJFhwEjgT8Ab1e0cF5EU4HsgSwM+XESWASeo6g8i0hKY\npKodRWSU+/ydw88r7zPsjsPUNTv3FfLVki18sWgL36zYysHiUhonxXNSp+ac2rU5g9unW2kTU6lg\n3nGU3feeDryuqoukZvfC7YCtwKsi0hNngaibgeYByWAz0Nx9ngGsD3j/BvfYzxKHiFwHXAdw1FFH\n1SA8YyJPkwYJnJ/ThvNz2lBQWMyU5Vv5YtEWvly8mffnbqB+fCyDOzRlWNcWDO3UzMqbmBrxkjjm\niMh4nH/w7xKRRkBpDT+zD3CTqs4SkSdxmqUOUVUVkSqV7VXVF4EXwbnjqEF8xkS0pIQ4TuvWktO6\ntaSopJRZeTv4YtFmxi/ezBeLthAbI/TLSmVY1xac2qUFLVJstUJTNV6aqmKAXkCequ4SkTQgQ1Xn\nV+sDnWG+M1U1090fhJM4jsaaqowJmdJSZd6GXYxfvIUvFm0mb+s+AHq2aXyoX+ToZrZiYV0WjOG4\nnVR1qYgcsS+jJisAisg3wDWqukxE7gMauC9tD+gcT1XV37sl3W/EaSrrCzylqsdVdH1LHMZUbmX+\nXr5Y5CSR+Rt2A5Cd3sC5E+nagp6tbYRWXROMxPGSql4bihUARaQXznDcBCAPuAqIAcYCRwFrcYbj\n7nD7U54BTsMZjnuVqlaYFSxxGFM1m3bt50v3TmTW6h2UlCotkhM5tetPI7Rstnr0sxUALXEYUy27\nCgqZsCSfLxZtZsqKrRwoKiWlfjwndWrGqV1bMKSDjdCKVsG44zinojfazHFjot/+whImL9/K+MWb\nmbAkn937i0iMj2Fw+3SGdW3BSZ1thFY0CcZw3DPcx2bAAOBrd/9EYDoQtonDGBMc9RNiOa1bC07r\n1oKiklJmr3ZHaC3awvjFzgitvu3cEVpdm9MyxRahqgu8jKoaD1xRNorJHfH0mqoOq4X4qsXuOIwJ\nLVVl/obdfLFoM18s2swqd4RWj9YpDOvagmFdm3N0M1uIKtIEs1bVElXtHLAfAywKPBZuLHEYU7tW\n5v94aJ7IvPW7AMgqG6HVpTk9Wze2tUQiQDATxzM4daLecQ9dCKxU1ZtqHGWIWOIwxj+bdx/gSzeJ\nzMzbTrE7QusUd65I3ywboRWugjqqSkTOBga7u1NU9cMaxhdSljiMCQ+7C4qYsHQL4xdtYfLyrewv\nKiE5MY6hnZpxUufmDO6QTkr9eL/DNC4bjmuJw5iwsr+whG9WODW0Ji7LZ8e+QmJjhJy2TTipczOG\ndmpOdnoDm3ToI0scljiMCVslpcr363fx9dItTFiSz9LNewFom5bk3I10as5x7VJJiLMmrdpkicMS\nhzERY+Ou/Uxcms/XS/OZtnIbB4tLaVgvjkHtmzK0UzNO6NiM9Eb1/A4z6lnisMRhTETaX1jC9FXb\nmLA0n6+X5LN5zwFEoEfrxpzUqRlDOzWja6tka9IKgWCOqlqAs7zr/7yEU7OqR/VCDB1LHMZEB1Vl\n8Q97+HpJPhOW5jNvwy5UoUVyIid2asZJnZpx/NFNrQRKkAQzcTzsPn3DfbzUfXweQFXXVjfIULHE\nYUx02rr3IJOW5TNxWT5Tlm/jx4PF1IuLYUB2GkM7N2dop2ZkNLbZ69UVzMTxnar2PuzY3OouHVsb\nLHEYE/0Ki0v5ds0OJizJZ8LSLazdXgBApxaN3OG+zejVpgmxNvHQs2Amju+BG1R1mrs/AHhOVXsF\nJdIQsMRhTN2iquRt2+c2aW3h2zU7KSlVmiTFc2LHZgzt3IzBHdJJTrQ5IxUJZuI4BngFSMHp19gJ\n/LomCzmFmiUOY+q23fuLmLJ8K18vzWfSsnx2FhQRFyMcm5nqzhlpRla6rXZ4uKCPqhKRFABV3V3D\n2ELOEocxpowzZ2QnE5Y4w33L5oxkpiUxtFNzTurcjGMzbc4IBPeO42bgVWAv8BLQB7hTVccHI9BQ\nsMRhjCnPhp0FTFzqjNKavmo7he6ckcEdmjK0U3NO6JhO04Z1c85IMBPHPFXtKSLDgJHAn4A3rHPc\nGBPpCgqLmbZyO18v3cLXS/PZsucgItCrTdmckeZ0btmozswZCcZCToeu5T6eDryuqoukrvwpGmOi\nWlJCHKd0ac4pXZqjqizatIev3buRR8cv59Hxy2mZksiA7Kb0z06jX1YqrZsk+R2277zccbwKZADt\ngJ5ALDBJVY8JfXjVY3ccxpiayt97gEnLtjJpWT4zVm1nZ0ERAG1S69M/K43+2Wn0z2pKi5REnyMN\nnmA2VcUAvYA8Vd0lImlAhqrOD06owWeJwxgTTKWlyvL8vcxYtZ0Zq7Yza/UOdu93Ekm7pg3ol5VK\nv6w0+mel0Sw5chOJ1aqyxGGMCZHSUqcUysy87czMcxLJ3gPFAGSnN3CSSHYa/bLSIqqj3RKHJQ5j\nTC0pKVUWbdrNjFVOIpm9egf7CksA6NC84aGmrb7t0mjSIMHnaMtX48QhIu1UdXXQI6sFljiMMX4q\nLillwcbdzMhzmrZy1+xkf5GTSDq1aOT2jziJJCUpfGazByNxzFHVY0RkgqqeFPQIQ8gShzEmnBSV\nlDJ/wy6njyTPSSQHi0sRgS4tkw/dkRzbLtXXsijBSBzfAf8GrgceP/x1VX2spkGGiiUOY0w4O1hc\nwrz1u91Eso2563ZRWFxKjED3jBT6uf0jx2am0rCel1kTwRGMeRwXASPccxoFKzBjjKnr6sXFcly7\nVI5rl8rNtOdAUQlz1+1k5qrtzMzbwStTVzNqch6xMUKP1in0z3ISSU5mE5ISai+RlMfLcNxfqOpn\ntRRPUNgdhzEmku0vLGHO2p3MyNvGzLwdzFu/i+JSJT5W6Nm68aE+kj5tm5AYH7xFrII5jyMFuBcY\n7B6aDDwQzsUOLXEYY6LJvoPF5K7deaiPZOHG3ZSUKgmxMfQ6qvGhPpLeRzWmXlz1E0kwE8f7wEJg\njHvoMqCnqp5T7ehCzBKHMSaa7T1QRO6anYdGbS3atJtShXpxMZzatQVPX9y78oscQTBrVWWr6rkB\n+/e7izsZY4zxQaPEeE7s1IwTOzUDnPVHZq/ewcy87STGh748vJfEsV9EBqrqVAAROR7YH9qwjDHG\neJVSP/5Qscba4CVxjAReL1vICWcFwCtCF5IxxphwVmniUNV5QE8RSXb394Q8KmOMMWHL84BgSxjG\nGGMAbJFdY4wxVWKJwxhjTJVU2lQlIrHAL4HMwPNrWqvKvW4usFFVh4tIO+BfQBowB7hMVQtFpB7w\nOnAMsB24UFXX1OSzjTHGVJ+XO46PgStx/kFvFLDV1M3AkoD9fwCPq+rROCO3rnaPXw3sdI8/7p5n\njDHGJ146x1urao9gfqiItMa5i3kQuE1EBBgKXOKeMga4D3geOMt9DvAe8IyIiEbjClTGGBMBvCSO\nz0TkVFUdH8TPfQL4PT/duaQBu1S12N3fAGS4zzOA9QCqWiwiu93ztwVeUESuA65zd38UkWU1iK/p\n4dcPY5EUK0RWvJEUK0RWvJEUK0RWvDWJta2Xk7wkjpnAhyISAxQBAqiqJlcnKhEZDuSr6hwROaE6\n1zgSVX0ReDEY1xKRXC/1WsJBJMUKkRVvJMUKkRVvJMUKkRVvbcTqJXE8BvQHFgSpeeh44EwROR1I\nBJKBJ4HGIhLn3nW0Bja6528E2gAbRCQOSMHpJDfGGOMDL53j64GFwepTUNW7VLW1qmbiLBb1tape\nCkwEznNPuwL4yH0+jp9KnJznnm/9G8YY4xMvdxx5wCQR+Qw4WHYwBEvH/gH4l4j8FfgOeNk9/jLw\nhoisBHbgJJtQC0qTVy2JpFghsuKNpFghsuKNpFghsuINeaxe1uO490jHVfX+kERkjDEmrFWaOIwx\nxphAXmaOTwT+J7uo6tCQRGSMMSaseenjuD3geSJwLlBczrkRS0Ta4JQ2aY6TKF9U1Sf9jap8IpII\nTAHq4fx3fE9Vj9isGC4OLzPjdzwVEZE1wF6gBCgO56GYItIYGA10w/l/99eqOsPfqI5MRDoC7wYc\nygLuUdUnfAqpQiJyK3ANzp/rAuAqVT3gb1RHJiI3A9fiTJl4KZR/ptVqqhKR2ap6XAji8Y2ItARa\nqupcEWmEUy9rhKou9jm0I3Jn2zdQ1R9FJB6YCtysqjN9Dq1cInIbkAMkR0jiyFHVsJ/0JSJjgG9U\ndbSIJABJqrrL77gq4/6Q2Aj0VdW1fsdzOBHJwPl71UVV94vIWOBTVX3N38j+l4h0w6n1dxxQCHwO\njFTVlaH4vEqH44pIasDWVESG4cyliCqq+oOqznWf78Wpo5VR8bv8o44f3d14dwvbDquAMjOj/Y4l\nmrgrcw7GHYWoqoWRkDRcJwGrwjFpBIgD6rtzyJKATT7HU57OwCxVLXDnwk0GzgnVh3mZxzEHp3lh\nDjAD+B0/FSCMSiKSCfQGZvkbScVEJFZEvgfygS9VNZzjLSszU+p3IB4pMF5E5rjlbMJVO2Ar8KqI\nfCcio0Wkgd9BeXQR8I7fQZRHVTcCjwLrgB+A3UEuvRRMC4FBIpImIknA6TgTp0Oi0sShqu1UNct9\nbK+qp6rq1FAF5DcRaQi8D9wS7qseqmqJqvbCmWl/nHu7GnYCy8z4HUsVDFTVPsAvgBtEZLDfAZUj\nDugDPK+qvYF9wJ3+hlQ5t0ntTODffsdSHhFpglNktR3QCmggIr/yN6ojU9UlOJXDx+M0U32P0z8X\nEuUmDhE5VkRaBOxfLiIfichTIpIaqoD85PYVvA+8paof+B2PV27TxETgNL9jKUdZmZk1OO2wQ0Xk\nTX9Dqpj7axNVzQc+xGk7DkcbgA0Bd5vv4SSScPcLYK6qbvE7kAqcDKxW1a2qWgR8AAzwOaZyqerL\nqnqMqg7GWZpieag+q6I7jlE4nSy4v7Yewhl1tJvImkXpidvZ/DKwJASz4oNORNLd0TSISH3gFGCp\nv1EdWTllZsLylxuAiDRwB0jgNvucitMUEHZUdTOw3h2tBE6/QVgO6DjMxYRxM5VrHdBPRJLcfx9O\n4udrCIUVEWnmPh6F07/xdqg+q6LhuLGqusN9fiHO8NT3gffddvVoczxwGbAg4Pvdraqf+hhTRVoC\nY9yRKTHAWFX9r88xRYvmOBWhwfk78raqfu5vSBW6CXjLbf7JA67yOZ4Kucn4FOA3fsdSEVWdJSLv\nAXNxpiB8R3j/aH5fRNJwqpjfEMpBEuUOxxWRhUAvdw2MpcB1qjql7DVVDcv2dGOMMaFV0R3HO8Bk\nEdkG7Ae+ARCRo3Gaq4wxxtRBFU4AFJF+OE0i41V1n3usA9CwbM6DMcaYusWKHBpjjKkSLxMAjTHG\nmEMscRhjjKkSSxwmIomIisg/A/ZvF5H7gnTt10TkvMrPrPHnnC8iS9ylC0IWl4hkisglVY/QmCOz\nxGEi1UHgHBFp6ncggdxieF5dDVyrqieGKh5XJlClxFHF72HqGEscJlIV40zGuvXwFw7/ZS4iP7qP\nJ4jIZLd0Tp6IPCQil4rIbBFZICLZAZc5WURyRWS5W2urrKjkIyLyrYjMF5HfBFz3GxEZxxFmbYvI\nxe71F4rIP9xj9wADgZdF5JEjvOcP7nvmichDR3h9TVnSFJEcEZnkPh8iIt+723fuDPiHcArgfS8i\nt3r9Hu4M+k/cGBaKyIVe/sOY6Ge/KkwkexaYLyIPV+E9PXFKUO/AmWU9WlWPE2cRnJuAW9zzMnHq\nU2UDE935S5fjVEg9VkTqAdNEpKxaah+gm6quDvwwEWmFU3zuGJz6QeNFZISqPiAiQ4HbVTX3sPf8\nAqe4Xl9VLahibbjbcWYNTxOnYOcBnKKHt5etgSJOtd9Kv4eInAtsUtVfuu+LuuUUTPXYHYeJWG71\n4teB31bhbd+6a68cBFbhVBMFZ3W3zIDzxqpqqaquwEkwnXBqVl3ulqSZBaQB7d3zZx+eNFzHApPc\nQnnFwFs462dU5GTgVVUtcL/njkrODzQNeExEfgs0dj/zcF6/xwLgFBH5h4gMUlWb+GsASxwm8j2B\n01cQuAZFMe7/2yISAyQEvHYw4HlpwH4pP78DP3yCk+IsyXmTqvZyt3YB6zPsq9G3qLpD3xFnSWcn\nSNWHcJY6rY9zJ9HpCO/19D1UdTnOHcgC4K9u85oxljhMZHN/jY/l54uLrcFpGgJnzYf4alz6fBGJ\ncfs9soBlwBfA9eKU30dEOkjliybNBoaIs3pmLE5V2MmVvOdL4CpxFuShnKaqNfz0Hc8tOygi2aq6\nQFX/AXyLc6e0F2gU8F5P38NtZitQ1TeBR4iMcu2mFlgfh4kG/wRuDNh/CfhIRObhLGpTnbuBdTj/\n6CfjrN18QERG4zRnzRURwVl5b0RFF1HVH0TkTpz1UgT4RFU/quQ9n4tILyBXRAqBT4G7DzvtfpyO\n9b8AkwKO3yIiJ+LcQS0CPnOfl7h/Hq8BT3r8Ht2BR0SkFKfi6vUVxW3qDis5YowxpkqsqcoYY0yV\nWOIwxhhTJZY4jDHGVIklDmOMMVViicMYY0yVWOIwxhhTJZY4jDHGVMn/A7JERMYC20F/AAAAAElF\nTkSuQmCC\n",
"text/plain": [
"
"
]
},
"metadata": {
"tags": []
}
}
]
},
{
"cell_type": "markdown",
"metadata": {
"id": "xTb4qztR9_ce",
"colab_type": "text"
},
"source": [
"###Unfortunatelly, the best number of clusters is not clear from the results of the elbow method. Possible candidates are 3, 5 and 6.\n",
"###Normaly, we would test all the above by profilling the clusters and making a business driven decission.\n",
"###For our purposes, we will examine only the slution with 3 groups."
]
},
{
"cell_type": "code",
"metadata": {
"id": "1GYVKms_JAfb",
"colab_type": "code",
"outputId": "75e13ec8-ddfc-4c0b-93ca-a401825658ef",
"colab": {
"base_uri": "https://localhost:8080/",
"height": 72
}
},
"source": [
"kmeans = KMeans(n_clusters=3, init='k-means++', max_iter=300, n_init=10, random_state=0)\n",
"kmeans.fit(df2)"
],
"execution_count": 0,
"outputs": [
{
"output_type": "execute_result",
"data": {
"text/plain": [
"KMeans(algorithm='auto', copy_x=True, init='k-means++', max_iter=300,\n",
" n_clusters=3, n_init=10, n_jobs=None, precompute_distances='auto',\n",
" random_state=0, tol=0.0001, verbose=0)"
]
},
"metadata": {
"tags": []
},
"execution_count": 17
}
]
},
{
"cell_type": "markdown",
"metadata": {
"id": "FJrQap7E_aAm",
"colab_type": "text"
},
"source": [
"###In ordder to best understand the clusters, we append a column indicating the cluster of each record to the non-ormalised data set."
]
},
{
"cell_type": "code",
"metadata": {
"id": "c5xfa9GiA1SE",
"colab_type": "code",
"outputId": "72aeb87f-6696-4370-8dc7-1d0c7d11cd7e",
"colab": {
"base_uri": "https://localhost:8080/",
"height": 226
}
},
"source": [
"dat['cluster']=kmeans.labels_\n",
"dat.head()"
],
"execution_count": 0,
"outputs": [
{
"output_type": "execute_result",
"data": {
"text/html": [
"