Mini-Batch Gradient Descent with Python
In machine learning, gradient descent is an optimization technique used for computing the model parameters (coefficients and bias) for algorithms like linear regression, logistic regression, neural networks, etc. In this technique, we repeatedly iterate through the training set and update the model parameters in accordance with the gradient of error with respect to the training set.
Depending on the number of training examples considered in updating the model parameters, we have 3-types of gradient descents:
- Batch Gradient Descent: Parameters are updated after computing the gradient of error with respect to the entire training set
- Stochastic Gradient Descent: Parameters are updated after computing the gradient of error with respect to a single training example
- Mini-Batch Gradient Descent: Parameters are updated after computing the gradient of error with respect to a subset of the training set
Batch Gradient Descent | Stochastic Gradient Descent | Mini-Batch Gradient Descent |
Since entire training data is considered before taking a step in the direction of gradient, therefore it takes a lot of time for making a single update. | Since only a single training example is considered before taking a step in the direction of gradient, we are forced to loop over the training set and thus cannot exploit the speed associated with vectorizing the code. | Since a subset of training examples is considered, it can make quick updates in the model parameters and can also exploit the speed associated with vectorizing the code. |
It makes smooth updates in the model parameters | It makes very noisy updates in the parameters | Depending upon the batch size, the updates can be made less noisy – greater the batch size less noisy is the update |
Thus, mini-batch gradient descent makes a compromise between the speedy convergence and the noise associated with gradient update which makes it a more flexible and robust algorithm.
Mini-Batch Gradient Descent:
Algorithm-
Let theta = model parameters and max_iters = number of epochs.
for itr = 1, 2, 3, …, max_iters:
for mini_batch (X_mini, y_mini):
- Forward Pass on the batch X_mini:
- Make predictions on the mini-batch
- Compute error in predictions (J(theta)) with the current values of the parameters
- Backward Pass:
- Compute gradient(theta) = partial derivative of J(theta) w.r.t. theta
- Update parameters:
- theta = theta – learning_rate*gradient(theta)
Below is the Python Implementation:
Step #1: First step is to import dependencies, generate data for linear regression and visualize the generated data. We have generated 8000 data examples, each having 2 attributes/features. These data examples are further divided into training set (X_train, y_train) and testing set (X_test, y_test) having 7200 and 800 examples respectively.
# importing dependencies import numpy as np import matplotlib.pyplot as plt # creating data mean = np.array([ 5.0 , 6.0 ]) cov = np.array([[ 1.0 , 0.95 ], [ 0.95 , 1.2 ]]) data = np.random.multivariate_normal(mean, cov, 8000 ) # visualising data plt.scatter(data[: 500 , 0 ], data[: 500 , 1 ], marker = '.' ) plt.show() # train-test-split data = np.hstack((np.ones((data.shape[ 0 ], 1 )), data)) split_factor = 0.90 split = int (split_factor * data.shape[ 0 ]) X_train = data[:split, : - 1 ] y_train = data[:split, - 1 ].reshape(( - 1 , 1 )) X_test = data[split:, : - 1 ] y_test = data[split:, - 1 ].reshape(( - 1 , 1 )) print ( "Number of examples in training set = % d" % (X_train.shape[ 0 ])) print ( "Number of examples in testing set = % d" % (X_test.shape[ 0 ])) |
Output:
Number of examples in training set = 7200
Number of examples in testing set = 800
Step #2: Next, we write the code for implementing linear regression using mini-batch gradient descent.
gradientDescent()
is the main driver function and other functions are helper functions used for making predictions – hypothesis()
, computing gradients – gradient()
, computing error – cost()
and creating mini-batches – create_mini_batches()
. The driver function initializes the parameters, computes the best set of parameters for the model and returns these parameters along with a list containing history of errors as the parameters get updated.
# linear regression using "mini-batch" gradient descent # function to compute hypothesis / predictions def hypothesis(X, theta): return np.dot(X, theta) # function to compute gradient of error function w.r.t. theta def gradient(X, y, theta): h = hypothesis(X, theta) grad = np.dot(X.transpose(), (h - y)) return grad # function to compute the error for current values of theta def cost(X, y, theta): h = hypothesis(X, theta) J = np.dot((h - y).transpose(), (h - y)) J / = 2 return J[ 0 ] # function to create a list containing mini-batches def create_mini_batches(X, y, batch_size): mini_batches = [] data = np.hstack((X, y)) np.random.shuffle(data) n_minibatches = data.shape[ 0 ] / / batch_size i = 0 for i in range (n_minibatches + 1 ): mini_batch = data[i * batch_size:(i + 1 ) * batch_size, :] X_mini = mini_batch[:, : - 1 ] Y_mini = mini_batch[:, - 1 ].reshape(( - 1 , 1 )) mini_batches.append((X_mini, Y_mini)) if data.shape[ 0 ] % batch_size ! = 0 : mini_batch = data[i * batch_size:data.shape[ 0 ]] X_mini = mini_batch[:, : - 1 ] Y_mini = mini_batch[:, - 1 ].reshape(( - 1 , 1 )) mini_batches.append((X_mini, Y_mini)) return mini_batches # function to perform mini-batch gradient descent def gradientDescent(X, y, learning_rate = 0.001 , batch_size = 32 ): theta = np.zeros((X.shape[ 1 ], 1 )) error_list = [] max_iters = 3 for itr in range (max_iters): mini_batches = create_mini_batches(X, y, batch_size) for mini_batch in mini_batches: X_mini, y_mini = mini_batch theta = theta - learning_rate * gradient(X_mini, y_mini, theta) error_list.append(cost(X_mini, y_mini, theta)) return theta, error_list |
Calling the gradientDescent()
function to compute the model parameters (theta) and visualizing the change in error function.
theta, error_list = gradientDescent(X_train, y_train) print ( "Bias = " , theta[ 0 ]) print ( "Coefficients = " , theta[ 1 :]) # visualising gradient descent plt.plot(error_list) plt.xlabel( "Number of iterations" ) plt.ylabel( "Cost" ) plt.show() |
Output:
Bias = [0.81830471]
Coefficients = [[1.04586595]]
Step #3 : Finally, we make predictions on the testing set and compute the mean absolute error in predictions.
# predicting output for X_test y_pred = hypothesis(X_test, theta) plt.scatter(X_test[:, 1 ], y_test[:, ], marker = '.' ) plt.plot(X_test[:, 1 ], y_pred, color = 'orange' ) plt.show() # calculating error in predictions error = np. sum (np. abs (y_test - y_pred) / y_test.shape[ 0 ]) print ( "Mean absolute error = " , error) |
Output:
Mean absolute error = 0.4366644295854125
The orange line represents the final hypothesis function: theta[0] + theta[1]*X_test[:, 1] + theta[2]*X_test[:, 2] = 0
Python Example for Beginners
Two Machine Learning Fields
There are two sides to machine learning:
- Practical Machine Learning:This is about querying databases, cleaning data, writing scripts to transform data and gluing algorithm and libraries together and writing custom code to squeeze reliable answers from data to satisfy difficult and ill defined questions. It’s the mess of reality.
- Theoretical Machine Learning: This is about math and abstraction and idealized scenarios and limits and beauty and informing what is possible. It is a whole lot neater and cleaner and removed from the mess of reality.
Data Science Resources: Data Science Recipes and Applied Machine Learning Recipes
Introduction to Applied Machine Learning & Data Science for Beginners, Business Analysts, Students, Researchers and Freelancers with Python & R Codes @ Western Australian Center for Applied Machine Learning & Data Science (WACAMLDS) !!!
Latest end-to-end Learn by Coding Recipes in Project-Based Learning:
Applied Statistics with R for Beginners and Business Professionals
Data Science and Machine Learning Projects in Python: Tabular Data Analytics
Data Science and Machine Learning Projects in R: Tabular Data Analytics
Python Machine Learning & Data Science Recipes: Learn by Coding
R Machine Learning & Data Science Recipes: Learn by Coding
Comparing Different Machine Learning Algorithms in Python for Classification (FREE)
Disclaimer: The information and code presented within this recipe/tutorial is only for educational and coaching purposes for beginners and developers. Anyone can practice and apply the recipe/tutorial presented here, but the reader is taking full responsibility for his/her actions. The author (content curator) of this recipe (code / program) has made every effort to ensure the accuracy of the information was correct at time of publication. The author (content curator) does not assume and hereby disclaims any liability to any party for any loss, damage, or disruption caused by errors or omissions, whether such errors or omissions result from accident, negligence, or any other cause. The information presented here could also be found in public knowledge domains.