Neural network prediction using Keras on a Data Frame.ipynb

{"cells":[{"cell_type":"markdown","id":"701048de-7a54-48c4-b441-9d3ddc412e2f","metadata":{},"source":["<a href=\"https://cognitiveclass.ai/?utm_medium=Exinfluencer&utm_source=Exinfluencer&utm_content=000026UJ&utm_term=10006555&utm_id=NA-SkillsNetwork-Channel-SkillsNetworkCoursesIBMDeveloperSkillsNetworkDL0101ENSkillsNetwork945-2022-01-01\"><img src=\"https://cf-courses-data.s3.us.cloud-object-storage.appdomain.cloud/IBMDeveloperSkillsNetwork-DL0101EN-SkillsNetwork/images/IDSN-logo.png\" width=\"400\"> </a>\n","\n","<h1 align=center><font size = 5>Regression Models with Keras</font></h1>\n"]},{"cell_type":"markdown","id":"bbaea2bd-fdfa-4cf2-a527-4a247a75ac0c","metadata":{},"source":["## Introduction\n"]},{"cell_type":"markdown","id":"aae0d793-e6d6-4e53-80fe-181ef4308c2a","metadata":{},"source":["As we discussed in the videos, despite the popularity of more powerful libraries such as PyToch and TensorFlow, they are not easy to use and have a steep learning curve. So, for people who are just starting to learn deep learning, there is no better library to use other than the Keras library. \n","\n","Keras is a high-level API for building deep learning models. It has gained favor for its ease of use and syntactic simplicity facilitating fast development. As you will see in this lab and the other labs in this course, building a very complex deep learning network can be achieved with Keras with only few lines of code. You will appreciate Keras even more, once you learn how to build deep models using PyTorch and TensorFlow in the other courses.\n","\n","So, in this lab, you will learn how to use the Keras library to build a regression model.\n"]},{"cell_type":"markdown","id":"022a0ec2-3cea-417e-8889-5c65caf971f0","metadata":{},"source":["<h2>Regression Models with Keras</h2>\n","\n","<h3>Objective for this Notebook<h3>    \n","<h5> 1. How to use the Keras library to build a regression model.</h5>\n","<h5> 2. Download and Clean dataset </h5>\n","<h5> 3. Build a Neural Network </h5>\n","<h5> 4. Train and Test the Network. </h5>     \n","\n"]},{"cell_type":"markdown","id":"ee476fb0-f7f9-4ee3-a8ca-5c9bb62056cd","metadata":{},"source":["## Table of Contents\n","\n","<div class=\"alert alert-block alert-info\" style=\"margin-top: 20px\">\n","\n","<font size = 3>\n","    \n","1. <a href=\"#item31\">Download and Clean Dataset</a>  \n","2. <a href=\"#item32\">Import Keras</a>  \n","3. <a href=\"#item33\">Build a Neural Network</a>  \n","4. <a href=\"#item34\">Train and Test the Network</a>  \n","\n","</font>\n","</div>\n"]},{"cell_type":"markdown","id":"7efcf931-5ad5-4324-ab17-31aa0e2b977a","metadata":{},"source":["<a id=\"item31\"></a>\n"]},{"cell_type":"markdown","id":"b1ea8f0d-afcb-48cf-90c1-cdd14bbebde5","metadata":{},"source":["## Download and Clean Dataset\n"]},{"cell_type":"markdown","id":"ec856cd2-190b-4644-92f7-77fce8ebb882","metadata":{},"source":["Let's start by importing the <em>pandas</em> and the Numpy libraries.\n"]},{"cell_type":"code","execution_count":null,"id":"37a30258-fcff-40db-96cb-789ca9060680","metadata":{"vscode":{"languageId":"python"}},"outputs":[],"source":["# All Libraries required for this lab are listed below. The libraries pre-installed on Skills Network Labs are commented. \n","# If you run this notebook on a different environment, e.g. your desktop, you may need to uncomment and install certain libraries.\n","\n","#!pip install numpy==1.21.4\n","#!pip install pandas==1.3.4\n","#!pip install keras==2.1.6"]},{"cell_type":"code","execution_count":null,"id":"d8c79c74-f278-4db8-94f5-c9101c2e9881","metadata":{"vscode":{"languageId":"python"}},"outputs":[],"source":["import pandas as pd\n","import numpy as np"]},{"cell_type":"markdown","id":"b4e0dcc2-99b1-466c-8086-6712cc5b40e4","metadata":{},"source":["We will be playing around with the same dataset that we used in the videos.\n","\n","<strong>The dataset is about the compressive strength of different samples of concrete based on the volumes of the different ingredients that were used to make them. Ingredients include:</strong>\n","\n","<strong>1. Cement</strong>\n","\n","<strong>2. Blast Furnace Slag</strong>\n","\n","<strong>3. Fly Ash</strong>\n","\n","<strong>4. Water</strong>\n","\n","<strong>5. Superplasticizer</strong>\n","\n","<strong>6. Coarse Aggregate</strong>\n","\n","<strong>7. Fine Aggregate</strong>\n"]},{"cell_type":"markdown","id":"6be3d86c-1d4d-43c2-88d0-0f4f681a331a","metadata":{},"source":["Let's download the data and read it into a <em>pandas</em> dataframe.\n"]},{"cell_type":"code","execution_count":null,"id":"d2a0e8cf-4872-406f-ae9f-7bbd13de4868","metadata":{"vscode":{"languageId":"python"}},"outputs":[],"source":["concrete_data = pd.read_csv('https://s3-api.us-geo.objectstorage.softlayer.net/cf-courses-data/CognitiveClass/DL0101EN/labs/data/concrete_data.csv')\n","concrete_data.head()"]},{"cell_type":"markdown","id":"2469361b-b36c-4ea0-9687-1d913f38ca0c","metadata":{},"source":["So the first concrete sample has 540 cubic meter of cement, 0 cubic meter of blast furnace slag, 0 cubic meter of fly ash, 162 cubic meter of water, 2.5 cubic meter of superplaticizer, 1040 cubic meter of coarse aggregate, 676 cubic meter of fine aggregate. Such a concrete mix which is 28 days old, has a compressive strength of 79.99 MPa. \n"]},{"cell_type":"markdown","id":"d709c0e6-fd74-4682-bd73-7512fc9d654f","metadata":{},"source":["#### Let's check how many data points we have.\n"]},{"cell_type":"code","execution_count":null,"id":"bb271395-98f0-42b8-8376-870dc5f4e0b6","metadata":{"vscode":{"languageId":"python"}},"outputs":[],"source":["concrete_data.shape"]},{"cell_type":"markdown","id":"34eae67d-3db4-4505-91fb-a7586a5e2688","metadata":{},"source":["So, there are approximately 1000 samples to train our model on. Because of the few samples, we have to be careful not to overfit the training data.\n"]},{"cell_type":"markdown","id":"fff37414-8422-49bd-a001-2106c5b88ac1","metadata":{},"source":["Let's check the dataset for any missing values.\n"]},{"cell_type":"code","execution_count":null,"id":"d2de8710-0378-4002-9a6f-c458b594e2d4","metadata":{"vscode":{"languageId":"python"}},"outputs":[],"source":["concrete_data.describe()"]},{"cell_type":"code","execution_count":null,"id":"691ea21d-a583-4183-9d7e-442da86aa485","metadata":{"vscode":{"languageId":"python"}},"outputs":[],"source":["concrete_data.isnull().sum()"]},{"cell_type":"markdown","id":"d31a4a08-b54a-4f2e-928b-676105401628","metadata":{},"source":["The data looks very clean and is ready to be used to build our model.\n"]},{"cell_type":"markdown","id":"bbe6524c-55c3-4fb9-a3e8-ed65f828c160","metadata":{},"source":["#### Split data into predictors and target\n"]},{"cell_type":"markdown","id":"fa62e25b-90cb-4a5c-b594-c172bacccbc8","metadata":{},"source":["The target variable in this problem is the concrete sample strength. Therefore, our predictors will be all the other columns.\n"]},{"cell_type":"code","execution_count":null,"id":"36b3e2ff-fe91-4749-a277-aa418ced9999","metadata":{"vscode":{"languageId":"python"}},"outputs":[],"source":["concrete_data_columns = concrete_data.columns\n","\n","predictors = concrete_data[concrete_data_columns[concrete_data_columns != 'Strength']] # all columns except Strength\n","target = concrete_data['Strength'] # Strength column"]},{"cell_type":"markdown","id":"6de05c92-ce6b-45b0-89d9-8247e82cef24","metadata":{},"source":["<a id=\"item2\"></a>\n"]},{"cell_type":"markdown","id":"c7fc5dcc-2f7e-49ef-8f4f-029bd6949f42","metadata":{},"source":["Let's do a quick sanity check of the predictors and the target dataframes.\n"]},{"cell_type":"code","execution_count":null,"id":"72fd40b0-b896-458e-8b01-5265e84b2f4f","metadata":{"vscode":{"languageId":"python"}},"outputs":[],"source":["predictors.head()"]},{"cell_type":"code","execution_count":null,"id":"aa50d52e-1237-4f54-9e74-3a6e53d6f492","metadata":{"vscode":{"languageId":"python"}},"outputs":[],"source":["target.head()"]},{"cell_type":"markdown","id":"54adb38c-defc-4156-98a2-0125354acc5f","metadata":{},"source":["Finally, the last step is to normalize the data by substracting the mean and dividing by the standard deviation.\n"]},{"cell_type":"code","execution_count":null,"id":"107277b0-9d5e-4ca1-8dcd-6ecea06ccfa7","metadata":{"vscode":{"languageId":"python"}},"outputs":[],"source":["predictors_norm = (predictors - predictors.mean()) / predictors.std()\n","predictors_norm.head()"]},{"cell_type":"markdown","id":"fdb24bf1-610f-4d16-99da-9aec9798b3e1","metadata":{},"source":["Let's save the number of predictors to *n_cols* since we will need this number when building our network.\n"]},{"cell_type":"code","execution_count":null,"id":"72c21bc6-ee8c-4d1e-9a8c-700d7f7a648a","metadata":{"vscode":{"languageId":"python"}},"outputs":[],"source":["n_cols = predictors_norm.shape[1] # number of predictors"]},{"cell_type":"markdown","id":"deaf62bf-fab4-47b4-b282-2f7e4ef0638a","metadata":{},"source":["<a id=\"item1\"></a>\n"]},{"cell_type":"markdown","id":"c7ad72d7-d37a-4cff-945b-e4424ede39f2","metadata":{},"source":["<a id='item32'></a>\n"]},{"cell_type":"markdown","id":"586f1914-04bc-4c2c-af43-63178902897e","metadata":{},"source":["## Import Keras\n"]},{"cell_type":"markdown","id":"0a777da9-1254-460c-a040-62cd694dba49","metadata":{},"source":["Recall from the videos that Keras normally runs on top of a low-level library such as TensorFlow. This means that to be able to use the Keras library, you will have to install TensorFlow first and when you import the Keras library, it will be explicitly displayed what backend was used to install the Keras library. In CC Labs, we used TensorFlow as the backend to install Keras, so it should clearly print that when we import Keras.\n"]},{"cell_type":"markdown","id":"fd147cc8-2e36-4bd6-a3a2-c4c9b9ed0c98","metadata":{},"source":["#### Let's go ahead and import the Keras library\n"]},{"cell_type":"code","execution_count":null,"id":"8dc60b83-f197-4a3e-a3bf-97843f6236c8","metadata":{"vscode":{"languageId":"python"}},"outputs":[],"source":["import keras"]},{"cell_type":"markdown","id":"dd5fb0cb-497f-48f8-b150-f272bc2fac4a","metadata":{},"source":["As you can see, the TensorFlow backend was used to install the Keras library.\n"]},{"cell_type":"markdown","id":"33ae3534-8421-48ef-97ae-5bd58dbce395","metadata":{},"source":["Let's import the rest of the packages from the Keras library that we will need to build our regressoin model.\n"]},{"cell_type":"code","execution_count":null,"id":"d57702ef-54f5-4485-a7b7-104f7c975f57","metadata":{"vscode":{"languageId":"python"}},"outputs":[],"source":["from keras.models import Sequential\n","from keras.layers import Dense"]},{"cell_type":"markdown","id":"f3169c0c-1f49-4ee4-ac0e-36d038a6690c","metadata":{},"source":["<a id='item33'></a>\n"]},{"cell_type":"markdown","id":"a49721c9-1ae0-43e1-80bb-7d87e21cb97a","metadata":{},"source":["## Build a Neural Network\n"]},{"cell_type":"markdown","id":"7f181e29-ad60-4a2a-8581-6f5c9e935d37","metadata":{},"source":["Let's define a function that defines our regression model for us so that we can conveniently call it to create our model.\n"]},{"cell_type":"code","execution_count":null,"id":"5e70c37c-ad25-4b70-b056-f021ccba6177","metadata":{"vscode":{"languageId":"python"}},"outputs":[],"source":["# define regression model\n","def regression_model():\n","    # create model\n","    model = Sequential()\n","    model.add(Dense(50, activation='relu', input_shape=(n_cols,)))\n","    model.add(Dense(50, activation='relu'))\n","    model.add(Dense(1))\n","    \n","    # compile model\n","    model.compile(optimizer='adam', loss='mean_squared_error')\n","    return model"]},{"cell_type":"markdown","id":"78496132-59ef-495f-b062-cbf6ef3cbc72","metadata":{},"source":["The above function create a model that has two hidden layers, each of 50 hidden units.\n"]},{"cell_type":"markdown","id":"b737d68b-a96b-48fe-9c5f-94d0d7e6b97b","metadata":{},"source":["<a id=\"item4\"></a>\n"]},{"cell_type":"markdown","id":"192d414a-951d-472e-ad0e-cf004a2f51ce","metadata":{},"source":["<a id='item34'></a>\n"]},{"cell_type":"markdown","id":"5961862e-19ff-4fd4-81a9-dfcce6cc9f8e","metadata":{},"source":["## Train and Test the Network\n"]},{"cell_type":"markdown","id":"9fd7a324-6b75-4689-aaa2-4a30ba85721d","metadata":{},"source":["Let's call the function now to create our model.\n"]},{"cell_type":"code","execution_count":null,"id":"d02d1ac4-d56b-4879-ab94-9c156f690bbf","metadata":{"vscode":{"languageId":"python"}},"outputs":[],"source":["# build the model\n","model = regression_model()"]},{"cell_type":"markdown","id":"5dcdc296-a32e-440f-afab-baa8d033f074","metadata":{},"source":["Next, we will train and test the model at the same time using the *fit* method. We will leave out 30% of the data for validation and we will train the model for 100 epochs.\n"]},{"cell_type":"code","execution_count":null,"id":"7ebe5f26-9a52-4506-a777-b651608e3840","metadata":{"vscode":{"languageId":"python"}},"outputs":[],"source":["# fit the model\n","model.fit(predictors_norm, target, validation_split=0.3, epochs=100, verbose=2)"]},{"cell_type":"markdown","id":"3d5114dd-f11c-4ca3-9b3a-ea440f164626","metadata":{},"source":["<strong>You can refer to this [link](https://keras.io/models/sequential/?utm_medium=Exinfluencer&utm_source=Exinfluencer&utm_content=000026UJ&utm_term=10006555&utm_id=NA-SkillsNetwork-Channel-SkillsNetworkCoursesIBMDeveloperSkillsNetworkDL0101ENSkillsNetwork945-2022-01-01) to learn about other functions that you can use for prediction or evaluation.</strong>\n"]},{"cell_type":"markdown","id":"dbeb9675-b233-4361-9c87-256854235211","metadata":{},"source":["Feel free to vary the following and note what impact each change has on the model's performance:\n","\n","1. Increase or decreate number of neurons in hidden layers\n","2. Add more hidden layers\n","3. Increase number of epochs\n"]},{"cell_type":"markdown","id":"63885dc9-f204-4ecc-a4bd-8aecd31246a8","metadata":{},"source":["### Thank you for completing this lab!\n","\n","This notebook was created by [Alex Aklson](https://www.linkedin.com/in/aklson/?utm_medium=Exinfluencer&utm_source=Exinfluencer&utm_content=000026UJ&utm_term=10006555&utm_id=NA-SkillsNetwork-Channel-SkillsNetworkCoursesIBMDeveloperSkillsNetworkDL0101ENSkillsNetwork945-2022-01-01). I hope you found this lab interesting and educational. Feel free to contact me if you have any questions!\n"]},{"cell_type":"markdown","id":"4cab296b-207d-4105-b544-37a12e0642fd","metadata":{},"source":["\n","## Change Log\n","\n","|  Date (YYYY-MM-DD) |  Version | Changed By  |  Change Description |\n","|---|---|---|---|\n","| 2020-09-21  | 2.0  | Srishti  |  Migrated Lab to Markdown and added to course repo in GitLab |\n","\n","\n","\n","<hr>\n","\n","## <h3 align=\"center\"> © IBM Corporation 2020. All rights reserved. <h3/>\n"]},{"cell_type":"markdown","id":"67bd639d-9eaa-4f5f-ad93-8ae631980288","metadata":{},"source":["This notebook is part of a course on **Coursera** called *Introduction to Deep Learning & Neural Networks with Keras*. If you accessed this notebook outside the course, you can take this course online by clicking [here](https://cocl.us/DL0101EN_Coursera_Week3_LAB1).\n"]},{"cell_type":"markdown","id":"bbc69e6a-816d-4314-836e-dc097d67d86a","metadata":{},"source":["<hr>\n","\n","Copyright &copy; 2019 [IBM Developer Skills Network](https://cognitiveclass.ai/?utm_medium=dswb&utm_source=bducopyrightlink&utm_content=000026UJ&utm_term=10006555&utm_id=NA-SkillsNetwork-Channel-SkillsNetworkCoursesIBMDeveloperSkillsNetworkDL0101ENSkillsNetwork945-2022-01-01&utm_campaign=bdu). This notebook and its source code are released under the terms of the [MIT License](https://bigdatauniversity.com/mit-license/?utm_medium=Exinfluencer&utm_source=Exinfluencer&utm_content=000026UJ&utm_term=10006555&utm_id=NA-SkillsNetwork-Channel-SkillsNetworkCoursesIBMDeveloperSkillsNetworkDL0101ENSkillsNetwork945-2022-01-01).\n"]}],"metadata":{"kernelspec":{"display_name":"Python","language":"python","name":"conda-env-python-py"},"language_info":{"name":""}},"nbformat":4,"nbformat_minor":4}