Machine Learning
Machine Learning
Machine Learning
and R Codes)
Sunil Ray, September 9, 2017
Login to Bookmark this article
Note: This article was originally published on August 10, 2015 and updated on Sept 9th, 2017
Overview
Major focus on commonly used machine learning algorithms
Algorithms covered- Linear regression, logistic regression, Naive Bayes, kNN, Random forest, etc.
Learn both theory and implementation of these algorithms in R and python
Introduction
Google’s self-driving cars and robots get a lot of press, but the company’s real future is in machine learning, the
technology that enables computers to get smarter and more personal.
We are probably living in the most defining period of human history. The period when computing moved from
large mainframes to PCs to cloud. But what makes it defining is not what has happened, but what is coming our
way in years to come.
What makes this period exciting and enthralling for someone like me is the democratization of the various tools
and techniques, which followed the boost in computing. Welcome to the world of data science!
Today, as a data scientist, I can build data-crunching machines with complex algorithms for a few dollars per
hour. But reaching here wasn’t easy! I had my dark days and nights.
Are you a beginner looking for a place to start your data science journey? Presenting two comprehensive
courses, full of knowledge and data science learning, curated just for you to learn data science (using Python)
from scratch:
The idea behind creating this guide is to simplify the journey of aspiring data scientists and machine learning
enthusiasts across the world. Through this guide, I will enable you to work on machine learning problems and
gain from experience. I am providing a high-level understanding of various machine learning algorithms
along with R & Python codes to run them. These should be sufficient to get your hands dirty.
Essentials of machine learning algorithms with implementation in R and Python
I have deliberately skipped the statistics behind these techniques, as you don’t need to understand them at the
start. So, if you are looking for statistical understanding of these algorithms, you should look elsewhere. But, if
you are looking to equip yourself to start building machine learning project, you are in for a treat.
How it works: This algorithm consist of a target / outcome variable (or dependent variable) which is to be
predicted from a given set of predictors (independent variables). Using these set of variables, we generate a
function that map inputs to desired outputs. The training process continues until the model achieves a desired
level of accuracy on the training data. Examples of Supervised Learning: Regression, Decision Tree, Random
Forest, KNN, Logistic Regression etc.
2. Unsupervised Learning
How it works: In this algorithm, we do not have any target or outcome variable to predict / estimate. It is used
for clustering population in different groups, which is widely used for segmenting customers in different groups
for specific intervention. Examples of Unsupervised Learning: Apriori algorithm, K-means.
3. Reinforcement Learning:
How it works: Using this algorithm, the machine is trained to make specific decisions. It works this way: the
machine is exposed to an environment where it trains itself continually using trial and error. This machine learns
from past experience and tries to capture the best possible knowledge to make accurate business decisions.
Example of Reinforcement Learning: Markov Decision Process
List of Common Machine Learning Algorithms
Here is the list of commonly used machine learning algorithms. These algorithms can be applied to almost any
data problem:
1. Linear Regression
2. Logistic Regression
3. Decision Tree
4. SVM
5. Naive Bayes
6. kNN
7. K-Means
8. Random Forest
9. Dimensionality Reduction Algorithms
10. Gradient Boosting algorithms
1. GBM
2. XGBoost
3. LightGBM
4. CatBoost
1. Linear Regression
It is used to estimate real values (cost of houses, number of calls, total sales etc.) based on continuous variable(s).
Here, we establish relationship between independent and dependent variables by fitting a best line. This best fit
line is known as regression line and represented by a linear equation Y= a *X + b.
The best way to understand linear regression is to relive this experience of childhood. Let us say, you ask a child
in fifth grade to arrange people in his class by increasing order of weight, without asking them their weights!
What do you think the child will do? He / she would likely look (visually analyze) at the height and build of
people and arrange them using a combination of these visible parameters. This is linear regression in real life!
The child has actually figured out that height and build would be correlated to the weight by a relationship, which
looks like the equation above.
In this equation:
Y – Dependent Variable
a – Slope
X – Independent variable
b – Intercept
These coefficients a and b are derived based on minimizing the sum of squared difference of distance between
data points and regression line.
Look at the below example. Here we have identified the best fit line having linear equation y=0.2811x+13.9. Now
using this equation, we can find the weight, knowing the height of a person.
Linear Regression is mainly of two types: Simple Linear Regression and Multiple Linear Regression. Simple
Linear Regression is characterized by one independent variable. And, Multiple Linear Regression(as the name
suggests) is characterized by multiple (more than 1) independent variables. While finding the best fit line, you
can fit a polynomial or curvilinear regression. And these are known as polynomial or curvilinear regression.
Here’s a coding window to try out your hand and build your own linear regression model in Python:
R Code
2. Logistic Regression
Don’t get confused by its name! It is a classification not a regression algorithm. It is used to estimate discrete
values ( Binary values like 0/1, yes/no, true/false ) based on given set of independent variable(s). In simple words,
it predicts the probability of occurrence of an event by fitting data to a logit function. Hence, it is also known as
logit regression. Since, it predicts the probability, its output values lies between 0 and 1 (as expected).
Let’s say your friend gives you a puzzle to solve. There are only 2 outcome scenarios – either you solve it or you
don’t. Now imagine, that you are being given wide range of puzzles / quizzes in an attempt to understand which
subjects you are good at. The outcome to this study would be something like this – if you are given a trignometry
based tenth grade problem, you are 70% likely to solve it. On the other hand, if it is grade fifth history question,
the probability of getting an answer is only 30%. This is what Logistic Regression provides you.
Coming to the math, the log odds of the outcome is modeled as a linear combination of the predictor variables.
odds= p/ (1-p) = probability of event occurrence / probability of not event occurrence
ln(odds) = ln(p/(1-p))
logit(p) = ln(p/(1-p)) = b0+b1X1+b2X2+b3X3....+bkXk
Above, p is the probability of presence of the characteristic of interest. It chooses parameters that maximize the
likelihood of observing the sample values rather than that minimize the sum of squared errors (like in ordinary
regression).
Now, you may ask, why take a log? For the sake of simplicity, let’s just say that this is one of the best mathematical
way to replicate a step function. I can go in more details, but that will beat the purpose of this article.
Build your own logistic regression model in Python here and check the accuracy:
R Code
x <- cbind(x_train,y_train)
# Train the model using the training sets and check score
logistic <- glm(y_train ~ ., data = x,family='binomial')
summary(logistic)
#Predict Output
predicted= predict(logistic,x_test)
Furthermore..
There are many different steps that could be tried in order to improve the model:
3. Decision Tree
This is one of my favorite algorithm and I use it quite frequently. It is a type of supervised learning algorithm that
is mostly used for classification problems. Surprisingly, it works for both categorical and continuous dependent
variables. In this algorithm, we split the population into two or more homogeneous sets. This is done based on
most significant attributes/ independent variables to make as distinct groups as possible. For more details, you
can read: Decision Tree Simplified.
source: statsexchange
In the image above, you can see that population is classified into four different groups based on multiple attributes
to identify ‘if they will play or not’. To split the population into different heterogeneous groups, it uses various
techniques like Gini, Information Gain, Chi-square, entropy.
The best way to understand how decision tree works, is to play Jezzball – a classic game from Microsoft (image
below). Essentially, you have a room with moving walls and you need to create walls such that maximum area
gets cleared off with out the balls.
So, every time you split the room with a wall, you are trying to create 2 different populations with in the same
room. Decision trees work in very similar fashion by dividing a population in as different groups as possible.
Let’s get our hands dirty and code our own decision tree in Python!
R Code
library(rpart)
x <- cbind(x_train,y_train)
# grow tree
fit <- rpart(y_train ~ ., data = x,method="class")
summary(fit)
#Predict Output
predicted= predict(fit,x_test)
For example, if we only had two features like Height and Hair length of an individual, we’d first plot these two
variables in two dimensional space where each point has two co-ordinates (these co-ordinates are known as
Support Vectors)
Now, we will find some line that splits the data between the two differently classified groups of data. This will
be the line such that the distances from the closest point in each of the two groups will be farthest away.
In the example shown above, the line which splits the data into two differently classified groups is the black line,
since the two closest points are the farthest apart from the line. This line is our classifier. Then, depending on
where the testing data lands on either side of the line, that’s what class we can classify the new data as.
Think of this algorithm as playing JezzBall in n-dimensional space. The tweaks in the game are:
You can draw lines/planes at any angles (rather than just horizontal or vertical as in the classic game)
The objective of the game is to segregate balls of different colors in different rooms.
And the balls are not moving.
Try your hand and design an SVM model in Python through this coding window:
R Code
library(e1071)
x <- cbind(x_train,y_train)
# Fitting model
fit <-svm(y_train ~ ., data = x)
summary(fit)
#Predict Output
predicted= predict(fit,x_test)
5. Naive Bayes
It is a classification technique based on Bayes’ theorem with an assumption of independence between predictors.
In simple terms, a Naive Bayes classifier assumes that the presence of a particular feature in a class is unrelated
to the presence of any other feature. For example, a fruit may be considered to be an apple if it is red, round, and
about 3 inches in diameter. Even if these features depend on each other or upon the existence of the other features,
a naive Bayes classifier would consider all of these properties to independently contribute to the probability that
this fruit is an apple.
Naive Bayesian model is easy to build and particularly useful for very large data sets. Along with simplicity,
Naive Bayes is known to outperform even highly sophisticated classification methods.
Bayes theorem provides a way of calculating posterior probability P(c|x) from P(c), P(x) and P(x|c). Look at the
equation below:
Here,
Example: Let’s understand it using an example. Below I have a training data set of weather and corresponding
target variable ‘Play’. Now, we need to classify whether players will play or not based on weather condition. Let’s
follow the below steps to perform it.
Step 2: Create Likelihood table by finding the probabilities like Overcast probability = 0.29 and probability of
playing is 0.64.
Step 3: Now, use Naive Bayesian equation to calculate the posterior probability for each class. The class with
the highest posterior probability is the outcome of prediction.
We can solve it using above discussed method, so P(Yes | Sunny) = P( Sunny | Yes) * P(Yes) / P (Sunny)
Here we have P (Sunny |Yes) = 3/9 = 0.33, P(Sunny) = 5/14 = 0.36, P( Yes)= 9/14 = 0.64
Now, P (Yes | Sunny) = 0.33 * 0.64 / 0.36 = 0.60, which has higher probability.
Naive Bayes uses a similar method to predict the probability of different class based on various attributes. This
algorithm is mostly used in text classification and with problems having multiple classes.
R Code
library(e1071)
x <- cbind(x_train,y_train)
# Fitting model
fit <-naiveBayes(y_train ~ ., data = x)
summary(fit)
#Predict Output
predicted= predict(fit,x_test)
These distance functions can be Euclidean, Manhattan, Minkowski and Hamming distance. First three functions
are used for continuous function and fourth one (Hamming) for categorical variables. If K = 1, then the case is
simply assigned to the class of its nearest neighbor. At times, choosing K turns out to be a challenge while
performing kNN modeling.
KNN can easily be mapped to our real lives. If you want to learn about a person, of whom you have no
information, you might like to find out about his close friends and the circles he moves in and gain access to
his/her information!
R Code
library(knn)
x <- cbind(x_train,y_train)
# Fitting model
fit <-knn(y_train ~ ., data = x,k=5)
summary(fit)
#Predict Output
predicted= predict(fit,x_test)
7. K-Means
It is a type of unsupervised algorithm which solves the clustering problem. Its procedure follows a simple and
easy way to classify a given data set through a certain number of clusters (assume k clusters). Data points inside
a cluster are homogeneous and heterogeneous to peer groups.
Remember figuring out shapes from ink blots? k means is somewhat similar this activity. You look at the shape
and spread to decipher how many different clusters / population are present!
In K-means, we have clusters and each cluster has its own centroid. Sum of square of difference between centroid
and the data points within a cluster constitutes within sum of square value for that cluster. Also, when the sum of
square values for all the clusters are added, it becomes total within sum of square value for the cluster solution.
We know that as the number of cluster increases, this value keeps on decreasing but if you plot the result you may
see that the sum of squared distance decreases sharply up to some value of k, and then much more slowly after
that. Here, we can find the optimum number of cluster.
Python Code
R Code
library(cluster)
fit <- kmeans(X, 3) # 5 cluster solution
8. Random Forest
Random Forest is a trademark term for an ensemble of decision trees. In Random Forest, we’ve collection of
decision trees (so known as “Forest”). To classify a new object based on attributes, each tree gives a classification
and we say the tree “votes” for that class. The forest chooses the classification having the most votes (over all the
trees in the forest).
1. If the number of cases in the training set is N, then sample of N cases is taken at random but with
replacement. This sample will be the training set for growing the tree.
2. If there are M input variables, a number m<<M is specified such that at each node, m variables are selected
at random out of the M and the best split on these m is used to split the node. The value of m is held
constant during the forest growing.
3. Each tree is grown to the largest extent possible. There is no pruning.
For more details on this algorithm, comparing with decision tree and tuning model parameters, I would suggest
you to read these articles:
R Code
library(randomForest)
x <- cbind(x_train,y_train)
# Fitting model
fit <- randomForest(Species ~ ., x,ntree=500)
summary(fit)
#Predict Output
predicted= predict(fit,x_test)
For example: E-commerce companies are capturing more details about customer like their demographics, web
crawling history, what they like or dislike, purchase history, feedback and many others to give them personalized
attention more than your nearest grocery shopkeeper.
As a data scientist, the data we are offered also consist of many features, this sounds good for building good
robust model but there is a challenge. How’d you identify highly significant variable(s) out 1000 or 2000? In such
cases, dimensionality reduction algorithm helps us along with various other algorithms like Decision Tree,
Random Forest, PCA, Factor Analysis, Identify based on correlation matrix, missing value ratio and others.
To know more about this algorithms, you can read “Beginners Guide To Learn Dimension Reduction
Techniques“.
Python Code
R Code
library(stats)
pca <- princomp(train, cor = TRUE)
train_reduced <- predict(pca,train)
test_reduced <- predict(pca,test)
GBM is a boosting algorithm used when we deal with plenty of data to make a prediction with high prediction
power. Boosting is actually an ensemble of learning algorithms which combines the prediction of several base
estimators in order to improve robustness over a single estimator. It combines multiple weak or average predictors
to a build strong predictor. These boosting algorithms always work well in data science competitions like Kaggle,
AV Hackathon, CrowdAnalytix.
R Code
library(caret)
x <- cbind(x_train,y_train)
# Fitting model
fitControl <- trainControl( method = "repeatedcv", number = 4, repeats = 4)
fit <- train(y ~ ., data = x, method = "gbm", trControl = fitControl,verbose = FALSE)
predicted= predict(fit,x_test,type= "prob")[,2]
GradientBoostingClassifier and Random Forest are two different boosting tree classifier and often people ask
about the difference between these two algorithms.
10.2. XGBoost
Another classic gradient boosting algorithm that’s known to be the decisive choice between winning and losing
in some Kaggle competitions.
The XGBoost has an immensely high predictive power which makes it the best choice for accuracy in events as
it possesses both linear model and the tree learning algorithm, making the algorithm almost 10x faster than
existing gradient booster techniques.
The support includes various objective functions, including regression, classification and ranking.
One of the most interesting things about the XGBoost is that it is also called a regularized boosting technique.
This helps to reduce overfit modelling and has a massive support for a range of languages such as Scala, Java,
R, Python, Julia and C++.
Supports distributed and widespread training on many machines that encompass GCE, AWS, Azure and Yarn
clusters. XGBoost can also be integrated with Spark, Flink and other cloud dataflow systems with a built in
cross validation at each iteration of the boosting process.
Python Code:
R Code:
require(caret)
x <- cbind(x_train,y_train)
# Fitting model
OR
10.3. LightGBM
LightGBM is a gradient boosting framework that uses tree based learning algorithms. It is designed to be
distributed and efficient with the following advantages:
The framework is a fast and high-performance gradient boosting one based on decision tree algorithms, used for
ranking, classification and many other machine learning tasks. It was developed under the Distributed Machine
Learning Toolkit Project of Microsoft.
Since the LightGBM is based on decision tree algorithms, it splits the tree leaf wise with the best fit whereas
other boosting algorithms split the tree depth wise or level wise rather than leaf-wise. So when growing on the
same leaf in Light GBM, the leaf-wise algorithm can reduce more loss than the level-wise algorithm and hence
results in much better accuracy which can rarely be achieved by any of the existing boosting algorithms.
Python Code:
num_round = 10
bst = lgb.train(param, train_data, num_round, valid_sets=[test_data])
bst.save_model('model.txt')
R Code:
library(RLightGBM)
data(example.binary)
#Parameters
lgbm.data.setField(handle.data, "label", y)
lgbm.booster.train(handle.booster, num_iterations, 5)
#Predict
pred <- lgbm.booster.predict(handle.booster, x.test)
#Test accuracy
sum(y.test == (y.pred > 0.5)) / length(y.test)
If you’re familiar with the Caret package in R, this is another way of implementing the LightGBM.
require(caret)
require(RLightGBM)
data(iris)
model <-caretModel.LGBM()
library(Matrix)
model.sparse <- caretModel.LGBM.sparse()
10.4. Catboost
CatBoost is a recently open-sourced machine learning algorithm from Yandex. It can easily integrate with deep
learning frameworks like Google’s TensorFlow and Apple’s Core ML.
The best part about CatBoost is that it does not require extensive data training like other ML models, and can
work on a variety of data formats; not undermining how robust it can be.
Make sure you handle missing data well before you proceed with the implementation.
Catboost can automatically deal with categorical variables without showing the type conversion error, which
helps you to focus on tuning your model better rather than sorting out trivial errors.
Learn more about Catboost from this article: https://www.analyticsvidhya.com/blog/2017/08/catboost-
automated-categorical-data/
Python Code:
import pandas as pd
import numpy as np
#Creating a training set for modeling and validation set to check model performance
X = train.drop(['Item_Outlet_Sales'], axis=1)
y = train.Item_Outlet_Sales
model.fit(X_train,
y_train,cat_features=categorical_features_indices,eval_set=(X_validation,
y_validation),plot=True)
submission = pd.DataFrame()
submission['Item_Identifier'] = test['Item_Identifier']
submission['Outlet_Identifier'] = test['Outlet_Identifier']
submission['Item_Outlet_Sales'] = model.predict(test)
R Code:
set.seed(1)
require(titanic)
require(caret)
require(catboost)
tt <- titanic::titanic_train[complete.cases(titanic::titanic_train),]
print(report)
print(importance)
Projects
Now, its time to take the plunge and actually play with some other real datasets. So are you ready to take on the
challenge? Accelerate your data science journey with the following Practice Problems:
End Notes
By now, I am sure, you would have an idea of commonly used machine learning algorithms. My sole intention
behind writing this article and providing the codes in R and Python is to get you started right away. If you are
keen to master machine learning, start right away. Take up problems, develop a physical understanding of the
process, apply these codes and see the fun!
Did you find this article useful ? Share your views and opinions in the comments section below.
8 Machine Learning Algorithms in Python – You Must Learn
by DataFlair Team · September 28, 2018
1. Objective
Previously, we discussed the techniques of machine learning with Python. Going deeper, today,
we will learn and implement 8 top Machine Learning Algorithms in Python.
a. Linear Regression
Linear regression is one of the supervised Machine learning algorithms in Python that observes
continuous features and predicts an outcome. Depending on whether it runs on a single variable or
on many features, we can call it simple linear regression or multiple linear regression.
This is one of the most popular Python ML algorithms and often under-appreciated. It assigns
optimal weights to variables to create a line ax+b to predict the output. We often use linear
regression to estimate real values like a number of calls and costs of houses based on continuous
variables. The regression line is the best line that fits Y=a*X+b to denote a relationship between
independent and dependent variables.
Do you know about Python Machine Learning Environment Setup
Let’s plot this for the diabetes dataset.
array([941.43097333])
1. >>> mean_squared_error(diabetes_y_test,diabetes_y_pred)
3035.0601152912695
0.410920728135835
1. >>> plt.plot(diabetes_X_test,diabetes_y_pred,color='pink',linewidth=3)
1. >>> plt.xticks(())
1. >>> plt.yticks(())
([], <a list of 0 Text yticklabel objects>)
1. >>> plt.show()
b. Logistic Regression
Logistic regression is a supervised classification is unique Machine Learning algorithms in Python
that finds its use in estimating discrete values like 0/1, yes/no, and true/false. This is based on a
given set of independent variables. We use a logistic function to predict the probability of an event
and this gives us an output between 0 and 1.
Although it says ‘regression’, this is actually a classification algorithm. Logistic regression fits data
into a logit function and is also called logit regression. Let’s plot this.
1. >>> x_test=np.linspace(-7,7,277)
2. >>> def model(x):
3. return 1/(1+np.exp(-x))
4. >>> loss=model(x_test*clf.coef_+clf.intercept_).ravel()
5. >>> plt.plot(x_test,loss,color='pink',linewidth=2.5)
1. >>> ols=linear_model.LinearRegression()
2. >>> ols.fit(x,y)
1. >>> plt.plot(x_test,ols.coef_*x_test+ols.intercept_,linewidth=1)
1. >>> plt.axhline(.4,color='.4')
1. >>> plt.ylabel('y')
Text(0,0.5,’y’)
1. >>> plt.xlabel('x')
Text(0.5,0,’x’)
1. >>> plt.xticks(range(-7,7))
2. >>> plt.yticks([0,0.4,1])
3. >>> plt.ylim(-.25,1.25)
(-0.25, 1.25)
1. >>> plt.xlim(-4,10)
(-4, 10)
1. >>> plt.show()
c. Decision Tree
A decision tree falls under supervised Machine Learning Algorithms in Python and comes of use
for both classification and regression- although mostly for classification. This model takes an
instance, traverses the tree, and compares important features with a determined conditional
statement. Whether it descends to the left child branch or the right depends on the result. Usually,
more important features are closer to the root.
Decision Tree, a Machine Learning algorithms in Python can work on both categorical and
continuous dependent variables. Here, we split a population into two or more homogeneous sets.
Let’s see the algorithm for this-
625
(625, 5)
01234
0B1111
1R1112
2R1113
3R1114
4R1115
1. >>> x,y,x_train,x_test,y_train,y_test=splitdataset(data)
2. >>> clf_gini=train_using_gini(x_train,x_test,y_train)
3. >>> clf_entropy=train_using_entropy(x_train,x_test,y_train)
4. >>> y_pred_gini=prediction(x_test,clf_gini)
1. >>> cal_accuracy(y_test,y_pred_gini)
1. >>> y_pred_entropy=prediction(x_test,clf_entropy)
1. >>> cal_accuracy(y_test,y_pred_entropy)
1. >>> plt.show()
Machine Learning Algorithms in Python – SVM
1. >>> for m, b, d in [(1, 0.65, 0.33), (0.5, 1.6, 0.55), (-0.2, 2.9, 0.2)]:
2. yfit = m * xfit + b
3. plt.plot(xfit, yfit, '-k')
4. plt.fill_between(xfit, yfit - d, yfit + d, edgecolor='none',
5. color='#AFFEDC', alpha=0.4)
1. >>> plt.xlim(-1,3.5)
(-1, 3.5)
1. >>> plt.show()
array([[16, 0, 0],
[ 0, 18, 0],
[ 0, 0, 11]], dtype=int64)
array([[16, 0, 0],
[ 0, 0, 18],
[ 0, 0, 11]], dtype=int64)
1. >>> logreg.predict(x)
array([0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 1, 2, 2,
2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2])
1. >>> y_pred=logreg.predict(x)
2. >>> len(y_pred)
150
0.96
1. >>> y_pred=knn.predict(x)
2. >>> metrics.accuracy_score(y,y_pred)
0.9666666666666667
1. >>> knn=KNeighborsClassifier(n_neighbors=1)
2. >>> knn.fit(x,y)
1. >>> y_pred=knn.predict(x)
2. >>> metrics.accuracy_score(y,y_pred)
1.0
(150, 4)
1. >>> y.shape
(150,)
(150, 4)
1. >>> y.shape
(150,)
(90, 4)
1. >>> x_test.shape
(60, 4)
1. >>> y_train.shape
(90,)
1. >>> y_test.shape
(60,)
1. >>> logreg=LogisticRegression()
2. >>> logreg.fit(x_train,y_train)
3. >>> y_pred=knn.predict(x_test)
4. >>> metrics.accuracy_score(y_test,y_pred)
0.9666666666666667
1. >>> knn=KNeighborsClassifier(n_neighbors=5)
2. >>> knn.fit(x_train,y_train)
1. >>> y_pred=knn.predict(x_test)
2. >>> metrics.accuracy_score(y_test,y_pred)
0.9666666666666667
1. >>> k_range=range(1,26)
2. >>> scores=[]
3. >>> for k in k_range:
4. knn = KNeighborsClassifier(n_neighbors=k)
5. knn.fit(x_train, y_train)
6. y_pred = knn.predict(x_test)
7. scores.append(metrics.accuracy_score(y_test, y_pred))
8. >>> scores
1. >>> plt.show()
g. k-Means
k-Means is an unsupervised algorithm that solves the problem of clustering. It classifies data using
a number of clusters. The data points inside a class are homogeneous and heterogeneous to peer
groups.
1. >>> centroids=kmeans.cluster_centers_
2. >>> labels=kmeans.labels_
3. >>> centroids
array([[1.16666667, 1.46666667],
[7.33333333, 9. ]])
1. >>> labels
array([0, 1, 0, 1, 0, 1])
1. >>> colors=['g.','r.','c.','y.']
2. >>> for i in range(len(x)):
3. print(x[i],labels[i])
4. plt.plot(x[i][0],x[i][1],colors[labels[i]],markersize=10)
[1. 2.] 0
[<matplotlib.lines.Line2D object at 0x0642AE10>] [5. 8.] 1
[<matplotlib.lines.Line2D object at 0x06438930>] [1.5 1.8] 0
[<matplotlib.lines.Line2D object at 0x06438BF0>] [8. 8.] 1
[<matplotlib.lines.Line2D object at 0x06438EB0>] [1. 0.6] 0
[<matplotlib.lines.Line2D object at 0x06438FB0>] [ 9. 11.] 1
[<matplotlib.lines.Line2D object at 0x043B1410>]
1. >>> plt.scatter(centroids[:,0],centroids[:,1],marker='x',s=150,linewidths=5,zorder=10)
1. >>> plt.show()
Machine Learning Algorithms in Python – K-Means
h. Random Forest
A random forest is an ensemble of decision trees. In order to classify every new object based on its
attributes, trees vote for class- each tree provides a classification. The classification with the most
votes wins in the forest.
1. >>> pl.xlabel('x')
Text(0.5,0,’x’)
1. >>> pl.ylabel('f(x)=log(x)')
Text(0,0.5,’f(x)=log(x)’)
1. >>> pl.legend(loc='best')
1. >>> pl.show()
1. >>> train,test=df[df['is_train']==True],df[df['is_train']==False]
2. >>> features=df.columns[:4]
3. >>> clf=RandomForestClassifier(n_jobs=2)
4. >>> y,_=pd.factorize(train['species'])
5. >>> clf.fit(train[features],y)
1. >>> preds=iris.target_names[clf.predict(test[features])]
2. >>> pd.crosstab(test['species'],preds,rownames=['actual'],colnames=['preds'])
So, this was all about Machine Learning Algorithms in Python Tutorial. Hope you like our
explanation.
3. Conclusion
Hence, today we discussed eight important Python Machine Learning Algorithms. Which one do
you think bears the most potential? Drop your suggestions in the comments below. We will surely
get back to you!
Related Topic- Python NumPy Tutorial
For reference
Statistical Data Analysis in Python
Share64
Tags: IPython, Jupyter, Pandas, Python, Statistical Analysis
This tutorial will introduce the use of Python for statistical data analysis, using data stored as Pandas DataFrame
objects, taking the form of a set of IPython notebooks.
Editor's note: This tutorial was originally published as course instructional material, and may contain out-of-
context references to other courses therein; this takes nothing away from the validity or usefulness of the
material.
Description
This tutorial will introduce the use of Python for statistical data analysis, using data stored as Pandas DataFrame
objects. Much of the work involved in analyzing data resides in importing, cleaning and transforming data in
preparation for analysis. Therefore, the first half of the course is comprised of a 2-part overview of basic and
intermediate Pandas usage that will show how to effectively manipulate datasets in memory. This includes tasks
like indexing, alignment, join/merge methods, date/time types, and handling of missing data. Next, we will cover
plotting and visualization using Pandas and Matplotlib, focusing on creating effective visual representations of
your data, while avoiding common pitfalls. Finally, participants will be introduced to methods for statistical data
modeling using some of the advanced functions in Numpy, Scipy and Pandas. This will include fitting your data
to probability distributions, estimating relationships among variables using linear and non-linear models, and a
brief introduction to bootstrapping methods. Each section of the tutorial will involve hands-on manipulation and
analysis of sample datasets, to be provided to attendees in advance.
The target audience for the tutorial includes all new Python users, though we recommend that users also attend
the NumPy and IPython session in the introductory track.
Student Instructions
For students familiar with Git, you may simply clone this repository to obtain all the materials (iPython
notebooks and data) for the tutorial. Alternatively, you may download a zip file containing the materials. A third
option is to simply view static notebooks by clicking on the titles of each section below.
Outline
Introduction to Pandas
Importing data
Series and DataFrame objects
Indexing, data selection and subsetting
Hierarchical indexing
Reading and writing files
Sorting and ranking
Missing data
Data summarization
Date/time types
Merging and joining DataFrame objects
Concatenation
Reshaping DataFrame objects
Pivoting
Data transformation
Permutation and sampling
Data aggregation and GroupBy operations
Statistical modeling
Fitting data to probability distributions
Fitting regression models
Model selection
Bootstrapping
Required Packages
For students running the latest version of Mac OS X (10.8), the easiest way to obtain all the packages is to
install the Scipy Superpack which works with Python 2.7.2 that ships with OS X.
Otherwise, another easy way to install all the necessary packages is to use Continuum Analytics' Anaconda.
The Ecological Detective: Confronting Models with Data, Ray Hilborn and Marc Mangel
Though targeted to ecologists, Mangel and Hilborn identify key methods that scientists can use to build useful
and credible models for their data. They don't shy away from the math, but the book is very readable and
example-laden.
Data Analysis Using Regression and Multilevel/Hierarchical Models, Andrew Gelman and Jennifer Hill
Frank Harrell's bag of tricks for regression modeling. I pull this off the shelf every week.
Bio: Christopher Fonnesbeck is an Assistant Professor in the Department of Biostatistics at the Vanderbilt
University School of Medicine. He specializes in computational statistics, Bayesian methods, meta-analysis,
and applied decision analysis. He originally hails from Vancouver, BC and received his Ph.D. from the
University of Georgia.
Getting started. Two of the most de-motivational words in the English language. The first step is often the
hardest to take, and when given too much choice in terms of direction it can often be debilitating.
Where to begin?
This post aims to take a newcomer from minimal knowledge of machine learning in Python all the way to
knowledgeable practitioner in 7 steps, all while using freely available materials and resources along the way.
The prime objective of this outline is to help you wade through the numerous free options that are available;
there are many, to be sure, but which are the best? Which complement one another? What is the best order in
which to use selected resources?
Moving forward, I make the assumption that you are not an expert in:
Machine learning
Python
Any of Python's machine learning, scientific computing, or data analysis libraries
It would probably be helpful to have some basic understanding of one or both of the first 2 topics, but even that
won't be necessary; some extra time spent on the earlier steps should help compensate.
Step 1: Basic Python Skills
If we intend to leverage Python in order to perform machine learning, having some base understanding of Python
is crucial. Fortunately, due to its widespread popularity as a general purpose programming language, as well as
its adoption in both scientific computing and machine learning, coming across beginner's tutorials is not very
difficult. Your level of experience in both Python and programming in general are crucial to choosing a starting
point.
First, you need Python installed. Since we will be using scientific computing and machine learning packages at
some point, I suggest that you install Anaconda. It is an industrial-strength Python implementation for Linux,
OSX, and Windows, complete with the required packages for machine learning, including numpy, scikit-learn,
and matplotlib. It also includes iPython Notebook, an interactive environment for many of our tutorials. I would
suggest Python 2.7, for no other reason than it is still the dominant installed version.
If you have no knowledge of programming, my suggestion is to start with the following free online book, then
move on to the subsequent materials:
If you have experience in programming but not with Python in particular, or if your Python is elementary, I
would suggest one or both of the following:
And for those looking for a 30 minute crash course in Python, here you go:
Of course, if you are an experienced Python programmer you will be able to skip this step. Even if so, I suggest
keeping the very readable Python documentation handy.
KDnuggets' own Zachary Lipton has pointed out that there is a lot of variation in what people consider a "data
scientist." This actually is a reflection of the field of machine learning, since much of what data scientists do
involves using machine learning algorithms to varying degrees. Is it necessary to intimately understand kernel
methods in order to efficiently create and gain insight from a support vector machine model? Of course not. Like
almost anything in life, required depth of theoretical understanding is relative to practical application. Gaining an
intimate understanding of machine learning algorithms is beyond the scope of this article, and generally requires
substantial amounts of time investment in a more academic setting, or via intense self-study at the very least.
The good news is that you don't need to possess a PhD-level understanding of the theoretical aspects of machine
learning in order to practice, in the same manner that not all programmers require a theoretical computer science
education in order to be effective coders.
Andrew Ng's Coursera course often gets rave reviews for its content; my suggestion, however, is to browse the
course notes compiled by a former student of the online course's previous incarnation. Skip over the Octave-
specific notes (a Matlab-like language unrelated to our Python pursuits). Be warned that these are not "official"
notes, but do seem to capture the relevant content from Andrew's course material. Of course, if you have the
time and interest, now would be the time to take Andrew Ng's Machine Learning course on Coursera.
There all sorts of video lectures out there if you prefer, alongside Ng's course mentioned above. I'm a fan of
Tom Mitchell, so here's a link to his recent lecture videos (along with Maria-Florina Balcan), which I find
particularly approachable:
You don't need all of the notes and videos at this point. A valid strategy involves moving forward to particular
exercises below, and referencing applicable sections of the above notes and videos when appropriate. For
example, when you come across an exercise implementing a regression model below, read the appropriate
regression section of Ng's notes and/or view Mitchell's regression videos at that time.
Alright. We have a handle on Python programming and understand a bit about machine learning. Beyond
Python there are a number of open source libraries generally used to facilitate practical machine learning. In
general, these are the main so-called scientific Python libraries we put to use when performing elementary
machine learning tasks (there is clearly subjectivity in this):
Scipy Lecture Notes, by Gaël Varoquaux, Emmanuelle Gouillart, and Olav Vahtras
10 Minutes to Pandas
You will see some other packages in the tutorials below, including, for example, Seaborn, which is a data
visualization library based on matplotlib. The aforementioned packages are (again, subjectively) the core of a
wide array of machine learning tasks in Python; however, understanding them should let you adapt to additional
and related packages without confusion when they are referenced in the following tutorials.
Python. Check.
Machine learning fundamentals. Check.
Numpy. Check.
Pandas. Check.
Matplotlib. Check.
The time has come. Let's start implementing machine learning algorithms with Python's de facto standard
machine learning library, scikit-learn.
Many of the following tutorials and exercises will be driven by the iPython (Jupyter) Notebook, which is an
interactive environment for executing Python. These iPython notebooks can optionally be viewed online or
downloaded and interacted with locally on your own computer.
Also note that the tutorials below are from a number of online sources. All Notebooks have been attributed to
the authors; if, for some reason, you find that someone has not been properly credited for their work, please let
me know and the situation will be rectified ASAP. In particular, I would like to tip my hat to Jake VanderPlas,
Randal Olson, Donne Martin, Kevin Markham, and Colin Raffel for their fantastic freely-available resources.
Our first tutorials for getting our feet wet with scikit-learn follow. I suggest doing all of these in order before
moving to the following steps.
A general introduction to scikit-learn, Python's most-used general purpose machine learning library, covering
the k-nearest neighbors algorithm:
A more in-depth and expanded introduction, including a starter project with a well-known dataset from start to
finish:
A focus on strategies for evaluating different models in scikit-learn, covering train/test dataset splits:
With a foundation having been laid in scikit-learn, we can move on to some more in-depth explorations of the
various common, and useful, algorithms. We start with k-means clustering, one of the most well-known
machine learning algorithms. It is a simple and often effective method for solving unsupervised learning
problems:
Next, we move back toward classification, and take a look at one of the most historically popular classification
methods:
We can then leverage regression for classification problems, via logistic regression:
We've gotten our feet wet with scikit-learn, and now we turn our attention to some more advanced topics. First
up are support vector machines, a not-necessarily-linear classifier relying on complex transformations of data
into higher dimensional space.
Next, random forests, an ensemble classifier, are examined via a Kaggle Titanic Competition walk-through:
Dimensionality reduction is a method for reducing the number of variables being considered in a problem.
Principal Component Analysis is a particular form of unsupervised dimensionality reduction:
Before moving on to the final step, we can take a moment to consider that we have come a long way in a
relatively short period of time.
Using Python and its machine learning libraries, we have covered some of the most common and well-known
machine learning algorithms (k-nearest neighbors, k-means clustering, support vector machines), investigated a
powerful ensemble technique (random forests), and examined some additional machine learning support tasks
(dimensionality reduction, model validation techniques). Along with some foundational machine learning skills,
we have started filling a useful toolkit for ourselves.
We will add one more in-demand tool to that kit before wrapping up.
Deep learning is everywhere! Deep learning builds on neural network research going back several decades, but
recent advances dating to the past several years have dramatically increased the perceived power of, and general
interest in, deep neural networks. If you are unfamiliar with deep learning, KDnuggets has many articles detailing
the numerous recent innovations, accomplishments, and accolades of the technology.
This final step does not purport to be a deep learning clinic of any sort; we will take a look at a few simple network
implementations in 2 of the leading contemporary Python deep learning libraries. For those interested in digging
deeper into deep learning, I recommend starting with the following free online book:
Theano
Theano is the first Python deep learning library we will look at. From the authors:
Theano is a Python library that allows you to define, optimize, and evaluate mathematical expressions involving
multi-dimensional arrays efficiently.
The following introductory tutorial on deep learning in Theano is lengthy, but it is quite good, very descriptive,
and heavily-commented:
Caffe
The other library we will test drive is Caffe. Again, from the authors:
Caffe is a deep learning framework made with expression, speed, and modularity in mind. It is developed by the
Berkeley Vision and Learning Center (BVLC) and by community contributors.
This tutorial is the cherry on the top of this article. While we have undertaken a few interesting examples above,
none likely compete with the following, which is implementing Google's #DeepDream using Caffe. Enjoy this
one! After understanding the tutorial, play around with it to get your processors dreaming on their own.
I didn't promise it would be quick or easy, but if you put the time in and follow the above 7 steps, there is no
reason that you won't be able to claim reasonable proficiency and understanding in a number of machine learning
algorithms and their implementation in Python using its popular libraries, including some of those on the cutting
edge of current deep learning research.
Bio: Matthew Mayo is a computer science graduate student currently working on his thesis parallelizing machine
learning algorithms. He is also a student of data mining, a data enthusiast, and an aspiring machine learning
scientist.
Related: