ML Cheatsheet

ML Cheatsheet Documentation
Team
Apr 18, 2021

Basics
1 Linear Regression 3
2 Gradient Descent 21
3 Logistic Regression 25
4 Glossary 39
5 Calculus 45
6 Linear Algebra 57
7 Probability 67
8 Statistics 69
9 Notation 71
10 Concepts 75
11 Forwardpropagation 81
12 Backpropagation 91
13 Activation Functions 97
14 Layers 105
15 Loss Functions 111
16 Optimizers 115
17 Regularization 121
18 Architectures 131
19 Classification Algorithms 145
20 Clustering Algorithms 151
i
21 Regression Algorithms 153
22 Reinforcement Learning 155
23 Datasets 157
24 Libraries 173
25 Papers 203
26 Other Content 209
27 Contribute 215
ii
Brief visual explanations of machine learning concepts with diagrams, code examples and links to resources for
learning more.
Warning: This document is under early stage development. If you find errors, please raise an issue or contribute
a better definition!
Basics 1
2 Basics
CHAPTER 1
Linear Regression
• Introduction
• Simple regression
– Making predictions
– Cost function
– Gradient descent
– Training
– Model evaluation
– Summary
• Multivariable regression
– Growing complexity
– Normalization
– Initialize weights
– Cost function
– Simplifying with matrices
– Bias term
3
1.1 Introduction
Linear Regression is a supervised machine learning algorithm where the predicted output is continuous and has a
constant slope. It’s used to predict values within a continuous range, (e.g. sales, price) rather than trying to classify
them into categories (e.g. cat, dog). There are two main types:
Simple regression
Simple linear regression uses traditional slope-intercept form, where 𝑚 and 𝑏 are the variables our algorithm will try
to “learn” to produce the most accurate predictions. 𝑥 represents our input data and 𝑦 represents our prediction.
𝑦 = 𝑚𝑥 + 𝑏
Multivariable regression
A more complex, multi-variable linear equation might look like this, where 𝑤 represents the coefficients, or weights,
our model will try to learn.
𝑓 (𝑥, 𝑦, 𝑧) = 𝑤1 𝑥 + 𝑤2 𝑦 + 𝑤3 𝑧
The variables 𝑥, 𝑦, 𝑧 represent the attributes, or distinct pieces of information, we have about each observation. For
sales predictions, these attributes might include a company’s advertising spend on radio, TV, and newspapers.
𝑆𝑎𝑙𝑒𝑠 = 𝑤1 𝑅𝑎𝑑𝑖𝑜 + 𝑤2 𝑇 𝑉 + 𝑤3 𝑁 𝑒𝑤𝑠
1.2 Simple regression
Let’s say we are given a dataset with the following columns (features): how much a company spends on Radio
advertising each year and its annual Sales in terms of units sold. We are trying to develop an equation that will let us
to predict units sold based on how much a company spends on radio advertising. The rows (observations) represent
companies.
Company Radio ($) Sales

Amazon 37.8 22.1
Google 39.3 10.4
Facebook 45.9 18.3
Apple 41.3 18.5
1.2.1 Making predictions
Our prediction function outputs an estimate of sales given a company’s radio advertising spend and our current values
for Weight and Bias.
𝑆𝑎𝑙𝑒𝑠 = 𝑊 𝑒𝑖𝑔ℎ𝑡 · 𝑅𝑎𝑑𝑖𝑜 + 𝐵𝑖𝑎𝑠
Weight the coefficient for the Radio independent variable. In machine learning we call coefficients weights.
Radio the independent variable. In machine learning we call these variables features.
Bias the intercept where our line intercepts the y-axis. In machine learning we can call intercepts bias. Bias offsets
all predictions that we make.
4 Chapter 1. Linear Regression

Our algorithm will try to learn the correct values for Weight and Bias. By the end of our training, our equation will
approximate the line of best fit.
Code
def predict_sales(radio, weight, bias):

return weight*radio + bias
1.2.2 Cost function
The prediction function is nice, but for our purposes we don’t really need it. What we need is a cost function so we
can start optimizing our weights.
Let’s use MSE (L2) as our cost function. MSE measures the average squared difference between an observation’s
actual and predicted values. The output is a single number representing the cost, or score, associated with our current
set of weights. Our goal is to minimize MSE to improve the accuracy of our model.
Math
Given our simple linear equation 𝑦 = 𝑚𝑥 + 𝑏, we can calculate MSE as:

𝑛
1 ∑︁
𝑀 𝑆𝐸 = (𝑦𝑖 − (𝑚𝑥𝑖 + 𝑏))2
𝑁 𝑖=1
1.2. Simple regression 5

Note:
• 𝑁 is the total number of observations (data points)
∑︀𝑛
• 𝑁1 𝑖=1 is the mean
• 𝑦𝑖 is the actual value of an observation and 𝑚𝑥𝑖 + 𝑏 is our prediction
Code
def cost_function(radio, sales, weight, bias):

companies = len(radio)
total_error = 0.0
for i in range(companies):
total_error += (sales[i] - (weight*radio[i] + bias))**2
return total_error / companies
1.2.3 Gradient descent
To minimize MSE we use Gradient Descent to calculate the gradient of our cost function. Gradient descent consists of
looking at the error that our weight currently gives us, using the derivative of the cost function to find the gradient (The
slope of the cost function using our current weight), and then changing our weight to move in the direction opposite of
the gradient. We need to move in the opposite direction of the gradient since the gradient points up the slope instead
of down it, so we move in the opposite direction to try to decrease our error.
Math
There are two parameters (coefficients) in our cost function we can control: weight 𝑚 and bias 𝑏. Since we need to
consider the impact each one has on the final prediction, we use partial derivatives. To find the partial derivatives, we
use the Chain rule. We need the chain rule because (𝑦 − (𝑚𝑥 + 𝑏))2 is really 2 nested functions: the inner function
𝑦 − (𝑚𝑥 + 𝑏) and the outer function 𝑥2 .
Returning to our cost function:
𝑛
1 ∑︁
𝑓 (𝑚, 𝑏) = (𝑦𝑖 − (𝑚𝑥𝑖 + 𝑏))2
𝑁 𝑖=1
Using the following:
(𝑦𝑖 − (𝑚𝑥𝑖 + 𝑏))2 = 𝐴(𝐵(𝑚, 𝑏))
We can split the derivative into

𝐴(𝑥) = 𝑥2
𝑑𝑓
= 𝐴′ (𝑥) = 2𝑥
𝑑𝑥
and
𝐵(𝑚, 𝑏) = 𝑦𝑖 − (𝑚𝑥𝑖 + 𝑏) = 𝑦𝑖 − 𝑚𝑥𝑖 − 𝑏
𝑑𝑥
= 𝐵 ′ (𝑚) = 0 − 𝑥𝑖 − 0 = −𝑥𝑖
𝑑𝑚
𝑑𝑥
= 𝐵 ′ (𝑏) = 0 − 0 − 1 = −1
𝑑𝑏

And then using the Chain rule which states:

𝑑𝑓 𝑑𝑓 𝑑𝑥
=
𝑑𝑚 𝑑𝑥 𝑑𝑚
𝑑𝑓 𝑑𝑓 𝑑𝑥
=
𝑑𝑏 𝑑𝑥 𝑑𝑏
We then plug in each of the parts to get the following derivatives
𝑑𝑓
= 𝐴′ (𝐵(𝑚, 𝑓 ))𝐵 ′ (𝑚) = 2(𝑦𝑖 − (𝑚𝑥𝑖 + 𝑏)) · −𝑥𝑖
𝑑𝑚
𝑑𝑓
= 𝐴′ (𝐵(𝑚, 𝑓 ))𝐵 ′ (𝑏) = 2(𝑦𝑖 − (𝑚𝑥𝑖 + 𝑏)) · −1
𝑑𝑏
We can calculate the gradient of this cost function as:
[︂ 𝑑𝑓 ]︂ [︂ 1 ∑︀ ]︂
𝑓 ′ (𝑚, 𝑏) = 𝑑𝑚 = 𝑁 ∑︀−𝑥𝑖 · 2(𝑦𝑖 − (𝑚𝑥𝑖 + 𝑏)) (1.1)
𝑑𝑓 1
𝑑𝑏 𝑁 −1 · 2(𝑦𝑖 − (𝑚𝑥𝑖 + 𝑏))
[︂ 1 ∑︀ ]︂
𝑁 −2𝑥𝑖 (𝑦𝑖 − (𝑚𝑥𝑖 + 𝑏))
= 1
∑︀ (1.2)
𝑁 −2(𝑦𝑖 − (𝑚𝑥𝑖 + 𝑏))
Code
To solve for the gradient, we iterate through our data points using our new weight and bias values and take the average
of the partial derivatives. The resulting gradient tells us the slope of our cost function at our current position (i.e.
weight and bias) and the direction we should update to reduce our cost function (we move in the direction opposite the
gradient). The size of our update is controlled by the learning rate.
def update_weights(radio, sales, weight, bias, learning_rate):
weight_deriv = 0
bias_deriv = 0
companies = len(radio)
for i in range(companies):
# Calculate partial derivatives
# -2x(y - (mx + b))
weight_deriv += -2*radio[i] * (sales[i] - (weight*radio[i] + bias))
# -2(y - (mx + b))

bias_deriv += -2*(sales[i] - (weight*radio[i] + bias))
# We subtract because the derivatives point in direction of steepest ascent

weight -= (weight_deriv / companies) * learning_rate
bias -= (bias_deriv / companies) * learning_rate
return weight, bias
1.2.4 Training
Training a model is the process of iteratively improving your prediction equation by looping through the dataset
multiple times, each time updating the weight and bias values in the direction indicated by the slope of the cost
function (gradient). Training is complete when we reach an acceptable error threshold, or when subsequent training
iterations fail to reduce our cost.
Before training we need to initialize our weights (set default values), set our hyperparameters (learning rate and
number of iterations), and prepare to log our progress over each iteration.

Code
def train(radio, sales, weight, bias, learning_rate, iters):

cost_history = []
for i in range(iters):
weight,bias = update_weights(radio, sales, weight, bias, learning_rate)
#Calculate cost for auditing purposes

cost = cost_function(radio, sales, weight, bias)
cost_history.append(cost)
# Log Progress
if i % 10 == 0:
print "iter={:d} weight={:.2f} bias={:.4f} cost={:.2}".format(i,
˓→weight, bias, cost)
return weight, bias, cost_history
1.2.5 Model evaluation
If our model is working, we should see our cost decrease after every iteration.
Logging
iter=1 weight=.03 bias=.0014 cost=197.25


Visualizing




Cost history
1.2.6 Summary
By learning the best values for weight (.46) and bias (.25), we now have an equation that predicts future sales based
on radio advertising investment.
𝑆𝑎𝑙𝑒𝑠 = .46𝑅𝑎𝑑𝑖𝑜 + .025
How would our model perform in the real world? I’ll let you think about it :)
1.3 Multivariable regression
Let’s say we are given data on TV, radio, and newspaper advertising spend for a list of companies, and our goal is to
predict sales in terms of units sold.
Company TV Radio News Units

Amazon 230.1 37.8 69.1 22.1
Google 44.5 39.3 23.1 10.4
Facebook 17.2 45.9 34.7 18.3
Apple 151.5 41.3 13.2 18.5
1.3. Multivariable regression 13

1.3.1 Growing complexity
As the number of features grows, the complexity of our model increases and it becomes increasingly difficult to
visualize, or even comprehend, our data.
One solution is to break the data apart and compare 1-2 features at a time. In this example we explore how Radio and
TV investment impacts Sales.
1.3.2 Normalization
As the number of features grows, calculating gradient takes longer to compute. We can speed this up by “normalizing”
our input data to ensure all values are within the same range. This is especially important for datasets with high
standard deviations or differences in the ranges of the attributes. Our goal now will be to normalize our features so
they are all in the range -1 to 1.
Code
For each feature column {

#1 Subtract the mean of the column (mean normalization)
#2 Divide by the range of the column (feature scaling)
}
Our input is a 200 x 3 matrix containing TV, Radio, and Newspaper data. Our output is a normalized matrix of the
same shape with all values between -1 and 1.

def normalize(features):
**
features - (200, 3)
features.T - (3, 200)
We transpose the input matrix, swapping

cols and rows to make vector math easier
**
for feature in features.T:

fmean = np.mean(feature)
frange = np.amax(feature) - np.amin(feature)
#Vector Subtraction
feature -= fmean
#Vector Division
feature /= frange
return features
Note: Matrix math. Before we continue, it’s important to understand basic Linear Algebra concepts as well as
numpy functions like numpy.dot().
Our predict function outputs an estimate of sales given our current weights (coefficients) and a company’s TV, radio,
and newspaper spend. Our model will try to identify weight values that most reduce our cost function.
𝑆𝑎𝑙𝑒𝑠 = 𝑊1 𝑇 𝑉 + 𝑊2 𝑅𝑎𝑑𝑖𝑜 + 𝑊3 𝑁 𝑒𝑤𝑠𝑝𝑎𝑝𝑒𝑟
def predict(features, weights):

**
features - (200, 3)
weights - (3, 1)
predictions - (200,1)
**
predictions = np.dot(features, weights)
return predictions
1.3.4 Initialize weights
W1 = 0.0
W2 = 0.0
W3 = 0.0
weights = np.array([
[W1],
[W2],
[W3]
])

1.3.5 Cost function
Now we need a cost function to audit how our model is performing. The math is the same, except we swap the 𝑚𝑥 + 𝑏
expression for 𝑊1 𝑥1 + 𝑊2 𝑥2 + 𝑊3 𝑥3 . We also divide the expression by 2 to make derivative calculations simpler.
𝑛
1 ∑︁
𝑀 𝑆𝐸 = (𝑦𝑖 − (𝑊1 𝑥1 + 𝑊2 𝑥2 + 𝑊3 𝑥3 ))2
2𝑁 𝑖=1
def cost_function(features, targets, weights):

**
features:(200,3)
targets: (200,1)
weights:(3,1)
returns average squared error among predictions
**
N = len(targets)
predictions = predict(features, weights)
# Matrix math lets use do this without looping

sq_error = (predictions - targets)**2
# Return average squared error among predictions

return 1.0/(2*N) * sq_error.sum()
Again using the Chain rule we can compute the gradient–a vector of partial derivatives describing the slope of the cost
function for each weight.
𝑓 ′ (𝑊1 ) = −𝑥1 (𝑦 − (𝑊1 𝑥1 + 𝑊2 𝑥2 + 𝑊3 𝑥3 )) (1.3)

𝑓 ′ (𝑊2 ) = −𝑥2 (𝑦 − (𝑊1 𝑥1 + 𝑊2 𝑥2 + 𝑊3(1.4)
𝑥3 ))
𝑓 ′ (𝑊3 ) = −𝑥3 (𝑦 − (𝑊1 𝑥1 + 𝑊2 𝑥2 + 𝑊3(1.5)
𝑥3 ))
def update_weights(features, targets, weights, lr):

'''
Features:(200, 3)
Targets: (200, 1)
Weights:(3, 1)
'''
#Extract our features

x1 = features[:,0]
x2 = features[:,1]
x3 = features[:,2]
# Use matrix cross product (*) to simultaneously

# calculate the derivative for each weight
d_w1 = -x1*(targets - predictions)
(continues on next page)

(continued from previous page)

# Multiply the mean derivative by the learning rate
# and subtract from our weights (remember gradient points in direction of
˓→steepest ASCENT)
weights[0][0] -= (lr * np.mean(d_w1))

return weights
And that’s it! Multivariate linear regression.
1.3.7 Simplifying with matrices
The gradient descent code above has a lot of duplication. Can we improve it somehow? One way to refactor would be
to loop through our features and weights–allowing our function to handle any number of features. However there is
another even better technique: vectorized gradient descent.
Math
We use the same formula as above, but instead of operating on a single feature at a time, we use matrix multiplication
to operative on all features and weights simultaneously. We replace the 𝑥𝑖 terms with a single feature matrix 𝑋.
𝑔𝑟𝑎𝑑𝑖𝑒𝑛𝑡 = −𝑋(𝑡𝑎𝑟𝑔𝑒𝑡𝑠 − 𝑝𝑟𝑒𝑑𝑖𝑐𝑡𝑖𝑜𝑛𝑠)
Code
X = [
[x1, x2, x3]
[x1, x2, x3]
.
.
.
[x1, x2, x3]
]
targets = [
[1],
[2],
[3]
]
def update_weights_vectorized(X, targets, weights, lr):

**
gradient = X.T * (predictions - targets) / N
X: (200, 3)
Targets: (200, 1)
Weights: (3, 1)
**
companies = len(X)
#1 - Get Predictions
predictions = predict(X, weights)

#2 - Calculate error/loss
error = targets - predictions
#3 Transpose features from (200, 3) to (3, 200)

# So we can multiply w the (200,1) error matrix.
# Returns a (3,1) matrix holding 3 partial derivatives --
# one for each feature -- representing the aggregate
# slope of the cost function across all observations
gradient = np.dot(-X.T, error)
#4 Take the average error derivative for each feature

gradient /= companies
#5 - Multiply the gradient by our learning rate

gradient *= lr
#6 - Subtract from our weights to minimize cost

weights -= gradient
return weights
1.3.8 Bias term
Our train function is the same as for simple linear regression, however we’re going to make one final tweak before
running: add a bias term to our feature matrix.
In our example, it’s very unlikely that sales would be zero if companies stopped advertising. Possible reasons for this
might include past advertising, existing customer relationships, retail locations, and salespeople. A bias term will help
us capture this base case.
Code
Below we add a constant 1 to our features matrix. By setting this value to 1, it turns our bias term into a constant.
bias = np.ones(shape=(len(features),1))
features = np.append(bias, features, axis=1)
After training our model through 1000 iterations with a learning rate of .0005, we finally arrive at a set of weights we
can use to make predictions:
𝑆𝑎𝑙𝑒𝑠 = 4.7𝑇 𝑉 + 3.5𝑅𝑎𝑑𝑖𝑜 + .81𝑁 𝑒𝑤𝑠𝑝𝑎𝑝𝑒𝑟 + 13.9
Our MSE cost dropped from 110.86 to 6.25.

References


CHAPTER 2
Gradient Descent
Gradient descent is an optimization algorithm used to minimize some function by iteratively moving in the direction
of steepest descent as defined by the negative of the gradient. In machine learning, we use gradient descent to update
the parameters of our model. Parameters refer to coefficients in Linear Regression and weights in neural networks.
2.1 Introduction
Consider the 3-dimensional graph below in the context of a cost function. Our goal is to move from the mountain in
the top right corner (high cost) to the dark blue sea in the bottom left (low cost). The arrows represent the direction of
steepest descent (negative gradient) from any given point–the direction that decreases the cost function as quickly as
possible. Source
Starting at the top of the mountain, we take our first step downhill in the direction specified by the negative gradient.
Next we recalculate the negative gradient (passing in the coordinates of our new point) and take another step in the
direction it specifies. We continue this process iteratively until we get to the bottom of our graph, or to a point where
we can no longer move downhill–a local minimum. image source.
21
2.2 Learning rate
The size of these steps is called the learning rate. With a high learning rate we can cover more ground each step, but
we risk overshooting the lowest point since the slope of the hill is constantly changing. With a very low learning rate,
we can confidently move in the direction of the negative gradient since we are recalculating it so frequently. A low
learning rate is more precise, but calculating the gradient is time-consuming, so it will take us a very long time to get
to the bottom.
2.3 Cost function
A Loss Functions tells us “how good” our model is at making predictions for a given set of parameters. The cost
function has its own curve and its own gradients. The slope of this curve tells us how to update our parameters to make
the model more accurate.
2.4 Step-by-step
Now let’s run gradient descent using our new cost function. There are two parameters in our cost function we can
control: 𝑚 (weight) and 𝑏 (bias). Since we need to consider the impact each one has on the final prediction, we need
to use partial derivatives. We calculate the partial derivatives of the cost function with respect to each parameter and
store the results in a gradient.
Math
Given the cost function:

𝑛
1 ∑︁
𝑓 (𝑚, 𝑏) = (𝑦𝑖 − (𝑚𝑥𝑖 + 𝑏))2
𝑁 𝑖=1
22 Chapter 2. Gradient Descent

The gradient can be calculated as:

𝑑𝑓
[︂ ]︂ [︂ 1 ∑︀ ]︂
𝑓 ′ (𝑚, 𝑏) = 𝑑𝑚 = 𝑁 ∑︀−2𝑥𝑖 (𝑦𝑖 − (𝑚𝑥𝑖 + 𝑏))
𝑑𝑓 1
𝑑𝑏 𝑁 −2(𝑦𝑖 − (𝑚𝑥𝑖 + 𝑏))
To solve for the gradient, we iterate through our data points using our new 𝑚 and 𝑏 values and compute the partial
derivatives. This new gradient tells us the slope of our cost function at our current position (current parameter values)
and the direction we should move to update our parameters. The size of our update is controlled by the learning rate.
Code
def update_weights(m, b, X, Y, learning_rate):

m_deriv = 0
b_deriv = 0
N = len(X)
for i in range(N):
# Calculate partial derivatives
# -2x(y - (mx + b))
m_deriv += -2*X[i] * (Y[i] - (m*X[i] + b))
# -2(y - (mx + b))

b_deriv += -2*(Y[i] - (m*X[i] + b))
# We subtract because the derivatives point in direction of steepest ascent

m -= (m_deriv / float(N)) * learning_rate
b -= (b_deriv / float(N)) * learning_rate
return m, b
References
2.4. Step-by-step 23
24 Chapter 2. Gradient Descent

CHAPTER 3
Logistic Regression
• Introduction
– Comparison to linear regression
– Types of logistic regression
• Binary logistic regression
– Sigmoid activation
– Decision boundary
– Cost function
– Mapping probabilities to classes
– Training
• Multiclass logistic regression
– Procedure
– Softmax activation
– Scikit-Learn example
3.1 Introduction
Logistic regression is a classification algorithm used to assign observations to a discrete set of classes. Unlike linear
regression which outputs continuous number values, logistic regression transforms its output using the logistic sigmoid
25
function to return a probability value which can then be mapped to two or more discrete classes.
3.1.1 Comparison to linear regression
Given data on time spent studying and exam scores. Linear Regression and logistic regression can predict different
things:
• Linear Regression could help us predict the student’s test score on a scale of 0 - 100. Linear regression
predictions are continuous (numbers in a range).
• Logistic Regression could help use predict whether the student passed or failed. Logistic regression predictions
are discrete (only specific values or categories are allowed). We can also view probability scores underlying the
model’s classifications.
3.1.2 Types of logistic regression
• Binary (Pass/Fail)
• Multi (Cats, Dogs, Sheep)
• Ordinal (Low, Medium, High)
3.2 Binary logistic regression
Say we’re given data on student exam results and our goal is to predict whether a student will pass or fail based on
number of hours slept and hours spent studying. We have two features (hours slept, hours studied) and two classes:
passed (1) and failed (0).
Studied Slept Passed

4.85 9.63 1
8.62 3.23 0
5.43 8.23 1
9.21 6.34 0
Graphically we could represent our data with a scatter plot.
26 Chapter 3. Logistic Regression

3.2.1 Sigmoid activation
In order to map predicted values to probabilities, we use the sigmoid function. The function maps any real value into
another value between 0 and 1. In machine learning, we use sigmoid to map predictions to probabilities.
Math
1
𝑆(𝑧) =
1 + 𝑒−𝑧
Note:
• 𝑠(𝑧) = output between 0 and 1 (probability estimate)
• 𝑧 = input to the function (your algorithm’s prediction e.g. mx + b)
• 𝑒 = base of natural log
3.2. Binary logistic regression 27

Graph
Code
def sigmoid(z):
return 1.0 / (1 + np.exp(-z))
3.2.2 Decision boundary
Our current prediction function returns a probability score between 0 and 1. In order to map this to a discrete class
(true/false, cat/dog), we select a threshold value or tipping point above which we will classify values into class 1 and
below which we classify values into class 2.
𝑝 ≥ 0.5, 𝑐𝑙𝑎𝑠𝑠 = 1
𝑝 < 0.5, 𝑐𝑙𝑎𝑠𝑠 = 0

For example, if our threshold was .5 and our prediction function returned .7, we would classify this observation as
positive. If our prediction was .2 we would classify the observation as negative. For logistic regression with multiple
classes we could select the class with the highest predicted probability.
Using our knowledge of sigmoid functions and decision boundaries, we can now write a prediction function. A
prediction function in logistic regression returns the probability of our observation being positive, True, or “Yes”. We
call this class 1 and its notation is 𝑃 (𝑐𝑙𝑎𝑠𝑠 = 1). As the probability gets closer to 1, our model is more confident that
the observation is in class 1.
Math
Let’s use the same multiple linear regression equation from our linear regression tutorial.
𝑧 = 𝑊0 + 𝑊1 𝑆𝑡𝑢𝑑𝑖𝑒𝑑 + 𝑊2 𝑆𝑙𝑒𝑝𝑡
This time however we will transform the output using the sigmoid function to return a probability value between 0 and
1.
1
𝑃 (𝑐𝑙𝑎𝑠𝑠 = 1) =
1 + 𝑒−𝑧
If the model returns .4 it believes there is only a 40% chance of passing. If our decision boundary was .5, we would
categorize this observation as “Fail.””
Code
We wrap the sigmoid function over the same prediction function we used in multiple linear regression

def predict(features, weights):

'''
Returns 1D array of probabilities
that the class label == 1
'''
z = np.dot(features, weights)
return sigmoid(z)
3.2.4 Cost function
Unfortunately we can’t (or at least shouldn’t) use the same cost function MSE (L2) as we did for linear regression.
Why? There is a great math explanation in chapter 3 of Michael Neilson’s deep learning book5 , but for now I’ll simply
say it’s because our prediction function is non-linear (due to sigmoid transform). Squaring this prediction as we do
in MSE results in a non-convex function with many local minimums. If our cost function has many local minimums,
gradient descent may not find the optimal global minimum.
Math
Instead of Mean Squared Error, we use a cost function called Cross-Entropy, also known as Log Loss. Cross-entropy
loss can be divided into two separate cost functions: one for 𝑦 = 1 and one for 𝑦 = 0.
The benefits of taking the logarithm reveal themselves when you look at the cost function graphs for y=1 and y=0.
These smooth monotonic functions7 (always increasing or always decreasing) make it easy to calculate the gradient
and minimize cost. Image from Andrew Ng’s slides on logistic regression1 .
5 http://neuralnetworksanddeeplearning.com/chap3.html
7 https://en.wikipedia.org/wiki/Monotoniconotonic_function
1 http://www.holehouse.org/mlclass/06_Logistic_Regression.html

The key thing to note is the cost function penalizes confident and wrong predictions more than it rewards confident
and right predictions! The corollary is increasing prediction accuracy (closer to 0 or 1) has diminishing returns on
reducing cost due to the logistic nature of our cost function.
Above functions compressed into one
Multiplying by 𝑦 and (1 − 𝑦) in the above equation is a sneaky trick that let’s us use the same equation to solve for
both y=1 and y=0 cases. If y=0, the first side cancels out. If y=1, the second side cancels out. In both cases we only
perform the operation we need to perform.
Vectorized cost function
Code
def cost_function(features, labels, weights):

'''
Using Mean Absolute Error
Features:(100,3)
Labels: (100,1)
Weights:(3,1)


Returns 1D matrix of predictions
Cost = (labels*log(predictions) + (1-labels)*log(1-predictions) ) / len(labels)
'''
observations = len(labels)
#Take the error when label=1

class1_cost = -labels*np.log(predictions)
#Take the error when label=0

class2_cost = (1-labels)*np.log(1-predictions)
#Take the sum of both costs

cost = class1_cost - class2_cost
#Take the average cost

cost = cost.sum() / observations
return cost
To minimize our cost, we use Gradient Descent just like before in Linear Regression. There are other more sophisti-
cated optimization algorithms out there such as conjugate gradient like BFGS, but you don’t have to worry about these.
Machine learning libraries like Scikit-learn hide their implementations so you can focus on more interesting things!
Math
One of the neat properties of the sigmoid function is its derivative is easy to calculate. If you’re curious, there is a
good walk-through derivation on stack overflow6 . Michael Neilson also covers the topic in chapter 3 of his book.
𝑠′ (𝑧) = 𝑠(𝑧)(1 − 𝑠(𝑧)) (3.1)
Which leads to an equally beautiful and convenient cost function derivative:
𝐶 ′ = 𝑥(𝑠(𝑧) − 𝑦)
Note:
• 𝐶 ′ is the derivative of cost with respect to weights
• 𝑦 is the actual class label (0 or 1)
• 𝑠(𝑧) is your model’s prediction
• 𝑥 is your feature or feature vector.
Notice how this gradient is the same as the MSE (L2) gradient, the only difference is the hypothesis function.
6 http://math.stackexchange.com/questions/78575/derivative-of-sigmoid-function-sigma-x-frac11e-x

