{ "metadata": { "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.0" }, "orig_nbformat": 2, "kernelspec": { "name": "python37364bit8d3b438fb5fc4430a93ac2cb74d693a7", "display_name": "Python 3.7.0 64-bit ('3.7')" }, "metadata": { "interpreter": { "hash": "70b38d7a306a849643e446cd70466270a13445e5987dfa1344ef2b127438fa4d" } } }, "nbformat": 4, "nbformat_minor": 2, "cells": [ { "source": [ "## Linear and Polynomial Regression for Pumpkin Pricing - Lesson 3\n", "\n", "Load up required libraries and dataset. Convert the data to a dataframe containing a subset of the data: \n", "\n", "- Only get pumpkins priced by the bushel\n", "- Convert the date to a month\n", "- Calculate the price to be an average of high and low prices\n", "- Convert the price to reflect the pricing by bushel quantity" ], "cell_type": "markdown", "metadata": {} }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "output_type": "execute_result", "data": { "text/plain": [ " City Name Type Package Variety Sub Variety Grade Date \\\n", "0 BALTIMORE NaN 24 inch bins NaN NaN NaN 4/29/17 \n", "1 BALTIMORE NaN 24 inch bins NaN NaN NaN 5/6/17 \n", "2 BALTIMORE NaN 24 inch bins HOWDEN TYPE NaN NaN 9/24/16 \n", "3 BALTIMORE NaN 24 inch bins HOWDEN TYPE NaN NaN 9/24/16 \n", "4 BALTIMORE NaN 24 inch bins HOWDEN TYPE NaN NaN 11/5/16 \n", "\n", " Low Price High Price Mostly Low ... Unit of Sale Quality Condition \\\n", "0 270.0 280.0 270.0 ... NaN NaN NaN \n", "1 270.0 280.0 270.0 ... NaN NaN NaN \n", "2 160.0 160.0 160.0 ... NaN NaN NaN \n", "3 160.0 160.0 160.0 ... NaN NaN NaN \n", "4 90.0 100.0 90.0 ... NaN NaN NaN \n", "\n", " Appearance Storage Crop Repack Trans Mode Unnamed: 24 Unnamed: 25 \n", "0 NaN NaN NaN E NaN NaN NaN \n", "1 NaN NaN NaN E NaN NaN NaN \n", "2 NaN NaN NaN N NaN NaN NaN \n", "3 NaN NaN NaN N NaN NaN NaN \n", "4 NaN NaN NaN N NaN NaN NaN \n", "\n", "[5 rows x 26 columns]" ], "text/html": "
\n\n\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n
City NameTypePackageVarietySub VarietyGradeDateLow PriceHigh PriceMostly Low...Unit of SaleQualityConditionAppearanceStorageCropRepackTrans ModeUnnamed: 24Unnamed: 25
0BALTIMORENaN24 inch binsNaNNaNNaN4/29/17270.0280.0270.0...NaNNaNNaNNaNNaNNaNENaNNaNNaN
1BALTIMORENaN24 inch binsNaNNaNNaN5/6/17270.0280.0270.0...NaNNaNNaNNaNNaNNaNENaNNaNNaN
2BALTIMORENaN24 inch binsHOWDEN TYPENaNNaN9/24/16160.0160.0160.0...NaNNaNNaNNaNNaNNaNNNaNNaNNaN
3BALTIMORENaN24 inch binsHOWDEN TYPENaNNaN9/24/16160.0160.0160.0...NaNNaNNaNNaNNaNNaNNNaNNaNNaN
4BALTIMORENaN24 inch binsHOWDEN TYPENaNNaN11/5/1690.0100.090.0...NaNNaNNaNNaNNaNNaNNNaNNaNNaN
\n

5 rows × 26 columns

\n
" }, "metadata": {}, "execution_count": 2 } ], "source": [ "import pandas as pd\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "\n", "pumpkins = pd.read_csv('../../data/US-pumpkins.csv')\n", "\n", "pumpkins.head()\n" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "output_type": "execute_result", "data": { "text/plain": [ " Month Variety City Package Low Price High Price Price\n", "70 1 3 1 0 5 3 13.636364\n", "71 1 3 1 0 10 7 16.363636\n", "72 2 3 1 0 10 7 16.363636\n", "73 2 3 1 0 9 6 15.454545\n", "74 2 3 1 0 5 3 13.636364" ], "text/html": "
\n\n\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n
MonthVarietyCityPackageLow PriceHigh PricePrice
7013105313.636364
71131010716.363636
72231010716.363636
7323109615.454545
7423105313.636364
\n
" }, "metadata": {}, "execution_count": 3 } ], "source": [ "from sklearn.preprocessing import LabelEncoder\n", "\n", "pumpkins = pumpkins[pumpkins['Package'].str.contains('bushel', case=True, regex=True)]\n", "\n", "new_columns = ['Package', 'Variety', 'City Name', 'Month', 'Low Price', 'High Price', 'Date']\n", "\n", "pumpkins = pumpkins.drop([c for c in pumpkins.columns if c not in new_columns], axis=1)\n", "\n", "price = (pumpkins['Low Price'] + pumpkins['High Price']) / 2\n", "\n", "month = pd.DatetimeIndex(pumpkins['Date']).month\n", "\n", "new_pumpkins = pd.DataFrame({'Month': month, 'Variety': pumpkins['Variety'], 'City': pumpkins['City Name'], 'Package': pumpkins['Package'], 'Low Price': pumpkins['Low Price'],'High Price': pumpkins['High Price'], 'Price': price})\n", "\n", "new_pumpkins.loc[new_pumpkins['Package'].str.contains('1 1/9'), 'Price'] = price/1.1\n", "\n", "new_pumpkins.loc[new_pumpkins['Package'].str.contains('1/2'), 'Price'] = price*2\n", "\n", "new_pumpkins.iloc[:, 0:-1] = new_pumpkins.iloc[:, 0:-1].apply(LabelEncoder().fit_transform)\n", "\n", "new_pumpkins.head()\n" ] }, { "source": [ "A scatterplot reminds us that we only have month data from August through December. We probably need more data to be able to draw conclusions in a linear fashion." ], "cell_type": "markdown", "metadata": {} }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "output_type": "execute_result", "data": { "text/plain": [ "" ] }, "metadata": {}, "execution_count": 4 }, { "output_type": "display_data", "data": { "text/plain": "
", "image/svg+xml": "\n\n\n\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n\n", "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXAAAAD4CAYAAAD1jb0+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAahUlEQVR4nO3df5Bd9Xnf8fdHawFrTCLLLIrQDwtjjZjWqiV7g8So08o4qhjw4A11W6vg4o5rJW08xUMqGwVNYndgIKME3DYZZ8B2rBSC7bHVhQKpqjHSZOIBOStLZsGgymB+LQraGMtgsgaxevrHPStWd+/VnrN77z0/7uc1s7P3PPdczjMHnWfP/Z7vD0UEZmZWPnPyTsDMzGbGBdzMrKRcwM3MSsoF3MyspFzAzcxK6m2dPNi5554by5Yt6+QhzcxKb//+/X8fEX318Y4W8GXLljE0NNTJQ5qZlZ6kZxvF3YRiZlZSLuBmZiXlAm5mVlIu4GZmJeUCbmZWUh3thWKdMXhghO27DvHisTHOn9fLlo0rGFi9KO+0zKzFXMArZvDACFt3DjN2fByAkWNjbN05DOAiblYxbkKpmO27Dp0s3hPGjo+zfdehnDIys3ZxAa+YF4+NZYqbWXm5gFfM+fN6M8XNrLxcwCtmy8YVzJ2jU2Jz54gtG1fklJGZtYsLeBVpmm0zq4RUBVzSM5KGJR2UNJTE5kvaLelw8vud7U3V0ti+6xDHx09d5/T4ePghplkFZbkD/1BErIqI/mT7BuC7EbEc+G6ybTnzQ0yz7jGbJpSPAjuS1zuAgdmnY7Plh5hm3SNtAQ/g/0raL2lzElsQEUeS138HLGh5dpbZlo0r6J3bc0qsd26PH2KaVVDakZj/NCJGJJ0H7Jb05OQ3IyIkRaMPJgV/M8DSpUtnlaxNb2K0pYfSm1WfIhrW3eYfkL4A/AL4NLA+Io5IWgjsjYjT3ub19/eHV+QxM8tG0v5Jzx9PmrYJRdLZks6ZeA38C+Ax4D7g2mS3a4F7W5eumZlNJ00TygLgf0ma2P8vI+L/SPpb4FuSPgU8C/zr9qVpZmb1pi3gEfE08P4G8Z8CH25HUmZmNj2PxDQzKynPB15BV9/5MN976uWT2+sunM/dn74kx4zMrB18B14x9cUb4HtPvczVdz6cU0Zm1i4u4BVTX7yni5tZebmAm5mVlAu4mVlJuYBXzLoL52eKm1l5uYBXzL/qX0rdgjzMUS1uZtXiAl4x23cd4kTd9DYnAi/oYFZB7gdeMV7QIZsNt+3l8NHXTm4vP+9sdl+/Pr+EzDLwHXjFeEGH9OqLN8Dho6+x4ba9+SRklpELeMVs2biiYRu4F3SYqr54Txc3KxoX8IoZevblhm3gQ896II9Z1biAV8w9+57PFDez8nIBr5jxJissNYt3s+XnnZ0pblY0LuAV0yNlinez3devn1Ks3QvFyiR1N0JJPcAQMBIRH5H0deCfAz9PdvlkRBxsfYqWxaY1S7jrkecaxm0qF+tsBg+MeMHsAsnSD/w64AngVybFtkTEt1ubks3GTQMrgVqb93gEPRKb1iw5GTebqcEDI2zdOczY8XEARo6NsXXnMICLeE5SFXBJi4ErgJuB69uakc3aTQMrXbCt5bbvOnSyeE8YOz7O9l2HXMBzkrYN/EvA54ATdfGbJT0q6XZJZzb6oKTNkoYkDY2Ojs4mVzPLkUf5Fs+0BVzSR4CjEbG/7q2twEXArwPzgc83+nxE3BER/RHR39fXN9t8zSwnHuVbPGnuwNcBV0p6BvgGcKmkuyLiSNS8Dvw5cHEb8zSznG3ZuILeuT2nxHrn9niUb46mLeARsTUiFkfEMuDjwEMRcY2khQCSBAwAj7U1UzPL1cDqRdxy1UoWzetFwKJ5vdxy1Uq3f+doNrMR3i2pDxBwEPjt1qRkZkU1sHqRC3aBZCrgEbEX2Ju8vrQN+ZiZWUoeiWlmVlIu4GZmJeUCbmZWUi7gZmYl5QJuZlZSLuBmZiXlVekryFN+prfshgemxJ659YocMimHbYPDnumyQHwHXjETU36OHBsjeGvKz8EDI3mnVjiNivfp4t1u2+Awdz3y3MnVncYjuOuR59g2OJxzZt3LBbxiTjflp9lseL3V4nEBrxhP+Wnt4vVWi8cFvGI85ae1i9dbLR4X8IrxlJ/WLs3WVfV6q/lxAa8YT/mZXrPeJu6F0thNAyu5Zu3Sk3fcPRLXrF3qXig5UnSw/aq/vz+GhoY6djwzsyqQtD8i+uvjvgM3Myup1AN5JPUAQ8BIRHxE0gXUllh7F7Af+EREvNGeNC0LD+RJ74IbHmDyd1ABP3ETSlNrbt7NS6++dZkvOOcM9t24IceMuluWO/DrgCcmbf8hcHtEvBf4GfCpViZmM+OBPOnVF2+ASOI2VX3xBnjp1TdYc/PunDKyVAVc0mLgCuArybaAS4FvJ7vsoLYupuXMA3nSa/b0x72aG6sv3tPFrf3S3oF/CfgccCLZfhdwLCLeTLZfABp+R5e0WdKQpKHR0dFZJWvT80Aes+4xbQGX9BHgaETsn8kBIuKOiOiPiP6+vr6Z/CcsAw/kMeseae7A1wFXSnqG2kPLS4H/BsyTNPEQdDHgRtYC8ECe9JqNH/S4wsYWnHNGpri137QFPCK2RsTiiFgGfBx4KCKuBvYAH0t2uxa4t21ZWmoeyJPeT269Ykqxdi+U5vbduGFKsXYvlHxlGsgjaT3wX5JuhO+hdkc+HzgAXBMRr5/u8x7IY2aWXbOBPJkWdIiIvcDe5PXTwMWtSM7MzLLzSEwzs5IqxZJqHlmYzUU3Psgvx99qGjurRzx58+U5ZlRcXiLMyqzwd+AeWZhNffEG+OV4cNGND+aUUXF5iTAru8IXcI8szKa+eE8X72ZeIszKrvAF3CMLrV28RJiVXeELuEcWWrt4iTAru8IXcI8szOasnsbFp1m8m3mJMCu7whdwjyzM5smbL59SrN0LpTEvEWZlV4puhAOrF7lgZ+BinV7/u+ez58lRXjw2xq/96ln0v3t+3imZpVaKAm7WDhNdVCd6OU10UQV8w2ClUPgmFLN2cRdVKzvfgVfQsgZLgj3jGfamcBfV7Dwqulh8B14xjYr36eLdzEuqZeNR0cXjAm5mqbjJqXhcwM0sFTc5FU+aNTHPkvR9ST+U9LikLybxr0v6iaSDyc+q9qdrZnnxqOjiSXMH/jpwaUS8H1gFXCZpbfLelohYlfwcbFuWZpa7LRtX0DPn1EFiPXPkUdE5SrMmZkTEL5LNucmPn/MU1LoLGw9EaRY3S2vo2ZcZP3HqpT9+Ihh69uWcMrJUbeCSeiQdBI4CuyNiX/LWzZIelXS7pDPblqWl9r2nGl9MzeJmaXn63eJJVcAjYjwiVgGLgYslvQ/YClwE/Dq1hY0/3+izkjZLGpI0NDo62qK0zazTPP1u8WTqhRIRx4A9wGURcSRpXnkd+HOaLHAcEXdERH9E9Pf19c0+YzPLhaffLZ40vVD6JM1LXvcCG4AnJS1MYgIGgMfamaiZ5ct34MWTZij9QmCHpB5qBf9bEXG/pIck9QECDgK/3cY8zcyszrQFPCIeBVY3iF/alozMzCwVj8Q0MyspF3Azs5JyAa+YXzmzJ1PcLK0F55yRKW7t5wJeMa+8Pp4pbpbWS6++kSlu7ecCbmZWUi7gZmYl5QJeMW4Dt3ZxG3jxuIBXzJVN1idsFu9mX/o3jaewbxbvdvtu3DClWC845wz23bghp4zMixpXzOlmjLtpYGWHsym2ZkuBbd91yAv1NuFiXSy+A68Yz1eRnpcIs7JzAa+YZvPCeb64qbxEmJWdm1Aq5u1n9PDaG1P7fL/9DD/ErLdl4wo++82pKwF6ibDmtg0Oc8++5xmPoEdi05olbprLke/AK6ZR8T5dvJv9/uBwpni32zY4zF2PPHeyOW48grseeY5tPl+5cQG3ruVRq9l4SbXicQE3s1T8gLx4XMDNLBUvqVY8aZZUO0vS9yX9UNLjkr6YxC+QtE/SjyV9U5KHYxWAe6Gk51Gr2WxasyRT3NovzR3468ClEfF+YBVwmaS1wB8Ct0fEe4GfAZ9qX5qWVrMvs/6SO5XbwLO5aWAl16xdevKOu0fimrVL3QslR2mWVAvgF8nm3OQngEuBf5vEdwBfAL7c+hTNrChuGljpgl0gqdrAJfVIOggcBXYDTwHHIuLNZJcXgIZjjyVtljQkaWh0dLQVOZuZGSkLeESMR8QqYDFwMXBR2gNExB0R0R8R/X19fTNM08zM6mXqhRIRx4A9wCXAPEkTTTCLgZEW52ZmZqeRphdKn6R5yeteYAPwBLVC/rFkt2uBe9uVpKX3zK1XZIp3M58rK7s0c6EsBHZI6qFW8L8VEfdL+hHwDUk3AQeAr7YxT0up2bDmbYPDfvjUgIu1lVmaXiiPAqsbxJ+m1h5uBeL5wM26h0diVoyHO5t1D08nWzE9UsNi7eHOjS274YEpMTerWFn4Drxi3tP39kzxbtaoeJ8ublY0LuAV8/ToP2SKm1l5uYBXjNvAzbqHC7iZWUm5gJuZlZQLeMV40v30PBLTys7dCCtm05ol3PXIcw3jNpWLtZWZC3jFTIy2vGff84xH0COxac0Sj8I0qyBFB3sn9Pf3x9DQUMeOZ2ZWBZL2R0R/fdxt4GZmJeUmlAraNjjsJhSzLuACXjHbBodPeYg5HnFy20XcrFrchFIxp5tO1syqJc2KPEsk7ZH0I0mPS7ouiX9B0oikg8nP5e1P16bjofRm3SNNE8qbwO9GxA8knQPsl7Q7ee/2iPij9qVnWc0RnGhQq+d4HI9Z5aRZkecIcCR5/aqkJ4BF7U7MZubMt81h7PiJhnEzq5ZMV7WkZdSWV9uXhD4j6VFJX5P0ziaf2SxpSNLQ6OjorJK16f2yQfE+XdzMyit1AZf0DuA7wGcj4hXgy8CFwCpqd+h/3OhzEXFHRPRHRH9fX18LUrbTOX9eb6a4mZVXqgIuaS614n13ROwEiIiXImI8Ik4Ad+IFjgthy8YVmeJmVl5peqEI+CrwRETcNim+cNJuvwk81vr0LKs/3XM4U9zMyitNL5R1wCeAYUkHk9jvAZskrQICeAb4rbZkaJkcPvpapriZlVeaXih/AzTqhPZg69MxM7O03LfMzKykXMArZvl5Z2eKm1l5uYBXzO7r108p1svPO5vd16/PJyEzaxvPRlhBLtZm3cF34GZmJeUCbmZWUm5CqSCvyGPWHVzAK8Yr8ph1DzehVIxX5DHrHi7gFeMVecy6hwt4xfSo8dI7zeJmVl4u4BWzac2STHEzKy8/xKyYiQeV7oViVn2KDraN9vf3x9DQUMeOZzadDbftPWWqXU87cHo+X9kMHhhh+65DvHhsjPPn9bJl4woGVmdfUljS/ojor4+7CcW6Vn0xgtq86Rtu25tPQgXn85XN4IERtu4cZuTYGAGMHBtj685hBg+MtOwYaVbkWSJpj6QfSXpc0nVJfL6k3ZIOJ78bLmpsVlRe/CIbn69stu86xNjx8VNiY8fH2b7rUMuOkaYN/E3gdyPiB5LOAfZL2g18EvhuRNwq6QbgBuDzLcvMZuzqOx/me0+9fHJ73YXzufvTl+SYkVn3efHYWKb4TEx7Bx4RRyLiB8nrV4EngEXAR4EdyW47gIGWZWUzVl+8Ab731MtcfefDOWVk1p3On9ebKT4TmdrAJS0DVgP7gAURcSR56++ABS3LymasvnhPF+9mXvwiG5+vbLZsXEHv3J5TYr1ze9iycUXLjpG6gEt6B/Ad4LMR8crk96LWlaVhdxZJmyUNSRoaHR2dVbJmreTFL7Lx+cpmYPUibrlqJYvm9SJg0bxebrlq5Yx6oTSTqh+4pLnUivfdEbEzCb8kaWFEHJG0EDja6LMRcQdwB9S6EbYgZ7OWeWXs+Gm37VS/86Hlp3SL+50PLc87pUIbWL2opQW7XppeKAK+CjwREbdNeus+4Nrk9bXAva1Pz7J6W5MR883i3WzNzbt56dU3Tom99OobrLl5d04ZFVsnusVZNmmaUNYBnwAulXQw+bkcuBXYIOkw8BvJtuXszSbfcZrFu1l98Z4u3u060S3Ospm2CSUi/gZodv/24damY2ZF1YlucZZNKeZCadVwVDObufPn9TLSoFi3slucZVP4ofRud7N2WXDOGZni3a4T3eIsm8IXcLe7Wbvsu3HDlGK94Jwz2HfjhpwyKrZOdIuzbArfhOJ2N2snF+ts2t0tzrIp/B14J4ajmpmVUeEL+JaNK5g759ROMHPnyO1uTZzV07jDULO4mZVX4ZtQgKmdGF2Lmvrgsnc2nPfkg8s8228jnrkxm22Dw17tqUAKfwe+fdchjo+fOgrl+Hj4IWYTnswqPc/cmM22wWHueuQ5xpNVvMYjuOuR59g2OJxzZt2r8AXcDzGtXfzHLpt79j2fKW7tV/gC7oeYZsUw3mT93GZxa7/CF3APHjArhh41fvjULG7tV/gC7sEDZsWwac2STHFrv1L0QvHgAbP8TfQ2cS+U4ihFATezYrhpYKULdoG4gJtZastueGBK7Jlbr8ghk3LYcNteDh997eR2q5egK3wbuJkVQ6Pifbp4t6sv3gCHj77Ghtv2tuwYaZZU+5qko5IemxT7gqSRuhV6zMwsUV+8p4vPRJo78K8DlzWI3x4Rq5KfB1uWkZmZpTJtAY+IvwY8NM3MrGBm0wb+GUmPJk0sTWdKkrRZ0pCkodHR0VkcztLwbIRmxbD8vLMzxWdipgX8y8CFwCrgCPDHzXaMiDsioj8i+vv6+mZ4OEvrl+ONhzU3i3ezRU2mY2gW73bNepu4F0pju69fP+XG6awe5d8LJSJeiojxiDgB3Alc3LKMzDqk0QK9p4t3u/dubdzbpFm8211958NTbpx+OR4tne1yRgVc0sJJm78JPNZsXzOrhjebfIlrFu92nZjtctqBPJLuAdYD50p6AfgDYL2kVUAAzwC/1bKMzMwslWkLeERsahD+ahtyMTOzDDwSs2Ka9TVxHxSzzlp34fxM8ZlwAa+YZs2RbqY066y7P33JlGLd6jVXPZlVxfRIDVdI8aT7Zp3X7gWyfQdeMZ50P71ODLQwaycX8Iq5aWAl16xdevKOu0fimrVLPYdzA7uvXz+lWLd6us8q8R+84lF0cEHS/v7+GBoa6tjxzKy12j2/tTUmaX9E9NfH3QZuXW3wwAjbdx3ixWNjnD+vly0bV3j5vtNwsS4WF3DrWoMHRti6c5ix4+NAbQj91p3DAC7iVgpuA7eutX3XoZPFe8LY8XG27zqUU0Zm2biAW9d6scmkVc3iZkXjAm5d6/wm08Y2i5sVjQu4da0tG1dkilvtucG6Wx/ighseYN2tDzF4YCTvlLqaC7h1rT/dczhTvNtNPPQdOTZG8NZDXxfx/LiAW9fqxKrhVeKHvsXjAm5mqfihb/FMW8CTRYuPSnpsUmy+pN2SDie/my5qbGbV4Ie+xZPmDvzrwGV1sRuA70bEcuC7ybZZqXhuj2y2bFxB79yeU2K9c3v80DdH0xbwiPhroH4Rt48CO5LXO4CBFudl1naezCqbgdWLuOWqlSya14uARfN6ueWqlR61mqNUk1lJWgbcHxHvS7aPRcS85LWAn01sN/jsZmAzwNKlSz/47LPPtiZzM7Mu0Wwyq1k/xIzaX4CmfwUi4o6I6I+I/r6+vtkezszMEjMt4C9JWgiQ/D7aupTMzCyNmRbw+4Brk9fXAve2Jh0zM0srTTfCe4CHgRWSXpD0KeBWYIOkw8BvJNtmZtZB084HHhGbmrz14RbnYmZmGXR0STVJo8BsuqGcC/x9i9JpJeeVXhFzAueVlfPKZrZ5vTsipvQC6WgBny1JQ4260uTNeaVXxJzAeWXlvLJpV16eC8XMrKRcwM3MSqpsBfyOvBNownmlV8ScwHll5byyaUtepWoDNzOzt5TtDtzMzBIu4GZmJVXIAi7pMkmHJP1Y0pS5xiWdKembyfv7ktkS887pk5JGJR1Mfv5Du3NKjjtlwY269yXpvyd5PyrpAwXJa72kn086X7/fgZyWSNoj6UeSHpd0XYN9On6+UuaVx/k6S9L3Jf0wyeuLDfbJ41pMk1de12OPpAOS7m/wXuvPVUQU6gfoAZ4C3gOcAfwQ+Ed1+/wn4M+S1x8HvlmAnD4J/EkO5+ufAR8AHmvy/uXAXwEC1gL7CpLXempTFHfyXC0EPpC8Pgf4fw3+P3b8fKXMK4/zJeAdyeu5wD5gbd0+Hb0WM+SV1/V4PfCXjf5fteNcFfEO/GLgxxHxdES8AXyD2gISk01eUOLbwIeTecnzzCkX0XjBjck+CvxF1DwCzJuYSTLnvDouIo5ExA+S168CTwD1qxF0/HylzKvjknPwi2RzbvJT3+uh09di2rw6TtJi4ArgK012afm5KmIBXwQ8P2n7Bab+Yz65T0S8CfwceFfOOQH8y+Rr97clLWljPlmkzT0PlyRfg/9K0j/u5IGTr6+rqd29TZbr+TpNXpDD+UqaBA5SmzJ6d0Q0PV8duhbT5gWdvx6/BHwOONHk/ZafqyIW8LL638CyiPgnwG7e+ktrjf2A2vwO7wf+BzDYqQNLegfwHeCzEfFKp447nWnyyuV8RcR4RKwCFgMXS3pfJ447nRR5dfR6lPQR4GhE7G/nceoVsYCPAJP/Wi5OYg33kfQ24FeBn+aZU0T8NCJeTza/AnywjflkkeZ8dlxEvDLxNTgiHgTmSjq33ceVNJdakbw7InY22CWX8zVdXnmdr0nHPwbsYeoC552+FlPllcP1uA64UtIz1JpYL5V0V90+LT9XRSzgfwssl3SBpDOoNfbfV7fP5AUlPgY8FMmTgbxyqmsnvZJaO2YR3Af8u6R3xVrg5xFxJO+kJP3aRPufpIup/Vts64WfHO+rwBMRcVuT3Tp+vtLkldP56pM0sfZtL7ABeLJut05fi6ny6vT1GBFbI2JxRCyjVh8eiohr6nZr+bmadj7wTouINyV9BthFrffH1yLicUn/FRiKiPuo/WP/n5J+TO1B2ccLkNN/lnQl8GaS0yfbmdME1RbcWA+cK+kF4A+oPdQhIv4MeJBaz4ofA/8A/PuC5PUx4D9KehMYAz7e7guf2l3SJ4DhpP0U4PeApZPyyuN8pckrj/O1ENghqYfaH4xvRcT9eV6LGfLK5Xqs1+5z5aH0ZmYlVcQmFDMzS8EF3MyspFzAzcxKygXczKykXMDNzErKBdzMrKRcwM3MSur/A6xFhQLt6CFWAAAAAElFTkSuQmCC\n" }, "metadata": { "needs_background": "light" } } ], "source": [ "import matplotlib.pyplot as plt\n", "plt.scatter('Month','Price',data=new_pumpkins)" ] }, { "source": [ "Try some different correlations" ], "cell_type": "markdown", "metadata": {} }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "output_type": "stream", "name": "stdout", "text": [ "0.3236397181608923\n0.6061712937226016\n" ] } ], "source": [ "print(new_pumpkins['City'].corr(new_pumpkins['Price']))\n", "\n", "print(new_pumpkins['Package'].corr(new_pumpkins['Price']))\n" ] }, { "source": [ "Drop unused columns" ], "cell_type": "markdown", "metadata": {} }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "output_type": "stream", "name": "stdout", "text": [ "\nInt64Index: 415 entries, 70 to 1742\nData columns (total 7 columns):\n # Column Non-Null Count Dtype \n--- ------ -------------- ----- \n 0 Month 415 non-null int64 \n 1 Variety 415 non-null int64 \n 2 City 415 non-null int64 \n 3 Package 415 non-null int64 \n 4 Low Price 415 non-null int64 \n 5 High Price 415 non-null int64 \n 6 Price 415 non-null float64\ndtypes: float64(1), int64(6)\nmemory usage: 25.9 KB\n" ] } ], "source": [ "\n", "new_pumpkins.dropna(inplace=True)\n", "new_pumpkins.info()\n", "\n", "\n" ] }, { "source": [ "Create a new dataframe" ], "cell_type": "markdown", "metadata": {} }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "output_type": "execute_result", "data": { "text/plain": [ " Package Price\n", "70 0 13.636364\n", "71 0 16.363636\n", "72 0 16.363636\n", "73 0 15.454545\n", "74 0 13.636364\n", "... ... ...\n", "1738 2 30.000000\n", "1739 2 28.750000\n", "1740 2 25.750000\n", "1741 2 24.000000\n", "1742 2 24.000000\n", "\n", "[415 rows x 2 columns]" ], "text/html": "
\n\n\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n
PackagePrice
70013.636364
71016.363636
72016.363636
73015.454545
74013.636364
.........
1738230.000000
1739228.750000
1740225.750000
1741224.000000
1742224.000000
\n

415 rows × 2 columns

\n
" }, "metadata": {}, "execution_count": 7 } ], "source": [ "new_columns = ['Package', 'Price']\n", "lin_pumpkins = new_pumpkins.drop([c for c in new_pumpkins.columns if c not in new_columns], axis='columns')\n", "\n", "lin_pumpkins\n" ] }, { "source": [ "Set X and y arrays to correspond to Package and Price" ], "cell_type": "markdown", "metadata": {} }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "X = lin_pumpkins.values[:, :1]\n", "y = lin_pumpkins.values[:, 1:2]\n" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "output_type": "stream", "name": "stdout", "text": [ "Model Accuracy: 0.3315342327998989\n" ] } ], "source": [ "from sklearn.linear_model import LinearRegression\n", "from sklearn.metrics import r2_score, mean_squared_error, mean_absolute_error\n", "from sklearn.model_selection import train_test_split\n", "\n", "\n", "X_train, X_test, y_train, y_test = train_test_split(X, y, test_size=0.2, random_state=0)\n", "lin_reg = LinearRegression()\n", "lin_reg.fit(X_train,y_train)\n", "\n", "pred = lin_reg.predict(X_test)\n", "\n", "accuracy_score = lin_reg.score(X_train,y_train)\n", "print('Model Accuracy: ', accuracy_score)" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "output_type": "display_data", "data": { "text/plain": "
", "image/svg+xml": "\n\n\n\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n\n", "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEGCAYAAABiq/5QAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAffUlEQVR4nO3de5BkZZnn8e+vq1qlhOFaMr1gVSIYS3Tj0ErJ4rgq0o3DKKE4ugZGobigKV2wq+HgtWJXd2NqZowYRdaJrjHlzqQDDurIos5INzCAoyzV0HJp3BW1q4Ro6eYuFgFW9rN/5Kmu6uqsqjzZefL6+0RkdOab5+R56kA+9dbznvO+igjMzKx7rGh2AGZm1lhO/GZmXcaJ38ysyzjxm5l1GSd+M7Mu09vsAKpxxBFHRC6Xa3YYZmZtZcuWLY9HRP/C9rZI/LlcjomJiWaHYWbWViRNVmp3qcfMrMs48ZuZdRknfjOzLuPEb2bWZZz4zcy6jBO/mWWuWCySy+VYsWIFuVyOYrHY7JC6Wltczmlm7atYLJLP55mengZgcnKSfD4PwPDwcDND61ru8ZtZpkZHR/ck/VnT09OMjo42KSLLtMcvaTvwW6AEzETEkKTDgOuBHLAdeF9EPJVlHGbWPFNTU6naLXuN6PG/NSLWRsRQ8vozwOaIeDWwOXltZh1qYGAgVbtlrxmlnncBVyfPrwbOakIMZtYgY2Nj9PX17dXW19fH2NhYkyKyrBN/AD+UtEVSPmk7MiJ2JM9/AxxZaUdJeUkTkiZ27dqVcZhmlpXh4WEKhQKDg4NIYnBwkEKh4IHdJlKWa+5KOioiHpX0CuBm4L8AN0bEIfO2eSoiDl3qc4aGhsKTtJmZpSNpy7wy+x6Z9vgj4tHk353Ad4CTgcckrUqCWgXszDIGMzPbW2aJX9LLJR00+xx4G/AAcCNwbrLZucB3s4rBzMz2leXlnEcC35E0e5xvRMQ/S7ob+Kak84FJ4H0ZxmBmZgtklvgj4pfAiRXanwDWZXVcMzNbmu/cNTPrMk78ZjXwpGPWzjxJm1lKnnTM2p17/GYpedIxa3dO/GYpedIxa3dO/GYpedIxa3dO/GYpedIxa3dO/GYpedIxa3dO/GY1GB4eZvv27ezevZvt27c76S9jZGSE3t5eJNHb28vIyEizQ+pqvpzTzDI1MjLC+Pj4ntelUmnP640bNzYrrK6W6bTM9eJpmc3aV29vL6VSaZ/2np4eZmZmmhBR92jKtMxmZpWS/lLtlj0nfjPLVE9PT6p2y54Tv5llanY6i2rbLXse3DWzTM0O4BYKBUqlEj09PeTzeQ/sNpEHd83MOlTTBncl9Ui6V9JNyeurJP1K0tbksTbrGMzMbE4jSj0fAx4C/mBe2ycj4oYGHNvMzBbItMcv6WjgHcBlWR7HzMyql3Wp5yvAp4DdC9rHJN0n6RJJL620o6S8pAlJE7t27co4TDOz7pFZ4pd0JrAzIrYseOuzwPHA64HDgE9X2j8iChExFBFD/f39WYVpZtZ1suzxvxF4p6TtwHXAaZL+PiJ2RNkLwJXAyRnGYGZmC2SW+CPisxFxdETkgLOBWyLiHEmrACQJOAt4IKsYzMxsX824gasoqR8QsBW4oAkxmJl1rYYk/oi4DbgteX5aI45pZmaVea4eM7Mu48RvZtZlnPjNzLqME7+ZWZdx4jcz6zJO/GZmXcaJ36wGa9asQdKex5o1a5odklnVnPjNUlqzZg3btm3bq23btm1O/tY2nPjNUlqY9JdrN2s1TvxmZl3Gid/MrMs48ZultHr16lTtZq3Gid8spQcffHCfJL969WoefPDBJkVklk4zpmU2a3tO8tbO3OM3M+syTvxmZl0m88QvqUfSvZJuSl4fI+kuSQ9Lul7SS7KOwazeRkZG6O3tRRK9vb2MjIw0OySzqjWix/8x4KF5r78IXBIRxwFPAec3IAazuhkZGWF8fJxSqQRAqVRifHzcyd/aRqaJX9LRwDuAy5LXAk4Dbkg2uZrygutmbaNQKKRqN2s1Wff4vwJ8CtidvD4ceDoiZpLXjwBHVdpRUl7ShKSJXbt2ZRymWfVme/rVtpu1mswSv6QzgZ0RsaWW/SOiEBFDETHU399f5+jMatfT05Oq3azVZNnjfyPwTknbgesol3guBQ6RNHv/wNHAoxnGYFZ3+Xw+VbtZq8ks8UfEZyPi6IjIAWcDt0TEMHAr8N5ks3OB72YVg1kWNm7cyIYNG/b08Ht6etiwYQMbN25scmTW7oaGQNr7sXlz/Y+jiKj/py48iHQqcHFEnCnpVZT/AjgMuBc4JyJeWGr/oaGhmJiYyDxOM7NGk5bfptY0LWlLRAwtbG/IlA0RcRtwW/L8l8DJjTiumVmr+d3v4MADq9/+5pvrH4Pn6jEzy9idd8Kb3pR+v6wKMp6ywawGvnPXlnPhhXN1+mqT/qteVU72s4+suMdvltLsnbuzZu/cBTzA2+Wqqdcv9NWvwkUX1T+WpTRkcHd/eXDXWklvb2/Fm7V6enqYmZmpsId1qt27oZbbNx5+GI49tv7xLNTUwV2zTuI7d7vb5CTkcun3m5mp7ZdEFlzjN0vJd+52n8sum6vXp0n68+v1rfS/hxO/WUq+c7c7nHjiXLL/yEeq2+f88xszOLu/XOoxS2l2ALdQKFAqlejp6SGfz3tgtwPUMjh7yy3w1rfWP5YsucdvVoONGzcyMzNDRDAzM+Okv4xisUgul2PFihXkcjmKxWKzQwLguef2nh6hWs8+O9erb7ekD+7xm1nGisUi+Xye6elpACYnJ/eUxYaHhxsezx13wJvfnH6/Vi7dpOUev5llanR0dE/SnzU9Pc3o6GjDYhgZmevVV5v0jz++Per1tXCP36wGqlAXaId7YpphamoqVXu91FKvHx+HCy6ofyytxonfLKVKSX+23cl/XwMDA0xOTlZsr6dab6b6xS/KUyV0E5d6zCxTY2Nj9PX17dXW19fH2NjYfn/29u1zJZw0Sb9UmivhdFvSByd+M8vY8PAwhUKBwcFBJDE4OEihUKh5YPfrX59L9sccU/1+8+v1K7o882U2V4+klwG3Ay+lXFK6ISI+L+kq4C3AM8mmH4qIrUt9lufqsVayWKkHXOfPymteAw88kG6fj3wECoVs4mkXzZir5wXgtIh4TtJK4E5JP0je+2RE3JDhsc2szdUyOHvbbfCWt9Q9lI6T5Zq7ERHPJS9XJg93h6ztrVu3LlW7VaceN1M56Vcn00qXpB5JW4GdwM0RcVfy1pik+yRdIumlWcZgVm+bF1n9erF2W9y//utcoj/ooOr3m1+vT7OflWWa+COiFBFrgaOBkyWdAHwWOB54PeUF1z9daV9JeUkTkiZ27dqVZZhm1kAf/ehcsj/11Or2Wb26c2+maoaGjG1HxNPArcAZEbEjKQO9AFzJIguvR0QhIoYiYqi/v78RYZpZRuaXcKodcP3a1+YS/YMPZhtft8ks8Uvql3RI8vwA4HTgZ5JWJW0CzgJSjtWbNZdr/MsrlWqr1//qV3PJ3rNcZyfLq3pWAVdL6qH8C+abEXGTpFsk9QMCtgJdcIO0dZIdO3akau8W99wDJ52Ufr9SydfVN1pmiT8i7gNeW6H9tKyOadYI27ZtS9XeyT78Ybj88vT7uU7fXJ6rx8xSqeX6+qOOgkceqX8sVhv/gWVmy6qlXn/ttXP1eif91uLEb5bS6tWrU7W3oyeeqC3ZP/nkXLI/55zs4rP948RvltKqVatStbeLq66aS/RHHFH9fvOvrz/00MzCszpyjd8spU66c/fII2HnzvT7eXC2vbnHb9Zl5pdwqk36F1zgO2c7Saoev6S+iJhefkszaxWlEvTW8Lf91q1w4on1j8ear6oev6Q/lrQN+Fny+kRJGzONzKxFtcOdu1u2zPXq0yT9+StTOel3rmpLPZcAfwI8ARARPwWqXKverLNs2rRpnyS/bt06Nm3a1KSIys47by7ZD+2z9MbivDJV96m6LxARv16w8lCp/uGYtYdmJ/lZtdxMdcwx8Mtf1j8Wax/V/n7/taQ/BkLSSkkXAw9lGJeZLaKW6+uLxblevZO+VdvjvwC4FDgKeBT4IXBhVkGZ2ZzHH4daZiZ/8klfV2+VVdXjj4jHI2I4Io6MiFdExDkR8UTWwZm1qpGREXp7e5FEb28vIyMjdf38K6+c69WnSfq+mcqqUe1VPVfPzq2fvD5U0hXZhWXWukZGRhgfH6dUKg9zlUolxsfH9zv59/fPJfvzzqt+P19fb2lVW+P/o2QVLQAi4ikqTLls1g0KiywhtVj7UubX6x9/vLp9LrrIyd72T7U1/hWSDk0SPpIOS7GvWUeZ7elX2z7fzAysXJn+mPfdB695Tfr9zCqpNnl/CfixpH+kvHLWe4GxpXaQ9DLgduClyXFuiIjPSzoGuA44HNgCfCAiXqwxfrOG6+npqZjke3p6Km4/MQGvf33643hlKstKtYO71wB/BjwG/Ab4s4i4dpndXgBOi4gTgbXAGZJOAb4IXBIRxwFPAefXGrxZM+QXWQx2fvu5586VcNIkfd9MZQDFYpFcLseKFSvI5XIUi8W6fr5iiSKhpD+IiGeT0s4+IuLJqg4i9QF3AhuA7wF/GBEzkt4AfCEi/mSp/YeGhmJiYqKaQ5ll7tBDD+Xpp5+u8E76gvtxx8HPf77/MVnnKBaL5PN5pqfnpkXr6+ujUCgwPDyc6rMkbYmIfe7jXq5P8Y3k3y3AxLzH7OvlDtojaSuwE7gZ+AXwdETMJJs8QvneALO2sXfSj3mP6lx33Vyv3knfFhodHd0r6QNMT08zOjpat2MsWeOPiDNVnqfhLRExlfbDI6IErE0uBf0OcHy1+0rKA3mAgYGBtIc2y8SWLVBLz/7pp+Hgg+sejnWgqanKqXax9losW0WMci3oe/tzkORS0FuBNwCHSJr9hXM05TuBK+1TiIihiBjqr+W2RbM6Ofjg/Z/8zEnfqrVYR7eeHeBqh4/ukZTqugRJ/bM3fUk6ADid8vw+t1K+KgjgXOC7aT7XrBHmX1//7LPV7+fr621/jY2N0dfXt1dbX18fY2NLXkiZSrWJ/z8AP5H0C0n3Sbpf0n3L7LMKuDXZ7m7g5oi4Cfg08AlJD1O+pPPyWoM3q6daJj9L9uSQQw51sre6GB4eplAoMDg4iCQGBwdrGthdypJX9ezZSBqs1B4Rk3WLZAm+qsey8PzzsKBjVZXx8fJShGatbrGrepYc3E1uwroAOA64H7h83hU5Zm1ndBT+8i/T7zczA4vcn2XWdpa7c/dq4PfAHcCfAquBj2UdlFk91bJYCbhOb51ruRr/6mQK5q9RHpB9UwNiMttvtdbrPThrrWD9+vVI2vNYv359XT9/ucT/+9knLvFYq6sl2Z90kpO9tZb169ezefPmvdo2b95c1+S/3JQNJeB3sy+BA4Dp5HlExB/ULZIleHDXKrn7bjj55PT7PfwwHHts/eMxqwct0XOp5mKcBZ+VfnA3IjycZS2lt7c8a2Va7s2bzfGc+tbyPDhrVl+e+NVakgdnrVutW7cuVXstnPitJfzud7Ul+8suc7K3zrJp06Z9kvy6devYtGlT3Y7hUo81zWc+A1/8Yvr9fDOVdbp6JvlKnPitoVyvN2s+J37LnJO9WWtxjd8yUUu9/uSTXa83awQnfquLf/u32pL99u1zif6uuzILz8zmcanHauYSjll7cuK3VJzszdpfZqUeSa+UdKukbZIelPSxpP0Lkh6VtDV5vD2rGKw+fDOVWWfJssY/A/x5RKwGTgEulLQ6ee+SiFibPL6fYQxWg+eeqy3ZX3GFk71ZO8is1BMRO4AdyfPfSnoIOCqr49n+ufhi+NKX0u9XKsEKXyJg1lYaUuOXlANeC9wFvBG4SNIHgQnKfxU81Yg4bG+u15t1p8z7apIOBL4FfDwingXGgWOBtZT/IqjYz5SUlzQhaWLXrl1Zh9k1XK83s0wTv6SVlJN+MSK+DRARj0VEKSJ2A18HKi6lERGFiBiKiKH+/v4sw+x4tST7oSEne7NOleVVPQIuBx6KiC/Pa181b7N3Aw9kFUO3+tGPakv2k5Nzif7uu7OLrxMUi0VyuRwrVqwgl8tRLBabHZJZ1bKs8b8R+ABwv6StSdvngPdLWgsEsB34aIYxdA3X6xunWCySz+eZnp4GYHJyknw+D8Dw8HAzQzOrypJr7rYKr7lbmZN9c+RyOSYnJ/dpHxwcZPv27Y0PyGwRNa25a63Hyb75pqamUrWbtRpfgd3inn22tnr9Ndd4cDYrAwMDqdrNWo0Tfwu68MK5RH/wwdXvVyrNJfoPfCC7+Lrd2NgYfX19e7X19fUxNjbWpIjM0nGpp0W4hNM+ZgdwR0dHmZqaYmBggLGxMQ/sWtvw4G4TOdmbWZYWG9x1qafBaqnXv+UtrtebWf048WfstttqS/ZTU3OJ/rbbsorOzLqRa/wZcAnHzFqZE3+dONmbWbtw4t8PTvZm1o5c40/hmWdqq9dfe60HZ82sdbjHv4yRERgfT7+fV6Yys1blxF+BSzhm1smc+BNO9mbWLTq6GDEyMkJvby+S6O3tZWRkZK/3a6nXn3GG6/Vm1t46NvGPjIwwPj5OqVQCoFQqMT7+cE3J/te/nkv0P/hBRgGbmTVIZqUeSa8ErgGOpLzaViEiLpV0GHA9kKO8Atf7IuKpeh+/UCgkz2rrlrs3b2adKsse/wzw5xGxGjgFuFDSauAzwOaIeDWwOXldd+Wefrrs7RKOmXWDzBJ/ROyIiHuS578FHgKOAt4FXJ1sdjVwVhbH7+npqWo7J3sz6zYNqfFLygGvBe4CjoyIHclbv6FcCqq0T17ShKSJXbt2pT7m7OLXC51++mVO9mbW1TJP/JIOBL4FfDwinp3/XpQXA6iYfiOiEBFDETHU39+f+rgbN25kw4YRenp6AdHT08uGDSP88IcfruGnMDPrHJkmfkkrKSf9YkR8O2l+TNKq5P1VwM6sjr9x40ZmZmaICGZmZti4cWNWhzIzq5tisUgul2PFihXkcjmKxWJdPz+zxC9JwOXAQxHx5Xlv3Qicmzw/F/huVjGYmbWbYrFIPp9ncnKSiGBycpJ8Pl/X5J/Z0ouS/iNwB3A/sDtp/hzlOv83gQFgkvLlnE8u9VmduvSimdlCuVyOycnJfdoHBwfZvn17qs9abOnFzK7jj4g7gcVuk1qX1XHNzNrZ1NRUqvZadOydu2Zm7WhgYCBVey06OvFnPUBiZlZvY2Nj9PX17dXW19fH2NhY3Y7RsYm/EQMkZmb1Njw8TKFQYHBwEEkMDg5SKBQYHh6u2zEyG9ytp1oGd+s5QGJm1o4WG9zt2B5/IwZIzMzaUccm/kYMkJiZtaOOTfyNGCAxM8tC296522yNGCAxM6u3tr5zt558566ZdYtG3LnbsT1+M7N25Dt399OaNWuQtOexZs2aZofU0pZbnN7Msuc7d/fDmjVr2LZt215t27Ztc/JfROXF6ced/M0arBEXpnRsjb88K3Rl7fAzN1pvb++epD9fT08PMzMzTYjIrHsVi0VGR0eZmppiYGCAsbGxmi5MWazG78RvgM+XWSfy4K4tabHF6atdtN7M2ocTvwGLL06/WLuZta8sl168QtJOSQ/Ma/uCpEclbU0eb8/q+JZOeXH6DXt6+D09PWzYsMHrFJt1oCyXXnwz8BxwTUSckLR9AXguIv4mzWe5xm9mll7Da/wRcTuw5Fq6ZmbWeM2o8V8k6b6kFHToYhtJykuakDSxa9euRsZnZtbRGp34x4FjgbXADuBLi20YEYWIGIqIof7+/tQHevnLX56q3cysWzQ08UfEYxFRiojdwNeBk7M61vPPP5+q3WD9+vV7TXGxfv36ZodkZhloaOKXtGrey3cDDyy27f7avXt3qvZut379ejZv3rxX2+bNm538zTpQb1YfLOkfgFOBIyQ9AnweOFXSWiCA7cBHszq+pbMw6S/XbmbtK7PEHxHvr9B8eVbHMzOz6vjOXTOzLuPEbwCsW7cuVbuZtS8nfgNg06ZN+yT5devWsWnTpiZFZGZZyazG32yHH344TzzxRMV2q8xJ3qw7dGyP/9JLL2XlypV7ta1cuZJLL720SRGZmbWGjk38w8PDXHnllQwODiKJwcFBrrzyyppWsTEz6yQduwKXmVm368oVuIrFIrlcjhUrVpDL5SgWi80Oycys6Tp2cLdYLJLP55mengZgcnJyz2pSLveYWTfr2B7/6OjonqQ/a3p6mtHR0SZFZGbWGjo28U9NTaVqNzPrFh2b+AcGBlK1m5l1i45N/GNjYxWv4x8bG2tSRGZmraFjEz/su+D6Uguwm5l1i45N/KOjo7z44ot7tb344ose3DWzrtexid+Du2ZmlWWW+CVdIWmnpAfmtR0m6WZJP0/+PTSr43tw18yssix7/FcBZyxo+wywOSJeDWxOXmdibGyMvr6+vdr6+vo8uGtmXS+zxB8RtwNPLmh+F3B18vxq4Kysjj88PEyhUNhrkrZCoeC7ds2s62U6SZukHHBTRJyQvH46Ig5Jngt4avZ1hX3zQB5gYGDgpMnJycziNDPrRC03SVuUf+Ms+lsnIgoRMRQRQ/39/Q2MzMysszU68T8maRVA8u/OBh/fzKzrNTrx3wicmzw/F/hug49vZtb1sryc8x+AHwP/XtIjks4H/ho4XdLPgfXJazMza6DM5uOPiPcv8ta6rI5pZmbLa4ulFyXtAvbnsp4jgMfrFE49Oa7qtWJM4LjSclzp7G9cgxGxz9UxbZH495ekiUqXNDWb46peK8YEjistx5VOVnF17Fw9ZmZWmRO/mVmX6ZbEX2h2AItwXNVrxZjAcaXluNLJJK6uqPGbmdmcbunxm5lZwonfzKzLdFTil3SGpP8r6WFJ+8z1L+mlkq5P3r8rmT202TF9SNIuSVuTx4ezjik57j4L5Sx4X5L+VxL3fZJe1yJxnSrpmXnn6783IKZXSrpV0jZJD0r6WIVtGn6+qoyrGefrZZL+j6SfJnH9jwrbNOO7WE1czfo+9ki6V9JNFd6r/7mKiI54AD3AL4BXAS8BfgqsXrDNCPB3yfOzgetbIKYPAX/bhPP1ZuB1wAOLvP924AeAgFOAu1okrlMpT/XdyHO1Cnhd8vwg4P9V+O/Y8PNVZVzNOF8CDkyerwTuAk5ZsE1Dv4sp4mrW9/ETwDcq/bfK4lx1Uo//ZODhiPhlRLwIXEd54Zf55i8EcwOwLlkXoJkxNUVUXihnvncB10TZT4BDZmdWbXJcDRcROyLinuT5b4GHgKMWbNbw81VlXA2XnIPnkpcrk8fCq0ga/V2sNq6Gk3Q08A7gskU2qfu56qTEfxTw63mvH2HfL8GebSJiBngGOLzJMQG8JykP3CDplRnGk0a1sTfDG5I/138gaU0jD5z8mf1ayr3F+Zp6vpaIC5pwvpLSxVbKU6/fHBGLnq8GfRerjQsa/338CvApYPci79f9XHVS4m9X/xvIRcQfATcz95vdKruH8vwjJwJfBf6pUQeWdCDwLeDjEfFso467nGXiasr5iohSRKwFjgZOlnRCI467nCriauj3UdKZwM6I2JLlcRbqpMT/KDD/t/PRSVvFbST1AgcDTzQzpoh4IiJeSF5eBpyUYTxpVHM+Gy4inp39cz0ivg+slHRE1seVtJJyci1GxLcrbNKU87VcXM06X/OO/zRwK3DGgrca/V2sKq4mfB/fCLxT0nbKpeDTJP39gm3qfq46KfHfDbxa0jGSXkJ5EOTGBdvMXwjmvcAtkYyYNCumBXXgd1Ku07aCG4EPJlernAI8ExE7mh2UpD+crW9KOpny/8OZJozkeJcDD0XElxfZrOHnq5q4mnS++iXNrq19AHA68LMFmzX6u1hVXI3+PkbEZyPi6IjIUc4Pt0TEOQs2q/u5ymw+/kaLiBlJFwH/Qvlqmisi4kFJ/xOYiIgbKX9JrpX0MOUBxLNbIKb/KumdwEwS04eyjGmWygvlnAocIekR4POUB7uIiL8Dvk/5SpWHgWngP7dIXO8FNkiaAZ4Hzs46YVDulX0AuD+pDwN8DhiYF1czzlc1cTXjfK0CrpbUQ/kXzTcj4qZmfhdTxNWU7+NCWZ8rT9lgZtZlOqnUY2ZmVXDiNzPrMk78ZmZdxonfzKzLOPGbmXUZJ37rKpJKyayLD0j6R0l9NXzGFyRdnEV8Zo3gxG/d5vmIWBsRJwAvAhc0OyCzRnPit252B3AcgKR/krQlmac9P7uByusp3JNMcrZ54QdI+kgy+dkByfO7k22/NfvXhKRjJf1E0v2S/kLSc/P2/2Syz32qMD+8WRac+K0rJXOe/Clwf9J0XkScBAxRvnvzcEn9wNeB9ySTnP2nBZ9xEXAmcFZEPA98OyJen2z7EHB+sumlwKUR8RrKs3bO7v824NWUp+9eC5wk6c3Z/MRmczpmygazKh0wb3qDOyjfDg/lZP/u5PkrKSfkfuD2iPgVQETMXyfgg5Snyj0rIn6ftJ0g6S+AQ4ADKU/VAfAG4Kzk+TeAv0mevy153Ju8PjA57u37+0OaLcWJ37rN88m0vHtIOhVYD7whIqYl3Qa8bJnPuZ9yL/1o4FdJ21WUfxH8VNKHKM85tBQBfxURX0sRv9l+c6nHrDzN7VNJ0j+e8tKJAD8B3izpGABJh83b517go8CNkv5d0nYQsCOZKnl43rY/Ad6TPJ8/wda/AOcl8+kj6ShJr6jjz2VWkRO/Gfwz0CvpIeCvKSdqImIXkAe+LemnwPXzd4qIO4GLge8lc9z/N8orYP2Ivaf7/TjwCUn3UR5MfibZ/4eUSz8/lnQ/5WX1DsrqhzSb5dk5zTKWXN3zfESEpLOB90dES6y9bN3JNX6z7J0E/G2yIMrTwHlNjse6nHv8ZmZdxjV+M7Mu48RvZtZlnPjNzLqME7+ZWZdx4jcz6zL/H1RP0V64LOVqAAAAAElFTkSuQmCC\n" }, "metadata": { "needs_background": "light" } } ], "source": [ "\n", "plt.scatter(X_test, y_test, color='black')\n", "plt.plot(X_test, pred, color='blue', linewidth=3)\n", "\n", "plt.xlabel('Package')\n", "plt.ylabel('Price')\n", "\n", "plt.show()\n", "\n" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "output_type": "execute_result", "data": { "text/plain": [ "array([[33.15655975]])" ] }, "metadata": {}, "execution_count": 12 } ], "source": [ "lin_reg.predict( np.array([ [2.75] ]) )" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "output_type": "execute_result", "data": { "text/plain": [ " Month Variety City Package Price\n", "70 1 3 1 0 13.636364\n", "71 1 3 1 0 16.363636\n", "72 2 3 1 0 16.363636\n", "73 2 3 1 0 15.454545\n", "74 2 3 1 0 13.636364\n", "... ... ... ... ... ...\n", "1738 1 1 9 2 30.000000\n", "1739 1 1 9 2 28.750000\n", "1740 1 1 9 2 25.750000\n", "1741 1 1 9 2 24.000000\n", "1742 1 1 9 2 24.000000\n", "\n", "[415 rows x 5 columns]" ], "text/html": "
\n\n\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n
MonthVarietyCityPackagePrice
70131013.636364
71131016.363636
72231016.363636
73231015.454545
74231013.636364
..................
1738119230.000000
1739119228.750000
1740119225.750000
1741119224.000000
1742119224.000000
\n

415 rows × 5 columns

\n
" }, "metadata": {}, "execution_count": 13 } ], "source": [ "new_columns = ['Variety', 'Package', 'City', 'Month', 'Price']\n", "poly_pumpkins = new_pumpkins.drop([c for c in new_pumpkins.columns if c not in new_columns], axis='columns')\n", "\n", "poly_pumpkins" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [ { "output_type": "execute_result", "data": { "text/plain": [ "" ], "text/html": "\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n
Month Variety City Package Price
Month1.0000000.171330-0.188728-0.144847-0.148783
Variety0.1713301.000000-0.248441-0.614855-0.863479
City-0.188728-0.2484411.0000000.3016040.323640
Package-0.144847-0.6148550.3016041.0000000.606171
Price-0.148783-0.8634790.3236400.6061711.000000
" }, "metadata": {}, "execution_count": 14 } ], "source": [ "corr = poly_pumpkins.corr()\n", "corr.style.background_gradient(cmap='coolwarm')" ] }, { "source": [ "Select the Package/Price columns" ], "cell_type": "markdown", "metadata": {} }, { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [], "source": [ "X=poly_pumpkins.iloc[:,3:4].values\n", "y=poly_pumpkins.iloc[:,4:5].values\n" ] }, { "source": [ "Create Polynomial Regression model" ], "cell_type": "markdown", "metadata": {} }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [ { "output_type": "display_data", "data": { "text/plain": "
", "image/svg+xml": "\n\n\n\n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n\n", "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEGCAYAAABiq/5QAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO3deXxU1f3/8dcnCVsUCwpaiCXRilawLhVRq1VL61J3rTutC2qUYN2KG6lbK1+xRcXaokalxRrcKK61tVapWn9uuFAFW8Ul1CCLCm5BIMn5/XHucJOQkJkk997JzPv5eMwj95zMzP1kYD5z56zmnENERPJHQdIBiIhIvJT4RUTyjBK/iEieUeIXEckzSvwiInmmKOkA0jFgwABXVlaWdBgiIt3Kyy+//JFzbmDL+m6R+MvKypgzZ07SYYiIdCtmVtNavZp6RETyjBK/iEieUeIXEckzSvwiInlGiV9EJM8o8YuI5BklfhGRPKPELyKSZ5T4RSRS1dXVlJWVUVBQQFlZGdXV1UmHlPe6xcxdEemeqqurKS8vp66uDoCamhrKy8sBGD16dJKh5TVd8YtIZCorK9cm/ZS6ujoqKysTikhAiV9EIrRw4cKM6iUeSvwiEpkhQ4ZkVC/xUOIXkcgceOCBGdVLPJT4RSQyjz76aEb1Eo9IR/WY2fvA50ADUO+cG2FmGwP3AGXA+8AxzrnlUcYhIslQG392iuOK//vOuR2dcyOC8sXAE865ocATQVlEcpDa+LNTEk09hwHTg+PpwOEJxCAiMZg4cSLFxcXN6oqLi5k4cWJCEQlEn/gd8Hcze9nMyoO6zZxzHwbHi4HNWnugmZWb2Rwzm7Ns2bKIwxSRKIwePZqqqipKS0sxM0pLS6mqqtLkrYSZcy66Jzcrcc7VmtmmwOPAz4CHnHP9mtxnuXOu//qeZ8SIEU577oqIZMbMXm7SzL5WpFf8zrna4OdS4H5gJLDEzAYFQQ0ClkYZg4iINBdZ4jezDcysb+oY2A94A3gIOCm420nAg1HFICIi64pyOOdmwP1mljrPDOfc38zsJeBeMzsVqAGOiTAGERFpIbLE75x7F9ihlfqPgR9EdV4REVk/zdwV6YCKigqKioowM4qKiqioqEg6pKyl9fizj9bjF8lQRUUFN91009pyQ0PD2vLUqVOTCisraT3+7BTpcM6uouGckk2KiopoaGhYp76wsJD6+voEIspeZWVl1NTUrFNfWlrK+++/H39AeSaR4Zwiuai1pL+++nymtXqykxK/SIYKCwszqs9nfk2e3sApwKgW9ZIUJX6RDKXaqNOtz2fjx1+H2TPANPyajIfQo0cPrdWTMHXuimQo1YFbVVVFQ0MDhYWFlJeXq2O3hZdfhssuOwDnmi7SdjOwU1IhSUCduyLS5WbOhBNPhJUrW/vtNEpLf6nO3Rioc1dEIuccXHUVHH1006S/HJjU5F5jqKnZOv7gZC019YhIl1i5Ek49Fe66q2ntW8Ahwc9vAkcH9bfy+efQt2/MQQqgK36RDtHM3eYWL4bvf7950h81CmA3fNIHOAv4ODgu5ZJL4oxQmlLiF8lQauZuatx+auZuvib/116DkSPhhRfCujPPhL/9DXwzT8pS4Jy1pd//Hp55JqYgpRl17opkSDN3Qw88AKNHQ7AiAwUFcMMNMG4cmMHw4cOZP39+i0c9DBwMwNChMHcu9OkTa9h5Q527Il1EM3d9J+6kSXDEEWHS32gjePRROOssn/QB5s2bx7Bhw5o9dujQ69hoI3/89ttw+eUxBi6AEr9IxvJ95u5XX8FJJ9GsjX7LLeH552H//de9/4QJE5rtuXv55acyeXL4+2uvhZdeij5uCSnxi2Qon2fuLlniO23/9Kewbu+94cUXYdtt171/anXOmpoanHNrV+fs06c66PyFxkY45RRYtSqev0EA51zW33beeWcnkk3Gjh3rCgsLHeAKCwvd2LFjkw4pcnPnOjdkiHO+ocffTjvNuVWr2n5MaWmpA9a5lZaWunfeca64OHyuyy6L72/JF8Ac10pOVeeuiLTr4YfhhBPgiy98uaAAJk+Gc88N2/NbU1BQQGs5xsxobGzkhhv8cwAUFfllHrbfPoI/IE+pc1dEMuYc/OY3cNhhYdLv2xceegjOO2/9SR/aXoUzVX/WWfDd7/q6+noYM8b/lGgp8YtIq1at8on4wgv9BwDAFlvAc8/BQQel9xwTJ06kuLi4WV1xcfHa1TkLC+H226FXL/+7l1/2nb0SLSV+EVnHsmWw777wxz+GdXvu6SdpDR+e/vOMHj2aqqqqZqN6qqqqmm27+K1vNR/Sefnl8N//dv5vkLapjV9EmnnjDTjkEGi6eObJJ8PNN4dX5l1tzRrYbTd45RVf3mMPePpp35cgHac2fhFp16OP+jb3VNI3g1//GqZNiy7pA/To4c9RFCwb+eyzfkkHiYYSv4jgHFx/vb/S//xzX7fhhvDgg3DBBe134naFHXZoPins4ovhvfeiP28+UuIXyXOrV8MZZ8D55/vJVABDhvir7kMOiTeWykpIrfBQVwennx52LEvXUeIXyWMffwz77Qe33hrW7b67n4mbxHj6Xr18k0+qbf+JJ/yoH+laSvwieerNN2HXXeGpp8K6n/4UnnwSNtssubh23dXPEUj5+c+htja5eHKREr9IHnrsMT+K5p13wrqrr4bp06F37+TiSvnlL2GrrfzxZ5/59f3V5NN1lPhF8ohzcOONcOCBPqECFBfDrFm+MzWOTtx0FBc3b+J55JGWWzpKZyjxi3SAma1zy3Zr1kBFBZx9dtiJu/nm8K9/+XX1s81ee8HYsWH57LNh6dLk4sklSvwiGWoryWdz8v/kEzjgAD8JK2XXXX0n7k47JRdXe665xo8wAt8R/bOfJRtPrlDiF8lxb73l2/OffDKsO/54mD0bBg1KLq509O0LVVVh+d574f77k4snVyjxi+Swf/zDX9m//XZY96tfQXV199nndv/9/ZIRKRUVsHx5m3eXNCjxi+Som27yzTsrVvhynz5w333wi19kTyduuq67Dr7+dX+8eLGfbCYdp8QvkmPq631beEUFpPZ/HzwYnnkGjjoq2dg6qn9//0GW8sc/+iGp0jFK/CIZamtF22xY6XbFCr9W/u9+F9aNGOE3M9955+Ti6gqHHw7HHBOWy8vDdYUkM0r8Ih3Q2j6mSVuwwHfi/v3vYd3RR/uZuYMHJxdXV7rxRthkE3+8cKGfeyCZU+IXyQGzZ8PIkc03MLn8crj7bj8ZKldsuin89rdheerU5ktOSHoiT/xmVmhmr5rZI0F5CzN7wcwWmNk9ZtYz6hhEulo2TeC69Va/0FpqpEvv3n6W6xVX5OZGJscfDwcfHJZPO82v5Cnpi+O/xTnAm03K1wDXO+e2ApYDp8YQg0iXyZYJXPX1fjGz8vJwg/Kvf91fAR93XKyhxMrMT0TbaCNfXrAALrss2Zi6m0gTv5ltDhwE3BaUDRgFzAzuMh04PMoYRHLRp5/CoYfClClh3U47+U7ckSOTiysuJSXNN2W//nq/H7CkJ+or/inAhUCwMgibACucc8H1CR8AJa090MzKzWyOmc1ZtmxZxGGKdB/vvuu3R/zrX8O6I4/0wzU33zy5uOJ26qnwgx/448ZGGDMGVq1KNqbuIrLEb2YHA0udcy935PHOuSrn3Ajn3IiBAwd2cXQi3dPTT/sr+vnzw7rKSj8xa4MNkosrCWa+fyPVeT1/PkycmGxM3UWUV/x7AIea2fvA3fgmnhuAfmYWbKnM5oC2WBBJw7Rp8MMf+sXKwO9WdeedcNVVudmJm44ttoBJk8Ly1VfD3LnJxdNdRPbfxTl3iXNuc+dcGXAc8KRzbjQwG0jNHzwJeDCqGESiEPcEroYGv+H5qaf6pZXBD2ucPRtGj47klN3KuHGwxx7+uL7eN/mkOruldUlcJ1wEnG9mC/Bt/tpRU7qduCZwff65n7E6eXJYt/32fjnl3XeP5JTdTkGB37SlVy9ffuWV5q+XrCuWxO+c+6dz7uDg+F3n3Ejn3FbOuaOdc+qOEWnF++/7TtxHHgnrDj0Unn0WSksTCysrbbMNXHllWL7iCvjPfxILJ+vlacugSHZ79lnfifvGG2HdRRf5teg33DC5uLLZz38erke0apVv8kktUifNKfGLdECUM3fvuANGjYLUKOYePfxqlJMm5W8nbjqKinwHeFEwdOS555ovVich/TcSyVBUM3cbG+GSS+Ckk2D1al83YIDfOeukkzr11Hlj++1hwoSwPGGCn/cgzSnxi2SBL77wk7CaDk3cbjs/E3fPPZOLqzuqrPSvHfg1fE4/HbJg8dSsosQvkrCFC31yf7DJwOaDDvLt/GVliYXVbfXs6Zt8Us1iTz4Jt92WbEzZRolfJEHPP+87cZtOOvr5z/2HQGoRMsncLrv41zFl/Hj44IPk4sk2SvwiCZkxA/bZB5Ys8eUePfyV6eTJUFiYaGg54corYehQf/zZZ3DmmWrySVHiF8lQZ2fuNjb6Dc9Hjw4XFdtkE3j8cT87V7pGnz5+YlfKX/7iP2xFiV+kQzo6c/fLL/12iE0XE9t2W7+k8N57RxRswqqrqykrK6OgoICysjKqq6tjO/f3vueXdEg5++zwG1Y+U+IXickHH8Bee8GsWWHdAQf48ebf/GZycUWpurqa8vJyampqcM5RU1NDeXl5rMn/6qthyBB//Mkn8LOfxXbqrKXEL9IBmU7gevFF34n7yith3bnnwsMPw9e+FnGwCaqsrKSuxb6IdXV1VFZWxhZD375++eaU++5r/uGbj5T4RTKU6QSue+7xzTgffujLRUV+68Drrw9nmeaqhQsXZlQflf32g1NOCcsVFf7qP18p8YtExDm/WNhxx8FXX/m6/v3hscfgjDMSDS02Q1JtLGnWR+naa2HQIH+8ZAmcf37sIWQNJX6RCKxc6RN+0xUjt9nGd+KOGpVcXHGbOHEixaktsgLFxcVMTGCrrP794aabwvL06c23r8wnSvwiXWzRIt+Je++9Yd2++/rJWqlx5fli9OjRVFVVUVpaiplRWlpKVVUVoxPaQeaww/wHcsoZZ/gx/vnGotpAoiuNGDHCzZkzJ+kwRID1L8Y2Z47j0EN98k8ZNw6mTMn99vzuYtkyGDYMPvrIl888s/k3gVxiZi8750a0rNcVv0iX+THf+16Y9AsL/bLAv/udkn42GTgQbrwxLN98M/zzn4mFkwglfpEuUQnMZOVKX+rXz7cfN508JNnj2GP9bmYpp53mV/LMF0r8Ip3SG6gGrlpbM3Sob8/fd9/EgpJ2mPnmndQcinfegUsvTTamOCnxi3TYZsBs4IS1Nd//vk/622yTWFCSpsGD4brrwvL11/t/u3ygxC/SIYOA54HdmtTdwmOPwcYbJxSSZOyUU8JvZs75fXpTC+flMiV+kYz1BGYBZUG5ATgbOJMePZKKSTrCDKqqYIMNfPnNN+Gqq9b/mFygxC+SAT/6eSrhlX49cDhwY1sPkSxXVtZ8y8urr4bXXkssnFgo8YtkYOpUgKaL5o8HHkkmGOkyFRXh3sYNDb7JZ82aZGOKkhK/SJqeesqvqBm6A7ghoWikKxUU+E1bevf25Vdfhd/8JtmYoqTEL5KGmho46iior0/VvATkyUpreWLrreGXvwzLV17p2/xzUUaJ38yK27+XSG6pq4Mjjgin+MMS4Ejgq+SCkkicdx6MCBY4WL3aN/k0NCQbUxTSSvxm9l0zmw/8JyjvYGZTI41MJAs452d1vvqqL/tROz8GPkgwKolKURFMm8ba0VnPP998eYdcke4V//XA/sDHAM65ucBeUQUlki2uvRbuuiss+yTwbFLhSAy+/W1oukHYhAl+Zm8uSbupxzn3vxZVOfgFSCT02GNw0UVhubw8fzZQyXeXXOI/AMDvrXD66amhvNGJc1P6dBP//8zsu4Azsx5mNh7I0W4PEViwwK/b3tjoy3vskZtf+aV1PXv6Jp+CIEPOnt18396uFvem9Okm/jOBcUAJUAvsGJRFcs7nn8Phh8OKFb5cUgIzZ/pkAFBYWNjq49qql+5pxAgYPz4sjx8P/2vZ7tFF4t6UXhuxiDTR2OiHbd5/vy/36gXPPAO77BLeZ30bsXSH95Okb+VK2HFHeOstXz7wQHjkEb/UQ1cqKCho9f+OmdGY+trZAZ3aiMXMpptZvybl/mY2rcPRiGSpiRPDpA9+HZemSR/aTvzr+0CQ7qlPHz+xK/VP++ijcOedXX+euDelT7epZ3vn3IpUwTm3HNgpkohEEvLQQ3DZZWH5nHPgxBPXvd8GqRW90qyX7m3PPeGss8LyOefA4sVde464N6VPN/EXmFn/VMHMNga0mZzkjDffhJ/8JCyPGgWTJ7d+3y+++CKjeun+/u///GJuAMuXN/8g6Apxb0qfVhu/mZ0ITADuAww4CpjonPtTJFG1oDZ+idKKFTByJLz9ti+XlcFLL8GAAa3fX238+enxx2G//cLyzJnw4x8nF086OtXG75y7Az9HfQmwGDgyrqQvEqWGBjjhhDDpFxfDAw+0nfQlf+27L5zaZGHWcePg44+Ti6cz1pv4zWyj4OfG+IQ/I7gtDupEurVLL/Wboqf84Q+www7JxSPZbfJkv2UjwJIlfm2f7qi9K/4Zwc+XgTlNbqlym8yst5m9aGZzzWyemV0Z1G9hZi+Y2QIzu8fMenbybxDpkHvv9ZtupFx8MRxzTHLxSPbr189v0p7ypz/5kT5doaSkBDNbeyspKemaJ27FehO/c+5g8w2aezvntmxy28I5t2U7z70KGOWc2wE/4esAM9sNuAa43jm3FbCc5rtaiMRi7ly/32rKj36UH1vuSecdeigcf3xYPuMM+Oyzzj1nSUkJixYtala3aNGiyJJ/u238zvdW/SXTJ3ZeaphDj+DmgFHAzKB+On7fOpHYfPSRn5mbmig5dCjMmAGaeCvp+u1vYeBAf/zBB3DhhZ17vpZJv736zkp3OOcrZrZL+3drzswKzew1YCnwOPAOsMI5l9rO4gP8MhCtPbbczOaY2Zxly5ZlemqRVtXXw7HHwvvv+3LfvvDgg/4rvEi6BgxovnbTLbf49Xy6i3QT/67A82b2jpn928xeN7N/t/cg51yDc25HYHNgJPCtdANzzlU550Y450YMTH20inTSBRfAk0+G5TvvhG23TS4e6b6OOcZ/c0w57TT48svk4slEuol/f2BLfDPNIcDBwc+0BLN+ZwO7A/3MLDX5a3P8om8ikbvjDpgyJSxfeaVvrxXpCDOYOjX8tvjuu36UWEcMTg0VSrO+s9obztnbzM4FLgAOAGqdczWpWzuPHZha38fM+gD74pdyno2fAAZwEvBgJ/8GkXa99JJfTz/liCPgF7/o2HO1NUlLk7fyz6BBcN11YXnKFHjuucyfp7a2dp0kP3jwYGpro7kubu+KfzowAngd+BFwbQbPPQiYHTQJvQQ87px7BLgION/MFgCbALdnHLVIBhYv9ol+1SpfHj4cpk8P11rPVEVFRUb1kttOPjmc0eucn+T1VQe2Y66trcU5t/YWVdKHdpZsMLPXnXPfDo6LgBedc9+JLJo2aMkG6ajVq/26O88GuyX26+ev/rfaquPPWVRUREMrO3AXFhZSX1/fyiMk19XUwHbbQWq5pgkT/EqvSevokg1rUgdNRuKIdBtnnx0m/YICuPvuziV9oNWkv756yX2lpXDNNWH5mmvg1Vcze46mk7dSt6i0l/h3MLPPgtvnwPapYzPr5JQFkWjdcou/pUyaBPvv3/nn1Q5c0pozz4S99vLHDQ0wZgysWbP+x6TEvcdDezN3C51zGwW3vs65oibHG0USkUgX+Ne/4Gc/C8vHH998G73OKG/aS5xGveSHggK47Tbo3duXX3sNfv3rZGNqi7ZelJzzwQd+v9QlS3x5p538B0GLfS46Zfjw4cyfP39tediwYcybN6/rTiDd1uTJfr4I+H2aX30Vhg1b/2OiWuq7U8syi3QXK1f6ETyppD9ggN9KsSuTfkVFRbOkDzB//nyN6hHAr9g5cqQ/Xr3aN/lkW/ePEr/kDOd8O2vqy2Fhod8so7S0a89TVVWVUb3kl8JCmDYNevTw5Rde8Gv7ZBMlfskZN9zgZ+emTJkCe+/d9efRqB5pz/DhzScIVlbCggVt3z/uSYFK/JITnniieeftmDF+h6QoaFSPpOPii2H77f3xypVw+unQ2Nj2/ZtO3krdoqLEL93ee+/5BbNSF9y77urXUIlqGLRG9Ug6evb0TT6p64F//hOypTVQiV+6tS+/9CskfvKJLw8aBLNmQa9e0Z1z6tSpjB07du0VfmFhIWPHjmXq1KnRnVS6pZ13Dkf4gF+3f+HC5OJJ0XBO6bac82vr33efL/fsCU89BbvtlmxcIk199RXsuCP897++fMABfrvGCCfmrqXhnJJzJk0Kkz745p24kn5FRQVFRUWYGUVFRRrKKW3q3Rtuvz1M9H/7m9+rN0lK/NItPfqoHymRMm6cXxUxDhUVFdx0001rR/E0NDRw0003KflLm/bYo/lM8nPP9avGJkVNPdLtvPWWnyDz6ae+vNde8I9/hOOmo6bVOaUjvvwSvv1tPxgB4Mgj4c9/jvacauqRnPDZZ3DYYWHS/8Y3fHNPXEkfNI5fOmaDDeDWW8PyrFl+gmESlPil22hshJ/8BP7zH1/u3RseeAA23TTeODSOXzrqBz/we/OmjBsHH38cfxxK/NJtXHEFPPxwWL79dvhO7NsCwT777JNRvUhTkydDSYk/XrrUt/fHTYlfuoVZs+BXvwrL48fDCSckE8uCNubet1Uv0tTXvgY33xyW77wT/vKXeGNQ4pes98YbcOKJYXnffeHqq5OLZ2EbM3Daqhdp6eCDYfTosHzGGWG/VRyU+CWrffKJ78z98ktf3nJLv31iUVFyMQ0ZMiSjepHWTJkCAwf649ra5jN8o6bEL1mrvh6OOw7efdeXN9gAHnwQNt442bjWtLGfXlv1Iq0ZMAB+//uwfOut8OST8ZxbiV+y1iWXwOOPh+U77oDttksunpRFixZlVC/SlqOO8hsHpZx2WvjtNkpK/JKVZszwox9SLr3UT3gRySVm/qq/f39ffu+95jPSo6LEL1nnlVeaL79wyCF+KKdILho0CK6/Piz/9rfw//5ftOdU4pessnSpX2b5q698+Vvf8sPdCrLof+qwNnbObqtepD0nnuhX7QS/6uyYMeF7IApZ9HaSfLdmDRx9NPzvf7680Ua+M3ejjZKNq6V58+atk+SHDRvGvHnzEopIujszuOUW2HBDX37rLb+rXFSU+CVrnHcePP20PzaDu+6CrbdONqa2zJs3r9kWeUr60llDhsCvfw3bbuubeg46KLpzKfFLVrj99uZD2yZOhAMPTC4ekSSccQa8+mr0+0oo8Uvinn8emi5lf/TRfqNqkXxTUBDttqFrzxP9KUTatmiRH6a5erUvb789/OEP8WxLJ5JN4tzVLcGJ75LvVq3ySf/DD3154439MssbbJBsXCJxS+3qlpLa1Q1g6tSpXX4+7cAliXDOz1KcNs2XCwvhscf8euUi+SaqXd20A5dklalTw6QPfpaukr7kq7h3dVPil9g99VTzzSd++lM455zk4hFJWty7uinxS6xqavzCVKlvryNG+Ikr6syVfFZeXp5RfWepc1diU1fnVyL86CNf3nRTuP9+6NMn2bhEkpbqwK2qqqKhoYHCwkLKy8sj6diFPOncLSwspLGxcW25oKAgsrYzaZ1zfsehu+7y5R49/Nrje+6ZbFwiuSxvO3dbJn2AxsbGyNrOpHXXXhsmfYAbb1TSF0lKZInfzL5hZrPNbL6ZzTOzc4L6jc3scTN7O/jZP6oYgHWSfnv10vUeewwuuigsl5f7qekiEqqurqasrIyCggLKysqorq6O7FxRXvHXAz93zg0DdgPGmdkw4GLgCefcUOCJoCw5asECv31i6nN2jz381b6IhKqrqykvL6empgbnHDU1NZSXl0eW/CNL/M65D51zrwTHnwNvAiXAYcD04G7TgcOjikGS9fnnfm39FSt8uaQEZs6Enj2TjUsk21RWVlJXV9esrq6ujsqItuOKpY3fzMqAnYAXgM2cc8EkfRYDm7XxmHIzm2Nmc5YtWxZHmNKFGhvhpJMgtVpxr15+BM/Xv55sXCLZaOHChRnVd1bkid/MNgT+DJzrnPus6e+cH1LU6rAi51yVc26Ec27EwIEDO3z+uCdGiDdxok/0KVVVsMsuycUjks2GDBmSUX1nRZr4zawHPulXO+dmBdVLzGxQ8PtBwNIoY4h7YoTAQw/BZZeF5XPO8VvLiUjrJk6cmFF9Z0U2jt/MDN+G/4lz7twm9b8BPnbOTTKzi4GNnXMXru+5OjuOv6KiIraJEfnuzTdh1119+z7AqFF+VE+RpgqKtMnWM3W9Mzm6rXH8USb+PYFngNeB1NjJCfh2/nuBIUANcIxz7pP1PZdW5+weVqyAkSPh7bd9uawMXnoJBgxINCyRrBd34o/sOsw59y+grb9G6zDmmIYGOOGEMOkXF/u19ZX0RbJPzs/cBejZsydmtvbWU+MJu9yll8Jf/xqW//AH2GGH5OIRkbblfOLv2bMna9asaVa3Zs0aJf8udO+9cPXVYfnii+GYY5KLR0TWL+cTf8uk3169ZGbuXDjllLD8ox/BVVclF49Id9RWO35UfbAaayEd9tFHfmZuasLh0KEwY4bfRlFEMhPnSsk5f8Uv0aivh2OPhfff9+W+feHBB6Ffv0TDEpE0KPFLh1xwgV9PP+XOO2HbbZOLR0TSp8QvGbvjDpgyJSxfeSUcemhy8YjkgpKSkmajD0tKSiI7lxK/ZOTFF/16+ilHHAG/+EVy8YjkgpKSEhYtWtSsbtGiRZElfyV+SdvixXDkkbBqlS8PHw7Tp0OB/heJdErLpN9efWfl/Ft22LBhGdXnu/79+zf7utm/v98gbfVqOOooqK319+vXz8/M7ds3wWBFpENyfjjn/PnzM6rPZ/3792dFateUwIoVK+jfvz/HHrucZ5/1dQUFcPfdsNVWCQQpIp2W81f8kr6WST+sP4ZbbgnLkybB/vvHFJRIHhg8eLZD9UIAAAnYSURBVHBG9Z2V81f83nbADsCKZrcvvoANNoD1LIwn7AH8bm3p+ONh/PjkohHJRbW1tet08A4ePJjaVNtqF8uTxH8I8H/r1Pbt62eZ9uvX8Vtuf3CU4PfR6QHATjvBbbfl8t8rkpyoknxr8iTxtz2dtKEBPv7Y3zoilz44+vXr16S5pzdwP6ktkQcM8FspFhcnFZ2IdJU8SfxzgRn4D4Dw1qfPYFau7Nwz59IHx/Lly5t08N4M7LI2xpkzobS0a84jIsnKk8Q/I7g1V1fnWLUKPv3U7x7VkVuufXAsX76cKVPgvPPCuilTYO+9O/d3ikj2yJPE37ZevWDTTf2tI3Ltg+NrX4Onngp/P2YMjBvXuRhFJLvkfeLvrFz+4Nh1V5g6NXv6IESkayjxJyyqD450P0xSa+mv63+88MKu9O79YazrhItI9JT4u7nOfnCsXh1+SGy99S74ju9i4FnAfw0wMyV/kRyixJ/nevaEgQP9DeYkHY6IxEBLNoiI5BklfhGRPKPELyKSZ5T4Za22OnDVsSuSW9S5K80oyYvkPl3xi4jkmZxP/Np6UUSkuZxP/PPmzVsnyQ8bNox58+YlFJGISLJyPvEDTJgwgdLSUsyM0tJSJkyYkHRIIiKJyfnO3erqasrLy6kLFqWpqamhvLwcgNGjRycZmohIInL+ir+ysnJt0k+pq6ujsrIyoYhERJKV84l/4cKFGdWLiOS6nE/8Q4YMyaheRCTX5XzinzhxYkb1IiK5LucT/4UXXphRvYhIrsv5xL9o0aKM6kVEcl3OJ34REWkussRvZtPMbKmZvdGkbmMze9zM3g5+9o/q/CIi0roor/j/CBzQou5i4Ann3FDgiaAcqcGDB2dULyKS6yJL/M65p4FPWlQfBkwPjqcDh0d1/pTa2tp1kvzgwYOpra2N+tQiIlkp7iUbNnPOfRgcLwY2a+uOZlYOlEPnx9wryYuIhBLr3HV+x482d/1wzlU550Y450YMHDgwxshERHJb3Il/iZkNAgh+Lo35/CIieS/uxP8QcFJwfBLwYMznFxHJe1EO57wLeA7Yxsw+MLNTgUnAvmb2NvDDoCwiIjGKrHPXOXd8G7/6QVTnFBGR9mnmrohInlHiFxHJM+ZHVWY3M1sG1HTBUw0APuqC5+lK2RgTZGdciil92RiXYkpPV8ZU6pxbZzx8t0j8XcXM5jjnRiQdR1PZGBNkZ1yKKX3ZGJdiSk8cMampR0Qkzyjxi4jkmXxL/FVJB9CKbIwJsjMuxZS+bIxLMaUn8pjyqo1fRETy74pfRCTvKfGLiOSZnEz8ZnaAmf3XzBaY2Tq7fJlZLzO7J/j9C2ZWlgUxnWxmy8zsteB2WgwxrbM9Zovfm5n9Noj532b2nSyIaR8z+7TJ63RZDDF9w8xmm9l8M5tnZue0cp9YX6s0Y0riteptZi+a2dwgritbuU+s7780Y4r9/Rect9DMXjWzR1r5XXSvk3Mup25AIfAOsCXQE5gLDGtxnwrg5uD4OOCeLIjpZOB3Mb9WewHfAd5o4/cHAn8FDNgNeCELYtoHeCTm12kQ8J3guC/wViv/frG+VmnGlMRrZcCGwXEP4AVgtxb3ifv9l05Msb//gvOeD8xo7d8pytcpF6/4RwILnHPvOudWA3fjt3xsqukWkDOBH5iZJRxT7Fzr22M2dRhwh/OeB/ql9lNIMKbYOec+dM69Ehx/DrwJlLS4W6yvVZoxxS74+78Iij2CW8sRJLG+/9KMKXZmtjlwEHBbG3eJ7HXKxcRfAvyvSfkD1n1DrL2Pc64e+BTYJOGYAH4cNBPMNLNvRBhPutKNO267B1/b/2pmw+M8cfB1eyf8VWNTib1W64kJEnitguaL1/AbLT3unGvztYrp/ZdOTBD/+28KcCHQ2MbvI3udcjHxd1cPA2XOue2Bxwk/6aW5V/Drj+wA3Ag8ENeJzWxD4M/Auc65z+I67/q0E1Mir5VzrsE5tyOwOTDSzLaL47ydjCnW95+ZHQwsdc69HOV52pKLib8WaPppvXlQ1+p9zKwI+BrwcZIxOec+ds6tCoq3ATtHGE+60nktY+Wc+yz1td059yjQw8wGRH1eM+uBT7DVzrlZrdwl9teqvZiSeq2anH8FMBs4oMWv4n7/tRtTAu+/PYBDzex9fNPvKDO7s8V9InudcjHxvwQMNbMtzKwnvlPkoRb3aboF5FHAky7oQUkqphbtwYfi22yT9hBwYjBiZTfgU+fch0kGZGZfT7VzmtlI/P/hSJNGcL7bgTedc9e1cbdYX6t0YkrotRpoZv2C4z7AvsB/Wtwt1vdfOjHF/f5zzl3inNvcOVeGzwdPOud+0uJukb1Oke3AlRTnXL2ZnQU8hh9NM805N8/MfgnMcc49hH/D/MnMFuA7Eo/LgpjONrNDgfogppOjjAnWbo+5DzDAzD4ALsd3fOGcuxl4FD9aZQFQB5ySBTEdBYw1s3pgJXBcxB/a4K/Ofgq8HrQTA0wAhjSJK+7XKp2YknitBgHTzawQ/0Fzr3PukSTff2nGFPv7rzVxvU5askFEJM/kYlOPiIishxK/iEieUeIXEckzSvwiInlGiV9EJM8o8UteMbOGYPXFN8zsPjMr7sBzXGFm46OITyQOSvySb1Y653Z0zm0HrAbOTDogkbgp8Us+ewbYCsDMHjCzl4P12stTdzC/j8IrwUJnT7R8AjM7PVgArU9w/FJw3z+nvk2Y2TfN7Hkze93MrjKzL5o8/oLgMf+2VtaJF4mCEr/kpWDtkx8BrwdVY5xzOwMj8LM4NzGzgcCtwI+Dhc6ObvEcZwEHA4c751YCs5xzuwT3fRM4NbjrDcANzrlv41ftTD1+P2AoftnuHYGdzWyvaP5ikVDOLdkg0o4+TZY4eAY/LR58sj8iOP4GPiEPBJ52zr0H4Jxruk/Aifglcw93zq0J6rYzs6uAfsCG+CU6AHYHDg+OZwCTg+P9gturQXnD4LxPd/aPFFkfJX7JNyuD5XnXMrN9gB8Cuzvn6szsn0Dvdp7ndfxV+ubAe0HdH/EfBHPN7GT8mkPrY8DVzrlbMohfpNPU1CPil7tdHiT9b+G3TgR4HtjLzLYAMLONmzzmVeAM4CEzGxzU9QU+DJZLHt3kvs8DPw6Omy609RgwJlhTHzMrMbNNu/DvEmmVEr8I/A0oMrM3gUn4RI1zbhlQDswys7nAPU0f5Jz7FzAe+Euwzv2l+F2wnqX5sr/nAueb2b/xncmfBo//O77p5zkzex2/vV7fqP5IkRStzikSsWB0z0rnnDOz44DjnXOJ77ks+Utt/CLR2xn4XbApygpgTMLxSJ7TFb+ISJ5RG7+ISJ5R4hcRyTNK/CIieUaJX0Qkzyjxi4jkmf8PbHQbzxSKCCoAAAAASUVORK5CYII=\n" }, "metadata": { "needs_background": "light" } } ], "source": [ "from sklearn.preprocessing import PolynomialFeatures\n", "from sklearn.pipeline import make_pipeline\n", "\n", "pipeline = make_pipeline(PolynomialFeatures(4), LinearRegression())\n", "\n", "X_train, X_test, y_train, y_test = train_test_split(X, y, test_size=0.2, random_state=0)\n", "\n", "pipeline.fit(np.array(X_train), y_train)\n", "\n", "y_pred=pipeline.predict(X_test)\n", "\n", "df = pd.DataFrame({'x': X_test[:,0], 'y': y_pred[:,0]})\n", "df.sort_values(by='x',inplace = True)\n", "points = pd.DataFrame(df).to_numpy()\n", "\n", "plt.plot(points[:, 0], points[:, 1],color=\"blue\", linewidth=3)\n", "plt.xlabel('Package')\n", "plt.ylabel('Price')\n", "\n", "plt.scatter(X,y, color=\"black\")\n", "\n", "plt.show()\n", "\n" ] }, { "cell_type": "code", "execution_count": 17, "metadata": {}, "outputs": [ { "output_type": "stream", "name": "stdout", "text": [ "Model Accuracy: 0.8537946517073784\n" ] } ], "source": [ "accuracy_score = pipeline.score(X_train,y_train)\n", "print('Model Accuracy: ', accuracy_score)\n", "\n" ] }, { "cell_type": "code", "execution_count": 18, "metadata": {}, "outputs": [ { "output_type": "execute_result", "data": { "text/plain": [ "array([[46.34509342]])" ] }, "metadata": {}, "execution_count": 18 } ], "source": [ "pipeline.predict( np.array([ [2.75] ]) )" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ] }