Pseudocode
Repeat {
1. Calculate gradient average

2. Multiply by learning rate
3. Subtract from weights
Code
def update_weights(features, labels, weights, lr):

'''
Vectorized Gradient Descent
Features:(200, 3)
Labels: (200, 1)
Weights:(3, 1)
'''
N = len(features)
#1 - Get Predictions
#2 Transpose features from (200, 3) to (3, 200)

# So we can multiply w the (200,1) cost matrix.
# Returns a (3,1) matrix holding 3 partial derivatives --
# one for each feature -- representing the aggregate
# slope of the cost function across all observations
gradient = np.dot(features.T, predictions - labels)
#3 Take the average cost derivative for each feature

gradient /= N
#4 - Multiply the gradient by our learning rate

gradient *= lr
#5 - Subtract from our weights to minimize cost

weights -= gradient
return weights
3.2.6 Mapping probabilities to classes
The final step is assign class labels (0 or 1) to our predicted probabilities.
Decision boundary
def decision_boundary(prob):
return 1 if prob >= .5 else 0

Convert probabilities to classes
def classify(predictions):
'''
input - N element array of predictions between 0 and 1
output - N element array of 0s (False) and 1s (True)
'''
decision_boundary = np.vectorize(decision_boundary)
return decision_boundary(predictions).flatten()
Example output
Probabilities = [ 0.967, 0.448, 0.015, 0.780, 0.978, 0.004]

Classifications = [1, 0, 0, 1, 1, 0]
3.2.7 Training
Our training code is the same as we used for linear regression.
def train(features, labels, weights, lr, iters):

cost_history = []
for i in range(iters):
weights = update_weights(features, labels, weights, lr)
#Calculate error for auditing purposes

cost = cost_function(features, labels, weights)
cost_history.append(cost)
# Log Progress
if i % 1000 == 0:
print "iter: "+str(i) + " cost: "+str(cost)
return weights, cost_history
If our model is working, we should see our cost decrease after every iteration.
iter: 0 cost: 0.635

iter: 1000 cost: 0.302
iter: 2000 cost: 0.264
Final cost: 0.2487. Final weights: [-8.197, .921, .738]

Cost history
Accuracy
Accuracy measures how correct our predictions were. In this case we simply compare predicted labels to true labels
and divide by the total.
def accuracy(predicted_labels, actual_labels):

diff = predicted_labels - actual_labels
return 1.0 - (float(np.count_nonzero(diff)) / len(diff))
Decision boundary
Another helpful technique is to plot the decision boundary on top of our predictions to see how our labels compare to
the actual labels. This involves plotting our predicted probabilities and coloring them with their true labels.

Code to plot the decision boundary
def plot_decision_boundary(trues, falses):

fig = plt.figure()
ax = fig.add_subplot(111)
no_of_preds = len(trues) + len(falses)
ax.scatter([i for i in range(len(trues))], trues, s=25, c='b', marker="o", label=

˓→'Trues')
ax.scatter([i for i in range(len(falses))], falses, s=25, c='r', marker="s",
˓→label='Falses')
plt.legend(loc='upper right');
ax.set_title("Decision Boundary")
ax.set_xlabel('N/2')
ax.set_ylabel('Predicted Probability')
plt.axhline(.5, color='black')
plt.show()
3.3 Multiclass logistic regression
Instead of 𝑦 = 0, 1 we will expand our definition so that 𝑦 = 0, 1...𝑛. Basically we re-run binary classification multiple
times, once for each class.

3.3.1 Procedure
1. Divide the problem into n+1 binary classification problems (+1 because the index starts at 0?).
2. For each class. . .
3. Predict the probability the observations are in that single class.
4. prediction = <math>max(probability of the classes)
For each sub-problem, we select one class (YES) and lump all the others into a second class (NO). Then we take the
class with the highest predicted value.
3.3.2 Softmax activation
The softmax function (softargmax or normalized exponential function) is a function that takes as input a vector of
K real numbers, and normalizes it into a probability distribution consisting of K probabilities proportional to the
exponentials of the input numbers. That is, prior to applying softmax, some vector components could be negative, or
greater than one; and might not sum to 1; but after applying softmax, each component will be in the interval [ 0 , 1
] , and the components will add up to 1, so that they can be interpreted as probabilities. The standard (unit) softmax
function is defined by the formula
𝑒𝑧(𝑖)
(𝑧𝑖 ) = ∑︀𝐾 𝑓 𝑜𝑟 𝑖 = 1, ., ., ., 𝐾 𝑎𝑛𝑑 𝑧 = 𝑧1 , ., ., ., 𝑧𝐾 (3.2)
𝑧(𝑗)
𝑗=1 𝑒
In words: we apply the standard exponential function to each element 𝑧𝑖 of the input vector 𝑧 and normalize these
values by dividing by the sum of all these exponentials; this normalization ensures that the sum of the components of
the output vector (𝑧) is 1.9
3.3.3 Scikit-Learn example
Let’s compare our performance to the LogisticRegression model provided by scikit-learn8 .
import sklearn
from sklearn.linear_model import LogisticRegression
from sklearn.cross_validation import train_test_split
# Normalize grades to values between 0 and 1 for more efficient computation

normalized_range = sklearn.preprocessing.MinMaxScaler(feature_range=(-1,1))
# Extract Features + Labels

labels.shape = (100,) #scikit expects this
features = normalized_range.fit_transform(features)
# Create Test/Train
features_train,features_test,labels_train,labels_test = train_test_split(features,
˓→labels,test_size=0.4)
# Scikit Logistic Regression

scikit_log_reg = LogisticRegression()
scikit_log_reg.fit(features_train,labels_train)
#Score is Mean Accuracy

9 https://en.wikipedia.org/wiki/Softmax_function
8 http://scikit-learn.org/stable/modules/linear_model.html#logistic-regression>
3.3. Multiclass logistic regression 37


scikit_score = clf.score(features_test,labels_test)
print 'Scikit score: ', scikit_score
#Our Mean Accuracy

observations, features, labels, weights = run()
probabilities = predict(features, weights).flatten()
classifications = classifier(probabilities)
our_acc = accuracy(classifications,labels.flatten())
print 'Our score: ',our_acc
Scikit score: 0.88. Our score: 0.89
References

CHAPTER 4
Glossary
Definitions of common machine learning terms.

Accuracy Percentage of correct predictions made by the model.
Algorithm A method, function, or series of instructions used to generate a machine learning model. Examples include
linear regression, decision trees, support vector machines, and neural networks.
Attribute A quality describing an observation (e.g. color, size, weight). In Excel terms, these are column headers.
Bias metric What is the average difference between your predictions and the correct value for that observation?
• Low bias could mean every prediction is correct. It could also mean half of your predictions are above
their actual values and half are below, in equal proportion, resulting in low average difference.
• High bias (with low variance) suggests your model may be underfitting and you’re using the wrong archi-
tecture for the job.
Bias term Allow models to represent patterns that do not pass through the origin. For example, if all my features
were 0, would my output also be zero? Is it possible there is some base value upon which my features have an
effect? Bias terms typically accompany weights and are attached to neurons or filters.
Categorical Variables Variables with a discrete set of possible values. Can be ordinal (order matters) or nominal
(order doesn’t matter).
Classification Predicting a categorical output.
• Binary classification predicts one of two possible outcomes (e.g. is the email spam or not spam?)
• Multi-class classification predicts one of multiple possible outcomes (e.g. is this a photo of a cat, dog,
horse or human?)
Classification Threshold The lowest probability value at which we’re comfortable asserting a positive classification.
For example, if the predicted probability of being diabetic is > 50%, return True, otherwise return False.
Clustering Unsupervised grouping of data into buckets.
Confusion Matrix Table that describes the performance of a classification model by grouping predictions into 4
categories.
• True Positives: we correctly predicted they do have diabetes
39
• True Negatives: we correctly predicted they don’t have diabetes

• False Positives: we incorrectly predicted they do have diabetes (Type I error)
• False Negatives: we incorrectly predicted they don’t have diabetes (Type II error)
Continuous Variables Variables with a range of possible values defined by a number scale (e.g. sales, lifespan).
Convergence A state reached during the training of a model when the loss changes very little between each iteration.
Deduction A top-down approach to answering questions or solving problems. A logic technique that starts with a
theory and tests that theory with observations to derive a conclusion. E.g. We suspect X, but we need to test our
hypothesis before coming to any conclusions.
Deep Learning Deep Learning is derived from one machine learning algorithm called perceptron or multi layer per-
ceptron that gain more and more attention nowadays because of its success in different fields like, computer
vision to signal processing and medical diagnosis to self-driving cars. As all other AI algorithms deep learn-
ing is from decades, but now today we have more and more data and cheap computing power that make this
algorithm really powerful to achieve state of the art accuracy. In modern world this algorithm knowns as arti-
ficial neural network. deep learning is much more than traditional artificial neural network. But it was highly
influenced by machine learning’s neural network and perceptron network.
Dimension Dimension for machine learning and data scientist is differ from physics, here Dimension of data means
how much feature you have in you data ocean(data-set). e.g in case of object detection application, flatten image
size and color channel(e.g 28*28*3) is a feature of the input set. In case of house price prediction (maybe) house
size is the data-set so we call it 1 dimentional data.
Epoch An epoch describes the number of times the algorithm sees the entire data set.
Extrapolation Making predictions outside the range of a dataset. E.g. My dog barks, so all dogs must bark. In
machine learning we often run into trouble when we extrapolate outside the range of our training data.
False Positive Rate Defined as
𝐹 𝑎𝑙𝑠𝑒𝑃 𝑜𝑠𝑖𝑡𝑖𝑣𝑒𝑠
𝐹 𝑃 𝑅 = 1 − 𝑆𝑝𝑒𝑐𝑖𝑓 𝑖𝑐𝑖𝑡𝑦 =
𝐹 𝑎𝑙𝑠𝑒𝑃 𝑜𝑠𝑖𝑡𝑖𝑣𝑒𝑠 + 𝑇 𝑟𝑢𝑒𝑁 𝑒𝑔𝑎𝑡𝑖𝑣𝑒𝑠
The False Positive Rate forms the x-axis of the ROC curve.
Feature With respect to a dataset, a feature represents an attribute and value combination. Color is an attribute.
“Color is blue” is a feature. In Excel terms, features are similar to cells. The term feature has other definitions
in different contexts.
Feature Selection Feature selection is the process of selecting relevant features from a data-set for creating a Machine
Learning model.
Feature Vector A list of features describing an observation with multiple attributes. In Excel we call this a row.
Gradient Accumulation A mechanism to split the batch of samples—used for training a neural network—into sev-
eral mini-batches of samples that will be run sequentially. This is used to enable using large batch sizes that
require more GPU memory than available.
Hyperparameters Hyperparameters are higher-level properties of a model such as how fast it can learn (learning
rate) or complexity of a model. The depth of trees in a Decision Tree or number of hidden layers in a Neural
Networks are examples of hyper parameters.
Induction A bottoms-up approach to answering questions or solving problems. A logic technique that goes from
observations to theory. E.g. We keep observing X, so we infer that Y must be True.
Instance A data point, row, or sample in a dataset. Another term for observation.
Label The “answer” portion of an observation in supervised learning. For example, in a dataset used to classify
flowers into different species, the features might include the petal length and petal width, while the label would
be the flower’s species.
40 Chapter 4. Glossary
Learning Rate The size of the update steps to take during optimization loops like Gradient Descent. With a high
learning rate we can cover more ground each step, but we risk overshooting the lowest point since the slope of
the hill is constantly changing. With a very low learning rate, we can confidently move in the direction of the
negative gradient since we are recalculating it so frequently. A low learning rate is more precise, but calculating
the gradient is time-consuming, so it will take us a very long time to get to the bottom.
Loss Loss = true_value(from data-set)- predicted value(from ML-model) The lower the loss, the better a model (un-
less the model has over-fitted to the training data). The loss is calculated on training and validation and its
interpretation is how well the model is doing for these two sets. Unlike accuracy, loss is not a percentage. It is a
summation of the errors made for each example in training or validation sets.
Machine Learning Mitchell (1997) provides a succinct definition: “A computer program is said to learn from expe-
rience E with respect to some class of tasks T and performance measure P, if its performance at tasks in T, as
measured by P, improves with experience E.” In simple language machine learning is a field in which human
made algorithms have an ability learn by itself or predict future for unseen data.
Model A data structure that stores a representation of a dataset (weights and biases). Models are created/learned when
you train an algorithm on a dataset.
Neural Networks Neural Networks are mathematical algorithms modeled after the brain’s architecture, designed to
recognize patterns and relationships in data.
Normalization Restriction of the values of weights in regression to avoid overfitting and improving computation
speed.
Noise Any irrelevant information or randomness in a dataset which obscures the underlying pattern.
Null Accuracy Baseline accuracy that can be achieved by always predicting the most frequent class (“B has the
highest frequency, so lets guess B every time”).
Observation A data point, row, or sample in a dataset. Another term for instance.
Outlier An observation that deviates significantly from other observations in the dataset.
Overfitting Overfitting occurs when your model learns the training data too well and incorporates details and noise
specific to your dataset. You can tell a model is overfitting when it performs great on your training/validation
set, but poorly on your test set (or new real-world data).
Parameters Parameters are properties of training data learned by training a machine learning model or classifier.
They are adjusted using optimization algorithms and unique to each experiment.
Examples of parameters include:
• weights in an artificial neural network
• support vectors in a support vector machine
• coefficients in a linear or logistic regression
Precision In the context of binary classification (Yes/No), precision measures the model’s performance at classifying
positive observations (i.e. “Yes”). In other words, when a positive value is predicted, how often is the prediction
correct? We could game this metric by only returning positive for the single observation we are most confident
in.
𝑇 𝑟𝑢𝑒𝑃 𝑜𝑠𝑖𝑡𝑖𝑣𝑒𝑠
𝑃 =
𝑇 𝑟𝑢𝑒𝑃 𝑜𝑠𝑖𝑡𝑖𝑣𝑒𝑠 + 𝐹 𝑎𝑙𝑠𝑒𝑃 𝑜𝑠𝑖𝑡𝑖𝑣𝑒𝑠
Recall Also called sensitivity. In the context of binary classification (Yes/No), recall measures how “sensitive” the
classifier is at detecting positive instances. In other words, for all the true observations in our sample, how many
did we “catch.” We could game this metric by always classifying observations as positive.
𝑅=
𝑇 𝑟𝑢𝑒𝑃 𝑜𝑠𝑖𝑡𝑖𝑣𝑒𝑠 + 𝐹 𝑎𝑙𝑠𝑒𝑁 𝑒𝑔𝑎𝑡𝑖𝑣𝑒𝑠
41
Recall vs Precision Say we are analyzing Brain scans and trying to predict whether a person has a tumor (True) or
not (False). We feed it into our model and our model starts guessing.
• Precision is the % of True guesses that were actually correct! If we guess 1 image is True out of 100
images and that image is actually True, then our precision is 100%! Our results aren’t helpful however
because we missed 10 brain tumors! We were super precise when we tried, but we didn’t try hard enough.
• Recall, or Sensitivity, provides another lens which with to view how good our model is. Again let’s say
there are 100 images, 10 with brain tumors, and we correctly guessed 1 had a brain tumor. Precision is
100%, but recall is 10%. Perfect recall requires that we catch all 10 tumors!
Regression Predicting a continuous output (e.g. price, sales).
Regularization Regularization is a technique utilized to combat the overfitting problem. This is achieved by adding
a complexity term to the loss function that gives a bigger loss for more complex models
Reinforcement Learning Training a model to maximize a reward via iterative trial and error.
ROC (Receiver Operating Characteristic) Curve A plot of the true positive rate against the false positive rate at
all classification thresholds. This is used to evaluate the performance of a classification model at different
classification thresholds. The area under the ROC curve can be interpreted as the probability that the model
correctly distinguishes between a randomly chosen positive observation (e.g. “spam”) and a randomly chosen
negative observation (e.g. “not spam”).
Segmentation It is the process of partitioning a data set into multiple distinct sets. This separation is done such that
the members of the same set are similar to each otherand different from the members of other sets.
Specificity In the context of binary classification (Yes/No), specificity measures the model’s performance at classi-
fying negative observations (i.e. “No”). In other words, when the correct label is negative, how often is the
prediction correct? We could game this metric if we predict everything as negative.
𝑇 𝑟𝑢𝑒𝑁 𝑒𝑔𝑎𝑡𝑖𝑣𝑒𝑠
𝑆=
𝑇 𝑟𝑢𝑒𝑁 𝑒𝑔𝑎𝑡𝑖𝑣𝑒𝑠 + 𝐹 𝑎𝑙𝑠𝑒𝑃 𝑜𝑠𝑖𝑡𝑖𝑣𝑒𝑠
Supervised Learning Training a model using a labeled dataset.

Test Set A set of observations used at the end of model training and validation to assess the predictive power of your
model. How generalizable is your model to unseen data?
Training Set A set of observations used to generate machine learning models.
Transfer Learning A machine learning method where a model developed for a task is reused as the starting point
for a model on a second task. In transfer learning, we take the pre-trained weights of an already trained model
(one that has been trained on millions of images belonging to 1000’s of classes, on several high power GPU’s
for several days) and use these already learned features to predict new classes.
True Positive Rate Another term for recall, i.e.
𝑇𝑃𝑅 =
𝑇 𝑟𝑢𝑒𝑃 𝑜𝑠𝑖𝑡𝑖𝑣𝑒𝑠 + 𝐹 𝑎𝑙𝑠𝑒𝑁 𝑒𝑔𝑎𝑡𝑖𝑣𝑒𝑠
The True Positive Rate forms the y-axis of the ROC curve.
Type 1 Error False Positives. Consider a company optimizing hiring practices to reduce false positives in job offers.
A type 1 error occurs when candidate seems good and they hire him, but he is actually bad.
Type 2 Error False Negatives. The candidate was great but the company passed on him.
Underfitting Underfitting occurs when your model over-generalizes and fails to incorporate relevant variations in
your data that would give your model more predictive power. You can tell a model is underfitting when it
performs poorly on both training and test sets.
Universal Approximation Theorem A neural network with one hidden layer can approximate any continuous func-
tion but only for inputs in a specific range. If you train a network on inputs between -2 and 2, then it will work
well for inputs in the same range, but you can’t expect it to generalize to other inputs without retraining the
model or adding more hidden neurons.
Unsupervised Learning Training a model to find patterns in an unlabeled dataset (e.g. clustering).
Validation Set A set of observations used during model training to provide feedback on how well the current param-
eters generalize beyond the training set. If training error decreases but validation error increases, your model is
likely overfitting and you should pause training.
Variance How tightly packed are your predictions for a particular observation relative to each other?
• Low variance suggests your model is internally consistent, with predictions varying little from each other
after every iteration.
• High variance (with low bias) suggests your model may be overfitting and reading too deeply into the
noise found in every training set.
References
43
CHAPTER 5
Calculus
• Introduction
• Derivatives
– Geometric definition
– Taking the derivative
– Step-by-step
– Machine learning use cases
• Chain rule
– How It Works
– Step-by-step
– Multiple functions
• Gradients
– Partial derivatives
– Step-by-step
– Directional derivatives
– Useful properties
• Integrals
– Computing integrals
– Applications of integration
* Computing probabilities
* Expected value
45
* Variance
5.1 Introduction
You need to know some basic calculus in order to understand how functions change over time (derivatives), and to
calculate the total amount of a quantity that accumulates over a time period (integrals). The language of calculus will
allow you to speak precisely about the properties of functions and better understand their behaviour.
Normally taking a calculus course involves doing lots of tedious calculations by hand, but having the power of com-
puters on your side can make the process much more fun. This section describes the key ideas of calculus which you’ll
need to know to understand machine learning concepts.
5.2 Derivatives
A derivative can be defined in two ways:

1. Instantaneous rate of change (Physics)
2. Slope of a line at a specific point (Geometry)
Both represent the same principle, but for our purposes it’s easier to explain using the geometric definition.
5.2.1 Geometric definition
In geometry slope represents the steepness of a line. It answers the question: how much does 𝑦 or 𝑓 (𝑥) change given
a specific change in 𝑥?
Using this definition we can easily calculate the slope between two points. But what if I asked you, instead of the slope
between two points, what is the slope at a single point on the line? In this case there isn’t any obvious “rise-over-run”
to calculate. Derivatives help us answer this question.
A derivative outputs an expression we can use to calculate the instantaneous rate of change, or slope, at a single point
on a line. After solving for the derivative you can use it to calculate the slope at every other point on the line.
5.2.2 Taking the derivative
Consider the graph below, where 𝑓 (𝑥) = 𝑥2 + 3.
46 Chapter 5. Calculus
The slope between (1,4) and (3,12) would be:

𝑦2 − 𝑦1 12 − 4
𝑠𝑙𝑜𝑝𝑒 = = =4
𝑥2 − 𝑥1 3−1
But how do we calculate the slope at point (1,4) to reveal the change in slope at that specific point?
One way would be to find the two nearest points, calculate their slopes relative to 𝑥 and take the average. But calculus
provides an easier, more precise way: compute the derivative. Computing the derivative of a function is essentially
the same as our original proposal, but instead of finding the two closest points, we make up an imaginary point an
infinitesimally small distance away from 𝑥 and compute the slope between 𝑥 and the new point.
In this way, derivatives help us answer the question: how does 𝑓 (𝑥) change if we make a very very tiny increase to x?
In other words, derivatives help estimate the slope between two points that are an infinitesimally small distance away
from each other. A very, very, very small distance, but large enough to calculate the slope.
In math language we represent this infinitesimally small increase using a limit. A limit is defined as the output value
a function approaches as the input value approaches another value. In our case the target value is the specific point at
which we want to calculate slope.
5.2.3 Step-by-step
Calculating the derivative is the same as calculating normal slope, however in this case we calculate the slope between
our point and a point infinitesimally close to it. We use the variable ℎ to represent this infinitesimally distance. Here
are the steps:
1. Given the function:
𝑓 (𝑥) = 𝑥2
2. Increment 𝑥 by a very small value ℎ(ℎ = 𝑥)
𝑓 (𝑥 + ℎ) = (𝑥 + ℎ)2
3. Apply the slope formula
5.2. Derivatives 47
𝑓 (𝑥 + ℎ) − 𝑓 (𝑥)
ℎ
4. Simplify the equation
𝑥2 + 2𝑥ℎ + ℎ2 − 𝑥2
ℎ
2𝑥ℎ + ℎ2
= 2𝑥 + ℎ
ℎ
5. Set ℎ to 0 (the limit as ℎ heads toward 0)
2𝑥 + 0 = 2𝑥
So what does this mean? It means for the function 𝑓 (𝑥) = 𝑥2 , the slope at any point equals 2𝑥. The formula is defined
as:
𝑓 (𝑥 + ℎ) − 𝑓 (𝑥)
lim
ℎ→0 ℎ
Code
Let’s write code to calculate the derivative of any function 𝑓 (𝑥). We test our function works as expected on the input
𝑓 (𝑥) = 𝑥2 producing a value close to the actual derivative 2𝑥.
def get_derivative(func, x):

"""Compute the derivative of `func` at the location `x`."""
h = 0.0001 # step size
return (func(x+h) - func(x)) / h # rise-over-run
def f(x): return x**2 # some test function f(x)=x^2

x = 3 # the location of interest
computed = get_derivative(f, x)
actual = 2*x
computed, actual # = 6.0001, 6 # pretty close if you ask me...
In general it’s preferable to use the math to obtain exact derivative formulas, but keep in mind you can always compute
derivatives numerically by computing the rise-over-run for a “small step” ℎ.
5.2.4 Machine learning use cases
Machine learning uses derivatives in optimization problems. Optimization algorithms like gradient descent use deriva-
tives to decide whether to increase or decrease weights in order to maximize or minimize some objective (e.g. a
model’s accuracy or error functions). Derivatives also help us approximate nonlinear functions as linear functions
(tangent lines), which have constant slopes. With a constant slope we can decide whether to move up or down the
slope (increase or decrease our weights) to get closer to the target value (class label).
5.3 Chain rule
The chain rule is a formula for calculating the derivatives of composite functions. Composite functions are functions
composed of functions inside other function(s).
5.3.1 How It Works
Given a composite function 𝑓 (𝑥) = 𝐴(𝐵(𝑥)), the derivative of 𝑓 (𝑥) equals the product of the derivative of 𝐴 with
respect to 𝐵(𝑥) and the derivative of 𝐵 with respect to 𝑥.
composite function derivative = outer function derivative * inner function derivative
For example, given a composite function 𝑓 (𝑥), where:
𝑓 (𝑥) = ℎ(𝑔(𝑥))
The chain rule tells us that the derivative of 𝑓 (𝑥) equals:
𝑑𝑓 𝑑ℎ 𝑑𝑔
= ·
𝑑𝑥 𝑑𝑔 𝑑𝑥
5.3.2 Step-by-step
Say 𝑓 (𝑥) is composed of two functions ℎ(𝑥) = 𝑥3 and 𝑔(𝑥) = 𝑥2 . And that:
𝑓 (𝑥) = ℎ(𝑔(𝑥)) (5.1)

2 3
= (𝑥 )
(5.2)
(5.3)
The derivative of 𝑓 (𝑥) would equal:
𝑑𝑓 𝑑ℎ 𝑑𝑔
= (5.4)
𝑑𝑥 𝑑𝑔 𝑑𝑥
𝑑ℎ 𝑑𝑔
= (5.5)
𝑑(𝑥2 ) 𝑑𝑥
Steps
1. Solve for the inner derivative of 𝑔(𝑥) = 𝑥2

𝑑𝑔
= 2𝑥
𝑑𝑥
2. Solve for the outer derivative of ℎ(𝑥) = 𝑥3 , using a placeholder 𝑏 to represent the inner function 𝑥2
𝑑ℎ
= 3𝑏2
𝑑𝑏
3. Swap out the placeholder variable (b) for the inner function (g(x))
3(𝑥2 )2
3𝑥4
4. Return the product of the two derivatives
3𝑥4 · 2𝑥 = 6𝑥5
5.3. Chain rule 49

5.3.3 Multiple functions
In the above example we assumed a composite function containing a single inner function. But the chain rule can also
be applied to higher-order functions like:
𝑓 (𝑥) = 𝐴(𝐵(𝐶(𝑥)))
The chain rule tells us that the derivative of this function equals:
𝑑𝑓 𝑑𝐴 𝑑𝐵 𝑑𝐶
=
𝑑𝑥 𝑑𝐵 𝑑𝐶 𝑑𝑥
We can also write this derivative equation 𝑓 ′ notation:
𝑓 ′ = 𝐴′ (𝐵(𝐶(𝑥)) · 𝐵 ′ (𝐶(𝑥)) · 𝐶 ′ (𝑥)
Steps
Given the function 𝑓 (𝑥) = 𝐴(𝐵(𝐶(𝑥))), lets assume:

𝐴(𝑥) = 𝑠𝑖𝑛(𝑥) (5.6)
2
𝐵(𝑥) = 𝑥
(5.7)
𝐶(𝑥) = 4𝑥
(5.8)
The derivatives of these functions would be:
𝐴′ (𝑥) = 𝑐𝑜𝑠(𝑥) (5.9)
′
𝐵 (𝑥)(5.10)
= 2𝑥
𝐶 ′ (𝑥) =4
(5.11)
We can calculate the derivative of 𝑓 (𝑥) using the following formula:
𝑓 ′ (𝑥) = 𝐴′ ((4𝑥)2 ) · 𝐵 ′ (4𝑥) · 𝐶 ′ (𝑥)
We then input the derivatives and simplify the expression:

𝑓 ′ (𝑥) = 𝑐𝑜𝑠((4𝑥)2 ) · 2(4𝑥) · 4 (5.12)
2
= 𝑐𝑜𝑠(16𝑥 ) · (5.13)
8𝑥 · 4
2
= 𝑐𝑜𝑠(16𝑥(5.14)
)32𝑥
5.4 Gradients
A gradient is a vector that stores the partial derivatives of multivariable functions. It helps us calculate the slope at a
specific point on a curve for functions with multiple independent variables. In order to calculate this more complex
slope, we need to isolate each variable to determine how it impacts the output on its own. To do this we iterate
through each of the variables and calculate the derivative of the function after holding all other variables constant.
Each iteration produces a partial derivative which we store in the gradient.
5.4.1 Partial derivatives
In functions with 2 or more variables, the partial derivative is the derivative of one variable with respect to the others.
If we change 𝑥, but hold all other variables constant, how does 𝑓 (𝑥, 𝑧) change? That’s one partial derivative. The next
variable is 𝑧. If we change 𝑧 but hold 𝑥 constant, how does 𝑓 (𝑥, 𝑧) change? We store partial derivatives in a gradient,
which represents the full derivative of the multivariable function.
5.4.2 Step-by-step
Here are the steps to calculate the gradient for a multivariable function:
1. Given a multivariable function
𝑓 (𝑥, 𝑧) = 2𝑧 3 𝑥2
2. Calculate the derivative with respect to 𝑥
𝑑𝑓
(𝑥, 𝑧)
𝑑𝑥
3. Swap 2𝑧 3 with a constant value 𝑏
𝑓 (𝑥, 𝑧) = 𝑏𝑥2
4. Calculate the derivative with 𝑏 constant
𝑑𝑓 𝑓 (𝑥 + ℎ) − 𝑓 (𝑥)
= lim (5.15)
𝑑𝑥 ℎ→0 ℎ
𝑏(𝑥 + ℎ)2 − 𝑏(𝑥2 )
= lim (5.16)
ℎ→0 ℎ
𝑏((𝑥 + ℎ)(𝑥 + ℎ)) − 𝑏𝑥2
= lim (5.17)
ℎ→0 ℎ
𝑏((𝑥2 + 𝑥ℎ + ℎ𝑥 + ℎ2 )) − 𝑏𝑥2
= lim (5.18)
ℎ→0 ℎ
𝑏𝑥2 + 2𝑏𝑥ℎ + 𝑏ℎ2 − 𝑏𝑥2
= lim (5.19)
ℎ→0 ℎ
2𝑏𝑥ℎ + 𝑏ℎ2
= lim (5.20)
ℎ→0 ℎ
= lim 2𝑏𝑥(5.21)
+ 𝑏ℎ
ℎ→0
(5.22)
As ℎ¯ > 0. . .
2bx + 0
5. Swap 2𝑧 3 back into the equation, to find the derivative with respect to 𝑥.
𝑑𝑓
(𝑥, 𝑧) = 2(2𝑧 3 )𝑥 (5.23)
𝑑𝑥
4𝑧 3 𝑥
=(5.24)
6. Repeat the above steps to calculate the derivative with respect to 𝑧
𝑑𝑓
(𝑥, 𝑧) = 6𝑥2 𝑧 2
𝑑𝑧
7. Store the partial derivatives in a gradient
[︂ 𝑑𝑓 ]︂
4𝑧 3 𝑥
[︂ ]︂
∇𝑓 (𝑥, 𝑧) = 𝑑𝑥 =
𝑑𝑓
𝑑𝑧
6𝑥2 𝑧 2
5.4.3 Directional derivatives
Another important concept is directional derivatives. When calculating the partial derivatives of multivariable func-
tions we use our old technique of analyzing the impact of infinitesimally small increases to each of our independent
variables. By increasing each variable we alter the function output in the direction of the slope.
5.4. Gradients 51
But what if we want to change directions? For example, imagine we’re traveling north through mountainous terrain
on a 3-dimensional plane. The gradient we calculated above tells us we’re traveling north at our current location. But
what if we wanted to travel southwest? How can we determine the steepness of the hills in the southwest direction?
Directional derivatives help us find the slope if we move in a direction different from the one specified by the gradient.
Math
The directional derivative is computed by taking the dot product11 of the gradient of 𝑓 and a unit vector ⃗𝑣 of “tiny
nudges” representing the direction. The unit vector describes the proportions we want to move in each direction. The
output of this calculation is a scalar number representing how much 𝑓 will change if the current input moves with
vector ⃗𝑣 .
Let’s say you have the function 𝑓 (𝑥, 𝑦, 𝑧) and you want to compute its directional derivative along the following
vector2 :
⎡ ⎤
2
⃗𝑣 = ⎣ 3 ⎦
−1
As described above, we take the dot product of the gradient and the directional vector:
⎡ ⎤ ⎡ ⎤
𝑑𝑓
2
⎢ 𝑑𝑥
𝑑𝑓 ⎥
⎣ 𝑑𝑦 ⎦ · ⎣ 3 ⎦
𝑑𝑓 −1
𝑑𝑧
We can rewrite the dot product as:
∇⃗𝑣𝑓 =2 𝑑𝑓 +3 𝑑𝑓 −1 𝑑𝑓
𝑑𝑥 𝑑𝑦 𝑑𝑧
This should make sense because a tiny nudge along ⃗𝑣 can be broken down into two tiny nudges in the x-direction,
three tiny nudges in the y-direction, and a tiny nudge backwards, by 1 in the z-direction.
5.4.4 Useful properties
There are two additional properties of gradients that are especially useful in deep learning. The gradient of a function:
1. Always points in the direction of greatest increase of a function (explained here)
2. Is zero at a local maximum or local minimum
5.5 Integrals
The integral of 𝑓 (𝑥) corresponds to the computation of the area under the graph of 𝑓 (𝑥). The area under 𝑓 (𝑥) between
the points 𝑥 = 𝑎 and 𝑥 = 𝑏 is denoted as follows:
∫︁ 𝑏
𝐴(𝑎, 𝑏) = 𝑓 (𝑥) 𝑑𝑥.
𝑎
11 https://en.wikipedia.org/wiki/Dot_product
2 https://www.khanacademy.org/math/multivariable-calculus/multivariable-derivatives/partial-derivative-and-gradient-articles/a/
directional-derivative-introduction
The area 𝐴(𝑎, 𝑏) is bounded by the function 𝑓 (𝑥) from above, by the 𝑥-axis from below,∫︀ and by two vertical lines at
𝑥 = 𝑎 and 𝑥 = 𝑏. The points 𝑥 = 𝑎 and 𝑥 = 𝑏 are called the limits of integration. The sign comes from the Latin
word summa. The integral is the “sum” of the values of 𝑓 (𝑥) between the two limits of integration.
The integral function 𝐹 (𝑐) corresponds to the area calculation as a function of the upper limit of integration:
∫︁ 𝑐
𝐹 (𝑐) ≡ 𝑓 (𝑥) 𝑑𝑥 .
0
There are two variables and one constant in this formula. The input variable 𝑐 describes the upper limit of integration.
The integration variable 𝑥 performs a sweep from 𝑥 = 0 until 𝑥 = 𝑐. The constant 0 describes the lower limit of
integration. Note that choosing 𝑥 = 0 for the starting point of the integral function was an arbitrary choice.
The integral function 𝐹 (𝑐) contains the “precomputed” information about the area under the graph of 𝑓 (𝑥). The
derivative function 𝑓 ′ (𝑥) tells us the “slope of the graph” property of the function 𝑓 (𝑥) for all values of 𝑥. Similarly,
the integral function 𝐹 (𝑐) tells us the “area under the graph” property of the function 𝑓 (𝑥) for all possible limits of
integration.
The area under 𝑓 (𝑥) between 𝑥 = 𝑎 and 𝑥 = 𝑏 is obtained by calculating the change in the integral function as follows:
∫︁ 𝑏
𝐴(𝑎, 𝑏) = 𝑓 (𝑥) 𝑑𝑥 = 𝐹 (𝑏) − 𝐹 (𝑎).
𝑎
5.5. Integrals 53
5.5.1 Computing integrals
We can approximate the total area under the function 𝑓 (𝑥) between 𝑥 = 𝑎 and 𝑥 = 𝑏 by splitting the region into tiny
vertical strips of width ℎ, then adding up the areas of the rectangular strips. The figure below shows how to compute
the area under 𝑓 (𝑥) = 𝑥2 between 𝑥 = 1 and 𝑥 = 3 by approximating it as four rectangular strips of width ℎ = 0.5.
Usually we want to choose ℎ to be a small number so that the approximation is accurate. Here is some sample code
that performs integration.
def get_integral(func, a, b):

"""Compute the area under `func` between x=a and x=b."""
h = 0.0001 # width of small rectangle
x = a # start at x=a
total = 0
while x <= b: # continue until x=b
total += h*func(x) # area of rect is base*height
x += h
return total
def f(x): return x**2 # some test function f(x)=x^2

computed = get_integral(f, 1, 3)
def actualF(x): return 1.0/3.0*x**3
actual = actualF(3) - actualF(1)
computed, actual # = 8.6662, 8.6666 # pretty close if you ask me...
You can find integral functions using the derivative formulas and some reverse engineering. To find an integral function
of the function 𝑓 (𝑥), we must find a function 𝐹 (𝑥) such that 𝐹 ′ (𝑥) = 𝑓 (𝑥). Suppose you’re given a function 𝑓 (𝑥)
and asked to find its integral function 𝐹 (𝑥):
∫︁
𝐹 (𝑥) = 𝑓 (𝑥) 𝑑𝑥.
This problem is equivalent to finding a function 𝐹 (𝑥) whose derivative is 𝑓 (𝑥):
𝐹 ′ (𝑥) = 𝑓 (𝑥).
For example, suppose you want to find the indefinite integral 𝑥2 𝑑𝑥. We can rephrase this problem as the search for
∫︀
some function 𝐹 (𝑥) such that
𝐹 ′ (𝑥) = 𝑥2 .
Remembering the derivative formulas we saw above, you guess that 𝐹 (𝑥) must contain an 𝑥3 term. Taking the
derivative of a cubic term results in a quadratic term. Therefore, the function you are looking for has the form
𝐹 (𝑥) = 𝑐𝑥3 , for some constant 𝑐. Pick the constant 𝑐 that makes this equation true:
𝐹 ′ (𝑥) = 3𝑐𝑥2 = 𝑥2 .
1
Solving 3𝑐 = 1, we find 𝑐 = 3and so the integral function is
∫︁
1
𝐹 (𝑥) = 𝑥2 𝑑𝑥 = 𝑥3 + 𝐶.
3
𝑑 1 3
[︀ ]︀ 2
You can verify that 𝑑𝑥 3𝑥 + 𝐶 = 𝑥 .
You can also verify Integrals using maths. Here is a set of formulas for your reference
5.5.2 Applications of integration
Integral calculations have widespread applications to more areas of science than are practical to list here. Let’s explore
a few examples related to probabilities.
Computing probabilities
A continuous random variable 𝑋 is described by its probability density function 𝑝(𝑥). A probability density function
𝑝(𝑥) is a positive function for which the total area under the curve is 1:
∫︁ ∞
𝑝(𝑥) ≥ 0, ∀𝑥 and 𝑝(𝑥) 𝑑𝑥 = 1.
−∞
The probability of observing a value of 𝑋 between 𝑎 and 𝑏 is given by the integral

∫︁ 𝑏
Pr(𝑎 ≤ 𝑋 ≤ 𝑏) = 𝑝(𝑥) 𝑑𝑥.
𝑎
Thus, the notion of integration is central to probability theory with continuous random variables.
We also use integration to compute certain characteristic properties of the random variable. The expected value and
the variance are two properties of any random variable 𝑋 that capture important aspects of its behaviour.
Expected value
The expected value of the random variable 𝑋 is computed using the formula
∫︁ ∞
𝜇= 𝑥 𝑝(𝑥).
−∞
The expected value is a single number that tells us what value of 𝑋 we can expect to obtain on average from the
random variable 𝑋. The expected value is also called the average or the mean of the random variable 𝑋.
5.5. Integrals 55
Variance
The variance of the random variable 𝑋 is defined as follows:

∫︁ ∞
𝜎2 = (𝑥 − 𝜇)2 𝑝(𝑥).
−∞
The variance formula computes the expectation of the squared distance of the random variable 𝑋 from its expected
value. The variance 𝜎 2 , also denoted var(𝑋), gives us an indication of how clustered or spread the values of 𝑋 are.
A small variance indicates the outcomes of 𝑋 are tightly clustered near the expected value 𝜇, while a large variance
indicates the outcomes of 𝑋 are widely spread. The square root of the variance is called the standard deviation and is
usually denoted 𝜎.
The expected value 𝜇 and the variance 𝜎 2 are two central concepts in probability theory and statistics because they
allow us to characterize any random variable. The expected value is a measure of the central tendency of the random
variable, while the variance 𝜎 2 measures its dispersion. Readers familiar with concepts from physics can think of the
expected value as the centre of mass of the distribution, and the variance as the moment of inertia of the distribution.
References
CHAPTER 6
Linear Algebra
• Vectors
– Notation
– Vectors in geometry
– Scalar operations
– Elementwise operations
– Dot product
– Hadamard product
– Vector fields
• Matrices
– Dimensions
– Scalar operations
– Elementwise operations
– Hadamard product
– Matrix transpose
– Matrix multiplication
– Test yourself
• Numpy
– Dot product
– Broadcasting
Linear algebra is a mathematical toolbox that offers helpful techniques for manipulating groups of numbers simulta-
57
neously. It provides structures like vectors and matrices (spreadsheets) to hold these numbers and new rules for how
to add, subtract, multiply, and divide them. Here is a brief overview of basic linear algebra concepts taken from my
linear algebra post on Medium.
6.1 Vectors
Vectors are 1-dimensional arrays of numbers or terms. In geometry, vectors store the magnitude and direction of a
potential change to a point. The vector [3, -2] says go right 3 and down 2. A vector with more than one dimension is
called a matrix.
6.1.1 Notation
There are a variety of ways to represent vectors. Here are a few you might come across in your reading.
⎡ ⎤ ⎛ ⎞
1 1 [︀ ]︀
𝑣 = ⎣2⎦ = ⎝2⎠ = 1 2 3
3 3
6.1.2 Vectors in geometry
Vectors typically represent movement from a point. They store both the magnitude and direction of potential changes
to a point. The vector [-2,5] says move left 2 units and up 5 units1 .
A vector can be applied to any point in space. The vector’s direction equals the slope of the hypotenuse created moving
up 5 and left 2. Its magnitude equals the length of the hypotenuse.
1 http://mathinsight.org/vector_introduction
58 Chapter 6. Linear Algebra

6.1.3 Scalar operations
Scalar operations involve a vector and a number. You modify the vector in-place by adding, subtracting, or multiplying
the number from all the values in the vector.
⎡ ⎤ ⎡ ⎤
2 3
⎣2⎦ + 1 = ⎣3⎦
2 3
6.1.4 Elementwise operations
In elementwise operations like addition, subtraction, and division, values that correspond positionally are combined to
produce a new vector. The 1st value in vector A is paired with the 1st value in vector B. The 2nd value is paired with
the 2nd, and so on. This means the vectors must have equal dimensions to complete the operation.*
[︂ ]︂ [︂ ]︂ [︂ ]︂
𝑎1 𝑏 𝑎 + 𝑏1
+ 1 = 1
𝑎2 𝑏2 𝑎2 + 𝑏2
y = np.array([1,2,3])
x = np.array([2,3,4])
y + x = [3, 5, 7]
y - x = [-1, -1, -1]
y / x = [.5, .67, .75]
See below for details on broadcasting in numpy.
6.1.5 Dot product
The dot product of two vectors is a scalar. Dot product of vectors and matrices (matrix multiplication) is one of the
most important operations in deep learning.
[︂ ]︂ [︂ ]︂
𝑎1 𝑏
· 1 = 𝑎1 𝑏1 + 𝑎2 𝑏2
𝑎2 𝑏2
y = np.array([1,2,3])
x = np.array([2,3,4])
np.dot(y,x) = 20
6.1.6 Hadamard product
Hadamard Product is elementwise multiplication and it outputs a vector.

[︂ ]︂ [︂ ]︂ [︂ ]︂
𝑎1 𝑏1 𝑎1 · 𝑏1
⊙ =
𝑎2 𝑏2 𝑎2 · 𝑏2
y = np.array([1,2,3])
x = np.array([2,3,4])
y * x = [2, 6, 12]
6.1. Vectors 59
6.1.7 Vector fields
A vector field shows how far the point (x,y) would hypothetically move if we applied a vector function to it like
addition or multiplication. Given a point in space, a vector field shows the power and direction of our proposed change
at a variety of points in a graph2 .
This vector field is an interesting one since it moves in different directions depending the starting point. The reason is
that the vector behind this field stores terms like 2𝑥 or 𝑥2 instead of scalar values like -2 and 5. For each point on the
graph, we plug the x-coordinate into 2𝑥 or 𝑥2 and draw an arrow from the starting point to the new location. Vector
fields are extremely useful for visualizing machine learning techniques like Gradient Descent.
6.2 Matrices
A matrix is a rectangular grid of numbers or terms (like an Excel spreadsheet) with special rules for addition, subtrac-
tion, and multiplication.
6.2.1 Dimensions
We describe the dimensions of a matrix in terms of rows by columns.

⎡ ⎤
2 4 [︂ ]︂
⎣ 5 −7⎦ 𝑎 2𝑎 8
18 7𝑎 − 4 10
12 5
The first has dimensions (3,2). The second (2,3).
a = np.array([
[1,2,3],
[4,5,6]
2 https://en.wikipedia.org/wiki/Vector_field


])
a.shape == (2,3)
b = np.array([
[1,2,3]
])
b.shape == (1,3)
6.2.2 Scalar operations
Scalar operations with matrices work the same way as they do for vectors. Simply apply the scalar to every element
in the matrix—add, subtract, divide, multiply, etc.
⎡ ⎤ ⎡ ⎤
2 3 3 4
⎣2 3⎦ + 1 = ⎣3 4⎦
2 3 3 4
# Addition
a = np.array(
[[1,2],
[3,4]])
a + 1
[[2,3],
[4,5]]
6.2.3 Elementwise operations
In order to add, subtract, or divide two matrices they must have equal dimensions. We combine corresponding values
in an elementwise fashion to produce a new matrix.
[︂ ]︂ [︂ ]︂ [︂ ]︂
𝑎 𝑏 1 2 𝑎+1 𝑏+2
+ =
𝑐 𝑑 3 4 𝑐+3 𝑑+4
a = np.array([
[1,2],
[3,4]])
b = np.array([
[1,2],
[3,4]])
a + b
[[2, 4],
[6, 8]]
a -- b
[[0, 0],
[0, 0]]
6.2. Matrices 61
6.2.4 Hadamard product
Hadamard product of matrices is an elementwise operation. Values that correspond positionally are multiplied to
produce a new matrix.
[︂ ]︂ [︂ ]︂ [︂ ]︂
𝑎1 𝑎2 𝑏1 𝑏2 𝑎1 · 𝑏1 𝑎2 · 𝑏2
⊙ =
𝑎3 𝑎4 𝑏3 𝑏4 𝑎3 · 𝑏3 𝑎4 · 𝑏4
a = np.array(
[[2,3],
[2,3]])
b = np.array(
[[3,4],
[5,6]])
# Uses python's multiply operator

a * b
[[ 6, 12],
[10, 18]]
In numpy you can take the Hadamard product of a matrix and vector as long as their dimensions meet the requirements
of broadcasting.
[︂ ]︂ [︂ ]︂ [︂ ]︂
𝑎1 𝑏1 𝑏2 𝑎1 · 𝑏1 𝑎1 · 𝑏2
⊙ =
𝑎2 𝑏3 𝑏4 𝑎2 · 𝑏3 𝑎2 · 𝑏4
6.2.5 Matrix transpose
Neural networks frequently process weights and inputs of different sizes where the dimensions do not meet the re-
quirements of matrix multiplication. Matrix transposition (often denoted by a superscript ‘T’ e.g. M^T) provides a
way to “rotate” one of the matrices so that the operation complies with multiplication requirements and can continue.
There are two steps to transpose a matrix:
1. Rotate the matrix right 90°
2. Reverse the order of elements in each row (e.g. [a b c] becomes [c b a])
As an example, transpose matrix M into T:
⎡ ⎤
𝑎 𝑏 [︂ ]︂
𝑎 𝑐 𝑒
⎣𝑐 𝑑⎦ ⇒
𝑏 𝑑 𝑓
𝑒 𝑓
a = np.array([
[1, 2],
[3, 4]])
a.T
[[1, 3],
[2, 4]]
6.2.6 Matrix multiplication
Matrix multiplication specifies a set of rules for multiplying matrices together to produce a new matrix.

Rules
Not all matrices are eligible for multiplication. In addition, there is a requirement on the dimensions of the resulting
matrix output. Source.
1. The number of columns of the 1st matrix must equal the number of rows of the 2nd
2. The product of an M x N matrix and an N x K matrix is an M x K matrix. The new matrix takes the rows of the
1st and columns of the 2nd
Steps
Matrix multiplication relies on dot product to multiply various combinations of rows and columns. In the image
below, taken from Khan Academy’s excellent linear algebra course, each entry in Matrix C is the dot product of a row
in matrix A and a column in matrix B3 .
The operation a1 · b1 means we take the dot product of the 1st row in matrix A (1, 7) and the 1st column in matrix B
(3, 5).
[︂ ]︂ [︂ ]︂
1 3
𝑎1 · 𝑏1 = · = (1 · 3) + (7 · 5) = 38
7 5
Here’s another way to look at it:

⎡ ⎤ ⎡ ⎤
𝑎 𝑏 [︂ ]︂ 1𝑎 + 3𝑏 2𝑎 + 4𝑏
1 2
⎣𝑐 𝑑⎦ · = ⎣ 1𝑐 + 3𝑑 2𝑐 + 4𝑑 ⎦
3 4
𝑒 𝑓 1𝑒 + 3𝑓 2𝑒 + 4𝑓
6.2.7 Test yourself
1. What are the dimensions of the matrix product?

[︂ ]︂ [︂ ]︂
1 2 1 2 3
· =2x3
5 6 5 6 7
2. What are the dimensions of the matrix product?
⎡ ⎤
⎡ ⎤ 1 2
1 2 3 4 ⎢5 6⎥
⎣5 6 7 8⎦·⎢⎣3
⎥=3x2
0⎦
9 10 11 12
2 1
3 https://www.khanacademy.org/math/precalculus/precalc-matrices/properties-of-matrix-multiplication/a/properties-of-matrix-multiplication
6.2. Matrices 63
3. What is the matrix product?

[︂ ]︂ [︂ ]︂ [︂ ]︂
2 3 5 4 19 23
· =
1 4 3 5 17 24
4. What is the matrix product?}
[︂ ]︂ [︂ ]︂
3 [︀ ]︀ 3 6 9
· 1 2 3 =
5 5 10 15
5. What is the matrix product?
⎡ ⎤
[︀ ]︀ 4 [︀ ]︀
1 2 3 · ⎣5⎦ = 32
6
6.3 Numpy
6.3.1 Dot product
Numpy uses the function np.dot(A,B) for both vector and matrix multiplication. It has some other interesting features
and gotchas so I encourage you to read the documentation here before use.
a = np.array([
[1, 2]
])
a.shape == (1,2)
b = np.array([
[3, 4],
[5, 6]
])
b.shape == (2,2)
# Multiply
mm = np.dot(a,b)
mm == [13, 16]
mm.shape == (1,2)
6.3.2 Broadcasting
In numpy the dimension requirements for elementwise operations are relaxed via a mechanism called broadcasting.
Two matrices are compatible if the corresponding dimensions in each matrix (rows vs rows, columns vs columns)
meet the following requirements:
1. The dimensions are equal, or
2. One dimension is of size 1
a = np.array([
[1],
[2]
])
b = np.array([
[3,4],
[5,6]
])


c = np.array([
[1,2]
])
# Same no. of rows

# Different no. of columns
# but a has one column so this works
a * b
[[ 3, 4],
[10, 12]]
# Same no. of columns

# Different no. of rows
# but c has one row so this works
b * c
[[ 3, 8],
[5, 12]]
# Different no. of columns

# Different no. of rows
# but both a and c meet the
# size 1 requirement rule
a + c
[[2, 3],
[3, 4]]
Tutorials
• Khan Academy Linear Algebra

• Deep Learning Book Math
• Andrew Ng Course Notes
• Linear Algebra Better Explained
• Understanding Matrices Intuitively
• Intro To Linear Algebra
• Immersive Math
References
6.3. Numpy 65

CHAPTER 7
Probability
• Links
• Screenshots
• License
Basic concepts in probability for machine learning.

This cheatsheet is a 10-page reference in probability that covers a semester’s worth of introductory probability.
The cheatsheet is based off of Harvard’s introductory probability course, Stat 110. It is co-authored by former Stat 110
Teaching Fellow William Chen and Stat 110 Professor Joe Blitzstein.
7.1 Links
• [Probability Cheatsheet PDF](http://www.wzchen.com/probability-cheatsheet/)
7.2 Screenshots
![First Page](http://i.imgur.com/Oa73huL.jpg) ![Second Page](http://i.imgur.com/dyvW2rB.jpg)
7.3 License
This work is licensed under a [Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International

License.][by-nc-sa].
[![Creative Commons License][by-nc-sa-img]][by-nc-sa]
67
References
68 Chapter 7. Probability
CHAPTER 8
Statistics
Basic concepts in statistics for machine learning.
References
69
70 Chapter 8. Statistics
CHAPTER 9
Notation
Commonly used math symbols in machine learning texts.
• Algebra
• Calculus
• Linear algebra
• Probability
• Set theory
• Statistics
Note: Use the table generator to quickly add new symbols. Import current tables into tablesgenerator from
figures/*.tgn. Export and save your changes. Also see helpful multiline editing in Sublime.
9.1 Algebra
Symbol Name Description Example

(𝑓 𝑔) composite function a nested function (f g)(x) = f(g(x))
delta change / difference x = x_1 - x_0
𝑒∑︀ Euler’s number e = 2.718281828 s = frac{1}{1+e^{-z}}
∏︀ summation sum of all values x_i = x_1 + x_2 + x_3
capital pi product of all values x_i = x_1x_2x_3
𝜖 epsilon tiny number near 0 lr = 1e-4
71
9.2 Calculus

𝑥′ derivative first derivative (x^2)’ = 2x
𝑥′′ second derivative second derivative (x^2)” = 2
lim limit function value as x approaches 0
nabla gradient f(a,b,c)
9.3 Linear algebra

[] brackets matrix or vector 𝑀 = [135]
· dot dot product (𝑍 = 𝑋 · 𝑊
⊙ hadamard hadamard product 𝐴=𝐵⊙𝐶
𝑋𝑇 transpose matrix transpose 𝑊𝑇 · 𝑋
⃗𝑥 vector vector 𝑣 = [123]
𝑋 matrix capitalized variables are matrices 𝑋, 𝑊, 𝐵
𝑥
ˆ unit vector vector of magnitude 1 𝑥
ˆ = [0.20.50.3]
9.4 Probability

𝑃 (𝐴) probability probability of event A P(x=1) = 0.5
9.5 Set theory

set list of distinct elements S = {1, 5, 7, 9}
9.6 Statistics

population mean mean of population values
𝑥
¯ sample mean mean of subset of population
2
population variance variance of population value
𝑠2 sample variance variance of subset of population
𝑋 standard deviation population standard deviation
𝑠 sample std dev standard deviation of sample
𝑋 correlation correlation of variables X and Y
𝑥
˜ median median value of variable x
72 Chapter 9. Notation
References
9.6. Statistics 73
74 Chapter 9. Notation
CHAPTER 10
Concepts
• Neural Network
• Neuron
• Synapse
• Weights
• Bias
• Layers
• Weighted Input
• Activation Functions
• Loss Functions
• Optimization Algorithms
• Gradient Accumulation
10.1 Neural Network
Neural networks are a class of machine learning algorithms used to model complex patterns in datasets using multiple
hidden layers and non-linear activation functions. A neural network takes an input, passes it through multiple layers of
hidden neurons (mini-functions with unique coefficients that must be learned), and outputs a prediction representing
the combined input of all the neurons.
75
Neural networks are trained iteratively using optimization techniques like gradient descent. After each cycle of train-
ing, an error metric is calculated based on the difference between prediction and target. The derivatives of this error
metric are calculated and propagated back through the network using a technique called backpropagation. Each neu-
ron’s coefficients (weights) are then adjusted relative to how much they contributed to the total error. This process is
repeated iteratively until the network error drops below an acceptable threshold.
10.2 Neuron
A neuron takes a group of weighted inputs, applies an activation function, and returns an output.
Inputs to a neuron can either be features from a training set or outputs from a previous layer’s neurons. Weights are
applied to the inputs as they travel along synapses to reach the neuron. The neuron then applies an activation function
to the “sum of weighted inputs” from each incoming synapse and passes the result on to all the neurons in the next
layer.
76 Chapter 10. Concepts

10.3 Synapse
Synapses are like roads in a neural network. They connect inputs to neurons, neurons to neurons, and neurons to
outputs. In order to get from one neuron to another, you have to travel along the synapse paying the “toll” (weight)
along the way. Each connection between two neurons has a unique synapse with a unique weight attached to it. When
we talk about updating weights in a network, we’re really talking about adjusting the weights on these synapses.
10.4 Weights
Weights are values that control the strength of the connection between two neurons. That is, inputs are typically
multiplied by weights, and that defines how much influence the input will have on the output. In other words: when
the inputs are transmitted between neurons, the weights are applied to the inputs along with an additional value (the
bias)
10.5 Bias
Bias terms are additional constants attached to neurons and added to the weighted input before the activation function
is applied. Bias terms help models represent patterns that do not necessarily pass through the origin. For example,
if all your features were 0, would your output also be zero? Is it possible there is some base value upon which your
features have an effect? Bias terms typically accompany weights and must also be learned by your model.
10.6 Layers
Input Layer
Holds the data your model will train on. Each neuron in the input layer represents a unique attribute in your dataset
(e.g. height, hair color, etc.).
Hidden Layer
Sits between the input and output layers and applies an activation function before passing on the results. There are often
multiple hidden layers in a network. In traditional networks, hidden layers are typically fully-connected layers—each
neuron receives input from all the previous layer’s neurons and sends its output to every neuron in the next layer. This
contrasts with how convolutional layers work where the neurons send their output to only some of the neurons in the
next layer.
10.3. Synapse 77
Output Layer
The final layer in a network. It receives input from the previous hidden layer, optionally applies an activation function,
and returns an output representing your model’s prediction.
10.7 Weighted Input
A neuron’s input equals the sum of weighted outputs from all neurons in the previous layer. Each input is multiplied
by the weight associated with the synapse connecting the input to the current neuron. If there are 3 inputs or neurons
in the previous layer, each neuron in the current layer will have 3 distinct weights—one for each each synapse.
Single Input
𝑍 = 𝐼𝑛𝑝𝑢𝑡 · 𝑊 𝑒𝑖𝑔ℎ𝑡
= 𝑋𝑊
Multiple Inputs
𝑛
∑︁
𝑍= 𝑥𝑖 𝑤𝑖
𝑖=1
= 𝑥1 𝑤1 + 𝑥2 𝑤2 + 𝑥3 𝑤3
Notice, it’s exactly the same equation we use with linear regression! In fact, a neural network with a single neuron
is the same as linear regression! The only difference is the neural network post-processes the weighted input with an
activation function.
10.8 Activation Functions
Activation functions live inside neural network layers and modify the data they receive before passing it to the next
layer. Activation functions give neural networks their power—allowing them to model complex non-linear relation-
ships. By modifying inputs with non-linear functions neural networks can model highly complex relationships between
features. Popular activation functions include relu and sigmoid.
Activation functions typically have the following properties:
• Non-linear - In linear regression we’re limited to a prediction equation that looks like a straight line. This is
nice for simple datasets with a one-to-one relationship between inputs and outputs, but what if the patterns in
our dataset were non-linear? (e.g. 𝑥2 , sin, log). To model these relationships we need a non-linear prediction
equation.1 Activation functions provide this non-linearity.
• Continuously differentiable—To improve our model with gradient descent, we need our output to have a nice
slope so we can compute error derivatives with respect to weights. If our neuron instead outputted 0 or 1
(perceptron), we wouldn’t know in which direction to update our weights to reduce our error.
• Fixed Range—Activation functions typically squash the input data into a narrow range that makes training the
model more stable and efficient.
10.9 Loss Functions
A loss function, or cost function, is a wrapper around our model’s predict function that tells us “how good” the model
is at making predictions for a given set of parameters. The loss function has its own curve and its own derivatives. The

slope of this curve tells us how to change our parameters to make the model more accurate! We use the model to make
predictions. We use the cost function to update our parameters. Our cost function can take a variety of forms as there
are many different cost functions available. Popular loss functions include: MSE (L2) and Cross-entropy Loss.
10.10 Optimization Algorithms
Be the first to contribute!
10.11 Gradient Accumulation
Gradient accumulation is a mechanism to split the batch of samples—used for training a neural network—into several
mini-batches of samples that will be run sequentially.
This is used to enable using large batch sizes that require more GPU memory than available. Gradient accumulation
helps in doing so by using mini-batches that require an amount of GPU memory that can be satisfied.
Gradient accumulation means running all mini-batches sequentially (generally on the same GPU) while accumulating
their calculated gradients and not updating the model variables - the weights and biases of the model. The model
variables must not be updated during the accumulation in order to ensure all mini-batches use the same model variable
values to calculate their gradients. Only after accumulating the gradients of all those mini-batches will we generate
and apply the updates for the model variables.
This results in the same updates for the model parameters as if we were to use the global batch.
10.10. Optimization Algorithms 79

More details, a technical and algorithmical deep-dive, how-to tutorials, and examples can be found at [2].
References

CHAPTER 11
Forwardpropagation
• Simple Network
– Steps
– Code
• Larger Network
– Architecture
– Weight Initialization
– Bias Terms
– Working with Matrices
– Dynamic Resizing
– Refactoring Our Code
– Final Result
81
11.1 Simple Network
Forward propagation is how neural networks make predictions. Input data is “forward propagated” through the net-
work layer by layer to the final layer which outputs a prediction. For the toy neural network above, a single pass of
forward propagation translates mathematically to:
𝑃 𝑟𝑒𝑑𝑖𝑐𝑡𝑖𝑜𝑛 = 𝐴( 𝐴( 𝑋𝑊ℎ )𝑊𝑜 )
Where 𝐴 is an activation function like ReLU, 𝑋 is the input and 𝑊ℎ and 𝑊𝑜 are weights.
11.1.1 Steps
1. Calculate the weighted input to the hidden layer by multiplying 𝑋 by the hidden weight 𝑊ℎ
2. Apply the activation function and pass the result to the final layer
3. Repeat step 2 except this time 𝑋 is replaced by the hidden layer’s output, 𝐻
11.1.2 Code
Let’s write a method feed_forward() to propagate input data through our simple network of 1 hidden layer. The output
of this method represents our model’s prediction.
def relu(z):
return max(0,z)
def feed_forward(x, Wh, Wo):

# Hidden layer
Zh = x * Wh
H = relu(Zh)
# Output layer
Zo = H * Wo
output = relu(Zo)
return output
x is the input to the network, Zo and Zh are the weighted inputs and Wo and Wh are the weights.
11.2 Larger Network
The simple network above is helpful for learning purposes, but in reality neural networks are much larger and more
complex. Modern neural networks have many more hidden layers, more neurons per layer, more variables per input,
82 Chapter 11. Forwardpropagation

more inputs per training set, and more output variables to predict. Here is a slightly larger network that will introduce
us to matrices and the matrix operations used to train arbitrarily large neural networks.
11.2.1 Architecture
To accomodate arbitrarily large inputs or outputs, we need to make our code more extensible by adding a few parame-
ters to our network’s __init__ method: inputLayerSize, hiddenLayerSize, outputLayerSize. We’ll still limit ourselves
to using one hidden layer, but now we can create layers of different sizes to respond to the different inputs or outputs.
INPUT_LAYER_SIZE = 1
HIDDEN_LAYER_SIZE = 2
OUTPUT_LAYER_SIZE = 2
11.2.2 Weight Initialization
Unlike last time where Wh and Wo were scalar numbers, our new weight variables will be numpy arrays. Each array
will hold all the weights for its own layer—one weight for each synapse. Below we initialize each array with the
numpy’s np.random.randn(rows, cols) method, which returns a matrix of random numbers drawn from a
normal distribution with mean 0 and variance 1.
def init_weights():
Wh = np.random.randn(INPUT_LAYER_SIZE, HIDDEN_LAYER_SIZE) * \
np.sqrt(2.0/INPUT_LAYER_SIZE)
Wo = np.random.randn(HIDDEN_LAYER_SIZE, OUTPUT_LAYER_SIZE) * \
np.sqrt(2.0/HIDDEN_LAYER_SIZE)
Here’s an example calling random.randn():
arr = np.random.randn(1, 2)
print(arr)
>> [[-0.36094661 -1.30447338]]
11.2. Larger Network 83

print(arr.shape)
>> (1,2)
As you’ll soon see, there are strict requirements on the dimensions of these weight matrices. The number of rows must
equal the number of neurons in the previous layer. The number of columns must match the number of neurons in the
next layer.
A good explanation of random weight initalization can be found in the Stanford CS231 course notes1 chapter on neural
networks.
11.2.3 Bias Terms
Bias terms allow us to shift our neuron’s activation outputs left and right. This helps us model datasets that do not
necessarily pass through the origin.
Using the numpy method np.full() below, we create two 1-dimensional bias arrays filled with the default value
0.2. The first argument to np.full is a tuple of array dimensions. The second is the default value for cells in the
array.
def init_bias():
Bh = np.full((1, HIDDEN_LAYER_SIZE), 0.1)
Bo = np.full((1, OUTPUT_LAYER_SIZE), 0.1)
return Bh, Bo
11.2.4 Working with Matrices
To take advantage of fast linear algebra techniques and GPUs, we need to store our inputs, weights, and biases in
matrices. Here is our neural network diagram again with its underlying matrix representation.
1 http://cs231n.github.io/neural-networks-2/#init

What’s happening here? To better understand, let’s walk through each of the matrices in the diagram with an emphasis
on their dimensions and why the dimensions are what they are. The matrix dimensions above flow naturally from the
architecture of our network and the number of samples in our training set.

Matrix dimensions
Var Name Di- Explanation

men-
sions
X Input (3, Includes 3 rows of training data, and each row has 1 attribute (height, price, etc.)
1)
Wh Hid- (1, These dimensions are based on number of rows equals the number of attributes for the obser-
den 2) vations in our training set. The number columns equals the number of neurons in the hidden
weights layer. The dimensions of the weights matrix between two layers is determined by the sizes of
the two layers it connects. There is one weight for every input-to-neuron connection between
the layers.
Bh Hid- (1, Each neuron in the hidden layer has is own bias constant. This bias matrix is added to the
den 2) weighted input matrix before the hidden layer applies ReLU.
bias
Zh Hid- (1, Computed by taking the dot product of X and Wh. The dimensions (1,2) are required by the
den 2) rules of matrix multiplication. Zh takes the rows of in the inputs matrix and the columns of
weighted weights matrix. We then add the hidden layer bias matrix Bh.
input
H Hid- (3, Computed by applying the Relu function to Zh. The dimensions are (3,2)—the number of
den 2) rows matches the number of training samples and the number of columns equals the number
acti- of neurons. Each column holds all the activations for a specific neuron.
va-
tions
Wo Out- (2, The number of rows matches the number of hidden layer neurons and the number of columns
put 2) equals the number of output layer neurons. There is one weight for every hidden-neuron-to-
weights output-neuron connection between the layers.
Bo Out- (1, There is one column for every neuron in the output layer.
put 2)
bias
Zo Out- (3, Computed by taking the dot product of H and Wo and then adding the output layer bias Bo.
put 2) The dimensions are (3,2) representing the rows of in the hidden layer matrix and the columns
weighted of output layer weights matrix.
input
O Out- (3, Each row represents a prediction for a single observation in our training set. Each column is
put 2) a unique attribute we want to predict. Examples of two-column output predictions could be
acti- a company’s sales and units sold, or a person’s height and weight.
va-
tions
11.2.5 Dynamic Resizing
Before we continue I want to point out how the matrix dimensions change with changes to the network architecture or
size of the training set. For example, let’s build a network with 2 input neurons, 3 hidden neurons, 2 output neurons,
and 4 observations in our training set.

Now let’s use same number of layers and neurons but reduce the number of observations in our dataset to 1 instance:
As you can see, the number of columns in all matrices remains the same. The only thing that changes is the number
of rows the layer matrices, which fluctuate with the size of the training set. The dimensions of the weight matrices
remain unchanged. This shows us we can use the same network, the same lines of code, to process any number of
observations.
11.2.6 Refactoring Our Code
Here is our new feed forward code which accepts matrices instead of scalar inputs.
def feed_forward(X):
'''
X - input matrix
Zh - hidden layer weighted input
Zo - output layer weighted input
H - hidden layer activation


y - output layer
yHat - output layer predictions
'''
# Hidden layer
Zh = np.dot(X, Wh) + Bh
H = relu(Zh)
# Output layer
Zo = np.dot(H, Wo) + Bo
yHat = relu(Zo)
return yHat
Weighted input
The first change is to update our weighted input calculation to handle matrices. Using dot product, we multiply the
input matrix by the weights connecting them to the neurons in the next layer. Next we add the bias vector using matrix
addition.
The first column in Bh is added to all the rows in the first column of resulting dot product of X and Wh. The second
value in Bh is added to all the elements in the second column. The result is a new matrix, Zh which has a column for
every neuron in the hidden layer and a row for every observation in our dataset. Given all the layers in our network are
fully-connected, there is one weight for every neuron-to-neuron connection between the layers.
The same process is repeated for the output layer, except the input is now the hidden layer activation H and the weights
Wo.
ReLU activation
The second change is to refactor ReLU to use elementwise multiplication on matrices. It’s only a small change, but its
necessary if we want to work with matrices. np.maximum() is actually extensible and can handle both scalar and
array inputs.

def relu(Z):
return np.maximum(0, Z)
In the hidden layer activation step, we apply the ReLU activation function np.maximum(0,Z) to every cell in the
new matrix. The result is a matrix where all negative values have been replaced by 0. The same process is repeated
for the output layer, except the input is Zo.
11.2.7 Final Result
Putting it all together we have the following code for forward propagation with matrices.
INPUT_LAYER_SIZE = 1
HIDDEN_LAYER_SIZE = 2
OUTPUT_LAYER_SIZE = 2
def init_weights():
Wh = np.random.randn(INPUT_LAYER_SIZE, HIDDEN_LAYER_SIZE) * \
np.sqrt(2.0/INPUT_LAYER_SIZE)
Wo = np.random.randn(HIDDEN_LAYER_SIZE, OUTPUT_LAYER_SIZE) * \
np.sqrt(2.0/HIDDEN_LAYER_SIZE)
def init_bias():
Bh = np.full((1, HIDDEN_LAYER_SIZE), 0.1)
Bo = np.full((1, OUTPUT_LAYER_SIZE), 0.1)
return Bh, Bo
def relu(Z):
return np.maximum(0, Z)
def relu_prime(Z):
'''
Z - weighted input matrix
Returns gradient of Z where all

negative values are set to 0 and
all positive values set to 1
'''
Z[Z < 0] = 0
Z[Z > 0] = 1
return Z
def cost(yHat, y):

cost = np.sum((yHat - y)**2) / 2.0
return cost
def cost_prime(yHat, y):

return yHat - y
def feed_forward(X):
'''
X - input matrix
Zh - hidden layer weighted input
Zo - output layer weighted input
H - hidden layer activation
y - output layer


yHat - output layer predictions
'''
# Hidden layer
H = relu(Zh)
# Output layer
Zo = np.dot(H, Wo) + Bo
yHat = relu(Zo)
References

CHAPTER 12
Backpropagation
• Chain rule refresher

• Applying the chain rule
• Saving work with memoization
• Code example
The goals of backpropagation are straightforward: adjust each weight in the network in proportion to how much it
contributes to overall error. If we iteratively reduce each weight’s error, eventually we’ll have a series of weights that
produce good predictions.
12.1 Chain rule refresher
As seen above, foward propagation can be viewed as a long series of nested equations. If you think of feed forward
this way, then backpropagation is merely an application of Chain rule to find the Derivatives of cost with respect to
any variable in the nested equation. Given a forward propagation function:
𝑓 (𝑥) = 𝐴(𝐵(𝐶(𝑥)))
A, B, and C are activation functions at different layers. Using the chain rule we easily calculate the derivative of 𝑓 (𝑥)
with respect to 𝑥:
𝑓 ′ (𝑥) = 𝑓 ′ (𝐴) · 𝐴′ (𝐵) · 𝐵 ′ (𝐶) · 𝐶 ′ (𝑥)
How about the derivative with respect to B? To find the derivative with respect to B you can pretend 𝐵(𝐶(𝑥)) is a
constant, replace it with a placeholder variable B, and proceed to find the derivative normally with respect to B.
𝑓 ′ (𝐵) = 𝑓 ′ (𝐴) · 𝐴′ (𝐵)
This simple technique extends to any variable within a function and allows us to precisely pinpoint the exact impact
each variable has on the total output.
91
12.2 Applying the chain rule
Let’s use the chain rule to calculate the derivative of cost with respect to any weight in the network. The chain rule
will help us identify how much each weight contributes to our overall error and the direction to update each weight to
reduce our error. Here are the equations we need to make a prediction and calculate total error, or cost:
Given a network consisting of a single neuron, total cost could be calculated as:
𝐶𝑜𝑠𝑡 = 𝐶(𝑅(𝑍(𝑋𝑊 )))
Using the chain rule we can easily find the derivative of Cost with respect to weight W.
𝐶 ′ (𝑊 ) = 𝐶 ′ (𝑅) · 𝑅′ (𝑍) · 𝑍 ′ (𝑊 )
𝑦 − 𝑦) · 𝑅′ (𝑍) · 𝑋
= (ˆ
Now that we have an equation to calculate the derivative of cost with respect to any weight, let’s go back to our toy
neural network example above
What is the derivative of cost with respect to 𝑊𝑜 ?

𝐶 ′ (𝑊𝑂 ) = 𝐶 ′ (ˆ
𝑦 ) · 𝑦ˆ′ (𝑍𝑂 ) · 𝑍𝑂
′
(𝑊𝑂 )
𝑦 − 𝑦) · 𝑅′ (𝑍𝑂 ) · 𝐻
= (ˆ
And how about with respect to 𝑊ℎ ? To find out we just keep going further back in our function applying the chain
rule recursively until we get to the function that has the Wh term.
𝐶 ′ (𝑊ℎ ) = 𝐶 ′ (ˆ
𝑦 ) · 𝑂′ (𝑍𝑜 ) · 𝑍𝑜′ (𝐻) · 𝐻 ′ (𝑍ℎ ) · 𝑍ℎ′ (𝑊ℎ )
𝑦 − 𝑦) · 𝑅′ (𝑍𝑜 ) · 𝑊𝑜 · 𝑅′ (𝑍ℎ ) · 𝑋
= (ˆ
92 Chapter 12. Backpropagation

And just for fun, what if our network had 10 hidden layers. What is the derivative of cost for the first weight 𝑤1 ?
𝑑𝐶 𝑑ˆ
𝑦 𝑑𝑍11
𝐶 ′ (𝑤1 ) =
· · ·
𝑑ˆ
𝑦 𝑑𝑍11 𝑑𝐻10
𝑑𝐻10 𝑑𝑍10 𝑑𝐻9 𝑑𝑍9 𝑑𝐻8 𝑑𝑍8 𝑑𝐻7
· · · · · · ·
𝑑𝑍10 𝑑𝐻9 𝑑𝑍9 𝑑𝐻8 𝑑𝑍8 𝑑𝐻7 𝑑𝑍7
𝑑𝑍7 𝑑𝐻6 𝑑𝑍6 𝑑𝐻5 𝑑𝑍5 𝑑𝐻4 𝑑𝑍4
· · · · · · ·
𝑑𝐻6 𝑑𝑍6 𝑑𝐻5 𝑑𝑍5 𝑑𝐻4 𝑑𝑍4 𝑑𝐻3
𝑑𝐻3 𝑑𝑍3 𝑑𝐻2 𝑑𝑍2 𝑑𝐻1 𝑑𝑍1
· · · · ·
𝑑𝑍3 𝑑𝐻2 𝑑𝑍2 𝑑𝐻1 𝑑𝑍1 𝑑𝑊1
See the pattern? The number of calculations required to compute cost derivatives increases as our network grows
deeper. Notice also the redundancy in our derivative calculations. Each layer’s cost derivative appends two new terms
to the terms that have already been calculated by the layers above it. What if there was a way to save our work
somehow and avoid these duplicate calculations?
12.3 Saving work with memoization
Memoization is a computer science term which simply means: don’t recompute the same thing over and over. In
memoization we store previously computed results to avoid recalculating the same function. It’s handy for speeding
up recursive functions of which backpropagation is one. Notice the pattern in the derivative equations below.
Each of these layers is recomputing the same derivatives! Instead of writing out long derivative equations for every
weight, we can use memoization to save our work as we backprop error through the network. To do this, we define 3
equations (below), which together encapsulate all the calculations needed for backpropagation. The math is the same,
but the equations provide a nice shorthand we can use to track which calculations we’ve already performed and save
our work as we move backwards through the network.
12.3. Saving work with memoization 93

We first calculate the output layer error and pass the result to the hidden layer before it. After calculating the hidden
layer error, we pass its error value back to the previous hidden layer before it. And so on and so forth. As we move
back through the network we apply the 3rd formula at every layer to calculate the derivative of cost with respect that
layer’s weights. This resulting derivative tells us in which direction to adjust our weights to reduce overall cost.
Note: The term layer error refers to the derivative of cost with respect to a layer’s input. It answers the question: how
does the cost function output change when the input to that layer changes?
Output layer error
To calculate output layer error we need to find the derivative of cost with respect to the output layer input, 𝑍𝑜 . It
answers the question—how are the final layer’s weights impacting overall error in the network? The derivative is then:
𝐶 ′ (𝑍𝑜 ) = (ˆ
𝑦 − 𝑦) · 𝑅′ (𝑍𝑜 )
𝑦 −𝑦)*𝑅′ (𝑍𝑜) sequence with the term 𝐸𝑜 . So our formula

To simplify notation, ml practitioners typically replace the (ˆ
for output layer error equals:
𝑦 − 𝑦) · 𝑅′ (𝑍𝑜 )
𝐸𝑜 = (ˆ
Hidden layer error
To calculate hidden layer error we need to find the derivative of cost with respect to the hidden layer input, Zh.
𝐶 ′ (𝑍ℎ ) = (ˆ
𝑦 − 𝑦) · 𝑅′ (𝑍𝑜 ) · 𝑊𝑜 · 𝑅′ (𝑍ℎ )
Next we can swap in the 𝐸𝑜 term above to avoid duplication and create a new simplified equation for Hidden layer
error:
𝐸ℎ = 𝐸𝑜 · 𝑊𝑜 · 𝑅′ (𝑍ℎ )
This formula is at the core of backpropagation. We calculate the current layer’s error, and pass the weighted error back
to the previous layer, continuing the process until we arrive at our first hidden layer. Along the way we update the
weights using the derivative of cost with respect to each weight.
Derivative of cost with respect to any weight
Let’s return to our formula for the derivative of cost with respect to the output layer weight 𝑊𝑜 .
𝐶 ′ (𝑊𝑂 ) = (ˆ
𝑦 − 𝑦) · 𝑅′ (𝑍𝑂 ) · 𝐻
We know we can replace the first part with our equation for output layer error 𝐸𝑜 . H represents the hidden layer
activation.
𝐶 ′ (𝑊𝑜 ) = 𝐸𝑜 · 𝐻
So to find the derivative of cost with respect to any weight in our network, we simply multiply the corresponding
layer’s error times its input (the previous layer’s output).
𝐶 ′ (𝑤) = 𝐶𝑢𝑟𝑟𝑒𝑛𝑡𝐿𝑎𝑦𝑒𝑟𝐸𝑟𝑟𝑜𝑟 · 𝐶𝑢𝑟𝑟𝑒𝑛𝑡𝐿𝑎𝑦𝑒𝑟𝐼𝑛𝑝𝑢𝑡
Note: Input refers to the activation from the previous layer, not the weighted input, Z.

Summary
Here are the final 3 equations that together form the foundation of backpropagation.
Here is the process visualized using our toy neural network example above.
12.4 Code example
def relu_prime(z):
if z > 0:
return 1
return 0
def cost(yHat, y):

return 0.5 * (yHat - y)**2
def cost_prime(yHat, y):

return yHat - y
def backprop(x, y, Wh, Wo, lr):

yHat = feed_forward(x, Wh, Wo)
# Layer Error
Eo = (yHat - y) * relu_prime(Zo)
Eh = Eo * Wo * relu_prime(Zh)
# Cost derivative for weights

dWo = Eo * H
dWh = Eh * x
# Update weights
12.4. Code example 95


Wh -= lr * dWh
Wo -= lr * dWo
References

CHAPTER 13
Activation Functions
• Linear
• ELU
• ReLU
• LeakyReLU
• Sigmoid
• Tanh
• Softmax
13.1 Linear
A straight line function where activation is proportional to input ( which is the weighted sum from neuron ).
97
Function Derivative
𝑅′ (𝑧, 𝑚) = 𝑚
{︀ }︀ {︀ }︀
𝑅(𝑧, 𝑚) = 𝑧 * 𝑚
def linear(z,m): def linear_prime(z,m):

return m*z return m
Pros
• It gives a range of activations, so it is not binary activation.

• We can definitely connect a few neurons together and if more than 1 fires, we could take the max ( or softmax)
and decide based on that.
Cons
• For this function, derivative is a constant. That means, the gradient has no relationship with X.
• It is a constant gradient and the descent is going to be on constant gradient.
• If there is an error in prediction, the changes made by back propagation is constant and not depending on the
change in input delta(x) !
13.2 ELU
Exponential Linear Unit or its widely known name ELU is a function that tend to converge cost to zero faster and
produce more accurate results. Different to other activation functions, ELU has a extra alpha constant which should
be positive number.
ELU is very similiar to RELU except negative inputs. They are both in identity function form for non-negative inputs.
On the other hand, ELU becomes smooth slowly until its output equal to -𝛼 whereas RELU sharply smoothes.
98 Chapter 13. Activation Functions

Function Derivative
z>0
.( e𝑧 ˘1) z <= 0
{︂ }︂
′ 1 𝑧>0
𝑅 (𝑧) =
.𝑒𝑧 𝑧<0
def elu(z,alpha):
return z if z >= 0 else alpha*(e^z -1)
def elu_prime(z,alpha):
return 1 if z > 0 else alpha*np.exp(z)
Pros
• ELU becomes smooth slowly until its output equal to -𝛼 whereas RELU sharply smoothes.
• ELU is a strong alternative to ReLU.
• Unlike to ReLU, ELU can produce negative outputs.
Cons
• For x > 0, it can blow up the activation with the output range of [0, inf].
13.3 ReLU
A recent invention which stands for Rectified Linear Units. The formula is deceptively simple: 𝑚𝑎𝑥(0, 𝑧). Despite its
name and appearance, it’s not linear and provides the same benefits as Sigmoid but with better performance.
13.3. ReLU 99
Function Derivative
z>0
0 z <= 0
{︂ }︂
′ 1 𝑧>0
𝑅 (𝑧) =
0 𝑧<0
def relu(z):
return max(0, z)
def relu_prime(z):
return 1 if z > 0 else 0
Pros
• It avoids and rectifies vanishing gradient problem.

• ReLu is less computationally expensive than tanh and sigmoid because it involves simpler mathematical opera-
tions.
Cons
• One of its limitation is that it should only be used within Hidden layers of a Neural Network Model.
• Some gradients can be fragile during training and can die. It can cause a weight update which will makes it
never activate on any data point again. Simply saying that ReLu could result in Dead Neurons.
• In another words, For activations in the region (x<0) of ReLu, gradient will be 0 because of which the weights
will not get adjusted during descent. That means, those neurons which go into that state will stop responding to
variations in error/ input ( simply because gradient is 0, nothing changes ). This is called dying ReLu problem.
• The range of ReLu is [0, inf). This means it can blow up the activation.

Further reading
• Deep Sparse Rectifier Neural Networks Glorot et al., (2011)

• Yes You Should Understand Backprop, Karpathy (2016)
13.4 LeakyReLU
LeakyRelu is a variant of ReLU. Instead of being 0 when 𝑧 < 0, a leaky ReLU allows a small, non-zero, constant
gradient 𝛼 (Normally, 𝛼 = 0.01). However, the consistency of the benefit across tasks is presently unclear.1
Function Derivative
z>0
𝛼𝑧 z <= 0
{︂ }︂
1 𝑧>0
𝑅′ (𝑧) =
𝛼 𝑧<0
def leakyrelu(z, alpha):

return max(alpha * z, z)
def leakyrelu_prime(z, alpha):

return 1 if z > 0 else alpha
Pros
• Leaky ReLUs are one attempt to fix the “dying ReLU” problem by having a small negative slope (of 0.01, or
so).
1 http://cs231n.github.io/neural-networks-1/
13.4. LeakyReLU 101

Cons
• As it possess linearity, it can’t be used for the complex Classification. It lags behind the Sigmoid and Tanh for
some of the use cases.
Further reading
• Delving Deep into Rectifiers: Surpassing Human-Level Performance on ImageNet Classification, Kaiming He
et al. (2015)
13.5 Sigmoid
Sigmoid takes a real value as input and outputs another value between 0 and 1. It’s easy to work with and has all the
nice properties of activation functions: it’s non-linear, continuously differentiable, monotonic, and has a fixed output
range.
Function Derivative
1
𝑆(𝑧) = 𝑆 ′ (𝑧) = 𝑆(𝑧) · (1 − 𝑆(𝑧))
1 + 𝑒−𝑧
def sigmoid(z): def sigmoid_prime(z):

return 1.0 / (1 + np.exp(-z)) return sigmoid(z) * (1-sigmoid(z))
Pros
• It is nonlinear in nature. Combinations of this function are also nonlinear!

• It will give an analog activation unlike step function.
• It has a smooth gradient too.
• It’s good for a classifier.

• The output of the activation function is always going to be in range (0,1) compared to (-inf, inf) of linear
function. So we have our activations bound in a range. Nice, it won’t blow up the activations then.
Cons
• Towards either end of the sigmoid function, the Y values tend to respond very less to changes in X.
• It gives rise to a problem of “vanishing gradients”.
• Its output isn’t zero centered. It makes the gradient updates go too far in different directions. 0 < output < 1,
and it makes optimization harder.
• Sigmoids saturate and kill gradients.
• The network refuses to learn further or is drastically slow ( depending on use case and until gradient /computa-
tion gets hit by floating point value limits ).
Further reading
• Yes You Should Understand Backprop, Karpathy (2016)
13.6 Tanh
Tanh squashes a real-valued number to the range [-1, 1]. It’s non-linear. But unlike Sigmoid, its output is zero-centered.
Therefore, in practice the tanh non-linearity is always preferred to the sigmoid nonlinearity.1
Function Derivative
𝑒𝑧 − 𝑒−𝑧 𝑡𝑎𝑛ℎ′ (𝑧) = 1 − 𝑡𝑎𝑛ℎ(𝑧)2

𝑡𝑎𝑛ℎ(𝑧) =
𝑒𝑧 + 𝑒−𝑧
def tanh(z): def tanh_prime(z):

return (np.exp(z) - np.exp(-z)) / return 1 - np.power(tanh(z), 2)
˓→ (np.exp(z) + np.exp(-z))
13.6. Tanh 103

Pros
• The gradient is stronger for tanh than sigmoid ( derivatives are steeper).
Cons
• Tanh also has the vanishing gradient problem.
13.7 Softmax
Softmax function calculates the probabilities distribution of the event over ‘n’ different events. In general way of
saying, this function will calculate the probabilities of each target class over all possible target classes. Later the
calculated probabilities will be helpful for determining the target class for the given inputs.
References

CHAPTER 14
Layers
• BatchNorm
• Convolution
• Dropout
• Linear
• LSTM
• Pooling
• RNN
14.1 BatchNorm
BatchNorm accelerates convergence by reducing internal covariate shift inside each batch. If the individual observa-
tions in the batch are widely different, the gradient updates will be choppy and take longer to converge.
The batch norm layer normalizes the incoming activations and outputs a new batch where the mean equals 0 and
standard deviation equals 1. It subtracts the mean and divides by the standard deviation of the batch.
Code
Code example from Agustinus Kristiadi

def BatchNorm():
# From https://wiseodd.github.io/techblog/2016/07/04/batchnorm/
# TODO: Add doctring for variable names. Add momentum to init.
def __init__(self):
pass
105
def forward(self, X, gamma, beta):

mu = np.mean(X, axis=0)
var = np.var(X, axis=0)
X_norm = (X - mu) / np.sqrt(var + 1e-8)

out = gamma * X_norm + beta
cache = (X, X_norm, mu, var, gamma, beta)
return out, cache, mu, var
def backward(self, dout, cache):

X, X_norm, mu, var, gamma, beta = cache
N, D = X.shape
X_mu = X - mu
std_inv = 1. / np.sqrt(var + 1e-8)
dX_norm = dout * gamma

dvar = np.sum(dX_norm * X_mu, axis=0) * -.5 * std_inv**3
dmu = np.sum(dX_norm * -std_inv, axis=0) + dvar * np.mean(-2. * X_mu, axis=0)
dX = (dX_norm * std_inv) + (dvar * 2 * X_mu / N) + (dmu / N)

dgamma = np.sum(dout * X_norm, axis=0)
dbeta = np.sum(dout, axis=0)
return dX, dgamma, dbeta
Further reading
• Original Paper
• Implementing BatchNorm in Neural Net
• Understanding the backward pass through Batch Norm
14.2 Convolution
In CNN, a convolution is a linear operation that involves multiplication of weight (kernel/filter) with the input and it
does most of the heavy lifting job.
Convolution layer consists of 2 major component 1. Kernel(Filter) 2. Stride
1. Kernel (Filter): A convolution layer can have more than one filter. The size of the filter should be smaller than
the size of input dimension. It is intentional as it allows filter to be applied multiple times at difference point
(position) on the input.Filters are helpful in understanding and identifying important features from given input.
By applying different filters (more than one filter) on the same input helps in extracting different features from
given input. Output from multiplying filter with the input gives Two dimensional array. As such, the output
array from this operation is called “Feature Map”.
2. Stride: This property controls the movement of filter over input. when the value is set to 1, then filter moves 1
column at a time over input. When the value is set to 2 then the filer jump 2 columns at a time as filter moves
over the input.
106 Chapter 14. Layers

Code
# this code demonstate on how Convolution works

# Assume we have a image of 4 X 4 and a filter fo 2 X 2 and Stride = 1
def conv_filter_ouput(input_img_section,filter_value):
# this method perfromas the multiplication of input and filter
# returns singular value
value = 0
for i in range(len(filter_value)):
for j in range(len(filter_value[0])):
value = value + (input_img_section[i][j]*filter_value[i][j])
return value
img_input = [[260.745, 261.332, 112.27 , 262.351],

[260.302, 208.802, 139.05 , 230.709],
[261.775, 93.73 , 166.118, 122.847],
[259.56 , 232.038, 262.351, 228.937]]
filter = [[1,0],
[0,1]]
filterX,filterY = len(filter),len(filter[0])
filtered_result = []
for i in range(0,len(img_mx)-filterX+1):
clm = []
for j in range(0,len(img_mx[0])-filterY+1):
clm.append(conv_filter_ouput(img_mx[i:i+filterX,j:j+filterY],filter))
filtered_result.append(clm)
print(filtered_result)
Further reading
• cs231n reference
14.2. Convolution 107

14.3 Dropout
A dropout layer takes the output of the previous layer’s activations and randomly sets a certain fraction (dropout rate)
of the activatons to 0, cancelling or ‘dropping’ them out.
It is a common regularization technique used to prevent overfitting in Neural Networks.
The dropout rate is the tunable hyperparameter that is adjusted to measure performance with different values. It is
typically set between 0.2 and 0.5 (but may be arbitrarily set).
Dropout is only used during training; At test time, no activations are dropped, but scaled down by a factor of dropout
rate. This is to account for more units being active during test time than training time.
For example:
• A layer in a neural net outputs a tensor (matrix) A of shape (batch_size, num_features).
• The dropout rate of the layer is set to 0.5 (50%).
• A random 50% of the values in A will be set to 0.
• These will then be multiplied with the weight matrix to form the inputs to the next layer.
The premise behind dropout is to introduce noise into a layer in order to disrupt any interdependent learning or
coincidental patterns that may occur between units in the layer, that aren’t significant.
Code
# layer_output is a 2D numpy matrix of activations
layer_output *= np.random.randint(0, high=2, size=layer_output.shape) # dropping out

˓→values

# scaling up by dropout rate during TRAINING time, so no scaling needs to be done at

˓→test time
layer_output /= 0.5
# OR
layer_output *= 0.5 # Scaling down during TEST time.
This results in the following operation.
All reference, images and code examples, unless mentioned otherwise, are from section 4.4.3 of Deep Learning for
Python by François Chollet.
14.4 Linear
14.5 LSTM
14.6 Pooling
Pooling layers often take convolution layers as input. A complicated dataset with many object will require a large
number of filters, each responsible finding pattern in an image so the dimensionally of convolutional layer can get
large. It will cause an increase of parameters, which can lead to over-fitting. Pooling layers are methods for reducing
this high dimensionally. Just like the convolution layer, there is kernel size and stride. The size of the kernel is smaller
than the feature map. For most of the cases the size of the kernel will be 2X2 and the stride of 2. There are mainly two
types of pooling layers.
The first type is max pooling layer. Max pooling layer will take a stack of feature maps (convolution layer) as input.
The value of the node in the max pooling layer is calculated by just the maximum of the pixels contained in the
window.
The other type of pooling layer is the Average Pooling layer. Average pooling layer calculates the average of pixels
contained in the window. Its not used often but you may see this used in applications for which smoothing an image
is preferable.
14.4. Linear 109

Code
def max_pooling(feature_map, size=2, stride=2):

"""
:param feature_map: Feature matrix of shape (height, width, layers)
:param size: size of kernal
:param stride: movement speed of kernal
:return: max-pooled feature vector
"""
pool_shape = (feature_map.shape[0]//stride, feature_map.shape[1]//stride, feature_
˓→map.shape[-1]) #shape of output
pool_out = numpy.zeros(pool_shape)
for layer in range(feature_map.shape[-1]):
#for each layer
row = 0
for r in numpy.arange(0,feature_map.shape[0], stride):
col = 0
for c in numpy.arange(0, feature_map.shape[1], stride):
pool_out[row, col, layer] = numpy.max([feature_map[c:c+size,
˓→r:r+size, layer]])
col = col + 1
row = row +1
return pool_out
14.7 RNN
References

CHAPTER 15
Loss Functions
• Cross-Entropy
• Hinge
• Huber
• Kullback-Leibler
• MAE (L1)
• MSE (L2)
15.1 Cross-Entropy
Cross-entropy loss, or log loss, measures the performance of a classification model whose output is a probability value
between 0 and 1. Cross-entropy loss increases as the predicted probability diverges from the actual label. So predicting
a probability of .012 when the actual observation label is 1 would be bad and result in a high loss value. A perfect
model would have a log loss of 0.
111
The graph above shows the range of possible loss values given a true observation (isDog = 1). As the predicted
probability approaches 1, log loss slowly decreases. As the predicted probability decreases, however, the log loss
increases rapidly. Log loss penalizes both types of errors, but especially those predictions that are confident and
wrong!
Cross-entropy and log loss are slightly different depending on context, but in machine learning when calculating error
rates between 0 and 1 they resolve to the same thing.
Code
def CrossEntropy(yHat, y):

if y == 1:
return -log(yHat)
else:
return -log(1 - yHat)
Math
In binary classification, where the number of classes 𝑀 equals 2, cross-entropy can be calculated as:
−(𝑦 log(𝑝) + (1 − 𝑦) log(1 − 𝑝))
112 Chapter 15. Loss Functions

If 𝑀 > 2 (i.e. multiclass classification), we calculate a separate loss for each class label per observation and sum the
result.
𝑀
∑︁
− 𝑦𝑜,𝑐 log(𝑝𝑜,𝑐 )
𝑐=1
Note:
• M - number of classes (dog, cat, fish)
• log - the natural log
• y - binary indicator (0 or 1) if class label 𝑐 is the correct classification for observation 𝑜
• p - predicted probability observation 𝑜 is of class 𝑐
15.2 Hinge
Used for classification.
Code
def Hinge(yHat, y):

return np.max(0, y - (1-2*y)*yHat)
15.3 Huber
Typically used for regression. It’s less sensitive to outliers than the MSE as it treats error as square only inside an
interval.
{︂ 1
ˆ)2
2 (𝑦 − 𝑦 𝑖𝑓 |(𝑦 − 𝑦ˆ)| < 𝛿
𝐿𝛿 =
𝛿((𝑦 − 𝑦ˆ) − 12 𝛿) 𝑜𝑡ℎ𝑒𝑟𝑤𝑖𝑠𝑒
Code
def Huber(yHat, y, delta=1.):

return np.where(np.abs(y-yHat) < delta,.5*(y-yHat)**2 , delta*(np.abs(y-yHat)-0.
˓→5*delta))
Further information can be found at Huber Loss in Wikipedia.
15.4 Kullback-Leibler
Code
15.2. Hinge 113

def KLDivergence(yHat, y):

"""
:param yHat:
:param y:
:return: KLDiv(yHat || y)
"""
return np.sum(yHat * np.log((yHat / y)))
15.5 MAE (L1)
Mean Absolute Error, or L1 loss. Excellent overview below [6] and [10].
Code
def L1(yHat, y):

return np.sum(np.absolute(yHat - y)) / y.size
15.6 MSE (L2)
Mean Squared Error, or L2 loss. Excellent overview below [6] and [10].
def MSE(yHat, y):

return np.sum((yHat - y)**2) / y.size
def MSE_prime(yHat, y):

return yHat - y
References
114 Chapter 15. Loss Functions

CHAPTER 16
Optimizers
What is Optimizer ?
It is very important to tweak the weights of the model during the training process, to make our predictions as correct
and optimized as possible. But how exactly do you do that? How do you change the parameters of your model, by
how much, and when?
Best answer to all above question is optimizers. They tie together the loss function and model parameters by updating
the model in response to the output of the loss function. In simpler terms, optimizers shape and mold your model into
its most accurate possible form by futzing with the weights. The loss function is the guide to the terrain, telling the
optimizer when it’s moving in the right or wrong direction.
Below are list of example optimizers
• Adagrad
• Adadelta
• Adam
• Conjugate Gradients
• BFGS
• Momentum
• Nesterov Momentum
• Newton’s Method
• RMSProp
• SGD
Image Credit: CS231n
115
16.1 Adagrad
Adagrad (short for adaptive gradient) adaptively sets the learning rate according to a parameter.
• Parameters that have higher gradients or frequent updates should have slower learning rate so that we do not
overshoot the minimum value.
• Parameters that have low gradients or infrequent updates should faster learning rate so that they get trained
quickly.
• It divides the learning rate by the sum of squares of all previous gradients of the parameter.
• When the sum of the squared past gradients has a high value, it basically divides the learning rate by a high
value, so the learning rate will become less.
• Similarly, if the sum of the squared past gradients has a low value, it divides the learning rate by a lower value,
so the learning rate value will become high.
• This implies that the learning rate is inversely proportional to the sum of the squares of all the previous gradients
of the parameter.
𝜕𝒥 (𝑤𝑡𝑖 )
𝑔𝑡𝑖 =
𝜕𝑊
𝑖
𝜕𝒥 (𝑤𝑡 )
𝑊 = 𝑊 − 𝛼 √︁∑︀
𝑡 2
𝑟=1 (𝑔𝑟𝑖 ) + 𝜀
Note:
• 𝑔𝑡𝑖 - the gradient of a parameter, :math: ‘Theta ‘ at an iteration t.
• 𝛼 - the learning rate
• 𝜖 - very small value to avoid dividing by zero
def Adagrad(data):
gradient_sums = np.zeros(theta.shape[0])
for t in range(num_iterations):
gradients = compute_gradients(data, weights)
gradient_sums += gradients ** 2
gradient_update = gradients / (np.sqrt(gradient_sums + epsilon))
weights = weights - lr * gradient_update
return weights
16.2 Adadelta
AdaDelta belongs to the family of stochastic gradient descent algorithms, that provide adaptive techniques for hyper-
parameter tuning. Adadelta is probably short for ‘adaptive delta’, where delta here refers to the difference between the
current weight and the newly updated weight.
The main disadvantage in Adagrad is its accumulation of the squared gradients. During the training process, the
accumulated sum keeps growing. From the above formala we can see that, As the accumulated sum increases learning
rate to shrink and eventually become infinitesimally small, at which point the algorithm is no longer able to acquire
additional knowledge.
116 Chapter 16. Optimizers

Adadelta is a more robust extension of Adagrad that adapts learning rates based on a moving window of gradient
updates, instead of accumulating all past gradients. This way, Adadelta continues learning even when many updates
have been done.
With Adadelta, we do not even need to set a default learning rate, as it has been eliminated from the update rule.
Implementation is something like this,
𝑣𝑡 = 𝜌𝑣𝑡−1 + (1 − 𝜌)∇2𝜃 𝐽(𝜃)

√
𝑤𝑡 + 𝜖
∆𝜃 = √ ∇𝜃 𝐽(𝜃)
𝑣𝑡 + 𝜖
𝜃 = 𝜃 − 𝜂∆𝜃
𝑤𝑡 = 𝜌𝑤𝑡−1 + (1 − 𝜌)∆𝜃2
def Adadelta(weights, sqrs, deltas, rho, batch_size):

eps_stable = 1e-5
for weight, sqr, delta in zip(weights, sqrs, deltas):
g = weight.grad / batch_size
sqr[:] = rho * sqr + (1. - rho) * nd.square(g)
cur_delta = nd.sqrt(delta + eps_stable) / nd.sqrt(sqr + eps_stable) * g
delta[:] = rho * delta + (1. - rho) * cur_delta * cur_delta
# update weight in place.
weight[:] -= cur_delta
16.3 Adam
Adaptive Moment Estimation (Adam) combines ideas from both RMSProp and Momentum. It computes adaptive
learning rates for each parameter and works as follows.
• First, it computes the exponentially weighted average of past gradients (𝑣𝑑𝑊 ).
• Second, it computes the exponentially weighted average of the squares of past gradients (𝑠𝑑𝑊 ).
𝑐𝑜𝑟𝑟𝑒𝑐𝑡𝑒𝑑
• Third, these averages have a bias towards zero and to counteract this a bias correction is applied (𝑣𝑑𝑊 ,
𝑐𝑜𝑟𝑟𝑒𝑐𝑡𝑒𝑑
𝑠𝑑𝑊 ).
• Lastly, the parameters are updated using the information from the calculated averages.
𝜕𝒥
𝑣𝑑𝑊 = 𝛽1 𝑣𝑑𝑊 + (1 − 𝛽1 )
𝜕𝑊
𝜕𝒥 2
𝑠𝑑𝑊 = 𝛽2 𝑠𝑑𝑊 + (1 − 𝛽2 )( )
𝜕𝑊
𝑐𝑜𝑟𝑟𝑒𝑐𝑡𝑒𝑑 𝑣𝑑𝑊
𝑣𝑑𝑊 =
1 − (𝛽1 )𝑡
𝑠𝑑𝑊
𝑠𝑐𝑜𝑟𝑟𝑒𝑐𝑡𝑒𝑑
𝑑𝑊 =
1 − (𝛽1 )𝑡
𝑣 𝑐𝑜𝑟𝑟𝑒𝑐𝑡𝑒𝑑
𝑊 = 𝑊 − 𝛼 √︁ 𝑑𝑊
𝑑𝑊 +𝜀
Note:
• 𝑣𝑑𝑊 - the exponentially weighted average of past gradients
• 𝑠𝑑𝑊 - the exponentially weighted average of past squares of gradients
16.3. Adam 117

• 𝛽1 - hyperparameter to be tuned
• 𝛽2 - hyperparameter to be tuned
𝜕𝒥
• 𝜕𝑊 - cost gradient with respect to current layer
• 𝑊 - the weight matrix (parameter to be updated)
16.4 Conjugate Gradients
16.5 BFGS
16.6 Momentum
Used in conjunction Stochastic Gradient Descent (sgd) or Mini-Batch Gradient Descent, Momentum takes into account
past gradients to smooth out the update. This is seen in variable 𝑣 which is an exponentially weighted average of the
gradient on previous steps. This results in minimizing oscillations and faster convergence.
𝜕𝒥
𝑣𝑑𝑊 = 𝛽𝑣𝑑𝑊 + (1 − 𝛽)
𝜕𝑊
𝑊 = 𝑊 − 𝛼𝑣𝑑𝑊
Note:
• 𝑣 - the exponentially weighted average of past gradients
𝜕𝒥
• 𝜕𝑊 - cost gradient with respect to current layer weight tensor
• 𝑊 - weight tensor
• 𝛽 - hyperparameter to be tuned
16.7 Nesterov Momentum
16.8 Newton’s Method

16.9 RMSProp
Another adaptive learning rate optimization algorithm, Root Mean Square Prop (RMSProp) works by keeping an
exponentially weighted average of the squares of past gradients. RMSProp then divides the learning rate by this
average to speed up convergence.
𝜕𝒥 2
𝑠𝑑𝑊 = 𝛽𝑠𝑑𝑊 + (1 − 𝛽)( )
𝜕𝑊
𝜕𝒥
𝜕𝑊
𝑊 = 𝑊 − 𝛼 √︁
𝑑𝑊 +𝜀
Note:
• 𝑠 - the exponentially weighted average of past squares of gradients
𝜕𝒥
• 𝜕𝑊 - cost gradient with respect to current layer weight tensor
• 𝑊 - weight tensor
• 𝛽 - hyperparameter to be tuned
16.10 SGD
SGD stands for Stochastic Gradient Descent.In Stochastic Gradient Descent, a few samples are selected randomly
instead of the whole data set for each iteration. In Gradient Descent, there is a term called “batch” which denotes the
total number of samples from a dataset that is used for calculating the gradient for each iteration. In typical Gradient
Descent optimization, like Batch Gradient Descent, the batch is taken to be the whole dataset. Although, using the
whole dataset is really useful for getting to the minima in a less noisy or less random manner, but the problem arises
when our datasets get really huge.
This problem is solved by Stochastic Gradient Descent. In SGD, it uses only a single sample to perform each iteration.
The sample is randomly shuffled and selected for performing the iteration.
Since only one sample from the dataset is chosen at random for each iteration, the path taken by the algorithm to reach
the minima is usually noisier than your typical Gradient Descent algorithm. But that doesn’t matter all that much
because the path taken by the algorithm does not matter, as long as we reach the minima and with significantly shorter
training time.
def SGD(data, batch_size, lr):

N = len(data)
np.random.shuffle(data)
mini_batches = np.array([data[i:i+batch_size]
for i in range(0, N, batch_size)])
for X,y in mini_batches:
backprop(X, y, lr)
References
16.9. RMSProp 119


CHAPTER 17
Regularization
• Data Augmentation
• Dropout
• Early Stopping
• Ensembling
• Injecting Noise
• L1 Regularization
• L2 Regularization
What is overfitting?
From Wikipedia overfitting is,

The production of an analysis that corresponds too closely or exactly to a particular set of data, and may therefore fail
to fit additional data or predict future observations reliably
What is Regularization?
It is a Techniques for combating overfitting and improving training.
17.1 Data Augmentation
Having more data (dataset / samples) is a best way to get better consistent estimators (ML model). In the real world
getting a large volume of useful data for training a model is cumbersome and labelling is an extremely tedious task.
121
Either labelling requires more manual annotation, example - For creating a better image classifier we can use Mturk
and involve more man power to generate dataset or doing survey in social media and asking people to participate and
generate dataset. Above process can yield good dataset however those are difficult to carry and expensive. Having
small dataset will lead to the well know Over fitting problem.
Data Augmentation is one of the interesting regularization technique to resolve the above problem. The concept is
very simple, this technique generates new training data from given original dataset. Dataset Augmentation provides a
cheap and easy way to increase the amount of your training data.
This technique can be used for both NLP and CV.
In CV we can use the techniques like Jitter, PCA and Flipping. Similarly in NLP we can use the techniques like
Synonym Replacement,Random Insertion, Random Deletion and Word Embeddings.
It is worth knowing that Keras’ provided ImageDataGenerator for generating Data Augmentation.
Sample code for random deletion
def random_deletion(words, p):

"""
Randomly delete words from the sentence with probability p
"""
#obviously, if there's only one word, don't delete it

if len(words) == 1:
return words
#randomly delete words with probability p

new_words = []
for word in words:
r = random.uniform(0, 1)
if r > p:
new_words.append(word)
#if you end up deleting all words, just return a random word
if len(new_words) == 0:
rand_int = random.randint(0, len(words)-1)
return [words[rand_int]]
return new_words
Furthermore, when comparing two machine learning algorithms train both with either augmented or non-augmented
dataset. Otherwise, no subjective decision can be made on which algorithm performed better
Further reading
• NLP Data Augmentation

• CV Data Augmentation
• Regularization
17.2 Dropout
122 Chapter 17. Regularization

What is Dropout?
Dropout is a regularization technique for reducing overfitting in neural networks by preventing complex co-adaptations
on training data
Dropout is a technique where randomly selected neurons are ignored during training. They are “dropped-out” ran-
domly. This means that their contribution to the activation of downstream neurons is temporally removed on the
forward pass and any weight updates are not applied to the neuron on the backward pass.
Simply put, It is the process of ignoring some of the neurons in particular forward or backward pass.
Dropout can be easily implemented by randomly selecting nodes to be dropped-out with a given probability (e.g. .1%)
each weight update cycle.
Most importantly Dropout is only used during the training of a model and is not used when evaluating the model.
image from https://www.cs.toronto.edu/~hinton/absps/JMLRdropout.pdf
import numpy as np
A = np.arange(20).reshape((5,4))
print("Given input: ")

print(A)
def dropout(X, drop_probability):

keep_probability = 1 - drop_probability
mask = np.random.uniform(0, 1.0, X.shape) < keep_probability
if keep_probability > 0.0:
scale = (1/keep_probability)
else:
scale = 0.0
return mask * X * scale
print("\n After Dropout: ")

print(dropout(A,0.5))
output from above code
17.2. Dropout 123

Given input:
[[ 0 1 2 3]
[ 4 5 6 7]
[ 8 9 10 11]
[12 13 14 15]
[16 17 18 19]]
After Dropout:
[[ 0. 2. 0. 0.]
[ 8. 0. 0. 14.]
[16. 18. 0. 22.]
[24. 0. 0. 0.]
[32. 34. 36. 0.]]
Further reading
• Dropout https://www.cs.toronto.edu/~hinton/absps/JMLRdropout.pdf
17.3 Early Stopping
One of the biggest problem in training neural network is how long to train the model.
Training too little will lead to underfit in train and test sets. Traning too much will have the overfit in training set and
poor result in test sets.
Here the challenge is to train the network long enough that it is capable of learning the mapping from inputs to outputs,
but not training the model so long that it overfits the training data.
One possible solution to solve this problem is to treat the number of training epochs as a hyperparameter and train the
model multiple times with different values, then select the number of epochs that result in the best accuracy on the
train or a holdout test dataset, But the problem is it requires multiple models to be trained and discarded.
Clearly, after ‘t’ epochs, the model starts overfitting. This is clear by the increasing gap between the train and the
validation error in the above plot.
One alternative technique to prevent overfitting is use validation error to decide when to stop. This approach is called
Early Stopping.
While building the model, it is evaluated on the holdout validation dataset after each epoch. If the accuracy of the
model on the validation dataset starts to degrade (e.g. loss begins to increase or accuracy begins to decrease), then the
training process is stopped. This process is called Early stopping.

Python implementation for Early stopping,
def early_stopping(theta0, (x_train, y_train), (x_valid, y_valid), n = 1, p = 100):

""" The early stopping meta-algorithm for determining the best amount of time to
˓→train.
REF: Algorithm 7.1 in deep learning book.
Parameters:
n: int; Number of steps between evaluations.
p: int; "patience", the number of evaluations to observe worsening validataion
˓→set.
theta0: Network; initial network.

x_train: iterable; The training input set.
y_train: iterable; The training output set.
x_valid: iterable; The validation input set.
y_valid: iterable; The validation output set.
Returns:
theta_prime: Network object; The output network.
i_prime: int; The number of iterations for the output network.
v: float; The validation error for the output network.
"""
# Initialize variables
theta = theta0.clone() # The active network
i = 0 # The number of training steps taken
j = 0 # The number of evaluations steps since last update of
˓→theta_prime
v = np.inf # The best evaluation error observed thusfar

theta_prime = theta.clone() # The best network found thusfar
i_prime = i # The index of theta_prime
while j < p:
# Update theta by running the training algorithm for n steps
for _ in range(n):
theta.train(x_train, y_train)
# Update Values
i += n
v_new = theta.error(x_valid, y_valid)
# If better validation error, then reset waiting time, save the network, and
˓→ update the best error value
if v_new < v:
j = 0
theta_prime = theta.clone()
i_prime = i
v = v_new
# Otherwise, update the waiting time

else:
j += 1
return theta_prime, i_prime, v
Further reading
• Regularization
17.3. Early Stopping 125

17.4 Ensembling
Ensemble methods combine several machine learning techniques into one predictive model. There are a few different
methods for ensembling, but the two most common are:
Bagging
• Bagging stands for bootstrap aggregation. One way to reduce the variance of an estimate is to average together
multiple estimates.
• It trains a large number of “strong” learners in parallel.
• A strong learner is a model that’s relatively unconstrained.
• Bagging then combines all the strong learners together in order to “smooth out” their predictions.
Boosting
• Boosting refers to a family of algorithms that are able to convert weak learners to strong learners.
• Each one in the sequence focuses on learning from the mistakes of the one before it.
• Boosting then combines all the weak learners into a single strong learner.
Bagging uses complex base models and tries to “smooth out” their predictions, while boosting uses simple base models
and tries to “boost” their aggregate complexity.
17.5 Injecting Noise
Noise is often introduced to the inputs as a dataset augmentation strategy. When we have a small dataset the network
may effectively memorize the training dataset. Instead of learning a general mapping from inputs to outputs, the model
may learn the specific input examples and their associated outputs. One approach for improving generalization error
and improving the structure of the mapping problem is to add random noise.
Adding noise means that the network is less able to memorize training samples because they are changing all of the
time, resulting in smaller network weights and a more robust network that has lower generalization error.
Noise is only added during training. No noise is added during the evaluation of the model or when the model is used
to make predictions on new data.
Random noise can be added to other parts of the network during training. Some examples include:
Noise Injection on Weights
• Noise added to weights can be interpreted as a more traditional form of regularization.

• In other words, it pushes the model to be relatively insensitive to small variations in the weights, finding points
that are not merely minima, but minima surrounded by flat regions.
Noise Injection on Outputs
• In the real world dataset, We can expect some amount of mistakes in the output labels. One way to remedy this
is to explicitly model the noise on labels.

• An example for Noise Injection on Outputs is label smoothing
Further reading
• Regularization
17.6 L1 Regularization
A regression model that uses L1 regularization technique is called Lasso Regression.
Mathematical formula for L1 Regularization.
Let’s define a model to see how L1 Regularization works. For simplicity, We define a simple linear regression model
Y with one independent variable.
In this model, W represent Weight, b represent Bias.
𝑊 = 𝑤1 , 𝑤2 ...𝑤𝑛
𝑋 = 𝑥1 , 𝑥2 ...𝑥𝑛
and the predicted result is 𝑌̂︀
𝑌̂︀ = 𝑤1 𝑥1 + 𝑤2 𝑥2 + ...𝑤𝑛 𝑥𝑛 + 𝑏
Following formula calculates the error without Regularization function
𝐿𝑜𝑠𝑠 = 𝐸𝑟𝑟𝑜𝑟(𝑌, 𝑌̂︀ )

Following formula calculates the error With L1 Regularization function
𝑛
∑︁
𝐿𝑜𝑠𝑠 = 𝐸𝑟𝑟𝑜𝑟(𝑌 − 𝑌̂︀ ) + 𝜆 |𝑤𝑖 |
1
Note: Here, If the value of lambda is Zero then above Loss function becomes Ordinary Least Square whereas very
large value makes the coefficients (weights) zero hence it under-fits.
One thing to note is that |𝑤| is differentiable when w!=0 as shown below,
{︃
d|𝑤| 1 𝑤>0
=
d𝑤 −1 𝑤 < 0
To understand the Note above,
Let’s substitute the formula in finding new weights using Gradient Descent optimizer.
𝜕𝐿1
𝑤𝑛𝑒𝑤 = 𝑤 − 𝜂
𝜕𝑤
When we apply the L1 in above formula it becomes,
d|𝑤|
𝑤𝑛𝑒𝑤 = 𝑤 − 𝜂.(𝐸𝑟𝑟𝑜𝑟(𝑌, 𝑌̂︀ ) + 𝜆 )
{︃ d𝑤
𝑤 − 𝜂.(𝐸𝑟𝑟𝑜𝑟(𝑌, 𝑌̂︀ ) + 𝜆) 𝑤 > 0
=
𝑤 − 𝜂.(𝐸𝑟𝑟𝑜𝑟(𝑌, 𝑌̂︀ ) − 𝜆) 𝑤 < 0
From the above formula,
17.6. L1 Regularization 127

• If w is positive, the regularization parameter 𝜆 > 0 will push w to be less positive, by subtracting 𝜆 from w.
• If w is negative, the regularization parameter 𝜆 < 0 will push w to be less negative, by adding 𝜆 to w. hence this
has the effect of pushing w towards 0.
Simple python implementation
def update_weights_with_l1_regularization(features, targets, weights, lr,lambda):

'''
Features:(200, 3)
Targets: (200, 1)
Weights:(3, 1)
'''

x1 = features[:,0]
x2 = features[:,1]
x3 = features[:,2]


weights[0][0] = (weights[0][0] - lr * np.mean(d_w1) - lambda) if weights[0][0] >

˓→0 else (weights[0][0] - lr * np.mean(d_w1) + lambda)

return weights
Use Case
L1 Regularization (or varient of this concept) is a model of choice when the number of features are high, Since it
provides sparse solutions. We can get computational advantage as the features with zero coefficients can simply be
ignored.
Further reading
• Linear Regression
17.7 L2 Regularization
A regression model that uses L2 regularization technique is called Ridge Regression. Main difference between L1 and
L2 regularization is, L2 regularization uses “squared magnitude” of coefficient as penalty term to the loss function.

Mathematical formula for L2 Regularization.
Let’s define a model to see how L2 Regularization works. For simplicity, We define a simple linear regression model
Y with one independent variable.
In this model, W represent Weight, b represent Bias.
𝑊 = 𝑤1 , 𝑤2 ...𝑤𝑛
𝑋 = 𝑥1 , 𝑥2 ...𝑥𝑛
and the predicted result is 𝑌̂︀
𝑌̂︀ = 𝑤1 𝑥1 + 𝑤2 𝑥2 + ...𝑤𝑛 𝑥𝑛 + 𝑏
Following formula calculates the error without Regularization function
𝐿𝑜𝑠𝑠 = 𝐸𝑟𝑟𝑜𝑟(𝑌, 𝑌̂︀ )
Following formula calculates the error With L2 Regularization function

𝑛
∑︁
𝐿𝑜𝑠𝑠 = 𝐸𝑟𝑟𝑜𝑟(𝑌 − 𝑌̂︀ ) + 𝜆 𝑤𝑖2
1
Note: Here, if lambda is zero then you can imagine we get back OLS. However, if lambda is very large then it will
add too much weight and it leads to under-fitting.
To understand the Note above,

Let’s substitute the formula in finding new weights using Gradient Descent optimizer.
𝜕𝐿2
𝑤𝑛𝑒𝑤 = 𝑤 − 𝜂
𝜕𝑤
When we apply the L2 in above formula it becomes,
𝜕𝐿2
𝑤𝑛𝑒𝑤 = 𝑤 − 𝜂.(𝐸𝑟𝑟𝑜𝑟(𝑌, 𝑌̂︀ ) + 𝜆 )
𝜕𝑤
= 𝑤 − 𝜂.(𝐸𝑟𝑟𝑜𝑟(𝑌, 𝑌̂︀ ) + 2𝜆𝑤)
Simple python implementation

def update_weights_with_l2_regularization(features, targets, weights, lr,lambda):
'''
Features:(200, 3)
Targets: (200, 1)
Weights:(3, 1)
'''

x1 = features[:,0]
x2 = features[:,1]
x3 = features[:,2]

17.7. L2 Regularization 129



weights[0][0] = weights[0][0] - lr * np.mean(d_w1) - 2 * lambda * weights[0][0]

return weights
Use Case
L2 regularization can address the multicollinearity problem by constraining the coefficient norm and keeping all the
variables. L2 regression can be used to estimate the predictor importance and penalize predictors that are not important.
One issue with co-linearity is that the variance of the parameter estimate is huge. In cases where the number of features
are greater than the number of observations, the matrix used in the OLS may not be invertible but Ridge Regression
enables this matrix to be inverted.
Further reading
• Ridge Regression
References

CHAPTER 18
Architectures
• Autoencoder
• CNN
• GAN
• MLP
• RNN
• VAE
18.1 Autoencoder
TODO: Description of Autoencoder use case and basic architecture. Figure from [1].
131
Model
An example implementation in PyTorch.

class Autoencoder(nn.Module):
def __init__(self, in_shape):
super().__init__()
c,h,w = in_shape
self.encoder = nn.Sequential(
nn.Linear(c*h*w, 128),
nn.ReLU(),
nn.Linear(128, 64),
nn.ReLU(),
nn.Linear(64, 12),
nn.ReLU()
)
self.decoder = nn.Sequential(
nn.Linear(12, 64),
nn.ReLU(),
nn.Linear(64, 128),
nn.ReLU(),
nn.Linear(128, c*h*w),
nn.Sigmoid()
)
def forward(self, x):

bs,c,h,w = x.size()
x = x.view(bs, -1)
x = self.encoder(x)
x = self.decoder(x)
x = x.view(bs, c, h, w)
return x
Training
def train(net, loader, loss_func, optimizer):

net.train()
132 Chapter 18. Architectures


for inputs, _ in loader:
inputs = Variable(inputs)
output = net(inputs)
loss = loss_func(output, inputs)
optimizer.zero_grad()
loss.backward()
optimizer.step()
Further reading
• Convolutional Autoencoders
• Deep Learning Book
18.2 CNN
TODO: Description of CNN use case and basic architecture. Figure from [2].
Model
Training
Further reading
• CS231 Convolutional Networks

18.2. CNN 133

18.3 GAN
A Generative Adversarial Network (GAN) is a type of network which creates novel tensors (often images, voices,
etc.). The generative portion of the architecture competes with the discriminator part of the architecture in a zero-sum
game. The goal of the generative network is to create novel tensors which the adversarial network attempts to classify
as real or fake. The goal of the generative network is generate tensors where the discriminator network determines that
the tensor has a 50% chance of being fake and a 50% chance of being real.
Figure from [3].
Model
Generator
class Generator(nn.Module):
def __init__(self):
super()
self.net = nn.Sequential(
nn.ConvTranspose2d( 200, 32 * 8, 4, 1, 0, bias=False),
nn.BatchNorm2d(32 * 8),
nn.ReLU(),
nn.ConvTranspose2d(32 * 8, 32 * 4, 4, 2, 1, bias=False),
nn.ReLU(),
nn.ConvTranspose2d( 32 * 4, 32 * 2, 4, 2, 1, bias=False),
nn.ReLU(),
nn.ConvTranspose2d( 32 * 2, 32, 4, 2, 1, bias=False),
nn.BatchNorm2d(32),
nn.ReLU(),
nn.ConvTranspose2d( 32, 1, 4, 2, 1, bias=False),
nn.Tanh()
)
def forward(self, tens):
return self.net(tens)

Discriminator
class Discriminator(nn.Module):
def __init__(self):
super()
self.net = nn.Sequential(
nn.Conv2d(1, 32, 4, 2, 1, bias=False),
nn.LeakyReLU(0.2),
nn.Conv2d(32, 32 * 2, 4, 2, 1, bias=False),
nn.LeakyReLU(0.2),
nn.Conv2d(32 * 2, 32 * 4, 4, 2, 1, bias=False),
nn.LeakyReLU(0.2),
# state size. (32*4) x 8 x 8
nn.Conv2d(32 * 4, 32 * 8, 4, 2, 1, bias=False),
nn.LeakyReLU(0.2),
# state size. (32*8) x 4 x 4
nn.Conv2d(32 * 8, 1, 4, 1, 0, bias=False),
nn.Sigmoid()
)
def forward(self, tens):

return self.net(tens)
Training
def train(netD, netG, loader, loss_func, optimizerD, optimizerG, num_epochs):

netD.train()
netG.train()
device = "cuda:0" if torch.cuda.is_available() else "cpu"
for epoch in range(num_epochs):
for i, data in enumerate(loader, 0):
netD.zero_grad()
realtens = data[0].to(device)
b_size = realtens.size(0)
label = torch.full((b_size,), 1, dtype=torch.float, device=device) # gen
˓→labels
output = netD(realtens)
errD_real = loss_func(output, label)
errD_real.backward() # backprop discriminator fake and real based on label
noise = torch.randn(b_size, 200, 1, 1, device=device)
fake = netG(noise)
label.fill_(0)
output = netD(fake.detach()).view(-1)
errD_fake = loss_func(output, label)
errD_fake.backward() # backprop discriminator fake and real based on label
errD = errD_real + errD_fake # discriminator error
optimizerD.step()
netG.zero_grad()
label.fill_(1)
output = netD(fake).view(-1)
errG = loss_func(output, label) # generator error
errG.backward()
optimizerG.step()
18.3. GAN 135

Further reading
• Generative Adversarial Networks

• PyTorch DCGAN Example
• Original Paper
18.4 MLP
A Multi Layer Perceptron (MLP) is a neural network with only fully connected layers. Figure from [5].
Model
An example implementation on FMNIST dataset in PyTorch. Full Code

1. The input to the network is a vector of size 28*28 i.e.(image from FashionMNIST dataset of dimension 28*28
pixels flattened to sigle dimension vector).
2. 2 fully connected hidden layers.
3. Output layer with 10 outputs.(10 classes)
class MLP(nn.Module):
def __init__(self):
super(MLP,self).__init__()
# define layers
self.fc1 = nn.Linear(in_features=28*28, out_features=500)
self.fc2 = nn.Linear(in_features=500, out_features=200)
self.fc3 = nn.Linear(in_features=200, out_features=100)
self.out = nn.Linear(in_features=100, out_features=10)

def forward(self, t):

# fc1 make input 1 dimentional
t = t.view(-1,28*28)
t = self.fc1(t)
t = F.relu(t)
# fc2
t = self.fc2(t)
t = F.relu(t)
# fc3
t = self.fc3(t)
t = F.relu(t)
# output
t = self.out(t)
return t
Training
def train(net, loader, loss_func, optimizer):

net.train()
n_batches = len(loader)
for inputs, targets in loader:
targets = Variable(targets)
output = net(inputs)
loss = loss_func(output, targets)
loss.backward()
optimizer.step()
# print statistics
running_loss = loss.item()
print('Training loss: %.3f' %( running_loss))
Evaluating
def main():
train_set = torchvision.datasets.FashionMNIST(
root = './FMNIST',
train = True,
download = False,
transform = transforms.Compose([
transforms.ToTensor()
])
)
mlp = MLP()
loader = torch.utils.data.DataLoader(train_set, batch_size = 1000)
optimizer = optim.Adam(mlp.parameters(), lr=0.01)
loss_func=nn.CrossEntropyLoss()
for i in range(0,15):
train(mlp,loader,loss_func,optimizer)
18.4. MLP 137


print("Finished Training")
torch.save(mlp.state_dict(), "./mlpmodel.pt")
test_set = torchvision.datasets.FashionMNIST(
root = './FMNIST',
train = False,
download = False,
transform = transforms.Compose([
transforms.ToTensor()
])
)
testloader = torch.utils.data.DataLoader(test_set, batch_size=4)
correct = 0
total = 0
with torch.no_grad():
for data in testloader:
images, labels = data
outputs = mlp(images)
_, predicted = torch.max(outputs.data, 1)
total += labels.size(0)
correct += (predicted == labels).sum().item()
print('Accuracy of the network on the 10000 test images: %d %%' % (
100 * correct / total))
Further reading
TODO
18.5 RNN
Description of RNN use case and basic architecture.

Model
class RNN(nn.Module):
def __init__(self, n_classes):
super().__init__()
self.hid_fc = nn.Linear(185, 128)
self.out_fc = nn.Linear(185, n_classes)
self.softmax = nn.LogSoftmax()
def forward(self, inputs, hidden):

inputs = inputs.view(1,-1)
combined = torch.cat([inputs, hidden], dim=1)
hid_out = self.hid_fc(combined)
out = self.out_fc(combined)
out = self.softmax(out)
return out, hid_out
Training
In this example, our input is a list of last names, where each name is a variable length array of one-hot encoded
characters. Our target is is a list of indices representing the class (language) of the name.
18.5. RNN 139

1. For each input name..

2. Initialize the hidden vector
3. Loop through the characters and predict the class
4. Pass the final character’s prediction to the loss function
5. Backprop and update the weights
def train(model, inputs, targets):

for i in range(len(inputs)):
target = Variable(targets[i])
name = inputs[i]
hidden = Variable(torch.zeros(1,128))
model.zero_grad()
for char in name:

input_ = Variable(torch.FloatTensor(char))
pred, hidden = model(input_, hidden)
loss = criterion(pred, target)

loss.backward()
for p in model.parameters():
p.data.add_(-.001, p.grad.data)
Further reading
• Jupyter notebook
18.6 VAE
Autoencoders can encode an input image to a latent vector and decode it, but they can’t generate novel images.
Variational Autoencoders (VAE) solve this problem by adding a constraint: the latent vector representation should
model a unit gaussian distribution. The Encoder returns the mean and variance of the learned gaussian. To generate a
new image, we pass a new mean and variance to the Decoder. In other words, we “sample a latent vector” from the
gaussian and pass it to the Decoder. It also improves network generalization and avoids memorization. Figure from
[4].

Loss Function
The VAE loss function combines reconstruction loss (e.g. Cross Entropy, MSE) with KL divergence.
def vae_loss(output, input, mean, logvar, loss_func):

recon_loss = loss_func(output, input)
kl_loss = torch.mean(0.5 * torch.sum(
torch.exp(logvar) + mean**2 - 1. - logvar, 1))
return recon_loss + kl_loss
Model
An example implementation in PyTorch of a Convolutional Variational Autoencoder.
class VAE(nn.Module):
def __init__(self, in_shape, n_latent):
super().__init__()
self.in_shape = in_shape
self.n_latent = n_latent
c,h,w = in_shape
self.z_dim = h//2**2 # receptive field downsampled 2 times
self.encoder = nn.Sequential(
nn.BatchNorm2d(c),
nn.Conv2d(c, 32, kernel_size=4, stride=2, padding=1), # 32, 16, 16
nn.BatchNorm2d(32),
nn.LeakyReLU(),
nn.Conv2d(32, 64, kernel_size=4, stride=2, padding=1), # 32, 8, 8
nn.BatchNorm2d(64),
nn.LeakyReLU(),
)
self.z_mean = nn.Linear(64 * self.z_dim**2, n_latent)
self.z_var = nn.Linear(64 * self.z_dim**2, n_latent)
self.z_develop = nn.Linear(n_latent, 64 * self.z_dim**2)
self.decoder = nn.Sequential(
18.6. VAE 141


nn.ConvTranspose2d(64, 32, kernel_size=3, stride=2, padding=0),
nn.BatchNorm2d(32),
nn.ReLU(),
nn.ConvTranspose2d(32, 1, kernel_size=3, stride=2, padding=1),
CenterCrop(h,w),
nn.Sigmoid()
)
def sample_z(self, mean, logvar):

stddev = torch.exp(0.5 * logvar)
noise = Variable(torch.randn(stddev.size()))
return (noise * stddev) + mean
def encode(self, x):

x = self.encoder(x)
x = x.view(x.size(0), -1)
mean = self.z_mean(x)
var = self.z_var(x)
return mean, var
def decode(self, z):

out = self.z_develop(z)
out = out.view(z.size(0), 64, self.z_dim, self.z_dim)
out = self.decoder(out)
return out
def forward(self, x):

mean, logvar = self.encode(x)
z = self.sample_z(mean, logvar)
out = self.decode(z)
return out, mean, logvar
Training
def train(model, loader, loss_func, optimizer):

model.train()
for inputs, _ in loader:
output, mean, logvar = model(inputs)

loss = vae_loss(output, inputs, mean, logvar, loss_func)
loss.backward()
optimizer.step()
Further reading
• Original Paper
• VAE Explained

References
18.6. VAE 143


CHAPTER 19
Classification Algorithms
Classification problems is when our output Y is always in categories like positive vs negative in terms of sentiment
analysis, dog vs cat in terms of image classification and disease vs no disease in terms of medical diagnosis.
19.1 Bayesian
Overlaps..
19.2 Decision Trees
Intuitions
Decision tree works by successively splitting the dataset into small segments until the target variable are the same or
until the dataset can no longer be split. It’s a greedy algorithm which make the best decision at the given time without
concern for the global optimality2 .
The concept behind decision tree is straightforward. The following flowchart show a simple email classification system
based on decision tree. If the address is “myEmployer.com”, it will classify it to “Email to read when bored”. Then
if the email contains the word “hockey”, this email will be classified as “Email from friends”. Otherwise, it will be
identified as “Spam: don’t read”. Image source2 .
2 Machine Learning in Action by Peter Harrington
145
Algorithm Explained
There are various kinds of decision tree algorithms such as ID3 (Iterative Dichotomiser 3), C4.5 and CART (Classifi-
cation and Regression Trees). The constructions of decision tree are similar5 :
1. Assign all training instances to the root of the tree. Set current node to root node.
2. Find the split feature and split value based on the split criterion such as information gain, information gain ratio
or gini coefficient.
3. Partition all data instances at the node based on the split feature and threshold value.
4. Denote each partition as a child node of the current node.
5. For each child node:
1. If the child node is “pure” (has instances from only one class), tag it as a leaf and return.
2. Else, set the child node as the current node and recurse to step 2.
ID3 creates a multiway tree. For each node, it trys to find the categorical feature that will yield the largest information
gain for the target variable.
C4.5 is the successor of ID3 and remove the restriction that the feature must be categorical by dynamically define a
discrete attribute that partitions the continuous attribute in the discrete set of intervals.
CART is similar to C4.5. But it differs in that it constructs binary tree and support regression problem3 .
The main differences are shown in the following table:
Dimensions ID3 C4.5 CART

Split Crite- Information Information gain ratio (Normalized infor- Gini coefficient for classification
rion gain mation gain) problems
Types of Fea- Categorical Categorical & numerical features Categorical & numerical features
tures feature
Type of Prob- Classification Classification Classification & regression
lem
Type of Tree Mltiway tree Mltiway tree Binary tree
5 Decision Trees
3 Scikit-learn Documentations: Tree algorithms: ID3, C4.5, C5.0 and CART
146 Chapter 19. Classification Algorithms

Code Implementation
We used object-oriented patterns to create the code for ID3, C4.5 and CART. We will first introduce the base class for
these three algorithms, then we explain the code of CART in details.
First, we create the base class TreeNode class and DecisionTree
class TreeNode:
def __init__(self, data_idx, depth, child_lst=[]):
self.data_idx = data_idx
self.depth = depth
self.child = child_lst
self.label = None
self.split_col = None
self.child_cate_order = None
def set_attribute(self, split_col, child_cate_order=None):

self.split_col = split_col
self.child_cate_order = child_cate_order
def set_label(self, label):

self.label = label
class DecisionTree()
def fit(self, X, y):
"""
X: train data, dimensition [num_sample, num_feature]
y: label, dimension [num_sample, ]
"""
self.data = X
self.labels = y
num_sample, num_feature = X.shape
self.feature_num = num_feature
data_idx = list(range(num_sample))
# Set the root of the tree
self.root = TreeNode(data_idx=data_idx, depth=0, child_lst=[])
queue = [self.root]
while queue:
node = queue.pop(0)
# Check if the terminate criterion has been met
if node.depth>self.max_depth or len(node.data_idx)==1:
# Set the label for the leaf node
self.set_label(node)
else:
# Split the node
child_nodes = self.split_node(node)
if not child_nodes:
self.set_label(node)
else:
queue.extend(child_nodes)
The CART algorithm, when constructing the binary tree, will try searching for the feature and threshold that will
yield the largest gain or the least impurity. The split criterion is a combination of the child nodes’ impurity. For the
child nodes’ impurity, gini coefficient or information gain are adopted in classification. For regression problem, mean-
square-error or mean-absolute-error are used. Example codes are showed below. For more details about the formulas,
please refer to Mathematical formulation for decision tree in scikit-learn documentation
19.2. Decision Trees 147

class CART(DecisionTree):
def get_split_criterion(self, node, child_node_lst):

total = len(node.data_idx)
split_criterion = 0
for child_node in child_node_lst:
impurity = self.get_impurity(child_node.data_idx)
split_criterion += len(child_node.data_idx) / float(total) * impurity
return split_criterion
def get_impurity(self, data_ids):

target_y = self.labels[data_ids]
total = len(target_y)
if self.tree_type == "regression":
res = 0
mean_y = np.mean(target_y)
for y in target_y:
res += (y - mean_y) ** 2 / total
elif self.tree_type == "classification":
if self.split_criterion == "gini":
res = 1
unique_y = np.unique(target_y)
for y in unique_y:
num = len(np.where(target_y==y)[0])
res -= (num/float(total))**2
elif self.split_criterion == "entropy":
unique, count = np.unique(target_y, return_counts=True)
res = 0
for c in count:
p = float(c) / total
res -= p * np.log(p)
return res
19.3 K-Nearest Neighbor
Introduction
K-Nearest Neighbor is a supervised learning algorithm both for classification and regression. The principle is to find
the predefined number of training samples closest to the new point, and predict the label from these training samples1 .
For example, when a new point comes, the algorithm will follow these steps:
1. Calculate the Euclidean distance between the new point and all training data
2. Pick the top-K closest training data
3. For regression problem, take the average of the labels as the result; for classification problem, take the most
common label of these labels as the result.
Code
Below is the Numpy implementation of K-Nearest Neighbor function. Refer to code example for details.
1 https://scikit-learn.org/stable/modules/neighbors.html#nearest-neighbors-classification

def KNN(training_data, target, k, func):

"""
training_data: all training data point
target: new point
k: user-defined constant, number of closest training data
func: functions used to get the the target label
"""
# Step one: calculate the Euclidean distance between the new point and all
˓→training data
neighbors= []
for index, data in enumerate(training_data):
# distance between the target data and the current example from the data.
distance = euclidean_distance(data[:-1], target)
neighbors.append((distance, index))
# Step two: pick the top-K closest training data

sorted_neighbors = sorted(neighbors)
k_nearest = sorted_neighbors[:k]
k_nearest_labels = [training_data[i][1] for distance, i in k_nearest]
# Step three: For regression problem, take the average of the labels as the
˓→result;
# for classification problem, take the most common label of these
˓→labels as the result.
return k_nearest, func(k_nearest_labels)
19.4 Logistic Regression
19.5 Random Forests
Random Forest Classifier using ID3 Tree: code example
19.6 Boosting
19.7 Support Vector Machines
References
19.4. Logistic Regression 149


CHAPTER 20
Clustering Algorithms
20.1 Centroid
20.2 Density
20.3 Distribution
20.4 Hierarchical
20.5 K-Means
20.6 Mean shift
151
References
152 Chapter 20. Clustering Algorithms

CHAPTER 21
Regression Algorithms
• Ordinary Least Squares

• Polynomial
• Lasso
• Ridge
• Stepwise
21.1 Ordinary Least Squares
OLS is the method with which linear regression is performed. The square of the difference from the mean is taken for
every data point, and the summed loss function is to be minimized.
𝑛
∑︁
𝑙= (𝑦𝑖 − 𝑦¯)2
𝑖=1
21.2 Polynomial
Polynomial regression is a modification of linear regression where the existing features are mapped to a polynomial
form. The problem is still a linear regression problem, but the input vector is now mapped to a higher dimensional
vector which serves as a pseudo-input vector of sorts.
x = (𝑥0 , 𝑥1 ) → x’ = (𝑥0 , 𝑥20 , 𝑥1 , 𝑥21 , 𝑥0 𝑥1 )
153
21.3 Lasso
Lasso Regression tries to reduce the ordinary least squares error similar to vanilla regression, but adds an extra term.
The sum of the 𝐿1 norm for every data point multiplied by a hyperparameter 𝛼 is used. This reduces model complexity
and prevents overfitting.
𝑛
∑︁ 𝑝
∑︁
𝑙= (𝑦𝑖 − 𝑦˜)2 + 𝛼 |𝑤𝑗 |
𝑖=1 𝑗=1
21.4 Ridge
Ridge regression is similar to lasso regression, but the regularization term uses the 𝐿2 norm instead.
𝑛
∑︁ 𝑝
∑︁
𝑙= (𝑦𝑖 − 𝑦˜)2 + 𝛼 𝑤𝑗2
𝑖=1 𝑗=1
21.5 Stepwise
Stepwise regression or spline regression helps us fit a piece wise function to the data. It is usually used with linear
models, but it can be generalized to higher degrees as well. The regression equation takes the form of
𝑦 = 𝑎𝑥 + 𝑏(𝑥 − 𝑥
¯)𝐻𝛼 + 𝑐
where 𝐻𝛼 is the shifted Heaviside step function, having its discontinuity at 𝛼.
References
154 Chapter 21. Regression Algorithms

CHAPTER 22
Reinforcement Learning
Some links to have a brief about Reinforcemnt Learning.

For Practical Application of Reinforcement Learning:https://towardsdatascience.com/applications-of-reinforcement-
learning-in-real-world-1a94955bcd12
Reinforcement Learning concept:https://www.geeksforgeeks.org/what-is-reinforcement-learning/
Reinforcement Learning Algorithms:https://medium.com/@SmartLabAI/reinforcement-learning-algorithms-an-
intuitive-overview-904e2dff5bbc
Programmatic and a theoretical introduction to reinforcement learning:https://spinningup.openai.com/
155
156 Chapter 22. Reinforcement Learning

CHAPTER 23
Datasets
Public datasets in vision, nlp and more forked from caesar0301’s awesome datasets wiki.
• Agriculture
• Art
• Biology
• Chemistry/Materials Science
• Climate/Weather
• Complex Networks
• Computer Networks
• Data Challenges
• Earth Science
• Economics
• Education
• Energy
• Finance
• GIS
• Government
• Healthcare
• Image Processing
• Machine Learning
• Museums
157
• Music
• Natural Language
• Neuroscience
• Physics
• Psychology/Cognition
• Public Domains
• Search Engines
• Social Networks
• Social Sciences
• Software
• Sports
• Time Series
• Transportation
23.1 Agriculture
• U.S. Department of Agriculture’s PLANTS Database

• U.S. Department of Agriculture’s Nutrient Database
23.2 Art
• Google’s Quick Draw Sketch Dataset
23.3 Biology
• 1000 Genomes
• American Gut (Microbiome Project)
• Broad Bioimage Benchmark Collection (BBBC)
• Broad Cancer Cell Line Encyclopedia (CCLE)
• Cell Image Library
• Complete Genomics Public Data
• EBI ArrayExpress
• EBI Protein Data Bank in Europe
• Electron Microscopy Pilot Image Archive (EMPIAR)
• ENCODE project
• Ensembl Genomes
158 Chapter 23. Datasets

• Gene Expression Omnibus (GEO)

• Gene Ontology (GO)
• Global Biotic Interactions (GloBI)
• Harvard Medical School (HMS) LINCS Project
• Human Genome Diversity Project
• Human Microbiome Project (HMP)
• ICOS PSP Benchmark
• International HapMap Project
• Journal of Cell Biology DataViewer
• MIT Cancer Genomics Data
• NCBI Proteins
• NCBI Taxonomy
• NCI Genomic Data Commons
• NIH Microarray data or FTP (see FTP link on RAW)
• OpenSNP genotypes data
• Pathguid - Protein-Protein Interactions Catalog
• Protein Data Bank
• Psychiatric Genomics Consortium
• PubChem Project
• PubGene (now Coremine Medical)
• Sanger Catalogue of Somatic Mutations in Cancer (COSMIC)
• Sanger Genomics of Drug Sensitivity in Cancer Project (GDSC)
• Sequence Read Archive(SRA)
• Stanford Microarray Data
• Stowers Institute Original Data Repository
• Systems Science of Biological Dynamics (SSBD) Database
• The Cancer Genome Atlas (TCGA), available via Broad GDAC
• The Catalogue of Life
• The Personal Genome Project or PGP
• UCSC Public Data
• UniGene
• Universal Protein Resource (UnitProt)
23.3. Biology 159

23.4 Chemistry/Materials Science
• NIST Computational Chemistry Comparison and Benchmark Database - SRD 101

• Open Quantum Materials Database
• Citrination Public Datasets
• Khazana Project
23.5 Climate/Weather
• Actuaries Climate Index

• Australian Weather
• Aviation Weather Center - Consistent, timely and accurate weather information for the world airspace system
• Brazilian Weather - Historical data (In Portuguese)
• Canadian Meteorological Centre
• Climate Data from UEA (updated monthly)
• European Climate Assessment & Dataset
• Global Climate Data Since 1929
• NASA Global Imagery Browse Services
• NOAA Bering Sea Climate
• NOAA Climate Datasets
• NOAA Realtime Weather Models
• NOAA SURFRAD Meteorology and Radiation Datasets
• The World Bank Open Data Resources for Climate Change
• UEA Climatic Research Unit
• WorldClim - Global Climate Data
• WU Historical Weather Worldwide
23.6 Complex Networks
• AMiner Citation Network Dataset

• CrossRef DOI URLs
• DBLP Citation dataset
• DIMACS Road Networks Collection
• NBER Patent Citations
• Network Repository with Interactive Exploratory Analysis Tools
• NIST complex networks data collection
• Protein-protein interaction network

• PyPI and Maven Dependency Network

• Scopus Citation Database
• Small Network Data
• Stanford GraphBase (Steven Skiena)
• Stanford Large Network Dataset Collection
• Stanford Longitudinal Network Data Sources
• The Koblenz Network Collection
• The Laboratory for Web Algorithmics (UNIMI)
• The Nexus Network Repository
• UCI Network Data Repository
• UFL sparse matrix collection
• WSU Graph Database
23.7 Computer Networks
• 3.5B Web Pages from CommonCrawl 2012

• 53.5B Web clicks of 100K users in Indiana Univ.
• CAIDA Internet Datasets
• ClueWeb09 - 1B web pages
• ClueWeb12 - 733M web pages
• CommonCrawl Web Data over 7 years
• CRAWDAD Wireless datasets from Dartmouth Univ.
• Criteo click-through data
• OONI: Open Observatory of Network Interference - Internet censorship data
• Open Mobile Data by MobiPerf
• Rapid7 Sonar Internet Scans
• UCSD Network Telescope, IPv4 /8 net
23.8 Data Challenges
• Bruteforce Database
• Challenges in Machine Learning
• CrowdANALYTIX dataX
• D4D Challenge of Orange
• DrivenData Competitions for Social Good
• ICWSM Data Challenge (since 2009)
• Kaggle Competition Data
23.7. Computer Networks 161

• KDD Cup by Tencent 2012

• Localytics Data Visualization Challenge
• Netflix Prize
• Space Apps Challenge
• Telecom Italia Big Data Challenge
• TravisTorrent Dataset - MSR’2017 Mining Challenge
• Yelp Dataset Challenge
23.9 Earth Science
• AQUASTAT - Global water resources and uses

• BODC - marine data of ~22K vars
• Earth Models
• EOSDIS - NASA’s earth observing system data
• Integrated Marine Observing System (IMOS) - roughly 30TB of ocean measurements or on S3
• Marinexplore - Open Oceanographic Data
• Smithsonian Institution Global Volcano and Eruption Database
• USGS Earthquake Archives
23.10 Economics
• American Economic Association (AEA)

• EconData from UMD
• Economic Freedom of the World Data
• Historical MacroEconomc Statistics
• International Economics Database and various data tools
• International Trade Statistics
• Internet Product Code Database
• Joint External Debt Data Hub
• Jon Haveman International Trade Data Links
• OpenCorporates Database of Companies in the World
• Our World in Data
• SciencesPo World Trade Gravity Datasets
• The Atlas of Economic Complexity
• The Center for International Data
• The Observatory of Economic Complexity
• UN Commodity Trade Statistics

• UN Human Development Reports
23.11 Education
• College Scorecard Data

• Student Data from Free Code Camp
23.12 Energy
• AMPds
• BLUEd
• COMBED
• Dataport
• DRED
• ECO
• EIA
• HES - Household Electricity Study, UK
• HFED
• iAWE
• PLAID - the Plug Load Appliance Identification Dataset
• REDD
• Tracebase
• UK-DALE - UK Domestic Appliance-Level Electricity
• WHITED
23.13 Finance
• CBOE Futures Exchange

• Google Finance
• Google Trends
• NASDAQ
• NYSE Market Data (see FTP link on RAW)
• OANDA
• OSU Financial data
• Quandl
• St Louis Federal
• Yahoo Finance
23.11. Education 163

23.14 GIS
• ArcGIS Open Data portal

• Cambridge, MA, US, GIS data on GitHub
• Factual Global Location Data
• Geo Spatial Data from ASU
• Geo Wiki Project - Citizen-driven Environmental Monitoring
• GeoFabrik - OSM data extracted to a variety of formats and areas
• GeoNames Worldwide
• Global Administrative Areas Database (GADM)
• Homeland Infrastructure Foundation-Level Data
• Landsat 8 on AWS
• List of all countries in all languages
• National Weather Service GIS Data Portal
• Natural Earth - vectors and rasters of the world
• OpenAddresses
• OpenStreetMap (OSM)
• Pleiades - Gazetteer and graph of ancient places
• Reverse Geocoder using OSM data & additional high-resolution data files
• TIGER/Line - U.S. boundaries and roads
• TwoFishes - Foursquare’s coarse geocoder
• TZ Timezones shapfiles
• UN Environmental Data
• World boundaries from the U.S. Department of State
• World countries in multiple formats
23.15 Government
• A list of cities and countries contributed by community

• Open Data for Africa
• OpenDataSoft’s list of 1,600 open data
23.16 Healthcare
• EHDP Large Health Data Sets

• Gapminder World demographic databases
• Medicare Coverage Database (MCD), U.S.

• Medicare Data Engine of medicare.gov Data

• Medicare Data File
• MeSH, the vocabulary thesaurus used for indexing articles for PubMed
• Number of Ebola Cases and Deaths in Affected Countries (2014)
• Open-ODS (structure of the UK NHS)
• OpenPaymentsData, Healthcare financial relationship data
• The Cancer Genome Atlas project (TCGA) and BigQuery table
• World Health Organization Global Health Observatory
23.17 Image Processing
• 10k US Adult Faces Database

• 2GB of Photos of Cats or Archive version
• Adience Unfiltered faces for gender and age classification
• Affective Image Classification
• Animals with attributes
• Caltech Pedestrian Detection Benchmark
• Chars74K dataset, Character Recognition in Natural Images (both English and Kannada are available)
• Face Recognition Benchmark
• GDXray: X-ray images for X-ray testing and Computer Vision
• ImageNet (in WordNet hierarchy)
• Indoor Scene Recognition
• International Affective Picture System, UFL
• Massive Visual Memory Stimuli, MIT
• MNIST database of handwritten digits, near 1 million examples
• Several Shape-from-Silhouette Datasets
• Stanford Dogs Dataset
• SUN database, MIT
• The Action Similarity Labeling (ASLAN) Challenge
• The Oxford-IIIT Pet Dataset
• Violent-Flows - Crowd Violence Non-violence Database and benchmark
• Visual genome
• YouTube Faces Database
23.17. Image Processing 165

23.18 Machine Learning
• Context-aware data sets from five domains

• Delve Datasets for classification and regression (Univ. of Toronto)
• Discogs Monthly Data
• eBay Online Auctions (2012)
• IMDb Database
• Keel Repository for classification, regression and time series
• Labeled Faces in the Wild (LFW)
• Lending Club Loan Data
• Machine Learning Data Set Repository
• Million Song Dataset
• More Song Datasets
• MovieLens Data Sets
• New Yorker caption contest ratings
• RDataMining - “R and Data Mining” ebook data
• Registered Meteorites on Earth
• Restaurants Health Score Data in San Francisco
• UCI Machine Learning Repository
• Yahoo! Ratings and Classification Data
• Youtube 8m
23.19 Museums
• Canada Science and Technology Museums Corporation’s Open Data

• Cooper-Hewitt’s Collection Database
• Minneapolis Institute of Arts metadata
• Natural History Museum (London) Data Portal
• Rijksmuseum Historical Art Collection
• Tate Collection metadata
• The Getty vocabularies
23.20 Music
• Nottingham Folk Songs

• Bach 10

23.21 Natural Language
• Automatic Keyphrase Extracttion

• Blogger Corpus
• CLiPS Stylometry Investigation Corpus
• ClueWeb09 FACC
• ClueWeb12 FACC
• DBpedia - 4.58M things with 583M facts
• Flickr Personal Taxonomies
• Freebase.com of people, places, and things
• Google Books Ngrams (2.2TB)
• Google MC-AFP, generated based on the public available Gigaword dataset using Paragraph Vectors
• Google Web 5gram (1TB, 2006)
• Gutenberg eBooks List
• Hansards text chunks of Canadian Parliament
• Machine Comprehension Test (MCTest) of text from Microsoft Research
• Machine Translation of European languages
• Microsoft MAchine Reading COmprehension Dataset (or MS MARCO)
• Multi-Domain Sentiment Dataset (version 2.0)
• Open Multilingual Wordnet
• Personae Corpus
• SaudiNewsNet Collection of Saudi Newspaper Articles (Arabic, 30K articles)
• SMS Spam Collection in English
• Universal Dependencies
• USENET postings corpus of 2005~2011
• Webhose - News/Blogs in multiple languages
• Wikidata - Wikipedia databases
• Wikipedia Links data - 40 Million Entities in Context
• WordNet databases and tools
23.22 Neuroscience
• Allen Institute Datasets

• Brain Catalogue
• Brainomics
• CodeNeuro Datasets
• Collaborative Research in Computational Neuroscience (CRCNS)
23.21. Natural Language 167

• FCP-INDI
• Human Connectome Project
• NDAR
• NeuroData
• Neuroelectro
• NIMH Data Archive
• OASIS
• OpenfMRI
• Study Forrest
23.23 Physics
• CERN Open Data Portal

• Crystallography Open Database
• NASA Exoplanet Archive
• NSSDC (NASA) data of 550 space spacecraft
• Sloan Digital Sky Survey (SDSS) - Mapping the Universe
23.24 Psychology/Cognition
• OSU Cognitive Modeling Repository Datasets
23.25 Public Domains
• Amazon
• Archive-it from Internet Archive
• Archive.org Datasets
• CMU JASA data archive
• CMU StatLab collections
• Data.World
• Data360
• Datamob.org
• Google
• Infochimps
• KDNuggets Data Collections
• Microsoft Azure Data Market Free DataSets
• Microsoft Data Science for Research

• Numbray
• Open Library Data Dumps
• Reddit Datasets
• RevolutionAnalytics Collection
• Sample R data sets
• Stats4Stem R data sets
• StatSci.org
• The Washington Post List
• UCLA SOCR data collection
• UFO Reports
• Wikileaks 911 pager intercepts
• Yahoo Webscope
23.26 Search Engines
• Academic Torrents of data sharing from UMB

• Datahub.io
• DataMarket (Qlik)
• Harvard Dataverse Network of scientific data
• ICPSR (UMICH)
• Institute of Education Sciences
• National Technical Reports Library
• Open Data Certificates (beta)
• OpenDataNetwork - A search engine of all Socrata powered data portals
• Statista.com - statistics and Studies
• Zenodo - An open dependable home for the long-tail of science
23.27 Social Networks
• 72 hours #gamergate Twitter Scrape

• Ancestry.com Forum Dataset over 10 years
• Cheng-Caverlee-Lee September 2009 - January 2010 Twitter Scrape
• CMU Enron Email of 150 users
• EDRM Enron EMail of 151 users, hosted on S3
• Facebook Data Scrape (2005)
• Facebook Social Networks from LAW (since 2007)
• Foursquare from UMN/Sarwat (2013)
23.26. Search Engines 169

• GitHub Collaboration Archive

• Google Scholar citation relations
• High-Resolution Contact Networks from Wearable Sensors
• Mobile Social Networks from UMASS
• Network Twitter Data
• Reddit Comments
• Skytrax’ Air Travel Reviews Dataset
• Social Twitter Data
• SourceForge.net Research Data
• Twitter Data for Online Reputation Management
• Twitter Data for Sentiment Analysis
• Twitter Graph of entire Twitter site
• Twitter Scrape Calufa May 2011
• UNIMI/LAW Social Network Datasets
• Yahoo! Graph and Social Data
• Youtube Video Social Graph in 2007,2008
23.28 Social Sciences
• ACLED (Armed Conflict Location & Event Data Project)

• Canadian Legal Information Institute
• Center for Systemic Peace Datasets - Conflict Trends, Polities, State Fragility, etc
• Correlates of War Project
• Cryptome Conspiracy Theory Items
• Datacards
• European Social Survey
• FBI Hate Crime 2013 - aggregated data
• Fragile States Index
• GDELT Global Events Database
• General Social Survey (GSS) since 1972
• German Social Survey
• Global Religious Futures Project
• Humanitarian Data Exchange
• INFORM Index for Risk Management
• Institute for Demographic Studies
• International Networks Archive

• International Social Survey Program ISSP

• International Studies Compendium Project
• James McGuire Cross National Data
• MacroData Guide by Norsk samfunnsvitenskapelig datatjeneste
• Minnesota Population Center
• MIT Reality Mining Dataset
• Notre Dame Global Adaptation Index (NG-DAIN)
• Open Crime and Policing Data in England, Wales and Northern Ireland
• Paul Hensel General International Data Page
• PewResearch Internet Survey Project
• PewResearch Society Data Collection
• Political Polarity Data
• StackExchange Data Explorer
• Terrorism Research and Analysis Consortium
• Texas Inmates Executed Since 1984
• Titanic Survival Data Set or on Kaggle
• UCB’s Archive of Social Science Data (D-Lab)
• UCLA Social Sciences Data Archive
• UN Civil Society Database
• Universities Worldwide
• UPJOHN for Labor Employment Research
• Uppsala Conflict Data Program
• World Bank Open Data
• WorldPop project - Worldwide human population distributions
23.29 Software
• FLOSSmole data about free, libre, and open source software development
23.30 Sports
• Basketball (NBA/NCAA/Euro) Player Database and Statistics

• Betfair Historical Exchange Data
• Cricsheet Matches (cricket)
• Ergast Formula 1, from 1950 up to date (API)
• Football/Soccer resources (data and APIs)
• Lahman’s Baseball Database
23.29. Software 171

• Pinhooker: Thoroughbred Bloodstock Sale Data

• Retrosheet Baseball Statistics
• Tennis database of rankings, results, and stats for ATP, WTA, Grand Slams and Match Charting Project
23.31 Time Series
• Databanks International Cross National Time Series Data Archive

• Hard Drive Failure Rates
• Heart Rate Time Series from MIT
• Time Series Data Library (TSDL) from MU
• UC Riverside Time Series Dataset
23.32 Transportation
• Airlines OD Data 1987-2008

• Bay Area Bike Share Data
• Bike Share Systems (BSS) collection
• GeoLife GPS Trajectory from Microsoft Research
• German train system by Deutsche Bahn
• Hubway Million Rides in MA
• Marine Traffic - ship tracks, port calls and more
• Montreal BIXI Bike Share
• NYC Taxi Trip Data 2009-
• NYC Taxi Trip Data 2013 (FOIA/FOILed)
• NYC Uber trip data April 2014 to September 2014
• Open Traffic collection
• OpenFlights - airport, airline and route data
• Philadelphia Bike Share Stations (JSON)
• Plane Crash Database, since 1920
• RITA Airline On-Time Performance data
• RITA/BTS transport data collection (TranStat)
• Toronto Bike Share Stations (XML file)
• Transport for London (TFL)
• Travel Tracker Survey (TTS) for Chicago
• U.S. Bureau of Transportation Statistics (BTS)
• U.S. Domestic Flights 1990 to 2009
• U.S. Freight Analysis Framework since 2007

CHAPTER 24
Libraries
Machine learning libraries and frameworks forked from josephmisti’s awesome machine learning.
• APL
• C
• C++
• Common Lisp
• Clojure
• Elixir
• Erlang
• Go
• Haskell
• Java
• Javascript
• Julia
• Lua
• Matlab
• .NET
• Objective C
• OCaml
• PHP
• Python
173
• Ruby
• Rust
• R
• SAS
• Scala
• Swift
24.1 APL
General-Purpose Machine Learning
• naive-apl - Naive Bayesian Classifier implementation in APL
24.2 C
• Darknet - Darknet is an open source neural network framework written in C and CUDA. It is fast, easy to install,
and supports CPU and GPU computation.
• Recommender - A C library for product recommendations/suggestions using collaborative filtering (CF).
• Hybrid Recommender System - A hybrid recomender system based upon scikit-learn algorithms.
Computer Vision
• CCV - C-based/Cached/Core Computer Vision Library, A Modern Computer Vision Library

• VLFeat - VLFeat is an open and portable library of computer vision algorithms, which has Matlab toolbox
Speech Recognition
• HTK -The Hidden Markov Model Toolkit. HTK is a portable toolkit for building and manipulating hidden
Markov models.
24.3 C++
Computer Vision
• DLib - DLib has C++ and Python interfaces for face detection and training general object detectors.
• EBLearn - Eblearn is an object-oriented C++ library that implements various machine learning models
• OpenCV - OpenCV has C++, C, Python, Java and MATLAB interfaces and supports Windows, Linux, Android
and Mac OS.
174 Chapter 24. Libraries

• VIGRA - VIGRA is a generic cross-platform C++ computer vision and machine learning library for volumes of
arbitrary dimensionality with Python bindings.
• BanditLib - A simple Multi-armed Bandit library.

• Caffe - A deep learning framework developed with cleanliness, readability, and speed in mind. [DEEP LEARN-
ING]
• CNTK by Microsoft Research, is a unified deep-learning toolkit that describes neural networks as a series of
computational steps via a directed graph.
• CUDA - This is a fast C++/CUDA implementation of convolutional [DEEP LEARNING]
• CXXNET - Yet another deep learning framework with less than 1000 lines core code [DEEP LEARNING]
• DeepDetect - A machine learning API and server written in C++11. It makes state of the art machine learning
easy to work with and integrate into existing applications.
• Disrtibuted Machine learning Tool Kit (DMTK) Word Embedding.
• DLib - A suite of ML tools designed to be easy to imbed in other applications
• DSSTNE - A software library created by Amazon for training and deploying deep neural networks using GPUs
which emphasizes speed and scale over experimental flexibility.
• DyNet - A dynamic neural network library working well with networks that have dynamic structures that change
for every training instance. Written in C++ with bindings in Python.
• encog-cpp
• Fido - A highly-modular C++ machine learning library for embedded electronics and robotics.
• igraph - General purpose graph library
• Intel(R) DAAL - A high performance software library developed by Intel and optimized for Intel’s architectures.
Library provides algorithmic building blocks for all stages of data analytics and allows to process data in batch,
online and distributed modes.
• LightGBM framework based on decision tree algorithms, used for ranking, classification and many other ma-
chine learning tasks.
• MLDB - The Machine Learning Database is a database designed for machine learning. Send it commands over
a RESTful API to store data, explore it using SQL, then train machine learning models and expose them as
APIs.
• mlpack - A scalable C++ machine learning library
• ROOT - A modular scientific software framework. It provides all the functionalities needed to deal with big data
processing, statistical analysis, visualization and storage.
• shark - A fast, modular, feature-rich open-source C++ machine learning library.
• Shogun - The Shogun Machine Learning Toolbox
• sofia-ml - Suite of fast incremental algorithms.
• Stan - A probabilistic programming language implementing full Bayesian statistical inference with Hamiltonian
Monte Carlo sampling
• Timbl - A software package/C++ library implementing several memory-based learning algorithms, among which
IB1-IG, an implementation of k-nearest neighbor classification, and IGTree, a decision-tree approximation of
IB1-IG. Commonly used for NLP.
24.3. C++ 175

• Vowpal Wabbit (VW) - A fast out-of-core learning system.

• Warp-CTC, on both CPU and GPU.
• XGBoost - A parallelized optimized general purpose gradient boosting library.
Natural Language Processing
• BLLIP Parser
• colibri-core - C++ library, command line tools, and Python binding for extracting and working with basic lin-
guistic constructions such as n-grams and skipgrams in a quick and memory-efficient way.
• CRF++ for segmenting/labeling sequential data & other Natural Language Processing tasks.
• CRFsuite for labeling sequential data.
• frog - Memory-based NLP suite developed for Dutch: PoS tagger, lemmatiser, dependency parser, NER, shallow
parser, morphological analyzer.
• libfolia](https://github.com/LanguageMachines/libfolia) - C++ library for the [FoLiA format
• MeTA](https://github.com/meta-toolkit/meta) - [MeTA : ModErn Text Analysis is a C++ Data Sciences Toolkit
that facilitates mining big text data.
• MIT Information Extraction Toolkit - C, C++, and Python tools for named entity recognition and relation ex-
traction
• ucto - Unicode-aware regular-expression based tokenizer for various languages. Tool and C++ library. Supports
FoLiA format.
Speech Recognition
• Kaldi - Kaldi is a toolkit for speech recognition written in C++ and licensed under the Apache License v2.0.
Kaldi is intended for use by speech recognition researchers.
Sequence Analysis
• ToPS - This is an objected-oriented framework that facilitates the integration of probabilistic models for se-
quences over a user defined alphabet.
Gesture Detection
• grt - The Gesture Recognition Toolkit. GRT is a cross-platform, open-source, C++ machine learning library
designed for real-time gesture recognition.
24.4 Common Lisp
• mgl, Gaussian Processes

• mgl-gpr - Evolutionary algorithms
• cl-libsvm - Wrapper for the libsvm support vector machine library

24.5 Clojure
• Clojure-openNLP - Natural Language Processing in Clojure (opennlp)

• Infections-clj - Rails-like inflection library for Clojure and ClojureScript
• Touchstone - Clojure A/B testing library

• Clojush - The Push programming language and the PushGP genetic programming system implemented in Clo-
jure
• Infer - Inference and machine learning in clojure
• Clj-ML - A machine learning library for Clojure built on top of Weka and friends
• DL4CLJ - Clojure wrapper for Deeplearning4j
• Encog
• Fungp - A genetic programming library for Clojure
• Statistiker - Basic Machine Learning algorithms in Clojure.
• clortex - General Machine Learning library using Numenta’s Cortical Learning Algorithm
• comportex - Functionally composable Machine Learning library using Numenta’s Cortical Learning Algorithm
• cortex - Neural networks, regression and feature learning in Clojure.
• lambda-ml - Simple, concise implementations of machine learning techniques and utilities in Clojure.
Data Analysis / Data Visualization
• Incanter - Incanter is a Clojure-based, R-like platform for statistical computing and graphics.
• PigPen - Map-Reduce for Clojure.
• Envision - Clojure Data Visualisation library, based on Statistiker and D3
24.6 Elixir
• Simple Bayes - A Simple Bayes / Naive Bayes implementation in Elixir.
• Stemmer stemming implementation in Elixir.
24.5. Clojure 177

24.7 Erlang
• Disco - Map Reduce in Erlang
24.8 Go
• go-porterstemmer - A native Go clean room implementation of the Porter Stemming algorithm.

• paicehusk - Golang implementation of the Paice/Husk Stemming Algorithm.
• snowball - Snowball Stemmer for Go.
• go-ngram - In-memory n-gram index with compression.
• gago - Multi-population, flexible, parallel genetic algorithm.

• Go Learn - Machine Learning for Go
• go-pr - Pattern recognition package in Go lang.
• go-ml - Linear / Logistic regression, Neural Networks, Collaborative Filtering and Gaussian Multivariate Dis-
tribution
• bayesian - Naive Bayesian Classification for Golang.
• go-galib - Genetic Algorithms library written in Go / golang
• Cloudforest - Ensembles of decision trees in go/golang.
• gobrain - Neural Networks written in go
• GoNN - GoNN is an implementation of Neural Network in Go Language, which includes BPNN, RBF, PCN
• MXNet - Lightweight, Portable, Flexible Distributed/Mobile Deep Learning with Dynamic, Mutation-aware
Dataflow Dep Scheduler; for Python, R, Julia, Go, Javascript and more.
• go-mxnet-predictor - Go binding for MXNet c_predict_api to do inference with pre-trained model
• go-graph - Graph library for Go/golang language.

• SVGo - The Go Language library for SVG generation
• RF - Random forests implementation in Go

24.9 Haskell
• haskell-ml - Haskell implementations of various ML algorithms.

• HLearn - a suite of libraries for interpreting machine learning models according to their algebraic structure.
• hnn - Haskell Neural Network library.
• hopfield-networks - Hopfield Networks for unsupervised learning in Haskell.
• caffegraph - A DSL for deep neural networks
• LambdaNet - Configurable Neural Networks in Haskell
24.10 Java
• Cortical.io as quickly and intuitively as the brain.

• CoreNLP - Stanford CoreNLP provides a set of natural language analysis tools which can take raw English
language text input and give the base forms of words
• Stanford Parser - A natural language parser is a program that works out the grammatical structure of sentences
• Stanford POS Tagger - A Part-Of-Speech Tagger (POS Tagger
• Stanford Name Entity Recognizer - Stanford NER is a Java implementation of a Named Entity Recognizer.
• Stanford Word Segmenter - Tokenization of raw text is a standard pre-processing step for many NLP tasks.
• Tregex, Tsurgeon and Semgrex.
• Stanford Phrasal: A Phrase-Based Translation System
• Stanford English Tokenizer - Stanford Phrasal is a state-of-the-art statistical phrase-based machine translation
system, written in Java.
• Stanford Tokens Regex - A tokenizer divides text into a sequence of tokens, which roughly correspond to
“words”
• Stanford Temporal Tagger - SUTime is a library for recognizing and normalizing time expressions.
• Stanford SPIED - Learning entities from unlabeled text starting with seed sets using patterns in an iterative
fashion
• Stanford Topic Modeling Toolbox - Topic modeling tools to social scientists and others who wish to perform
analysis on datasets
• Twitter Text Java - A Java implementation of Twitter’s text processing library
• MALLET - A Java-based package for statistical natural language processing, document classification, cluster-
ing, topic modeling, information extraction, and other machine learning applications to text.
• OpenNLP - a machine learning based toolkit for the processing of natural language text.
• LingPipe - A tool kit for processing text using computational linguistics.
• ClearTK components in Java and is built on top of Apache UIMA.
24.9. Haskell 179

• Apache cTAKES is an open-source natural language processing system for information extraction from elec-
tronic medical record clinical free-text.
• ClearNLP - The ClearNLP project provides software and resources for natural language processing. The project
started at the Center for Computational Language and EducAtion Research, and is currently developed by the
Center for Language and Information Research at Emory University. This project is under the Apache 2 license.
• CogcompNLP developed in the University of Illinois’ Cognitive Computation Group, for example illinois-core-
utilities which provides a set of NLP-friendly data structures and a number of NLP-related utilities that support
writing NLP applications, running experiments, etc, illinois-edison a library for feature extraction from illinois-
core-utilities data structures and many other packages.
• aerosolve - A machine learning library by Airbnb designed from the ground up to be human friendly.
• Datumbox - Machine Learning framework for rapid development of Machine Learning and Statistical applica-
tions
• ELKI
• Encog - An advanced neural network and machine learning framework. Encog contains classes to create a wide
variety of networks, as well as support classes to normalize and process data for these neural networks. Encog
trains using multithreaded resilient propagation. Encog can also make use of a GPU to further speed processing
time. A GUI based workbench is also provided to help model and train neural networks.
• FlinkML in Apache Flink - Distributed machine learning library in Flink
• H2O - ML engine that supports distributed learning on Hadoop, Spark or your laptop via APIs in R, Python,
Scala, REST/JSON.
• htm.java - General Machine Learning library using Numenta’s Cortical Learning Algorithm
• java-deeplearning - Distributed Deep Learning Platform for Java, Clojure,Scala
• Mahout - Distributed machine learning
• Meka.
• MLlib in Apache Spark - Distributed machine learning library in Spark
• Hydrosphere Mist - a service for deployment Apache Spark MLLib machine learning models as realtime, batch
or reactive web services.
• Neuroph - Neuroph is lightweight Java neural network framework
• ORYX - Lambda Architecture Framework using Apache Spark and Apache Kafka with a specialization for
real-time large-scale machine learning.
• Samoa SAMOA is a framework that includes distributed machine learning for data streams with an interface to
plug-in different stream processing platforms.
• RankLib - RankLib is a library of learning to rank algorithms
• rapaio - statistics, data mining and machine learning toolbox in Java
• RapidMiner - RapidMiner integration into Java code
• Stanford Classifier - A classifier is a machine learning tool that will take data items and place them into one of
k classes.
• SmileMiner - Statistical Machine Intelligence & Learning Engine
• SystemML language.

• WalnutiQ - object oriented model of the human brain

• Weka - Weka is a collection of machine learning algorithms for data mining tasks
• LBJava - Learning Based Java is a modeling language for the rapid development of software systems, offers
a convenient, declarative syntax for classifier and constraint definition directly in terms of the objects in the
programmer’s application.
Speech Recognition
• CMU Sphinx - Open Source Toolkit For Speech Recognition purely based on Java speech recognition library.
• Flink - Open source platform for distributed stream and batch data processing.
• Hadoop - Hadoop/HDFS
• Spark - Spark is a fast and general engine for large-scale data processing.
• Storm - Storm is a distributed realtime computation system.
• Impala - Real-time Query for Hadoop
• DataMelt - Mathematics software for numeric computation, statistics, symbolic calculations, data analysis and
data visualization.
• Dr. Michael Thomas Flanagan’s Java Scientific Library
Deep Learning
• Deeplearning4j - Scalable deep learning for industry with parallel GPUs
24.11 Javascript
• Twitter-text - A JavaScript implementation of Twitter’s text processing library

• NLP.js - NLP utilities in javascript and coffeescript
• natural - General natural language facilities for node
• Knwl.js - A Natural Language Processor in JS
• Retext - Extensible system for analyzing and manipulating natural language
• TextProcessing - Sentiment analysis, stemming and lemmatization, part-of-speech tagging and chunking, phrase
extraction and named entity recognition.
• NLP Compromise - Natural Language processing in the browser
24.11. Javascript 181

• D3.js
• High Charts
• NVD3.js
• dc.js
• chartjs
• dimple
• amCharts
• D3xter - Straight forward plotting built on D3
• statkit - Statistics kit for JavaScript
• datakit - A lightweight framework for data analysis in JavaScript
• science.js - Scientific and statistical computing in JavaScript.
• Z3d - Easily make interactive 3d plots built on Three.js
• Sigma.js - JavaScript library dedicated to graph drawing.
• C3.js- customizable library based on D3.js for easy chart drawing.
• Datamaps- Customizable SVG map/geo visualizations using D3.js.
• ZingChart- library written on Vanilla JS for big data visualization.
• cheminfo - Platform for data visualization and analysis, using the visualizer project.
• Convnet.js - ConvNetJS is a Javascript library for training Deep Learning models[DEEP LEARNING]
• Clusterfck - Agglomerative hierarchical clustering implemented in Javascript for Node.js and the browser
• Clustering.js - Clustering algorithms implemented in Javascript for Node.js and the browser
• Decision Trees - NodeJS Implementation of Decision Tree using ID3 Algorithm
• DN2A - Digital Neural Networks Architecture
• figue - K-means, fuzzy c-means and agglomerative clustering
• Node-fann bindings for Node.js
• Kmeans.js - Simple Javascript implementation of the k-means algorithm, for node.js and the browser
• LDA.js - LDA topic modeling for node.js
• Learning.js - Javascript implementation of logistic regression/c4.5 decision tree
• Machine Learning - Machine learning library for Node.js
• machineJS - Automated machine learning, data formatting, ensembling, and hyperparameter optimization for
competitions and exploration- just give it a .csv file!
• mil-tokyo - List of several machine learning libraries
• Node-SVM - Support Vector Machine for nodejs
• Brain - Neural networks in JavaScript [Deprecated]

• Bayesian-Bandit - Bayesian bandit implementation for Node and the browser.

• Synaptic - Architecture-free neural network library for node.js and the browser
• kNear - JavaScript implementation of the k nearest neighbors algorithm for supervised learning
• NeuralN - C++ Neural Network library for Node.js. It has advantage on large dataset and multi-threaded train-
ing.
• kalman - Kalman filter for Javascript.
• shaman - node.js library with support for both simple and multiple linear regression.
• ml.js - Machine learning and numerical analysis tools for Node.js and the Browser!
• Pavlov.js - Reinforcement learning using Markov Decision Processes
Misc
• sylvester - Vector and Matrix math for JavaScript.

• simple-statistics as well as in node.js.
• regression-js - A javascript library containing a collection of least squares fitting methods for finding a trend in
a set of data.
• Lyric - Linear Regression library.
• GreatCircle - Library for calculating great circle distance.
24.12 Julia
• MachineLearning - Julia Machine Learning library

• MLBase - A set of functions to support the development of machine learning algorithms
• PGM - A Julia framework for probabilistic graphical models.
• DA - Julia package for Regularized Discriminant Analysis
• Regression
• Local Regression - Local regression, so smooooth!
• Naive Bayes - Simple Naive Bayes implementation in Julia
• Mixed Models mixed-effects models
• Simple MCMC - basic mcmc sampler implemented in Julia
• Distance - Julia module for Distance evaluation
• Decision Tree - Decision Tree Classifier and Regressor
• Neural - A neural network in Julia
• MCMC - MCMC tools for Julia
• Mamba for Bayesian analysis in Julia
24.12. Julia 183

• GLM - Generalized linear models in Julia

• Online Learning
• GLMNet - Julia wrapper for fitting Lasso/ElasticNet GLM models using glmnet
• Clustering - Basic functions for clustering data: k-means, dp-means, etc.
• SVM - SVM’s for Julia
• Kernal Density - Kernel density estimators for julia
• Dimensionality Reduction - Methods for dimensionality reduction
• NMF - A Julia package for non-negative matrix factorization
• ANN - Julia artificial neural networks
• Mocha - Deep Learning framework for Julia inspired by Caffe
• XGBoost - eXtreme Gradient Boosting Package in Julia
• ManifoldLearning - A Julia package for manifold learning and nonlinear dimensionality reduction
• Merlin - Flexible Deep Learning Framework in Julia
• ROCAnalysis - Receiver Operating Characteristics and functions for evaluation probabilistic binary classifiers
• GaussianMixtures - Large scale Gaussian Mixture Models
• ScikitLearn - Julia implementation of the scikit-learn API
• Knet - Koç University Deep Learning Framework
• Topic Models - TopicModels for Julia

• Text Analysis - Julia package for text analysis
• Graph Layout - Graph layout algorithms in pure Julia

• Data Frames Meta - Metaprogramming tools for DataFrames
• Julia Data - library for working with tabular data in Julia
• Data Read - Read files from Stata, SAS, and SPSS
• Hypothesis Tests - Hypothesis tests for Julia
• Gadfly - Crafty statistical graphics for Julia.
• Stats - Statistical tests for Julia
• RDataSets - Julia package for loading many of the data sets available in R
• DataFrames - library for working with tabular data in Julia
• Distributions - A Julia package for probability distributions and associated functions.
• Data Arrays - Data structures that allow missing values

• Time Series - Time series toolkit for Julia

• Sampling - Basic sampling algorithms for Julia
Misc Stuff / Presentations
• DSP.
• JuliaCon Presentations - Presentations for JuliaCon
• SignalProcessing - Signal Processing tools for Julia
• Images - An image library for Julia
24.13 Lua
• Torch7
• cephes - Cephes mathematical functions library, wrapped for Torch. Provides and wraps the 180+ special
mathematical functions from the Cephes mathematical library, developed by Stephen L. Moshier. It is used,
among many other places, at the heart of SciPy.
• autograd - Autograd automatically differentiates native Torch code. Inspired by the original Python version.
• graph - Graph package for Torch
• randomkit - Numpy’s randomkit, wrapped for Torch
• signal - A signal processing toolbox for Torch-7. FFT, DCT, Hilbert, cepstrums, stft
• nn - Neural Network package for Torch
• torchnet - framework for torch which provides a set of abstractions aiming at encouraging code re-use as well
as encouraging modular programming
• nngraph - This package provides graphical computation for nn library in Torch7.
• nnx - A completely unstable and experimental package that extends Torch’s builtin nn library
• rnn - A Recurrent Neural Network library that extends Torch’s nn. RNNs, LSTMs, GRUs, BRNNs, BLSTMs,
etc.
• dpnn - Many useful features that aren’t part of the main nn package.
• dp - A deep learning library designed for streamlining research and development using the Torch7 distribution.
It emphasizes flexibility through the elegant use of object-oriented design patterns.
• optim - An optimization library for Torch. SGD, Adagrad, Conjugate-Gradient, LBFGS, RProp and more.
• unsup.
• manifold - A package to manipulate manifolds
• svm - Torch-SVM library
• lbfgs - FFI Wrapper for liblbfgs
• vowpalwabbit - An old vowpalwabbit interface to torch.
• OpenGM - OpenGM is a C++ library for graphical modeling, and inference. The Lua bindings provide a simple
way of describing graphs, from Lua, and then optimizing them with OpenGM.
24.13. Lua 185

• sphagetti module for torch7 by @MichaelMathieu

• LuaSHKit - A lua wrapper around the Locality sensitive hashing library SHKit
• kernel smoothing - KNN, kernel-weighted average, local linear regression smoothers
• cutorch - Torch CUDA Implementation
• cunn - Torch CUDA Neural Network Implementation
• imgraph - An image/graph library for Torch. This package provides routines to construct graphs on images,
segment them, build trees out of them, and convert them back to images.
• videograph - A video/graph library for Torch. This package provides routines to construct graphs on videos,
segment them, build trees out of them, and convert them back to videos.
• saliency - code and tools around integral images. A library for finding interest points based on fast integral
histograms.
• stitch - allows us to use hugin to stitch images and apply same stitching to a video sequence
• sfm - A bundle adjustment/structure from motion package
• fex - A package for feature extraction in Torch. Provides SIFT and dSIFT modules.
• OverFeat - A state-of-the-art generic dense feature extractor
• Numeric Lua
• Lunatic Python
• SciLua
• Lua - Numerical Algorithms
• Lunum
Demos and Scripts
• Core torch7 demos repository. * linear-regression, logistic-regression * face detector (training and detection
as separate demos) * mst-based-segmenter * train-a-digit-classifier * train-autoencoder * optical flow demo *
train-on-housenumbers * train-on-cifar * tracking with deep nets * kinect demo * filter-bank visualization *
saliency-networks
• Training a Convnet for the Galaxy-Zoo Kaggle challenge(CUDA demo)
• Music Tagging - Music Tagging scripts for torch7
• torch-datasets - Scripts to load several popular datasets including: * BSR 500 * CIFAR-10 * COIL * Street
View House Numbers * MNIST * NORB
• Atari2600 - Scripts to generate a dataset with static frames from the Arcade Learning Environment
24.14 Matlab
Computer Vision
• Contourlets - MATLAB source code that implements the contourlet transform and its utility functions.
• Shearlets - MATLAB code for shearlet transform

• Curvelets - The Curvelet transform is a higher dimensional generalization of the Wavelet transform designed to
represent images at different scales and different angles.
• Bandlets - MATLAB code for bandlet transform
• mexopencv - Collection and a development kit of MATLAB mex functions for OpenCV library
• NLP - An NLP library for Matlab
• Training a deep autoencoder or a classifier on MNIST

• Convolutional-Recursive Deep Learning for 3D Object Classification - Convolutional-Recursive Deep Learning
for 3D Object Classification[DEEP LEARNING]
• t-Distributed Stochastic Neighbor Embedding technique for dimensionality reduction that is particularly well
suited for the visualization of high-dimensional datasets.
• Spider - The spider is intended to be a complete object orientated environment for machine learning in Matlab.
• LibSVM - A Library for Support Vector Machines
• LibLinear - A Library for Large Linear Classification
• Machine Learning Module - Class on machine w/ PDF,lectures,code
• Caffe - A deep learning framework developed with cleanliness, readability, and speed in mind.
• Pattern Recognition Toolbox - A complete object-oriented environment for machine learning in Matlab.
• Pattern Recognition and Machine Learning - This package contains the matlab implementation of the algorithms
described in the book Pattern Recognition and Machine Learning by C. Bishop.
• Optunity - A library dedicated to automated hyperparameter optimization with a simple, lightweight API to
facilitate drop-in replacement of grid search. Optunity is written in Python but interfaces seamlessly with MAT-
LAB.
• matlab_gbl - MatlabBGL is a Matlab package for working with graphs.

• gamic - Efficient pure-Matlab implementations of graph algorithms to complement MatlabBGL’s mex functions.
24.15 .NET
Computer Vision
• OpenCVDotNet - A wrapper for the OpenCV project to be used with .NET applications.
• Emgu CV - Cross platform wrapper of OpenCV which can be compiled in Mono to e run on Windows, Linus,
Mac OS X, iOS, and Android.
• AForge.NET - Open source C# framework for developers and researchers in the fields of Computer Vision and
Artificial Intelligence. Development has now shifted to GitHub.
24.15. .NET 187

• Accord.NET - Together with AForge.NET, this library can provide image processing and computer vision al-
gorithms to Windows, Windows RT and Windows Phone. Some components are also available for Java and
Android.
• Stanford.NLP for .NET - A full port of Stanford NLP packages to .NET and also available precompiled as a
NuGet package.
• Accord-Framework -The Accord.NET Framework is a complete framework for building machine learning, com-
puter vision, computer audition, signal processing and statistical applications.
• Accord.MachineLearning - Support Vector Machines, Decision Trees, Naive Bayesian models, K-means, Gaus-
sian Mixture models and general algorithms such as Ransac, Cross-validation and Grid-Search for machine-
learning applications. This package is part of the Accord.NET Framework.
• DiffSharp for machine learning and optimization applications. Operations can be nested to any level, meaning
that you can compute exact higher-order derivatives and differentiate functions that are internally making use of
differentiation, for applications such as hyperparameter optimization.
• Vulpes - Deep belief and deep learning implementation written in F# and leverages CUDA GPU execution with
Alea.cuBase.
• Encog - An advanced neural network and machine learning framework. Encog contains classes to create a wide
variety of networks, as well as support classes to normalize and process data for these neural networks. Encog
trains using multithreaded resilient propagation. Encog can also make use of a GPU to further speed processing
time. A GUI based workbench is also provided to help model and train neural networks.
• Neural Network Designer - DBMS management system and designer for neural networks. The designer appli-
cation is developed using WPF, and is a user interface which allows you to design your neural network, query
the network, create and configure chat bots that are capable of asking questions and learning from your feed
back. The chat bots can even scrape the internet for information to return in their output as well as to use for
learning.
• Infer.NET - Infer.NET is a framework for running Bayesian inference in graphical models. One can use In-
fer.NET to solve many different kinds of machine learning problems, from standard problems like classifica-
tion, recommendation or clustering through to customised solutions to domain-specific problems. Infer.NET has
been used in a wide variety of domains including information retrieval, bioinformatics, epidemiology, vision,
and many others.
• numl - numl is a machine learning library intended to ease the use of using standard modeling techniques for
both prediction and clustering.
• Math.NET Numerics - Numerical foundation of the Math.NET project, aiming to provide methods and algo-
rithms for numerical computations in science, engineering and every day use. Supports .Net 4.0, .Net 3.5 and
Mono on Windows, Linux and Mac; Silverlight 5, WindowsPhone/SL 8, WindowsPhone 8.1 and Windows 8
with PCL Portable Profiles 47 and 344; Android/iOS with Xamarin.
• Sho to enable fast and flexible prototyping. The environment includes powerful and efficient libraries for lin-
ear algebra as well as data visualization that can be used from any .NET language, as well as a feature-rich
interactive shell for rapid development.

24.16 Objective C
• YCML.
• MLPNeuralNet - Fast multilayer perceptron neural network library for iOS and Mac OS X. MLPNeuralNet
predicts new examples by trained neural network. It is built on top of the Apple’s Accelerate Framework, using
vectorized operations and hardware acceleration if available.
• MAChineLearning - An Objective-C multilayer perceptron library, with full support for training through back-
propagation. Implemented using vDSP and vecLib, it’s 20 times faster than its Java equivalent. Includes sample
code for use from Swift.
• BPN-NeuralNetwork. This network can be used in products recommendation, user behavior analysis, data
mining and data analysis.
• Multi-Perceptron-NeuralNetwork and designed unlimited-hidden-layers.
• KRHebbian-Algorithm in neural network of Machine Learning.
• KRKmeans-Algorithm - It implemented K-Means the clustering and classification algorithm. It could be used
in data mining and image compression.
• KRFuzzyCMeans-Algorithm the fuzzy clustering / classification algorithm on Machine Learning. It could be
used in data mining and image compression.
24.17 OCaml
• Oml - A general statistics and machine learning library.

• GPR - Efficient Gaussian Process Regression in OCaml.
• Libra-Tk - Algorithms for learning and inference with discrete probabilistic models.
• TensorFlow - OCaml bindings for TensorFlow.
24.18 PHP
• jieba-php - Chinese Words Segmentation Utilities.
• PHP-ML - Machine Learning library for PHP. Algorithms, Cross Validation, Neural Network, Preprocessing,
Feature Extraction and much more in one library.
• PredictionBuilder - A library for machine learning that builds predictions using a linear regression.
• Rubix ML - A high-level machine learning and deep learning library for the PHP language.
24.16. Objective C 189

24.19 Python
Computer Vision
• Scikit-Image - A collection of algorithms for image processing in Python.

• SimpleCV - An open source computer vision framework that gives access to several high-powered computer
vision libraries, such as OpenCV. Written on Python and runs on Mac, Windows, and Ubuntu Linux.
• Vigranumpy - Python bindings for the VIGRA C++ computer vision library.
• OpenFace - Free and open source face recognition with deep neural networks.
• PCV - Open source Python module for computer vision
• NLTK - A leading platform for building Python programs to work with human language data.
• Pattern - A web mining module for the Python programming language. It has tools for natural language pro-
cessing, machine learning, among others.
• Quepy - A python framework to transform natural language questions to queries in a database query language
• TextBlob tasks. Stands on the giant shoulders of NLTK and Pattern, and plays nicely with both.
• YAlign - A sentence aligner, a friendly tool for extracting parallel sentences from comparable corpora.
• jieba - Chinese Words Segmentation Utilities.
• SnowNLP - A library for processing Chinese text.
• spammy - A library for email Spam filtering built on top of nltk
• loso - Another Chinese segmentation library.
• genius - A Chinese segment base on Conditional Random Field.
• KoNLPy - A Python package for Korean natural language processing.
• nut - Natural language Understanding Toolkit
• Rosetta
• BLLIP Parser
• PyNLPl](https://github.com/proycon/pynlpl) - Python Natural Language Processing Library. General purpose
NLP library for Python. Also contains some specific modules for parsing common NLP formats, most notably
for [FoLiA, but also ARPA language models, Moses phrasetables, GIZA++ alignments.
• python-ucto
• python-frog
• python-zpar](https://github.com/EducationalTestingService/python-zpar) - Python bindings for [ZPar, a statisti-
cal part-of-speech-tagger, constiuency parser, and dependency parser for English.
• colibri-core - Python binding to C++ library for extracting and working with with basic linguistic constructions
such as n-grams and skipgrams in a quick and memory-efficient way.
• spaCy - Industrial strength NLP with Python and Cython.
• PyStanfordDependencies - Python interface for converting Penn Treebank trees to Stanford Dependencies.
• Distance - Levenshtein and Hamming distance computation

• Fuzzy Wuzzy - Fuzzy String Matching in Python

• jellyfish - a python library for doing approximate and phonetic matching of strings.
• editdistance - fast implementation of edit distance
• textacy - higher-level NLP built on Spacy
• stanford-corenlp-python](https://github.com/dasmith/stanford-corenlp-python) - Python wrapper for [Stanford
CoreNLP
• auto_ml - Automated machine learning for production and analytics. Lets you focus on the fun parts of ML,
while outputting production-ready code, and detailed analytics of your dataset and results. Includes support for
NLP, XGBoost, LightGBM, and soon, deep learning.
• machine learning](https://github.com/jeff1evesque/machine-learning) - automated build consisting of a
[web-interface](https://github.com/jeff1evesque/machine-learning#web-interface), and set of [programmatic-
interface, are stored into a NoSQL datastore.
• XGBoost Library
• Bayesian Methods for Hackers - Book/iPython notebooks on Probabilistic Programming in Python
• Featureforge A set of tools for creating and testing machine learning features, with a scikit-learn compatible
API
• scikit-learn - A Python module for machine learning built on top of SciPy.
• metric-learn - A Python module for metric learning.
• SimpleAI Python implementation of many of the artificial intelligence algorithms described on the book “Ar-
tificial Intelligence, a Modern Approach”. It focuses on providing an easy to use, well documented and tested
library.
• astroML - Machine Learning and Data Mining for Astronomy.
• graphlab-create implemented on top of a disk-backed DataFrame.
• BigML - A library that contacts external servers.
• pattern - Web mining module for Python.
• NuPIC - Numenta Platform for Intelligent Computing.
• Pylearn2](https://github.com/lisa-lab/pylearn2) - A Machine Learning library based on [Theano.
• keras](https://github.com/fchollet/keras) - Modular neural network library based on [Theano.
• Lasagne - Lightweight library to build and train neural networks in Theano.
• hebel - GPU-Accelerated Deep Learning Library in Python.
• Chainer - Flexible neural network framework
• prohpet - Fast and automated time series forecasting framework by Facebook.
• gensim - Topic Modelling for Humans.
• topik - Topic modelling toolkit
24.19. Python 191

• PyBrain - Another Python Machine Learning Library.

• Brainstorm - Fast, flexible and fun neural networks. This is the successor of PyBrain.
• Crab - A exible, fast recommender engine.
• python-recsys - A Python library for implementing a Recommender System.
• thinking bayes - Book on Bayesian Analysis
• Image-to-Image Translation with Conditional Adversarial Networks](https://github.com/williamFalcon/pix2pix-
keras) - Implementation of image to image (pix2pix) translation from the paper by [isola et al.[DEEP
LEARNING]
• Restricted Boltzmann Machines -Restricted Boltzmann Machines in Python. [DEEP LEARNING]
• Bolt - Bolt Online Learning Toolbox
• CoverTree - Python implementation of cover trees, near-drop-in replacement for scipy.spatial.kdtree
• nilearn - Machine learning for NeuroImaging in Python
• imbalanced-learn - Python module to perform under sampling and over sampling with various techniques.
• Shogun - The Shogun Machine Learning Toolbox
• Pyevolve - Genetic algorithm framework.
• Caffe - A deep learning framework developed with cleanliness, readability, and speed in mind.
• breze - Theano based library for deep and recurrent neural networks
• pyhsmm, focusing on the Bayesian Nonparametric extensions, the HDP-HMM and HDP-HSMM, mostly with
weak-limit approximations.
• mrjob - A library to let Python program run on Hadoop.
• SKLL - A wrapper around scikit-learn that makes it simpler to conduct experiments.
• neurolab - https://github.com/zueve/neurolab
• Spearmint - Spearmint is a package to perform Bayesian optimization according to the algorithms outlined in
the paper: Practical Bayesian Optimization of Machine Learning Algorithms. Jasper Snoek, Hugo Larochelle
and Ryan P. Adams. Advances in Neural Information Processing Systems, 2012.
• Pebl - Python Environment for Bayesian Learning
• Theano - Optimizing GPU-meta-programming code generating array oriented optimizing math compiler in
Python
• TensorFlow - Open source software library for numerical computation using data flow graphs
• yahmm - Hidden Markov Models for Python, implemented in Cython for speed and efficiency.
• python-timbl - A Python extension module wrapping the full TiMBL C++ programming interface. Timbl is an
elaborate k-Nearest Neighbours machine learning toolkit.
• deap - Evolutionary algorithm framework.
• pydeep - Deep Learning In Python
• mlxtend - A library consisting of useful tools for data science and machine learning tasks.
• neon](https://github.com/NervanaSystems/neon) - Nervana’s [high-performance Python-based Deep Learning
framework [DEEP LEARNING]
facilitate drop-in replacement of grid search.

• Neural Networks and Deep Learning - Code samples for my book “Neural Networks and Deep Learning” [DEEP
LEARNING]
• Annoy - Approximate nearest neighbours implementation
• skflow - Simplified interface for TensorFlow, mimicking Scikit Learn.
• TPOT - Tool that automatically creates and optimizes machine learning pipelines using genetic programming.
Consider it your personal data science assistant, automating a tedious part of machine learning.
• pgmpy A python library for working with Probabilistic Graphical Models.
• DIGITS is a web application for training deep learning models.
• Orange - Open source data visualization and data analysis for novices and experts.
• milk - Machine learning toolkit focused on supervised classification.
• TFLearn - Deep learning library featuring a higher-level API for TensorFlow.
• REP - an IPython-based environment for conducting data-driven research in a consistent and reproducible way.
REP is not trying to substitute scikit-learn, but extends it and provides better user experience.
• rgf_python Library.
• gym - OpenAI Gym is a toolkit for developing and comparing reinforcement learning algorithms.
• skbayes - Python package for Bayesian Machine Learning with scikit-learn API
• fuku-ml - Simple machine learning library, including Perceptron, Regression, Support Vector Machine, Decision
Tree and more, it’s easy to use and easy to learn for beginners.
• SciPy - A Python-based ecosystem of open-source software for mathematics, science, and engineering.
• NumPy - A fundamental package for scientific computing with Python.
• Numba complier to LLVM aimed at scientific Python by the developers of Cython and NumPy.
• NetworkX - A high-productivity software for complex networks.
• igraph - binding to igraph library - General purpose graph library
• Pandas - A library providing high-performance, easy-to-use data structures and data analysis tools.
• Open Mining
• PyMC - Markov Chain Monte Carlo sampling toolkit.
• zipline - A Pythonic algorithmic trading library.
• PyDy - Short for Python Dynamics, used to assist with workflow in the modeling of dynamic motion based
around NumPy, SciPy, IPython, and matplotlib.
• SymPy - A Python library for symbolic mathematics.
• statsmodels - Statistical modeling and econometrics in Python.
• astropy - A community Python library for Astronomy.
• matplotlib - A Python 2D plotting library.
• bokeh - Interactive Web Plotting for Python.
24.19. Python 193

• plotly - Collaborative web plotting for Python and matplotlib.

• vincent - A Python to Vega translator.
• d3py](https://github.com/mikedewar/d3py) - A plotting library for Python, based on [D3.js.
• PyDexter - Simple plotting for Python. Wrapper for D3xterjs; easily render charts in-browser.
• ggplot - Same API as ggplot2 for R.
• ggfortify - Unified interface to ggplot2 popular R packages.
• Kartograph.py - Rendering beautiful SVG maps in Python.
• pygal - A Python SVG Charts Creator.
• PyQtGraph - A pure-python graphics and GUI library built on PyQt4 / PySide and NumPy.
• pycascading
• Petrel - Tools for writing, submitting, debugging, and monitoring Storm topologies in pure Python.
• Blaze - NumPy and Pandas interface to Big Data.
• emcee - The Python ensemble sampling toolkit for affine-invariant MCMC.
• windML - A Python Framework for Wind Energy Analysis and Prediction
• vispy - GPU-based high-performance interactive OpenGL 2D/3D data visualization library
• cerebro2 A web-based visualization and debugging platform for NuPIC.
• NuPIC Studio An all-in-one NuPIC Hierarchical Temporal Memory visualization and debugging super-tool!
• SparklingPandas
• Seaborn - A python visualization library based on matplotlib
• bqplot
• pastalog - Simple, realtime visualization of neural network training performance.
• caravel - A data exploration platform designed to be visual, intuitive, and interactive.
• Dora - Tools for exploratory data analysis in Python.
• Ruffus - Computation Pipeline library for python.
• SOMPY.
• somoclu Massively parallel self-organizing maps: accelerate training on multicore CPUs, GPUs, and clusters,
has python API.
• HDBScan - implementation of the hdbscan algorithm in Python - used for clustering
• visualize_ML - A python package for data exploration and data analysis.
• scikit-plot - A visualization library for quick and easy generation of common plots in data analysis and machine
learning.
Neural networks
• Neural networks - NeuralTalk is a Python+numpy project for learning Multimodal Recurrent Neural Networks
that describe images with sentences.
• Neuron neural networks learned with Gradient descent or LeLevenberg–Marquardt algorithm.

• Data Driven Code - Very simple implementation of neural networks for dummies in python without using any
libraries, with detailed comments.
24.20 Ruby
• Treat - Text REtrieval and Annotation Toolkit, definitely the most comprehensive toolkit I’ve encountered so far
for Ruby
• Ruby Linguistics - Linguistics is a framework for building linguistic utilities for Ruby objects in any language. It
includes a generic language-independent front end, a module for mapping language codes into language names,
and a module which contains various English-language utilities.
• Stemmer - Expose libstemmer_c to Ruby
• Ruby Wordnet - This library is a Ruby interface to WordNet
• Raspel - raspell is an interface binding for ruby
• UEA Stemmer - Ruby port of UEALite Stemmer - a conservative stemmer for search and indexing
• Twitter-text-rb - A library that does auto linking and extraction of usernames, lists and hashtags in tweets
• Ruby Machine Learning - Some Machine Learning algorithms, implemented in Ruby

• Machine Learning Ruby
• jRuby Mahout - JRuby Mahout is a gem that unleashes the power of Apache Mahout in the world of JRuby.
• CardMagic-Classifier - A general classifier module to allow Bayesian and other types of classifications.
• rb-libsvm - Ruby language bindings for LIBSVM which is a Library for Support Vector Machines
• Random Forester - Creates Random Forest classifiers from PMML files
• rsruby - Ruby - R bridge

• data-visualization-ruby - Source code and supporting content for my Ruby Manor presentation on Data Visual-
isation with Ruby
• ruby-plot - gnuplot wrapper for ruby, especially for plotting roc curves into svg files
• plot-rb - A plotting library in Ruby built on top of Vega and D3.
• scruffy - A beautiful graphing toolkit for Ruby
• SciRuby
• Glean - A data management tool for humans
• Bioruby
• Arel
24.20. Ruby 195

Misc
• Big Data For Chimps

• Listof](https://github.com/kevincobain2000/listof) - Community based data collection, packed in gem. Get list
of pretty much anything (stop words, countries, non words) in txt, json or hash. [Demo/Search for a list
24.21 Rust
• deeplearn-rs - deeplearn-rs provides simple networks that use matrix multiplication, addition, and ReLU under
the MIT license.
• rustlearn - a machine learning framework featuring logistic regression, support vector machines, decision trees
and random forests.
• rusty-machine - a pure-rust machine learning library.
• leaf](https://github.com/autumnai/leaf) - open source framework for machine intelligence, sharing concepts
from TensorFlow and Caffe. Available under the MIT license. [**[Deprecated]**
• RustNN - RustNN is a feedforward neural network library.
24.22 R
• ahaz - ahaz: Regularization for semiparametric additive hazards regression

• arules - arules: Mining Association Rules and Frequent Itemsets
• biglasso - biglasso: Extending Lasso Model Fitting to Big Data in R
• bigrf - bigrf: Big Random Forests: Classification and Regression Forests for Large Data Sets
• ‘bigRR <http://cran.r-project.org/web/packages/bigRR/index.html) - bigRR: Generalized Ridge Regres-
sion (with special advantage for p >> n cases>‘__
• bmrm - bmrm: Bundle Methods for Regularized Risk Minimization Package
• Boruta - Boruta: A wrapper algorithm for all-relevant feature selection
• bst - bst: Gradient Boosting
• C50 - C50: C5.0 Decision Trees and Rule-Based Models
• caret - Classification and Regression Training: Unified interface to ~150 ML algorithms in R.
• caretEnsemble - caretEnsemble: Framework for fitting multiple caret models as well as creating ensembles of
such models.
• Clever Algorithms For Machine Learning
• CORElearn - CORElearn: Classification, regression, feature evaluation and ordinal evaluation
• CoxBoost - CoxBoost: Cox models by likelihood based boosting for a single survival endpoint or competing
risks
• Cubist - Cubist: Rule- and Instance-Based Regression Modeling

• e1071, TU Wien
• earth - earth: Multivariate Adaptive Regression Spline Models
• elasticnet - elasticnet: Elastic-Net for Sparse Estimation and Sparse PCA
• ElemStatLearn - ElemStatLearn: Data sets, functions and examples from the book: “The Elements of Statistical
Learning, Data Mining, Inference, and Prediction” by Trevor Hastie, Robert Tibshirani and Jerome Friedman
Prediction” by Trevor Hastie, Robert Tibshirani and Jerome Friedman
• evtree - evtree: Evolutionary Learning of Globally Optimal Trees
• forecast - forecast: Timeseries forecasting using ARIMA, ETS, STLM, TBATS, and neural network models
• forecastHybrid - forecastHybrid: Automatic ensemble and cross validation of ARIMA, ETS, STLM, TBATS,
and neural network models from the “forecast” package
• fpc - fpc: Flexible procedures for clustering
• frbs - frbs: Fuzzy Rule-based Systems for Classification and Regression Tasks
• GAMBoost - GAMBoost: Generalized linear and additive models by likelihood based boosting
• gamboostLSS - gamboostLSS: Boosting Methods for GAMLSS
• gbm - gbm: Generalized Boosted Regression Models
• glmnet - glmnet: Lasso and elastic-net regularized generalized linear models
• glmpath - glmpath: L1 Regularization Path for Generalized Linear Models and Cox Proportional Hazards Model
• GMMBoost - GMMBoost: Likelihood-based Boosting for Generalized mixed models
• grplasso - grplasso: Fitting user specified models with Group Lasso penalty
• grpreg - grpreg: Regularization paths for regression models with grouped covariates
• h2o - A framework for fast, parallel, and distributed machine learning algorithms at scale – Deeplearning,
Random forests, GBM, KMeans, PCA, GLM
• hda - hda: Heteroscedastic Discriminant Analysis
• Introduction to Statistical Learning
• ipred - ipred: Improved Predictors
• kernlab - kernlab: Kernel-based Machine Learning Lab
• klaR - klaR: Classification and visualization
• lars - lars: Least Angle Regression, Lasso and Forward Stagewise
• lasso2 - lasso2: L1 constrained estimation aka ‘lasso’
• LiblineaR - LiblineaR: Linear Predictive Models Based On The Liblinear C/C++ Library
• LogicReg - LogicReg: Logic Regression
• Machine Learning For Hackers
• maptree - maptree: Mapping, pruning, and graphing tree models
• mboost - mboost: Model-Based Boosting
• medley - medley: Blending regression models, using a greedy stepwise approach
• mlr - mlr: Machine Learning in R
• mvpart - mvpart: Multivariate partitioning
24.22. R 197
• ncvreg - ncvreg: Regularization paths for SCAD- and MCP-penalized regression models
• nnet - nnet: Feed-forward Neural Networks and Multinomial Log-Linear Models
• oblique.tree - oblique.tree: Oblique Trees for Classification Data
• pamr - pamr: Pam: prediction analysis for microarrays
• party - party: A Laboratory for Recursive Partytioning
• partykit - partykit: A Toolkit for Recursive Partytioning
• penalized penalized estimation in GLMs and in the Cox model
• penalizedLDA - penalizedLDA: Penalized classification using Fisher’s linear discriminant
• penalizedSVM - penalizedSVM: Feature Selection SVM using penalty functions
• quantregForest - quantregForest: Quantile Regression Forests
• randomForest - randomForest: Breiman and Cutler’s random forests for classification and regression
• randomForestSRC
• rattle - rattle: Graphical user interface for data mining in R
• rda - rda: Shrunken Centroids Regularized Discriminant Analysis
• rdetools in Feature Spaces
• REEMtree Data
• relaxo - relaxo: Relaxed Lasso
• rgenoud - rgenoud: R version of GENetic Optimization Using Derivatives
• rgp - rgp: R genetic programming framework
• Rmalschains in R
• rminer in classification and regression
• ROCR - ROCR: Visualizing the performance of scoring classifiers
• RoughSets - RoughSets: Data Analysis Using Rough Set and Fuzzy Rough Set Theories
• rpart - rpart: Recursive Partitioning and Regression Trees
• RPMM - RPMM: Recursively Partitioned Mixture Model
• RSNNS
• RWeka - RWeka: R/Weka interface
• RXshrink - RXshrink: Maximum Likelihood Shrinkage via Generalized Ridge or Least Angle Regression
• sda - sda: Shrinkage Discriminant Analysis and CAT Score Variable Selection
• SDDA - SDDA: Stepwise Diagonal Discriminant Analysis
• SuperLearner](https://github.com/ecpolley/SuperLearner) and [subsemble - Multi-algorithm ensemble learning
packages.
• svmpath - svmpath: svmpath: the SVM Path algorithm
• tgp - tgp: Bayesian treed Gaussian process models
• tree - tree: Classification and regression trees
• varSelRF - varSelRF: Variable selection using random forests

• XGBoost.R Library
facilitate drop-in replacement of grid search. Optunity is written in Python but interfaces seamlessly to R.
• igraph - binding to igraph library - General purpose graph library
• TDSP-Utilities.
• ggplot2 - A data visualization package based on the grammar of graphics.
24.23 SAS
• Enterprise Miner - Data mining and machine learning that creates deployable models using a GUI or code.
• Factory Miner - Automatically creates deployable machine learning models across numerous market or customer
segments using a GUI.
• SAS/STAT - For conducting advanced statistical analysis.

• University Edition - FREE! Includes all SAS packages necessary for data analysis and visualization, and in-
cludes online SAS courses.
High Performance Machine Learning
• High Performance Data Mining - Data mining and machine learning that creates deployable models using a GUI
or code in an MPP environment, including Hadoop.
• High Performance Text Mining - Text mining using a GUI or code in an MPP environment, including Hadoop.
• Contextual Analysis - Add structure to unstructured text using a GUI.

• Sentiment Analysis - Extract sentiment from text using a GUI.
• Text Miner - Text mining using a GUI or code.
Demos and Scripts
• ML_Tables - Concise cheat sheets containing machine learning best practices.

• enlighten-apply - Example code and materials that illustrate applications of SAS machine learning techniques.
24.23. SAS 199

• enlighten-integration - Example code and materials that illustrate techniques for integrating SAS with other
analytics technologies in Java, PMML, Python and R.
• enlighten-deep - Example code and materials that illustrate using neural networks with several hidden layers in
SAS.
• dm-flow - Library of SAS Enterprise Miner process flow diagrams to help you learn by example about specific
data mining topics.
24.24 Scala
• ScalaNLP - ScalaNLP is a suite of machine learning and numerical computing libraries.

• Breeze - Breeze is a numerical processing library for Scala.
• Chalk - Chalk is a natural language processing library.
• FACTORIE - FACTORIE is a toolkit for deployable probabilistic modeling, implemented as a software library in
Scala. It provides its users with a succinct language for creating relational factor graphs, estimating parameters
and performing inference.

• Scalding - A Scala API for Cascading
• Summing Bird - Streaming MapReduce with Scalding and Storm
• Algebird - Abstract Algebra for Scala
• xerial - Data management utilities for Scala
• simmer - Reduce your data. A unix filter for algebird-powered aggregation.
• PredictionIO - PredictionIO, a machine learning server for software developers and data engineers.
• BIDMat - CPU and GPU-accelerated matrix library intended to support large-scale exploratory data analysis.
• Wolfe Declarative Machine Learning
• Flink - Open source platform for distributed stream and batch data processing.
• Spark Notebook - Interactive and Reactive Data Science using Scala and Spark.
• Conjecture - Scalable Machine Learning in Scalding

• brushfire - Distributed decision tree ensemble learning in Scala
• ganitha - scalding powered machine learning
• adam - A genomics processing engine and specialized file format built using Apache Avro, Apache Spark and
Parquet. Apache 2 licensed.

• bioscala - Bioinformatics for the Scala programming language

• BIDMach - CPU and GPU-accelerated Machine Learning Library.
• Figaro - a Scala library for constructing probabilistic models.
• H2O Sparkling Water - H2O and Spark interoperability.
• FlinkML in Apache Flink - Distributed machine learning library in Flink
• DynaML - Scala Library/REPL for Machine Learning Research
• Saul - Flexible Declarative Learning-Based Programming.
• SwiftLearner - Simply written algorithms to help study ML or write your own implementations.
24.25 Swift
• Swift AI - Highly optimized artificial intelligence and machine learning library written in Swift.
• BrainCore - The iOS and OS X neural network framework
• swix - A bare bones library that includes a general matrix language and wraps some OpenCV for iOS develop-
ment.
• DeepLearningKit an Open Source Deep Learning Framework for Apple’s iOS, OS X and tvOS. It currently
allows using deep convolutional neural network models trained in Caffe on Apple operating systems.
• AIToolbox - A toolbox framework of AI modules written in Swift: Graphs/Trees, Linear Regression, Support
Vector Machines, Neural Networks, PCA, KMeans, Genetic Algorithms, MDP, Mixture of Gaussians.
• MLKit - A simple Machine Learning Framework written in Swift. Currently features Simple Linear Regression,
Polynomial Regression, and Ridge Regression.
• Swift Brain - The first neural network / machine learning library written in Swift. This is a project for AI
algorithms in Swift for iOS and OS X development. This project includes algorithms focused on Bayes theorem,
neural networks, SVMs, Matrices, etc..
24.25. Swift 201


CHAPTER 25
Papers
• Deep Learning
– Understanding
– Optimization / Training Techniques
– Unsupervised / Generative Models
– Image Segmentation / Object Detection
– Image / Video
– Natural Language Processing
– Speech / Other
– Reinforcement Learning
– New papers
– Classic Papers
25.1 Machine Learning
25.2 Deep Learning
Forked from terryum’s awesome deep learning papers.
203
25.2.1 Understanding
• Distilling the knowledge in a neural network (2015), G. Hinton et al. [pdf]

• Deep neural networks are easily fooled: High confidence predictions for unrecognizable images (2015), A.
Nguyen et al. [pdf]
• How transferable are features in deep neural networks? (2014), J. Yosinski et al. [pdf]
• CNN features off-the-Shelf: An astounding baseline for recognition (2014), A. Razavian et al. [pdf]
• Learning and transferring mid-Level image representations using convolutional neural networks (2014), M.
Oquab et al. [pdf]
• Visualizing and understanding convolutional networks (2014), M. Zeiler and R. Fergus [pdf]
• Decaf: A deep convolutional activation feature for generic visual recognition (2014), J. Donahue et al. [pdf]
25.2.2 Optimization / Training Techniques
• Batch normalization: Accelerating deep network training by reducing internal covariate shift (2015), S. Loffe
and C. Szegedy [pdf]
• Delving deep into rectifiers: Surpassing human-level performance on imagenet classification (2015), K. He et
al. [pdf]
• Dropout: A simple way to prevent neural networks from overfitting (2014), N. Srivastava et al. [pdf]
• Adam: A method for stochastic optimization (2014), D. Kingma and J. Ba [pdf]
• Improving neural networks by preventing co-adaptation of feature detectors (2012), G. Hinton et al. [pdf]
• Random search for hyper-parameter optimization (2012) J. Bergstra and Y. Bengio [pdf]
25.2.3 Unsupervised / Generative Models
• Pixel recurrent neural networks (2016), A. Oord et al. [pdf]

• Improved techniques for training GANs (2016), T. Salimans et al. [pdf]
• Unsupervised representation learning with deep convolutional generative adversarial networks (2015), A. Rad-
ford et al. [pdf]
• DRAW: A recurrent neural network for image generation (2015), K. Gregor et al. [pdf]
• Generative adversarial nets (2014), I. Goodfellow et al. [pdf]
• Auto-encoding variational Bayes (2013), D. Kingma and M. Welling [pdf]
• Building high-level features using large scale unsupervised learning (2013), Q. Le et al. [pdf]
25.2.4 Image Segmentation / Object Detection
• You only look once: Unified, real-time object detection (2016), J. Redmon et al. [pdf]
• Fully convolutional networks for semantic segmentation (2015), J. Long et al. [pdf]
• Faster R-CNN: Towards Real-Time Object Detection with Region Proposal Networks (2015), S. Ren et al. [pdf]
• Fast R-CNN (2015), R. Girshick [pdf]
• Rich feature hierarchies for accurate object detection and semantic segmentation (2014), R. Girshick et al. [pdf]
204 Chapter 25. Papers

• Semantic image segmentation with deep convolutional nets and fully connected CRFs, L. Chen et al. [pdf]
• Learning hierarchical features for scene labeling (2013), C. Farabet et al. [pdf]
25.2.5 Image / Video
• Image Super-Resolution Using Deep Convolutional Networks (2016), C. Dong et al. [pdf]
• A neural algorithm of artistic style (2015), L. Gatys et al. [pdf]
• Deep visual-semantic alignments for generating image descriptions (2015), A. Karpathy and L. Fei-Fei [pdf]
• Show, attend and tell: Neural image caption generation with visual attention (2015), K. Xu et al. [pdf]
• Show and tell: A neural image caption generator (2015), O. Vinyals et al. [pdf]
• Long-term recurrent convolutional networks for visual recognition and description (2015), J. Donahue et al.
[pdf]
• VQA: Visual question answering (2015), S. Antol et al. [pdf]
• DeepFace: Closing the gap to human-level performance in face verification (2014), Y. Taigman et al. [pdf]:
• Large-scale video classification with convolutional neural networks (2014), A. Karpathy et al. [pdf]
• DeepPose: Human pose estimation via deep neural networks (2014), A. Toshev and C. Szegedy [pdf]
• Two-stream convolutional networks for action recognition in videos (2014), K. Simonyan et al. [pdf]
• 3D convolutional neural networks for human action recognition (2013), S. Ji et al. [pdf]
25.2.6 Natural Language Processing
• Neural Architectures for Named Entity Recognition (2016), G. Lample et al. [pdf]
• Exploring the limits of language modeling (2016), R. Jozefowicz et al. [pdf]
• Teaching machines to read and comprehend (2015), K. Hermann et al. [pdf]
• Effective approaches to attention-based neural machine translation (2015), M. Luong et al. [pdf]
• Conditional random fields as recurrent neural networks (2015), S. Zheng and S. Jayasumana. [pdf]
• Memory networks (2014), J. Weston et al. [pdf]
• Neural turing machines (2014), A. Graves et al. [pdf]
• Neural machine translation by jointly learning to align and translate (2014), D. Bahdanau et al. [pdf]
• Sequence to sequence learning with neural networks (2014), I. Sutskever et al. [pdf]
• Learning phrase representations using RNN encoder-decoder for statistical machine translation (2014), K. Cho
et al. [pdf]
• A convolutional neural network for modeling sentences (2014), N. Kalchbrenner et al. [pdf]
• Convolutional neural networks for sentence classification (2014), Y. Kim [pdf]
• Glove: Global vectors for word representation (2014), J. Pennington et al. [pdf]
• Distributed representations of sentences and documents (2014), Q. Le and T. Mikolov [pdf]
• Distributed representations of words and phrases and their compositionality (2013), T. Mikolov et al. [pdf]
• Efficient estimation of word representations in vector space (2013), T. Mikolov et al. [pdf]
25.2. Deep Learning 205

• Recursive deep models for semantic compositionality over a sentiment treebank (2013), R. Socher et al. [pdf]
• Generating sequences with recurrent neural networks (2013), A. Graves. [pdf]
25.2.7 Speech / Other
• End-to-end attention-based large vocabulary speech recognition (2016), D. Bahdanau et al. [pdf]
• Deep speech 2: End-to-end speech recognition in English and Mandarin (2015), D. Amodei et al. [pdf]
• Speech recognition with deep recurrent neural networks (2013), A. Graves [pdf]
• Deep neural networks for acoustic modeling in speech recognition: The shared views of four research groups
(2012), G. Hinton et al. [pdf]
• Context-dependent pre-trained deep neural networks for large-vocabulary speech recognition (2012) G. Dahl et
al. [pdf]
• Acoustic modeling using deep belief networks (2012), A. Mohamed et al. [pdf]
25.2.8 Reinforcement Learning
• End-to-end training of deep visuomotor policies (2016), S. Levine et al. [pdf]

• Learning Hand-Eye Coordination for Robotic Grasping with Deep Learning and Large-Scale Data Collection
(2016), S. Levine et al. [pdf]
• Asynchronous methods for deep reinforcement learning (2016), V. Mnih et al. [pdf]
• Deep Reinforcement Learning with Double Q-Learning (2016), H. Hasselt et al. [pdf]
• Mastering the game of Go with deep neural networks and tree search (2016), D. Silver et al. [pdf]
• Continuous control with deep reinforcement learning (2015), T. Lillicrap et al. [pdf]
• Human-level control through deep reinforcement learning (2015), V. Mnih et al. [pdf]
• Deep learning for detecting robotic grasps (2015), I. Lenz et al. [pdf]
• Playing atari with deep reinforcement learning (2013), V. Mnih et al. [pdf]
25.2.9 New papers
• Deep Photo Style Transfer (2017), F. Luan et al. [pdf]

• Evolution Strategies as a Scalable Alternative to Reinforcement Learning (2017), T. Salimans et al. [pdf]
• Deformable Convolutional Networks (2017), J. Dai et al. [pdf]
• Mask R-CNN (2017), K. He et al. [pdf]
• Learning to discover cross-domain relations with generative adversarial networks (2017), T. Kim et al. [pdf]
• Deep voice: Real-time neural text-to-speech (2017), S. Arik et al., [pdf]
• PixelNet: Representation of the pixels, by the pixels, and for the pixels (2017), A. Bansal et al. [pdf]
• Batch renormalization: Towards reducing minibatch dependence in batch-normalized models (2017), S. Ioffe.
[pdf]
• Wasserstein GAN (2017), M. Arjovsky et al. [pdf]
• Understanding deep learning requires rethinking generalization (2017), C. Zhang et al. [pdf]

• Least squares generative adversarial networks (2016), X. Mao et al. [pdf]
25.2.10 Classic Papers
• An analysis of single-layer networks in unsupervised feature learning (2011), A. Coates et al. [pdf]
• Deep sparse rectifier neural networks (2011), X. Glorot et al. [pdf]
• Natural language processing (almost) from scratch (2011), R. Collobert et al. [pdf]
• Recurrent neural network based language model (2010), T. Mikolov et al. [pdf]
• Stacked denoising autoencoders: Learning useful representations in a deep network with a local denoising
criterion (2010), P. Vincent et al. [pdf]
• Learning mid-level features for recognition (2010), Y. Boureau [pdf]
• A practical guide to training restricted boltzmann machines (2010), G. Hinton [pdf]
• Understanding the difficulty of training deep feedforward neural networks (2010), X. Glorot and Y. Bengio [pdf]
• Why does unsupervised pre-training help deep learning (2010), D. Erhan et al. [pdf]
• Learning deep architectures for AI (2009), Y. Bengio. [pdf]
• Convolutional deep belief networks for scalable unsupervised learning of hierarchical representations (2009),
H. Lee et al. [pdf]
• Greedy layer-wise training of deep networks (2007), Y. Bengio et al. [pdf]
• A fast learning algorithm for deep belief nets (2006), G. Hinton et al. [pdf]
• Gradient-based learning applied to document recognition (1998), Y. LeCun et al. [pdf]
• Long short-term memory (1997), S. Hochreiter and J. Schmidhuber. [pdf]
25.2. Deep Learning 207


CHAPTER 26
Other Content
Books, blogs, courses and more forked from josephmisiti’s awesome machine learning
• Blogs
– Data Science
– Machine learning
– Math
• Books
– Machine learning
– Deep learning
– Probability & Statistics
– Linear Algebra
• Courses
• Podcasts
• Tutorials
26.1 Blogs
26.1.1 Data Science
• https://jeremykun.com/
• http://iamtrask.github.io/
• http://blog.explainmydata.com/
209
• http://andrewgelman.com/
• http://simplystatistics.org/
• http://www.evanmiller.org/
• http://jakevdp.github.io/
• http://blog.yhat.com/
• http://wesmckinney.com
• http://www.overkillanalytics.net/
• http://newton.cx/~peter/
• http://mbakker7.github.io/exploratory_computing_with_python/
• https://sebastianraschka.com/blog/index.html
• http://camdavidsonpilon.github.io/Probabilistic-Programming-and-Bayesian-Methods-for-Hackers/
• http://colah.github.io/
• http://www.thomasdimson.com/
• http://blog.smellthedata.com/
• https://sebastianraschka.com/
• http://dogdogfish.com/
• http://www.johnmyleswhite.com/
• http://drewconway.com/zia/
• http://bugra.github.io/
• http://opendata.cern.ch/
• https://alexanderetz.com/
• http://www.sumsar.net/
• https://www.countbayesie.com
• http://blog.kaggle.com/
• http://www.danvk.org/
• http://hunch.net/
• http://www.randalolson.com/blog/
• https://www.johndcook.com/blog/r_language_for_programmers/
• http://www.dataschool.io/
26.1.2 Machine learning
• OpenAI
• Distill
• Andrej Karpathy Blog
• Colah’s Blog
• WildML
210 Chapter 26. Other Content

• FastML
• TheMorningPaper
26.1.3 Math
• http://www.sumsar.net/
• http://allendowney.blogspot.ca/
• https://healthyalgorithms.com/
• https://petewarden.com/
• http://mrtz.org/blog/
26.2 Books
26.2.1 Machine learning
• Real World Machine Learning [Free Chapters]

• An Introduction To Statistical Learning - Book + R Code
• Elements of Statistical Learning - Book
• Probabilistic Programming & Bayesian Methods for Hackers - Book + IPython Notebooks
• Think Bayes - Book + Python Code
• Information Theory, Inference, and Learning Algorithms
• Gaussian Processes for Machine Learning
• Data Intensive Text Processing w/ MapReduce
• Reinforcement Learning: - An Introduction
• Mining Massive Datasets
• A First Encounter with Machine Learning
• Pattern Recognition and Machine Learning
• Machine Learning & Bayesian Reasoning
• Introduction to Machine Learning - Alex Smola and S.V.N. Vishwanathan
• A Probabilistic Theory of Pattern Recognition
• Introduction to Information Retrieval
• Forecasting: principles and practice
• Practical Artificial Intelligence Programming in Java
• Introduction to Machine Learning - Amnon Shashua
• Reinforcement Learning
• A Quest for AI
• Introduction to Applied Bayesian Statistics and Estimation for Social Scientists - Scott M. Lynch
26.2. Books 211

• Bayesian Modeling, Inference and Prediction

• A Course in Machine Learning
• Machine Learning, Neural and Statistical Classification
• Bayesian Reasoning and Machine Learning Book+MatlabToolBox
• R Programming for Data Science
• Data Mining - Practical Machine Learning Tools and Techniques Book
26.2.2 Deep learning
• Deep Learning - An MIT Press book

• Coursera Course Book on NLP
• NLTK
• NLP w/ Python
• Foundations of Statistical Natural Language Processing
• An Introduction to Information Retrieval
• A Brief Introduction to Neural Networks
• Neural Networks and Deep Learning
26.2.3 Probability & Statistics
• Think Stats - Book + Python Code

• From Algorithms to Z-Scores - Book
• The Art of R Programming
• Introduction to statistical thought
• Basic Probability Theory
• Introduction to probability - By Dartmouth College
• Principle of Uncertainty
• Probability & Statistics Cookbook
• Advanced Data Analysis From An Elementary Point of View
• Introduction to Probability - Book and course by MIT
• The Elements of Statistical Learning: Data Mining, Inference, and Prediction. -Book
• An Introduction to Statistical Learning with Applications in R - Book
• Learning Statistics Using R
• Introduction to Probability and Statistics Using R - Book
• Advanced R Programming - Book
• Practical Regression and Anova using R - Book
• R practicals - Book
• The R Inferno - Book

26.2.4 Linear Algebra
• Linear Algebra Done Wrong

• Linear Algebra, Theory, and Applications
• Convex Optimization
• Applied Numerical Computing
• Applied Numerical Linear Algebra
26.3 Courses
• CS231n, Convolutional Neural Networks for Visual Recognition, Stanford University

• CS224d, Deep Learning for Natural Language Processing, Stanford University
• Oxford Deep NLP 2017, Deep Learning for Natural Language Processing, University of Oxford
• Artificial Intelligence (Columbia University) - free
• Machine Learning (Columbia University) - free
• Machine Learning (Stanford University) - free
• Neural Networks for Machine Learning (University of Toronto) - free
• Machine Learning Specialization (University of Washington) - Courses: Machine Learning Foundations: A
Case Study Approach, Machine Learning: Regression, Machine Learning: Classification, Machine Learn-
ing: Clustering & Retrieval, Machine Learning: Recommender Systems & Dimensionality Reduction,Machine
Learning Capstone: An Intelligent Application with Deep Learning; free
• Machine Learning Course (2014-15 session) (by Nando de Freitas, University of Oxford) - Lecture slides and
video recordings.
• Learning from Data (by Yaser S. Abu-Mostafa, Caltech) - Lecture videos available
26.4 Podcasts
• The O’Reilly Data Show

• Partially Derivative
• The Talking Machines
• The Data Skeptic
• Linear Digressions
• Data Stories
• Learning Machines 101
• Not So Standard Deviations
• TWIMLAI
-_‘Machine Learning Guide <http://ocdevel.com/mlg>‘_
26.3. Courses 213

26.5 Tutorials

CHAPTER 27
Contribute
Become a contributor! Check out our github for more information.
215

ML Cheatsheet

Uploaded by

Document Informationclick to expand document information

Copyright:

Available Formats

ML Cheatsheet

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

ML Cheatsheet

Uploaded by

Copyright:

Available Formats

ML Cheatsheet Documentation

Apr 18, 2021

15 Loss Functions 111

19 Classification Algorithms 145

20 Clustering Algorithms 151

22 Reinforcement Learning 155

26 Other Content 209

𝑆𝑎𝑙𝑒𝑠 = 𝑤1 𝑅𝑎𝑑𝑖𝑜 + 𝑤2 𝑇 𝑉 + 𝑤3 𝑁 𝑒𝑤𝑠

1.2 Simple regression

Company Radio ($) Sales

1.2.1 Making predictions

𝑆𝑎𝑙𝑒𝑠 = 𝑊 𝑒𝑖𝑔ℎ𝑡 · 𝑅𝑎𝑑𝑖𝑜 + 𝐵𝑖𝑎𝑠

4 Chapter 1. Linear Regression

def predict_sales(radio, weight, bias):

1.2.2 Cost function

Given our simple linear equation 𝑦 = 𝑚𝑥 + 𝑏, we can calculate MSE as:

1.2. Simple regression 5

def cost_function(radio, sales, weight, bias):

1.2.3 Gradient descent

Using the following:

(𝑦𝑖 − (𝑚𝑥𝑖 + 𝑏))2 = 𝐴(𝐵(𝑚, 𝑏))

We can split the derivative into

6 Chapter 1. Linear Regression

And then using the Chain rule which states:

# -2(y - (mx + b))

# We subtract because the derivatives point in direction of steepest ascent

return weight, bias

1.2. Simple regression 7

def train(radio, sales, weight, bias, learning_rate, iters):

#Calculate cost for auditing purposes

return weight, bias, cost_history

1.2.5 Model evaluation

iter=1 weight=.03 bias=.0014 cost=197.25

8 Chapter 1. Linear Regression

1.2. Simple regression 9

10 Chapter 1. Linear Regression

1.2. Simple regression 11

12 Chapter 1. Linear Regression

𝑆𝑎𝑙𝑒𝑠 = .46𝑅𝑎𝑑𝑖𝑜 + .025

1.3 Multivariable regression

Company TV Radio News Units

1.3. Multivariable regression 13

1.3.1 Growing complexity

For each feature column {

14 Chapter 1. Linear Regression

We transpose the input matrix, swapping

for feature in features.T:

1.3.3 Making predictions

𝑆𝑎𝑙𝑒𝑠 = 𝑊1 𝑇 𝑉 + 𝑊2 𝑅𝑎𝑑𝑖𝑜 + 𝑊3 𝑁 𝑒𝑤𝑠𝑝𝑎𝑝𝑒𝑟

def predict(features, weights):

1.3.4 Initialize weights

1.3. Multivariable regression 15

1.3.5 Cost function

def cost_function(features, targets, weights):

predictions = predict(features, weights)

# Matrix math lets use do this without looping

# Return average squared error among predictions

1.3.6 Gradient descent